Transcriptome comparative analysis of immune tissues from asymptomatic and diseased Epinephelus moara naturally infected with nervous necrosis virus

Decoding the fish genome opens a new era in important trait research and molecular breeding in China

Aquaculture represents the fastest-growing global food production sector, as it has become an essential component of the global food supply. China has the world's largest aquaculture industry in terms of production volume. However, the sustainable development of fish culture is hindered by several concerns, including germplasm degradation and disease outbreaks. The practice of genomic breeding, which relies heavily on genome information and genotypephenotype relationships, has significant potential for increasing the efficiency of aquaculture production. In 2014, the completion of the genome sequencing and annotation of the Chinese tongue sole signified the beginning of the fish genomics era in China. Since then, domestic researchers have made dramatic progress in functional genomic studies. To date, the genomes of more than 60 species of fish in China have been assembled and annotated. Based on these reference genomes, evolutionary, comparative, and functional genomic studies have revolutionized our understanding of a wide range of biologically and economically important traits of fishes, including growth and development, sex determination, disease resistance, metamorphosis, and pigmentation. Furthermore, genomic tools and breeding techniques such as SNP arrays, genomic selection, and genome editing have greatly accelerated genetic improvement through the incorporation of functional genomic information into breeding activities. This review aims to summarize the current status, advances, and perspectives of the genome resources, genomic study of important traits, and genomic breeding techniques of fish in China. The review will provide aquaculture researchers, fish breeders, and farmers with updated information concerning fish genomic research and breeding technology. The summary will help to promote the genetic improvement of production traits and thus will support the sustainable development of fish aquaculture.

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.

Transcriptome analysis of the response of four immune related organs of tilapia (Oreochromis niloticus) to the addition of resveratrol in feed

Resveratrol (RES) has been found to have immunological enhancement effects on Oreochromis niloticus. In O. nilocticus, the liver, spleen and kidney act as immune target tissues, while intestine works for nutrition sensing organ. In the present study, we determined RES administration on these immune tissues transcriptomic response in genetically improved farmed tilapia (GIFT), and further analyzed the relationship between transcriptomic response and intestinal microbiota. As results, hepatic hemosiderin and intestinal goblet cells significantly increased with RES addition. Kyoto encyclopedia of genes and genomes (KEGG) pathways associated with herpes simplex virus 1 infection, calcium signaling pathway, cell adhesion molecules, apoptosis, and mitogen-activated protein kinase (MAPK)/peroxisome proliferators-activated receptors (PPAR) signaling pathways were enriched. In particular, the differentially enriched genes (DEGs) associated pathways were present in different sampling tissues, times, and comparisons, interestingly, the PPAR signaling pathway was enriched with increasing time of RES addition. The assembled DEGs presented verified expression in the kidney, liver, spleen, and intestine tissues, and fabp6 was highly expressed in the intestine. Serial DEGs of fatty acid-binding proteins (fabp7, fabp7a, fabp10a) decreased in the liver and kidney, and fabp6 significantly increased in the spleen. With time, the pathways of energy metabolism, glycan biosynthesis, and metabolism decreased and increased in the intestinal metagenome. Some Candidatus branches significantly increased (C. cerribacteria and C. harrisonbacteria) and while others decreased (C. glodbacteria, etc.), whereas C. verstraetearchaeota fluctuated with RES addition. slc27a6 and dbi were negatively correlated with bacteria involved in the lipid, energy, and carbohydrate metabolism pathways. The present study suggests that RES supplementation affected lipid metabolism in immune-related organs may be related to the PPAR signaling pathway.

Transcriptome and 16S rRNA analysis revealed the response of largemouth bass (Micropterus salmoides) to Rhabdovirus infection

To better understand the response of largemouth bass (Micropterus salmoides) to Micropterus salmoides rhabdovirus (MSRV) infection, we investigated the intestinal bacterial flora and transcriptome profile of fish at 72 hours post-infection (hpi). Total of 1574 differentially expressed genes (DEGs) were identified in largemouth bass spleen following MSRV infection, including 573 upregulated and 1001 downregulated genes. KEGG and GO enrichment analysis revealed that upregulated genes were enriched in certain antiviral related signaling pathway, including NOD-like receptor (NLR), RIG-I like receptors (RLR) and regulation of the interferon (IFN)-γ-mediated signaling pathway, whereas some immune-related DEGs enriched in focal adhesion (FA) and ECM-receptor interaction(ECM-RI) were downregulated, as well as genes associated with metabolic processes, such as peroxisome proliferator-activated receptors (PPAR), adipocytokine signaling pathway, Glycerolipid and Retinol metabolism. Furthermore, the principal component analysis (PCA) and phylogenetic analysis revealed that MSRV infection significantly affected the microbiota of largemouth bass intestine; the LEfSe analysis showed that relative abundances of Streptococcus were significantly increased, while the content of Akkermansia, Enterococcus and Lactobacillus were remarkably decreased in the fish intestine following MSRV infection. Additionally, a high correlation was determined between the expressions of interferon-related upregulated genes and the relative abundance of Streptococcus by redundancy analysis (RDA). These results collectively illustrated that intestinal microbiota composition might be associated with the immune-related gene expression in largemouth bass in response to MSRV infection.

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.

The dynamic immune responses of Mandarin fish (Siniperca chuatsi) to ISKNV in early infection based on full-length transcriptome analysis and weighted gene co-expression network analysis

Mandarin fish (Siniperca chuatsi) been seriously harmed by infectious spleen and kidney necrosis virus (ISKNV) in recent years, but the early immune response mechanism of infection is still unknown. Here, we performed RNA sequencing on the spleens of mandarin fish infected with ISKNV at 0, 12, 24, 48, and 72 h post-infection (hpi) using short-read Illumina RNA sequencing and long-read Pacific Biosciences isoform sequencing to generate a full-length transcriptome. The immune responses of mandarin fish infected with ISKNV at the molecular level were characterized by RNA-seq analysis and weighted gene co-expression network analysis (WGCNA). A total of 26,528 full-length transcript sequences were obtained. There were 2,729 (1,680 up-regulated and 1,112 down-regulated), 1,874 (1,136 up-regulated and 738 down-regulated), 2,032 (1,158 up-regulated and 847 down-regulated), and 4,176 (2,233 up-regulated and 1,943 down-regulated) differentially expressed genes (DEGs) in mandarin fish at 12, 24, 48, and 72 hpi, compared with uninfected fish, respectively. A total of four modules of co-expressed DEGs identified by WGCNA were significantly positively correlated to the four time points after infection, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the immune-related DEGs in all these modules were mainly enriched in Phagosome, Endocytosis, Herpes simplex infection, and Cytokine-cytokine receptor interaction pathways. Further analysis showed that oher signaling pathways, including CAMs, NOD-like receptor and ER protein processing, Intestinal immune network for IgA production, TLR pathway, and Apoptosis significantly enriched in four modules corresponding to 12, 24, 48, and 72 hpi respectively, had specifically participated in the immune response. Hub genes identified based on the high-degree nodes in the WGCN, including CAM3, IL-8, CCL21, STING, SNX1, PFR and TBK1, and some DEGs such as MHCI, MHCII, TfR, STING, TNF α, TBK1, IRF1, and NF-kB, BCR, IgA and Bcl-XL had involved in dynamic molecular response of mandarin fish to ISKNV infection. In sum, this study provides a set of full-length transcriptome of the spleen tissue of mandarin fish for the first time and revealed a group of immune genes and pathways involved in different temporal responses to ISKNV infection, which has implications for resource conservation and aiding the development of strategies to prevent virus early infection for mandarin fish.

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.

Combining of transcriptomic and proteomic data to mine immune-related genes and proteins in the liver of Cynoglossus semilaevis challenged with Vibrio anguillarum

The liver is a multi-functional organ including metabolism, substance synthesis, detoxification, and various immune functions, and its role in immunity has attracted more and more attention. However, research on the liver immune response of fish infected by pathogenic bacteria is currently lacking. In this study, the transcriptomics and proteomics of the liver of Cynoglossus semilaevis infected with Vibrio anguillarum were analyzed. A total of 1470 genes and 497 proteins were differentially expressed in the pairwise comparison of obvious symptoms of infection (HOSG), no obvious symptoms of infection (NOSG) and PBS treatment (CG). Gene ontology and KEGG enrichment pathways analysis showed that differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were mainly enriched in toll-like receptors (TLRs), complement and coagulation cascades, nucleotide oligomerization domain (NOD)-like receptors (NLRs), mitogen-activated protein kinase (MAPK) and phagosome signaling pathways, which suggested the combined action of the five pathways were significant to enhance the liver immune defense. The combination of transcriptomic and proteomic analysis showed that ITGβ1, C3, C5 and MRC1 were significantly up-regulated, which might play an important role in the liver immune response to the recognition of V. anguillarum, inflammatory response and phagocytosis. The transcriptome and proteome data we obtained provide information on some key genes and proteins for further study of the mechanism of liver immune response.

Molecular characterization of the T cell costimulatory receptors CD28 and CTLA4 in the European sea bass

For the activation of T cells, it is necessary the specific recognition of the peptide by the T cell receptors (TCR) in the surface of antigen-presenting cells (APCs) and additional signals delivered by costimulatory receptors. In fish, knowledge about the presence of these costimulatory signals is limited and functional evidence almost absent. Thus, in this study, we have identified the stimulatory CD28 and the inhibitory cytotoxic T-lymphocyte-associated protein 4 (CTLA4) coreceptors in the European sea bass (Dicentrarchus labrax), and evaluated their transcription. In parallel, the transcription encoding for the T cell markers CD8α and CD4 was also evaluated. Both coreceptors showed the canonical architecture including a signal peptide, an immunoglobulin domain, a transmembrane region and a cytosolic tail. Protein predictions and phylogenetic tree identify them as true mammalian orthologues of CD28 and CTLA4. We found these genes constitutively expressed in all studied organs of European sea bass with high expression in lymphoid organs (thymus, spleen and head-kidney) and liver. The molecular expression pattern of these genes was up-regulated in head-kidney leucocytes stimulated with T mitogens as concanavalin A and phytohemagglutinin (PHA), but not with the B cell mitogen lipopolysaccharide (LPS). Fish challenged with nodavirus (NNV) evidenced a differential and opposing regulation of the cd28 and ctla4 transcription levels in the brain, the target organ for viral replication, and head-kidney. While cd28 transcription tends to decrease over the infection time in both organs the expression of the ctla4 gene tends to increase. Interestingly, the coreceptor expression is highly and significantly correlated to the transcription of the T cell markers. Our results highlight the important role of CD28 and CTLA4 as costimulatory receptors of T cells in European sea bass but further studies are deserved.

A Chromosomal-scale Reference Genome of the Kelp Grouper Epinephelus moara

The kelp grouper, Epinephelus moara, is a carnivorous coral reef fish widely distributing in the Indo-Pacific Oceans. Compared to other grouper species, E. moara is featured for its wide tolerance to environmental stresses, such as temperature and salinity. In addition, it is popularly used as a breeding species for hybrid grouper production. Here, we report the de novo sequencing and assembly of E. moara genome using a combination of Illumina pair-ended, PacBio and Hi-C technologies. We generated a 1.08 Gb genome assembly with a scaffold N50 of 44.93 Mb and contig N50 of 2.22 Mb. The scaffolds and contigs were clustered and oriented into 24 chromosomes, and 24,135 protein-coding genes were predicted, among which 96.8% were annotated. More than 95.6% of the conserved complete genes were successfully retrieved by BUSCO analysis. Comparative genomic analyses showed that some expanded gene families in the E. moara genome were significantly enriched in innate immune pathways. The E. moara genome provides a valuable resource for genetic improvement and genomic breeding of groupers, as well as evolutionary and comparative study with other grouper species.

Transcriptome analysis reveals immune pathways underlying resistance in the common carp Cyprinus carpio against the oomycete Aphanomyces invadans

Infection with Aphanomyces invadans is a serious fish disease with major global impacts. Despite affecting over 160 fish species, some of the species like the common carp Cyprinus carpio are resistant to A. invadans infection. In the present study, we investigated the transcriptomes of head kidney of common carp experimentally infected with A. invadans. In time course analysis, 5288 genes were found to be differentially expressed (DEGs), of which 731 were involved in 21 immune pathways. The analysis of immune-related DEGs suggested that efficient processing and presentation of A. invadans antigens, enhanced phagocytosis, recognition of pathogen-associated molecular patterns, and increased recruitment of leukocytes to the sites of infection contribute to resistance of common carp against A. invadans. Herein, we provide a systematic understanding of the disease resistance mechanisms in common carp at molecular level as a valuable resource for developing disease management strategies for this devastating fish-pathogenic oomycete.

Prevalence of Infectious Spleen and Kidney Necrosis Virus (ISKNV), Nervous Necrosis Virus (NNV) and Ectoparasites in Juvenile Epinephelus spp. Farmed in Aceh, Indonesia

A cross-sectional survey was used to estimate the prevalence of infections with the Infectious spleen and kidney necrosis virus (ISKNV, Megalocytivirus), nervous necrosis virus (NNV, Betanodavirus), and infestations with ectoparasites during the rainy season in juvenile grouper (Epinephelus spp.) farmed in Aceh, Indonesia. The survey was intended to detect aquatic pathogens present at 10% prevalence with 95% confidence, assuming 100% sensitivity and specificity using a sample size of 30 for each diagnostic test. Eight populations of grouper from seven farms were sampled. Additional targeted sampling was conducted for populations experiencing high mortality. Infection with NNV was detected at all farms with seven of the eight populations being positive. The apparent prevalence for NNV ranged from 0% (95% CI: 0–12) to 73% (95% CI: 54–88). All of the fish tested from the targeted samples (Populations 9 and 10) were positive for NNV and all had vacuolation of the brain and retina consistent with viral nervous necrosis (VNN). Coinfections with ISKNV were detected in five populations, with the highest apparent prevalence being 13% (95% CI: 4–31%). Trichodina sp., Cryptocaryonirritans and Gyrodactylus sp. were detected at three farms, with 66% to 100% of fish being infested. Hybrid grouper sourced from a hatchery were 5.4 and 24.9 times more likely to have a NNV infection and a higher parasite load compared to orange-spotted grouper collected from the wild (p < 0.001). This study found that VNN remains a high-impact disease in grouper nurseries in Aceh, Indonesia.

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.

Integrated transcriptomic and proteomic analyses of grass carp intestines after vaccination with a double-targeted DNA vaccine of Vibrio mimicus

Intestinal mucosal immunity plays a vital role against Vibrio mimicus infection because it is an enteric pathogen causing serious vibriosis in fish. In the previous studies, we developed an oral double-targeted DNA vaccine of V. mimicus and demonstrated that the vaccine could elicit significantly higher intestinal mucosal immune response than did naked DNA vaccine. But, little is known underlying regulatory molecular mechanisms of the enhanced intestinal mucosal immunity. Here the transcriptome and proteome in the intestines of the grass carps immunized or not with the double-targeted DNA vaccine were investigated by using RNA-seq and iTRAQ-coupled LC-MS/MS. Compared with the control group, a total of 5339 differentially expressed genes (DEGs) and 1173 differentially expressed proteins (DEPs) were identified in the immunized fish intestines. Subsequently, the integrated analysis between transcriptome and proteome data revealed that 250 DEPs were matched with the corresponding DEGs (named associated DEPs/DEGs) at both transcriptome and proteome levels. Fifty of all the associated DEPs/DEGs were immune-related and mainly enriched in phagosome, antigen-processing and presentation, complement and coagulation cascades, NLRs and MAPK signaling pathways via Gene Ontology and KEGG pathway analyses, which suggested the coordination of the five activated pathways was essential to the enhanced intestinal mucosal immune response in the immunized fish. The protein-protein interaction analysis showed that 60 of the 63 immune-related DEPs to form an integrated network. Additionally, randomly selected DEGs and DEPs were respectively validated by quantitative real-time RT-PCR and multiple reaction monitoring (MRM) assay, indicating that the both RNA-Seq and iTRAQ results in the study were reliable. Overall, our comprehensive transcriptome and proteome data provide some key genes and their protein products for further research on the regulatory molecular mechanisms underlying the enhanced intestinal mucosal immunity.

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.