Discovery of miRNAs and Their Corresponding miRNA Genes in Atlantic Cod (Gadus morhua): Use of Stable miRNAs as Reference Genes Reveals Subgroups of miRNAs That Are Highly Expressed in Particular Organs

miRNA-seq analysis of liver tissue from largemouth bass (Micropterus salmoides) in response to oxytetracycline and enzyme-treated soy protein

The specific miRNA regulation triggered by enzyme-treated soybean protein in response to well-known stressors, such as the prophylactic use of the antimicrobial oxytetracycline, remains unknown. Hence, this study aimed to evaluate the regulatory changes of hepatic miRNAs induced by oxytetracycline and enzyme-treated soybean protein in largemouth bass dietary formulations. The experiment was designed with three groups: the normal control (NC), the oxytetracycline exposure treatment group (OTC), and the pre-treatment with enzyme-treated soybean protein before oxytetracycline exposure group (ETSP). miRNA sequencing was employed to characterize the differences between these groups. In conclusion, the NC group exhibited up-regulation of 13 host miRNAs and down-regulation of 1 miRNA compared to the OTC group, whereas the ETSP group showed an increasing trend of 36 host miRNAs and a decreasing trend of 13 host miRNAs compared to the OTC group. Nine miRNAs were identified as prudential targets for enzyme-treated soy protein, protecting the largemouth bass liver from oxytetracycline. Furthermore, gene ontology analysis revealed nine key miRNAs that mediate signaling pathways with significant differences. The cellular lipid metabolic process was identified as the most important biological process, and the propanoate metabolism pathway was highlighted as significant. These results will facilitate further exploration of the mechanism by which enzyme-treated soy protein alleviates the effects of oxytetracycline on largemouth bass in water environments.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

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.

LncRNA-adm2 targets adm2 via cid-miR-n3 and negatively regulates the inflammatory response in grass carp (Ctenopharyngodon idella)

Long non-coding RNAs (lncRNAs), which impact gene expression following pathogen infections, have garnered significant attention in recent years. Recent discoveries have revealed that lncRNAs play a crucial role in fish immune responses to pathogen infections. We investigated the influence of lncRNA-adm2 on the antibacterial immune response generated by Aeromonas hydrophila in grass carp (Ctenopharyngodon idella) through the adsorption of cid-miR-n3. Furthermore, we found that cid-miR-n3 interacts with lncRNA-adm2 and targets the 3' UTR of adm2. The upregulation of lncRNA-adm2 expression led to the suppression of pro-inflammatory cytokines (il-1β and il-6) in CIK cells, while anti-inflammatory cytokines (il-10) increased. Our research provides evidence that lncRNAs are involved in the antibacterial immune response of fish, expanding our understanding of the function of lncRNAs in teleosts.

MicroRNA qPCR normalization in Nile tilapia (Oreochromis niloticus): Effects of acute cold stress on potential reference targets

The Nile tilapia (Oreochromis niloticus) is one of the most important cultured fish worldwide, but tilapia culture is largely affected by low temperatures. Recent studies suggest that microRNAs (miRNAs) regulate cold tolerance traits in fish. In general, qPCR-based methods are the simplest and most accurate forms of miRNA quantification. However, qPCR data heavily depends on appropriate normalization. Therefore, the aim of the present study is to determine whether the expression of previously tested, stably expressed miRNAs are affected by acute cold stress in Nile tilapia. For this purpose, one small nuclear RNA (U6) and six candidate reference miRNAs (miR-23a, miR-25-3, Let-7a, miR-103, miR-99-5, and miR-455) were evaluated in four tissues (blood, brain, liver, and gills) under two experimental conditions (acute cold stress and control) in O. niloticus. The stability of the expression of each candidate reference miRNA was analyzed by four independent methods (the delta Ct method, geNorm, NormFinder, and BestKeeper). Further, consensual comprehensive ranking of stability was built with RefFinder. Overall, miR-103 was the most stable reference miRNA in this study, and miR-103 and Let-7a were the best combination of reference targets. Equally important, Let-7a, miR-23a, and miR-25-3 remained consistently stable across different tissues and experimental groups. Considering all variables, U6, miR-99-5, and miR-455 were the least stable candidates under acute cold stress. Most important, suitable reference miRNAs were validated in O. niloticus, facilitating further accurate miRNA quantification in this species.

Multiomics analysis revealed miRNAs as potential regulators of the immune response in Carassius auratus gills to Aeromonas hydrophila infection

The gill of fish is an important immune organ for pathogen defense, but its microRNA (miRNA) expression and regulatory mechanism remain unclear. In this study, we report on the histopathological and immunohistochemical features of the gills of the crucian carp Carassius auratus challenged with Aeromonas hydrophila. Small RNA libraries of the gills were constructed and sequenced on the Illumina HiSeq 2000 platform. A total of 1,165 differentially expressed miRNAs (DEMs) were identified in gills, of which 539 known and 7 unknown DEMs were significantly screened (p < 0.05). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the potential target genes/proteins were primarily involved in 33 immune-related pathways, in which the inflammatory responses were focused on the Toll-like receptor (TLR), mitogen-activated protein kinase (MAPK), and nuclear factor kappa B (NF-κB) signaling pathways. Moreover, the expression levels of 14 key miRNAs (e.g., miR-10, miR-17, miR-26a, miR-144, miR-145, and miR-146a) and their target genes (e.g., TNFα, TLR4, NF-κB, TAB1, PI3K, and IRAK1) were verified. In addition, the protein levels based on isobaric tags for relative and absolute quantification (iTRAQ) were significantly associated with the results of the quantitative real-time PCR (qRT-PCR) analysis (p < 0.01). miR-17/pre-miR-17 were identified in the regulation expression of the NF-κB target gene, and the phylogenetic tree analysis showed that the pre-miR-17 of C. auratus with the closest similarity to the zebrafish Danio rerio is highly conserved in teleosts. This is the first report of the multi-omics analysis of the miRNAs and proteins in the gills of C. auratus infected with A. hydrophila, thus enriching knowledge on the regulation mechanism of the local immune response in Cyprinidae fish.

Expression Analysis in Atlantic Salmon Liver Reveals miRNAs Associated with Smoltification and Seawater Adaptation

Simple Summary

Smoltification is a developmental process that preadapts Atlantic salmon for a life in seawater. Suboptimal smoltification and poor timing of transfer to seawater is associated with increased mortality. MicroRNAs (miRNAs) are small non-coding genes. They regulate gene expression post-transcriptionally as part of the miRNA induce silencing complex (miRISC) where they guide miRISC to particular mRNAs (target genes). The aim of this study was to identify Atlantic salmon miRNAs expressed in liver that are associated with smoltification and adaptation to seawater as well as to predict their target genes. In total, 62 guide miRNAs were identified, and by their expression patterns they were clustered into three groups. Target gene predictions followed by gene enrichment analysis of the predicted targets indicated that the guide miRNAs were involved in post-transcriptional regulation of important smoltification associated biological processes. Some of these were energy metabolism, protein metabolism and transport, circadian rhythm, stress and immune response. Together, the results indicate that certain miRNAs are involved in the regulation of many of the important changes occurring in the liver during this developmental transition.

Abstract

Optimal smoltification is crucial for normal development, growth, and health of farmed Atlantic salmon in seawater. Here, we characterize miRNA expression in liver to reveal whether miRNAs regulate gene expression during this developmental transition. Expression changes of miRNAs and mRNAs was studied by small-RNA sequencing and microarray analysis, respectively. This revealed 62 differentially expressed guide miRNAs (gDE-miRNAs) that could be divided into three groups with characteristic dynamic expression patterns. Three of miRNA families are known as highly expressed in liver. A rare arm shift was observed during smoltification in the Atlantic salmon-specific novel-ssa-miR-16. The gDE-miRNAs were predicted to target 2804 of the genes revealing expression changes in the microarray analysis. Enrichment analysis revealed that targets were significantly enriched in smoltification-associated biological process groups. These included lipid and cholesterol synthesis, carbohydrate metabolism, protein metabolism and protein transport, immune system genes, circadian rhythm and stress response. The results indicate that gDE-miRNAs may regulate many of the changes associated with this developmental transition in liver. The results pave the way for validation of the predicted target genes and further study of gDE-miRNA and their targets by functional assays.

Selection of References for microRNA Quantification in Japanese Flounder (Paralichthys olivaceus) Normal Tissues and Edwardsiella tarda-Infected Livers

MicroRNA (miRNA) plays essential roles in post-transcriptional regulation of protein coding genes, and the quantitative real-time polymerase chain reaction (qRT-PCR) is the powerful and broadly employed tool to conduct studies of miRNA expression. Identifying appropriate references to normalize quantitative data is a prerequisite to ensure the qRT-PCR accuracy. Until now, there has been no report about miRNA reference for qRT-PCR in Japanese flounder (Paralichthys olivaceus), one important marine cultured fish along the coast of Northern Asia. In this study, combined with miRNA-Seq analysis and literature search, 10 candidates (miR-34a-5p, miR-205-5p, miR-101a-3p, miR-22-3p, miR-23a-3p, miR-210-5p, miR-30c-5p, U6, 5S rRNA, and 18S rRNA) were chosen as potential references to test their expression stability among P. olivaceus tissues, and in livers of P. olivaceus infected with Edwardsiella tarda at different time points. The expression stability of these candidates was analyzed by qRT-PCR and evaluated with Delta CT, BestKeeper, geNorm, as well as NormFinder methods, and RefFinder was employed to estimate the comprehensive ranking according to the four methods. As the result, miR-22-3p and miR-23a-3p were proved to be the suitable combination as reference miRNAs for both P. olivaceus normal tissues and livers infected with E. tarda, and they were successfully applied to normalize miR-7a and miR-221-5p expression in P. olivaceus livers in response to E. tarda infection. All these results provide valuable information for P. olivaceus miRNA quantitative expression analysis in the future.

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.

Evolution after Whole-Genome Duplication: Teleost MicroRNAs

MicroRNAs (miRNAs) are important gene expression regulators implicated in many biological processes, but we lack a global understanding of how miRNA genes evolve and contribute to developmental canalization and phenotypic diversification. Whole-genome duplication events likely provide a substrate for species divergence and phenotypic change by increasing gene numbers and relaxing evolutionary pressures. To understand the consequences of genome duplication on miRNA evolution, we studied miRNA genes following the teleost genome duplication (TGD). Analysis of miRNA genes in four teleosts and in spotted gar, whose lineage diverged before the TGD, revealed that miRNA genes were retained in ohnologous pairs more frequently than protein-coding genes, and that gene losses occurred rapidly after the TGD. Genomic context influenced retention rates, with clustered miRNA genes retained more often than nonclustered miRNA genes and intergenic miRNA genes retained more frequently than intragenic miRNA genes, which often shared the evolutionary fate of their protein-coding host. Expression analyses revealed both conserved and divergent expression patterns across species in line with miRNA functions in phenotypic canalization and diversification, respectively. Finally, major strands of miRNA genes experienced stronger purifying selection, especially in their seeds and 3'-complementary regions, compared with minor strands, which nonetheless also displayed evolutionary features compatible with constrained function. This study provides the first genome-wide, multispecies analysis of the mechanisms influencing metazoan miRNA evolution after whole-genome duplication.

Modulation of hepatic miRNA expression in Atlantic salmon (Salmo salar) by family background and dietary fatty acid composition

This study finds significant differences in hepatic fatty acid composition between four groups of Atlantic salmon (Salmo salar) consisting of offspring from families selected for high and low capacities to express the delta 6 desaturase isomer b and fed diets with 10% or 75% fish oil. The results demonstrated that hepatic lipid metabolism was affected by experimental conditions (diet/family). The fatty acid composition in the four groups mirrored the differences in dietary composition, but it was also associated with the family groups. Small RNA sequencing followed by RT-qPCR identified 12 differentially expressed microRNAs (DE miRNAs), with expression associated with family groups (miR-146 family members, miR-200b, miR-214, miR-221, miR-125, miR-135, miR-137, miR_nov_1), diets (miR-203, miR-462) or both conditions. All the conserved DE miRNAs have been reported as associated with lipid metabolism in other vertebrates. In silico predictions revealed 37 lipid metabolism pathway genes, including desaturases, transcription factors and key enzymes in the synthesis pathways as putative targets (e.g., srebp-1 and 2, Δ6fad_b and c, hmdh, elovl4 and 5b, cdc42). RT-qPCR analysis of selected target genes showed expression changes that were associated with diet and with family groups (d5fad, d6fad_a, srebp-1). There was a reciprocal difference in the abundance of ssa-miR-203a-3p and srebp-1 in one group comparison, whereas other predicted targets did not reveal any evidence of being negatively regulated by degradation. More experimental studies are needed to validate and fully understand the predicted interactions and how the DE miRNAs may participate in the regulation of hepatic lipid metabolism.

MicroRNA-29b modulates the innate immune response by suppressing IFNγs production in orange-spotted grouper (Epinephelus coioides)

Interferon-γ (IFNγ), a type II interferon, is essential to host resistance against various infections. Unlike other vertebrates, fish have two types of IFNγs, IFNγ1 (also named IFNγ-rel) and IFNγ2. MicroRNAs (miRNAs) regulate multiple biological processes by suppressing mRNA translation or inducing mRNA degradation. Among them, miR-29 can directly target IFNγ and affact innate and adaptive immune responses in mice. There are five members of the miR-29 family in orange-spotted grouper (Epinephelus coioides), which share the same miRNA seed region. However, whether miR-29 directly targets E. coioides IFNγs and regulate IFNγ production is still unknown. In the present study, the negative correlation between miR-29b and both IFNγs in immune tissues of healthy E. coioides and grouper spleen cells (GS cells) stimulated with LPS or poly I:C was demonstrated. Moreover, dual-luciferase reporter assays and western blotting were performed to demonstrate that miR-29b suppressed E. coioides IFNγ production. Studies of NO production in GS cells after miR-29b transfection revealed that miR-29b overexpression affected NO production through the downregulation of IFNγ expression. Taken together, our results suggest that miR-29b may directly target E. coioides IFNγs and modulate IFNγ-mediated innate immune responses by suppressing E. coioides IFNγs production.

Selection of reference genes for miRNA quantitative PCR and its application in miR-34a/Sirtuin-1 mediated energy metabolism in Megalobrama amblycephala

MiRNAs are small, non-coding RNAs that downregulate gene expression at post-transcriptional levels. They have emerged as important regulators involved in metabolism, immunity, and cancer. Real-time quantitative PCR (RT-qPCR) is an effective and main method for quantifying target miRNA. For robust RT-qPCR method, suitable reference genes play crucial roles in data normalization. Blunt snout bream (Megalobrama amblycephala) is an economically important aquaculture species; however, no reference genes dedicated for qPCR method has been identified for this species so far. The objective of this study was to screen stable reference genes for miRNA RT-qPCR and demonstrated its application in energy metabolism in blunt snout bream. The stabilities of ten potential reference genes (miR-21-1-5p, miR-107a-3p, miR-222a-3p, miR-146a-5p, miR-101a-3p, miR-22a-3p, miR-103-3p, miR-456-3p, miR-221-3p, and U6 (RNU6A)) were evaluated in nine tissues (brain, muscle, liver, skin, spleen, heart, gill, intestine, and eye) under normal condition and in three tissues (liver, intestine, and spleen) under four stresses (heat stress, ammonia stress, bacterial challenge, and glycolipid stress). Using GeNorm, NormFinder, and RefFinder softwares, we discovered that different tissues and stresses are both important variability factors for the expression stability of miRNAs. After verifying miR-34a/Sirtuin-1 expressions in high-carbohydrate diet-induced blunt snout bream, we eventually identified that the most stable reference gene in this species was miR-221-3p, and the best combination of reference genes were miR-221-3p and miR-103-3p.

miRNA analysis with Prost! reveals evolutionary conservation of organ-enriched expression and post-transcriptional modifications in three-spined stickleback and zebrafish

MicroRNAs (miRNAs) can have organ-specific expression and functions; they can originate from dedicated miRNA genes, from non-canonical miRNA genes, or from mirror-miRNA genes and can also experience post-transcriptional variation. It remains unclear, however, which mechanisms of miRNA production or modification are organ-specific and the extent of their evolutionary conservation. To address these issues, we developed the software Prost! (PRocessing Of Short Transcripts), which, among other features, helps quantify mature miRNAs, accounts for post-transcriptional processing, such as nucleotide editing, and identifies mirror-miRNAs. Here, we applied Prost! to annotate and analyze miRNAs in three-spined stickleback (Gasterosteus aculeatus), a model fish for evolutionary biology reported to have a miRNome larger than most teleost fish. Zebrafish (Danio rerio), a distantly related teleost with a well-known miRNome, served as comparator. Our results provided evidence for the existence of 286 miRNA genes and 382 unique mature miRNAs (excluding mir430 gene duplicates and the vaultRNA-derived mir733), which doesn't represent a miRNAome larger than other teleost miRNomes. In addition, small RNA sequencing data from brain, heart, testis, and ovary in both stickleback and zebrafish identified suites of mature miRNAs that display organ-specific enrichment, many of which are evolutionarily-conserved in the brain and heart in both species. These data also supported the hypothesis that evolutionarily-conserved, organ-specific mechanisms may regulate post-transcriptional variations in miRNA sequence. In both stickleback and zebrafish, miR2188-5p was edited frequently with similar nucleotide changes in the seed sequence with organ specific editing rates, highest in the brain. In summary, Prost! is a new tool to identify and understand small RNAs, to help clarify a species' miRNA biology as shown here for an important model for the evolution of developmental mechanisms, and to provide insight into organ-enriched expression and the evolutionary conservation of miRNA post-transcriptional modifications.

Expanding the miRNA Repertoire in Atlantic Salmon; Discovery of IsomiRs and miRNAs Highly Expressed in Different Tissues and Developmental Stages

MicroRNAs (miRNAs) are important post-transcriptional gene expression regulators. Here, 448 different miRNA genes, including 17 novel miRNAs, encoding for 589 mature Atlantic salmon miRNAs were identified after sequencing 111 samples (fry, pathogen challenged fry, various developmental and adult tissues). This increased the reference miRNAome with almost one hundred genes. Prior to isomiR characterization (mature miRNA variants), the proportion of erroneous sequence variants (ESVs) arising in the analysis pipeline was assessed. The ESVs were biased towards 5’ and 3’ end of reads in unexpectedly high proportions indicating that measurements of ESVs rather than Phred score should be used to avoid misinterpreting ESVs as isomiRs. Forty-three isomiRs were subsequently discovered. The biological effect of the isomiRs measured as increases in target diversity was small (<3%). Five miRNA genes showed allelic variation that had a large impact on target gene diversity if present in the seed. Twenty-one miRNAs were ubiquitously expressed while 31 miRNAs showed predominant expression in one or few tissues, indicating housekeeping or tissue specific functions, respectively. The miR-10 family, known to target Hox genes, were highly expressed in the developmental stages. The proportion of miR-430 family members, participating in maternal RNA clearance, was high at the earliest developmental stage.

MicroRNAs profiles of Chinese Perch Brain (CPB) cells infected with Siniperca chuatsi rhabdovirus (SCRV)

MicroRNAs are non-coding RNAs, which widely participate in biological processes. In recent years, Siniperca chuatsi rhabdovirus (SCRV) has caused mass mortality in Chinese perch (Siniperca chuatsi). To identify specific miRNAs involved in SCRV infection, deep sequencing of microRNA on Chinese perch brain cell line (CPB) with or without SCRV infection were performed at 6 and 12 h post of infection (hpi). Totally 382 miRNAs were identified, including 217 known miRNA aligned with zebrafish miRNAs and 165 novel miRNAs by MiRDeep2 program. Of which 15 and 35 differentially-expressed miRNAs were determined respectively to 6 and 12 hpi. Nine miRNAs were selected randomly from the differentially-expressed miRNAs and validated by quantitative real-time PCR (qRT-PCR). These results were consistent with the microRNA sequencing results. Besides, target genes of 98 differentially-expressed miRNAs were predicted. Three of miRNAs (miR-122, miR-214, miR-135a) were selected, and its effects were analyzed in CPC cells transfected with appropriate miRNA mimics/inhibitors to evaluate its regulation effects by qRT-PCR and western blot. The results demonstrated that miR-214 inhibited the replication of SCRV, while miR-122 promoted the replication of SCRV and there was no correlation between the miR-135a and SCRV replication. These results will pave a new way for the development of effective strategies against the SCRV infection.

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

Non-coding RNAs Function as Immune Regulators in Teleost Fish

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

Identification 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.

Discovery of microRNAs associated with the antiviral immune response of Atlantic cod macrophages

MicroRNAs (miRNAs) are known to play important immunoregulatory roles in teleosts, although miRNAs involved in the antiviral immune response of Atlantic cod (Gadus morhua) were previously uncharacterised. Using deep sequencing and qPCR, the present study was conducted to identify miRNAs responsive to the viral mimic, polyriboinosinic polyribocytidylic acid (pIC) in Atlantic cod macrophages. Macrophage samples isolated from Atlantic cod (n=3) and treated with pIC or phosphate buffered saline (PBS control) for 24 and 72h were used for miRNA profiling. Following deep sequencing, DESeq2 analyses identified four (miR-731-3p, miR-125b-3-3p, miR-150-3p and miR-462-3p) and two (miR-2188-3p and miR-462-3p) significantly differentially expressed miRNAs at 24 and 72h post-stimulation (HPS), respectively. Sequencing-identified miRNAs were subjected to qPCR validation using a larger number of biological replicates (n=6) exposed to pIC or PBS over time (i.e. 12, 24, 48 and 72 HPS). As in sequencing, miR-731-3p, miR-462-3p and miR-2188-3p showed significant up-regulation by pIC. The sequencing results were not qPCR-validated for miR-125b-3-3p and miR-150-3p as up- and down-regulated miRNAs at 24 HPS, respectively; however, qPCR results showed significant up-regulation in response to pIC stimulation at later time points (i.e. 48 and/or 72 HPS). We also used qPCR to assess the expression of other miRNAs that were previously shown as immune responsive in other vertebrates. qPCR results at 48 and/or 72 HPS revealed that miR-128-3-5p, miR-214-1-5p and miR-451-3p were induced by pIC, whereas miR-30b-3p and miR-199-1-3p expression were repressed in response to pIC. The present study identified ten pIC-stimulated miRNAs, suggesting them as important in antiviral immune responses of Atlantic cod macrophages. Some pIC-responsive miRNAs identified in this study were predicted to target putative immune-related genes of Atlantic cod (e.g. miR-30b-3p targeting herc4), although the regulatory functions of these miRNAs need to be validated by future studies.

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.