Transcriptome analysis of genes responding to NNV infection in Asian seabass epithelial cells

Genomic characterization and transcription analysis of European sea bass (Dicentrarchus labrax) rtp3 genes

Fish RTP3, belonging to the receptor-transporting protein family, display several functions, including a putative antiviral role as virus-responsive gene. In this work, we have identified and characterized two different European sea bass rtp3 genes. In addition, an in vivo transcription analysis in response to LPS, poly I:C and betanodavirus infection (RGNNV genotype) has been performed. The sequence analysis showed that European sea bass displays two rtp3 genes, X1 and X2, composed of two exons and a single intron (1007-bp and 888-bp long, respectively), located within the ORF sequence. The full-length cDNA is 1969 bp for rtp3 X1, and 1491 bp for rtp3 X2. Several ATTTA motifs have been found in the intron sequence of both genes, whereas rtp3 X1 also contains this motif in both untranslated regions. The transcription analyses revealed significant level of rtp3 X2 mRNA in brain and head kidney after LPS and poly I:C inoculation; however, the induction elicited by RGNNV infection was much higher, suggesting an essential role for this protein in controlling NNV infections.

Mapping of a major QTL for increased robustness and detection of genome assembly errors in Asian seabass (Lates calcarifer)

BACKGROUND

For Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either 'sensitive' or 'robust' individuals during the 28-day observation period.

RESULTS

We constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly.

CONCLUSION

According to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

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.

Cytokine-like activity of European sea bass ISG15 protein on RGNNV-infected E-11 cells

IFN-I generates an antiviral state by inducing the expression of numerous genes, called IFN-stimulated genes, ISGs, including ISG15, which is the only ISG with cytokine-like activity. In a previous study, we developed the Dl_ISG15_E11 cell line, which consisted of E11 cells able to express and secrete sea bass ISG15. The current study is a step forward, analysing the effect of secreted sea bass ISG15 on RGNNV replication in E11 cells, and looking into its immunomodulatory activity in order to corroborate its cytokine-like activity. The medium from ISG15-produccing cells compromised RGNNV replication, as it has been demonstrated both, by reduction in the viral genome synthesis and, specially, in the yield of infective viral particles. The implication of sea bass ISG15 in this protection has been demonstrated by ISG15 removal, which decreased the percentage of surviving cells upon viral infection, and by incubation of RGNNV-infected cells with a recombinant sea bass ISG15 protein, which resulted in almost full protection. Furthermore, the immunomodulatory activity of extracellular sea bass ISG15 has been demonstrated, which reaffirms a cytokine-like role for this protein.

The Insertion in the 3' UTR of Pmel17 Is the Causal Variant for Golden Skin Color in Tilapia

Understanding of the relationships between genotypes and phenotypes is a central problem in biology. Although teleosts have colorful phenotypes, not much is known about their underlying mechanisms. Our previous study showed that golden skin color in Mozambique tilapia was mapped in the major locus containing the Pmel gene, and an insertion in 3' UTR of Pmel17 was fully correlated with the golden color. However, the molecular mechanism of how Pmel17 determines the golden skin color is unknown. In this study, knockout of Pmel17 with CRISPR/Cas9 in blackish tilapias resulted in golden coloration, and rescue of Pmel17 in golden tilapias recovered the wild-type blackish color, indicating that Pmel17 is the gene determining the golden and blackish color. Functional analysis in vitro showed that the insertion in the 3' UTR of Pmel17 reduced the transcripts of Pmel17. Our data supplies more evidence to support that Pmel17 is the gene for blackish and golden colors, and highlights that the insertion in the 3' UTR of Pmel17 is the causative mutation for the golden coloration.

Genetic Basis for Resistance Against Viral Nervous Necrosis: GWAS and Potential of Genomic Prediction Explored in Farmed European Sea Bass (Dicentrarchus labrax)

Viral nervous necrosis (VNN) is an infectious disease caused by the red-spotted grouper nervous necrosis virus (RGNNV) in European sea bass and is considered a serious concern for the aquaculture industry with fry and juveniles being highly susceptible. To understand the genetic basis for resistance against VNN, a survival phenotype through the challenge test against the RGNNV was recorded in populations from multiple year classes (YC2016 and YC2017). A total of 4,851 individuals from 181 families were tested, and a subset (n∼1,535) belonging to 122 families was genotyped using a ∼57K Affymetrix Axiom array. The survival against the RGNNV showed low to moderate heritability with observed scale estimates of 0.18 and 0.25 obtained using pedigree vs. genomic information, respectively. The genome-wide association analysis showed a strong signal of quantitative trait loci (QTL) at LG12 which explained ∼33% of the genetic variance. The QTL region contained multiple genes (ITPK1, PLK4, HSPA4L, REEP1, CHMP2, MRPL35, and SCUBE) with HSPA4L and/or REEP1 genes being highly relevant with a likely effect on host response in managing disease-associated symptoms. The results on the accuracy of predicting breeding values presented 20–43% advantage in accuracy using genomic over pedigree-based information which varied across model types and applied validation schemes.

Zebrafish as a Vertebrate Model for Studying Nodavirus Infections

Nervous necrosis virus (NNV) is a neurotropic pathogenic virus affecting a multitude of marine and freshwater fish species that has a high economic impact on aquaculture farms worldwide. Therefore, the development of new tools and strategies aimed at reducing the mortality caused by this virus is a pivotal need. Although zebrafish is not considered a natural host for NNV, the numerous experimental advantages of this species make zebrafish an attractive model for studying different aspects of the disease caused by NNV, viral encephalopathy and retinopathy (VER). In this work, we established the best way and age to infect zebrafish larvae with NNV, obtaining significant mortalities in 3-day-postfertilization larvae when the virus was inoculated directly into the brain or by intramuscular microinjection. As occurs in naturally susceptible fish species, we confirmed that after intramuscular injection the virus was able to migrate to the central nervous system (CNS). As expected, due to the severe damage that this virus causes to the CNS, alterations in the swimming behavior of the zebrafish larvae were also observed. Taking advantage of the existence of transgenic fluorescent zebrafish lines, we were able to track the migration of different innate immune cells, mainly neutrophils, to the site of infection with NNV via the brain. However, we did not observe colocalization between the viral particles and neutrophils. RNA-Seq analysis of NNV-infected and uninfected larvae at 1, 3 and 5 days postinfection (dpi) revealed a powerful modulation of the antiviral immune response, especially at 5 dpi. We found that this response was dominated by, though not restricted to, the type I interferon system, the major defence mechanism in the innate immune response against viral pathogens. Therefore, as zebrafish larvae are able to develop the main characteristic of NNV infection and respond with an efficient immune arsenal, we confirmed the suitability of zebrafish larvae for modelling VER disease and studying different aspects of NNV pathogenesis, immune response and screening of antiviral drugs.

Transcriptomic Analysis of Fish Hosts Responses to Nervous Necrosis Virus

Nervous necrosis virus (NNV) has been responsible for mass mortalities in the aquaculture industry worldwide, with great economic and environmental impact. The present review aims to summarize the current knowledge of gene expression responses to nervous necrosis virus infection in different fish species based on transcriptomic analysis data. Four electronic databases, including PubMed, Web of Science, and SCOPUS were searched, and more than 500 publications on the subject were identified. Following the application of the appropriate testing, a total of 24 articles proved eligible for this review. NNV infection of different host species, in different developmental stages and tissues, presented in the eligible publications, are described in detail, revealing and highlighting genes and pathways that are most affected by the viral infection. Those transcriptome studies of NNV infected fish are oriented in elucidating the roles of genes/biomarkers for functions of special interest, depending on each study’s specific emphasis. This review presents a first attempt to provide an overview of universal host reaction mechanisms to viral infections, which will provide us with new perspectives to overcome NNV infection to build healthier and sustainable aquaculture systems.

Nanoplastics Increase Fish Susceptibility to Nodavirus Infection and Reduce Antiviral Immune Responses

Nanoplastics (NPs) might cause different negative effects on aquatic organisms at different biological levels, ranging from single cells to whole organisms, including cytotoxicity, reproduction, behavior or oxidative stress. However, the impact of NPs on disease resistance is almost unknown. The objective of this study was to assess whether exposure to 50 nm functionalized polystyrene NPs impacts fish susceptibility to viral diseases both in vitro and in vivo. In particular, we focused on the nervous necrosis virus (NNV), which affects many fish species, producing viral encephalopathy and retinopathy (VER), and causes great economic losses in marine aquaculture. In vitro and in vivo approaches were used. A brain cell line (SaB-1) was exposed to 1 μg mL⁻¹ of functionalized polystyrene NPs (PS-NH₂, PS-COOH) and then infected with NNV. Viral titers were increased in NP-exposed cells whilst the transcription of inflammatory and antiviral markers was lowered when compared to those cells only infected with NNV. In addition, European sea bass (Dicentrarchus labrax) juveniles were intraperitoneally injected with the same NPs and then challenged with NNV. Our results indicated that NPs increased the viral replication and clinical signs under which the fish died although the cumulate mortality was unaltered. Again, exposure to NPs produced a lowered inflammatory and antiviral response. Our results highlight that the presence of NPs might impact the infection process of NNV and fish resistance to the disease, posing an additional risk to marine organisms.

Effect of rearing density on nervous necrosis virus infection in Senegalese sole (Solea senegalensis)

Intensive fish farming at high densities results in a wide range of adverse consequences on fish welfare, including pathogen spreading, stress and increased mortality rates. In this work, we have assessed whether the survival of Senegalese sole infected with the nervous necrosis virus (NNV), a pathogen responsible for severe disease outbreaks, is affected by rearing density. Based on the different fish ratios per surface area (g cm^-2 ) and water volume (g L^-1 ), our research showed an earlier mortality onset in the tanks containing NNV-infected fish reared at medium density (MD: 0.071 g cm^-2 /5 g L^-1 ) and high density (HD: 0.142 g cm^-2 /10 g L^-1 ), as well as higher cumulative mortality values. However, transcription analysis of hsp70, gr1 and pepck genes, well-known stress biomarkers, seems to indicate that none of the challenged fish were under high stress conditions. NNV load was slightly higher both in dead and in sampled fish from MD and HD groups, and especially in the rearing water from these groups, where peaks in mortality seemed to correlate with increasing NNV load in the water. In conclusion, our results suggest that rearing NNV-infected Senegalese sole at high densities resulted in an earlier mortality onset and higher cumulative values and viral load.

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.

Linking Pedigree Information to the Gene Expression Phenotype to Understand Differential Family Survival Mechanisms in Highly Fecund Fish: A Case Study in the Larviculture of Pacific Bluefin Tuna

Mass spawning in fish culture often brings about a marked variance in family size, which can cause a reduction in effective population sizes in seed production for stock enhancement. This study reports an example of combined pedigree information and gene expression phenotypes to understand differential family survival mechanisms in early stages of Pacific bluefin tuna, Thunnus orientalis, in a mass culture tank. Initially, parentage was determined using the partial mitochondrial DNA control region sequence and 11 microsatellite loci at 1, 10, 15, and 40 days post-hatch (DPH). A dramatic proportional change in the families was observed at around 15 DPH; therefore, transcriptome analysis was conducted for the 15 DPH larvae using a previously developed oligonucleotide microarray. This analysis successfully addressed the family-specific gene expression phenotypes with 5739 differentially expressed genes and highlighted the importance of expression levels of gastric-function-related genes at the developmental stage for subsequent survival. This strategy demonstrated herein can be broadly applicable to species of interest in aquaculture to comprehend the molecular mechanism of parental effects on offspring survival, which will contribute to the optimization of breeding technologies.

Omics Strategies in Current Advancements of Infectious Fish Disease Management

Aquaculture is an important industry globally as it remains one of the significant alternatives of animal protein source supplies for humankind. Yet, the progression of this industry is being dampened by the increasing rate of fish mortality, mainly the outbreak of infectious diseases. Consequently, the regress in aquaculture ultimately results in the economy of multiple countries being affected due to the decline of product yields and marketability. By 2025, aquaculture is expected to contribute approximately 57% of fish consumption worldwide. Without a strategic approach to curb infectious diseases, the increasing demands of the aquaculture industry may not be sustainable and hence contributing to the over-fishing of wild fish. Recently, a new holistic approach that utilizes multi-omics platforms including transcriptomics, proteomics, and metabolomics is unraveling the intricate molecular mechanisms of host-pathogen interaction. This approach aims to provide a better understanding of how to improve the resistance of host species. However, no comprehensive review has been published on multi-omics strategies in deciphering fish disease etiology and molecular regulation. Most publications have only covered particular omics and no constructive reviews on various omics findings across fish species, particularly on their immune systems, have been described elsewhere. Our previous publication reviewed the integration of omics application for understanding the mechanism of fish immune response due to microbial infection. Hence, this review provides a thorough compilation of current advancements in omics strategies for fish disease management in the aquaculture industry. The discovery of biomarkers in various fish diseases and their potential advancement to complement the recent progress in combatting fish disease is also discussed in this review.

Transcriptome analysis reveals a complex response to the RGNNV/SJNNV reassortant Nervous Necrosis Virus strain in sea bream larvae

The gilthead sea bream (Sparus aurata) is a marine fish of great importance for Mediterranean aquaculture. This species has long been considered resistant to Nervous Necrosis Virus (NNV), an RNA virus that causes massive mortalities in several farmed fish animals. However, the recent appearance of RGNNV/SJNNV reassortant strains started to pose a serious threat to sea bream hatcheries, as it is able to infect larvae and juveniles of this species. While host response to NNV has been extensively studied in adult fish, little attention has been devoted to early life history stages, which are generally the most sensitive ones. Here we report for the first time a time-course RNA-seq analysis on 21-day old fish gilthead sea bream larvae experimentally infected with a RGNNV/SJNNV strain. NNV-infected and mock-infected samples were collected at four time points (6 h, 12 h, 24 h, and 48 h post infection). Four biological replicates, each consisting of five pooled larvae, were analysed for each time point and group. A large set of genes were found to be significantly regulated, especially at early time points (6 h and 12 h), with several heat shock protein encoding transcripts being up-regulated (e.g. hspa5, dnaj4, hspa9, hsc70), while many immune genes were down-regulated (e.g. myd88 and irf5 at T06, pik3r1, stat3, jak1, il12b and il6st at T12). A gene set enrichment analysis (GSEA) identified several altered pathways/processes. For instance, the formation of peroxisomes, which are important anti-viral components as well as essential for nervous system homeostasis, and the autophagy pathway were down-regulated at 6 h and 24 h post infection (hpi). Finally, two custom "reactomes" (i.e. significant gene sets observed in other studies) were defined and used. The first reactome integrated the transcriptomic response to NNV in different fish species, while the second one included all genes found to be stimulated either by interferon (IFN) or by IFN and Chikungunya virus in zebrafish. Genes in both reactomes showed predominant up-regulation at 6hpi and 12hpi and a general down-regulation at 24hpi. Such evidence suggest a certain degree of similarity between the response of sea bream and that of other fish species to NNV, while the observed down-regulation of IFN- and viral-stimulated pathways argues for a possible interference of NNV against the host response.

Scale Drop Disease Virus (SDDV) and Lates calcarifer Herpes Virus (LCHV) Coinfection Downregulate Immune-Relevant Pathways and Cause Splenic and Kidney Necrosis in Barramundi Under Commercial Farming Conditions

Marine farming of barramundi (Lates calcarifer) in Southeast Asia is currently severely affected by viral diseases. To better understand the biological implications and gene expression response of barramundi in commercial farming conditions during a disease outbreak, the presence of pathogens, comparative RNAseq, and histopathology targeting multiple organs of clinically “sick” and “healthy” juveniles were investigated. Coinfection of scale drop disease virus (SDDV) and L. calcarifer herpes virus (LCHV) were detected in all sampled fish, with higher SDDV viral loads in sick than in healthy fish. Histopathology showed that livers in sick fish often had moderate to severe abnormal fat accumulation (hepatic lipidosis), whereas the predominant pathology in the kidneys shows moderate to severe inflammation and glomerular necrosis. The spleen was the most severely affected organ, with sick fish presenting severe multifocal and coalescing necrosis. Principal component analysis (PC1 and PC2) explained 70.3% of the observed variance and strongly associated the above histopathological findings with SDDV loads and with the sick phenotypes, supporting a primary diagnosis of the fish being impacted by scale drop disease (SDD). Extracted RNA from kidney and spleen of the sick fish were also severely degraded likely due to severe inflammation and tissue necrosis, indicating failure of these organs in advanced stages of SDD. RNAseq of sick vs. healthy barramundi identified 2,810 and 556 differentially expressed genes (DEGs) in the liver and muscle, respectively. Eleven significantly enriched pathways (e.g., phagosome, cytokine-cytokine-receptor interaction, ECM-receptor interaction, neuroactive ligand-receptor interaction, calcium signaling, MAPK, CAMs, etc.) and gene families (e.g., tool-like receptor, TNF, lectin, complement, interleukin, chemokine, MHC, B and T cells, CD molecules, etc.) relevant to homeostasis and innate and adaptive immunity were mostly downregulated in sick fish. These DEGs and pathways, also previously identified in L. calcarifer as general immune responses to other pathogens and environmental stressors, suggest a failure of the clinically sick fish to cope and overcome the systemic inflammatory responses and tissue degeneration caused by SDD.

BEI Inactivated Vaccine Induces Innate and Adaptive Responses and Elicits Partial Protection upon Reassortant Betanodavirus Infection in Senegalese Sole

Nervous necrosis virus (NNV), the causative agent of viral encephalopathy and retinopathy (VER), is one of the most threatening viruses affecting marine and freshwater fish species worldwide. Senegalese sole is a promising fish species in Mediterranean aquaculture but also highly susceptible to NNV and VER outbreaks, that puts its farming at risk. The development of vaccines for aquaculture is one of best tools to prevent viral spread and sudden outbreaks, and virus inactivation is the simplest and most cost-effective method available. In this work, we have designed two inactivated vaccines based on the use of formalin or binary ethylenimine (BEI) to inactivate a reassortant NNV strain. After vaccination, the BEI-inactivated vaccine triggered the production of specific IgM-NNV antibodies and stimulated innate and adaptive immune responses at transcriptional level (rtp3, mx, mhcii and tcrb coding genes). Moreover, it partially improved survival after an NNV in vivo challenge, reducing the mid-term viral load and avoiding the down-regulation of immune response post-challenge. On the other hand, the formalin-inactivated vaccine improved the survival of fish upon infection without inducing the production of IgM-NNV antibodies and only stimulating the expression of herc4 and mhcii genes (in head-kidney and brain, respectively) during the vaccination period; this suggests that other immune-related pathways may be involved in the partial protection provoked. Although these vaccines against NNV showed encouraging results, further studies are needed to improve sole protection and to fully understand the underlying immune mechanism.

Functional-genomic analysis reveals intraspecies diversification of antiviral receptor transporter proteins in Xenopus laevis

The Receptor Transporter Protein (RTP) family is present in most, if not all jawed vertebrates. Most of our knowledge of this protein family comes from studies on mammalian RTPs, which are multi-function proteins that regulate cell-surface G-protein coupled receptor levels, influence olfactory system development, regulate immune signaling, and directly inhibit viral infection. However, mammals comprise less than one-tenth of extant vertebrate species, and our knowledge about the expression, function, and evolution of non-mammalian RTPs is limited. Here, we explore the evolutionary history of RTPs in vertebrates. We identify signatures of positive selection in many vertebrate RTP clades and characterize multiple, independent expansions of the RTP family outside of what has been described in mammals. We find a striking expansion of RTPs in the African clawed frog, Xenopus laevis, with 11 RTPs in this species as opposed to 1 to 4 in most other species. RNA sequencing revealed that most X. laevis RTPs are upregulated following immune stimulation. In functional assays, we demonstrate that at least three of these X. laevis RTPs inhibit infection by RNA viruses, suggesting that RTP homologs may serve as antiviral effectors outside of Mammalia.

Immunogene expression analysis in betanodavirus infected-Senegalese sole using an OpenArray® platform

The transcriptomic response of Senegalese sole (Solea senegalensis) triggered by two betanodaviruses with different virulence to that fish species has been assessed using an OpenArray® platform based on TaqMan™ quantitative PCR. The transcription of 112 genes per sample has been evaluated at two sampling times in two organs (head kidney and eye/brain-pooled samples). Those genes were involved in several roles or pathways, such as viral recognition, regulation of type I (IFN-1)-dependent immune responses, JAK-STAT cascade, interferon stimulated genes, protein ubiquitination, virus responsive genes, complement system, inflammatory response, other immune system effectors, regulation of T-cell proliferation, and proteolysis and apoptosis. The highly virulent isolate, wSs160.3, a wild type reassortant containing a RGNNV-type RNA1 and a SJNNV-type RNA2 segments, induced the expression of a higher number of genes in both tested organs than the moderately virulent strain, a recombinant harbouring mutations in the protruding domain of the capsid protein. The number of differentially expressed genes was higher 2 days after the infection with the wild type isolate than at 3 days post-inoculation. The wild type isolate also elicited an exacerbated interferon 1 response, which, instead of protecting sole against the infection, increases the disease severity by the induction of apoptosis and inflammation-derived immunopathology, although inflammation seems to be modulated by the complement system. Furthermore, results derived from this study suggest a potential important role for some genes with high expression after infection with the highly virulent virus, such as rtp3, sacs and isg15. On the other hand, the infection with the mutant does not induce immune response, probably due to an altered recognition by the host, which is supported by a different viral recognition pathway, involving myd88 and tbkbp1.

Differential immunogene expression profile of European sea bass (Dicentrarchus labrax, L.) in response to highly and low virulent NNV

European sea bass is highly susceptible to the nervous necrosis virus, RGNNV genotype, whereas natural outbreaks caused by the SJNNV genotype have not been recorded. The onset and severity of an infectious disease depend on pathogen virulence factors and the host immune response. The importance of RGNNV capsid protein amino acids 247 and 270 as virulence factors has been previously demonstrated in European sea bass; however, sea bass immune response against nodaviruses with different levels of virulence has been poorly characterized. Knowing the differences between the immune response against both kinds of isolates may be key to get more insight into the host mechanisms responsible for NNV virulence. For this reason, this study analyses the transcription of immunogenes differentially expressed in European sea bass inoculated with nodaviruses with different virulence: a RGNNV virus obtained by reverse genetics (rDl956), highly virulent to sea bass, and a mutated virus (Mut_247+270Dl956, RGNNV virus displaying SJNNV-type amino acids at positions 247 and 270 of the capsid protein), presenting lower virulence. This study has been performed in brain and head kidney, and the main differences between the immunogene responses triggered by both viruses have been observed in brain. The immunogene response in this organ is stronger after inoculation with the most virulent virus, and the main differences involved genes related with IFN I system, inflammatory response, cell-mediated response, and apoptosis. The lower virulence of Mut_247+270Dl956 to European sea bass can be associated with a delayed IFN I response, as well as an early and transitory inflammation and cell-mediated responses, suggesting that those can be pivotal elements in controlling the viral infection, and therefore, their functional activity could be analysed in future studies. In addition, this study supports the role of capsid amino acids at positions 247 and 270 as important determinants of RGNNV virulence to European sea bass.

Characterization of GAB3 and its association with NNV resistance in the Asian seabass

Understanding the functions of genes related to disease resistance and identifying polymorphisms in these genes are essential in molecular breeding for disease resistance. Viral nervous necrosis (VNN) is one of the major diseases in the Asian seabass, Lates calcarifer. Our previous works on QTL mapping, GWAS and cell-line transcriptome analysis of the Asian seabass after NNV challenge revealed that the gene GAB3 might be a candidate gene for VNN resistance. In this study, we cloned and characterized GAB3, and identified SNPs in the gene of the Asian seabass. The cDNA of the gene was 2165 bp, containing an ORF of 1674 bp encoding 557 amino acids. The gene consisted of 10 exons and nine introns. It was ubiquitously expressed in normal fish. An analysis of the association between two SNPs in the second intron and NNV resistance in 1035 fish descended from 43 families revealed that the two SNPs were significantly associated with VNN resistance. After NNV infection, the expression of GAB3 was significantly increased in the brain, spleen, muscle and gut, and was suppressed in the liver. The GAB3 protein was localized in the nucleus. Overexpression of GAB3 with specific GAB3-pcDNA was positively correlated to increased viral RNA and titer in NNV-infected Asian seabass cells. Our study provides new evidence to support that GAB3 may be an important gene related to NNV resistance. In addition, the SNPs provide DNA markers for the selection of candidate genes resistance to NNV at the juvenile stage of Asian seabass.

Metabolic profiles of fish nodavirus infection in vitro: RGNNV induced and exploited cellular fatty acid synthesis for virus infection

Red-spotted grouper nervous necrosis virus (RGNNV), the causative agent of viral nervous necrosis disease, has caused high mortality and heavy economic losses in marine aquaculture worldwide. However, changes in host cell metabolism during RGNNV infection remain largely unknown. Here, the global metabolic profiling during RGNNV infection and the roles of cellular fatty acid synthesis in RGNNV infection were investigated. As the infection progressed, 71 intracellular metabolites were significantly altered in RGNNV-infected cells compared with mock-infected cells. The levels of metabolites involved in amino acid biosynthesis and metabolism were significantly decreased, whereas those that correlated with fatty acid synthesis were significantly up-regulated during RGNNV infection. Among them, tryptophan and oleic acid were assessed as the most crucial biomarkers for RGNNV infection. In addition, RGNNV infection induced the formation of lipid droplets and re-localization of fatty acid synthase (FASN), indicating that RGNNV induced and required lipogenesis for viral infection. The exogenous addition of palmitic acid (PA) enhanced RGNNV infection, and the inhibition of FASN and acetyl-CoA carboxylase (ACC) significantly decreased RGNNV replication. Additionally, not only inhibition of palmitoylation and phospholipid synthesis, but also destruction of fatty acid β-oxidation significantly decreased viral replication. These data suggest that cellular fatty acid synthesis and mitochondrial β-oxidation are essential for RGNNV to complete the viral life cycle. Thus, it has been demonstrated for the first time that RGNNV infection in vitro overtook host cell metabolism and, in that process, cellular fatty acid synthesis was an essential component for RGNNV replication.

RNA-Seq analysis of European sea bass (Dicentrarchus labrax L.) infected with nodavirus reveals powerful modulation of the stress response

Nodavirus, or nervous necrosis virus (NNV), is the causative agent of viral encephalopathy and retinopathy (VER), a severe disease affecting numerous fish species worldwide. European sea bass, a cultured species of great economic importance, is highly susceptible to the disease. To better understand the response of this organism to NNV, we conducted RNA-Seq analysis of the brain and head kidney from experimentally infected and uninfected sea bass juveniles at 24 and 72 hours post-infection (hpi). Contrary to what was expected, we observed modest modulation of immune-related genes in the brain, the target organ of this virus, and some of these genes were even downregulated. However, genes involved in the stress response showed extremely high modulation. Accordingly, the genes encoding the enzymes implicated in the synthesis of cortisol were almost the only overexpressed genes in the head kidney at 24 hpi. This stress response was attenuated after 72 h in both tissues, and a progressive immune response against the virus was mounted. Moreover, experiments were conducted to determine how stress activation could impact NNV replication. Our results show the complex interplay between viral activity, the stress reaction and the immune response.

Establishment of a brain cell line (FuB-1) from mummichog (Fundulus heteroclitus) and its application to fish virology, immunity and nanoplastics toxicology

The marine fish mummichog (Fundulus heteroclitus), extensively used as research model, including in ecotoxicology, for over a century has been surpassed by other fish species. This fact may be associated with the lack of cell lines from this species, excellent models for the comprehension of fish physiology, immunology, toxicology and virology, that contribute to the reduction in the number of animals used in research. We have generated, for the first time, a brain-derived cell line from mummichog, FuB-1, and evaluated its application to the fields of fish virology, immunity and toxicology. First, FuB-1 cells show epithelial morphology and neural stem/astroglial origin. Secondly, FuB-1 cells effectively supports the replication of both spring viremia carp (SVCV) and infectious pancreatic necrosis (IPNV) viruses, but not nodavirus (NNV), indicating its potential use for fish virology. Related to this, FuB-1 cells infected with NNV up-regulate the transcription of genes related to the antiviral immune response, leading to cell resistance; while they are unaltered when infected with IPNV and SVCV, facilitating viral replication. Finally, FuB-1 cells were used for toxicological purposes and we demonstrated that exposure to either polystyrene nanoplastics (PS-100) or several human-usage pharmaceuticals are cytotoxic. Additionally, PS-100 particles increase the antioxidant catalase and glutathione S-transferase activities and decrease the total non-protein thiols in FuB-1 cells. However, PS-100 particles are able to reduce the cytotoxic effects induced by the pharmaceuticals. In conclusion, we have generated a cell line from mummichog, which might represent a valuable model for fish studies in the fields of virology, immunology and toxicology.

Betanodavirus and VER Disease: A 30-year Research Review

The outbreaks of viral encephalopathy and retinopathy (VER), caused by nervous necrosis virus (NNV), represent one of the main infectious threats for marine aquaculture worldwide. Since the first description of the disease at the end of the 1980s, a considerable amount of research has gone into understanding the mechanisms involved in fish infection, developing reliable diagnostic methods, and control measures, and several comprehensive reviews have been published to date. This review focuses on host–virus interaction and epidemiological aspects, comprising viral distribution and transmission as well as the continuously increasing host range (177 susceptible marine species and epizootic outbreaks reported in 62 of them), with special emphasis on genotypes and the effect of global warming on NNV infection, but also including the latest findings in the NNV life cycle and virulence as well as diagnostic methods and VER disease control.

The HIF1αn gene and its association with hypoxia tolerance in the Asian seabass

Hypoxia is one of the major challenges in aquaculture industry. Breeding of fish tolerant to hypoxia is an important task in genetic improvement of aquaculture species. The Asian seabass, Lates calcarifer, is an important foodfish species. We identified and characterized the hypoxia-inducible factor inhibitor (HIF1αn) gene in the Asian seabass. The full-length cDNA sequence of the HIF1αn was 3425 bp, with an ORF of 1065 bp, encoding 354 amino acids. The genomic sequence of the gene was 8667 bp in length, and contained eight exons and seven introns. Phylogenetic analysis of the gene in fish and tetrapods revealed that the HIF1αn in the Asian seabass was closely related to that of tilapia (Oreochromis niloticus). The HIF1αn was highly up-regulated in the gill, spleen and heart after 3.5-hours hypoxia treatment. We identified three SNPs in the third and fourth introns of the HIF1αn gene. The SNP (i.e. SNP 9332241 (C/T)) in the fourth intron was significantly (P < 0.01) associated with hypoxia tolerance. This SNP might be useful in selecting Asian seabass for improved tolerance to hypoxia. Our data also provide useful information for further detailed analysis of the function of the HIF1αn gene in hypoxia tolerance.

Transcriptomic analysis of Baltic cod (Gadus morhua) liver infected with Contracaecum osculatum third stage larvae indicates parasitic effects on growth and immune response

High infection levels due to third-stage larvae of the anisakid nematode Contracaecum osculatum have been documented in cod from the eastern part of the Baltic sea during the latest decades. The nematode larvae mainly infect the liver of Baltic cod and prevalence of infection has reached 100% with a mean intensity up to 80 parasites per host in certain areas and size classes. Low condition factors of the cod have been observed concomitant with the rise in parasite abundance suggesting a parasitic effect on growth parameters. To investigate any association between parasite infection and physiological status of the host we performed a comparative transcriptomic analysis of liver obtained from C. osculatum infected and non-infected cod. A total of 47,025 predicted gene models showed expression in cod liver and sequences corresponding to 2084 (4.43%) unigenes were differentially expressed in infected liver when compared to non-infected liver. Of the differentially expressed unigenes (DEGs) 1240 unigenes were up-regulated while 844 unigenes were down-regulated. The Gene Ontology (GO) enrichment analysis showed that 1304 DEGs were represented in cellular process and single-organism process, cell and cell part, binding and catalytic activity. As determined by the Kyoto Encyclopedia of Gene and Genomes (KEGG) Pathways analysis, 454 DEGs were involved in 138 pathways. Ninety-seven genes were related to metabolic pathways including carbohydrate, lipid, and amino acid metabolism. Thirteen regulated genes were playing a role in immune response such as Toll-like receptor signaling, NOD-like receptor signaling, RIG-I-like receptor signalling and thirty-six genes were associated with growth processes. This indicates that the nematode infection in Baltic cod may affect on molecular mechanisms involving metabolism, immune function and growth.

Comparative transcriptome analysis reveals the role of p53 signalling pathway during red-spotted grouper nervous necrosis virus infection in Lateolabrax japonicus brain cells

Nervous necrosis virus (NNV) is one of the fish pathogens that have caused mass mortalities of many marine and freshwater fishes in the world. To better comprehend the molecular immune mechanism of sea perch (Lateolabrax japonicus) against NNV infection, the comparative transcriptome analysis of red-spotted grouper nervous necrosis virus (RGNNV)-infected or mock-infected L. japonicus brain (LJB) cells was performed via RNA sequencing technology. Here, 1,969 up-regulated genes and 9,858 down-regulated genes, which were widely implicated in immune response pathways, were identified. Furthermore, we confirmed that p53 signalling pathway was repressed at 48 hr post-RGNNV infection, as indicated by up-regulation of Mdm2 and down-regulation of p53 and its downstream target genes, including Bax, Casp8 and CytC. Overexpression of L. japonicus p53 (Ljp53) significantly inhibited RGNNV replication and up-regulated the expression of apoptosis-related genes, whereas the down-regulation caused by pifithrin-α led to the opposite effect, suggesting Ljp53 might promote cell apoptosis to repress virus replication. Luciferase assay indicated that Ljp53 could enhance the promoter activities of zebrafish interferon (IFN)1, indicating that Ljp53 could exert its anti-RGNNV activities by enforcing the type I IFN response. This study revealed the potential antiviral role of p53 during NNV infection.

European sea bass brain DLB-1 cell line is susceptible to nodavirus: A transcriptomic study

Viral diseases are responsible for high rates of mortality and subsequent economic losses in modern aquaculture. The nervous necrosis virus (NNV) produces viral encephalopathy and retinopathy (VER), which affects the fish central nervous system. It is considered one of the most serious viral diseases in marine aquaculture, the European sea bass (Dicentrarchus labrax) being amongst the most susceptible. We have evaluated the European sea bass brain derived cell line (DLB-1) susceptibility to NNV genotypes and evaluated its transcriptomic profile. DLB-1 cells supported NNV gene transcription and replication since strains belonging to the four NNV genotypes produce cytopathic effects. Afterwards, DLB-1 cells were infected with an RGNNV strain, the one which showed the highest replication, for 12 and 72 h and an RNA-seq analysis was performed to identify potential genes involved in the host-NNV interactions. Differential expression analysis showed the up-regulation of many genes related to immunity, heat-shock proteins or apoptosis but not to proteasome or autophagy processes. These data suggest that the immune response, mainly the interferon (IFN) pathway, is not powerful enough to abrogate the infection, and cells finally suffer stress and die by apoptosis liberating infective particles. GO enrichment also revealed, for the first time, the down-regulation of terms related to brain/neuron biology indicating molecular mechanisms causing the pathogenic effect of NNV. This study opens the way to understand key elements in sea bass brain and NNV interactions.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

Transcriptome analysis is a powerful tool that enables a deep understanding of complicated physiological pathways, including immune responses. RNA sequencing (RNA-Seq)-based transcriptome analysis and various bioinformatics tools have also been used to study non-model animals, including aquaculture species for which reference genomes are not available. Rapid developments in these techniques have not only accelerated investigations into the process of pathogenic infection and defense strategies in fish, but also used to identify immunity-related genes in fish. These findings will contribute to fish immunotherapy for the prevention and treatment of bacterial infections through the design of more specific and effective immune stimulants, adjuvants, and vaccines. Until now, there has been little information regarding the universality and diversity of immune reactions against pathogenic infection in fish. Therefore, one of the aims of this paper is to introduce the RNA-Seq technique for examination of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of immune responses against bacteria, based on our previous findings in largemouth bass (Micropterus salmoides) against Nocardia seriolae, gray mullet (Mugil cephalus) against Lactococcus garvieae, orange-spotted grouper (Epinephelus coioides) against Vibrio harveyi, and koi carp (Cyprinus carpio) against Aeromonas sobria, using RNA-seq techniques. We demonstrated that only 39 differentially expressed genes (DEGs) were present in all species. However, the number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs in largemouth bass against N. seriolae, 819 DEGs in mullets against L. garvieae, 909 in groupers against V. harveyi, and 1471 in carps against A. sobria. The DEGs in different fish species were also representative of specific immune-related pathways. The results of this study will enhance our understanding of the immune responses of fish, and will aid in the development of effective vaccines, therapies, and disease-resistant strains.

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.

Comparative transcriptome analysis between the short-term stress and long-term adaptation of the Ruditapes philippinarum in response to benzo[a]pyrene

In order to monitor the pollution of polycyclic aromatic hydrocarbons (PAHs) in the seawater environment, screening biomarkers capable of monitoring PAHs is the focus of many studies. The transcriptomic profiles of the digestive gland tissue from the R. philippinarum groups after the exposure to BaP (4 μg/L) at four time points (0, 0.5, 6 and 15 days) were investigated to globally screen the key genes and pathways involved in the responses to short-term stress and long-term adaptation of BaP resistance. By comparative transcriptome analysis, 233, 282 and 58 differentially expressed genes (DEGs) were identified at 0.5 day, 6 day and 15 day (vs 0 day). The differential expression genes were related to stress response, detoxification metabolic process and innate immunity. DEGs of each group at different stages were clustered in six profiles based on gene expression pattern. Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis were used on all genes to determine the biological functions and processes. We selected Multidrug resistance protein 3 (MRP3), transcriptional regulator ATRX-like isoform X2 (ATRX) as biomarker indicator genes for short-term pollution monitoring and NADH dehydrogenase [ubiquinone] 1 (NQO1), Complement C1q-like protein 4 (C1q), Glutathione-S-transferase theta (GST), E3 ubiquitin-protein ligase (E3) for long-term pollution monitoring based on the different expression patterns and the function in detoxification and antioxidant defense system. Besides, the expression of seven genes was measured through Quantitative real-time PCR (qPCR) according to their gene expression patterns which was confirmed by the DGE analysis. Taken together, adoption of transcriptomic analysis to explore the bivalves' mRNA abundance changes and detoxification metabolic mechanism under the BaP stress at different time points can aid the development of sensitive and informed molecular endpoints for application towards ecotoxicogenomic monitoring of bivalves.

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.

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.

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

Characterization of a novel disease resistance gene rtp3 and its association with VNN disease resistance in Asian seabass

Asian seabass, an important food fish in Southeast Asia, has suffered from nervous necrosis virus (NNV) infection, resulting in massive mortality of Asian seabass larvae and enormous economic losses. Identification and characterization of disease resistance genes is important. Previous transcriptome analysis of Asians seabass epithelial cells after NNV infection revealed a highly inducible gene, receptor-transporting protein 3 (rtp3), indicating it could play an important role in Asian seabass - NNV interaction. To characterize this gene, we determined its expression pattern and subcellular localization. The rtp3 was highly induced in most examined tissues and organs of Asian seabass after NNV infection, and protein Rtp3 was localized in cytoplasm. Further association study in multiple families revealed that a microsatellite marker, (GT)ntt(GT)n, in the 3' UTR of rtp3 was significantly associated with VNN disease resistance in Asian seabass. Our results imply that rtp3 may be a novel disease resistance gene in Asian seabass. This data could improve our understanding of molecular interaction between Asian seabass and NNV, and has the potential to be applied in marker-assisted selection for disease resistance breeding in Asian seabass.

Characterization of the Transcriptome and Gene Expression of Brain Tissue in Sevenband Grouper (Hyporthodus septemfasciatus) in Response to NNV Infection

Grouper is one of the favorite sea food resources in Southeast Asia. However, the outbreaks of the viral nervous necrosis (VNN) disease due to nervous necrosis virus (NNV) infection have caused mass mortality of grouper larvae. Many aqua-farms have suffered substantial financial loss due to the occurrence of VNN. To better understand the infection mechanism of NNV, we performed the transcriptome analysis of sevenband grouper brain tissue, the main target of NNV infection. After artificial NNV challenge, transcriptome of brain tissues of sevenband grouper was subjected to next generation sequencing (NGS) using an Illumina Hi-seq 2500 system. Both mRNAs from pooled samples of mock and NNV-infected sevenband grouper brains were sequenced. Clean reads of mock and NNV-infected samples were de novo assembled and obtained 104,348 unigenes. In addition, 628 differentially expressed genes (DEGs) in response to NNV infection were identified. This result could provide critical information not only for the identification of genes involved in NNV infection, but for the understanding of the response of sevenband groupers to NNV infection.

Transcriptomic profiles of striped snakehead fish cells (SSN-1) infected with red-spotted grouper nervous necrosis virus (RGNNV) with an emphasis on apoptosis pathway

Nervous necrosis virus (NNV), the causative agent of viral nervous necrosis (VNN) disease, has caused mass mortality of cultured marine and freshwater fish worldwide, resulting in enormous economic losses in the aquaculture industry. However, the molecular mechanisms underlying the pathogenicity of NNV are still poorly understood. In this study, the transcriptomic profiles of striped snakehead fish (Channa striatus) cells (SSN-1) infected with red-spotted grouper NNV (RGNNV) were investigated using deep RNA sequencing technique. From 254,955,234 raw reads, a total of 253,338,544 clean reads were obtained and they were assembled into 93,372 unigenes. Differentially expressed genes (DEGs) were identified from RGNNV-infected or mock-infected SSN-1 cells, including 1184 up-regulated and 1456 down-regulated genes at 3 h (h) post of infection (poi), and 1138 up-regulated and 2073 down-regulated genes at 24 h poi, respectively. These DEGs were involved in many pathways related to viral pathogenesis, including retinoic acid-inducible gene I (RIG-I) like receptors pathway, apoptosis pathway, oxidative phosphorylation, PI3K-Akt signaling pathway, and MAPK signaling pathway. Subsequent analysis focusing on the apoptosis pathway showed that the expression of Endonuclease G (EndoG) was up-regulated upon RGNNV infection at both 3 and 24 h poi. Therefore, EndoG gene was cloned and its function was further characterized. The results showed that over-expression of EndoG could also induce cellular apoptosis in SSN-1 cells, indicating that RGNNV infection might induce apoptosis of SSN-1 cells via EndoG-associated mitochondrial pathway. These results will shed a new light on the pathogenesis of NNV.

MKK7 confers different activities to viral infection of Singapore grouper iridovirus (SGIV) and nervous necrosis virus (NNV) in grouper

Mitogen-activated protein kinase 7 (MKK7) is one of the major stress-activated protein kinase (SAPK)-activating kinases in response to environmental or physiological stimuli. Here a MKK7 named as Ec-MKK7 was identified from orange-spotted grouper, Epinephelus coioides. The full-length cDNA of Ec-MKK7 was 1853 bp, with an open reading frame (ORF) of 1272 bp encoding a putative protein of 423 amino acids. A characteristic S-K-A-K-T motif was contained in the domain of dual-specificity protein kinase, mitogen-activated protein kinase kinase 7 (PKc_MKK7). Intracellular localization showed that Ec-MKK7 was localized in both the cytoplasm and the nucleus of grouper spleen (GS) and/or grouper brain (EAGB) cells. Moreover, Ec-MKK7 was universally expressed in all examined tissues and showed expression modulation to challenges of lipopolysacchride (LPS), Singapore grouper iridovirus (SGIV) and polyriboinosinic polyribocytidylic acid (poly I:C) in vivo. A gene targeting strategy over-expressing Ec-MKK7 was performed to examine the activities of MKK7 to viral infection in vitro. Our data showed that Ec-MKK7 was involved in the evasion and replication of SGIV but played an antiviral role to the infection of nervous necrosis virus (NNV). All results demonstrated that Ec-MKK7 could play important roles in grouper innate immunity and show distinct functions on virus infection.

The expression and function of hsp30-like small heat shock protein genes in amphibians, birds, fish, and reptiles

Small heat shock proteins (sHSPs) are a superfamily of molecular chaperones with important roles in protein homeostasis and other cellular functions. Amphibians, reptiles, fish and birds have a shsp gene called hsp30, which was also referred to as hspb11 or hsp25 in some fish and bird species. Hsp30 genes, which are not found in mammals, are transcribed in response to heat shock or other stresses by means of the heat shock factor that is activated in response to an accumulation of unfolded protein. Amino acid sequence analysis revealed that representative HSP30s from different classes of non-mammalian vertebrates were distinct from other sHSPs including HSPB1/HSP27. Studies with amphibian and fish recombinant HSP30 determined that they were molecular chaperones since they inhibited heat- or chemically-induced aggregation of unfolded protein. During non-mammalian vertebrate development, hsp30 genes were differentially expressed in selected tissues. Also, heat shock-induced stage-specific expression of hsp30 genes in frog embryos was regulated at the level of chromatin structure. In adults and/or tissue culture cells, hsp30 gene expression was induced by heat shock, arsenite, cadmium or proteasomal inhibitors, all of which enhanced the production of unfolded/damaged protein. Finally, immunocytochemical analysis of frog and chicken tissue culture cells revealed that proteotoxic stress-induced HSP30 accumulation co-localized with aggresome-like inclusion bodies. The congregation of damaged protein in aggresomes minimizes the toxic effect of aggregated protein dispersed throughout the cell. The current availability of probes to detect the presence of hsp30 mRNA or encoded protein has resulted in the increased use of hsp30 gene expression as a marker of proteotoxic stress in non-mammalian vertebrates.