Expression profiling analysis of the microRNA response of Cynoglossus semilaevis to Vibrio anguillarum and other stimuli

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.

Integrated analysis of microRNA and mRNA expression profiles in Preeclampsia

BACKGROUND

Preeclampsia (PE), a pregnancy specific syndrome, is one kind of common gestational hypertension disease, which can cause maternal and perinatal mortality and morbidity. This study was conducted to identify key microRNAs (miRNAs), mRNAs and related signaling pathways in the pathogenesis of PE.

METHODS

Whole transcriptome sequencing and small RNA sequencing of the peripheral blood from 3 PE patients and 3 normal pregnant women were performed. Differential expressed (DE) miRNAs were identified using the DEseq2 package. Target genes of the selected upregulated and downregulated DE miRNAs were predicted. Based on the hypergeometric distribution of DE miRNA target genes, we analyzed GO enrichment and KEGG pathway enrichment using R.

RESULTS

Total 1291 and 1281 novel RNAs were obtained from the preeclampsia patients and healthy individuals. 70 miRNAs were screened out with significant levels with 51 significantly upregulated and 19 significantly downregulated. 44,306 genes were predicted as the targets of these miRNAs. Besides, KEGG pathway analysis revealed that the upregulated miRNAs were enriched in Glycosaminoglycan biosynthesis-chondroitin sulfate / dermatan sulfate, Base excision repair and the downregulated miRNAs were enriched in Tuberculosis, Phagosome.

CONCLUSION

We constructed regulatory networks of miRNAs and target genes, there were 2208 negative miRNA-mRNA interactions in total. The network and pathway information illustrate the potential functions of mRNAs and miRNAs in PE pathogenesis.

Acute cold stress leads to zebrafish ovarian dysfunction by regulating miRNA and mRNA

Temperature is a critical factor that regulates the reproduction processes in teleost. However, the gonadal response mechanism to cold stress in fish remains largely unknown. In the present study, female zebrafish were exposed to different extents of low temperatures at 18 °C and 10 °C for 48 h. The ovarian histology was remarkably damaged after cold stress exposure. Integrated analysis of miRNA-mRNA was used to investigate the ovarian response to acute cold stress. A large number of mRNAs and miRNAs were altered by cold stress, which are involved in extensive biological processes. It is indicated that the signal transduction of MAPK and Calcium signaling pathway is highly engaged in zebrafish ovary to adapt to cold stress. The immune system was dysregulated by cold stress while the ovarian autophagy was activated. Remarkably increased gene number related to reproductive functions was identified in the cold stress at 10 °C compared to the control. The cold stress-induced dysregulated reproductive genes include star, hsd3b1, hsd17b1, inha, insl3, amh, nanos1 and foxl2. Combined with the dysregulated insulin, IGF and progesterone signaling, it is suggested that cold stress affects ovarian function in multiple aspects, including oocyte meiosis, folliculogenesis, final maturation and ovarian maintenance. On the other hand, the ovarian miRNA-mRNA regulatory network response to cold stress was also constructed. Overall, our result revealed the ovarian response to cold stress in zebrafish and provided insight into the fish adaptation mechanism to acute temperature change.

Transcriptomic Analysis Reveals Functional Interaction of mRNA-lncRNA-miRNA in Trachinotus ovatus Infected by Cryptocaryon irritans

The skin of Trachinotus ovatus is a crucial component of the mucosal immune system and serves as the primary site of infection by Cryptocaryon irritans. In order to investigate the significant role of skin in C. irritans infection, a comprehensive transcriptome analysis was conducted on skin tissues from the infection group, infection-adjacent group, and infection group compared with the infection-adjacent group (ATT_vs_PER, ADJ_vs_PER, ATT_vs_ADJ). This study identified differentially expressed long non-coding RNAs (DE lncRNAs), microRNAs (DE miRNAs), and differentially expressed genes (DEGs). The prediction of lncRNA target genes was accomplished by utilizing positional relationship (co-location) and expression correlation (co-expression) with protein-coding genes. Subsequently, functional enrichment analysis was conducted on the target genes of differentially expressed lncRNAs, revealing their involvement in signaling pathways such as tight junction, MAPK, and cell adhesion molecules. This study describes the regulatory network of lncRNA-miRNA-mRNA in T. ovatus skin tissue infected with C. irritans. Functional prediction analysis showed that differentially expressed lncRNA and miRNA may regulate the expression of immune genes such as interleukin-8 (il8) to resist the infection of C. irritans. Conducting additional research on these non-coding RNAs will facilitate a deeper understanding of their immune regulatory function in T. ovatus during C. irritans infection. The study of non-coding RNA in this study laid a foundation for revealing the molecular mechanism of the immune system of T. ovatus to respond to the infection of C. irritans. It provided a choice for the molecular breeding of Trachinotus ovatus against C. irritans.

The immune-related circRNA-miRNA-mRNA ceRNA regulatory network in the liver of turbot (Scophthalmus maximus L.) induced by Vibrio anguillarum

Recently, Vibrio anguillarum, a Gram-negative pathogenic bacterium, has been becoming a major constraint on the development of the turbot aquaculture industry because of its characteristics of worldwide distribution, broad host range and potentially devastating impacts. Although the functions of protein-coding mRNAs in the immune response against bacterial infection have been reported, as well as several non-coding RNAs (ncRNAs), such as circular RNAs (circRNAs) and microRNAs (miRNAs), the relationships between mRNAs and ncRNAs in the immune system of turbot liver are still limited during bacterial infection. In present study, the comprehensive analyses of whole-transcriptome sequencing were conducted in turbot liver infected by V. anguillarum. The differential expression was analyzed in the data of circRNAs, miRNAs, and mRNAs. The interactions of miRNA-circRNA pairs and miRNA-mRNA pairs were predicted basing on the negative regulatory relationships between miRNAs and their target circRNAs\mRNAs. The circRNA-related ceRNA regulatory networks were constructed for the analyses of regulated mechanism in turbot immune system. Subsequently, the RT-qPCR was carried out to verify the results of sequencing. Finally, we identified 31 circRNAs, 53 miRNAs and 948 mRNAs with differential expression. Gene set enrichment analyses using Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways showed that innate immunity was principally activated at the early stages of infection, while adaptive immunity was activated after 24 h. Finally, 65 circRNA-miRNA-mRNA pathways were constructed, based on the hypothesis of ceRNA regulatory networks. In conclusion, our findings provide new insights on the underlying immune response to bacterial infection and identify novel target genes for the prevention and control of disease in turbot.

Integrative analysis of miRNA and mRNA expression associated with the immune response in the intestine of rainbow trout (Oncorhynchus mykiss) infected with infectious hematopoietic necrosis virus

Rainbow trout (Oncorhynchus mykiss), an economically important cold-water fish cultured worldwide, suffers from infectious hematopoietic necrosis virus (IHNV) infection, resulting in huge financial losses. In order to understand the immune response of rainbow trout during virus infection, we explored trout intestine transcriptome profiles following IHNV challenge, and identified 3355 differentially expressed genes (DEGs) and 80 differentially expressed miRNAs (DEMs). Transcriptome analysis revealed numerous DEGs involved in immune responses, such as toll-like receptor 3 (TLR3), toll-like receptor 7/8 (TLR7/8), tripartite motif-containing 25 (TRIM25), DExH-Box helicase 58 (DHX58), interferon-induced with helicase C domain 1 (IFIH1), interferon regulatory factor 3 (IRF3/7), signal transducer and activator of transcription 1 (STAT1) and heat shock protein 90-alpha 1 (HSP90A1). Integrated analysis identified five key miRNAs (miR-19-y, miR-181-z, miR-203-y, miR-143-z and miR-206-y) targeting at least two important immune genes (TRIM25, DHX58, STAT1, TLR7/8 and HSP90A1). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that DEGs and target genes were significantly enriched in various immune-related terms including immune system process, binding, cell part and pathways of Toll-like receptor signalling, RIG-I-like receptor signalling, NOD-like receptor signalling, JAK-STAT signalling, PI3K-Akt signalling, NF-kappa B signalling, IL-17 signalling and AGE-RAGE signalling. In addition, protein-protein interaction networks (PPI) was used to display highly interactive DEG networks involving eight immune-related pathways. The expression trends of 12 DEGs and 10 DEMs were further verified by quantitative real-time PCR, which confirmed the reliability of the transcriptome sequencing results. This study expands our understanding of the immune response of rainbow trout infected with IHNV, and provides valuable resources for future studies on the immune molecular mechanism and disease resistance breeding.

Integrated analysis of the miRNA and mRNA expression profiles in Leiocassis longirostris at gonadal maturation

Leiocassis longirostris is a commercially important fish species that shows a sexually dimorphic growth pattern. A lack of molecular data from the gonads of this species has hindered research and selective breeding efforts. In this study, we conducted a comprehensive analysis of the expression profile of miRNA and mRNA to explore their regulatory roles in the gonadal maturation stage of L. longirostris. We identified 60 differentially expressed miRNAs and 20,752 differentially expressed genes by sequencing. A total of 90 miRNAs and 21 target genes involved in gonad development and sex determination were identified. Overall, the results of this study enhance our understanding of the molecular mechanisms underlying sex determination and differentiation and provide valuable genomic information for the selective breeding of L. longirostris.

Comparative analysis of the miRNA-mRNA regulation networks in turbot (Scophthalmus maximus L.) following Vibrio anguillarum infection

MicroRNAs could not only regulate posttranscriptional silencing of target genes in eukaryotic organisms, but also have positive effect on their target genes as well. These microRNAs have been reported to be involved in mucosal immune responses to pathogen infection in teleost. Therefore, we constructed the immune-related miRNA-mRNA networks in turbot intestine following Vibrio anguillarum infection. In our results, 1550 differentially expressed (DE) genes and 167 DE miRNAs were identified. 113 DE miRNAs targeting 89 DE mRNAs related to immune response were used to construct miRNA-mRNA interaction networks. Functional analysis showed that target genes were associated with synthesis and degradation of ketone bodies, mucin type O-Glycan biosynthesis, homologous recombination, biotin metabolism, and intestinal immune network for IgA production that were equivalent to the function of IgT and IgM in fish intestine. Finally, 10 DE miRNAs and 7 DE mRNAs were selected for validating the accuracy of high-throughput sequencing results by qRT-PCR. The results of this study will provide valuable information for the elucidation of the regulation mechanisms of miRNA-mRNA interactions involved in disease resistance in teleost mucosal immune system.

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.

Research progress on the regulation of nutrition and immunity by microRNAs in fish

MicroRNAs (miRNAs) are a class of highly conserved, endogenous non-coding single-stranded small RNA molecules with a length of 18-25 nucleotides. MiRNAs can negatively regulate the target gene through complementary pairing with the mRNA. It has been more than 20 years since the discovery of miRNA molecules, and many achievements have been made in fish research. This paper reviews the research progress in the regulation of fish nutrition and immunity by miRNAs in recent years. MiRNAs regulate the synthesis of long-chain polyunsaturated fatty acids, and are involved in the metabolism of glucose, lipids, as well as cholesterol in fish. Moreover, miRNAs play various roles in antibacterial and antiviral immunity of fish. They can promote the immune response of fish, but may also participate in the immune escape mechanism of bacteria or viruses. One important aspect of miRNAs regulation on fish immunity is mediated by targeting pattern recognition receptors and downstream signaling factors. Together, current results indicate that miRNAs are widely involved in the complex regulatory network of fish. Further studies on fish miRNAs may deepen our understanding of the regulatory network of fish nutrition and immunity, and have the potential to promote the development of microRNA-based products and detection reagents that can be applied in aquaculture industry.

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.

Integrated Analysis of circRNA-miRNA-mRNA Regulatory Networks in the Intestine of Sebastes schlegelii Following Edwardsiella tarda Challenge

Sebastes schlegelii, an important aquaculture species, has been widely cultured in East Asian countries. With the increase in the cultivation scale, various diseases have become major threats to the industry. Evidence has shown that non-coding RNAs (ncRNAs) have remarkable functions in the interactions between pathogens and their hosts. However, little is known about the mechanisms of circular RNAs (circRNAs) and coding RNAs in the process of preventing pathogen infection in the intestine in teleosts. In this study, we aimed to uncover the global landscape of mRNAs, circRNAs, and microRNAs (miRNAs) in response to Edwardsiella tarda infection at different time points (0, 2, 6, 12, and 24 h) and to construct regulatory networks for exploring the immune regulatory mechanism in the intestine of S. schlegelii. In total, 1,794 mRNAs, 87 circRNAs, and 79 miRNAs were differentially expressed. The differentially expressed RNAs were quantitatively validated using qRT-PCR. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that most of the differentially expressed mRNA genes and the target genes of ncRNAs were related to immune signaling pathways, such as the NF-κB signal pathway, pathogen recognition receptors related to signaling pathways (Toll-like receptors and Nod-like receptors), and the chemokine signaling pathway. Based on these differentially expressed genes, 624 circRNA-miRNA pairs and 2,694 miRNA-mRNA pairs were predicted using the miRanda software. Integrated analyses generated 25 circRNA-miRNA-mRNA interaction networks. In a novel_circ_0004195/novel-530/IκB interaction network, novel_530 was upregulated, while its two targets, novel_circ_0004195 and IκB, were downregulated after E. tarda infection. In addition, two circRNA-miRNA-mRNA networks related to apoptosis (novel_circ_0003210/novel_152/apoptosis-stimulating of p53 protein 1) and interleukin (novel_circ_0001907/novel_127/interleukin-1 receptor type 2) were also identified in our study. We thus speculated that the downstream NF-κB signaling pathway, p53 signaling pathway, and apoptosis pathway might play vital roles in the immune response in the intestine of S. schlegelii. This study revealed a landscape of RNAs in the intestine of S. schlegelii during E. tarda infection and provided clues for further study on the immune mechanisms and signaling networks based on the circRNA-miRNA-mRNA axis in S. schlegelii.

Systematic Identification and Analysis of Circular RNAs of Japanese Flounder (Paralichthys olivaceus) in Response to Vibrio anguillarum Infection

Circular RNA (circRNA) is a new class of non-coding RNA that is structured into a closed loop without polyadenylation. Recent studies showed that circRNAs are involved in the host immune response to pathogen infection. Japanese flounder (Paralichthys olivaceus), an important economical marine fish cultured in north Asia, is affected by Vibrio anguillarum, a pathogenic bacterium that can infect a large number of fish. In this study, we systematically explored the circRNAs in the spleen of V. anguillarum-infected flounder at different infection time points. A total of 6581 circRNAs were identified, 148 of which showed differential expression patterns after V. anguillarum infection and were named DEcirs. Most of the DEcirs were strongly time-specific. The parental genes of the DEcirs were identified and functionally classified into diverse pathways, including immune-related pathways. Among the immune-related DEcirs, seven were predicted to sponge 18 targeted miRNAs that were differentially expressed during V. anguillarum infection (named DETmiRs). Further analysis showed that the DEcirs and their corresponding DETmiRs intertwined into complicated immune related networks. These results indicate that in flounder, circRNAs are regulated by V. anguillarum and form interactive networks with mRNAs and miRNAs that likely play important roles in the immune defense against pathogen infection.

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.

Pol-miR-363-3p plays a significant role in the immune defense of Japanese flounder Paralichthys olivaceus against bacterial and viral infection

In this study, we examined the function of a Japanese flounder (Paralichthys olivaceus) microRNA (miRNA), pol-miR-363-3p. We found that pol-miR-363-3p targets an ubiquitin-specific protease (USP), USP32. USP is a family of deubiquitinating enzymes essential to the functioning of the ubiquitin proteasome system. In mammals, USP32 is known to be associated with cancer and immunity. In fish, the function of USP32 is unknown. We found that flounder USP32 (PoUSP32) expression was detected in the major tissues of flounder, particularly intestine. In vitro and in vivo studies showed that pol-miR-363-3p directly regulated PoUSP32 in a negative manner by interaction with the 3'UTR of PoUSP32. Overexpression of pol-miR-363-3p or interference with PoUSP32 expression in flounder cells significantly blocked Streptococcus iniae infection. Consistently, in vivo knockdown of pol-miR-363-3p or overexpression of PoUSP32 enhanced dissemination of S. iniae in flounder tissues, whereas in vivo knockdown of PoUSP32 inhibited S. iniae dissemination. In addition, pol-miR-363-3p knockdown also significantly promoted the tissue dissemination of the viral pathogen megalocytivirus, which, as well as S. iniae, regulated pol-miR-363-3p expression. Together these results revealed an important role of pol-miR-363-3p in flounder immune defense against bacterial and viral infection.

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.

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.

MicroRNAs Involved in the Regulation of LC-PUFA Biosynthesis in Teleosts: miR-33 Enhances LC-PUFA Biosynthesis in Siganus canaliculatus by Targeting insig1 which in Turn Upregulates srebp1

Post-transcriptional regulatory mechanisms play important roles in the regulation of LC-PUFA biosynthesis. Our previous study revealed that miR-33 could increase the expression of fatty acyl desaturases (fads2) in the rabbitfish Siganus canaliculatus, but the specific mechanism is unknown. Here, we confirmed that miR-33 could target the 3'UTR of insulin-induced gene 1 (insig1), resulting in downregulation of its protein level in the rabbitfish hepatocyte line (SCHL). In vitro overexpression of miR-33 inhibited the mRNA level of insig1 and increased the mRNA levels of Δ6Δ5 fads2 and elovl5, as well as srebp1. In SCHL cells, proteolytic activation of sterol-regulatory-element-binding protein-1 (Srebp1) was blocked by Insig1, with overexpression of insig1 decreasing mature Srebp1 level, while inhibition of insig1 led to the opposite effect. Srebp1 could enhance the promoter activity of Δ6Δ5 fads2 and elovl5, whose expression levels decreased with knockdown of srebp1 in SCHL. Overexpression of miR-33 also resulted in a higher conversion of 18:3n-3 to 18:4n-3 and 20:5n-3 to 22:5n-3, linked to desaturation and elongation via Δ6Δ5 Fads2 and Elovl5, respectively. The results suggested that the mechanism by which miR-33 regulates LC-PUFA biosynthesis in rabbitfish is through enhancing the expression of srebp1 by targeting insig1. The findings here provide more insight to the mechanism of miRNAs involvement in the regulation of LC-PUFA biosynthesis in teleosts.

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.

Transcriptomic analysis of microRNAs-mRNAs regulating innate immune response of zebrafish larvae against Vibrio parahaemolyticus infection

In recent years, microRNAs (miRNAs) have been shown to play important roles in immunity. Analyses of the functions of miRNAs and their targets are useful in understanding the regulation of the immune response. To understand the relationships between miRNAs and their targets during infection, we used zebrafish as an infection model in which to characterize the miRNA and mRNA transcriptomes of zebrafish larvae infected with Vibrio parahaemolyticus. We identified the differentially expressed miRNAs and mRNAs. Overall, 37 known zebrafish miRNAs were differentially expressed in the infection group and 107 predicted target genes of 26 miRNAs were differentially expressed in the mRNA transcriptome. These targets with specific Gene Ontology (GO) terms, such as peripheral nervous system neuron axonogenesis, organophosphate metabolic process, heme binding, protein binding, tetrapyrrole binding, protein dimerization activity, and aromatase activity, which regulate nerve conduction, energy metabolism, hematopoiesis, and protein synthesis. They were also associated with Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways such as phototransduction, tryptophan metabolism, notch signaling, and purine metabolism. Our findings indicate that miRNAs regulate the innate immune response via complex networks, and zebrafish (Danio rerio, dre)-miR-205-3p, dre-miR-141-5p, dre-miR-200a-5p, dre-miR-92a-2-5p, dre-miR-192, and dre-miR-1788 may play important roles in the innate immune response by regulating target genes.

Network of microRNA-transcriptional factor-mRNA in cold response of turbot Scophthalmus maximus

The aim of this study is to understand fish cold-tolerant mechanism. We analyzed the transcriptional reactions to the cold condition in turbot Scophthalmus maximus by using RNA-seq and microRNA (miRNA)-seq. Meio-gynogenetic diploid turbots were treated at 0 °C to distinguish the cold-tolerant (CT) and cold-sensitive (CS) groups. The results showed that there were quite different responses at both mRNA and miRNA levels, with more up-regulated mRNAs (1069 vs. 194) and less down-regulated miRNAs (4 vs. 1) in CT versus CS relative to the control group. The network of miRNA-transcription factor-mRNA, regulating turbot different response to cold stress, was constructed, which involved in cell cycle, component of cell membrane, signal transduction, and circadian rhythm pathways. The above information demonstrates mechanisms by which cold tolerance is increased in fish.

Identification and characteristics of muscle growth-related microRNA in the Pacific abalone, Haliotis discus hannai

BACKGROUND

The Pacific abalone, Haliotis discus hannai, is the most important cultivated abalone in China. Improving abalone muscle growth and increasing the rate of growth are important genetic improvement programs in this industry. MicroRNAs are important small noncoding RNA molecules that regulate post-transcription gene expression. However, no miRNAs have been reported to regulate muscle growth in H. discus hannai.

RESULTS

we profiled six small RNA libraries for three large abalone individuals (L_HD group) and three small individuals (S_HD group) using RNA sequencing technology. A total of 205 miRNAs, including 200 novel and 5 known miRNAs, were identified. In the L_HD group, 3 miRNAs were up-regulated and 7 were down-regulated compared to the S_HD specimens. Bioinformatics analysis of miRNA target genes revealed that miRNAs participated in the regulation of cellular metabolic processes, the regulation of biological processes, the Wnt signaling pathway, ECM-receptor interaction, and the MAPK signaling pathway, which are associated with regulating growth. Bone morphogenetic protein 7 (BMP7) was verified as a target gene of hdh-miR-1984 by a luciferase reporter assay and we examined the expression pattern in different developmental stages.

CONCLUSION

This is the first study to demonstrate that miRNAs are related to the muscle growth of H. discus hannai. This information could be used to study the mechanisms of abalone muscle growth. These DE-miRNAs may be useful as molecular markers for functional genomics and breeding research in abalone and closely related species.

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.

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.

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila)

Background

The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper.

Results

A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data.

Conclusions

The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3995-9) contains supplementary material, which is available to authorized users.

Long Non-Coding RNAs (lncRNAs) of Sea Cucumber: Large-Scale Prediction, Expression Profiling, Non-Coding Network Construction, and lncRNA-microRNA-Gene Interaction Analysis of lncRNAs in Apostichopus japonicus and Holothuria glaberrima During LPS Challenge and Radial Organ Complex Regeneration

Long non-coding RNA (lncRNA) structurally resembles mRNA but cannot be translated into protein. Although the systematic identification and characterization of lncRNAs have been increasingly reported in model species, information concerning non-model species is still lacking. Here, we report the first systematic identification and characterization of lncRNAs in two sea cucumber species: (1) Apostichopus japonicus during lipopolysaccharide (LPS) challenge and in heathy tissues and (2) Holothuria glaberrima during radial organ complex regeneration, using RNA-seq datasets and bioinformatics analysis. We identified A. japonicus and H. glaberrima lncRNAs that were differentially expressed during LPS challenge and radial organ complex regeneration, respectively. Notably, the predicted lncRNA-microRNA-gene trinities revealed that, in addition to targeting protein-coding transcripts, miRNAs might also target lncRNAs, thereby participating in a potential novel layer of regulatory interactions among non-coding RNA classes in echinoderms. Furthermore, the constructed coding-non-coding network implied the potential involvement of lncRNA-gene interactions during the regulation of several important genes (e.g., Toll-like receptor 1 [TLR1] and transglutaminase-1 [TGM1]) in response to LPS challenge and radial organ complex regeneration in sea cucumbers. Overall, this pioneer systematic identification, annotation, and characterization of lncRNAs in echinoderm pave the way for similar studies and future genetic, genomic, and evolutionary research in non-model species.

High throughput sequencing of small RNAs transcriptomes in two Crassostrea oysters identifies microRNAs involved in osmotic stress response

Increasing evidence suggests that microRNAs post-transcriptionally regulate gene expression and are involved in responses to biotic and abiotic stress. However, the role of miRNAs involved in osmotic plasticity remains largely unknown in marine bivalves. In the present study, we performed low salinity challenge with two Crassostrea species (C. gigas and C. hongkongensis), and conducted high-throughput sequencing of four small RNA libraries constructed from the gill tissues. A total of 202 and 87 miRNAs were identified from C. gigas and C. hongkongensis, respectively. Six miRNAs in C. gigas and two in C. hongkongensis were differentially expressed in response to osmotic stress. The expression profiles of these eight miRNAs were validated by qRT-PCR. Based on GO enrichment and KEGG pathway analysis, genes associated with microtubule-based process and cellular component movement were enriched in both species. In addition, five miRNA-mRNA interaction pairs that showed opposite expression patterns were identified in the C. hongkongensis, Differential expression analysis identified the miRNAs that play important regulatory roles in response to low salinity stress, providing insights into molecular mechanisms that are essential for salinity tolerance in marine bivalves.

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.

High-Throughput Profiling of Caenorhabditis elegans Starvation-Responsive microRNAs

MicroRNAs (miRNAs) are non-coding RNAs of ~22 nucleotides in length that regulate gene expression by interfering with the stability and translation of mRNAs. Their expression is regulated during development, under a wide variety of stress conditions and in several pathological processes. In nature, animals often face feast or famine conditions. We observed that subjecting early L4 larvae from Caenorhabditis elegans to a 12-hr starvation period produced worms that are thinner and shorter than well-fed animals, with a decreased lipid accumulation, diminished progeny, reduced gonad size, and an increased lifespan. Our objective was to identify which of the 302 known miRNAs of C. elegans changed their expression under starvation conditions as compared to well-fed worms by means of deep sequencing in early L4 larvae. Our results indicate that 13 miRNAs (miR-34-3p, the family of miR-35-3p to miR-41-3p, miR-39-5p, miR-41-5p, miR-240-5p, miR-246-3p and miR-4813-5p) were upregulated, while 2 miRNAs (let-7-3p and miR-85-5p) were downregulated in 12-hr starved vs. well-fed early L4 larvae. Some of the predicted targets of the miRNAs that changed their expression in starvation conditions are involved in metabolic or developmental process. In particular, miRNAs of the miR-35 family were upregulated 6–20 fold upon starvation. Additionally, we showed that the expression of gld-1, important in oogenesis, a validated target of miR-35-3p, was downregulated when the expression of miR-35-3p was upregulated. The expression of another reported target, the cell cycle regulator lin-23, was unchanged during starvation. This study represents a starting point for a more comprehensive understanding of the role of miRNAs during starvation in C. elegans.