Discovery and characterization of miRNA genes in Atlantic salmon (Salmo salar) by use of a deep sequencing approach

Differential Digestive Stability of Food-Derived microRNAs: The Case of miR-30c-5p and miR-92a-3p in Polyfloral Honey

Dietary microRNAs (miRs) represent a new area in food science. Although they have been found in many foods, including honey, more research is needed about their stability and fate during digestion. Hence, this study aimed to analyze the digestive stability of two selected miRs in honey. We extracted miR-92a-3p and miR-30c-5p from pasteurized and unpasteurized forms of polyfloral honey using two different methods and kits: a column-based manual method and a phenol-free semi-automated magnetic-bead-based method. The latter option was used for the subsequent analysis of samples according to the INFOGEST static in vitro digestion protocol. Also, the honey samples were examined for exosome-like particles using dynamic light scattering. Although the expression levels of both miRs were significantly lower following intestinal digestion, we found a difference in the resilience of the miRs to gastrointestinal conditions, with miR-30c-5p being relatively stable compared to miR-92a-3p following digestion, regardless of the honey's pasteurization treatment. Moreover, there was marked heterogeneity in the extracellular vesicle profile of the pasteurized sample. We identified the presence of two broadly conserved miRs in honey: miR-92a-3p and miR-30c-5p. Despite honey exhibiting high digestibility, miR-92a-3p was less resilient than miR-30c-5p, demonstrating considerable resistance under gastrointestinal conditions. Although further research is needed, the results obtained from this study may represent a starting point for utilizing honey as a source of exogenous miRNAs for preventive strategies and more "natural" treatments.

Expression Analysis of Moritella viscosa-Challenged Atlantic Salmon Identifies Disease-Responding Genes, MicroRNAs and Their Predicted Target Genes and Pathways

Moritella viscosa is a bacterial pathogen causing winter-ulcer disease in Atlantic salmon. The lesions on affected fish lead to increased mortality, decreased fish welfare, and inferior meat quality in farmed salmon. MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional regulation by guiding the miRNA-induced silencing complex to specific mRNA transcripts (target genes). The goal of this study was to identify miRNAs responding to Moritella viscosa in salmon by investigating miRNA expression in the head-kidney and the muscle/skin from lesion sites caused by the pathogen. Protein coding gene expression was investigated by microarray analysis in the same materials. Seventeen differentially expressed guide-miRNAs (gDE-miRNAs) were identified in the head-kidney, and thirty-nine in lesion sites, while the microarray analysis reproduced the differential expression signature of several thousand genes known as infection-responsive. In silico target prediction and enrichment analysis suggested that the gDE-miRNAs were predicted to target genes involved in immune responses, hemostasis, angiogenesis, stress responses, metabolism, cell growth, and apoptosis. The majority of the conserved gDE-miRNAs (e.g., miR-125, miR-132, miR-146, miR-152, miR-155, miR-223 and miR-2188) are known as infection-responsive in other vertebrates. Collectively, the findings indicate that gDE-miRNAs are important post-transcriptional gene regulators of the host response to bacterial infection.

Identification and expression analysis of sex biased miRNAs in chinese hook snout carp Opsariichthys bidens

As an economically important fish, Opsariichthys bidens has obvious sexual dimorphism and strong reproductive capacity, but no epigenetics study can well explain its phenotypic variations. In recent years, many microRNAs involved in the regulation of reproductive development have been explored. In this study, the small RNA libraries of O. bidens on the testis and ovary were constructed and sequenced. A total of 295 known miRNAs were obtained and 100 novel miRNAs were predicted. By comparing testis and ovary libraries, 115 differentially expressed (DE) miRNAs were selected, of which 53 were up-regulated and 62 were down-regulated. A total of 64 GO items (padj < 0.01) and 206 KEGG pathways (padj < 0.01) were enriched in the target gene of miRNA. After that, the expression levels of nine DE miRNAs, including let-7a, miR-146b, miR-18c, miR-202-5p, miR-135c, miR-9-5p, miR-34c-3p, miR-460-5p and miR-338 were verified by qRT-PCR. Furthermore, bidirectional prediction of DE miRNAs and sex-related genes was carried out and the targeting correlation between miR-9-5p and nanos1 was verified by Dual-Luciferase reporter assay. Our findings identified the differentially expressed miRNA and paved the way to new possibilities for the follow-up study on the mechanism of miRNA-mRNA interaction in the gonads of O. bidens.

Differential Expression of miRNAs and Their Predicted Target Genes Indicates That Gene Expression in Atlantic Salmon Gill Is Post-Transcriptionally Regulated by miRNAs in the Parr-Smolt Transformation and Adaptation to Sea Water

Smoltification (parr-smolt transformation) is a complex developmental process consisting of developmental changes that lead to remodeling of the Atlantic salmon gill. Here, the expression changes of miRNAs and mRNAs were studied by small-RNA sequencing and microarray analysis, respectively, to identify miRNAs and their predicted targets associated with smoltification and subsequent sea water adaptation (SWA). In total, 18 guide miRNAs were identified as differentially expressed (gDE miRNAs). Hierarchical clustering analysis of expression changes divided these into one cluster of 13 gDE miRNAs with decreasing expression during smoltification and SWA that included the miRNA-146, miRNA-30 and miRNA-7132 families. Another smaller cluster that showed increasing expression consisted of miR-101a-3p, miR-193b-5p, miR-499a-5p, miR-727a-3p and miR-8159-5p. The gDE miRNAs were predicted to target 747 of the genes (DE mRNAs), showing expression changes in the microarray analysis. The predicted targets included genes encoding NKA-subunits, aquaporin-subunits, cystic fibrosis transmembrane conductance regulator and the solute carrier family. Furthermore, the predicted target genes were enriched in biological processes associated with smoltification and SWA (e.g., immune system, reactive oxygen species, stress response and extracellular matrix organization). Collectively, the results indicate that remodeling of the gill involves the post-transcriptional regulation of gene expression by the characterized gDE miRNAs.

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.

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.

Whole-Genome Transcript Expression Profiling Reveals Novel Insights into Transposon Genes and Non-Coding RNAs during Atlantic Salmon Seawater Adaptation

The growing amount of genome information and transcriptomes data available allows for a better understanding of biological processes. However, analysis of complex transcriptomic experimental designs involving different conditions, tissues, or times is relevant. This study proposes a novel approach to analyze complex data sets combining transcriptomes and miRNAs at the chromosome-level genome. Atlantic salmon smolts were transferred to seawater under two strategies: (i) fish group exposed to gradual salinity changes (GSC) and (ii) fish group exposed to a salinity shock (SS). Gills, intestine, and head kidney samples were used for total RNA extraction, followed by mRNA and small RNA illumina sequencing. Different expression patterns among the tissues and treatments were observed through a whole-genome transcriptomic approach. Chromosome regions highly expressed between experimental conditions included a great abundance of transposable elements. In addition, differential expression analysis showed a greater number of transcripts modulated in response to SS in gills and head kidney. miRNA expression analysis suggested a small number of miRNAs involved in the smoltification process. However, target analysis of these miRNAs showed a regulatory role in growth, stress response, and immunity. This study is the first to evidence the interplaying among mRNAs and miRNAs and the structural relationship at the genome level during Atlantic salmon smoltification.

MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3'UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome

Complete 3′UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3′UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3′UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3′UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3′UTR motifs.

Aberration-corrected ultrafine analysis of miRNA reads at single-base resolution: a k-mer lattice approach

Raw sequencing reads of miRNAs contain machine-made substitution errors, or even insertions and deletions (indels). Although the error rate can be low at 0.1%, precise rectification of these errors is critically important because isoform variation analysis at single-base resolution such as novel isomiR discovery, editing events understanding, differential expression analysis, or tissue-specific isoform identification is very sensitive to base positions and copy counts of the reads. Existing error correction methods do not work for miRNA sequencing data attributed to miRNAs’ length and per-read-coverage properties distinct from DNA or mRNA sequencing reads. We present a novel lattice structure combining kmers, (k – 1)mers and (k + 1)mers to address this problem. The method is particularly effective for the correction of indel errors. Extensive tests on datasets having known ground truth of errors demonstrate that the method is able to remove almost all of the errors, without introducing any new error, to improve the data quality from every-50-reads containing one error to every-1300-reads containing one error. Studies on experimental miRNA sequencing datasets show that the errors are often rectified at the 5′ ends and the seed regions of the reads, and that there are remarkable changes after the correction in miRNA isoform abundance, volume of singleton reads, overall entropy, isomiR families, tissue-specific miRNAs, and rare-miRNA quantities.

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

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

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.

Characterization of miRNAs in Extracellular Vesicles Released From Atlantic Salmon Monocyte-Like and Macrophage-Like Cells

Cell-derived extracellular vesicles (EVs) participate in cell-cell communication via transfer of molecular cargo including genetic material like miRNAs. In mammals, it has previously been established that EV-mediated transfer of miRNAs can alter the development or function of immune cells, such as macrophages. Our previous research revealed that Atlantic salmon head kidney leukocytes (HKLs) change their morphology, phagocytic ability and miRNA profile from primarily “monocyte-like” at Day 1 to primarily “macrophage-like” at Day 5 of culture. Therefore, we aimed to characterize the miRNA cargo packaged in EVs released from these two cell populations. We successfully isolated EVs from Atlantic salmon HKL culture supernatants using the established Vn96 peptide-based pull-down. Isolation was validated using transmission electron microscopy, nanoparticle tracking analysis, and Western blotting. RNA-sequencing identified 19 differentially enriched (DE) miRNAs packaged in Day 1 versus Day 5 EVs. Several of the highly abundant miRNAs, including those that were DE (e.g. ssa-miR-146a, ssa-miR-155 and ssa-miR-731), were previously identified as DE in HKLs and are associated with macrophage differentiation and immune response in other species. Interestingly, the abundance relative of the miRNAs in EVs, including the most abundant miRNA (ssa-miR-125b), was different than the miRNA abundance in HKLs, indicating selective packaging of miRNAs in EVs. Further study of the miRNA cargo in EVs derived from fish immune cells will be an important next step in identifying EV biomarkers useful for evaluating immune cell function, fish health, or response to disease.

Characterization of Differentially Expressed miRNAs and Their Predicted Target Transcripts during Smoltification and Adaptation to Seawater in Head Kidney of Atlantic Salmon

Smoltification and early seawater phase are critical developmental periods with physiological and biochemical changes in Atlantic salmon that facilitates survival in saltwater. MicroRNAs (miRNAs) are known to have important roles in development, but whether any miRNAs are involved in regulation of gene expression during smoltification and the adaption to seawater is largely unknown. Here, small RNA sequencing of materials from head kidney before, during smoltification and post seawater transfer were used to study expression dynamics of miRNAs, while microarray analysis was applied to study mRNA expression dynamics. Comparing all timepoints, 71 miRNAs and 2709 mRNAs were identified as differentially expressed (DE). Hierarchical clustering analysis of the DE miRNAs showed three major clusters with characteristic expression changes. Eighty-one DE mRNAs revealed negatively correlated expression patterns to DE miRNAs in Cluster I and III. Furthermore, 42 of these mRNAs were predicted as DE miRNA targets. Gene enrichment analysis of negatively correlated target genes showed they were enriched in gene ontology groups hormone biosynthesis, stress management, immune response, and ion transport. The results strongly indicate that post-transcriptional regulation of gene expression by miRNAs is important in smoltification and sea water adaption, and this study identifies several putative miRNA-target pairs for further functional studies.

Heterogenic Origin of Micro RNAs in Atlantic Salmon (Salmo salar) Seminal Plasma

The origin and contribution of seminal plasma RNAs into the whole semen RNA repertoire are poorly known, frequently being overlooked or neglected. In this study, we used high-throughput sequencing and RT-qPCR to profile microRNA (miRNA) constituents in the whole semen, as well as in fractionated spermatozoa and seminal plasma of Atlantic salmon (Salmo salar). We found 85 differentially accumulated miRNAs between spermatozoa and the seminal plasma. We identified a number of seminal plasma-enriched and spermatozoa-enriched miRNAs. We localized the expression of some miRNAs in juvenile and mature testes. Two abundant miRNAs, miR-92a-3p and miR-202-5p, localized to both spermatogonia and somatic supporting cells in immature testis, and they were also highly abundant in somatic cells in mature testis. miR-15c-5p, miR-30d-5p, miR-93a-5p, and miR-730-5p were detected only in mature testis. miRs 92a-3p, 202-5p, 15c-5p, and 30d-5p were also detected in a juvenile ovary. The RT-qPCR experiment demonstrated lack of correlation in miRNA transcript levels in seminal plasma versus blood plasma. Our results indicate that salmon semen is rich in miRNAs, which are present in both spermatozoa and seminal plasma. Testicular-supporting somatic cells are likely the source of seminal plasma enrichment, whereas blood plasma is unlikely to contribute to the seminal plasma miRNA repertoire.

SalMotifDB: a tool for analyzing putative transcription factor binding sites in salmonid genomes

Background

Recently developed genome resources in Salmonid fish provides tools for studying the genomics underlying a wide range of properties including life history trait variation in the wild, economically important traits in aquaculture and the evolutionary consequences of whole genome duplications. Although genome assemblies now exist for a number of salmonid species, the lack of regulatory annotations are holding back our mechanistic understanding of how genetic variation in non-coding regulatory regions affect gene expression and the downstream phenotypic effects.

Results

We present SalMotifDB, a database and associated web and R interface for the analysis of transcription factors (TFs) and their cis-regulatory binding sites in five salmonid genomes. SalMotifDB integrates TF-binding site information for 3072 non-redundant DNA patterns (motifs) assembled from a large number of metazoan motif databases. Through motif matching and TF prediction, we have used these multi-species databases to construct putative regulatory networks in salmonid species. The utility of SalMotifDB is demonstrated by showing that key lipid metabolism regulators are predicted to regulate a set of genes affected by different lipid and fatty acid content in the feed, and by showing that our motif database explains a significant proportion of gene expression divergence in gene duplicates originating from the salmonid specific whole genome duplication.

Conclusions

SalMotifDB is an effective tool for analyzing transcription factors, their binding sites and the resulting gene regulatory networks in salmonid species, and will be an important tool for gaining a better mechanistic understanding of gene regulation and the associated phenotypes in salmonids. SalMotifDB is available at https://salmobase.org/apps/SalMotifDB.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6051-0) contains supplementary material, which is available to authorized users.

Identification and expression of microRNAs in european eels Anguilla anguilla from two natural sites with different pollution levels

MicroRNAs (miRNAs) are a class of small non-coding RNA that control multiple biological processes through negative post-transcriptional regulation of gene expression. Recently a role of miRNAs in the response of aquatic organisms to environmental toxicants emerged. Toxicant-induced changes in miRNA expression might then represent novel biomarkers to evaluate the health status of these organisms. In this study, we aimed to identify the miRNA repertoire in the liver of the European eel Anguilla anguilla and to compare their differential expression between a polluted site located in the Gironde Estuary and a pristine site in Arcachon Bay (France). A total of 299 mature miRNAs were identified. In polluted water, 19 miRNAs were up-regulated and 22 were down-regulated. We predicted that these differentially expressed miRNAs could target 490 genes that were involved in ribosome biogenesis, response to hormones, response to chemical and chromatin modification. Moreover, we observed only few examples (29) of negative correlation between the expression levels of miRNAs and their targets suggesting that, in the system studied, miRNAs might not only regulate gene expression directly by degrading mRNA but also by inhibiting protein translation or by regulating other epigenetic processes. This study is the first example of in situ investigation of the role of miRNAs in the response of a fish species to water quality. Our findings provide new insights into the involvement of epigenetic mechanisms in the response of animals chronically exposed to pollution and pave the way for the utilization of miRNAs in aquatic ecotoxicology.

Discovery of MicroRNAs from Batrachuperus yenyuanensis Using Deep Sequencing and Prediction of Their Targets

MicroRNAs (miRNAs), a family of ∼22-nucleotide non-coding single-stranded RNA molecules, are considered as key post-transcriptional regulators of gene expression that regulate various biological processes in living organism. Many miRNAs have been identified in animals; however, few have been reported in Hynobiidae species. The present study is aimed to identify a full repertoire of miRNAs in Batrachuperus yenyuanensis (Yenyuan stream salamander), which would significantly increase our knowledge of miRNAs in amphibians. A small RNA library was constructed from B. yenyuanensis and sequenced using deep sequencing. As a result, 1,717,751 clean reads were obtained, representing 356 known and 80 novel miRNAs. Additionally, expression levels of eight randomly selected miRNAs in B. yenyuanensis were confirmed using the stem-loop quantitative real-time reverse transcription PCR. In addition, 13,972 targets were predicted for these identified miRNAs, although the physiological functions of many of these targets remain unknown. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested that the predicted targets are involved in a variety of physiological regulatory functions in B. yenyuanensis. These results provide useful information for further research on the miRNAs involved in the growth and development of B. yenyuanensis, as well as adaptation of this species to its high-altitude habitats.

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

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

High-throughput sequencing reveals microRNAs in response to heat stress in the head kidney of rainbow trout (Oncorhynchus mykiss)

Recently, the research of animal microRNAs (miRNAs) has attracted wide attention for its regulatory effect in the development process and the response to abiotic stresses. Rainbow trout is a commercially and cold water fish species, and usually encounters heat stress, which affects its growth and leads to a huge economic loss. But there were few investigations about the roles of miRNAs in heat stress in rainbow trout. In this study, miRNAs of rainbow trout which were involved in heat stress were identified by high-throughput sequencing of six small RNA libraries from head kidney tissues under control (18 °C) and heat-treated (24 °C) conditions. A total of 392 conserved miRNAs and 989 novel miRNAs were identified, of which 78 miRNAs were expressed in different response to heat stress. Ten of these miRNAs were further validated by quantitative real-time PCR. In addition to, including 393 negative correlation miRNA-target gene pairs, several important regulatory pathways were involved in heat stress of the potential target genes, including protein processing in endoplasmic reticulum, NOD-like receptor signaling pathway, and phagosome. Our data significantly advance understanding of heat stress regulatory mechanism of miRNA in the head kidney of rainbow trout, which provide a useful resource for the cultivation of rainbow trout.

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.

Marine Environmental Epigenetics

Marine organisms' persistence hinges on the capacity for acclimatization and adaptation to the myriad of interacting environmental stressors associated with global climate change. In this context, epigenetics-mechanisms that facilitate phenotypic variation through genotype-environment interactions-are of great interest ecologically and evolutionarily. Our comprehensive review of marine environmental epigenetics guides our recommendations of four key areas for future research: the dynamics of wash-in and wash-out of epigenetic effects, the mechanistic understanding of the interplay of different epigenetic marks and the interaction with the microbiome, the capacity for and mechanisms of transgenerational epigenetic inheritance, and the evolutionary implications of the interaction of genetic and epigenetic features. Emerging insights in marine environmental epigenetics can be applied to critical issues such as aquaculture, biomonitoring, and biological invasions, thereby improving our ability to explain and predict the responses of marine taxa to global climate change.

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.

Deep sequencing of small RNAs from 11 tissues of grass carp Ctenopharyngodon idella and discovery of sex-related microRNAs

This research identified 169 known microRNAs (miRNAs), 380 novel miRNAs, and 30,538 targets in 11 tissues (blood, brain, derma, gill, heart, intestine, kidney, liver, muscle, pronephros, and spleen) from grass carp Ctenopharyngodon idella with high-throughput sequencing (HTS). Transcripts per million (TPM) expression analysis detected 41 brain-enriched miRNAs (accounting for 61.19% of all tissue-enriched miRNAs). Real-time quantitative PCR (RT-qPCR) confirmed that 21 of 24 randomly selected tissue-enriched miRNAs from the TPM analysis were indeed tissue-enriched (P < 0.05), suggesting the HTS and TPM analyses were reliable. Nine of the 41 brain-enriched miRNAs are complementary to members of the double-sex and mab-3 related transcription factor family (dmrt) involved in sex differentiation. RT-qPCR revealed that cid-miR-138 was more highly expressed in testis than in ovary (P < 0.01), while the reverse was true for target gene dmrt4a (P < 0.01). This opposite expression pattern suggested the direct participation of cid-miR-138-dmrt4a in neuroendocrine mechanisms related to brain-pituitary networks during sex development. The discovery of miRNAs from 11 C. idella tissues expands the available fish miRNA database, and enhances our understanding of the role of sex-related miRNAs in tissue differentiation and maintenance of specific tissue functions in fishes.

Involvement of selected cellular miRNAs in the in vitro and in vivo infection of infectious salmon anemia virus (ISAV)

Infectious salmon anemia virus (ISAV) is the causative agent of infectious salmon anemia (ISA), a relatively novel disease primarily affecting farmed salmon species, primarily in Salmo salar specimens, causing severe outbreaks in most producer countries. Although ISAV has been extensively studied at the molecular level, not much is known about the host/cell interaction at the small RNA level. MicroRNAs (miRNAs) are small, non-coding RNA that regulate mRNA expression at the post-transcriptional level. In recent years, the putative role of these molecules in host-pathogen interactions has drawn particular attention because of their pivotal involvement as regulatory elements in a number of eukaryotic organisms. Given the importance of the salmon industry in Chile, a deep understanding of the interaction between ISAV and its hosts is of importance. In the present work, we studied the kinetic expression of selected miRNAs during ISAV infection, both in vitro and in vivo. Based on initial experimental data derived from a small RNA-Seq analysis, a group of miRNAs that were differentially expressed in infected cells were selected for analysis. As a result, two miRNAs, miR-462a-5p and miR-125 b-5p, showed increased and decreased expression, respectively, during ISAV infection.

Genome-wide microRNA screening in Nile tilapia reveals pervasive isomiRs' transcription, sex-biased arm switching and increasing complexity of expression throughout development

MicroRNAs (miRNAs) are key regulators of gene expression in multicellular organisms. The elucidation of miRNA function and evolution depends on the identification and characterization of miRNA repertoire of strategic organisms, as the fast-evolving cichlid fishes. Using RNA-seq and comparative genomics we carried out an in-depth report of miRNAs in Nile tilapia (Oreochromis niloticus), an emergent model organism to investigate evo-devo mechanisms. Five hundred known miRNAs and almost one hundred putative novel vertebrate miRNAs have been identified, many of which seem to be teleost-specific, cichlid-specific or tilapia-specific. Abundant miRNA isoforms (isomiRs) were identified with modifications in both 5p and 3p miRNA transcripts. Changes in arm usage (arm switching) of nine miRNAs were detected in early development, adult stage and even between male and female samples. We found an increasing complexity of miRNA expression during ontogenetic development, revealing a remarkable synchronism between the rate of new miRNAs recruitment and morphological changes. Overall, our results enlarge vertebrate miRNA collection and reveal a notable differential ratio of miRNA arms and isoforms influenced by sex and developmental life stage, providing a better picture of the evolutionary and spatiotemporal dynamics of miRNAs.

MicroRNA regulation of Toll-like receptor signaling pathways in teleost fish

The innate immune system is the first line defense mechanism that recognizes, responds to, controls or eliminates invading pathogens. Toll-like receptors (TLRs) are a critical family of pattern recognition receptors (PRRs) tightly regulated by complex mechanisms involving many molecules to ensure a beneficial outcome in response to foreign invaders. MicroRNAs (miRNAs), a transcriptional and posttranscriptional regulator family in a wide range of biological processes, have been identified as new molecules related to the regulation of TLR-signaling pathways in immune responses. To date, at least 22 TLR types have been identified in more than a dozen different fish species. However, the functions and underlying mechanisms of miRNAs in the regulation of inflammatory responses related to the TLR-signaling pathway in fish is lacking. In this review, we summarize the regulation of miRNA expression profiles in the presence of TLR ligands or pathogen infections in teleost fish. We focus on the effects of miRNAs in regulating TLR-signaling pathways by targeting multiple molecules, including TLRs themselves, TLR-associated signaling proteins, and TLR-induced cytokines. An understanding of the relationship between the TLR-signaling pathways and miRNAs may provide new insights for drug intervention to manipulate immune responses in fish.

Fishing Into the MicroRNA Transcriptome

In the last decade, several studies have been focused on revealing the microRNA (miRNA) repertoire and determining their functions in farm animals such as poultry, pigs, cattle, and fish. These small non-protein coding RNA molecules (18-25 nucleotides) are capable of controlling gene expression by binding to messenger RNA (mRNA) targets, thus interfering in the final protein output. MiRNAs have been recognized as the main regulators of biological features of economic interest, including body growth, muscle development, fat deposition, and immunology, among other highly valuable traits, in aquatic livestock. Currently, the miRNA repertoire of some farmed fish species has been identified and characterized, bringing insights about miRNA functions, and novel perspectives for improving health and productivity. In this review, we summarize the current advances in miRNA research by examining available data on Neotropical and other key species exploited by fisheries and in aquaculture worldwide and discuss how future studies on Neotropical fish could benefit from this knowledge. We also make a horizontal comparison of major results and discuss forefront strategies for miRNA manipulation in aquaculture focusing on forward-looking ideas for forthcoming research.

Identification and differential expression of microRNAs in testis and ovary of Amur sturgeon (Acipenser schrenckii)

BACKGROUND

MicroRNAs (miRNAs) cooperate with sex-related genes in post-transcriptional regulation and play extremely important roles in the establishment of sexually dimorphic traits in animals. However, the gonad miRNAs and expression patterns of miRNAs in sturgeon have not been investigated.

METHODS

In the present study, we used high-throughput small RNA sequencing (RNA-Seq) to discover gonad miRNAs from the ovaries and testes of Amur sturgeons (Acipenser schrenckii). Further, microarray and real-time PCR assays were performed to identify the expression patterns of gonad miRNAs.

RESULTS

As a result, a total of 679 conserved and 51 novel miRNAs were successfully discovered in the gonads of A. schrenckii. Moreover, we found wide sequence variations (isomiRs) in gonad miRNAs, including 5' and 3' isomiRs. Our microarray analysis further characterized the 730 miRNAs expression profiles, which indicated that 117 differentially expressed miRNAs were detected with sex-biased patterns: 71 testis-biased and 46 ovary-biased miRNAs. Based on bioinformatics prediction, we found that there were functional differences between the testis-biased and ovary-biased miRNA targets involved in reproductive-related GO and KEGG pathways. Further, the association of the differentially expressed miRNAs and sex-related target mRNAs was uncovered. Finally, the expression patterns of 11 sex-biased miRNAs and 7 sex-related targets were validated in testes and ovaries using real-time PCR. Putative, negatively expressed miRNA-mRNA relationships were confirmed, such as Dmrt1 and asc-miR-2779, AR and asc-miR-203b-3p, foxl2 and asc-miR-30d.

CONCLUSION

This study provides information regarding the gonad miRNAs in sturgeon. The differential expression miRNAs in the gonads will help us to further understand the role of miRNA-mediated post-transcriptional regulation in the ovary and testis of sturgeon.

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.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.

Epigenetic considerations in aquaculture

Epigenetics has attracted considerable attention with respect to its potential value in many areas of agricultural production, particularly under conditions where the environment can be manipulated or natural variation exists. Here we introduce key concepts and definitions of epigenetic mechanisms, including DNA methylation, histone modifications and non-coding RNA, review the current understanding of epigenetics in both fish and shellfish, and propose key areas of aquaculture where epigenetics could be applied. The first key area is environmental manipulation, where the intention is to induce an ‘epigenetic memory’ either within or between generations to produce a desired phenotype. The second key area is epigenetic selection, which, alone or combined with genetic selection, may increase the reliability of producing animals with desired phenotypes. Based on aspects of life history and husbandry practices in aquaculture species, the application of epigenetic knowledge could significantly affect the productivity and sustainability of aquaculture practices. Conversely, clarifying the role of epigenetic mechanisms in aquaculture species may upend traditional assumptions about selection practices. Ultimately, there are still many unanswered questions regarding how epigenetic mechanisms might be leveraged in aquaculture.

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

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

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.

Modulation of Atlantic salmon miRNome response to sea louse infestation

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

Integrative testis transcriptome analysis reveals differentially expressed miRNAs and their mRNA targets during early puberty in Atlantic salmon

BACKGROUND

Our understanding of the molecular mechanisms implementing pubertal maturation of the testis in vertebrates is incomplete. This topic is relevant in Atlantic salmon aquaculture, since precocious male puberty negatively impacts animal welfare and growth. We hypothesize that certain miRNAs modulate mRNAs relevant for the initiation of puberty. To explore which miRNAs regulate mRNAs during initiation of puberty in salmon, we performed an integrated transcriptome analysis (miRNA and mRNA-seq) of salmon testis at three stages of development: an immature, long-term quiescent stage, a prepubertal stage just before, and a pubertal stage just after the onset of single cell proliferation activity in the testis.

RESULTS

Differentially expressed miRNAs clustered into 5 distinct expression profiles related to the immature, prepubertal and pubertal salmon testis. Potential mRNA targets of these miRNAs were predicted with miRmap and filtered for mRNAs displaying negatively correlated expression patterns. In summary, this analysis revealed miRNAs previously known to be regulated in immature vertebrate testis (miR-101, miR-137, miR-92b, miR-18a, miR-20a), but also miRNAs first reported here as regulated in the testis (miR-new289, miR-30c, miR-724, miR-26b, miR-new271, miR-217, miR-216a, miR-135a, miR-new194 and the novel predicted n268). By KEGG enrichment analysis, progesterone signaling and cell cycle pathway genes were found regulated by these differentially expressed miRNAs. During the transition into puberty we found differential expression of miRNAs previously associated (let7a/b/c), or newly associated (miR-15c, miR-2184, miR-145 and the novel predicted n7a and b) with this stage. KEGG enrichment analysis revealed that mRNAs of the Wnt, Hedgehog and Apelin signaling pathways were potential regulated targets during the transition into puberty. Likewise, several regulated miRNAs in the pubertal stage had earlier been associated (miR-20a, miR-25, miR-181a, miR-202, let7c/d/a, miR-125b, miR-222a/b, miR-190a) or have now been found connected (miR-2188, miR-144, miR-731, miR-8157 and the novel n2) to the initiation of puberty.

CONCLUSIONS

This study has - for the first time - linked testis maturation to specific miRNAs and their inversely correlated expressed targets in Atlantic salmon. The study indicates a broad functional conservation of already known miRNAs and associated pathways involved in the transition into puberty in vertebrates. The analysis also reveals miRNAs not previously associated with testis tissue or its maturation, which calls for further functional studies in the testis.

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.

Discovery and characterization of the feline miRNAome

The domestic cat is an important human companion animal that can also serve as a relevant model for ~250 genetic diseases, many metabolic and degenerative conditions, and forms of cancer that are analogous to human disorders. MicroRNAs (miRNAs) play a crucial role in many biological processes and their dysregulation has a significant impact on important cellular pathways and is linked to a variety of diseases. While many species already have a well-defined and characterized miRNAome, miRNAs have not been carefully studied in cats. As a result, there are no feline miRNAs present in the reference miRNA databases, diminishing the usefulness of medical research on spontaneous disease in cats for applicability to both feline and human disease. This study was undertaken to define and characterize the cat miRNAome in normal feline tissues. High-throughput sequencing was performed on 12 different normal cat tissues. 271 candidate feline miRNA precursors, encoding a total of 475 mature sequences, were identified, including several novel cat-specific miRNAs. Several analyses were performed to characterize the discovered miRNAs, including tissue distribution of the precursors and mature sequences, genomic distribution of miRNA genes and identification of clusters, and isomiR characterization. Many of the miRNAs were regulated in a tissue/organ-specific manner.

Selection of reference genes for microRNA analysis associated to early stress response to handling and confinement in Salmo salar

MicroRNAs are key non-coding RNA molecules that play a relevant role in the regulation of gene expression through translational repression and/or transcript cleavage during normal development and physiological adaptation processes like stress. Quantitative reverse transcription polymerase chain reaction (RT-qPCR) has become the approach normally used to determine the levels of microRNAs. However, this approach needs the use of endogenous reference. An improper selection of endogenous references can result in confusing interpretation of data. The aim of this study was to identify and validate appropriate endogenous reference miRNA genes for normalizing RT-qPCR survey of miRNAs expression in four different tissues of Atlantic salmon, under handling and confinement stress conditions associated to early or primary stress response. Nine candidate reference normalizers, including microRNAs and nuclear genes, normally used in vertebrate microRNA expression studies were selected from literature, validated by RT-qPCR and analyzed by the algorithms geNorm and NormFinder. The results revealed that the ssa-miR-99-5p gene was the most stable overall and that ssa-miR-99-5p and ssa-miR-23a-5p genes were the best combination. Moreover, the suitability of ssa-miR-99-5p and ssa-miR-23a-5p as endogeneuos reference genes was demostrated by the expression analysis of ssa-miR-193-5p gene.

Identification of differentially expressed Atlantic salmon miRNAs responding to salmonid alphavirus (SAV) infection

Background

MicroRNAs (miRNAs) control multiple biological processes including the innate immune responses by negative post-transcriptional regulation of gene expression. As there were no studies on the role(s) of miRNAs in viral diseases in Atlantic salmon, we aimed to identify miRNAs responding to salmonid alphavirus (SAV) infection. Their expression were studied at different time points post infection with SAV isolates associated with different mortalities. Furthermore, the genome sequences of the identified miRNAs were analysed to reveal putative cis-regulatory elements, and, finally, their putative target genes were predicted.

Results

Twenty differentially expressed miRNAs (DE miRNAs) were identified. The expression of the majority of these increased post infection with maximum levels reached after the viral load were stabilized or decreasing. On the other hand, some miRNAs (e.g. the miRNA-21 family) showed decreased expression at the early time points post infection. There were significant differences in the temporal expression of individual miRNA associated with different SAV isolates. Target gene prediction in SAV responsive immune network genes showed that seventeen of the DE miRNAs could target 24 genes (e.g. IRF3, IRF7). Applying the Atlantic salmon transcriptome as input 28 more immune network genes were revealed as putative targets (e.g. IRF5, IRF4). The majority of the predicted target genes promote inflammatory response. The upstream sequences of the miRNA genes revealed a high density of cis-regulatory sequences known as binding sites for immune network transcription factors (TFs). A high expression in the late phase could therefore be due to increased transcription promoted by immune response activated TFs. Based on the in silico target predictions, we discuss their putative roles as early promotors or late inhibitors of inflammation. We propose that the differences in expressions associated with different SAV isolates could contribute to their differences in mortality rates.

Conclusions

This study represents the first steps in exploring miRNAs important in viral-host interaction in Atlantic salmon. We identified several miRNAs responding to SAV infection. Some likely to prohibit harmful inflammation while other may promote an early immune response. Their predicted functions need to be validated and further studied in functional assays to fully understand their roles in immune homeostasis.

Electronic supplementary material

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

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

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.

Estrogen regulation of gene expression in the teleost fish immune system

Elucidating the mechanisms of estrogens-induced immunomodulation in teleost fish is of great importance due to the observed worldwide continuing decrease in pristine environments. However, little is know about the immunotoxicological consequences of exposure to these chemicals in fish, or of the mechanisms through which these effects are mediated. In this review, we summarize the results showing estrogens (natural or synthetic) acting through estrogen receptors and regulating specific target genes, also through microRNAs (miRNAs), leading to modulation of the immune functioning. The identification and characterization of miRNAs will provide new opportunities for functional genome research on teleost immune system and can also be useful when screening for novel molecule biomarkers for environmental pollution.

The miR-33 gene is identified in a marine teleost: a potential role in regulation of LC-PUFA biosynthesis in Siganus canaliculatus

As the first marine teleost demonstrated to have the ability to biosynthesize long-chain polyunsaturated fatty acids (LC-PUFA) from C₁₈ PUFA precursors, rabbitfish Siganus canaliculatus provides a good model for studying the regulatory mechanisms of LC-PUFA biosynthesis in teleosts. Here the potential roles of miR-33 in such regulation were investigated. The miR-33 gene was identified within intron 16 of the gene encoding sterol regulatory element-binding protein 1 (Srebp1), an activator of LC-PUFA biosynthesis. Expression of miR-33 in rabbitfish tissues correlated with that of srebp1, while its expression in liver was highly responsive to ambient salinities and PUFA components, factors affecting LC-PUFA biosynthesis. Srebp1 activation promoted the expression of Δ4 and Δ6 Δ5 fatty acyl desaturases (Fad), key enzymes for LC-PUFA biosynthesis, accompanied by elevated miR-33 abundance in rabbitfish hepatocytes. miR-33 overexpression induced the expression of the two fad, but suppressed that of insulin-induced gene 1 (insig1), which encodes a repressor blocking Srebp proteolytic activation and has targeting sites of miR-33. These results indicated that miR-33, cooperating with Srebp1, may be involved in regulation of LC-PUFA biosynthesis by facilitating fad expression, probably through targeting insig1. To our knowledge, this is the first report of the participation of miR-33 in LC-PUFA biosynthesis in vertebrates.

Illumina Sequencing Reveals Aberrant Expression of MicroRNAs and Their Variants in Whitefish (Coregonus lavaretus) Liver after Exposure to Microcystin-LR

Molecular analyses show that challenging fish with microcystin-LR (MC-LR) causes perturbations of microRNA (miRNA) signaling. However, the significance and scope of these alterations is currently unknown. To address this issue, we studied miRNA gene expression in the liver of juvenile whitefish, C. lavaretus, during 28 days of exposure to a subacute dose of MC-LR (100 μg·kg^-1 body mass). Using genomic resources of Atlantic salmon (AGKD03), the mature miRNA library of Atlantic salmon (miRBase-21) and bioinformatics tools (sRNAbench), we discovered and annotated a total of 377 distinct mature miRNAs belonging to 93 families of evolutionary conserved miRNAs, as well as 24 novel mature miRNA candidates that were mapped to 14 distinct S. salar miRNA precursors. miRNA-Seq transcriptome profiling of liver tissues revealed differential miRNA expression in control and treated fish at 14 days (73 miRNAs were modulated) and at 28 days (83 miRNAs) of the treatment, subsequently validated by qPCR for nine selected differentially expressed miRNAs. Additional qPCR study confirmed the miRNA-Seq data and revealed consistent, aberrant miRNAs expression profile in the later phase of MC-LR hepatotoxicity (7–28 d). Functional annotation analysis revealed that the aberrantly expressed miRNAs have target genes involved in cytoskeletal remodeling, cell metabolism, cell cycle regulation and apoptosis; dysregulation of these processes in liver cells leads to cirrhosis and hepatocellular carcinoma in humans. To enable deeper insight into the molecular responses of liver cells in fish exposed to MC-LR, we expanded the miRNAome analysis by inclusion of miRNA variants (isomiRs) profiles, and we showed that the isomiR profiles of liver specific MiR122, and a few other miRNAs, correlated with MC-LR treatment. Given the importance of isomiRs for disease biology in mammals, we believe that further research focused on the miRNA isoforms will bring us closer to better understanding the molecular mechanisms of MC-LR hepatotoxicity.

MiRNA expression profile of ionizing radiation-induced liver injury in mouse using deep sequencing

In order to investigate the potential regulatory roles of microRNAs (miRNAs) in mouse response to ionizing radiation (IR), the small RNA libraries from liver tissues of mice with or without ionizing radiation (IR) were sequenced by high-throughput deep sequencing technology. A total of 270 miRNAs including 212 known and 58 potentially novel miRNAs were identified. Within these miRNAs, there were 48 miRNAs that were differentially expressed, including 27 known and 21 novel miRNAs. The results of quantitative RT-polymerase chain reaction (qRT-PCR) were in consistent with the sequencing analysis. Target gene prediction, function annotation, and pathway of the identified miRNAs were analyzed using RNAhybrid, miRanda software and Swiss-Prot, Gene Ontology (GO), Clusters of Orthologous Groups (COG), Kyoto Encyclopedia of Genes, and Genomes (KEGG) and non-redundant (NR) databases. These results should be useful to investigate the biological function of miRNAs under IR-induced liver injury.

Evolution and Distribution of Teleost myomiRNAs: Functionally Diversified myomiRs in Teleosts

Myosin heavy chain (MYH) genes belong to a multigene family, and the regulated expression of each member determines the physiological and contractile muscle properties. Among these, MYH6, MYH7, and MYH14 occupy unique positions in the mammalian MYH gene family because of their specific expression in slow/cardiac muscles and the existence of intronic micro(mi) RNAs. MYH6, MYH7, and MYH14 encode miR-208a, miR-208b, and miR-499, respectively. These MYH encoded miRNAs are designated as myomiRs because of their muscle-specific expression and functions. In mammals, myomiRs and host MYHs form a transcription network involved in muscle fiber-type specification; thus, genomic positions and expression patterns of them are well conserved. However, our previous studies revealed divergent distribution and expression of MYH14/miR-499 among teleosts, suggesting the unique evolution of myomiRs and host MYHs in teleosts. Here, we examined distribution and expression of myomiRs and host MYHs in various teleost species. The major cardiac MYH isoforms in teleosts are an intronless gene, atrial myosin heavy chain (amhc), and ventricular myosin heavy chain (vmhc) gene that encodes an intronic miRNA, miR-736. Phylogenetic analysis revealed that vmhc/miR-736 is a teleost-specific myomiR that differed from tetrapoda MYH6/MYH7/miR-208s. Teleost genomes also contain species-specific orthologs in addition to vmhc and amhc, indicating complex gene duplication and gene loss events during teleost evolution. In medaka and torafugu, miR-499 was highly expressed in slow/cardiac muscles whereas the expression of miR-736 was quite low and not muscle specific. These results suggest functional diversification of myomiRs in teleost with the diversification of host MYHs.

Integrated analysis of miRNA and mRNA expression profiles in tilapia gonads at an early stage of sex differentiation

Background

MicroRNAs (miRNAs) represent a second regulatory network that has important effects on gene expression and protein translation during biological process. However, the possible role of miRNAs in the early stages of fish sex differentiation is not well understood. In this study, we carried an integrated analysis of miRNA and mRNA expression profiles to explore their possibly regulatory patterns at the critical stage of sex differentiation in tilapia.

Results

We identified 279 pre-miRNA genes in tilapia genome, which were highly conserved in other fish species. Based on small RNA library sequencing, we identified 635 mature miRNAs in tilapia gonads, in which 62 and 49 miRNAs showed higher expression in XX and XY gonads, respectively. The predicted targets of these sex-biased miRNAs (e.g., miR-9, miR-21, miR-30a, miR-96, miR-200b, miR-212 and miR-7977) included genes encoding key enzymes in steroidogenic pathways (Cyp11a1, Hsd3b, Cyp19a1a, Hsd11b) and key molecules involved in vertebrate sex differentiation (Foxl2, Amh, Star1, Sf1, Dmrt1, and Gsdf). These genes also showed sex-biased expression in tilapia gonads at 5 dah. Some miRNAs (e.g., miR-96 and miR-737) targeted multiple genes involved in steroid synthesis, suggesting a complex miRNA regulatory network during early sex differentiation in this fish.

Conclusions

The sequence and expression patterns of most miRNAs in tilapia are conserved in fishes, indicating the basic functions of vertebrate miRNAs might share a common evolutionary origin. This comprehensive analysis of miRNA and mRNA at the early stage of molecular sex differentiation in tilapia XX and XY gonads lead to the discovery of differentially expressed miRNAs and their putative targets, which will facilitate studies of the regulatory network of molecular sex determination and differentiation in fishes.

Electronic supplementary material

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

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

Background

Atlantic cod (Gadus morhua) is among the economically most important species in the northern Atlantic Ocean and a model species for studying development of the immune system in vertebrates. MicroRNAs (miRNAs) are an abundant class of small RNA molecules that regulate fundamental biological processes at the post-transcriptional level. Detailed knowledge about a species miRNA repertoire is necessary to study how the miRNA transcriptome modulate gene expression. We have therefore discovered and characterized mature miRNAs and their corresponding miRNA genes in Atlantic cod. We have also performed a validation study to identify suitable reference genes for RT-qPCR analysis of miRNA expression in Atlantic cod. Finally, we utilized the newly characterized miRNA repertoire and the dedicated RT-qPCR method to reveal miRNAs that are highly expressed in certain organs.

Results

The discovery analysis revealed 490 mature miRNAs (401 unique sequences) along with precursor sequences and genomic location of the miRNA genes. Twenty six of these were novel miRNA genes. Validation studies ranked gmo-miR-17-1—5p or the two-gene combination gmo-miR25-3p and gmo-miR210-5p as most suitable qPCR reference genes. Analysis by RT-qPCR revealed 45 miRNAs with significantly higher expression in tissues from one or a few organs. Comparisons to other vertebrates indicate that some of these miRNAs may regulate processes like growth, lipid metabolism, immune response to microbial infections and scar damage repair. Three teleost-specific and three novel Atlantic cod miRNAs were among the differentially expressed miRNAs.

Conclusions

The number of known mature miRNAs was considerably increased by our identification of miRNAs and miRNA genes in Atlantic cod. This will benefit further functional studies of miRNA expression using deep sequencing methods. The validation study showed that stable miRNAs are suitable reference genes for RT-qPCR analysis of miRNA expression. Applying RT-qPCR we have identified several miRNAs likely to have important regulatory functions in particular organs.

Characterization of an extensive rainbow trout miRNA transcriptome by next generation sequencing

Background

MicroRNAs (miRNAs) have emerged as important post-transcriptional regulators of gene expression in a wide variety of physiological processes. They can control both temporal and spatial gene expression and are believed to regulate 30 to 70 % of the genes. Data are however limited for fish species, with only 9 out of the 30,000 fish species present in miRBase. The aim of the current study was to discover and characterize rainbow trout (Oncorhynchus mykiss) miRNAs in a large number of tissues using next-generation sequencing in order to provide an extensive repertoire of rainbow trout miRNAs.

Results

A total of 38 different samples corresponding to 16 different tissues or organs were individually sequenced and analyzed independently in order to identify a large number of miRNAs with high confidence. This led to the identification of 2946 miRNA loci in the rainbow trout genome, including 445 already known miRNAs. Differential expression analysis was performed in order to identify miRNAs exhibiting specific or preferential expression among the 16 analyzed tissues. In most cases, miRNAs exhibit a specific pattern of expression in only a few tissues. The expression data from sRNA sequencing were confirmed by RT-qPCR. In addition, novel miRNAs are described in rainbow trout that had not been previously reported in other species.

Conclusion

This study represents the first characterization of rainbow trout miRNA transcriptome from a wide variety of tissue and sets an extensive repertoire of rainbow trout miRNAs. It provides a starting point for future studies aimed at understanding the roles of miRNAs in major physiological process such as growth, reproduction or adaptation to stress. These rainbow trout miRNAs repertoire provide a novel resource to advance genomic research in salmonid species.

Electronic supplementary material

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