Identification and characterization of Cynoglossus semilaevis microRNA response to Vibrio anguillarum infection through high-throughput sequencing

Characterizing and identifying of miRNAs involved in berberine modulating glucose metabolism of Megalobrama amblycephala

The present study, as one part of a larger project that aimed to investigate the effects of dietary berberine (BBR) on fish growth and glucose regulation, mainly focused on whether miRNAs involve in BBR's modulation of glucose metabolism in fish. Blunt snout bream Megalobrama amblycephala (average weight of 20.36 ± 1.44 g) were exposed to the control diet (NCD, 30% carbohydrate), the high-carbohydrate diet (HCD, 43% carbohydrate) and the berberine diet (HCB, HCD supplemented with 50 mg/kg BBR). After 10 weeks' feeding trial, intraperitoneal injection of glucose was conducted, and then, the plasma and liver were sampled at 0 h, 1 h, 2 h, 6 h, and 12 h. The results showed the plasma glucose levels in all groups rose sharply and peaked at 1 h after glucose injection. Unlike the NCD and HCB groups, the plasma glucose in the HCD group did not decrease after 1 h, while remained high level until at 2 h. The NCD group significantly increased liver glycogen content at times 0-2 h compared to the other two groups and then liver glycogen decreased sharply until at times 6-12 h. To investigate the role of BBR that may cause the changes in plasma glucose and liver glycogen, miRNA high-throughput sequencing was performed on three groups of liver tissues at 2 h time point. Eventually, 20 and 12 differentially expressed miRNAs (DEMs) were obtained in HCD vs NCD and HCB vs HCD, respectively. Through function analyzing, we found that HCD may affect liver metabolism under glucose loading through the NF-κB pathway; and miRNAs regulated by BBR mainly play roles in adipocyte lipolysis, niacin and nicotinamide metabolism, and amino acid transmembrane transport. In the functional exploration of newly discovered novel:Chr12_18892, we found its target gene, adenylate cyclase 3 (adcy3), was widely involved in lipid decomposition, amino acid metabolism, and other pathways. Furthermore, a targeting relationship of novel:Chr12_18892 and adcy3 was confirmed by double luciferase assay. Thus, BBR may promote novel:Chr12_18892 to regulate the expression of adcy3 and participate in glucose metabolism.

miR-144 affects the immune response and activation of inflammatory responses in Cynoglossus semilaevis by regulating the expression of CsMAPK6

MicroRNAs are increasingly recognized for their pivotal role in the immune system, yet the specific regulatory functions of fish-derived microRNAs remain largely unexplored. In this research, we discovered a novel miRNA, Cse-miR-144, in the Chinese tongue sole (Cynoglossus semilaevis), characterized by a 73-base pair precursor and a 21-nucleotide mature sequence. Our findings revealed that the expression of Cse-miR-144 was notably inhibited by various Vibrio species. Utilizing bioinformatics and dual-luciferase assay techniques, we established that the pro-inflammatory cytokine gene CsMAPK6 is a direct target of Cse-miR-144. Subsequent in vitro and in vivo western blotting analyses confirmed that Cse-miR-144 can effectively reduce the protein levels of CsMAPK6 post-transcriptionally. Moreover, CsMAPK6 is known to be involved in the activation of the Nuclear Factor kappa-light-chain-enhancer of activated B cells (NF-kB). Additional investigations using qPCR and ELISA demonstrated that suppression of Cse-miR-144 leads to an upsurge in the liver mRNA levels of various immune genes (including MYD88, TRAF6, NF-κB, TRAF2, TRAF3, and TNF), alongside a marked increase in the production and secretion of pro-inflammatory cytokines (IL-1β, IL-6, and IL-8) in the bloodstream of C. semilaevis. These findings collectively underscore the potential of Cse-miR-144 as a key inhibitor of CsMAPK and its crucial role in modulating the immune and inflammatory responses in teleost fish. Compared to the siRNA, miRNA is a better tool in controlling the expression of target gene with a lower cost.

Effects of miR-722 on gene expression and alternative splicing in the liver of half-smooth tongue sole after infection with Vibrio anguillarum

MicroRNAs play crucial roles in various biological processes, including but not limited to differentiation, development, disease, and immunity. However, their immunoregulatory roles in half-smooth tongue sole are lacking. Our previous studies indicated that miR-722 could target C5aR1 to modulate the complement pathway to alleviate inflammatory response and even affect the mortality after the bacterial infection with Vibrio anguillarum. Driven by the purpose of revealing the underlying mechanisms, in this study, we investigated the effects of miR-722 on the gene expression and alternative splicing (AS) in the liver of half-smooth tongue sole after Vibrio anguillarum infection, with the approach of miR-722 overexpression/silencing and subsequent RNA-seq. Among the different comparisons, the I group (miR-722 inhibitor and V. anguillarum) versus blank control (PBS) exhibited the highest number of differentially expressed genes (DEGs), suggesting that the immune response was overactivated after inhibiting the miR-722. In addition, enrichment analyses were performed to reveal the functions of DEGs and differential AS (DAS) genes, reflecting the enrichment of RNA splicing and immune-related pathways including NF-κB and T cell receptor signaling pathway. Comparing the M group (miR-722 mimic and V. anguillarum) with the negative control (random sequence and V. anguillarum), two immune-related genes, cd48 and mapk8, were differentially expressed, of which mapk8 was also differentially spliced, indicating their importance in the immune response. Furthermore, representative gene analysis was performed, suggesting their corresponding functional changes due to AS. To verify the RNA-seq data, quantitative real-time PCR was employed with twenty pairs of primers for DEGs and DAS events. Overall, our results demonstrated that miR-722 could mediate the transcriptome-wide changes of gene expression and AS in half-smooth tongue sole, and provided insights into the regulatory role of miR-722 in immune responses, laying the foundation for further functional analyses and practical applications in aquaculture.

Regulatory mechanism of miR-722 on C5aR1 and its functions against bacterial inflammation in half-smooth tongue sole (Cynoglossus semilaevis)

MicroRNAs (miRNAs) are small non-coding RNAs involved in various biological processes, including immunity. Previously, we investigated the miRNAs of half-smooth tongue sole (Cynoglossus semilaevis) and found that miR-722 (designated Cse-miR-722) was significantly differentially expressed after infection with Vibrio anguillarum, reflecting its importance in immune response. Our preliminary bioinformatic analysis suggested that Cse-miR-722 could target C5aR1 (designated CsC5aR1), which was known to play crucial roles in complement activation and inflammatory response, as a receptor of C5a. However, the underlying mechanisms of their interactions and specific functions in inflammatory and immune response are still enigmas. In this study, we successfully cloned the precursor sequence of Cse-miR-722 (94 bp) and the full length of CsC5aR1 (1541 bp, protein molecular weight 39 kDa). The target gene of Cse-miR-722 was verified as CsC5aR1 by a dual luciferase reporter assay, and Cse-miR-722 was confirmed to regulate CsC5aR1 at the protein level using quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting. The expression of CsC5aR1 and Cse-miR-722 in liver cells and four immune tissues of half-smooth tongue sole changed significantly after LPS stimulation and infection with V. anguillarum. To explore the functional role of Cse-miR-722 in half-smooth tongue sole, we performed both in vitro and in vivo experiments. Cse-miR-722 was observed to affect phagocytosis and respiratory burst activity of macrophages by regulating CsC5aR1 in half-smooth tongue sole. Furthermore, we found that Cse-miR-722 regulated the expression of CsC5aR1, CsC5a, and the inflammatory factors CsIL1-β, CsIL6, CsIL8, and CsTNF-α both in vitro and in vivo. In addition, Cse-miR-722 reduced mortality and pathological damage. This study clarified the regulatory mechanism of Cse-miR-722 on CsC5aR1 and provided insight into the regulatory roles of Cse-miR-722 in immune responses, laying a theoretical foundation for the feasibility of using miR-722 to prevent and control bacterial diseases in teleost.

Diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) have altered microRNA responses in immune tissues after infection with Vibrio anguillarum

Production of sterile fishes through artificial retention of a third set of chromosomes (triploidy) is a sustainable alternative for aquaculture since it reduces escapee pressure on wild populations. However, these fishes have reduced survival in stressful conditions and in response to infection. In this study, the impact of Vibrio anguillarum infection on diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) was investigated to identify if there was any significant immune regulation by microRNAs (miRNA). Small RNAs from hindgut, head kidney, and spleen were sequenced to determine if miRNA transcript abundance was altered due to ploidy and infection in nine-month old full-sibling diploids and triploids. All three tissues had differentially expressed miRNA prior to infection, indicating subtle changes in epigenetic regulation due to increased ploidy. Additionally, miRNA were altered by infection, but there was only a difference in spleen miRNA expression between diploids and triploids at three days of infection. Furthermore, one miRNA (ssa-miR-2188-3p) was confirmed as having an altered response to infection in triploids compared to diploids, implicating potential immune dysregulation due to increased ploidy. The miRNAs identified in this study are predicted to target immune pathways, providing evidence for their importance in regulating responses to pathogens. This study is the first to investigate how increased ploidy alters miRNA expression in response to infection. Additionally, it provides evidence for epigenetic dysregulation in triploid fishes, which may contribute to their poor performance in response to stress.

A novel teleost microRNA regulates autophagy and NF-κB activation during bacterial infection

MicroRNAs (miRNAs) are a class of non-coding RNAs with regulatory functions in many cellular processes, including immune defense. In this study, we identified novel-m0089-3p, a novel miRNA with unknown function, in the teleost fish Japanese flounder (Paralichthys olivaceus) and investigated its immune function. Novel-m0089-3p was found to target the autophagy-associated gene ATG7 and negatively regulate ATG7 expression via interaction with the 3' UTR of ATG7. During the infection of the bacterial pathogen Edwardsiella tarda, novel-m0089-3p expression was induced in flounder, which in turn repressed ATG7 expression. Overexpression of novel-m0089-3p or blocking ATG7 expression inhibited autophagy and promoted the intracellular replication of E. tarda. Novel-m0089-3p overexpression, as well as E. tarda infection, activated NF-κB and stimulated the expression of inflammatory cytokines. Together these results revealed an important role of novel-m0089-3p in response to bacterial infection.

Identification and evaluation of potential probiotics against skin-ulceration disease in the Chinese tongue sole (Cynoglossus semilaevis)

In this study, three highly pathogenic bacterial strains (Vibrio harveyi TB6, Vibrio alginolyticus TN1, and Vibrio parahaemolyticus TN3) were isolated from skin ulcers and intestines of diseased Chinese tongue sole (Cynoglossus semilaevis). The bacteria were investigated using hemolytic activity tests, in vitro co-culture with intestinal epithelial cells, and artificial infection of C. semilaevis. A further 126 strains were isolated from the intestines of healthy C. semilaevis. The three pathogens were used as indicator bacteria, and the antagonistic strains were identified from the 126 strains. The activities of exocrine digestive enzymes in the strains were also tested. Four strains with antibacterial and digestive enzyme activities were obtained and the best strains, Bacillus subtilis Y2 and Bacillus amyloliquefaciens Y9, were selected according to their ability to protect epithelial cells from infection. In addition, the effects of strains Y2 and Y9 at the individual level were investigated, finding that the activities of the immune-related enzymes superoxide dismutase, catalase, acid phosphatase, and peroxidase were significantly increased in the sera of the treatment group compared with the control group (p < 0.05). The specific growth rate (SGR, %) was also increased, especially in the Y2 group, and was significantly higher compared with the controls (p < 0.05). The result of the artificial infection test showed that the cumulative mortality within 72 h in the Y2 group was the lowest (50.5%), and in the Y9 group (68.5%) it was significantly lower than that in the control group (100%) (p < 0.05). Analysis of the intestinal microbial communities indicated that Y2 and Y9 could alter the composition of the intestinal flora, increasing both species richness and evenness, and inhibiting the growth of Vibrio in the intestine. These results suggested food supplemented with Y2 and Y9 could improve both immune function and disease resistance, as well as have a positive effect on the growth performance and the intestinal morphology of C. semilaevis.

Cromileptes altivelis microRNA Transcriptome Analysis upon Nervous Necrosis Virus (NNV) Infection and the Effect of cal-miR-155 on Cells Apoptosis and Virus Replication

MicroRNAs (miRNAs) could regulate various biological processes. Nervous necrosis virus (NNV) is one of the primary germs of the Humpback grouper (Cromileptes altivelis), a commercial fish of great importance for Asian aquaculture. However, there is limited available information on the host-virus interactions of C. altivelis. miRNAs have been shown to play key roles in the host response to infection by a variety of pathogens. To better understand the regulatory mechanism of miRNAs, we constructed miRNA transcriptomes and identified immune-related miRNAs of C. altivelis spleen in response to NNV infection. Reads from the three libraries were mapped onto the Danio rerio reference genome. As a result, a total of 942 mature miRNAs were determined, with 266 known miRNAs and 676 novel miRNAs. Among them, thirty-two differentially expressed miRNAs (DEmiRs) were identified compared to the PBS control. These DEmiRs were targeted on 895 genes, respectively, by using miRanda v3.3a. Then, 14 DEmiRs were validated by qRT-PCR and showed consistency with those obtained from high-throughput sequencing. In order to study the relationship between viral infection and host miRNA, a cell line from C. altivelis brain (CAB) was used to examine the expressions of five known DEmiRs (miR-132-3p, miR-194a, miR-155, miR-203b-5p, and miR-146) during NNV infection. The results showed that one miRNA, cal-miRNA-155, displayed significantly increased expression in response to the virus infection. Subsequently, it was proved that overexpression of cal-miR-155 enhanced cell apoptosis with or without NNV infection and inhibited virus replication in CAB cells. Oppositely, the cal-miRNA-155 inhibitor markedly suppressed apoptosis in CAB cells. The results of the apoptosis-related genes mRNA expression also showed the regulation of cal-miR-155 on the apoptosis process in CAB cells. These findings verify that miR-155 might exert a function as a pro-apoptotic factor in reply to NNV stimulation in CAB cells and help us further study the molecular mechanisms of the pathogenesis of NNV in C. altivelis.

Expression characteristics and interaction networks of microRNAs in spleen tissues of grass carp (Ctenopharyngodon idella)

The spleen is an important immune organ in fish. MicroRNAs (miRNAs) have been shown to play an important role in the regulation of immune function. However, miRNA expression profiles and their interaction networks associated with the postnatal late development of spleen tissue are still poorly understood in fish. The grass carp (Ctenopharyngodon idella) is an important economic aquaculture species in China. Here, two small RNA libraries were constructed from the spleen tissue of healthy grass carp at one-year-old and three-year-old. A total of 324 known conserved miRNAs and 9 novel miRNAs were identified by using bioinformatic analysis. Family analysis showed that 23 families such as let-7, mir-1, mir-10, mir-124, mir-8, mir-7, mir-9, and mir-153 were highly conserved between vertebrates and invertebrates. In addition, 14 families such as mir-459, mir-430, mir-462, mir-7147, mir-2187, and mir-722 were present only in fish. Expression analysis showed that the expression patterns of miRNAs in the spleen of one-year-old and three-year-old grass carp were highly consistent, and the percentage of miRNAs with TPM > 100 was above 39%. Twenty significant differentially expressed (SDE) miRNAs were identified. Gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that these SDE miRNAs were primarily involved in erythrocyte differentiation, lymphoid organ development, immune response, lipid metabolic process, the B cell receptor signaling pathway, the T cell receptor signaling pathway, and the PPAR signaling pathway. In addition, the following miRNA-mRNA interaction networks were constructed: immune and hematopoietic, cell proliferation and differentiation, and lipid metabolism. This study determined the miRNA transcriptome as well as miRNA-mRNA interaction networks in normal spleen tissue during the late development stages of grass carp. The results expand the number of known miRNAs in grass carp and are a valuable resource for better understanding the molecular biology of the spleen development in grass carp.

Mucus piRNAs profiles of Vibrio harveyi-infected Cynoglossus semilaevis: A hint for fish disease monitoring

The half-smooth tongue sole, Cynoglossus semilaevis, is an important cultured flatfish species. Vibrio harveyi is a common pathogen to this fish, which may result in great economic loss to C. semilaevis culture industry. piRNAs, a non-coding RNAs with 26-32 nt, have been regarded as promising biomarkers for cancer diagnosis and fish diseases. Here, we extracted the RNA from mucus of C. semilaevis and constructed the differential expression profiles of piRNAs between the sick fish (MS) and healthy fish (MC). We identified 45,696 differentially expressed piRNAs including 22,735 up-regulated piRNAs and 22,961 down-regulated piRNAs in MS group compared with MC group. The GO enrichment and KEGG pathway enrichment analyses of the differential piRNAs were carried out. The result showed immunity-related target genes mainly involved in immune system process, response to stimulus, cell killing, immune system, infectious diseases and cell growth and death. The 10 most differentially expressed piRNAs were chosen to perform the qRT-PCR, while only seven piRNAs were consistent with the sequence result. Compared with MC group, the expression levels of piR-mmu-72173>piR-rno-62831>piR-xtr-704880, piR-dme-15546979, piR-mmu-49941660, piR-mmu-29283297 and piR-mmu-1758399 were significantly lower, and piR-gga-10574 and piR-gga-134812 were significantly higher in MS group. These piRNAs may be potential biomarkers during the V. harveyi infection of C. semilaevis. This study could provide a new method to identify the infection status of C. semilaevis and understand better about the innate and adaptive immune system in C. semilaevis during bacterial infection.

Genome-wide integrated analysis reveals functions of lncRNA-miRNA-mRNA interactions in Atlantic salmon challenged by Aeromonas salmonicida

Aeromonas salmonicida (A. salmonicida) is a pathogenic bacterium that causes serious problems in the global Atlantic salmon aquaculture industry. In this study, we comprehensively analyzed the profiles of lncRNAs, miRNAs and mRNAs in gills of Atlantic salmon at high-dose A. salmonicida infection (3.06 × 10⁸ CFU/mL), low-dose A. salmonicida infection (3.06 × 10⁵ CFU/mL), and a PBS (100 μL) control. We identified 65 differentially expressed lncRNAs, 41 miRNAs, and 512 mRNAs between the control group and infection groups. Functional analysis showed that these genes were significantly enriched in the p53 signaling pathway, Wnt signaling pathway, mTOR signaling pathway, JAK-STAT signaling pathway, and Toll-like receptor signaling pathway. In addition, we predicted key genes in immune-related pathways and constructed a lncRNA-miRNA-mRNA network based on whole transcriptomic analysis. We further predicted three lncRNA-miRNA-mRNA axes as potential novel biomarkers in regulating the immune response of Atlantic salmon against A. salmonicida infection.

Immune role of the complement component 6 gene and its associated novel miRNA, miR-727, in half-smooth tongue sole (Cynoglossus semilaevis)

The complement component 6 (C6) gene is a component of the membrane attack complex (MAC), which causes rapid lytic destruction of bacteria. MicroRNAs (miRNAs) are small noncoding RNAs that regulate gene stability, including that of immune genes. However, current research on the function of C6 and its regulation by miRNAs is lacking. In the present study, we identified and characterized C6 and a novel miRNA, miR-727 (designated CsC6 and Cse-miR-727, respectively), of the half-smooth tongue sole (Cynoglossus semilaevis) that responded to infection with Vibrio anguillarum, a Gram-negative pathogen of marine fish. The full-length cDNA of CsC6 contained a 256 bp 5' untranslated region (5'-UTR), a 2820 bp open reading frame (ORF) encoding 939 amino acids, and a 205 bp 3'-UTR. SMART analysis showed that CsC6 contains typical C6 domains, including three TSP1 domains, one LDLa domain, one MACPF domain, two CCP domains and two FIMAC domains. CsC6 and Cse-miR-727 are widely expressed in the 13 tissues of half-smooth tongue sole, and their expression in immune tissues is significantly changed after V. anguillarum infection, generally showing an inverse trend. We confirmed that CsC6 was the target gene of Cse-miR-727 using the dual luciferase reporter assay and that Cse-miR-727 regulated CsC6 at the protein level using quantitative real-time polymerase chain reaction (qRT-PCR) and western blotting. The hepatic expression levels of not only the MAC components C7, C8α, C8β, C8γ and C9 but also the MAPKs, NF-κβ, AP-1, IL1β, IL6 and TNFα, which are involved in many signaling pathways, changed significantly in half-smooth tongue sole following stimulation with the Cse-miR-727 agomir and inhibitor. This evidence suggested that CsC6 could be mediated by Cse-miR-727 to affect MAC assembly and immune signaling pathways in half-smooth tongue soles. To our best knowledge, this study is the first to investigate the regulatory mechanism and immune response of complement genes mediated by miRNAs in fish.

PUF60 of Japanese flounder is regulated by pol-miR-novel_395 and involved in pathogen infection, autophagy, and apoptosis

MicroRNAs (miRNAs) are evolutionary conserved, non-coding small RNAs that have been shown to regulate diverse biological processes including immunity. In a previous study, a novel miRNA of Japanese flounder (Paralichthys olivaceus), pol-miR-novel_395, was found to be responsive in expression to the infection of the bacterial pathogen Edwardsiella tarda. In the present study, we examined the regulation and immune effect of pol-miR-novel_395 and its target gene. We found that pol-miR-novel_395 expression was regulated by E. tarda and megalocytivirus, and pol-miR-novel_395 targeted the gene of PUF60 (poly (U)-binding-splicing factor 60 kDa) of flounder (named PoPUF60). Constitutive expression of PoPUF60 occurred in relatively high levels in the heart and liver of flounder. Bacterial infection upregulated PoPUF60 expression, whereas viral infection downregulated PoPUF60 expression. Interference with PoPUF60 expression or overexpression of pol-miR-novel_395 in flounder cells strongly potentiated E. tarda infection. Consistently, in vivo knockdown of PoPUF60 enhanced bacterial dissemination in the tissues of flounder but blocked viral replication, whereas in vivo overexpression of PoPUF60 inhibited bacterial dissemination but facilitated viral replication. Additionally, pol-miR-novel_395 and PoPUF60 were involved in the process of autophagy and apoptosis. Collectively, these results indicated that PoPUF60 and pol-miR-novel_395 play an important role in pathogen infection, autophagy, and apoptosis.

Comparative Mucous miRomics in Cynoglossus semilaevis Related to Vibrio harveyi Caused Infection

Epidermal mucus is an important barrier and regulating mediator in fish. MicroRNAs (miRNAs) are proved to be involved in various biological processes, also as promising biomarkers for disease diagnosis. Vibrio harveyi has long been a noticeable bacterial pathogen in Cynoglossus semilaevis aquaculture. To find the evidence whether there are indicating miRNAs in mucus and whether the miRNAs are related to infections caused by V. harveyi, miRNA profiles of mucus from V. harveyi infected fish and healthy controls were screened by small RNA sequencing and verified by quantitative real-time PCR. This is the first report about miRNA profiling of flatfish mucus, aiming at illustrating the pathogenesis of V. harveyi caused infection and developing disease-related biomarkers. The results revealed significant differences in expression levels of some miRNAs between infected fish and healthy ones. Three hundred differentially expressed miRNAs were obtained after filtering through FC > 2 or FC < 0.5 and most of the differential miRNAs were downregulated. After verification through qRT-PCR, four unique miRNAs, dre-miR-451, dre-miR-184, dre-miR-205-5p > ssa-miR-205b-5p, and dre-miR-181a-5p > ssa-miR-181a-5p, were identified as V. harveyi infection-related signatures, consistent with sequencing trend. The expression levels of these four miRNAs in the infected fish were all significantly lower than controls. These miRNAs in mucus could be used to differentiate diseased and healthy fish in a non-invasive way with practical value for large-scale disease screening. They also provided new insights into the mechanism underlying the bacterial infections in fish.

Expression of ssa-miR-155 during ISAV infection in vitro: Putative role as a modulator of the immune response in Salmo salar

Multiple cellular components are involved in pathogen-host interaction during viral infection; in this context, the role of miRNAs have become highly relevant. We assessed the expression of selected miRNAs during an in vitro infection of a Salmo salar cell line with Infectious Salmon Anemia Virus (ISAV), the causative agent of a severe disease by the same name. Salmon orthologs for miRNAs that regulate antiviral responses were measured using RT-qPCR in an in vitro time-course assay. We observed a modulation of specific miRNAs expression, where ssa-miR-155-5p was differentially over-expressed. Using in silico analysis, we identified the putative mRNA targets for ssa-miR-155-5p, finding a high prevalence of hosts immune response-related genes; moreover, several mRNAs involved in the viral infective process were also identified as targets for this miRNA. Our results suggest a relevant role for miR-155-5p in Salmo salar during an ISAV infection as a regulator of the immune response to the virus.

Proteomics of mucosal exosomes of Cynoglossus semilaevis altered when infected by Vibrio harveyi

The cargo of exosomes contains proteins with various functions, which might be promising biomarkers for disease diagnosis and prognosis. To explore the impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis from a different perspective and develop promising biomarkers for infection, the exosomes from epidermal mucus of healthy controls（EC）and sick fish（ES）were extracted and identified, coupled with proteomic screening through iTRAQ followed with LC-MS/MS. 1531 credible proteins were obtained relating to structural, metabolic and immunological functions. 359 different expressed proteins (DEPs) (FC > 2 or FC < 0.5) were found, with 161 up-regulated and 198 down-regulated in ES. Based on the database of C. semilaevis on Uniprot, 71 proteins were characterized as concrete names, including 19 up-regulated proteins and 52 down-regulated proteins, and were selected as subjects for further studies. Ferritin, Toll-like receptor 5S protein and Calcium-transporting ATPase were upregulated, while Histone H2B and Eukaryotic translation initiation factor 5A were downregulated, consistent with the expression levels of related mRNAs in skin tissue verified by qRT-PCR. The integrated analysis between miRomics and proteomics also provided possible regulatory relationships mediated by mucous exosomes during infection. The signature proteins in mucosal exosomes could make sense in the explanation of the infection defending mechanism and the development of biomarkers which can differentiate diseased and healthy C. semilaevis individuals.

The emerging regulatory roles of noncoding RNAs in immune function of fish: MicroRNAs versus long noncoding RNAs

The genome could be considered as raw data expressed in proteins and various types of noncoding RNAs (ncRNAs). However, a large portion of the genome is dedicated to ncRNAs, which in turn represent a considerable amount of the transcriptome. ncRNAs are modulated on levels of type and amount whenever any physiological process occurs or as a response to external modulators. ncRNAs, typically forming complexes with other partners, are key molecules that influence diverse cellular processes. Based on the knowledge of mammalian biology, ncRNAs are known to regulate and control diverse trafficking pathways and cellular activities. Long noncoding RNAs (lncRNAs) notably have diverse and more regulatory roles than microRNAs. Expanding these studies on fish has derived the same conclusion with relevance to other species, including invertebrates, explored the potentials to harness such types of RNA to further understand the biology of such organisms, and opened gates for applying recent technologies, such as RNA interference and delivering micromolecules as microRNAs to living cells and possibly to target organs. These technologies should improve aquaculture productivity and fish health, as well as help understand fish biology.

Extracellular vesicles piwi-interacting RNAs from skin mucus for identification of infected Cynoglossus semilaevis with Vibrio harveyi

Extracellular vesicles play a regulatory role in intracellular and intercellular transmission through a variety of biological information molecules, including mRNA, small RNAs and proteins. piRNAs are one kind of regulatory small RNAs in the vesicles at the post transcriptional level. Hereby, we isolated the extracellular vesicles from skin mucus and screened the piRNA profiles of these vesicles, aiming at developing biomarkers related to bacterial infections in Cynoglossus semilaevis. The different profilings of piRNAs in mucous extracellular vesicles of C. semilaevis were compared through small RNA sequencing, between fish infected with Vibrio harveyi and healthy ones. The number of clean reads on the alignment of exosome sick (ES) group was 105, 345 and that of exosome control (EC) group was 455, 144. GO and KEGG pathway enrichment analysis showed that most of the target genes were involved in cellular process, response to stimulus, biological regulation, immune system process and signal transduction, signal molecular and interaction, transport and catabolism. The 45 final candidate piRNAs related to immunity or infectious diseases included 20 piRNAs with high expression in the ES group and 25 piRNAs with a low expression in the ES group. After verification by qRT-PCR, there was significant difference of five piRNAs expression level between infected fish and healthy fish, in line with the sequencing. The expression level of piR-mmu-16401212, piR-mmu-26829319 and piR-gga-244092 in infected fish were significantly lower than that of control group, while piR-gga-71717 and piR-gga-99034 were higher, which implying that these piRNAs in mucous extracellular vesicles can be used to identify diseased fish from normal ones. This work supplied a novel class of biomarker for infection diagnosis in fish, and it will be benefit for screening disease resistant breeding of C. semilaevis.

Detecting Cynoglossus semilaevis infected with Vibrio harveyi using micro RNAs from mucous exosomes

Exosomes are important mediators of vesicle transportation and contain microRNAs (miRNAs) that mediate transcriptional gene knockout and silencing in biological processes. Moreover, exosomic miRNAs are promising biomarkers for disease diagnosis and physiological status indication in many species, including fish. The impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis aquaculture is becoming more and more serious as the industry expands. To overcome this challenge, miRNAs in mucous exosomes were screened by small RNA sequencing and verified by quantitative real-time PCR to develop biomarkers. This is the first capture of exosomes from flatfish mucus coupled with miRNA profiling. The results revealed significant differences in expression levels of some miRNAs between infected and healthy fish. Three unique miRNAs were identified for V. harveyi infection diagnosis; expression levels of dre-miR-205-5p and dre-miR-205-5p in infected fish were significantly lower than controls, while dre-miR-100-5p expression was higher. These miRNAs in mucous exosomes could be used to differentiate diseased and healthy fish in an early screening method with practical value for breeding disease-resistant C. semilaevis.

Micro-Transcriptome Analysis Reveals Immune-Related MicroRNA Regulatory Networks of Paralichthys olivaceus Induced by Vibrio anguillarum Infection

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.

The novel fish miRNA pol-miR-novel_171 and its target gene FAM49B play a critical role in apoptosis and bacterial infection

MicroRNAs (miRNAs) are a type of small, non-coding RNAs that participate in many cellular and biological processes by regulating mRNA stability. In a previous study, we identified 96 Japanese flounder (Paralichthys olivaceus) miRNAs responsive to the infection of Edwardsiella tarda, a bacterial pathogen to fish as well as humans. In the current study, we examined the regulation and function of one novel miRNA, i.e., pol-miR-novel_171, from the above 96 miRNA pool. We found that pol-miR-novel_171 expression was regulated by E. tarda and megalocytivirus in a pathogen-specific manner, and that pol-miR-novel_171 targeted the gene of FAM49B (family with sequence similarity 49 member B) of flounder (named PoFAM49B) by negative interaction with the 3'-UTR of PoFAM49B. To date, the function fish FAM49B is unknown. We found that PoFAM49B expressed in multiple tissues of flounder, and recombinant PoFAM49B interacted with and inhibited the growth of Gram-negative bacterial pathogens. Interference with PoFAM49B expression in flounder cells promoted E. tarda infection. Similar effects on E. tarda infection were observed with pol-miR-novel_171 overexpression. Consistently, in vivo knockdown of PoFAM49B in flounder enhanced E. tarda dissemination in fish tissues. Furthermore, interference with PoFAM49B expression, or overexpression of pol-miR-novel_171, promoted apoptosis of flounder cells, while in vitro and in vivo knockdown of PoFAM49B augmented the expressions of key apoptosis-associated genes. These results revealed for the first time the immune function of fish FAM49B and the regulatory mechanism of a novel fish miRNA by demonstrating that pol-miR-novel_171, via PoFAM49B, played a critical role in apoptosis and anti-bacterial immunity.

Characterization of Streptococcus iniae-induced microRNA profiles in Paralichthys olivaceus and identification of pol-3p-10740_175 as a regulator of antibacterial immune response

MicroRNAs (miRNAs) are involved in many biological activities including immune defense against pathogens. In this study, we applied high-throughput sequencing technology to examine miRNAs in Japanese flounder (Paralichthys olivaceus) infected with Streptococcus iniae at different times. A total of 1038 miRNAs were identified, of which, 249 were novel miRNAs, and 81 showed differential expression (named DEmiRNAs) after S. iniae infection. Of the 81 DEmiRNAs identified, 34 and 58 occurred at 6 h and 24 h post-infection, respectively; most DEmiRNAs were strongly time-specific, and only 13.6% of the DEmiRNAs were shared between the two time points. A total of 9582 target genes were predicted for the 81 DEmiRNAs. The putative target genes were enriched in various GO and KEGG pathways of biological processes and molecular/cellular functions, in particular endocytosis, regulation of transcription, lysososme, and the signaling pathways of MAPK, ErbB, and AMPK. One of the DEmiRNAs, pol-3p-10740_175, was found to target dual specificity phosphatase 6 (Dusp6) and repress the expression of the latter. Transfection of flounder FG cells with pol-3p-10740_175 caused a significant inhibition on S. iniae invasion. The results of this study provided the first S. iniae-induced miRNA profile in Japanese flounder and indicated that flounder miRNAs play an important role in antibacterial immunity.

Effects of microRNA-731 on inflammation and apoptosis by targeting CiGadd45aa in grass carp

Ctenopharyngodon idellagrowth arrest and DNA damage-inducible 45aa (CiGadd45aa) is a member of the Gadd45 family of immune-related proteins in grass carp. There is increasing evidence that microRNAs (miRNAs) are involved in the regulation of inflammatory responses and apoptosis. However, little is known about the regulatory effects of miRNAs on CiGadd45aa expression in grass carp. In the present study, CiGadd45aa was identified as a target gene of miR-731 based on miRNA expression profiling and dual-luciferase reporter assay. Our study revealed that miR-731 targets CiGadd45aa and regulates the expression of proinflammatory factors, thereby regulating immune response in grass carp. In addition, CiGadd45aa and miR-731 were both found induced apoptosis. Hence, our findings provide a theoretical basis for exploring the molecular mechanism by which miR-731 regulates inflammation and apoptosis in grass carp.

Dietary Immunostimulant CpG Modulates MicroRNA Biomarkers Associated with Immune Responses in Atlantic Salmon (Salmo salar)

MicroRNAs (miRNAs) are key regulators in fish immune responses. However, no study has previously characterized the impact of polyriboinosinic polyribocytidylic acid (pIC) and formalin-killed typical Aeromonas salmonicida (ASAL) on miRNA expression in Atlantic salmon fed a commercial diet with and without immunostimulant CpG. To this end, first, we performed small RNA deep sequencing and qPCR analyses to identify and confirm pIC- and/or ASAL-responsive miRNAs in the head kidney of salmon fed a control diet. DESeq2 analyses identified 12 and 18 miRNAs differentially expressed in pIC and ASAL groups, respectively, compared to the controls. Fifteen of these miRNAs were studied by qPCR; nine remained significant by qPCR. Five miRNAs (miR-27d-1-2-5p, miR-29b-2-5p, miR-146a-5p, miR-146a-1-2-3p, miR-221-5p) were shown by qPCR to be significantly induced by both pIC and ASAL. Second, the effect of CpG-containing functional feed on miRNA expression was investigated by qPCR. In pre-injection samples, 6 of 15 miRNAs (e.g., miR-181a-5-3p, miR-462a-3p, miR-722-3p) had significantly lower expression in fish fed CpG diet than control diet. In contrast, several miRNAs (e.g., miR-146a-1-2-3p, miR-192a-5p, miR-194a-5p) in the PBS- and ASAL-injected groups had significantly higher expression in CpG-fed fish. Multivariate statistical analyses confirmed that the CpG diet had a greater impact on miRNA expression in ASAL-injected compared with pIC-injected fish. This study identified immune-relevant miRNA biomarkers that will be valuable in the development of diets to combat infectious diseases of salmon.

Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

The sustainable development of aquaculture has been impeded by infectious diseases worldwide. However, the genomic architecture and the genetic basis underlying the disease resistance remain poorly understood, which severely hampers both the understanding of the evolution of fish disease resistance traits and the prevention of these diseases in the aquaculture community. Cynoglossus semilaevis is a representative and commercially-important flatfish species. Here we combined genome-wide association study and Fst and nucleotide diversity filtration to identify loci important for the disease resistance. Based on 1,016,774 single-nucleotide polymorphisms (SNPs) identified from 650 Gb genome resequencing data of 505 individuals, we detected 33 SNPs significantly associated with disease resistance and 79 candidate regions after filtration steps. Both the allele frequencies and genotype frequencies of the associated loci were significantly different between the resistant and susceptible fish, suggesting a role in the genetic basis of disease resistance. The SNP with strongest association with disease resistance was located in Chr 17, at 145 bp upstream of fblx19 gene, and overlapped with the major quantitative trait locus previously identified. Several genes, such as plekha7, nucb2, and fgfr2, were also identified to potentially play roles in the disease resistance. Furthermore, the expression of some associating genes were likely under epigenetic regulations between the bacterial resistant and susceptible families. These results provide insights into the mechanism that enable variation of disease resistance to bacterial pathogen infection. The identified polymorphisms and genes are valuable targets and molecular resources for disease resistance and other traits, and for advanced breeding practice for superior germplasm in fish aquaculture.

Global profiling and characterization of Japanese flounder (Paralichthys olivaceus) kidney microRNAs regulated by Edwardsiella tarda infection in a time-dependent fashion

Japanese flounder (Paralichthys olivaceus) is an important economic fish species farmed in China and other countries. It is susceptible to infection by Edwardsiella tarda, a severe fish pathogen with a broad host range. In this study, we employed high-throughput deep sequencing technology to identify, in a global scale, flounder kidney microRNAs (miRNAs) induced by E. tarda at different stages of infection. Differentially expressed miRNAs (DEmiRNAs) and mRNAs (DEmRNAs) exhibiting significantly altered expression levels before and after E. tarda infection were examined. A total of 96 DEmiRNAs were identified, for which 2779 target genes were predicted. Eighty-seven miRNA-mRNA pairs, involving 29 DEmiRNAs and 86 DEmRNAs, showed negative correlations in their expression patterns. GO and KEGG enrichment analysis revealed that the putative target genes of the DEmiRNAs were associated with diverse biological processes, cellular components, and molecular functions. One of the DEmiRNAs, pol-miR-182-5p, was demonstrated to regulate sphingosine-1-phosphate receptor 1 (PoS1PR1) negatively in a manner that depended on the specific interaction between the seed sequence of pol-miR-182-5p and the 3'-UTR of PoS1PR1. Overexpression of pol-miR-182-5p in flounder cells promoted apoptosis and inhibited cellular viability. Knockdown of PoS1PR1 in flounder enhanced E. tarda invasion and dissemination in fish tissues. These results provide new insights into miRNA-mediated anti-bacterial immunity in flounder.

A transcriptome analysis focusing on splenic immune-related mciroRNAs of rainbow trout upon Aeromonas salmonicida subsp. salmonicida infection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that can regulate the immune responses during pathogen infection. Aeromonas salmonicida (A. salmonicida) subsp. salmonicida is the causative agent of furunculosis in salmon and trout. To identify the miRNAs and investigate the specific miRNAs in rainbow trout upon A. salmonicida subsp. salmonicida infection, we performed high throughput sequencing using the spleens of rainbow trout infected with and without an A. salmonicida subsp. salmonicida clinical isolate. A total of 381 known miRNAs and 926 novel miRNAs were identified. Eleven known and 16 novel miRNAs were found to be differentially expressed upon infection. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses indicated that the target genes of the differentially expressed miRNAs were closely associated with immune responses and biological regulations. Additionally, over- and suppressed expression of miR-155-5p significantly enhanced and reduced the IL-2 and IL-1β expressions in RTG-2 cells induced by A. salmonicida, respectively. To our knowledge, this is the first experimental study on the miRNAs of rainbow trout upon A. salmonicida infection. The results here might lay a foundation for the further understanding of the roles of miRNAs in the immune responses during A. salmonicida infection in rainbow trout.

Dynamics of MiRNA Transcriptome in Turbot (Scophthalmus maximus L.) Intestine Following Vibrio anguillarum Infection

MicroRNAs (miRNAs) are a group of small non-coding RNAs, which could bind to the 3'-untranslated regions of their target mRNAs to regulate gene expression in various biological processes, including immune-regulated signaling pathways. Turbot (Scophthalmus maximus L.), an important commercial fish species in China, has been suffering with Vibrio anguillarum infection resulted in dramatic economic loss. Therefore, we investigated the expression profiles of miRNAs, as well as the immune-related miRNA-mRNA pairs in turbot intestine at 1 h, 4 h, and 12 h following V. anguillarum infection. As a result, 266 predicted novel miRNAs and 283 conserved miRNAs belonging to 92 miRNA families were detected. A total of 44 miRNAs were differentially expressed in the intestine following V. anguillarum infection. Following prediction, the potential target genes of differentially expressed miRNAs were grouped into a wide range of functional categories, including immune defense/evasion, inflammatory responses, RIG-I signaling pathway, and Toll-like receptor signaling pathway. Moreover, we selected 15 differentially expressed immune genes and their related differentially expressed miRNAs to construct an interaction network for V. anguillarum infection in turbot. These results suggested that in teleost, as in higher vertebrates, miRNAs prominently contribute to immune responses, protecting the host against infection. In addition, this is the first report of comprehensive identification of turbot miRNAs being differentially regulated in the intestine related to V. anguillarum infection. Our results provided an opportunity for further understanding of the molecular mechanisms of miRNA regulation in turbot host-pathogen interactions.

Catching the complexity of salmon-louse interactions

The study of host-parasite relationships is an integral part of the immunology of aquatic species, where the complexity of both organisms has to be overlayed with the lifecycle stages of the parasite and immunological status of the host. A deep understanding of how the parasite survives in its host and how they display molecular mechanisms to face the immune system can be applied for novel parasite control strategies. This review highlights current knowledge about salmon and sea louse, two key aquatic animals for aquaculture research worldwide. With the aim to catch the complexity of the salmon-louse interactions, molecular information gleaned through genomic studies are presented. The host recognition system and the chemosensory receptors found in sea lice reveal complex molecular components, that in turn, can be disrupted through specific molecules such as non-coding RNAs.

Edwardsiella tarda-induced miR-7a functions as a suppressor in PI3K/AKT/GSK3β signaling pathway by targeting insulin receptor substrate-2 (IRS2a and IRS2b) in Paralichthys olivaceus

To study the effect of Edwardsiella tarda infection on miRNAs expression profile in Japanese flounder, fish were injected intraperitoneally with E. tarda. The miRNAs involved in regulating immune responses were analyzed by high-throughput sequencing. A total of 164 mature miRNAs were identified, of which 17 miRNAs were differentially expressed (DE miRNAs) after E. tarda infection, indicating that they were immune-related miRNAs. To further examine the relationship between the miRNAs and their predicted target mRNAs, a total of 22 predicted target mRNAs, mainly related to endocytic signaling pathway, NF-κB signaling pathway, and p53 signaling pathway, were detected with miRNA mimics in HEK-293T cells by dual-luciferase reporter experiments. Finally, we confirmed that insulin receptor substrate-2 (IRS2a and IRS2b) were regulated by miR-7a. And the target sites of the 3' untranslated region (UTR) of IRS2a and IRS2b were verified by dual-luciferase reporter experiments. Furthermore, we found that the E. tarda and LPS significantly increased host miR-7a expression. In vivo and in vitro studies revealed that IRS2-mediated PI3K/AKT/GSK3β signaling pathway was suppressed. Taken together, these results implied that miR-7a might be a key regulator of PI3K/AKT/GSK3β signaling pathway via suppressing the IRS2a and IRS2b genes.

Identification and expression profiling analysis of microRNAs in Nile tilapia (Oreochromis niloticus) in response to Streptococcus agalactiae infection

MicroRNAs (miRNAs) play vital regulatory roles in various biological processes, including in immune responses. Nile tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify immune-related miRNAs of O. niloticus, 4 libraries from liver during S. agalactiae infection (0 h, 5 h, 50 h, and 7 d) were sequenced by high-throughput sequencing technology in tilapia. We obtained 10,703,531, 11,507,163, 11,180,179 and 13,408,414 clean reads per library, respectively. In our results, a total of 482 miRNAs were identified through bioinformatic analysis, including 220 conserved miRNAs and 262 putative novel miRNAs. Moreover, 21 (4.36%), 50 (10.37%), and 46 (9.54%) miRNAs were significantly differentially expressed at 5 h, 50 h and 7 d, respectively. In addition, 6939 target genes regulated by these differentially expressed miRNAs were predicted, and their functional annotations were predicted by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, which revealed that a majority of differentially expressed miRNAs were involved in apoptotic process, metabolic process, and immune responses. Finally, Real-time quantitative PCR experiments were performed for 7 miRNAs by stem-loop RT-PCR, and a general agreement was confirmed between the sequencing and RT-qPCR data. To our understanding, this is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in liver related to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions, and genetic resources for molecular assistant selection for disease resistant breeding program.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

Non-coding RNAs (ncRNAs) are functional RNA molecules that are transcribed from DNA but not translated into proteins. ncRNAs function as key regulators of gene expression and chromatin modification. Recently, the functional role of ncRNAs in teleost fish has been extensively studied. Teleost fish are a highly diverse group among the vertebrate lineage. Fish are also important in terms of aquatic ecosystem, food production and human life, being the source of animal proteins worldwide and models of biomedical research. However, teleost fish always suffer from the invasion of infectious pathogens including viruses and bacteria, which has resulted in a tremendous economic loss to the fishing industry worldwide. Emerging evidence suggests that ncRNAs, especially miRNAs and lncRNAs, may serve as important regulators in cytokine and chemokine signaling, antigen presentation, complement and coagulation cascades, and T cell response in teleost fish. In this review, we summarize current knowledge and understanding of the roles of both miRNAs and lncRNAs in immune regulation in teleost fish. Molecular mechanism insights into the function of ncRNAs in fish immune response may contribute to the development of potential biomarkers and therapeutic targets for the prevention and treatment of fish diseases.

Identification of microRNAs with heat stress responsive and immune properties in Marsupenaeus japonicus based on next-generation sequencing and bioinformatics analysis: Essential regulators in the heat stress-host interactions

Summer mortality syndrome is one of the most serious issue for Marsupenaeus japonicus aquaculture in China. Since it causes massive economic loss and threatens sustainability of M. japonicus aquaculture industry, thus, there is an urgent desire to reveal the heat stress-host interactions mechanisms that lead to mass mortalities of M. japonicus in hot summer months. MicroRNAs (miRNAs) are small noncoding RNAs that involved in regulation of diverse biological processes, including stress and immune response, and might serve as potential regulators in the heat stress-host interactions. In the present study, miRNAs with heat stress responsive and immune properties were identified and characterized in M. japonicus by small RNA sequencing and bioinformatics analysis. In total, 79 host miRNAs were identified, among which 15 miRNAs were differentially expressed in response to heat stress. Target genes prediction and function annotation revealed that a variety of host cellular processes, such as signal transduction, transcription, anti-stress response, ribosomal biogenesis, lipid metabolism, cytoskeleton, etc, were potentially subject to miRNA-mediated regulation in response to heat stress. Furthermore, a total of 30 host miRNAs that potentially involved in interaction with white spot syndrome virus (WSSV) were obtained via predicting and analyzing the target genes from WSSV. The results showed that a batch of WSSV genes that code for structural proteins and enzymes that are essential for WSSV infection and proliferation, such as envelope proteins, capsid proteins, immediate-early proteins, collagen-like protein, protein kinase, thymidylate synthetase, TATA-box bind protein, etc, were predicted to be targeted by host miRNAs. Several of the host miRNAs with predicted antiviral capacity were down-regulated under heat stress, indicating a repression of host miRNA-mediated antiviral immune response. This study highlighted the essential roles of host miRNAs in the heat stress-host interactions and provided valuable information for further investigation on the mechanism of miRNA-mediated heat stress and immune response of shrimp.

Under control: The innate immunity of fish from the inhibitors' perspective

The innate immune response involves a concerted network of induced gene products, preformed immune effectors, biochemical signalling cascades and specialised cells. However, the multifaceted activation of these defensive measures can derail or overshoot and, if left unchecked, overwhelm the host. A plenty of regulatory devices therefore mediate the fragile equilibrium between pathogen defence and pathophysiological manifestations. Over the past decade in particular, an almost complete set of teleostean sequences orthologous to mammalian immunoregulatory factors has been identified in various fish species, which prove the remarkable conservation of innate immune-control concepts among vertebrates. This review will present the current knowledge on more than 50 teleostean regulatory factors (plus additional fish-specific paralogs) that are of paramount importance for controlling the clotting cascade, the complement system, pattern-recognition pathways and cytokine-signalling networks. A special focus lies on those immunoregulatory features that have emerged as potential biomarker genes in transcriptome-wide research studies. Moreover, we report on the latest progress in elucidating control elements that act directly with immune-gene-encoding nucleic acids, such as transcription factors, hormone receptors and micro- and long noncoding RNAs. Investigations into the function of teleostean inhibitory factors are still mainly based on gene-expression profiling or overexpression studies. However, in support of structural and in-vitro analyses, evidence from in-vivo trials is also available and revealed many biochemical details on piscine immune regulation. The presence of multiple gene copies in fish adds a degree of complexity, as it is so far hardly understood if they might play distinct roles during inflammation. The present review addresses this and other open questions that should be tackled by fish immunologists in future.

Acute Hepatopancreatic Necrosis Disease (AHPND) related microRNAs in Litopenaeus vannamei infected with AHPND-causing strain of Vibrio parahemolyticus

BACKGROUND

Acute hepatopancreatic necrosis disease (AHPND) has emerged as a major debilitating disease that causes massive shrimp death resulting in substantial economic losses in shrimp aquaculture. Given that several diseases and infections have been associated with microRNAs (miRNAs), we conducted a comparative transcriptomic analysis using the AHPND (VA) and non-AHPND (VN) strains of Vibrio parahemolyticus to identify miRNAs potentially involved in AHPND pathogenesis in Litopenaeus vannamei.

RESULTS

A total of 83 miRNAs (47 upregulated and 36 downregulated) were significantly differentially expressed between the VA and VN challenged groups, while 222 target genes of these miRNAs were predicted. Functional enrichment analysis revealed that the miRNAs target genes were involved in multiple biological processes including metabolic pathways, amoebiasis, Vibrio cholerae infection etc. Finally, interaction network and qPCR (Real-time Quantitative PCR) analysis of 12 potential key AHPND-related miRNAs and their predicted target genes, revealed their possible involvement in modulating several immune-related processes in the pathogenesis of AHPND.

CONCLUSIONS

We have shown using comparative transcriptomic analysis, miRNAs and their target genes that are responsive to AHPND V. parahemolyticus infection in shrimp, therefore suggesting their possible role in defense response to AHPND V. parahemolyticus infection.

MicroRNA-based transcriptomic responses of Atlantic salmon during infection by the intracellular bacterium Piscirickettsia salmonis

MicroRNAs (miRNAs) are small non-coding RNAs that have emerged as key regulators in diverse biological processes across taxa. However, despite the importance of these transcripts, little is known about their role during the immune response in salmonids. Because of this, we use deep sequencing technologies to explore the microRNA-based transcriptomic response of the Atlantic salmon (Salmo salar) to the intracellular bacteria Piscirickettsia salmonis, one of the main threats to salmon aquaculture in Chile. Hence, 594 different miRNAs were identified from head kidney and spleen transcriptomic data. Among them, miRNA families mir-181, mir-143 and mir-21 were the most abundant in control groups, while after infection with P. salmonis, mir-21, mir-181 and mir-30 were the most predominant families. Furthermore, transcriptional analysis revealed 84 and 25 differentially expressed miRNAs in head kidney and spleen respectively, with an overlapping response of 10 miRNAs between the analyzed tissues. Target prediction, coupled with GO enrichment analysis, revealed that the possible targets of the most regulated miRNAs were genes involved in the immune response, such as cortisol metabolism, chemokine-mediated signaling pathway and neutrophil chemotaxis genes. Among these, predicted putative target genes such as C-C motif chemokine 19-like, stromal cell-derived factor 1-like, myxovirus resistance protein 2 and hepcidin-1 were identified. Overall, our results suggest that miRNA expression in co-modulation with transcription activity of target genes is related to putative roles of non-coding RNAs in the immune response of Atlantic salmon against intracellular bacterial pathogens.

Modulation of Atlantic salmon miRNome response to sea louse infestation

MicroRNAs are non-coding RNA that plays a crucial role in post-transcriptional regulation and immune system regulation. On other hand, sea lice are prevalent parasites that affect salmon farming, generating different degrees of immune suppression depending on the salmon and sea louse species. Caligus rogercresseyi for example, which affects the salmon industry in Chile, decreases Th1 response, macrophage activation, TLR-mediated response and iron regulation in infected fish. In this study, we explore Atlantic salmon miRNome during infestation by C. rogercresseyi. Using small RNA sequencing, we annotated 1718 miRNAs for skin and head kidney from infected Atlantic salmon. The most abundant families identified were mir-10, mir-21, mir-30, mir-181 and let7. Significant differences were found between tissue, with 1404 annotated miRNA in head kidney and 529 in skin. Differential analysis of transcript expression indicated that at an early stage of infestation miRNA expression was higher in head kidney than in skin tissue, revealing tissue-specific expression patterns. In parallel, miRNA target prediction using 3'UTRs from highly regulated immune-related genes and iron metabolism showed that mir-140-4 and mir-181a-2-5 modulate the expression of TLR22 and Aminolevulinic acid synthase, respectively. This study contributes knowledge about the immune response of Atlantic salmon during infestation with sea lice.

miRNAs associated with immune response in teleost fish

MicroRNAs (miRNAs) have been identified as important post transcriptional regulators of gene expression. In higher vertebrates, a subset of miRNAs has been identified as important regulators of a number of key genes in immune system gene networks, and this paper review recent studies on miRNAs associated with immune response in teleost fish. Challenge studies conducted in several species have identified differently expressed miRNAs associated with viral or bacterial infection. The results from these studies point out several miRNAs that are likely to have evolutionary conserved functions that are related to immune response in teleost fish. Changed expression levels of mature miRNAs from the five miRNA genes miRNA-462, miRNA-731, miRNA-146, miRNA-181 and miRNA-223 are observed following viral as well as bacterial infection in several teleost fish. Furthermore, significant changes in expression of mature miRNAs from the five genes miRNA-21, miRNA-155, miRNA-1388, miRNA-99 and miRNA-100 are observed in multiple studies of virus infected fish while changes in expression of mature miRNA from the three genes miRNA-122, miRNA-192 and miRNA-451 are observed in several studies of fish with bacterial infections. Interestingly, some of these genes are not present in higher vertebrates. The function of the evolutionary conserved miRNAs responding to infection depends on the target gene(s) they regulate. A few target genes have been identified while a large number of target genes have been predicted by in silico analysis. The results suggest that many of the targets are genes from the host's immune response gene networks. We propose a model with expected temporal changes in miRNA expression if they target immune response activators/effector genes or immune response inhibitors, respectively. The best way to understand the function of a miRNA is to identify its target gene(s), but as the amount of genome resources for teleost fish is limited, with less well characterized genomes and transcriptomes, identifying the true target genes of the miRNAs associated with the immune response is a challenge. Identifying such target genes by applying new methods and approaches will likely be the next important step to understand the function of the miRNAs associated with immune response in teleost fish.

Characterization of conserved and novel miRNAs using deep sequencing and prediction of miRNA targets in Crucian carp (Carassius auratus)

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs of -22 nucleotides that can base pair with their target mRNAs, which represses their translation or induces their degradation in various biological processes. However, little is known about identification of miRNAs and their target genes in C. auratus. In the present study, a small RNA library from pooled tissue of C. auratus was constructed and sequenced using the deep sequencing. A total of 320 conserved miRNAs (belonging to 105 families) as well as 11 potentially novel miRNAs were identified. Stem-loop qRT-PCR analysis confirmed that both conserved and novel miRNAs were expressed in C. auratus, and some of them were preferentially expressed in certain tissues. Subsequently, a total of 1668 potential target genes were predicted for these identified miRNAs and GO analysis showed that most of the targets were involved in lots of physiological actions. This study represents a first large-scale identification and characterization of C. auratus miRNAs and their potential target genes. Taken together, our results add new information to existing data on C. auratus miRNAs and should be useful for investigating the biological functions of miRNAs in fishes and other aquatic species.

Transcriptome-wide analysis of microRNAs in Branchiostoma belcheri upon Vibrio parahemolyticus infection

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that participate in diverse biological processes via regulating expressions of target genes at post-transcriptional level. Amphioxus, as modern survivor of an ancient chordate lineage, is a model organism for comparative genomics study. However, miRNAs involved in regulating immune responses in Branchiostoma belcheri are largely unclear. Here, we systematically investigated the microRNAs (miRNAs) involved in regulating immune responses in the cephalochordate amphioxus (Branchiostoma belcheri) through next-generation deep sequencing of amphioxus samples infected with Vibrio parahemolyticus. We identified 198 novel amphioxus miRNAs, consisting of 12 conserved miRNAs, 33 candidate star miRNAs and 153 potential amphioxus-specific-miRNAs. Using microarray profiling, 14 miRNAs were differentially expressed post infection, suggesting they are immune-related miRNAs. Eight miRNAs (bbe-miR-92a-3p, bbe-miR-92c-3p, bbe-miR-210-5p, bbe-miR-22-3p, bbe-miR-1∼bbe-miR-133 and bbe-miR-217∼bbe-miR-216 clusters) were significantly increased at 12 h post-infection, while bbe-miR-2072-5p was downregulated at 6 h and 12 h. Three miRNAs, bbe-miR-1-3p, bbe-miR-22-3p and bbe-miR-92a-3p, were confirmed to be involved in immune responses to infection by qRT-PCR. Our findings further clarify important regulatory roles of miRNAs in the innate immune response to bacterial infection in amphioxus.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Identification and Characterization of MicroRNAs in the Liver of Blunt Snout Bream (Megalobrama amblycephala) Infected by Aeromonas hydrophila

MicroRNAs (miRNAs) are small RNA molecules that play key roles in regulation of various biological processes. In order to better understand the biological significance of miRNAs in the context of Aeromonas hydrophila infection in Megalobrama amblycephala, small RNA libraries obtained from fish liver at 0 (non-infection), 4, and 24 h post infection (poi) were sequenced using Illumina deep sequencing technology. A total of 11,244,207, 9,212,958, and 7,939,157 clean reads were obtained from these three RNA libraries, respectively. Bioinformatics analysis identified 171 conserved miRNAs and 62 putative novel miRNAs. The existence of ten randomly selected novel miRNAs was validated by RT-PCR. Pairwise comparison suggested that 61 and 44 miRNAs were differentially expressed at 4 and 24 h poi, respectively. Furthermore, the expression profiles of nine randomly selected miRNAs were validated by qRT-PCR. MicroRNA target prediction, gene ontology (GO) annotation, and Kyoto Encylopedia of Genes and Genomes (KEGG) analysis indicated that a variety of biological pathways could be affected by A. hydrophila infection. Additionally, transferrin (TF) and transferrin receptor (TFR) genes were confirmed to be direct targets of miR-375. These results will expand our knowledge of the role of miRNAs in the immune response of M. amblycephala to A. hydrophila infection, and facilitate the development of effective strategies against A. hydrophila infection in M. amblycephala.

Comprehensive identification and profiling of Nile tilapia (Oreochromis niloticus) microRNAs response to Streptococcus agalactiae infection through high-throughput sequencing

MicroRNAs are a kind of small non-coding RNAs that participate in various biological processes. Deregulated microRNA expression is associated with several types of diseases. Tilapia (Oreochromis niloticus) is an important commercial fish species in China. To identify miRNAs and investigate immune-related miRNAs of O. niloticus, we applied high-throughput sequencing technology to identify and analyze miRNAs from tilapia infected with Streptococcus agalactiae at a timescale of 72 h divided into six different time points. The results showed that a total of 3009 tilapia miRNAs were identified, including in 1121 miRNAs which have homologues in the currently available databases and 1878 novel miRNAs. The expression levels of 218 tilapia miRNAs were significantly altered at 6 h-72 h post-bacterial infection (pi), and these miRNAs were therefore classified as differentially expressed tilapia miRNAs. For the 1121 differentially expressed tilapia miRNAs target 41961 genes. GO and KEGG enrichment analysis revealed that some target genes of tilapia miRNAs were grouped mainly into the categories of apoptotic process, signal pathway, and immune response. This is the first report of comprehensive identification of O. niloticus miRNAs being differentially regulated in spleen in normal conditions relating to S. agalactiae infection. This work provides an opportunity for further understanding of the molecular mechanisms of miRNA regulation in O. niloticus host-pathogen interactions.

Identification and profiling of growth-related microRNAs of the swimming crab Portunus trituberculatus by using Solexa deep sequencing

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of mRNAs. The swimming crab Portunus trituberculatus is one of the most important crustacean species for aquaculture in China. However, to date no miRNAs have been reported to for modulating growth in P. trituberculatus. To investigate miRNAs involved in the growth of this species, we constructed six small RNA libraries for big individuals (BIs) and small individuals (SIs) from a highly inbred family. Six mixed RNA pools of five tissues (eyestalk, gill, heart, hepatopancreas, and muscle) were obtained. By aligning sequencing data with those for known miRNAs, a total of 404 miRNAs, including 339 known and 65 novel miRNAs, were identified from the six libraries. MiR-100 and miR-276a-3p were among the most prominent miRNA species. We identified seven differentially expressed miRNAs between the BIs and SIs, which were validated using real-time PCR. Preliminary analyzes of their putative target genes and GO and KEGG pathway analyzes showed that these differentially expressed miRNAs could play important roles in global transcriptional depression and cell differentiation of P. trituberculatus. This study reveals the first miRNA profile related to the body growth of P. trituberculatus, which would be particularly useful for crab breeding programs.

Identification and characterization of immune-related microRNAs in blunt snout bream, Megalobrama amblycephala

MicroRNAs (miRNAs) play vital roles in diverse biological processes, including in immune response. Blunt snout bream (Megalobrama amblycephala) is a prevalent and important commercial endemic freshwater fish species in China's intensive polyculture systems. To identify immune-related miRNAs of M. amblycephala, two small RNA (sRNA) libraries from immune tissues with or without lipopolysaccharide (LPS) stimulation were constructed and sequenced using the high-throughput sequencing technology. Totally, 16,425,543 and 15,076,813 raw reads, corresponding to 14,156,755 and 13,445,869 clean reads, were obtained in the normal and infected libraries, respectively. A total of 324 miRNAs, including 218 known miRNAs and 106 putative novel miRNAs were identified by bioinformatic analysis. We analyzed differentially expressed miRNAs between two libraries using pairwise comparison. 113 (34.88%) miRNAs were found to be significantly differentially expressed between two libraries, with 63 (55.75%) exhibiting elevated expression in LPS stimulation sample. Thereinto, a number of known miRNAs were identified immune-related. Real-time quantitative PCR (RT-qPCR) were implemented for 12 miRNAs of two samples, and agreement was confirmed between the sequencing and RT-qPCR data. Target genes likely regulated by these differentially expressed miRNAs were predicted using computational prediction. The functional annotation of target genes by Gene Ontology enrichment (GO) and Kyoto Encyclopedia of Genes and Genomes pathway analysis (KEGG) indicated that a majority of differential miRNAs might involved in immune response. To our knowledge, this is the first comprehensive study of miRNAs in response to LPS stimulation in M. amblycephala, even in fish. These results deepened our understanding of the role of miRNAs in the intricate host's immune system, and should be useful to develop new control strategies for host immune defense against various bacterial invasions in M. amblycephala.

Expression Profile Analysis of miR-221 and miR-222 in Different Tissues and Head Kidney Cells of Cynoglossus semilaevis, Following Pathogen Infection

Half-smooth tongue sole (Cynoglossus semilaevis) is an important marine commercial fish species in China, which suffers from widespread disease outbreaks. Recently, in this regard, our group identified immune-related microRNAs (miRNAs) of C. semilaevis following Vibrio anguillarum infection. Furthermore, miRNA microarray was utilized to characterize the immune roles of important miRNA candidates in response to bacterial infection. Therefore, in the present study, we characterized miR-221 and miR-222 and profiled their expression after challenge. Here, miR-221 and miR-222 precursors were predicted to have a typical hairpin structure. Both miRNAs were expressed in a broad range of tissues in C. semilaevis, while miR-221 and miR-222 were significantly differentially expressed in the immune tissues of C. semilaevis among three small RNA libraries [control group (CG), bacteria-challenged fish without obvious symptoms of infection (NOSG), and bacteria-challenged fish with obvious symptoms of infection (HOSG)]. In order to further characterize and understand the immune response of miR-221 and miR-222, therefore, we profiled miR-221 and miR-222 expression in selected immune tissues after challenge with V. anguillarum. Both miR-221 and miR-222 were upregulated in the liver and spleen, while different expression patterns were observed in the head kidney. In addition, in half-smooth tongue sole head kidney cell line after challenge with lipopolysaccharide (LPS), polyinosinic:polycytidylic acid (poly I:C), peptidoglycan (PGN), and red-spotted grouper nervous necrosis virus (RGNNV), both miR-221 and miR-222 showed significant difference in expression response to pathogen. Meanwhile, the target gene of miR-221 and miR-222 was predicted, which indicated that tumor necrosis factor receptor-associated factor 6 (TRAF6) and interleukin-1 beta (IL-1β) were the target genes of miR-221 and miR-222, respectively. Collectively, these findings indicated that miR-221 and miR-222 have putative roles in innate immune response during C. semilaevis exposure to pathogens. Our findings could expand the knowledge of immune function of C. semilaevis miRNA and guide future studies on C. semilaevis immunity.