Search

Disease outbreaks in aquaculture have recently intensified. In particular, visceral white-nodules disease, caused by Pseudomonas plecoglossicida , has severely hindered the small yellow croaker ( Larimichthys polyactis ) aquaculture industry. However, research on this disease is limited. To address this gap, the present study employed a 100K SNP chip to genotype individuals from an F1 full-sib family, identify single nucleotide polymorphisms (SNPs), and construct a genetic linkage map for this species. A high-density genetic linkage map spanning a total length of 1395.72 cM with an average interval of 0.08 cM distributed across 24 linkage groups was obtained. Employing post-infection survival time as an indicator of disease resistance, 13 disease resistance-related quantitative trait loci (QTLs) were detected, and these regions included 169 genes. Functional enrichment analyses pinpointed 11 candidate disease resistance-related genes. RT-qPCR analysis revealed that the genes of chmp1a and arg1 are significantly differentially expressed in response to P. plecoglossicida infection in spleen and liver tissues, indicating their pivotal functions in disease resistance. In summary, in addition to successfully constructing a high-density genetic linkage map, this study reports the first QTL mapping for visceral white-nodules disease resistance. These results provide insight into the intricate molecular mechanisms underlying disease resistance in the small yellow croaker.

Ras-related associated with diabetes (RRAD) is a member of the Ras GTPase superfamily that plays a role in several cellular functions, such as cell proliferation and differentiation. In particular, the superfamily acts as an NF-κB signaling pathway inhibitor and calcium regulator to participate in the immune response pathway. A recent transcriptome study revealed that rrad was expressed in the spleen of disease-resistant Japanese flounder ( Paralichthys olivaceus ) individuals compared with disease-susceptible individuals, and the results were also verified by qPCR. Thus, the present study aimed to explore how rrad regulates antimicrobial immunity via the NF-κB pathway. First, the coding sequence of P. olivaceus rrad was identified. The sequence was 1092 bp in length, encoding 364 amino acids. Based on phylogenetic and structural relationship analyses, P. olivaceus rrad appeared to be more closely related to teleosts. Next, rrad expression differences between disease-resistant and disease-susceptible individuals in immune-related tissues were evaluated, and the results revealed that rrad was expressed preferentially in the spleen of disease-resistant individuals. In response to Edwardsiella piscicida infection, rrad expression in the spleen changed. In vitro, co-culture was carried out to assess the hypo-methylated levels of the rrad promoter in the disease-resistant spleen, which was consistent with the high mRNA expression. The siRNA-mediated knockdown of rrad performed with the gill cell line of P. olivaceus affected many rrad -network-related genes, i.e., dcp1b , amagt , rus1 , rapgef1 , ralbp1 , plce1 , rasal1 , nckipsd , prkab2 , cytbc-1 , sh3 , and others, as well as some inflammation-related genes, such as bal2 and Il-1β . In addition, flow cytometry analysis showed that rrad overexpression was more likely to induce cell apoptosis, with establishing a link between rrad 's function and its potential roles in regulating the NF-κB pathway. Thus,. the current study provided some clarity in terms of understanding the immune response about rrad gene differences between disease-resistant and disease-susceptible P. olivaceus individuals. This study provides a molecular basis for fish rrad gene functional analysis and may serve as a reference for in-depth of bacterial disease resistance of teleost.

The disease caused by Largemouth bass ranavirus (LMBV) is one of the most severe viral diseases in largemouth bass ( Micropterus salmoides ). It is crucial to evaluate the genetic resistance of largemouth bass to LMBV and develop markers for disease-resistance breeding. In this study, 100 individuals (45 resistant and 55 susceptible) were sequenced and evaluated for resistance to LMBV and a total of 2,579,770 variant sites (SNPs-single-nucleotide polymorphisms (SNPs) and insertions-deletions (InDels)) were identified. A total of 2348 SNPs-InDels and 1018 putative candidate genes associated with LMBV resistance were identified by genome-wide association analyses (GWAS). Furthermore, GO and KEGG analyses revealed that the 10 candidate genes ( MHC II , p38 MAPK , AMPK , SGK1 , FOXO3 , FOXO6 , S1PR1 , IL7R , RBL2 , and GADD45 ) were related to intestinal immune network for IgA production pathway and FoxO signaling pathway. The acquisition of candidate genes related to resistance will help to explore the molecular mechanism of resistance to LMBV in largemouth bass. The potential polymorphic markers identified in this study are important molecular markers for disease resistance breeding in largemouth bass.

Hsp40-Hsp70 typically function in concert as molecular chaperones, and their roles in post-infection immune responses are increasingly recognized. However, in the economically important fish species Scophthalmus maximus (turbot), there is still a lack in the systematic identification, interaction models, and binding site analysis of these proteins. Herein, 62 Hsp40 genes and 16 Hsp70 genes were identified in the turbot at a genome-wide level and were unevenly distributed on 22 chromosomes through chromosomal distribution analysis. Phylogenetic and syntenic analysis provided strong evidence in supporting the orthologies and paralogies of these HSPs. Protein-protein interaction and expression analysis was conducted to predict the expression profile after challenging with Aeromonas salmonicida . dnajb1b and hspa1a were found to have a co-expression trend under infection stresses. Molecular docking was performed using Auto-Dock Tool and PyMOL for this pair of chaperone proteins. It was discovered that in addition to the interaction sites in the J domain, the carboxyl-terminal domain of Hsp40 also plays a crucial role in its interaction with Hsp70. This is important for the mechanistic understanding of the Hsp40-Hsp70 chaperone system, providing a theoretical basis for turbot disease resistance breeding, and effective value for the prevention of certain diseases in turbot.

E3 ubiquitin ligases, key components of the ubiquitin proteasome system, orchestrate protein degradation through ubiquitylation and profoundly impact cellular biology. Small HERC E3 ligases (HERC3-6) have diverse functions in mammals, including roles in spermatogenesis, protein degradation, and immunity. Until now, only mammals' HERC3, HERC5, and HERC6 are known to participate in immune responses, with major involvement in the antiviral response. Interestingly, an exclusive HERC7 has been characterized in fish showing great molecular conservation and antiviral roles. Thus, this study identifies and characterizes the herc7 gene in the European sea bass teleost. The European sea bass herc7 gene and the putative protein show good conservation of the promoter binding sites for interferons and the RCC1 and HECT domains characteristic of HERC proteins, respectively. The phylogenetic analysis shows a unique cluster with the fish-exclusive HERC7 orthologues. During ontogeny, the herc7 gene is expressed from 3 days post-fertilization onwards, being constitutively and widely distributed in adult tissues. In vitro, stimulated leucocytes up-regulate the herc7 gene in response to mitogens and viruses, pointing to a role in the immune response. Furthermore, sea bass herc7 expression is related to the interferon response intensity and viral load in different tissues upon in vivo infection with red-grouper betanodavirus (RGNNV), suggesting the potential involvement of fish HERC7 in ISGylation-based antiviral activity, similarly to mammalian HERC5. This study broadens the understanding of small HERC proteins in fish species and highlights HERC7 as a potential contributor to the immune response in European sea bass, with implications for antiviral defense mechanisms. Future research is needed to unravel the precise actions and functions of HERC7 in teleost fish immunity, providing insights into direct antiviral activity and viral evasion.

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., Competing Interests: The authors declare no conflict of interest.

The innate immune response in Salmo salar , mediated by pattern recognition receptors (PRRs), is crucial for defending against pathogens. This study examined DDX41 protein functions as a cytosolic/nuclear sensor for cyclic dinucleotides, RNA, and DNA from invasive intracellular bacteria. The investigation determined the existence, conservation, and functional expression of the ddx41 gene in S. salar . In silico predictions and experimental validations identified a single ddx41 gene on chromosome 5 in S. salar , showing 83.92% homology with its human counterpart. Transcriptomic analysis in salmon head kidney confirmed gene transcriptional integrity. Proteomic identification through mass spectrometry characterized three unique peptides with 99.99% statistical confidence. Phylogenetic analysis demonstrated significant evolutionary conservation across species. Functional gene expression analysis in SHK-1 cells infected by Piscirickettsia salmonis and Renibacterium salmoninarum indicated significant upregulation of DDX41, correlated with increased proinflammatory cytokine levels and activation of irf3 and interferon signaling pathways. In vivo studies corroborated DDX41 activation in immune responses, particularly when S. salar was challenged with P. salmonis , underscoring its potential in enhancing disease resistance. This is the first study to identify the DDX41 pathway as a key component in S. salar innate immune response to invading pathogens, establishing a basis for future research in salmonid disease resistance.

In aquaculture, viral diseases pose a significant threat and can lead to substantial economic losses. The primary defense against viral invasion is the innate immune system, with interferons (IFNs) playing a crucial role in mediating the immune response. With advancements in molecular biology, the role of non-coding RNA (ncRNA), particularly microRNAs (miRNAs), in gene expression has gained increasing attention. While the function of miRNAs in regulating the host immune response has been extensively studied, research on their immunomodulatory effects in teleost fish, including silver carp ( Hyphthalmichthys molitrix ), is limited. Therefore, this research aimed to investigate the immunomodulatory role of microRNA-30b-5p (miR-30b-5p) in the antiviral immune response of silver carp ( Hypophthalmichthys molitrix ) by targeting cytokine receptor family B5 (CRFB5) via the JAK/STAT signaling pathway. In this study, silver carp were stimulated with polyinosinic-polycytidylic acid (poly (I:C)), resulting in the identification of an up-regulated miRNA (miR-30b-5p). Through a dual luciferase assay, it was demonstrated that CRFB5, a receptor shared by fish type I interferon, is a novel target of miR-30b-5p. Furthermore, it was found that miR-30b-5p can suppress post-transcriptional CRFB5 expression. Importantly, this study revealed for the first time that miR-30b-5p negatively regulates the JAK/STAT signaling pathway, thereby mediating the antiviral immune response in silver carp by targeting CRFB5 and maintaining immune system stability. These findings not only contribute to the understanding of how miRNAs act as negative feedback regulators in teleost fish antiviral immunity but also suggest their potential therapeutic measures to prevent an excessive immune response.

YTH domain-containing genes are important readers of N 6 -methyladenosine (m 6 A) modifications with ability to directly affect the fates of distinct RNAs in organisms. Despite their importance, little is known about YTH domain-containing genes in teleosts until now. In the present study, a total of 10 YTH domain-containing genes have been systematically identified and functionally characterized in rainbow trout ( Oncorhynchus mykiss ). According to the phylogenetic tree, gene structure and syntenic analysis, these YTH domain-containing genes could be classified into three evolutionary subclades, including YTHDF, YTHDC1 and YTHDC2. Of them, the copy number of OmDF1 , OmDF2 , OmDF3 , and OmDC1 were duplicated or even triplicated in rainbow trout due to the salmonid-specific whole-genome duplication event. The three-dimensional protein structure analysis revealed that there were similar structures and the same amino acid residues that were associated with cage formation between humans and rainbow trout, implying their similar manners in binding to m 6 A modification. Additionally, the results of qPCR experiment indicated that the expression patterns of a few YTH domain-containing genes, especially OmDF1b , OmDF3a and OmDF3b , were significantly different in liver tissue of rainbow trout under four different temperatures (7 °C, 11 °C, 15 °C, and 19 °C). The expression levels of OmDF1a , OmDF1b and OmDC1a were obviously repressed in spleen tissue of rainbow trout at 24 h after Yersinia ruckeri infection, while increased expression was detected in OmDF3b . This study provides a systemic overview of YTH domain-containing genes in rainbow trout and reveals their biological roles in responses to temperature stress and bacterial infection.

Vasoactive intestinal peptide (VIP), a member of secretin/glucagon family, is involved in a variety of biological activities such as gut motility, immune responses, and carcinogenesis. In this study, the VIP precursor gene ( On-VIP ) and its receptor gene VIPR1 ( On-VIPR1 ) were identified from Nile tilapia ( Oreochromis niloticus ), and the functions of On-VIP in the immunomodulation of Nile tilapia against bacterial infection were investigated and characterized. On-VIP and On-VIPR1 contain a 450 bp and a 1326 bp open reading frame encoding deduced protein of 149 and 441 amino acids, respectively. Simultaneously, the transcript of both On-VIP and On-VIPR1 were highly expressed in the intestine and sharply induced by Streptococcus agalatiae . Moreover, the positive signals of On-VIP and On-VIPR1 were detected in the longitudinal muscle layer and mucosal epithelium of intestine, respectively. Furthermore, both in vitro and in vivo experiments indicated several immune functions of On-VIP, including reduction of P65 , P38 , MyD88 , STAT3 , and AP1 , upregulation of CREB and CBP , and suppression of inflammation. Additionally, in vivo experiments proved that On-VIP could protect Nile tilapia from bacterial infection and promote apoptosis and pyroptosis. These data lay a theoretical basis for further understanding of the mechanism of VIP guarding bony fish against bacterial infection.

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

A. hydrophila is an important pathogen that mainly harms aquatic animals and has exhibited resistance to a variety of antibiotics. Here, to seek an effective alternative for antibiotics, the effects of umbelliferone (UM) at sub-MICs on A. hydrophila virulence factors and the quorum-sensing system were studied. Subsequently, RNA sequencing was employed to explore the potential mechanisms for the antivirulence activity of umbelliferone. Meanwhile, the protective effect of umbelliferone on grass carp infected with A. hydrophila was studied in vivo. Our results indicated that umbelliferone could significantly inhibit A. hydrophila virulence such as hemolysis, biofilm formation, swimming and swarming motility, and their quorum-sensing signals AHL and AI-2. Transcriptomic analysis showed that umbelliferone downregulated expression levels of genes related to exotoxin, the secretory system (T2SS and T6SS), iron uptake, etc. Animal studies demonstrated that umbelliferone could significantly improve the survival of grass carps infected with A. hydrophila , reduce the bacterial load in the various tissues, and ameliorate cardiac, splenic, and hepatopancreas injury. Collectively, umbelliferone can reduce the pathogenicity of A. hydrophila and is a potential drug for treating A. hydrophila infection.

Yellow catfish ( Pelteobagrus fulvidraco ) is an important aquaculture fish susceptible to Aeromonas veronii infection, which causes acute death resulting in huge economic losses. Understanding the molecular processes of host immune defense is indispensable to disease control. Here, we conducted the integrated and comparative analyses of the transcriptome and metabolome of yellow catfish in response to A. veronii infection at the invaded stage and recovering stage. The crosstalk between A. veronii -induced genes and metabolites uncovered the key biomarkers and pathways that strongest contribute to different response strategies used by yellow catfish at corresponding defense stages. We found that at the A. veronii invading stage, the immune defense was strengthened by synthesizing lipids with energy consumption to repair the skin defense line and accumulate lipid droplets promoting intracellular defense line; triggering an inflammatory response by elevating cytokine IL-6, IL-10 and IL-1β following PAMP-elicited mitochondrial signaling, which was enhanced by ROS produced by impaired mitochondria; and activating apoptosis by up-regulating caspase 3, 7 and 8 and Prostaglandin F1α, meanwhile down-regulating FoxO3 and BCL6. Apoptosis was further potentiated via oxidative stress caused by mitochondrial dysfunction and exceeding inflammatory response. Additionally, cell cycle arrest was observed. At the fish recovering stage, survival strategies including sugar catabolism with D-mannose decreasing; energy generation through the TCA cycle and Oxidative phosphorylation pathways; antioxidant protection by enhancing Glutathione (oxidized), Anserine, and α-ketoglutarate; cell proliferation by inducing Cyclin G2 and CDKN1B; and autophagy initiated by FoxO3, ATG8 and ATP6V1A were highlighted. This study provides a comprehensive picture of yellow catfish coping with A. veronii infection, which adds new insights for deciphering molecular mechanisms underlying fish immunity and developing stage-specific disease control techniques in aquaculture.

The large yellow croaker ( Larimichthys crocea ) is one of the most economically valuable mariculture fish in China. Infection of Pseudomonas plecoglossicida can cause an outbreak of "internal organ white-spot disease", which seriously affects the aquaculture of the large yellow croaker. Ubiquitylation is closely related to the post-translation modification of proteins and plays a vital role in many hosts' immune defense pathways, while the E2-binding enzyme is a key factor in ubiquitination. Our previous genome-wide association study found that the ubiquitin-binding enzyme E2G1 (designed LcUbe2g1 ) was one of the candidate genes related to disease resistance in large yellow croaker. In this study, we analyzed the molecular characteristics, function, and immune mechanism of the LcUbe2g1 . The full-length cDNA is 812 bp, with an open reading frame of 513 bp, encoding 170 amino acid residues. The results of the RT-qPCR and immunohistochemistry analysis revealed that its transcription and translation were significantly activated by the infection of P. plecoglossicida in large yellow croaker. Immunocytochemistry experiments verified the co-localization of Lc UBE2G1 and the ubiquitin proteins in the head kidney cells of large yellow croaker. Through GST pull-down, we found that Lc UBE2G1 interacted with NEDD8 to co-regulate the ubiquitination process. The above results indicate that Lc UBE2G1 is essential in the regulation of ubiquitination against P. plecoglossicida infection in large yellow croaker, which lays a foundation for further study on the resistance mechanism of internal organ white-spot disease.

Phosphoribosyl pyrophosphate synthetases (EC 2.7.6.1) are key enzymes in the biological synthesis of phosphoribosyl pyrophosphate and are involved in diverse developmental processes. In our previous study, the PRPS1 gene was discovered as a key disease-resistance candidate gene in yellow drum, Nibea albiflora , in response to the infection of Vibrio harveyi , through genome-wide association analysis. This study mainly focused on the characteristics and its roles in immune responses of the PRPS1 gene in yellow drum. In the present study, the NaPRPS1 gene was cloned from yellow drum, encoding a protein of 320 amino acids. Bioinformatic analysis showed that Na PRPS1 was highly conserved during evolution. Quantitative RT-PCR demonstrated that NaPRPS1 was highly expressed in the head-kidney and brain, and its transcription and translation were significantly activated by V. harveyi infection examined by RT-qPCR and immunohistochemistry analysis, respectively. Subcellular localization revealed that Na PRPS1 was localized in cytoplasm. In addition, semi-in vivo pull-down assay coupled with mass spectrometry identified myeloid differentiation factor 88 (MyD88) as an Na PRPS1-interacting patterner, and their interaction was further supported by reciprocal pull-down assay and co-immunoprecipitation. The inducible expression of MyD88 by V. harveyi suggested that the linker molecule MyD88 in innate immune response may play together with Na PRPS1 to coordinate the immune signaling in yellow drum in response to the pathogenic infection. We provide new insights into important functions of PRPS1, especially PRPS1 in the innate immunity of teleost fishes, which will benefit the development of marine fish aquaculture.

Rainbow trout ( Oncorhynchus mykiss ) serves as one of the most important commercial fish with an annual production of around 800,000 tonnes. However, infectious diseases, such as furunculosis caused by Aeromonas salmonicida infection, results in great economic loss in trout culture. The brain and kidney are two important organs associated with "sickness behaviors" and immunomodulation in response to disease. Therefore, we worked with 60 trout and investigated transcriptional responses and enrichment pathways between healthy and infected trout. We observed that furunculosis resulted in the activation of toll-like receptors with neuroinflammation and neural dysfunction in the brain, which might cause the "sickness behaviors" of infected trout including anorexia and lethargy. We also showed the salmonid-specific whole genome duplication contributed to duplicated colony stimulating factor 1 ( csf-1 ) paralogs, which play an important role in modulating brain immunomodulation. Enrichment analyses of kidneys showed up-regulated immunomodulation and down-regulated neural functions, suggesting an immune-neural interaction between the brain and kidney. Moreover, the kidney endocrine network was activated in response to A. salmonicida infection, further convincing the communications between endocrine and immune systems in regulating internal homeostasis. Our study provided a foundation for pathophysiological responses of the brain and kidney in response to furunculosis and potentially offered a reference for generating disease-resistant trout strains.

The protozoan parasite Cryptocaryon irritans causes marine white spot disease in a wide range of fish hosts, including gilthead seabream, a very sensitive species with great economic importance in the Mediterranean area. Thus, we aimed to evaluate the immunity of gilthead seabream after a severe natural outbreak of C. irritans . Morphological alterations and immune cell appearance in the gills were studied by light microscopy and immunohistochemical staining. The expression of several immune-related genes in the gills and head kidney were studied by qPCR, including inflammatory and immune cell markers, antimicrobial peptides (AMP), and cell-mediated cytotoxicity (CMC) molecules. Serum humoral innate immune activities were also assayed. Fish mortality reached 100% 8 days after the appearance of the C. irritans episode. Gill filaments were engrossed and packed without any space between filaments and included parasites and large numbers of undifferentiated and immune cells, namely acidophilic granulocytes. Our data suggest leukocyte mobilization from the head kidney, while the gills show the up-regulated transcription of inflammatory, AMPs, and CMC-related molecules. Meanwhile, only serum bactericidal activity was increased upon infection. A potent local innate immune response in the gills, probably orchestrated by AMPs and CMC, is triggered by a severe natural outbreak of C. irritans .

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, Ci Df, 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 Ci Df was found to be 753 bp, encoding 250 amino acids with a molecular mass of 27.06 kDa. Ci Df 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 Ci Df shares the highest identity (96%) and similarity (98%) with Df from Anabarilius grahami . Phylogenetic analysis indicated that Ci Df and other fish Dfs formed a distinct evolutionary branch. Similar to most Dfs from other vertebrates, the Ci Df gene structure is characterized by four introns and five exons. The incubation of recombinant Ci Df protein with grass carp serum significantly increased the C3b content, demonstrating the conserved function of Ci Df in the AP in promoting C3 cleavage, similar to Dfs in mammals. Ci Df 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 Ci Df 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 Ci Df 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 .

Megalocytivirus is an important viral pathogen to many farmed fishes, including Japanese flounder ( Paralichthys olivaceus ). In this study, we examined megalocytivirus-induced RNA responses in the spleen of flounder by high-throughput sequencing and integrative analysis of various RNA-seq data. A total of 1327 microRNAs (miRNAs), including 368 novel miRNAs, were identified, among which, 171 (named DEmiRs) exhibited significantly differential expressions during viral infection in a time-dependent manner. For these DEmiRs, 805 differentially expressed target mRNAs (DETmRs) were predicted, whose expressions not only significantly changed after megalocytivirus infection but were also negatively correlated with their paired DEmiRs. Integrative analysis of immune-related DETmRs and their target DEmiRs identified 12 hub DEmiRs, which, together with their corresponding DETmRs, formed an interaction network containing 84 pairs of DEmiR and DETmR. In addition to DETmRs, 19 DEmiRs were also found to regulate six key immune genes (mRNAs) differentially expressed during megalocytivirus infection, and together they formed a network consisting of 21 interactive miRNA-messenger RNA (mRNA) pairs. Further analysis identified 9434 circular RNAs (circRNAs), 169 of which (named DEcircRs) showed time-specific and significantly altered expressions during megalocytivirus infection. Integrated analysis of the DETmR-DEmiR and DEcircR-DEmiR interactions led to the identification of a group of competing endogenous RNAs (ceRNAs) constituted by interacting triplets of circRNA, miRNA, and mRNA involved in antiviral immunity. Together these results indicate that complicated regulatory networks of different types of non-coding RNAs and coding RNAs are involved in megalocytivirus infection.

MAF1 is a global suppressor of RNA polymerase III-dependent transcription, and is conserved from yeast to human. Growing evidence supports the involvement of MAF1 in the immune response of mammals, but its biological functions in fish are unknown. We isolated and characterized Maf1 from the olive flounder Paralichthys olivaceus ( PoMaf1 ). The coding region of PoMaf1 comprised 738 bp encoding a 245-amino-acid protein. The deduced PoMAF1 amino acid sequence shared features with those of MAF1 orthologues from vertebrates. PoMaf1 mRNA was detected in all tissues examined, and the levels were highest in eye and muscle tissue. The PoMaf1 mRNA level increased during early development. In addition, the PoMaf1 transcript level decreased during viral hemorrhagic septicemia virus (VHSV) infection of flounder hirame natural embryo (HINAE) cells. To investigate the role of PoMaf1 in VHSV infection, single-cell-derived PoMaf1 knockout HINAE cells were generated using the clustered regularly interspaced short palindromic repeats/CRISPR-associated-9 (CRISPR/Cas9) system, and cell clones with complete disruption of PoMaf1 were selected. PoMaf1 disruption increased the VHSV glycoprotein (G) mRNA levels during VHSV infection of HINAE cells, implicating PoMAF1 in the immune response to VSHV infection. To our knowledge, this is the first study to characterize fish Maf1 , which may play a role in the response to viral infection.

Microcystin-LR (MC-LR) is a potent hepatotoxin for which a substantial gap in knowledge persists regarding the underlying molecular mechanisms of liver toxicity and injury. Although long non-coding RNAs (lncRNAs) have been extensively studied in model organisms, our knowledge concerning the role of lncRNAs in liver injury is limited. Given that lncRNAs show low levels of sequence conservation, their role becomes even more unclear in non-model organisms without an annotated genome, like whitefish ( Coregonus lavaretus ). The objective of this study was to discover and profile aberrantly expressed polyadenylated lncRNAs that are involved in MC-LR-induced liver injury in whitefish. Using RNA sequencing (RNA-Seq) data, we de novo assembled a high-quality whitefish liver transcriptome. This enabled us to find 94 differentially expressed (DE) putative evolutionary conserved lncRNAs, such as MALAT1, HOTTIP, HOTAIR or HULC, and 4429 DE putative novel whitefish lncRNAs, which differed from annotated protein-coding transcripts (PCTs) in terms of minimum free energy, guanine-cytosine (GC) base-pair content and length. Additionally, we identified DE non-coding transcripts that might be 3' autonomous untranslated regions (3'UTRs) of mRNAs. We found both evolutionary conserved lncRNAs as well as novel whitefish lncRNAs that could serve as biomarkers of liver injury.

Ozone is a strong oxidant, and its use in aquaculture has been shown to improve water quality and fish health. At present, it is predominantly used in freshwater systems due to the high risk of toxic residual oxidant exposure in brackish water and seawater. Here, we report the effects of ozone on Atlantic salmon ( Salmo salar ) post-smolts (~100 g), in a brackish water (12 ppt) flow-through system. Salmon were exposed to oxidation reduction potential concentrations of 250 mV (control), 280 mV (low), 350 mV (medium), 425 mV (high) and 500 mV (very high). The physiological impacts of ozone were characterized by blood biochemical profiling, histopathologic examination and gene expression analysis in skin and gills. Fish exposed to 425 mV and higher showed ≥33% cumulative mortality in less than 10 days. No significant mortalities were recorded in the remaining groups. The skin surface quality and the thickness of the dermal and epidermal layers were not significantly affected by the treatments. On the other hand, gill histopathology showed the adverse effects of increasing ozone doses and the changes were more pronounced in the group exposed to 350 mV and higher. Cases of gill damages such as necrosis, lamellar fusion and hypertrophy were prevalent in the high and very high groups. Expression profiling of key biomarkers for mucosal health supported the histology results, showing that gills were significantly more affected by higher ozone doses compared to the skin. Increasing ozone doses triggered anti-oxidative stress and inflammatory responses in the gills, where transcript levels of glutathione reductase , copper/zinc superoxide dismutase , interleukin 1β and interleukin were significantly elevated. Heat shock protein 70 was significantly upregulated in the skin of fish exposed to 350 mV and higher. Bcl-2 associated x protein was the only gene marker that was significantly upregulated by increasing ozone doses in both mucosal tissues. In conclusion, the study revealed that short-term exposure to ozone at concentrations higher than 350 mV in salmon in brackish water resulted in significant health and welfare consequences, including mortality and gill damages. The results of the study will be valuable in developing water treatment protocols for salmon farming.

MicroRNAs (miRNAs) are non-coding regulatory RNAs that play a vital part in the host immune response to pathogen infection. Japanese flounder ( Paralichthys olivaceus ) is an important aquaculture fish species that has suffered from bacterial diseases, including that caused by Vibrio anguillarum infection. In a previous study, we examined the messenger RNA (mRNA) expression profiles of flounder during V. anguillarum infection and identified 26 hub genes in the flounder immune response. In this study, we performed the micro-transcriptome analysis of flounder spleen in response to V. anguillarum infection at 3 different time points. Approximately 277 million reads were obtained, from which 1218 miRNAs were identified, including 740 known miRNAs and 478 novel miRNAs. Among the miRNAs, 206 were differentially expressed miRNAs (DEmiRs), and 104 of the 206 DEmiRs are novel miRNAs identified for the first time. Most of the DEmiRs were strongly time-dependent. A total of 1355 putative target genes of the DEmiRs (named DETGs) were identified based on integrated analysis of miRNA-mRNA expressions. The DETGs were enriched in multiple functional categories associated with immunity. Thirteen key DEmiRs and 66 immune DETGs formed an intricate regulatory network constituting 106 pairs of miRNAs and DETGs that span five immune pathways. Furthermore, seven of the previously identified hub genes were found to be targeted by 73 DEmiRs, and together they formed interlinking regulatory networks. These results indicate that V. anguillarum infection induces complicated miRNA response with extensive influences on immune gene expression in Japanese flounder.

Proliferative kidney disease is an emerging disease among salmonids in Europe and North America caused by the myxozoan parasite Tetracapsuloides bryosalmonae. The decline of endemic brown trout ( Salmo trutta ) in the Alpine streams of Europe is fostered by T. bryosalmonae infection. Toll-like receptors (TLRs) are a family of pattern recognition receptors that acts as sentinels of the immune system against the invading pathogens. However, little is known about the TLRs' response in salmonids against the myxozoan infection. In the present study, we identified and evaluated TLR1, TLR19, and TLR13-like genes of brown trout using data-mining and phylogenetic analysis. The expression pattern of TLRs was examined in the posterior kidney of brown trout infected with T. bryosalmonae at various time points. Typical Toll/interleukin-1 receptor protein domain was found in all tested TLRs. However, TLR13-like chr2 had a short amino acid sequence with no LRR domain. Phylogenetic analysis illustrated that TLR orthologs are conserved across vertebrates. Similarly, a conserved synteny gene block arrangement was observed in the case of TLR1 and TLR19 across fish species. Interestingly, all tested TLRs showed their maximal relative expression from 6 to 10 weeks post-exposure to the parasite. Our results suggest that these TLRs may play an important role in the innate defense mechanism of brown trout against the invading T. bryosalmonae ., Competing Interests: The authors declare no conflict of interest.

Parasitic sea lice (e.g., Lepeophtheirus salmonis ) cause costly outbreaks in salmon farming. Molecular insights into parasite-induced host responses will provide the basis for improved management strategies. We investigated the early transcriptomic responses in pelvic fins of Atlantic salmon parasitized with chalimus I stage sea lice. Fin samples collected from non-infected (i.e. pre-infected) control (PRE) and at chalimus-attachment sites (ATT) and adjacent to chalimus-attachment sites (ADJ) from infected fish were used in profiling global gene expression using 44 K microarrays. We identified 6568 differentially expressed probes (DEPs, FDR < 5%) that included 1928 shared DEPs between ATT and ADJ compared to PRE. The ATT versus ADJ comparison revealed 90 DEPs, all of which were upregulated in ATT samples. Gene ontology/pathway term network analyses revealed profound changes in physiological processes, including extracellular matrix (ECM) degradation, tissue repair/remodeling and wound healing, immunity and defense, chemotaxis and signaling, antiviral response, and redox homeostasis in infected fins. The QPCR analysis of 37 microarray-identified transcripts representing these functional themes served to confirm the microarray results with a significant positive correlation ( p < 0.0001). Most immune/defense-relevant transcripts were downregulated in both ATT and ADJ sites compared to PRE, suggesting that chalimus exerts immunosuppressive effects in the salmon's fins. The comparison between ATT and ADJ sites demonstrated the upregulation of a suite of immune-relevant transcripts, evidencing the salmon's attempt to mount an anti-lice response. We hypothesize that an imbalance between immunomodulation caused by chalimus during the early phase of infection and weak defense response manifested by Atlantic salmon makes it a susceptible host for L. salmonis ., Competing Interests: Authors S.S. and R.G.T. in the representation of Cargill Innovation, participated in the data interpretation and discussion, but had no role in the design of the microarray experiment and transcript expression experiment, the data collection and analysis, the preparation of this manuscript, and the decision to submit the manuscript for publication. All authors declare no conflict of interest.

Rock bream iridovirus (RBIV) is a notorious agent that causes high mortality in aquaculture of rock bream ( Oplegnathus fasciatus ). Despite severity of this virus, no transcriptomic studies on RBIV-infected rock bream that can provide fundamental information on protective mechanism against the virus have been reported so far. This study aimed to investigate physiological mechanisms between host and RBIV through transcriptomic changes in the spleen based on RNA-seq. Depending on infection intensity and sampling time point, fish were divided into five groups: uninfected healthy fish at week 0 as control (0C), heavy infected fish at week 0 (0H), heavy mixed RBIV and bacterial infected fish at week 0 (0MH), uninfected healthy fish at week 3 (3C), and light infected fish at week 3 (3L). We explored clusters from 35,861 genes with Fragments Per Kilo-base of exon per Million mapped fragments (FPKM) values of 0.01 or more through signed co-expression network analysis using WGCNA package. Nine of 22 modules were highly correlated with viral infection (|gene significance (GS) vs. module membership (MM) |> 0.5, p -value < 0.05). Expression patterns in selected modules were divided into two: heavy infected (0H and 0MH) and control and light-infected groups (0C, 3C, and 3L). In functional analysis, genes in two positive modules (5448 unigenes) were enriched in cell cycle, DNA replication, transcription, and translation, and increased glycolysis activity. Seven negative modules (3517 unigenes) built in this study showed significant decreases in the expression of genes in lymphocyte-mediated immune system, antigen presentation, and platelet activation, whereas there was significant increased expression of endogenous apoptosis-related genes. These changes lead to RBIV proliferation and failure of host defense, and suggests the importance of blood cells such as thrombocytes and B cells in rock bream in RBIV infection. Interestingly, a hub gene, pre-mRNA processing factor 19 (PRPF19) showing high connectivity (kME), and expression of this gene using qRT-PCR was increased in rock bream blood cells shortly after RBIV was added. It might be a potential biomarker for diagnosis and vaccine studies in rock bream against RBIV. This transcriptome approach and our findings provide new insight into the understanding of global rock bream-RBIV interactions including immune and pathogenesis mechanisms.

Hepcidin is a liver-derived peptide hormone that is related to iron balance and immunity in humans. However, its function in Siniperca chuatsi has not been well elucidated. In this study, we analyzed the expression and function of the S. chuatsi hepcidin ( Sc-hep ) gene. Sc-hep was specifically expressed in the liver and appeared to be one of the most highly expressed genes in the liver. After spleen and kidney necrosis virus (ISKNV) infection and lipopolysaccharide (LPS) and polyinosinic-polycytidylic acid (Poly I:C) stimulation, the expression of Sc-hep in the liver increased by approximately 110, 6500, and 225 times, respectively. After ferrous sulfate (FS) injection, the expression of Sc-hep in the liver increased approximately 520-fold. We found that miR-19c-5p could inhibit Sc-hep expression. Five CpG dinucleotides distributed in the promoter region showed no differential methylation between the liver and the stomach, both presenting high methylation rates. After FS or LPS injection, the expression of three iron balance-related genes ( FPN1 , TFR1 , and FTN ) and five immune-related cytokine genes ( IL-1 β, IL8 , TNF- α, TLR22 , and SOCS3 ) significantly changed. These results indicate that Sc-hep participates in the regulation of iron balance and plays an important role in the immune system. Sc-hep increased approximately 52-fold when mandarin fish were domesticated with artificial diets. Sc-hep might be used as a real-time biomarker of mandarin fish liver because its expression markedly varies under different physiological conditions.

Integrin β-1 (ITGB1) is a transmembrane protein belonging to the integrin family and it plays an important role in viral entry. In this study, the itgb1b gene of the rare minnow, Gobiocypris rarus , was cloned and analyzed. To investigate the possible role of itgb1b on grass carp reovirus (GCRV) infection, we generated an ITGB1b-deficient rare minnow (ITGB1b -/- ) using the CRISPR/Cas9 system. Following stimulation with GCRV, the survival time of the -ITGB1b -/- rare minnows was extended in comparison to the wild-type minnows. Moreover, the relative copy number of GCRV and the level of clathrin-mediated endocytosis-associated and apoptosis-related gene expression in the ITGB1b -/- rare minnows was significantly lower than that of the wild-type minnows. These results suggested that the absence of itgb1b reduced viral entry efficiency and the expression of apoptosis-related genes. Moreover, the data suggested that itgb1b played an important role in mediating the entry of viruses into the cells via clathrin. Therefore, these findings provide novel insight into the function of itgb1b in the process of GCRV infection.

Grass carp reovirus (GCRV) causes huge economic loss to the grass carp cultivation industry but the mechanism remains largely unknown. In this study, we investigated the global and complement gene-specific DNA methylation in grass carp after GCRV infection aimed to uncover the mechanism underlying GCRV infection. The global DNA methylation level was increased after GCRV infection. Expression levels of enzymes involved in DNA methylation including DNA methyltransferase (DNMT), ten-eleven translocation proteins (TETs), and glycine N -methyltransferase (GNMT) were significantly altered after GCRV infection. In order to investigate the relationship between the gene expression level and DNA methylation level, two representative complement genes, complement component 3 ( C3 ) and kininogen-1 ( KNG1 ), were selected for further analysis. mRNA expression levels of the two genes were significantly increased at 5 and 7 days after GCRV infection, whereas the DNA methylation level at the 5' flanking regions of the two genes were down-regulated at the same time-points. Moreover, a negative correlation was detected between gene expression levels and DNA methylation levels of the two genes. Therefore, the current data revealed a global and complement gene-specific DNA methylation profile after GCRV infection. Our study would provide new insights into understanding the mechanism underlying GCRV infection., Competing Interests: The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

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

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., Competing Interests: The authors declare no conflict of interest.

Lymphocystis disease virus (LCDV) infection may induce a variety of host gene expression changes associated with disease development; however, our understanding of the molecular mechanisms underlying host-virus interactions is limited. In this study, RNA sequencing (RNA-seq) was employed to investigate differentially expressed genes (DEGs) in the gill of the flounder ( Paralichthys olivaceus ) at one week post LCDV infection. Transcriptome sequencing of the gill with and without LCDV infection was performed using the Illumina HiSeq 2500 platform. In total, RNA-seq analysis generated 193,225,170 clean reads aligned with 106,293 unigenes. Among them, 1812 genes were up-regulated and 1626 genes were down-regulated after LCDV infection. The DEGs related to cellular process and metabolism occupied the dominant position involved in the LCDV infection. A further function analysis demonstrated that the genes related to inflammation, the ubiquitin-proteasome pathway, cell proliferation, apoptosis, tumor formation, and anti-viral defense showed a differential expression. Several DEGs including β actin , toll-like receptors, cytokine-related genes, antiviral related genes, and apoptosis related genes were involved in LCDV entry and immune response. In addition, RNA-seq data was validated by quantitative real-time PCR. For the first time, the comprehensive gene expression study provided valuable insights into the host-pathogen interaction between flounder and LCDV., Competing Interests: The authors declare no conflict of interest.

Grass carp hemorrhagic disease, caused by the grass carp reovirus (GCRV), is a major disease that hampers the development of grass carp aquaculture in China. The mechanism underlying GCRV infection is still largely unknown. Circular RNAs (circRNAs) are important regulators involved in various biological processes. In the present study, grass carp were infected with GCRV, and spleen samples were collected at 0 (control), 1, 3, 5, and 7 days post-infection (dpi). Samples were used to construct and sequence circRNA libraries, and a total of 5052 circRNAs were identified before and after GCRV infection, of which 41 exhibited differential expression compared with controls. Many parental genes of the differentially expressed circRNAs are involved in metal ion binding, protein ubiquitination, enzyme activity, and nucleotide binding. Moreover, 72 binding miRNAs were predicted from the differentially expressed circRNAs, of which eight targeted genes were predicted to be involved in immune responses, blood coagulation, hemostasis, and complement and coagulation cascades. Upregulation of these genes may lead to endothelial and blood cell damage and hemorrhagic symptoms. Our results indicate that an mRNA-miRNA-circRNA network may be present in grass carp infected with GCRV, providing new insight into the mechanism underlying grass carp reovirus infection., Competing Interests: The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

This study describes the development and validation of an enriched oligonucleotide-microarray platform for Sparus aurata (SAQ) to provide a platform for transcriptomic studies in this species. A transcriptome database was constructed by assembly of gilthead sea bream sequences derived from public repositories of mRNA together with reads from a large collection of expressed sequence tags (EST) from two extensive targeted cDNA libraries characterizing mRNA transcripts regulated by both bacterial and viral challenge. The developed microarray was further validated by analysing monocyte/macrophage activation profiles after challenge with two Gram-negative bacterial pathogen-associated molecular patterns (PAMPs; lipopolysaccharide (LPS) and peptidoglycan (PGN)). Of the approximately 10,000 EST sequenced, we obtained a total of 6837 EST longer than 100 nt, with 3778 and 3059 EST obtained from the bacterial-primed and from the viral-primed cDNA libraries, respectively. Functional classification of contigs from the bacterial- and viral-primed cDNA libraries by Gene Ontology (GO) showed that the top five represented categories were equally represented in the two libraries: metabolism (approximately 24% of the total number of contigs), carrier proteins/membrane transport (approximately 15%), effectors/modulators and cell communication (approximately 11%), nucleoside, nucleotide and nucleic acid metabolism (approximately 7.5%) and intracellular transducers/signal transduction (approximately 5%). Transcriptome analyses using this enriched oligonucleotide platform identified differential shifts in the response to PGN and LPS in macrophage-like cells, highlighting responsive gene-cassettes tightly related to PAMP host recognition. As observed in other fish species, PGN is a powerful activator of the inflammatory response in S. aurata macrophage-like cells. We have developed and validated an oligonucleotide microarray (SAQ) that provides a platform enriched for the study of gene expression in S. aurata with an emphasis upon immunity and the immune response.

In the present study, IL-1β cDNA was identified and analyzed from largemouth bass ( Micropterus salmoides ). Full length IL-1β mRNA was obtained using Rapid Amplification of cDNA Ends (RACE), which contains 78 bp 3'-UTR, a 455 bp 5'-UTR, and an open reading frame (ORF) of 702 bp coding for 233 amino acid residues. The molecular weight and theoretical isoelectric point of largemouth bass IL-1β protein was predicted to be 26.7 kDa and 6.08 respectively. A largemouth bass IL-1β phylogenetic analysis showed a close relation to the IL-1βs of striped trumpeter ( Latris lineata ), Chinese perch ( Siniperca chuatsi ), and Japanese sea bass ( Lateolabrax japonicus ). Peptidoglycan upregulated IL-1β in the spleen and head kidney, while lipopolysaccharide upregulated detectable levels of IL-1β in the spleen only. Largemouth bass, challenged with Nocardia seriolae (1.0 × 10⁶ cfu/mL), showed a significant increase in IL-1β at 3 and 5 days post infection (dpi) in the spleen, while in the head kidney significant expression was found at 2 and 3 dpi, peaking at 3 dpi. Furthermore, tumor necrosis factor α (TNF-α) showed significantly higher expression in the spleen at 3 and 5 dpi, and in the head kidney at 1 and 3 dpi, with expression decreasing at 5 dpi in both tissues., Competing Interests: The authors declare no conflict of interest.

MicroRNAs (miRNAs) play important roles in regulation of many biological processes in eukaryotes, including pathogen infection and host interactions. Flavobacterium columnare (FC) infection can cause great economic loss of common carp (Cyprinus carpio) which is one of the most important cultured fish in the world. However, miRNAs in response to FC infection in common carp has not been characterized. To identify specific miRNAs involved in common carp infected with FC, we performed microRNA sequencing using livers of common carp infected with and without FC. A total of 698 miRNAs were identified, including 142 which were identified and deposited in the miRbase database (Available online: http://www.mirbase.org/) and 556 had only predicted miRNAs. Among the deposited miRNAs, eight miRNAs were first identified in common carp. Thirty of the 698 miRNAs were differentially expressed miRNAs (DIE-miRNAs) between the FC infected and control samples. From the DIE-miRNAs, seven were selected randomly and their expression profiles were confirmed to be consistent with the microRNA sequencing results using RT-PCR and qRT-PCR. In addition, a total of 27,363 target genes of the 30 DIE-miRNAs were predicted. The target genes were enriched in five Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including focal adhesion, extracellular matrix (ECM)-receptor interaction, erythroblastic leukemia viral oncogene homolog (ErbB) signaling pathway, regulation of actin cytoskeleton, and adherent junction. The miRNA expression profile of the liver of common carp infected with FC will pave the way for the development of effective strategies to fight against FC infection.

MicroRNAs (miRNAs) play important roles in mediating multiple biological processes in eukaryotes and are being increasingly studied to evaluate their roles associated with cellular changes following viral infection. Snakehead fish Vesiculovirus (SHVV) has caused mass mortality in snakehead fish during the past few years. To identify specific miRNAs involved in SHVV infection, we performed microRNA deep sequencing on a snakehead fish cell line (SSN-1) with or without SHVV infection. A total of 205 known miRNAs were identified when they were aligned with the known zebrafish miRNAs, and nine novel miRNAs were identified using MiRDeep2 software. Eighteen and 143 of the 205 known miRNAs were differentially expressed at three and 24 h post-infection (poi), respectively. From the differentially-expressed miRNAs, five were randomly selected to validate their expression profiles using quantitative reverse transcription polymerase chain reaction (qRT-PCR), and their expression profiles were consistent with the microRNA sequencing results. In addition, the target gene prediction of the SHVV genome was performed for the differentially-expressed host miRNAs, and a total of 10 and 58 differentially-expressed miRNAs were predicted to bind to the SHVV genome at three and 24 h poi, respectively. The effects of three selected miRNAs (miR-130-5p, miR-214 and miR-216b) on SHVV multiplication were evaluated using their mimics and inhibitors via qRT-PCR and Western blotting. The results showed that all three miRNAs were able to inhibit the multiplication of SHVV; whereas the mechanisms underlying the SHVV multiplication inhibited by the specific miRNAs need to be further characterized in the future.

Channel catfish (Ictalurus punctatus) is an important agricultural fish that has been plagued by Streptococcus iniae (S. iniae) infections in recent years, some of them severe. C5a peptidase is an important virulent factor of S. iniae. In this study, the subunit vaccine containing the truncated part of C5a peptidase (pSCPI) was mixed with aluminum hydroxide gel (AH), propolis adjuvant (PA), and Freund's Incomplete Adjuvant (FIA). The immunogenicity of the pSCPI was detected by Western-blot in vitro. The relative percent survival (RPS), lysozyme activity, antibody titers, and the expression of the related immune genes were monitored in vivo to evaluate the immune effects of the three different adjuvants. The results showed that pSCPI exerted moderate immune protection (RPS = 46.43%), whereas each of the three adjuvants improved the immune protection of pSCPI. The immunoprotection of pSCPI + AH, pSCPI + PA, and pSCPI + FIA was characterized by RPS values of 67.86%, 75.00% and, 85.71%, respectively. Further, each of the three different adjuvanted pSCPIs stimulated higher levels of lysozyme activity and antibody titers than the unadjuvanted pSCPI and/or PBS buffer. In addition, pSCPI + FIA and pSCPI + PA induced expression of the related immune genes under investigation, which was substantially higher than the levels stimulated by PBS. pSCPI + AH significantly stimulated the induction of MHC II β, CD4-L2, and IFN-γ, while it induced slightly higher production of TNF-α and even led to a decrease in the levels of IL-1β, MHC I α, and CD8 α. Therefore, we conclude that compared with the other two adjuvants, FIA combined with pSCPI is a more promising candidate adjuvant against S. iniae in channel catfish.

The 27.8 kDa membrane protein expressed in flounder (Paralichthys olivaceus) gill cells was proved to be a receptor mediating lymphocystis disease virus (LCDV) infection. In this study, SDS-PAGE and Western blotting demonstrated that 27.8 kDa receptor (27.8R) was shared by flounder and turbot (Scophthalmus maximus). Indirect immunofluorescence assay (IIFA) and immunohistochemistry showed that 27.8R was widely expressed in tested tissues of healthy turbot. The indirect enzyme-linked immunosorbent assay indicated that 27.8R expression was relatively higher in stomach, gill, heart, and intestine, followed by skin, head kidney, spleen, blood cells, kidney and liver, and lower in ovary and brain in healthy turbot, and it was significantly up-regulated after LCDV infection. Meanwhile, real-time quantitative PCR demonstrated that LCDV was detected in heart, peripheral blood cells, and head kidney at 3 h post infection (p.i.), and then in other tested tissues at 12 h p.i. LCDV copies increased in a time-dependent manner, and were generally higher in the tissues with higher 27.8R expression. Additionally, IIFA showed that 27.8R and LCDV were detected at 3 h p.i. in some leukocytes. These results suggested that 27.8R also served as a receptor in turbot, and LCDV can infect some leukocytes which might result in LCDV spreading to different tissues in turbot.