Characterization of grass carp spleen transcriptome during GCRV infection

The role of MASP1 in the complement system and expression characteristics in response to GCRV infection in grass carp

The grass carp (Ctenopharyngodon idella) constitutes a significant economic resource within the aquaculture sector of our nation, yet it has been chronically afflicted by the Grass Carp Reovirus (GCRV) disease. The complement system, a vital component of fish's innate immunity, plays a crucial role in combating viral infections. This research investigates the potential role of MASP1, a key molecule in the lectin pathway of the complement system, in the GCRV infection in grass carp. An analysis of the molecular characteristics of MASP1 in grass carp revealed that its identity and similarity percentages range from 35.10 to 91.00 % and 35.30-91.00 %, respectively, in comparison to other species. Phylogenetically, MASP1 in C. idella aligns closely with species such as Danio rerio, Cyprinus carpio, and Carassius carassius, exhibiting chromosomal collinearity with the zebrafish. Subsequent tissue analysis in both healthy and GCRV-infected grass carp indicated that MASP1's basal expression was predominantly in the liver. Post-GCRV infection, MASP1 expression in various tissues exhibited temporal variations: peaking in the liver on day 5, spleen on day 7, and kidney on day 14. Furthermore, employing Complement Component 3 (C3) as a benchmark for complement system activation, it was observed that MASP1 could activate and cleave C3 to C3b. MASP1 also demonstrated an inhibitory effect on GCRV replication (compared with the control group, VP2 and VP7 decreased by 6.82-fold and 4.37-fold) and enhanced the expression of antiviral genes, namely IRF3, IRF7 and IFN1 (compared with the control group, increased 2.25-fold, 45.38-fold and 22.37-fold, respectively). In vivo protein injection experiments substantiated MASP1's influence on the relative mRNA expression levels of C3 in various tissues and its protein expression in serum. This study also verified that C3 could modulate the expression of antiviral genes such as IFN1 and IRF3.

Preliminary study of BF/C2 on immune mechanism of grass carp against GCRV infection

BF/C2 is a crucial molecule in the coagulation complement cascade pathway and plays a significant role in the immune response of grass carp through the classical, alternative, and lectin pathways during GCRV infection. In vivo experiments demonstrated that the mRNA expression levels of BF/C2 (A, B) in grass carp positively correlated with GCRV viral replication at various stages of infection. Excessive inflammation leading to death coincided with peak levels of BF/C2 (A, B) mRNA expression and GCRV viral replication. Correspondingly, BF/C2 (A, B) recombinant protein, CIK cells and GCRV co-incubation experiments yielded similar findings. Therefore, 3 h (incubation period) and 9 h (death period) were selected as critical points for this study. Transcriptome sequencing analysis revealed significant differences in the expression of BF/C2A and BF/C2B during different stages of CIK infection with GCRV and compared to the blank control group (PBS). Specifically, the BF/C2A_3 and BF/C2A_9 groups exhibited 2729 and 2228 differentially expressed genes (DEGs), respectively, with 1436 upregulated and 1293 downregulated in the former, and 1324 upregulated and 904 downregulated in the latter. The BF/C2B_3 and BF/C2B_9 groups showed 2303 and 1547 DEGs, respectively, with 1368 upregulated and 935 downregulated in the former, and 818 upregulated and 729 downregulated in the latter. KEGG functional enrichment analysis of these DEGs identified shared pathways between BF/C2A and PBS groups at 3 and 9 h, including the C-type lectin receptor signaling pathway, protein processing in the endoplasmic reticulum, Toll-like receptor signaling pathway, Salmonella infection, apoptosis, tight junction, and adipocytokine signaling pathway. Additionally, the BF/C2B groups at 3 and 9 h shared pathways related to protein processing in the endoplasmic reticulum, glycolysis/gluconeogenesis, and biosynthesis of amino acids. The mRNA levels of these DEGs were validated in cellular models, confirming consistency with the sequencing results. In addition, the mRNA expression levels of these candidate genes (mapk1, il1b, rela, nfkbiab, akt3a, hyou1, hsp90b1, dnajc3a et al.) in the head kidney, kidney, liver and spleen of grass carp immune tissue were significantly different from those of the control group by BF/C2 (A, B) protein injection in vivo. These candidate genes play an important role in the response of BF/C2 (A, B) to GCRV infection and it also further confirmed that BF/C2 (A, B) of grass carp plays an important role in coping with GCRV infection.

Comparative Analysis of mRNA, microRNA of Transcriptome, and Proteomics on CIK Cells Responses to GCRV and Aeromonas hydrophila

Grass Carp Reovirus (GCRV) and Aeromonas hydrophila (Ah) are the causative agents of haemorrhagic disease in grass carp. This study aimed to investigate the molecular mechanisms and immune responses at the miRNA, mRNA, and protein levels in grass carp kidney cells (CIK) infected by Grass Carp Reovirus (GCRV, NV) and Aeromonas hydrophilus (Bacteria, NB) to gain insight into their pathogenesis. Within 48 h of infection with Grass Carp Reovirus (GCRV), 99 differentially expressed microRNA (DEMs), 2132 differentially expressed genes (DEGs), and 627 differentially expressed proteins (DEPs) were identified by sequencing; a total of 92 DEMs, 3162 DEGs, and 712 DEPs were identified within 48 h of infection with Aeromonas hydrophila. It is worth noting that most of the DEGs in the NV group were primarily involved in cellular processes, while most of the DEGs in the NB group were associated with metabolic pathways based on KEGG enrichment analysis. This study revealed that the mechanism of a grass carp haemorrhage caused by GCRV infection differs from that caused by the Aeromonas hydrophila infection. An important miRNA-mRNA-protein regulatory network was established based on comprehensive transcriptome and proteome analysis. Furthermore, 14 DEGs and 6 DEMs were randomly selected for the verification of RNA/small RNA-seq data by RT-qPCR. Our study not only contributes to the understanding of the pathogenesis of grass carp CIK cells infected with GCRV and Aeromonas hydrophila, but also serves as a significant reference value for other aquatic animal haemorrhagic diseases.

Comparative Transcriptomics Analysis Reveals Unique Immune Response to Grass Carp Reovirus Infection in Barbel Chub (Squaliobarbus curriculus)

Grass carp (Ctenopharyngodon idella) and barbel chub (Squaliobarbus curriculus)-both Leuciscinae subfamily species-demonstrate differences in grass carp reovirus (GCRV) infection resistance. We infected barbel chubs with type II GCRV and subjected their liver, spleen, head kidney, and trunk kidney samples to investigate anti-GCRV immune mechanisms via RNA sequencing and quantitative real-time polymerase chain reaction (qRT-PCR). We identified 139, 970, 867, and 2374 differentially expressed genes (DEGs) in the liver, spleen, head kidney, and trunk kidney, respectively. Across all four tissues, gene ontology analysis revealed significant immune response-related DEG enrichment, and the Kyoto Encyclopedia of Genes and Genomes analysis revealed pattern recognition receptor (PRR) and cytokine-related pathway enrichment. We noted autophagy pathway enrichment in the spleen, head kidney, and trunk kidney; apoptosis pathway enrichment in the spleen and trunk kidney; and complement- and coagulation-cascade pathway enrichment in only the spleen. Comparative transcriptome analysis between GCRV-infected barbel chubs and uninfected barbel chubs comprehensively revealed that PRR, cytokine-related, complement- and coagulation-cascade, apoptosis, and autophagy pathways are potential key factors influencing barbel chub resistance to GCRV infection. qRT-PCR validation of 11 immune-related DEGs confirmed our RNA-seq data's accuracy. These findings provide a theoretical foundation and empirical evidence for the understanding of GCRV infection resistance in barbel chub and hybrid grass carp-barbel chub breeding.

The Teleost CXCL13-CXCR5 Axis Induces Inflammatory Cytokine Expression through the Akt-NF-κB, p38-AP-1, and p38-NF-κB Pathways

The ancestors of chemokines originate in the most primitive of vertebrates, which has recently attracted great interest in the immune functions and the underlying mechanisms of fish chemokines. In the current study, we identified an evolutionarily conserved chemokine, CiCXCL13, from a teleost fish, grass carp. CiCXCL13 was characterized by a typical SCY (small cytokine CXC) domain and four cysteine residues (C34, C36, C61, C77), with the first two cysteines separated by a random amino acid residue, although it shared 24.2-54.8% identity with the counterparts from other vertebrates. CiCXCL13 was an inducible chemokine, whose expression was significantly upregulated in the immune tissues of grass carps after grass carp reovirus infection. CiCXCL13 could bind to the membrane of grass carp head kidney leukocytes and promote cell migration, NO release, and the expression of >15 inflammatory cytokines, including IL-1β, TNF-α, IL-10 and TGF-β1, thus regulating the inflammatory response. Mechanistically, CiCXCL13 interacted with its evolutionarily conserved receptor CiCXCR5 and activated the Akt-NF-κB and p38-AP-1 pathways, as well as a previously unrevealed p38-NF-κB pathway, to efficiently induce inflammatory cytokine expression, which was distinct from that reported in mammals. Zebrafish CXCL13 induced inflammatory cytokine expression through Akt, p38, NF-κB, and AP-1 as CiCXCL13. Meanwhile, the CiCXCL13-CiCXCR5 axis-mediated inflammatory activity was negatively shaped by grass carp atypical chemokine receptor 2 (CiACKR2). The present study is, to our knowledge, the first to comprehensively define the immune function of CXCL13 in inflammatory regulation and the underlying mechanism in teleosts, and it provides a valuable perspective on the evolution and biology of fish chemokines.

Integrative transcriptome analysis reveals alternative polyadenylation potentially contributes to GCRV early infection

Introduction

Grass carp reovirus (GCRV), a member of the Aquareovirus genus in the Reoviridae family, is considered to be the most pathogenic aquareovirus. Productive viral infection requires extensive interactions between viruses and host cells. However, the molecular mechanisms underlying GCRV early infection remains elusive.

Methods

In this study we performed transcriptome and DNA methylome analyses with Ctenopharyngodon idellus kidney (CIK) cells infected with GCRV at 0, 4, and 8 h post infection (hpi), respectively.

Results

We found that at early infection stage the differentially expressed genes related to defense response and immune response in CIK cells are activated. Although DNA methylation pattern of CIK cells 8 hpi is similar to mock-infected cells, we identified a considerable number of genes that selectively utilize alternative polyadenylation sites. Particularly, we found that biological processes of cytoskeleton organization and regulation of microtubule polymerization are statistically enriched in the genes with altered 3'UTRs.

Discussion

Our results suggest that alternative polyadenylation potentially contributes to GCRV early infection.

Structural comparison and expression function analysis of BF/C2 in Ctenopharyngodon idella and Squaliobarbus curriculus

Ctenopharyngodon idella and Squaliobarbus curriculus, members of the Cyprinidae family and Yaroideae subfamily, have shown different levels of resistance to grass carp reo virus (GCRV), with S. curriculus exhibiting higher resilience. In the pursuit to explore the distinctions in the structural and expression traits of BF/C2 (A,B) between the two species, we conducted an analysis involving the cloning and examination of various coding sequences (CDS). We successfully cloned the CDS of ci-BF/C2A and ci-BF/C2B from C. idella, which spanned 2259 bp and 2514 bp respectively, encoding 752 and 837 amino acids. Similarly, the CDS of sc-BF/C2A and sc-BF/C2B from S. curriculus were cloned, featuring lengths of 1353 bp and 2517 bp and encoding 450 and 838 amino acids, respectively. A chromosome collinearity assessment revealed that ci-BF/C2A demonstrated collinearity with sc-BF/C2A, a finding not replicated with ci-BF/C2B and sc-BF/C2B. Delving into gene structure, we discerned that ci-BF/C2A harbored a greater number of Tryp_SPc domains compared to sc-BF/C2A. Following this, we engineered and purified six prokaryotic recombinant proteins: CI-BF/C2A, CI-BF/C2A1 (a variant resulting from the deletion of the Tryp_SPc domain of CI-BF/C2A), CI-BF/C2A2 (representing the Tryp_SPc domain of CI-BF/C2A), CI-BF/C2B, SC-BF/C2A, and SC-BF/C2B. Through serum co-incubation tests with these recombinant proteins, we established the activation of the complement marker C3 in each case. Utilizing fluorescence quantitative expression analysis, we observed ubiquitous expression of ci-BF/C2A and ci-BF/C2B across all grass carp tissues, predominantly in the liver. This pattern mirrored in S. curriculus, where sc-BF/C2A was highly expressed in the gills, and sc-BF/C2B manifested notably in the liver. Kidney cell infection experiments on both species revealed enhanced resistance to GCRV post-incubation with the recombinant proteins. Notably, cells treated with SC-BF/C2 (A, B) exhibited pronounced resilience compared to those treated with CI-BF/C2 (A, B, A1, A2). However, cells incubated with CI-BF/C2A1 and CI-BF/C2A2 showed strengthen resistance relative to cells treated with CI-BF/C2A and CI-BF/C2B. In GCRV infection trials on grass carp, ci-BF/C2A and ci-BF/C2B expressions reached a zenith on the seventh day post-infection, highlighting a distinctive functional mode in immune defense against GCRV infection orchestrated by BF/C2. The empirical data underscores the pivotal role of the Tryp_SPc domain in immune responses to GCRV infection, pinpointing its influence on ci-BF/C2A expression. Conclusively, this investigation provides a foundational understanding of the unique immune function characteristics of BF/C2 in grass carp, paving the way for further scholarly exploration in this realm.

Functional and expression analysis reveals the involvement of integrin αvβ3 in antiviral immunity of grass carp (Ctenopharyngodon idella)

Integrins are α-β heterodimeric cell receptors that can bind the protein components of pathogens, and play crucial roles in mammalian immune responses, but the immune functions mediated by integrins remains largely unknown in teleost fish. In this study, an integrin αvβ3 (GCαvβ3) originally assembled by αv (GCαv) and β3 (GCβ3) subunits, was identified from a teleost fish grass carp Ctenopharyngodon idella. The pairwise alignment analyses showed that the amino acid sequences of GCαv and GCβ3 shared high similarity (75.2-95.1%) and identity (58.6-90.7%) with their homologs from other vertebrates. Both GCαv and GCβ3 harbored the conserved protein domains and motifs, and were clustered in fish branch of the phylogenetic tree containing the counterparts from various vertebrates. Co-immunoprecipitation displayed that GCβ3 could interact with the grass carp reovirus (GCRV) outer capsid protein VP5. Two incubation experiments revealed that the interaction of GCRV or VP5 proteins with GCβ3 could induce the expressions of type I interferons (IFNs) including IFN2 and IFN3 in grass carp ovary cell line. The functional analysis demonstrated that GCαvβ3 served as a receptor of viral protein components to be involved in antiviral immunity as human integrin αvβ3 did. In addition, both GCαv and GCβ3 were significantly upregulated in various tissues of grass carp after GCRV infection. This study might provide fundamental basis for understanding the molecular characteristics and immune functions of GCαvβ3, and offer a new insight into the antiviral immune mechanism specific to the integrins in grass carp.

Functional and Expressional Analyses Reveal the Distinct Role of Complement Factor I in Regulating Complement System Activation during GCRV Infection in Ctenopharyngodon idella

Complement factor I (CFI), a complement inhibitor, is well known for regulating the complement system activation by degrading complement component 3b (C3b) in animal serum, thus becoming involved in innate defense. Nevertheless, the functional mechanisms of CFI in the complement system and in host-pathogen interactions are far from being clarified in teleost fish. In the present study, we cloned and characterized the CFI gene, CiCFI, from grass carp (Ctenopharyngodon idella) and analyzed its function in degrading serum C3b and expression changes after grass carp reovirus (GCRV) infection. The open reading frame of CiCFI was found to be 2121 bp, encoding 706 amino acids with a molecular mass of 79.06 kDa. The pairwise alignments showed that CiCFI shared the highest identity (66.9%) with CFI from Carassius gibelio and the highest similarity (78.7%) with CFI from Danio rerio. The CiCFI protein was characterized by a conserved functional core Tryp_SPc domain with the catalytic triad and substrate binding sites. Phylogenetic analysis indicated that CiCFI and the homologs CFIs from other teleost fish formed a distinct evolutionary branch. Similar with the CFIs reported in mammals, the recombinant CiCFI protein could significantly reduce the C3b content in the serum, demonstrating the conserved function of CiCFI in the complement system in the grass carp. CiCFI mRNA and protein showed the highest expression level in the liver. After GCRV infection, the mRNA expressions of CiCFI were first down-regulated, then up-regulated, and then down-regulated to the initial level, while the protein expression levels maintained an overall downward trend to the late stage of infection in the liver of grass carps. Unexpectedly, the protein levels of CiCFI were also continuously down-regulated in the serum of grass carps during GCRV infection, while the content of serum C3b proteins first increases and then returns to the initial level, suggesting a distinct role of CiCFI in regulating complement activation and fish-virus interaction. Combining our previous results that complement factor D, a complement enhancer, shows continuously up-regulated expression levels in grass carps during GCRV infection, and this study may provide the further essential data for the full picture of complex complement regulation mechanism mediated by Df and CFI of the grass carp during pathogen infection.

Cloning of Grass Carp Chemokine XC Receptor 1 (XCR1) Gene and Evaluation of Its Expression in Various Organs after GCRV Infection

Interaction between the chemokine receptor XCR1 and its ligand is closely related to the immune function in animals; however, there are only a few reports on role of XCR1 in the immune system of fish. We aimed to analyze the expression of XCR1 in various organs or tissues of grass carp before and after Grass Carp Reovirus (GCRV) infection to better understand the function of XCR1 in resistance to GCRV infection. We cloned and sequenced the cDNA of grass carp XCR1 and analyzed the molecular structure of XCR1 based its amino acid sequence. Further, we analyzed the relative expression levels of XCR1 in different organs or tissues of male parent grass carp with GCRV resistance (P1) and their first-generation offspring (F1) before and after GCRV infection. Our results show that the total length of cDNA of the grass carp XCR1 gene is 1659 bp and encodes 365 amino acids. XCR1 contains seven conserved transmembrane helical domains. The homologous tertiary structure of XCR1 is similar to its homologs in other species. After artificial GCRV infection, there were significant differences in the expression of the grass carp XCR1 gene in different tissues, at different time points, and between P1 and F1 fish. These results will contribute to our understanding of the role of XCR1 in fish immune responses and contribute to the development of GCRV-resistant grass carp.

Functional Identification of Complement Factor D and Analysis of Its Expression during GCRV Infection in Grass Carp (Ctenopharyngodon idella)

Complement factor D (Df) is a serine protease well known for activating the alternative pathway (AP) in mammals by promoting the cleavage of complement component 3 (C3), thus becoming involved in innate defense. In teleost fish, however, the functional mechanisms of Df in the AP and against pathogen infection are far from clear. In the present study, we cloned and characterized the Df gene, CiDf, from grass carp (Ctenopharyngodon idella) and analyzed its function in promoting C3 cleavage and expression changes after grass carp reovirus (GCRV) infection. The open reading frame of CiDf was found to be 753 bp, encoding 250 amino acids with a molecular mass of 27.06 kDa. CiDf harbors a conserved Tryp_SPc domain, with three conserved residues representing the catalytic triad and three conserved binding sites in the substrate specificity pocket. Pairwise alignment showed that CiDf shares the highest identity (96%) and similarity (98%) with Df from Anabarilius grahami. Phylogenetic analysis indicated that CiDf and other fish Dfs formed a distinct evolutionary branch. Similar to most Dfs from other vertebrates, the CiDf gene structure is characterized by four introns and five exons. The incubation of recombinant CiDf protein with grass carp serum significantly increased the C3b content, demonstrating the conserved function of CiDf in the AP in promoting C3 cleavage, similar to Dfs in mammals. CiDf mRNA expression was widely detected in various tissues and levels were relatively higher in the liver, spleen, and intestine of grass carp. During GCRV infection over a 168-hour period, a high level of CiDf mRNA expression in the liver, spleen, and intestine was maintained at 144 and 168 h, suggesting AP activity at the late stage of GCRV infection. Collectively, the above results reveal the conserved structure and function of CiDf and its distinct expression patterns after GCRV infection, which provide a key basis for studying the roles of Df and AP during GCRV infection in the grass carp C. idella.

Host-microbiota interactions and responses to grass carp reovirus infection in Ctenopharyngodon idellus

Gut microbiota could facilitate host to defense diseases, but fish-microbiota interactions during viral infection and the underlying mechanism are poorly understood. We examined interactions and responses of gut microbiota to grass carp reovirus (GCRV) infection in Ctenopharyngodon idellus, which is the most important aquaculture fish worldwide. We found that GCRV infection group with serious haemorrhagic symptoms (G7s) showed considerably different gut microbiota, especially with an abnormally high abundance of gram-negative anaerobic Cetobacterium somerae. It also showed the lowest (p < 0.05) alpha-diversity but with much higher ecological process of homogenizing dispersal (28.8%), confirming a dysbiosis of the gut microbiota after viral infection. Interestingly, signaling pathways of NOD-like receptors (NLRs), toll-like receptors (TLRs), and lipopolysaccharide (LPS) stimulation genes were significantly (q-value < 0.01) enriched in G7s, which also significantly (p < 0.01) correlated with the core gut microbial genera of Cetobacterium and Acinetobacter. The results suggested that an expansion of C. somerae initiated by GCRV could aggravate host inflammatory reactions through the LPS-related NLRs and TLRs pathways. This study advances our understanding of the interplay between fish immunity and gut microbiota challenged by viruses; it also sheds new insights for ecological defense of fish diseases with the help of gut microbiota.

Comparison of the blood parameters and histopathology between grass carp infected with a virulent and avirulent isolates of genotype II grass carp reovirus

Grass carp hemorrhagic disease caused by grass carp reovirus (GCRV) is the most important disease for grass carp aquaculture. Its typical clinical symptom is haemorrhaging, although the mechanism was remained unclear. In this study, we investigated the differences in blood parameters and histopathological features between grass carp infected with a virulent and avirulent isolates of genotype II GCRV. Infection with the virulent isolate resulted in increases in 8 routine blood and 2 serum biochemical parameters (P < 0.05); while 9 routine blood and 5 biochemical parameters were significantly decreased (P < 0.05) compared with fish infected with the avirulent isolate. The majority of these alterations were related to hemorrhage, inflammatory reactions and organic damage. The histopathologic changes were primarily vasodilation and hyperaemia in multiple organs, lymphocyte and macrophage infiltration as well as severe vacuolar degeneration in spleen, kidney and liver. The histopathology changes in fish infected with the avirulent isolate were minimal. These results indicated that the pathogenicity of GCRV was primarily reflected in destruction of the blood circulatory system and parenchymatous organs. This study lays the foundation for further research on the pathogenesis of bleeding caused by GCRV infection and the use of blood parameters and histopathology as tools for disease diagnosis.

Different responses in one-year-old and three-year-old grass carp reveal the mechanism of age restriction of GCRV infection

Grass carp is an important fish species in Chinese aquaculture, and can be afflicted by a hemorrhagic disease caused by the grass carp reovirus (GCRV). Interestingly, the affects of GCRV infection of grass carp are age-restricted, meaning that one-year-old grass carp can be infected and can suffer hemorrhagic disease, but three-year-old carp are not so afflicted. In this study, we investigated the mechanism responsible for this age-restricted pathology. We evaluated the relative copy number of GCRV RNA, the expression levels of proteins in blood, and changes in DNA methylation in carp from the two age groups after infection with GCRV. After GCRV infection, the relative copy number of GCRV RNA in three-year-old grass carp was significantly lower than in one-year-old carp. The differences in circulating protein levels mainly occurred in concentrated in complement and coagulation proteins, and the expression levels of these proteins were significantly higher in three-year-old grass carp than in one-year-old carp. Moreover, the expression levels of DNA methylation-related genes in the liver and spleen of one-year-old grass carp were significantly higher than those of three-year-old carp. These results suggested that as age of grass carp increases, faster and more efficient response of the immune system after viral infection, especially the complement system, and differences in DNA methylation may be important factors that affect the age restriction observed in GCRV infection. Our study provides new insights into the mechanisms underlying age restriction of GCRV infection.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Transcriptomics Sequencing Provides Insights into Understanding the Mechanism of Grass Carp Reovirus Infection

Grass carp is an important aquaculture fish species in China that is affected by severe diseases, especially haemorrhagic disease caused by grass carp reovirus (GCRV). However, the mechanisms of GCRV invasion and infection remain to be elucidated. In the present study, Ctenopharyngodon idellus kidney (CIK) cells were infected with GCRV, harvested at 0, 8, 24, and 72 h post infection, respectively, and then subjected to transcriptomics sequencing. Each sample yielded more than 6 Gb of clean data and 40 million clean reads. To better understand GCRV infection, the process was divided into three phases: the early (0-8 h post infection), middle (8-24 h post infection), and late (24-72 h) stages of infection. A total of 76 (35 up-regulated, 41 down-regulated), 553 (463 up-regulated, 90 down-regulated), and 284 (150 up-regulated, 134 down-regulated) differently expressed genes (DEGs) were identified during the early, middle, and late stages of infection, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that DEGs were mainly involved in carbohydrate biosynthesis, transport, and endocytosis in the early stage, phagocytosis and lysosome pathways were mainly enriched in the middle stage, and programmed cell death, apoptosis, and inflammation were largely associated with the late stage. These results suggest GCRV infection is a gradual process involving adsorption on the cell surface, followed by endocytosis into cells, transport by lysosomes, and eventually resulted in cell necrosis and/or apoptosis. Our findings provide insight into the mechanisms of grass carp reovirus infection.

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.

Transcriptome profiling of developing spleen tissue and discovery of immune-related genes in grass carp (Ctenopharyngodon idella)

Grass carp Ctenopharyngodon idella is an important freshwater aquaculture species. However, studies regarding transcriptomic profiling of developing spleen tissue in the grass carp are lacking. Here, the transcriptome sequencing from the spleen tissue of one-year-old (cis1) and three-year-old (cis3) grass carp was performed using Illumina paired-end sequencing technology. The de novo assemblies yielded 48,970 unigenes with average lengths of 1264.51 bp from the two libraries. The assembled unigenes were evaluated and functionally annotated by comparing with sequences in major public databases including Nr, COG, Swiss-Prot, KEGG, Pfam and GO. Comparative analysis of expression levels revealed that a total of 38,254 unigenes were expressed in both the cis1 and cis3 libraries, while 4356 unigenes were expressed only in the cis1 library, and 3312 unigenes were expressed only in the cis3 library. Meanwhile, 1782 unigenes (including 903 down-regulated and 879 up-regulated unigenes) were differentially expressed between the two developmental stages of the grass carp spleen. Based on GO and KEGG enrichment analysis, these differentially expressed genes widely participated in the regulation of immunity and response in the grass carp. Moreover, the main components of six immune-related pathways were identified, including complement and coagulation cascades, Toll-like receptor signaling, B-cell receptor signaling, T-cell receptor signaling, antigen processing and presentation, and chemokine signaling. Finally, two identified transcripts including TLR 8 and complement component C8 were validated for reliability by RT-PCR. Collectively, the results obtained in this study will provide a basis for the study of molecular mechanisms in grass carp spleen development.