Transcriptome characterization and gene expression of Epinephelus spp in endoplasmic reticulum stress-related pathway during betanodavirus infection in vitro

Functional characterization of Malabar grouper (Epinephelus malabaricus) interferon regulatory factor 9 involved in antiviral response

IRF9 is a crucial component in the JAK-STAT pathway. IRF9 interacts with STAT1 and STAT2 to form IFN-I-stimulated gene factor 3 (ISGF3) in response to type I IFN stimulation, which promotes ISG transcription. However, the mechanism by which IFN signaling regulates Malabar grouper (Epinephelus malabaricus) IRF9 is still elusive. Here, we explored the nd tissue-specific mRNA distribution of the MgIRF9 gene, as well as its antiviral function in E. malabaricus. MgIRF9 encodes a protein of 438 amino acids with an open reading frame of 1317 base pairs. MgIRF9 mRNA was detected in all tissues of a healthy M. grouper, with the highest concentrations in the muscle, gills, and brain. It was significantly up-regulated by nervous necrosis virus infection and poly (I:C) stimulation. The gel mobility shift test demonstrated a high-affinity association between MgIRF9 and the promoter of zfIFN in vitro. In GK cells, grouper recombinant IFN-treated samples showed a significant response in ISGs and exhibited antiviral function. Subsequently, overexpression of MgIRF9 resulted in a considerable increase in IFN and ISGs mRNA expression (ADAR1, ADAR1-Like, and ADAR2). Co-immunoprecipitation studies demonstrated that MgIRF9 and STAT2 can interact in vivo. According to the findings, M. grouper IRF9 may play a role in how IFN signaling induces ISG gene expression in grouper species.

Profile of Innate Immunity in Gilthead Seabream Larvae Reflects Mortality upon Betanodavirus Reassortant Infection and Replication

Historically, gilthead seabream (Sparus aurata) has been considered a fish species resistant to nervous necrosis virus (NNV) disease. Nevertheless, mortality in seabream hatcheries, associated with typical clinical signs of the viral encephalopathy and retinopathy (VER) disease has been confirmed to be caused by RGNNV/SJNNV reassortants. Because of this, seabream larvae at 37 and 86 days post-hatching (dph) were infected by immersion with RGNNV/SJNNV and SJNNV/RGNNV reassortants under laboratory conditions, and mortality, viral replication and immunity were evaluated. Our results show that gilthead seabream larvae, mainly those at 37 dph, are susceptible to infection with both NNV reassortant genotypes, with the highest impact from the RGNNV/SJNNV reassortant. In addition, viral replication occurs at both ages (37 and 86 dph) but the recovery of infective particles was only confirmed in 37 dph larvae,; this value was also highest with the RGNNV/SJNNV reassortant. Larvae immunity, including the expression of antiviral, inflammatory and cell-mediated cytotoxicity genes, was affected by NNV infection. Levels of the natural killer lysin (Nkl) peptide were increased in SJNNV/RGNNV-infected larvae of 37 dph, though hepcidin was not. Our results demonstrate that the seabream larvae are susceptible to both NNV reassortants, though mainly to RGNNV/SJNNV, in an age-dependent manner.

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.

Comparative transcriptome analysis of two contrasting resistant and susceptible Aegilops tauschii accessions to wheat leaf rust (Puccinia triticina) using RNA-sequencing

Leaf rust, caused by Puccinia triticina Eriks., is the most common rust disease of wheat (Triticum aestivum L.) worldwide. Owing to the rapid evolution of virulent pathotypes, new and effective leaf rust resistance sources must be found. Aegilops tauschii, an excellent source of resistance genes to a wide range of diseases and pests, may provide novel routes for resistance to this disease. In this study, we aimed to elucidate the transcriptome of leaf rust resistance in two contrasting resistant and susceptible Ae. tauschii accessions using RNA-sequencing. Gene ontology, analysis of pathway enrichment and transcription factors provided an apprehensible review of differentially expressed genes and highlighted biological mechanisms behind the Aegilops–P. triticina interaction. The results showed the resistant accession could uniquely recognize pathogen invasion and respond precisely via reducing galactosyltransferase and overexpressing chromatin remodeling, signaling pathways, cellular homeostasis regulation, alkaloid biosynthesis pathway and alpha-linolenic acid metabolism. However, the suppression of photosynthetic pathway and external stimulus responses were observed upon rust infection in the susceptible genotype. In particular, this first report of comparative transcriptome analysis offers an insight into the strength and weakness of Aegilops against leaf rust and exhibits a pipeline for future wheat breeding programs.

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.

The Alteration of Intestinal Microbiota Profile and Immune Response in Epinephelus coioides during Pathogen Infection

Epinephelus coioides, or grouper, is a high economic value fish species that plays an important role in the aquaculture industry in Asia. However, both viral and bacterial diseases have threatened grouper for many years, especially nervous necrosis virus, grouper iridovirus and Vibrio harveyi, which have caused a bottleneck in the grouper industry. Currently, intestinal microbiota can provide novel insights into the pathogenesis-related factors involved in pathogen infection. Hence, we investigated the comparison of intestinal microbiota communities in control group and pathogen-infected grouper through high-throughput sequencing of the 16S rRNA gene. Our results showed that microbial diversity was decreased, whereas microbial richness was increased during pathogen infection. The individuals in each group were distributed distinctly on the PLSDA diagram, especially the GIV group. Proteobacteria and Firmicutes were the most abundant bacterial phyla in all groups. Interestingly, beneficial genera, Faecalibacterium and Bifidobacterium, predominated in the intestines of the control group. In contrast, the intestines of pathogen-infected grouper had higher levels of harmful genera such as Sphingomonas, Atopostipes, Staphylococcus and Acinetobacter. Additionally, we investigated the expression levels of innate and adaptive immune-related genes after viral and bacterial infection. The results revealed that immunoglobulin T and proinflammatory cytokine levels in the intestine increased after pathogen infection. Through these unique bacterial compositions in diseased and uninfected fish, we could establish a novel therapeutic approach and bacterial marker for preventing and controlling these diseases.

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.

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.

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.

Blast2Fish: a reference-based annotation web tool for transcriptome analysis of non-model teleost fish

BACKGROUND

Transcriptome analysis by next-generation sequencing has become a popular technique in recent years. This approach is quite suitable for non-model organism study, as de novo assembly is independent of prior genomic sequences of organisms. De novo sequencing has benefited many studies on commercially important fish species. However, to understand the functions of these assembled sequences, they still need to be annotated with existing sequence databases. By combining Basic Local Alignment Search Tool (BLAST) and Gene Ontology analysis, we were able to identify homologous sequences of assembled sequences and describe their characteristics using pre-defined tags for each gene, though the above conventional annotation results obtained for non-model assembled sequences was still associated with a lack of pre-defined tags and poorly documented records in the database.

RESULTS

We introduced Blast2Fish, a novel approach for performing functional enrichment analysis on non-model teleost fish transcriptome data. The Blast2Fish pipeline was designed to be a reference-based enrichment method. Instead of annotating the BLAST single top hit by a pre-defined gene-to-tag database, we included 500 hits to search related PubMed articles and parse biological terms. These descriptive terms were then sorted and recorded as annotations for the query. The results showed that Blast2Fish was capable of providing meaningful annotations on immunology topics for non-model fish transcriptome analysis.

CONCLUSION

Blast2Fish provides a novel approach for annotating sequences of non-model fish. The reference-based strategy allows annotation to be performed without pre-defined tags for each gene. This method strongly benefits non-model teleost fish studies for gene functional enrichment analysis.

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.

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

Epinephelus moara is an economically important fish in Southeast Asian countries but is suffering from nervous necrosis virus (NNV) infection. A deeper understanding of the host-NNV interaction mechanisms makes sense for disease control, however, at present, the pathogenesis of natural NNV infection and the resistance mechanism in host remains poorly understood. In this study, asymptomatic and diseased E. moara with clinical symptoms of viral nervous necrosis (VNN) from a grouper farm were both detected with a positive RT-PCR signal of NNV, then transcriptome sequencing of their immune tissues (liver, spleen and kidney) were performed for comparation analysis. The de novo assemblies yielded 53,789 unigenes which had a length varied from 201 to 19,675 bp and a N50 length of 2115 bp, and 29,451 unigenes were functionally annotated, with 83, 250 and 5632 unigenes being differentially expressed in liver, spleen and kidney respectively. KEGG pathway enrichment analysis of the DEGs showed many DEGs were enriched in immune related pathways. Although the expression of class I major histocompatibility complex (MHC) was significantly higher in three immune tissues of the diseased grouper, many immune related genes, including humoral immune molecules (such as antibodies), the cellular mediated cytotoxic molecules (such as perforin) and some adhesion related genes were down regulated in the diseased grouper. Our results provided many unigenes that might play important roles in NNV resistance for further research. Furthermore, a total of 8666 unigenes containing 11,623 simple sequence repeats (SSRs) were identified, which provided useful information for screening molecular markers associated with NNV resistance in E. moara.

Comparative analysis of growth performance and liver transcriptome response of juvenile Ancherythroculter nigrocauda fed diets with different protein levels

This study aimed at investigating the effects of dietary protein levels on the growth and liver transcriptome in juvenile Ancherythroculter nigrocauda. Six semi-purified diets were formulated containing 25 (control), 30, 35, 40, 45, and 50% protein. Each diet was fed to three groups of 35 fish (mean initial weight: 5.86 ± 0.10 g) for 56 days. The rate of weight gain and specific growth rate increased with dietary protein levels from 25% to 40%, but remained unchanged when fed with 45 or 50% dietary protein. The feed conversion ratio was significantly influenced by the dietary protein levels, being the lowest in fish fed 40% protein. Illumina RNA-seq analysis was performed to investigate liver gene expression changes under different dietary protein treatments. A total of 367.78 million clean reads were obtained from the six libraries. Compared with 25% protein treatment library, there were 734, 1946, 1755, 2726, and 1523 upregulated genes, and 407, 1882, 1865, 2216 and 1624 downregulated genes in the 30, 35, 40, 45, and 50% protein treatment libraries, respectively. Trend analysis of these differentially expressed genes (DEGs) identified six statistically significant trends. A series of DEGs that related to protein metabolism, growth and development, lipid metabolism and immune and stress response were identified. Moreover, gene ontology enrichment analysis of the DEGs demonstrated that cellular process, single-organism process, metabolic process and biological regulation were the most highly overrepresented biological processes. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that protein processing in endoplasmic reticulum, PPAR signaling pathway, complement and coagulation cascades, and cytochrome P450 (CYP450s) were significantly enriched in the dietary protein treatment groups. Furthermore, qPCR results showed excellent agreement on those of RNA-seq for both up- and down-regulated genes (including fasn, accα, SCD, CPT-I, igf1, ST, AST, trdmt1, hsp70, cyp450, MHC-II, C4, tgfβ, ube4b, apoE and abcb7). Thus, our results provide the baseline information for the feed formulation and nutritional research for A. nigrocauda.

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.

Comparative transcriptome analysis of two selenium-accumulating genotypes of Aegilops tauschii Coss. in response to selenium

Background

Selenium (Se), an essential micronutrient in both animals and humans, has various biological functions, and its deficiency can lead to various diseases. The most common method for increasing Se uptake is the consumption of Se-rich plants, which transform inorganic Se into organic forms. Wheat is eaten daily by many people. The Se content of Aegilops tauschii (Ae. tauschii), one of the ancestors of hexaploid common wheat, is generally higher than that of wheat. In this study, two genotypes of Ae. tauschii with contrasting Se-accumulating abilities were subjected to different Se treatments followed by high-throughput transcriptome sequencing.

Results

Sequencing of 12 transcriptome libraries of Ae. tauschii grown under different Se treatments produced about a total of 47.72 GB of clean reads. After filtering out rRNA sequences, approximately 19.3 million high-quality clean reads were mapped to the reference genome (ta IWGSC_MIPSv2.1 genome DA). The total number of reference genome gene is 32,920 and about 26,407 known genes were detected in four groups. Functional annotation of these mapped genes revealed a large number of genes and some pivotal pathways that may participate in Se metabolism. The expressions of several genes potentially involved in Se metabolism were confirmed by quantitative real-time PCR.

Conclusions

Our study, the first to examine Se metabolism in Ae. tauschii, has provided a theoretical foundation for future elucidation of the mechanism of Se metabolism in this species.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0700-1) contains supplementary material, which is available to authorized users.

Gradual and Acute Temperature Rise Induces Crossing Endocrine, Metabolic, and Immunological Pathways in Maraena Whitefish (Coregonus maraena)

The complex and still poorly understood nature of thermoregulation in various fish species complicates the determination of the physiological status on the basis of diagnostic marker genes and indicative molecular pathways. The present study aimed to compare the physiological impacts of both gradual and acute temperature rise from 18 to 24°C on maraena whitefish in aquaculture. Microarray-based transcriptome profiles in the liver, spleen and kidney of heat-stressed maraena whitefish revealed the modulation of a significantly higher number of genes in those groups exposed to gradually rising temperatures compared with the acutely stressed groups, which might reflect early adaptation mechanisms. Moreover, we suggest a common set of 11 differentially expressed genes that indicate thermal stress induced by gradual or acute temperature rise in the three selected tissues. Besides the two pathways regulated in both data sets unfolded protein response and aldosterone signaling in epithelial cells, we identified unique tissue- and stress type-specific pathways reflecting the crossroads between signal transduction, metabolic and immunologic pathways to cope with thermal stress. In addition, comparing lists of differentially regulated genes with meta-analyzed published data sets revealed that “acute temperature rise”-responding genes that encode members of the HSP70, HSP90, and HSP40 families; their functional homologs; co-chaperones and stress-signal transducers are well-conserved across different species, tissues and/or cell types and experimental approaches.

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.

Transcriptome profiling analysis of grouper during nervous necrosis virus persistent infection

Nervous necrosis virus (NNV) infection has been considered a serious disease in farmed grouper. Particularly, the persistent infection model conducts the grouper into a carrier state that continues to spread the virus through spawning. This particular model makes disease control more difficult in the aquaculture industry. In the present study, we used RNA-Seq, a high-throughput method based on next-generation sequencing, to profile the expression of genes during the period of NNV persistent infection. We evaluated the transcriptomic changes in the brain tissue of grouper. The inactivated-NNV vaccine was used as a comparison group. Based on the differentially expressed genes, highly immune cell active signaling and surface receptor expression were triggered during persistent infection. The interferon-induced response was also highly expressed in the infected brain tissue. However, critical negative regulatory factors of T-cells, such as PD-L1 and LAG3, were up-regulated. The present transcriptome study revealed a comprehensive view of the state of NNV persistent infection and provided insights into the state of impaired NNV clearance in the grouper.

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

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

Cloning and expression of Malabar grouper (Epinephelus malabaricus) ADAR1 gene in response to immune stimulants and nervous necrosis virus

ADARs are RNA editing catalysts that bind double-stranded RNA and convert adenosine to inosine, a process that can lead to destabilization of dsRNA structures and suppression of mRNA translation. In mammals, ADAR1 genes are involved in various cellular pathways, including interferon (IFN)-mediated response. However, the function of fish ADAR1 remains unclear. We report here the cloning of ADAR1 in Malabar grouper (Epinephelus malabaricus) (MgADAR1) and its response to various immune stimulants. The MgADAR1 cDNA is 5371-bp long, consisting of an open reading frame encoding a putative protein of 1381 amino acids, a 235-nt 5'-terminal untranslated region (UTR), and a 990-nt 3'-UTR. The deduced amino acid sequence exhibits signature features of a chitin synthesis regulation domain, two Z-DNA-binding domains (Z alpha), three dsRNA binding motifs (DSRM) and one tRNA-specific and dsRNA adenosine deaminase domain (ADEAMc). MgADAR1 mRNA expressed ubiquitously in tissues of healthy Malabar grouper, with elevated levels in the brain, gills and eyes. In response to poly (I: C), the MgADAR1 mRNA level was significantly up-regulated in the brain and spleen, but not head kidney. Upon nervous necrosis virus (NNV) infection the level of MgADAR1 increased in the brain, whereas Mx increased in the brain, spleen and head kidney. Induction of MgADAR1 by poly (I: C) and NNV was also observed in vitro. Additionally, the expression of MgADAR1 was upregulated by recombinant grouper IFN in grouper cells. These data indicate an intricate interplay between ADAR1 and NNV infection in grouper as MgADAR1 might be regulated in a tissue-specific manner.

Grouper STAT1a is involved in antiviral immune response against iridovirus and nodavirus infection

Signal Transducer and Activator of Transcription 1 (STAT1) has been demonstrated to function as a critical mediator in multiple cell processes, such as cell proliferation, cell death, and innate immune response. Interestingly, two orthologues of human STAT1, including STAT1a and STAT1b genes have been identified in different fish. However, the detailed roles of fish STAT1a in virus replication still remained largely uncertain. Here, we cloned a STAT1a from orange-spotted grouper Epinephelus coioides (EcSTAT1a) and characterized its roles during fish virus infection. EcSTAT1a encoded a 751-aa peptide which shared 97% and 93% identity to STAT1 from mandarin fish (Siniperca chuatsi) and Malabar grouper (Epinephelus malabaricus), respectively. Amino acid alignment analysis showed that EcSTAT1a contained a STAT-int domain, a STAT-alpha domain, a STAT-bind domain (DNA binding domain), a SH2 domain and a STAT1-TAZ2 bind domain. In examined tissues from healthy grouper, the expression of EcSTAT1a was predominant in intestine, gill and liver. In grouper cells, the relative expression levels of EcSTAT1a was significantly increased during red-spotted grouper nervous necrosis virus (RGNNV) or Singapore grouper iridovirus (SGIV) infection. Under fluorescence microscopy, we found that EcSTAT1a mainly localized in the cytoplasm. The ectopic expression of EcSTAT1a in vitro significantly delayed the cytopathic effect (CPE) progression evoked by RGNNV and SGIV. Further studies showed that the expression levels of viral genes, including SGIV major capsid protein (MCP), VP19, ICP-18, LITAF and RGNNV coat protein (CP), RNA-dependent RNA polymerase (RdRp) were all significantly reduced in EcSTAT1a overexpressing cells compared to the control vector transfected cells, suggested that EcSTAT1a exerted antiviral activity against iridovirus and nodavirus. Furthermore, overexpression of EcSTAT1a significantly increased the expression of interferon related cytokines or effectors and pro-inflammatory factors. Together, our results elucidated that EcSTAT1a might function as a critical antiviral factor by regulating the host interferon immune and inflammation response.

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.

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

Asian seabass is an important food fish in Southeast Asia. Viral nervous necrosis (VNN) disease, triggered by nervous necrosis virus (NNV) infection, has caused mass mortality of Asian seabass larvae, resulting in enormous economic losses in the Asian seabass industry. In order to better understand the complex molecular interaction between Asian seabass and NNV, we investigated the transcriptome profiles of Asian seabass epithelial cells, which play an essential role in immune regulation, after NNV infection. Using the next generation sequencing (NGS) technology, we sequenced mRNA from eight samples (6, 12, 24, 48 h post-inoculation) of mock and NNV-infected Asian seabass epithelial cell line, respectively. Clean reads were de novo assembled into a transcriptome consisting of 89026 transcripts with a N50 of 2617 bp. Furthermore, 251 differentially expressed genes (DEGs) in response to NNV infection were identified. Top DEGs include protein asteroid homolog 1-like (ASTE1), receptor-transporting protein 3 (RTP3), heat shock proteins 30 (HSP30) and 70 (HSP70), Viperin, interferon regulatory factor 3 (IRF3) and other genes related to innate immunity. Our data suggest that abundant and diverse genes corresponding to NNV infection. The results of this study could also offer vital information not only for identification of novel genes involved in Asian seabass-NNV interaction, but also for our understanding of the molecular mechanism of Asian seabass' response to viral infection. In addition, 24807 simple sequence repeats (SSRs) were detected in the assembled transcriptome, providing valuable resources for studying genetic variations and accelerating quantitative trait loci (QTL) mapping for disease resistance in Asian seabass in the future.

In Silico Comparative Transcriptome Analysis of Two Color Morphs of the Common Coral Trout (Plectropomus Leopardus)

The common coral trout is one species of major importance in commercial fisheries and aquaculture. Recently, two different color morphs of Plectropomus leopardus were discovered and the biological importance of the color difference is unknown. Since coral trout species are poorly characterized at the molecular level, we undertook the transcriptomic characterization of the two color morphs, one black and one red coral trout, using Illumina next generation sequencing technologies. The study produced 55162966 and 54588952 paired-end reads, for black and red trout, respectively. De novo transcriptome assembly generated 95367 and 99424 unique sequences in black and red trout, respectively, with 88813 sequences shared between them. Approximately 50% of both trancriptomes were functionally annotated by BLAST searches against protein databases. The two trancriptomes were enriched into 25 functional categories and showed similar profiles of Gene Ontology category compositions. 34110 unigenes were grouped into 259 KEGG pathways. Moreover, we identified 14649 simple sequence repeats (SSRs) and designed primers for potential application. We also discovered 130524 putative single nucleotide polymorphisms (SNPs) in the two transcriptomes, supplying potential genomic resources for the coral trout species. In addition, we identified 936 fast-evolving genes and 165 candidate genes under positive selection between the two color morphs. Finally, 38 candidate genes underlying the mechanism of color and pigmentation were also isolated. This study presents the first transcriptome resources for the common coral trout and provides basic information for the development of genomic tools for the identification, conservation, and understanding of the speciation and local adaptation of coral reef fish species.

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

Infectious salmon anaemia virus (ISAV) is an orthomyxovirus causing high mortality in farmed Atlantic salmon (Salmo salar). The collective data from the Atlantic salmon-ISAV interactions, performed "in vitro" using various salmon cell lines and "in vivo" fish infected with different ISAV isolates, have shown a strong regulation of immune related transcripts during the infection. Despite this strong defence response, the majority of fish succumb to infections with ISAV. The deficient protection of the host against ISAV is in part due to virulence factors of the virus, which allow evade the host-defence machinery. As such, the viral replication is uninhibited and viral loads quickly spread to several tissues causing massive cellular damage before the host can develop an effective cell-mediated and humoral outcome. To interrogate the correlation of the viral replication with the host defence response, we used fish that have been infected by cohabitation with ISAV-injected salmons. Whole gene expression patterns were measured with RNA-seq using RNA extracted from Head-kidney, Liver and Gills. The results show divergent mRNA abundance of functional modules related to interferon pathway, adaptive/innate immune response and cellular proliferation/differentiation. Furthermore, gene regulation in distinct tissues during the infection process was independently controlled within the each tissue and the observed mRNA expression suggests high modulation of the ISAV-segment transcription. Importantly this is the first time that strong correlations between functional modules containing significant immune process with protein-protein affinities and viral-segment transcription have been made between different tissues of ISAV-infected fish.

De novo characterization of fall dormant and nondormant alfalfa (Medicago sativa L.) leaf transcriptome and identification of candidate genes related to fall dormancy

Alfalfa (Medicago sativa L.) is one of the most widely cultivated perennial forage legumes worldwide. Fall dormancy is an adaptive character related to the biomass production and winter survival in alfalfa. The physiological, biochemical and molecular mechanisms causing fall dormancy and the related genes have not been well studied. In this study, we sequenced two standard varieties of alfalfa (dormant and non-dormant) at two time points and generated approximately 160 million high quality paired-end sequence reads using sequencing by synthesis (SBS) technology. The de novo transcriptome assembly generated a set of 192,875 transcripts with an average length of 856 bp representing about 165.1 Mb of the alfalfa leaf transcriptome. After assembly, 111,062 (57.6%) transcripts were annotated against the NCBI non-redundant database. A total of 30,165 (15.6%) transcripts were mapped to 323 Kyoto Encyclopedia of Genes and Genomes pathways. We also identified 41,973 simple sequence repeats, which can be used to generate markers for alfalfa, and 1,541 transcription factors were identified across 1,350 transcripts. Gene expression between dormant and non-dormant alfalfa at different time points were performed, and we identified several differentially expressed genes potentially related to fall dormancy. The Gene Ontology and pathways information were also identified. We sequenced and assembled the leaf transcriptome of alfalfa related to fall dormancy, and also identified some genes of interest involved in the fall dormancy mechanism. Thus, our research focused on studying fall dormancy in alfalfa through transcriptome sequencing. The sequencing and gene expression data generated in this study may be used further to elucidate the complete mechanisms governing fall dormancy in alfalfa.

Characterization of the transcriptome of the Asian gypsy moth Lymantria dispar identifies numerous transcripts associated with insecticide resistance

Although the Asian gypsy moth Lymantria dispar causes extensive forest damage worldwide, little is known regarding the genes involved in its development or response to insecticides. Accordingly, characterization of the transcriptome of L. dispar larvae would promote the development of toxicological methods for its control. RNA-seq analysis of L. dispar larvae messenger RNA (mRNA) generated 62,063 unigenes with N50 of 993 bp, from which 23,975 unique sequences (E-value < 10(-5)) were identified using a BLASTx search of the NCBI non-redundant (nr) database. Using functional classification in the Gene Ontology (GO) and Clusters of Orthologous Groups (COG) databases, 7,309 indentified sequences were categorized into 51 functional groups and 8,079 sequences were categorized into 25 functional groups, respectively. Moreover, we identified a large number of transcripts encoding known insecticide targets, or proteins involved in the metabolism of insecticides. Reads per kilobase of unigene length per million mapped reads (RPKM) analysis identified 39 high abundance transcripts, of which 27 exhibited significantly altered expression patterns across the egg, larvae, pupae, male and female adult stages. Our study provides the most comprehensive transcriptomic sequence resource for L. dispar, which will form the basis for future identification of candidate insecticide resistance genes in L. dispar.

Expression analysis of heat shock protein genes during Aeromonas hydrophila infection in rohu, Labeo rohita, with special reference to molecular characterization of Grp78

Heat shock protein (Hsp) genes are stress-related genes that activate the host immune system during infection. Hsp genes expression in fish, studied during bacterial infections, is mostly confined to Hsp70 and Hsp90. The present study is an expression analysis of seven Hsp genes: Apg2, Hsp90, Hsp70, glucose-regulated protein 78 (Grp78), heat shock cognate 70 (Hsc70), Grp75, and Hsp30 during Aeromonas hydrophila infection in rohu, Labeo rohita. Forty-eight rohu juveniles were challenged with 3 × 10(7) cfu bacteria per 20 g of body weight intraperitoneally. The expression of these genes was studied in infected liver, anterior kidney, and spleen tissues of rohu at different time periods: 0, 1, 3, 6, 12, 24, 48, 72 h, 7, and 15 days post-infection by qPCR. The Hsp gene modulation was greater in liver as compared to spleen and kidney tissues. Induced expression of Apg2, Hsp90, Grp78, Grp75, and Hsc70 was noticed during peak periods (3 to 24 h post-challenge) of bacterial infectivity. Hsp70 was found to be down-regulated during the process of infection. In contrast to the other six genes, Hsp30 showed a variable expression pattern in all three tissues. Grp78 was found to be the most highly (1,587-fold) expressed gene in liver at 12 h post-challenge. Further, molecular characterization of Grp78 revealed it to be a highly conserved Hsp gene, essential not only during infection but also during early developmental stages of rohu, and its expression was noticed in all organs studied except in gill tissues, which indicated its potential immune regulatory role during infection process.

Transcriptome analysis of medaka following epinecidin-1 and TH1-5 treatment of NNV infection

Nervous necrosis virus (NNV) infects a wide range of larval and juvenile fish species, thereby causing enormous economic losses in the aquaculture industry. Possible solutions to this problem include the use of antimicrobial peptides (AMPs), which directly inhibit bacterial growth, and also modulate host signaling mechanisms. The AMPs epinecidin (Epi)-1 and Tilapia hepcidin (TH) 1-5 have been demonstrated to be effective against Nervous necrosis virus infection in medaka (Oryzias latipes). However, the underlying molecular mechanisms are yet to be explored. Here, microarray analyses were performed to examine how NNV infection and/or epinecidin-1 or TH1-5 treatment affects gene expression in medaka; such analyses enabled the prediction of host signaling pathways affected by virus infection and/or regulated by epinecidin-1 and TH1-5. Transcriptome analysis revealed altered expression of genes involved in B cell activation, T cell activation, adipocytokine signaling, and mast cell activation. We subsequently used real-time PCR to analyze expression of key genes involved in these signaling mechanisms. Medaka infected with NNV exhibited up-regulation of PVALB, CEBPA, IFIM, IFN, IL-6ST, NF-kB2, SOC3, SP1, and TGFB1, and such increases were prevented by pre-treatment with epinecidin-1 or TH1-5. Immunohistochemistry using the anti-NNV antibody to stain brain and eye sections revealed that epinecidin-1 treatment during or after infection clears viral load, while TH1-5 treatment only reduces viral numbers if applied during infection. These observations demonstrate that epinecidin-1 and TH1-5 modulate NNV-induced host signaling mechanisms, thereby preventing viral multiplication in host organisms.

Immunity to betanodavirus infections of marine fish

Betanodaviruses cause viral nervous necrosis in numerous fish species, but some species are resistant to infection by these viruses. It is essential to fully characterize the immune responses that underlie this protective response. Complete characterization of the immune responses against nodaviruses may allow the development of methods that stimulate fish immunity and of an effective betanodavirus vaccine. Such strategies could include stimulation of specific immune system responses or blockage of factors that decrease the immune response. The innate immune system clearly provides a front-line defense, and this includes the production of interferons and other cytokines. Interferons that are released inside infected cells and that suppress viral replication may be the most ancient form of innate immunity. This review focuses on the immune responses of fish to betanodavirus infection.

Contrasted innate responses to two viruses in zebrafish: insights into the ancestral repertoire of vertebrate IFN-stimulated genes

Ease of imaging and abundance of genetic tools make the zebrafish an attractive model host to understand host-pathogen interactions. However, basic knowledge regarding the identity of genes involved in antiviral immune responses is still lagging in this species. We conducted a microarray analysis of the larval zebrafish response to two models of RNA virus infections with very different outcomes. Chikungunya virus (CHIKV) induces a rapid and protective IFN response. Infection with infectious hematopoietic necrosis virus is lethal and is associated with a delayed and inefficient IFN response. A typical signature of IFN-stimulated genes (ISGs) was observed with both viruses, but was stronger for CHIKV. We further compared the zebrafish and human ISG repertoires and made a genomic and phylogenic characterization of the main gene families. We describe a core set of well-induced ISGs conserved across vertebrates, as well as multigenic families diversified independently in each taxon. The conservation of ISGs involved in antiviral signaling indicates conservation of the main feedback loops in these pathways. Whole-mount in situ hybridization of selected transcripts in infected larvae revealed a typical pattern of expression for ISGs in the liver, gut, and blood vessels with both viruses. We further show that some inflammatory genes were additionally induced through IFN-independent pathways by infectious hematopoietic necrosis virus and not by CHIKV. This study provides a useful reference set for the analysis of host-virus interactions in zebrafish and highlights the differences between protective and nonprotective antiviral innate responses.

Identification of the STAT1 gene and the characterisation of its immune response to immunostimulants, including nervous necrosis virus (NNV) infection, in Malabar grouper (Epinephelus malabaricus)

Signal Transducer and Activator of Transcription (STAT)-1 is an indispensable signal transduction protein that is involved in the interferon pathway. STAT-1 plays an important role in the innate immune response. The full-length cDNA of Malabar grouper (Epinephelus malabaricus) STAT-1, MgSTAT1, was cloned. Phylogenetic analysis was performed based on the amino acid sequence. Our results indicate that STAT1 is highly conserved with other vertebrates. We also report the expression of MgSTAT1 in different tissues treated with immune stimulants, including LPS, CpG ODN, and poly (I:C), in vivo. The expression of MgSTAT1 was significantly induced in the head kidney upon treatment with poly (I:C) compared to the control. Moreover, the results indicate that MgSTAT1 is up-regulated during nervous necrosis virus (NNV) infection. This study reveals that similar to the mammalian antiviral response, MgSTAT1 mediates the immune response in Malabar grouper.

Gene expression changes leading extreme alkaline tolerance in Amur ide (Leuciscus waleckii) inhabiting soda lake

Background

Amur ide (Leuciscus waleckii) is an economically and ecologically important cyprinid species in Northern Asia. The Dali Nor population living in the soda lake Dali Nor can adapt the extremely high alkalinity, providing us a valuable material to understand the adaptation mechanism against extreme environmental stress in teleost.

Results

In this study, we generated high-throughput RNA-Seq data from three tissues gill, liver and kidney of L. waleckii living in the soda lake Dali Nor and the fresh water lake Ganggeng Nor, then performed parallel comparisons of three tissues. Our results showed that out of assembled 64,603 transcript contigs, 28,391 contigs had been assigned with a known function, corresponding to 20,371 unique protein accessions. We found 477, 2,761 and 3,376 differentially expressed genes (DEGs) in the gill, kidney, and liver, respectively, of Dali Nor population compared to Ganggeng Nor population with FDR ≤ 0.01and fold-change ≥ 2. Further analysis revealed that well-known functional categories of genes and signaling pathway, which are associated with stress response and extreme environment adaptation, have been significantly enriched, including the functional categories of “response to stimulus”, “transferase activity”, “transporter activity” and “oxidoreductase activity”, and signaling pathways of “mTOR signaling”, “EIF2 signaling”, “superpathway of cholesterol biosynthesis”. We also identified significantly DEGs encoding important modulators on stress adaptation and tolerance, including carbonic anhydrases, heat shock proteins, superoxide dismutase, glutathione S-transferases, aminopeptidase N, and aminotransferases.

Conclusions

Overall, this study demonstrated that transcriptome changes in L. waleckii played a role in adaptation to complicated environmental stress in the highly alkalized Dali Nor lake. The results set a foundation for further analyses on alkaline-responsive candidate genes, which help us understand teleost adaptation under extreme environmental stress and ultimately benefit future breeding for alkaline-tolerant fish strains.