Identification and differential expression of microRNAs during metamorphosis of the Japanese flounder (Paralichthys olivaceus)

Isolation and identification of vapA-absent Aeromonas salmonicida in diseased snakehead Channa argus in China

Aeromonas salmonicida is the typical pathogen causing furunculosis, reported widely in salmonids. Because of multiple serotypes, the control of A. salmonicida-caused disease has increasingly received much attention. Recently, A. salmonicida infection was reported in non-salmonid fish species. Here, a pathogenic A. salmonicida, named as As-s, was isolated from cultured snakehead (Channa argus) in a local fish farm in Shandong, China. As-s displayed clear hemolysis, amylase, and positive catalase activities, and grew at a wide range of temperatures (10-37 °C) and pH values (5.5-8.5). As-s was highly sensitive to cefuroxime sodium, ceftriaxone, ceftazidime, piperacillin, and cefoperazone and also apparently sensitive to chloramphenicol, erythromycin, and 25% cinnamaldehyde. The Virulence array protein gene cloning' results suggested that As-s has this gene compared with the other two vapA-containing strains, despite a close relationship of these strains via phylogenetic analysis. Severe ulcers on skin, muscle, and abnormal liver, and hemorrhage in pectoral/ventral fins and anal region were observed, and exophthalmos were also noticed in infected juvenile snakehead, as well as necrosis and infiltration of blood cells emerged in the internal organs using pathological section. In addition, As-s caused high mortality in snakehead, consistently with its immune gene response. This study reports the first isolation of vapA-absent A. salmonicida in snakehead.

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

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

High-Temperature-Induced Differential Expression of miRNA Mediates Liver Inflammatory Response in Tsinling Lenok Trout (Brachymystax lenok tsinlingensis)

High-temperature stress poses a significant environmental challenge for aquatic organisms, including tsinling lenok trout (Brachymystax lenok tsinlingensis). This study aimed to investigate the role of microRNAs (miRNAs) in inducing liver inflammation in tsinling lenok trout under high-temperature stress. Tsinling lenok trout were exposed to high-temperature conditions (24 °C) for 8 h, and liver samples were collected for analysis. Through small RNA sequencing, we identified differentially expressed miRNAs in the liver of high-temperature-stressed tsinling lenok trout compared to the control group (maintained at 16 °C). Several miRNAs, including novel-m0105-5p and miR-8159-x, showed significant changes in expression levels. Additionally, we conducted bioinformatics analysis to explore the potential target genes of these differentially expressed miRNAs. Our findings revealed that these miRNA target genes are involved in inflammatory response pathways, such as NFKB1 and MAP3K5. The downregulation of novel-m0105-5p and miR-8159-x in the liver of high-temperature-stressed tsinling lenok trout suggests their role in regulating liver inflammatory responses. To validate this, we performed a dual-luciferase reporter assay to confirm the regulatory relationship between miRNAs and target genes. Our results demonstrated that novel-m0105-5p and miR-8159-x enhance the inflammatory response of hepatocytes by promoting the expression of NFKB1 and MAP3K5, respectively. In conclusion, our study provides evidence that high-temperature stress induces liver inflammation in tsinling lenok trout through dysregulation of miRNAs. Understanding the molecular mechanisms underlying the inflammatory response in tsinling lenok trout under high-temperature stress is crucial for developing strategies to mitigate the negative impacts of environmental stressors on fish health and aquaculture production.

miR-130a targets CiGadd45bb to modulate the inflammatory response to bacterial infection in Ctenopharyngodon idella kidney (CIK) cells

Septicemia is a systemic inflammatory response to bacterial infection that results in a hyper-inflammatory state, which could lead to septic shock and death in grass carp (Ctenopharyngodon idella). The aim of this study was to determine the underlying mechanism of microRNA (miR-130a) in bacteria-infected grass carp. Expression levels of miR-130a against Aeromonas hydrophila (A. hydrophila) infection in Ctenopharyngodon idella kidney cells (CIK) were analyzed. Luciferase reporter assay, quantitative reverse transcription-polymerase chain reaction were performed to explore whether Ctenopharyngodon idella growth arrest and DNA damage-inducible 45 (CiGadd45bb) was a target of miR-130a. MiR-130a mimic, inhibitor and miR-control were transfected to CIK respectively. After transfection, the expression levels of proinflammatory genes were determined. Here we show that CiGadd45bb as a target of miR-130a. We also confirmed that miR-130a levels were significantly higher after being stimulated for 4 h and lower after 12 h (P < 0.01). Overexpressing miR-130a strikingly inhibited p38, JNK, ERK and TNF-a genes (P < 0.01) and silencing miR-130a activated p38, JNK, ERK, TNF-a, IFN and IL-8 (P < 0.01). Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by miR-130a in grass carp.

Characterization of miRNAs in Embryonic, Larval, and Adult Lumpfish Provides a Reference miRNAome for Cyclopterus lumpus

Simple Summary

Lumpfish (Cyclopterus lumpus) is an emergent aquaculture species, and its miRNA repertoire is still unknown. miRNAs are critical post-transcriptional modulators of teleost gene expression. Therefore, a lumpfish reference miRNAome was characterized by small RNA sequencing and miRDeep analysis of samples from different organs and developmental stages. The resulting miRNAome, an essential reference for future expression analyses, consists of 443 unique mature miRNAs from 391 conserved and eight novel miRNA genes. Enrichment of specific miRNAs in particular organs and developmental stages indicates that some conserved lumpfish miRNAs regulate organ and developmental stage-specific functions reported in other teleosts.

Abstract

MicroRNAs (miRNAs) are endogenous small RNA molecules involved in the post-transcriptional regulation of protein expression by binding to the mRNA of target genes. They are key regulators in teleost development, maintenance of tissue-specific functions, and immune responses. Lumpfish (Cyclopterus lumpus) is becoming an emergent aquaculture species as it has been utilized as a cleaner fish to biocontrol sea lice (e.g., Lepeophtheirus salmonis) infestation in the Atlantic Salmon (Salmo salar) aquaculture. The lumpfish miRNAs repertoire is unknown. This study identified and characterized miRNA encoding genes in lumpfish from three developmental stages (adult, embryos, and larvae). A total of 16 samples from six different adult lumpfish organs (spleen, liver, head kidney, brain, muscle, and gill), embryos, and larvae were individually small RNA sequenced. Altogether, 391 conserved miRNA precursor sequences (discovered in the majority of teleost fish species reported in miRbase), eight novel miRNA precursor sequences (so far only discovered in lumpfish), and 443 unique mature miRNAs were identified. Transcriptomics analysis suggested organ-specific and age-specific expression of miRNAs (e.g., miR-122-1-5p specific of the liver). Most of the miRNAs found in lumpfish are conserved in teleost and higher vertebrates, suggesting an essential and common role across teleost and higher vertebrates. This study is the first miRNA characterization of lumpfish that provides the reference miRNAome for future functional studies.

Characterization of MDA5 and microRNA-203 negatively regulates the RLR signaling pathway via targeting MDA5 in miiuy croaker

MDA5 is a member of retinoic acid-inducible gene I (RIG-I)-like receptors (RLR receptors), which may play a crucial role in the immune regulation process. Recently, microRNAs (miRNAs) have been shown to act as an important regulator in the RLRs signaling pathway. Additionally, the MDA5 gene, as a significant cytosolic pathogen recognition receptor (PRR), its characteristics and functions have been extensively investigated, while less research has been done on the mechanisms of MDA5-miRNA mediated gene regulation. In this study, the evolution and functional characterization of MDA5 from miiuy croaker (mmiMDA5) were characterized. Comparative genomic analysis demonstrated that the ascidiacea and superclass do not have the MDA5 gene in the process of evolution. MDA5 contains four structural domains: CARD, ResIII, Helicase C, and RIG-I C-RD. The MDA5 was ubiquitously expressed in all tested miiuy croaker tissues. Moreover, the expressions were significantly up-regulated after stimulation with poly (I: C), which indicated that MDA5 might be involved in the antiviral immune response. The bioinformatics predicted programs have indicated that miR-203 has a direct negative regulatory effect on MDA5 in miiuy croaker. Furthermore, the dual-luciferase reporter assay have showed that miR-203 was the direct negative regulator of MDA5 in miiuy croaker. This study is the first to demonstrate that miRNA can suppress cytokines by regulating the RLR signaling pathway in teleost fish, providing some new ideas for studying miRNA-mediated regulation of immune responses in mammals.

Identification and expression analysis of cobia (Rachycentron canadum) liver-related miRNAs under hypoxia stress

At present, due to the influence of global warming, seasonal change, diurnal variation, and eutrophication of the water body, hypoxia has become one of the major factors limiting the stable development of cobia (Rachycentron canadum) culture. In this study, the miRNAs involved in hypoxia stress were screened, and the target genes of miRNAs were annotated and analyzed. The results showed that a total of 184 conservative microRNA (miRNA) and 121 newly predicted miRNA were obtained by sequencing the liver of control (C) and hypoxic (dissolved oxygen, DO (2.64 ± 0.25) mg/L; 3 h) (S) groups. The pathways involved in energy metabolism included starch and sucrose metabolism (ko00500), glycosaminoglycan degradation (ko00531), and galactose metabolism (ko00052). The results indicate that the body maintains physiological activities by regulating some important pathways at the transcriptional level under hypoxia stress, such as the conversion of aerobic metabolism and anaerobic metabolism, the reduction of energy consumption, and the promotion of red blood cell proliferation to maintain the homeostasis of the body.

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.

Megalocytivirus Induces Complicated Fish Immune Response at Multiple RNA Levels Involving mRNA, miRNA, and circRNA

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

Pol-miR-150 regulates anti-bacterial and viral infection in Japanese flounder (Paralichthys olivaceus) via the lysosomal protein LMP2L

MiR-150 is a microRNA (miRNA) present in a number of teleost species, but its target and regulation mechanism are unknown. Similarly, lysosome membrane protein 2-like (LMP2L) is a gene identified in fish but with unknown function. In this study, we examined the regulation mechanism and function of flounder miR-150 (named pol-miR-150) and its target gene LMP2L (named PoLMP2L) in association with bacterial and viral infection. We found that pol-miR-150 expression was not only modulated by the bacterial pathogen Streptococcus iniae but also by the viral pathogen megalocytivirus. Pol-miR-150 targeted PoLMP2L by binding to the 3'-untranslated region (3'-UTR) of PoLMP2L and inhibited PoLMP2L expression in vitro and in vivo. PoLMP2L is a member of the CD36 superfamily of scavenger receptors and homologous to but phylogenetically distinct from lysosomal integral membrane protein type 2 (LIMP2). PoLMP2L was localized mainly in the lysosomes and expressed in multiple organs of flounder. In vivo knockdown and overexpression of PoLMP2L enhanced and suppressed, respectively, S. iniae dissemination in flounder tissues, whereas in vivo knockdown and overexpression of pol-miR-150 produced the opposite effects on S. iniae dissemination. In addition, pol-miR-150 knockdown also significantly inhibited the replication of megalocytivirus. The results of this study revealed the regulation mechanism and immune functions of fish miR-150 and LMP2L, and indicated that LMP2L and miR-150 play an important role in the antimicrobial immunity of fish.

Molecular Mechanisms Regulating Muscle Plasticity in Fish

Growth rates in fish are largely dependent on genetic and environmental factors, of which the latter can be highly variable throughout development. For this reason, muscle growth in fish is particularly dynamic as muscle structure and function can be altered by environmental conditions, a concept referred to as muscle plasticity. Myogenic regulatory factors (MRFs) like Myogenin, MyoD, and Pax7 control the myogenic mechanisms regulating quiescent muscle cell maintenance, proliferation, and differentiation, critical processes central for muscle plasticity. This review focuses on recent advancements in molecular mechanisms involving microRNAs (miRNAs) and DNA methylation that regulate the expression and activity of MRFs in fish. Findings provide overwhelming support that these mechanisms are significant regulators of muscle plasticity, particularly in response to environmental factors like temperature and nutritional challenges. Genetic variation in DNA methylation and miRNA expression also correlate with variation in body weight and growth, suggesting that genetic markers related to these mechanisms may be useful for genomic selection strategies. Collectively, this knowledge improves the understanding of mechanisms regulating muscle plasticity and can contribute to the development of husbandry and breeding strategies that improve growth performance and the ability of the fish to respond to environmental challenges.

Identification and Characterization of lncRNAs Related to the Muscle Growth and Development of Japanese Flounder (Paralichthys olivaceus)

Long noncoding RNAs (lncRNAs) play an important role in many life activities, but the expression pattern and function of lncRNAs in Japanese flounder skeletal muscle are unclear. In this study, 751 lncRNAs were selected from skeletal muscle in different development stages of the Japanese flounder [stage A: larval 7 days post hatching (dph); stage B: juvenile about 90 dph; stage C (female) and stage D (male): adult about 24 months] using coding potential analysis methods. In total, 232, 211, 194, 28, 29, and 14 differentially expressed lncRNAs and 9549, 8673, 9181, 1821, 1080, and 557 differentially expressed mRNAs were identified in comparisons of A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. We identified the cis- and trans-regulatory target genes of differentially expressed lncRNAs, and lncRNA-gene interaction networks were constructed using the Cytoscape program. In total, there were 200, 200, 200, 93, 47, and 11 cis-regulation relationships, and 29, 19, 24, 38, 8, and 47 trans-regulation relationships in the comparisons between A versus B, A versus C, A versus D, B versus C, B versus D, and C versus D, respectively. These results indicate that lncRNA may participate in the development of Japanese flounder skeletal muscle through cis- or trans-acting mechanisms, thus providing a scientific basis for further study of the biological function of lncRNA in Japanese flounder skeletal muscle. Based on these relationships, functional annotation of the related lncRNAs was performed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Differentially expressed genes associated with muscle development were enriched in multiple pairs of comparisons (e.g., differentially expressed genes LOC109640370, LOC109634180, LOC109643555, rusc1, and LOC109626999 were enriched in the actin-binding term, and differentially expressed genes LOC109640370, was, LOC109644970, LOC109643555, and itga9 were enriched in the regulation of the actin cytoskeleton pathway in the KEGG pathway analysis in the comparison of stages C and D). We predicted lncRNA-mRNA, miRNA-mRNA, and lncRNA-miRNA regulatory relationships and constructed interactive networks using Cytoscape software. Co-expression networks show that most lncRNAs interact with one or two predicted miRNAs involved in muscle growth and development. These results provide a basis for further study of the function of lncRNAs on skeletal muscle in different developmental stages of Japanese flounder.

Distinct expression patterns of seven crucial microRNAs during early embryonic development in medaka (Oryzias latipes)

MicroRNAs (i.e. miRNAs) are small non-coding RNAs that play essential modulation roles in embryonic development in vertebrates. Paternal and maternal miRNAs contribute to the development of post-fertilization embryo and zygotic genome activation. The pattern of expression and their roles in embryonic development of medaka are not clearly understood. The present study, therefore, examined a temporal expression of seven miRNAs, ola-let-7a, ola-miR-202-3p, ola-miR-126-3p, ola-miR-122, ola-miR-92a, ola-miR-125a-3p and ola-miR-430a in sperm, oocytes, and embryos during early developmental stages. Three unique expression patterns of miRNAs were observed. ola-let7a, ola-miR-202-3p and ola-miR-126-3p showed both paternal and maternal expression, and ola-miR-122, ola-miR-92a, ola-miR-125a-3p showed maternal expression only. The expression of six out of seven miRNAs significantly decreased after maternal-zygotic transition (MZT), whereas ola-miR-430a expression initiated only after MZT. The temporal dynamic expression of these miRNAs suggests their potential roles in early embryogenesis and genome-zygotic activation in medaka.

Gonad-Specific Transcriptomes Reveal Differential Expression of Gene and miRNA Between Male and Female of the Discus Fish (Symphysodon aequifasciatus)

The discus fish (Symphysodon aequifasciatus) is an ornamental fish that is well-known around the world. In artificial reproduction, they must be matched by one male and one female, whereas phenotype investigation indicated that there are no significant differences in appearance between males and females, which causes great difficulties in the mating during artificial reproduction. So, it is of great importance to establish artificial sex identification methods for the discus fish. The molecular mechanism of the sexual dimorphism of the discus fish was previously unknown. In this study, we constructed six cDNA libraries from three adult testes and three adult ovaries and performed RNA sequencing for identifying sex-biased candidate genes and microRNAs (miRNAs). A total of 50,082 non-redundant genes (unigenes) were identified, of which 18,570 unigenes were significantly overexpressed in testes, and 11,182 unigenes were significantly overexpressed in ovaries. A total of 551 miRNAs were identified, of which 47 miRNAs were differentially expressed between testes and ovaries. Eight differentially expressed unigenes, seven differentially expressed miRNAs and one non-differential miRNA were validated by quantitative real-time polymerase chain reaction. Twenty-four of these differentially expressed miRNAs and their 15 predicted target genes constituted 41 miRNA–mRNA interaction pairs, and some of vital sex-related metabolic pathways were also identified. These results revealed these differentially expressed genes and miRNAs between ovary and testis might be involved in regulating gonadal development, sex determination, gametogenesis, and physiological function maintenance, and there are complex regulatory networks between genes and miRNAs. It can help us understand the molecular mechanism of the sexual dimorphism and obtain a high-efficiency sex identification method in the artificial reproduction process of the discus fish.

miR-148 targets CiGadd45ba and CiGadd45bb to modulate the inflammatory response to bacterial infection in grass carp

In grass carp (Ctenopharyngodon idella), septicemia is a systemic inflammatory response to bacterial infection. Once infected bacteria, a hyper-inflammatory state that could lead to septic shock and death. There is increasing evidence that microRNAs are involved in the regulation of the inflammatory response. Ctenopharyngodon idella growth arrest and DNA damage-inducible 45 ba and bb (CiGadd45ba and CiGadd45bb) are two subtypes of Gadd45b. In the present study, miR-148 was confirmed to be involved in the inflammatory response after infection with Aeromonas hydrophila. Dual-luciferase reporter assays and miRNA expression profiling confirmed that miR-148 targeted both CiGadd45ba and CiGadd45bb. Transfection with miR148 mimics and inhibitors altered the expression levels of proinflammatory genes, suggesting that miRNAs regulate the immune response in grass carp. Our results provide a theoretical basis for studying the molecular mechanism underlying the regulation of inflammation by miR-148 in grass carp.

miR-430a regulates the development of left-right asymmetry by targeting sqt in the teleost

microRNAs (miRNAs) are short, endogenous non-coding RNAs that contain approximately 18-22 nucleotides. miRNAs are involved in gene regulation by recognizing and binding the 3'UTR of target gene. In our previous data, miR-430 family showed significant differential expression modes through metamorphosis in Japanese flounder. It was speculated that miR-430a plays a key role in left-right patterning. We predicted the targets of miR-430a and gene ontology (GO) was performed. We speculated miR-430a is involved in the basal molecular function and organ development. In Japanese flounder, sqt as a target of miR-430a was enriched into heart development term. Sqt has been reported to participate in mesendoderm formation and organ development. Cardiac morphogenesis is the first asymmetric development process, which breaks left-right symmetry in bilateria. It was used as a marker to detect L-R asymmetric effects of miR-430a. Overexpression and suppression of miR-430a resulted in abnormal KV (Kupffer's vesicles) development and disordered in nodal-related expression with consequent cardiac laterality. Squint mRNA of Japanese flounder (Posqt) as a target of miR-430a was overexpressed and caused similar phenotype with miR-430a suppression group, such as longer cilia in KV and high range of clmc2 and spaw ectopic expression. Moreover, rescue experiments were performed and suggested that cardiac and KV defections, induced by overexpressing miR-430a, could be rescued by injecting Posqt mRNA. These results suggested that miR-430a regulates the development of left-right asymmetry by targeting sqt in the teleost.

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.

The role of miR-92 in regulating early development and metamorphosis of Japanese flounder Paralichthys olivaceus

MicroRNAs are a class of short non-coding RNAs that contain approximately 22 nucleotides and play a regulatory role in RNA silencing and translational repression. miR-92 belongs to the miR-17-92 family and has a regulatory effect on cell proliferation, apoptosis, and expression of proto-oncogenes and tumor suppressor genes. However, its function in flatfish is unclear. In this study, we used farmed Japanese flounder, Paralichthys olivaceus, and showed that gata5 is a target gene of miR-92. Experiments on miR-92 overexpression indicated that gata5 and sox17 were downregulated, while the transcription level of ntl increased. By contrast, depletion of miR-92 resulted in increased gata5 and sox17 levels and reduced ntl level. Moreover, thiourea treatment indicated that miR-92 may inhibit the metamorphic development of Japanese flounder. Our study suggests that miR-92 regulates the fate of endoderm and mesoderm by controlling gata5.

Genetic features of Haliotis discus hannai by infection of vibrio and virus

BACKGROUND

Haliotis discus hannai more commonly referred to as the Pacific Abalone is of significant commercial and economical value in South Korea, with it being the second largest producer in the world. Despite this significance there is a lack of genetic studies with regards to the species. Most existing studies focused mainly on environmental factors.

OBJECTIVE

To provide a comprehensive review describing the genetic feature of Haliotis discus hannai by infection of vibrio and virus.

METHODS

This review summarized the immune response in the Haliotis spp. with regards to immunological genes such as Cathepsin B, C-type lectin and Toll-like receptors. Genetic studies with regards to transposable elements and miRNAs are few and far between. A study identified LTR retrotransposon Ty3/gypsy in the species. As to miRNA, a single study identified numerous miRNAs in the Haliotis discus hannai.

CONCLUSION

This paper sought to provide an overview of genetic perspective with regards to immune response genes, transposable elements and miRNAs.

miR-21 targets jnk and ccr7 to modulate the inflammatory response of grass carp following bacterial infection

Grass carp septicemia is a systemic inflammatory response that develops following a bacterial infection. The hyperinflammatory state develops could lead to septic shock and lethality. There is increasing evidence that microRNAs are involved in the regulation of the inflammatory response. In the present study, miR-21 was confirmed to be involved in the inflammatory response following infection with Aeromonas hydrophila and LPS stimulation. Both jnk and ccr7 were identified as target gene of miR-21 by overexpression, inhibition, and dual luciferase reporter assays experiments. Meanwhile, miR-21 targets the jnk and ccr7 to modulate downstream pro-inflammatory factors tnf-α, il-1β, il-6, and il-12. Our results provide a theoretical basis for exploring the molecular mechanism of grass carp miR-21 regulating inflammation.

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.

Comparative expression analysis of let-7 microRNAs during ovary development in Megalobrama amblycephala

As a critical regulator of gene expression, let-7 family miRNAs have been reported to be involved in multiple physiological processes. In this study, in order to elucidate the putative regulatory effect of let-7 miRNAs during fish gonadal development and to identify which member is crucial for this regulation, the expression of ten members including let-7a/b/c/d/e/f/g/h/i/j were quantified in ovary, pituitary, and brain tissues during the different ovarian developmental stages of blunt snout bream Megalobrama amblycephala. According to the data from analysis of expression patterns, let-7a showed the highest expression value in almost all the tested ovaries, pituitaries, and brains, with let-7b and let-7d moderately expressed, following by other let-7 miRNAs. In terms of the differential expression levels of ten let-7 miRNAs at each developmental stage, the results showed that let-7a/b/d/f/h expression gradually increased during the ovary development from stage I to V and dropped significantly at the phase VI in ovary tissues. However, the expression of let-7a/b/e/f in pituitary increased during the ovary development from stage I to IV and declined at stage V. Among all the let-7 miRNAs, let-7a/b/d had the highest expression and their expression patterns were consistent with the gonad development of M. amblycephala. Furthermore, the mostly predicted target genes of let-7 miRNAs are significantly enriched in signaling pathways closely related to gonadal development through KEGG enrichment analysis. These results indicate that let-7 miRNA members, especially let-7a/b/d, may play important roles in the regulation of ovary development in M. amblycephala through negatively regulating expression of their target genes.

Identification and Characterization of MicroRNAs in Skin of Chinese Giant Salamander (Andrias davidianus) by the Deep Sequencing Approach

MicroRNAs (miRNA) play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. However, the miRNAs from skin of Andrias davidianus have not been reported. In this study, a small-RNA cDNA library was constructed and sequenced from skin of A. davidianus. A total of 513 conserved miRNAs belonging to 174 families were identified. The remaining 108 miRNAs we identified were novel and likely to be skin tissue-specific but were expressed at low levels. The presence of randomly selected 15 miRNAs identified and their expression in eight different tissues from A. davidianus were validated by stem-loop qRT-PCR. For better understanding the functions of miRNAs, 129,791 predicated target genes were analyzed by GO and their pathways illustrated by KEGG pathway analyses. The results show that these identified miRNAs from A. davidianus skin are involved in a broad range of physiological functions including metabolism, growth, development, and immune responses. This study exhaustively identifies miRNAs and their target genes, which will ultimately pave the way for understanding their role in skin of A. davidianus and other amphibians. Further studies are necessary to better understand miRNA-mediated gene regulation.

Identification and expression analysis of miRNA in hybrid snakehead by deep sequencing approach and their targets prediction

MicroRNAs (miRNAs) play important regulatory roles in numerous biological processes, but there is no report on miRNAs of hybrid snakehead. In this study, four independent small RNA libraries were constructed from the spleen, liver kidney and muscle of hybrid snakehead. These libraries were sequenced using deep sequencing technology, as result, a total of 1,067,172, 1,152,002, 1,653,941 and 970,866 clean reads from these four libraries were obtained. 252 known miRNAs and 63 putative novel miRNAs in these small RNA dataset were identified. The stem-loop RT-qPCR analysis indicated that eight known miRNAs and two novel miRNAs show different expression in eight different kinds of tissues. For better understanding the functions of miRNAs, 95,947 predicated target genes were analyzed by GO and their pathways, the results indicated that these targets of the identified miRNAs are involved in a broad range of physiological functions.

microRNA-210 participates in regulating RIG-I signaling pathway via targeting DUBA in miiuy croaker after poly(I:C) stimulation

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that participate in the regulation of various biological processes. A series of microRNAs have been shown to be important regulators of both innate and adaptive immune responses, including RIG-I signaling pathway. In this study, we evaluated the regulation role of miR-210 in the RLRs signaling pathway of miiuy croaker. Upon poly(I:C) stimulation, the expression of miR-210 in both miiuy croaker spleen tissues and macrophages were significantly upregulated. By means of the dual luciferase reporter assay, a direct interaction between miR-210 and the 3-untranslated region (UTR) of Deubiquitinating enzyme A (DUBA) was confirmed, and we found that miR-210 could reduce the luciferase levels of wild-type DUBA 3'UTR, whereas mutant-type led to a complete abrogation of the negative effect. Furthermore, the negative regulatory effects of pre-miR-210 on DUBA have been indicated in a dose- and time-dependent manners. As DUBA is an important regulator involved in the RLRs signaling pathway and could bind with and regulate TRAF3, we also examined the expression patterns of DUBA and TRAF3 in vivo and in vitro. We found that the expression of both DUBA and TRAF3 were significantly changed upon poly(I:C) stimulation in miiuy croaker. The expression patterns between miR-210 and DUBA showed a negative correlation, which indicated that miR-210 can target and downregulate the expression of DUBA. Overall, these results will enrich the knowledge of immune response related miRNAs in miiuy croaker, which will be useful for better understanding the complicated regulatory networks in fish species.

The DNA methylation status of MyoD and IGF-I genes are correlated with muscle growth during different developmental stages of Japanese flounder (Paralichthys olivaceus)

Many genes related to muscle growth modulate myoblast proliferation and differentiation and promote muscle hypertrophy. MyoD is a myogenic determinant that contributes to myoblast determination, and insulin-like growth factor 1 (IGF-I) interacts with MyoD to regulate muscle hypertrophy and muscle mass. In this study, we aimed to assess DNA methylation and mRNA expression patterns of MyoD and IGF-I during different developmental stages of Japanese flounder, and to examine the relationship between MyoD and IGF-I gene. DNA and RNA were extracted from muscles, and DNA methylation of MyoD and IGF-I promoter and exons was detected by bisulfite sequencing. The relative expression of MyoD and IGF-I was measured by quantitative polymerase chain reaction. IGF-I was measured by radioimmunoassay. Interestingly, the lowest expression of MyoD and IGF-I emerged at larva stage, and the mRNA expression was negatively associated with methylation. We hypothesized that many skeletal muscle were required to complete metamorphosis; thus, the expression levels of MyoD and IGF-I genes increased from larva stage and then decreased. The relative expression levels of MyoD and IGF-I exhibited similar patterns, suggesting that MyoD and IGF-I regulated muscle growth through combined effects. Changes in the concentrations of IGF-I hormone were similar to those of IGF-I gene expression. Our results the mechanism through which MyoD and IGF-I regulate muscle development and demonstrated that MyoD interacted with IGF-I to regulate muscle growth during different developmental stages.

Discovery and characterization of conserved and novel microRNAs from blunt snout bream (Megalobrama amblycephala) by deep sequencing

MicroRNAs (miRNAs) are short, single stranded RNA molecules with approximately 22 nts in length, which regulate the stability and translation of messenger RNAs in several organisms. To increase the repertoire of miRNAs characterized in M. amblycephala, we used the deep sequencing technology to sequence a small RNA library using pooled RNA sample isolated from the 4 different tissues of M. amblycephala. A total of 309 conserved miRNAs that originated from 131 miRNA families were detected. 15 novel candidates miRNA were identified. Randomly selected 6 miRNAs were analyzed by stem-loop qRT-PCR and differential expression patterns were observed in 6 different tissues of M. amblycephala. Furthermore, the potential targets were predicted. GO analysis showed that most of the targets were involved in a broad range of physiological functions including fish growth, development, metabolism, stress responses and so on. Overall, our results significantly increased the number of novel miRNAs in M. amblycephala, which should be useful for further investigation into the role of miRNAs in regulating diverse biological processes.

Inducible microRNA-122 modulates RIG-I signaling pathway via targeting DAK in miiuy croaker after poly(I:C) stimulation

MicroRNA-122 (miR-122) was originally identified in mouse and then lots of researches on miR-122 had been performed in mammals. However, the functional study of miR-122 were restricted in fish. In miiuy croaker, miR-122 is sensitive to poly(I:C) stimulation. In this study, a combination of bioinformatics and experimental techniques were used to investigate the functions of miR-122. DAK is a putative target gene of miR-122 which was predicted by bioinformatics, and further the luciferase reporter assays were used to confirm the target sites in DAK 3'untranslated region. The inhibiting effect of miR-122 mimics or pre-miR-122 on DAK presented the dose and time dependent manners, and the pre-miR-122 showed stronger inhibiting effect on DAK than the miR-122 mimics. Therefore, the miR-122 participate in regulating RIG-I-like receptors signaling pathway via inhibiting DAK which is the inhibitors of MDA5. The expression of miR-122 and DAK showed negative relationship in both miiuy croaker spleen and macrophages, which imply that miR-122 may regulate DAK at the post-transcriptional level. These results will enhance our understanding about the regulation of miRNAs on immune response in fish.

miR-1338-5p Modulates Growth Hormone Secretion and Glucose Utilization by Regulating ghitm in Genetically Improved Farmed Tilapia (GIFT, Oreochromis niloticus)

MicroRNAs (miRNAs) are endogenous, non-coding small RNA molecules about 22 nt in length, which could regulate the expressions of target genes and participate in growth and development of organisms. Genetically improved farmed tilapia (GIFT, Oreochromis niloticus) is an important economic freshwater species in China and the growth performance is one of the main breeding indicators. Growth hormone inducible transmembrane protein (ghitm) plays an important role in growth and development of both mammals and invertebrates; however, little studies have been reported on fish. Our previous experiments indicated that miR-1338-5p expression may be negatively correlated with ghitm expression. In this study, we firstly used qRT-PCR and northern blot to verify the expression of miR-1338-5p and ghitm, and determined the binding site of miR-1338-5p in the ghitm 3'-untranslated region (UTR) by luciferase reporter assay. Secondly, juveniles GIFT injected with miR-1338-5p antagomir were used to analyze the regulatory function of the miR-1338-5p-ghitm pair in vivo. The results showed that the ghitm 3'-UTR was complementary to the 5' 2-8-nt site of miR-1338-5p. Inhibition of miR-1338-5p promoted ghitm expression in the pituitary and liver of GIFT. ghitm could interfere in the growth hormone (Gh)-growth hormone receptor (Ghr)-insulin-like growth factor (Igf) signaling pathway by competing with the ghr1 for combination with Gh, and then reduce the growth of GIFT. Moreover, the reduction of Gh in serum may regulate insulin secretion and result in the increasing sugar and fat storage in serum and liver. Our results suggest that miR-1338-5p participates in the growth and development of GIFT through the regulation of ghitm, which provides theoretical support for the study of the fish growth mechanism.

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.

The integrated analysis of RNA-seq and microRNA-seq depicts miRNA-mRNA networks involved in Japanese flounder (Paralichthys olivaceus) albinism

Albinism, a phenomenon characterized by pigmentation deficiency on the ocular side of Japanese flounder (Paralichthys olivaceus), has caused significant damage. Limited mRNA and microRNA (miRNA) information is available on fish pigmentation deficiency. In this study, a high-throughput sequencing strategy was employed to identify the mRNA and miRNAs involved in P. olivaceus albinism. Based on P. olivaceus genome, RNA-seq identified 21,787 know genes and 711 new genes by transcripts assembly. Of those, 235 genes exhibited significantly different expression pattern (fold change ≥2 or ≤0.5 and q-value≤0.05), including 194 down-regulated genes and 41 up-regulated genes in albino versus normally pigmented individuals. These genes were enriched to 81 GO terms and 9 KEGG pathways (p≤0.05). Among those, the pigmentation related pathways-Melanogenesis and tyrosine metabolism were contained. High-throughput miRNA sequencing identified a total of 475 miRNAs, including 64 novel miRNAs. Furthermore, 33 differentially expressed miRNAs containing 13 up-regulated and 20 down-regulated miRNAs were identified in albino versus normally pigmented individuals (fold change ≥1.5 or ≤0.67 and p≤0.05). The next target prediction discovered a variety of putative target genes, of which, 134 genes including Tyrosinase (TYR), Tyrosinase-related protein 1 (TYRP1), Microphthalmia-associated transcription factor (MITF) were overlapped with differentially expressed genes derived from RNA-seq. These target genes were significantly enriched to 254 GO terms and 103 KEGG pathways (p<0.001). Of those, tyrosine metabolism, lysosomes, phototransduction pathways, etc., attracted considerable attention due to their involvement in regulating skin pigmentation. Expression patterns of differentially expressed mRNA and miRNAs were validated in 10 mRNA and 10 miRNAs by qRT-PCR. With high-throughput mRNA and miRNA sequencing and analysis, a series of interested mRNA and miRNAs involved in fish pigmentation are identified. And the miRNA-mRNA regulatory network also provides a solid starting point for further elucidation of fish pigmentation deficiency.

Evolutionary conservation and expression of miR-10a-3p in olive flounder and rock bream

MicroRNAs (miRNAs) are small non-coding RNAs (ncRNAs) that mainly bind to the seed sequences located within the 3' untranslated region (3' UTR) of target genes. They perform an important biological function as regulators of gene expression. Different genes can be regulated by the same miRNA, whilst different miRNAs can be regulated by the same genes. Here, the evolutionary conservation and expression pattern of miR-10a-3p in olive flounder and rock bream was examined. Binding sites (AAAUUC) to seed region of the 3' UTR of target genes were highly conserved in various species. The expression pattern of miR-10a-3p was ubiquitous in the examined tissues, whilst its expression level was decreased in gill tissues infected by viral hemorrhagic septicemia virus (VHSV) compared to the normal control. In the case of rock bream, the spleen, kidney, and liver tissues showed dominant expression levels of miR-10a-3p. Only the liver tissues in the rock bream samples infected by the iridovirus indicated a dominant miR-10a-3p expression. The gene ontology (GO) analysis of predicted target genes for miR-10a-3p revealed that multiple genes are related to binding activity, catalytic activity, cell components as well as cellular and metabolic process. Overall the results imply that the miR-10a-3p could be used as a biomarker to detect VHSV infection in olive flounder and iridovirus infection in rock bream. In addition, the data provides fundamental information for further study of the complex interaction between miR-10a-3p and gene expression.

Identification and profiling of growth-related microRNAs in Chinese perch (Siniperca chuatsi)

Background

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that play important roles in the regulation of diverse biological processes in eukaryotes. Chinese perch (Siniperca chuatsi) is one of the most economically important fish species widely cultured in China. Growth is an extremely important characteristic in fish. Individual differences in body size are common in Siniperca chuatsi, which significantly influence the aquaculture production of Siniperca chuatsi. However, the underline growth-related regulatory factors, such as miRNAs, are still unknown.

Results

To investigate the growth-related miRNAs in Siniperca chuatsi, two RNA libraries from four growth-related tissues (brain, pituitary, liver, and muscle) of Siniperca chuatsi at 6-month stage with relatively high or low growth rates (big-size group or small-size group) were obtained and sequenced using Solexa sequencing. A total of 252 known miRNAs and 12 novel miRNAs were identified. The expression patterns of these miRNAs in big-size group and small-size group were compared, and the results showed that 31 known and 5 novel miRNAs were differently expressed (DE). Furthermore, to verify the Solexa sequencing, five DE miRNAs were randomly selected and quantified by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The results showed that their expression patterns were consistent with those of Solexa sequencing. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of target genes of DE miRNAs was performed. It showed that the target genes were involved in multiple biological processes including metabolic process, suggesting that metabolic process played an important role in growth of fish.

Conclusions

Siniperca chuatsi is a popular and economically important species in aquaculture. In this study, miRNAs in Siniperca chuatsi with different growth rates were identified, and their expression profiles were compared. The data provides the first large-scale miRNA profiles related to growth of Siniperca chuatsi, which is expected to contribute to a better understanding of the role of miRNAs in regulating the biological processes of growth and possibly useful for Siniperca chuatsi breeding.

Electronic supplementary material

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

Identification of Differentially Expressed Micrornas Associate with Glucose Metabolism in Different Organs of Blunt Snout Bream (Megalobrama amblycephala)

Blunt snout bream (Megalobrama amblycephala) is a widely favored herbivorous fish species and is a frequentlyused fish model for studying the metabolism physiology. This study aimed to provide a comprehensive illustration of the mechanisms of a high-starch diet (HSD) induced lipid metabolic disorder by identifying microRNAs (miRNAs) controlled pathways in glucose and lipid metabolism in fish using high-throughput sequencing technologies. Small RNA libraries derived from intestines, livers, and brains of HSD and normal-starch diet (NSD) treated M. amblycephala were sequenced and 79, 124 and 77 differentially expressed miRNAs (DEMs) in intestines, livers, and brains of HSD treated fish were identified, respectively. Bioinformatics analyses showed that these DEMs targeted hundreds of predicted genes were enriched into metabolic pathways and biosynthetic processes, including peroxisome proliferator-activated receptor (PPAR), glycolysis/gluconeogenesis, and insulin signaling pathway. These analyses confirmed that miRNAs play crucial roles in glucose and lipid metabolism related to high wheat starch treatment. These results provide information on further investigation of a DEM-related mechanism dysregulated by a high carbohydrate diet.

Evolution of Fish Let-7 MicroRNAs and Their Expression Correlated to Growth Development in Blunt Snout Bream

The lethal-7 (let-7) miRNA, known as one of the first founding miRNAs, is present in multiple copies in a genome and has diverse functions in animals. In this study, comparative genomic analysis of let-7 miRNAs members in fish species indicated that let-7 miRNA is a sequence conserved family in fish, while different species have the variable gene copy numbers. Among the ten members including let-7a/b/c/d/e/f/g/h/i/j, the let-7a precursor sequence was more similar to ancestral sequences, whereas other let-7 miRNA members were separate from the late differentiation of let-7a. The mostly predicted target genes of let-7 miRNAs are involved in biological process, especially developmental process and growth through Gene Ontology (GO) enrichment analysis. In order to identify the possible different functions of these ten miRNAs in fish growth development, their expression levels were quantified in adult males and females of Megalobrama amblycephala, as well as in 3-, 6-, and 12-months-old individuals with relatively slow- and fast-growth rates. These ten miRNAs had similar tissue expression patterns between males and females, with higher expression levels in the brain and pituitary than that in other tissues (p < 0.05). Among these miRNAs, the relative expression level of let-7a was the highest among almost all the tested tissues, followed by let-7b, let-7d and let-7c/e/f/g/h/i/j. As to the groups with different growth rates, the expression levels of let-7 miRNAs in pituitary and brain from the slow-growth group were always significantly higher than that in the fast-growth group (p < 0.05). These results suggest that let-7 miRNA members could play an important role in the regulation of growth development in M. amblycephala through negatively regulating expression of their target genes.

microRNA-145 regulates the RLR signaling pathway in miiuy croaker after poly(I:C) stimulation via targeting MDA5

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that participate in diverse biological processes via degrading the target mRNAs or repressing translation. In this study, the regulation of miRNA to the RLR (RIG-I-like receptor) signaling pathway by degrading the target mRNAs was researched in miiuy croaker. MDA5, a microRNA-145-5p (miR-145-5p) putative target gene, was predicted by bioinformatics, and the target sites from the 3'untranslated region of MDA5 transcripts were confirmed using luciferase reporter assays. Pre-miR-145 was more effective in inhibiting MDA5 than miR-145-5p mimic, and the effect was dose- and time-dependent. The expression patterns of miR-145-5p and MDA5 were analyzed in liver and kidney from miiuy croaker. Results implied that miR-145-5p may function via degrading the MDA5 mRNAs, thereby regulating the RLR signaling pathway. Studies on miR-145-5p will enrich knowledge of its functions in immune response regulation in fish, as well as offer a basis for regulatory networks that are composed of numerous miRNAs.

Integrated analysis neurimmiRs of tilapia (Oreochromis niloticus) involved in immune response to Streptococcus agalactiae, a pathogen causing meningoencephalitis in teleosts

MicroRNAs (miRNAs) are a class of noncoding RNA molecules and play important roles in a wide spectrum of biological processes, including in immune response. Recent years have witnessed considerable amount of research interest in studies on miRNA-mediated modulation gene function during neuroinflammation. Here, we evaluated Streptococcus agalactiae infected tilapia (Oreochromis niloticus) brain for the expression profile of miRNAs, potential functions and their correlation with genes involved in inflammatory pathways. A total of 1981 miRNAs were identified, including in 486 miRNAs which have homologues in the currently available databases and 1945 novel miRNAs. The expression levels of 547 miRNAs were significantly altered at 6 h-48 h post-bacterial infection, and these miRNAs were therefore classified as differentially expressed tilapia miRNAs. Real-time PCR were implemented for 14 miRNAs co-expressed in five samples, and agreement was confirmed between the high-throughput sequencing and real-time PCR data. For the 486 differentially expressed miRNAs target 41,820 genes. GO and KEGG enrichment analysis revealed that some target genes of miRNAs were grouped mainly into the categories of apoptotic, signal pathwayand immune response. This is the first report of comprehensive identification of teleost miRNAs being differentially regulated in brain in normal conditions relating to bacterial infection.

Gene expression patterns regulating embryogenesis based on the integrated de novo transcriptome assembly of the Japanese flounder

The Japanese flounder (Paralichthys olivaceus) is one of the most important commercial and biological marine fishes. However, the molecular biology involved during embryogenesis and early development of the Japanese flounder remains largely unknown due to a lack of genomic resources. A comprehensive and integrated transcriptome is necessary to study the molecular mechanisms of early development and to allow for the detailed characterization of gene expression patterns during embryogenesis; this approach is critical to understanding the processes that occur prior to mesectoderm formation during early embryonic development. In this study, more than 117.8 million 100bp PE reads were generated from pooled RNA extracted from unfertilized eggs to 41dph (days post-hatching) embryos and were sequenced using Illumina pair-end sequencing technology. In total, 121,513 transcripts (≥200bp) were obtained using de novo assembly. A sequence similarity search indicated that 52,338 transcripts show significant similarity to 22,462 known proteins from the NCBI non-redundant database and the Swiss-Prot protein database and were annotated using Blast2GO. GO terms were assigned to 44,627 transcripts with 12,006 functional terms, and 10,024 transcripts were assigned to 133 KEGG pathways. Furthermore, gene expression differences between the unfertilized egg and the gastrula embryo were analysed using Illumina RNA-Seq with single-read sequencing technology, and 24,837 differentially and specifically expressed transcripts were identified and included 5,286 annotated transcripts and 19,569 non-annotated transcripts. All of the expressed transcripts in the unfertilized egg and gastrula embryo were further classified as maternal, zygotic, or maternal-zygotic transcripts, which may help us to understand the roles of these transcripts during the embryonic development of the Japanese flounder. Thus, the results will contribute to an improved understanding of the gene expression patterns and signalling pathways that control the molecular mechanisms of early embryonic development.

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.

Dynamic mRNA and miRNA expression analysis in response to intermuscular bone development of blunt snout bream (Megalobrama amblycephala)

Intermuscular bone (IB), which occurs only in the myosepta of lower teleosts, is attracting more attention because they are difficult to remove and make the fish unpleasant to eat. By gaining a better understanding of the genetic regulation of IB development, an integrated analysis of miRNAs and mRNAs expression profiling was performed on Megalobrama amblycephala. Four key development stages were selected for transcriptome and small RNA sequencing. A number of significantly differentially expressed miRNAs/genes associated with bone formation and differentiation were identified and the functional characteristics of these miRNAs/genes were revealed by GO function and KEGG pathway analysis. These were involved in TGF-β, ERK and osteoclast differentiation pathways known in the literature to affect bone formation and differentiation. MiRNA-mRNA interaction pairs were detected from comparison of expression between different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating bone development and differentiation. A negative regulation effect of two miRNAs was verified through dual luciferase reporter assay. As a unique public resource for gene expression and regulation during the IB development, this study is expected to provide forwards ideas and resources for further biological researches to understand the IBs’ development.

The miiuy croaker microRNA transcriptome and microRNA regulation of RIG-I like receptor signaling pathway after poly(I:C) stimulation

MicroRNAs (miRNAs) as endogenous small non-coding RNAs play key regulatory roles in diverse biological processes via degrading the target mRNAs or inhibiting protein translation. Previously many researchers have reported the identification, characteristic of miRNAs and the interaction with its target gene. But, the study on the regulation of miRNAs to biological processes via regulatory the key signaling pathway was still limited. In order to comprehend the regulatory mechanism of miRNAs, two small RNA libraries from the spleen of miiuy croaker individuals with or without poly(I:C) infection were constructed. The 197 conserved miRNAs and 75 novel miRNAs were identified, and 14 conserved and 8 novel miRNAs appeared significant variations. Those differently expressed miRNAs relate to immune regulation of miiuy croaker. Furthermore, expressions of four differently expressed miRNAs were validated by qRT-PCR, and the result was consistent with sequencing data. The target genes of the differently expressed miRNAs in the two libraries were predicted, and some candidate target genes were involved in the RIG-I-like receptor (RLR) signaling pathway. The negative regulation of miRNAs to target genes were confirmed by comparing the expression pattern of miRNAs and their target genes. The results of regulating target genes were that firstly directly or indirectly activating the downstream signaling cascades and subsequent inducting the type I interferon, inflammatory cytokines and apoptosis. These studies could help us to deeper understand the roles of miRNAs played in the fish immune system, and provide a new way to investigate the defense mechanism of fish.

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.

MicroRNAs of the mesothorax in Qinlingacris elaeodes, an alpine grasshopper showing a wing polymorphism with unilateral wing form

MicroRNAs (miRNAs) are now recognized as key post-transcriptional regulators in regulation of phenotypic diversity. Qinlingacris elaeodes is a species of the alpine grasshopper, which is endemic to China. Adult individuals have three wing forms: wingless, unilateral-winged and short-winged. This is an ideal species to investigate the phenotypic plasticity, development and evolution of insect wings because of its case of unilateral wing form in both the sexes. We sequenced a small RNA library prepared from mesothoraxes of the adult grasshoppers using the Illumina deep sequencing technology. Approximately 12,792,458 raw reads were generated, of which the 854,580 high-quality reads were used only for miRNA identification. In this study, we identified 49 conserved miRNAs belonging to 41 families and 69 species-specific miRNAs. Moreover, seven miRNA*s were detected both for conserved miRNAs and species-specific miRNAs, which were supported by hairpin forming precursors based on polymerase chain reaction. This is the first description of miRNAs in alpine grasshoppers. The results provide a useful resource for further studies on molecular regulation and evolution of miRNAs in grasshoppers. These findings not only enrich the miRNAs for insects but also lay the groundwork for the study of post-transcriptional regulation of wing forms.