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

A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis)

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

Exploring transcriptional and post-transcriptional epigenetic regulation of crf and 11βhsd2 in rainbow trout brain during chronic social stress

Using dominance hierarchies in juvenile rainbow trout (Oncorhynchus mykiss) as a model of chronic social stress in fish, we explored whether epigenetic transcriptional and post-transcriptional mechanisms are involved in the gene expression of corticotropin-releasing factor (crf) and 11β-hydroxysteroid dehydrogenase (11βhsd2), key factors involved in the regulation of the endocrine stress axis response. In juvenile rainbow trout pairs, subordinate individuals display sustained elevation of circulating cortisol concentrations. Cortisol production is controlled by the hypothalamic-pituitary-interrenal (HPI) axis in fish and initiated by CRF release from the preoptic area (POA). Given that crf is modulated during chronic social stress, and that such stress has been implicated in the epigenetic regulation of crf in other taxa, we probed a role for epigenetic regulation of crf transcript abundance in chronically stressed rainbow trout. We also investigated the regulation of the cortisol-metabolising enzyme 11βhsd2 in the POA, which is upregulated in subordinates. The potential involvement of DNA methylation and microRNAs (miRNAs) in the regulation of crf transcript abundance was investigated during social stress in the POA of fish, as was the potential involvement of miRNAs in 11βhsd2 regulation. Although transcript abundances of crf were elevated in subordinate fish after 4 days, DNA methylation profiles within putative promoter sequences upstream of the crf gene were not significantly affected by chronic stress. An inverse relationship between crf and its predicted posttranscriptional regulator miR-103a-3p in the POA suggests that miRNAs may be involved in mediating the effects of chronic social stress on key components of the endocrine stress axis.

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.

Circulating MicroRNAs Indicative of Sex and Stress in the European Seabass (Dicentrarchus labrax): Toward the Identification of New Biomarkers

MicroRNAs (miRNAs) constitute a new category of biomarkers. Studies on miRNAs in non-mammalian species have drastically increased in the last few years. Here, we explored the use of miRNAs as potential, poorly invasive markers, to identify sex and characterize acute stress in fish. The European seabass (Dicentrarchus labrax) was chosen as a model because of its rapid response to stress and its specific sex determination system, devoid of sexual chromosomes. We performed a small RNA-sequencing analysis in the blood plasma of male and female European seabass (mature and immature) as well as in the blood plasma of juveniles submitted to an acute stress and sampled throughout the recovery period (at 0 h, 0.5 h, 1.5 h and 6 h). In immature individuals, both miR-1388-3p and miR-7132a-5p were up-regulated in females, while miR-499a-5p was more abundant in males. However, no miRNAs were found to be differentially expressed between sexes in the blood plasma of mature individuals. For the acute stress analysis, five miRNAs (miR-155-5p, miR-200a-3p, miR-205-1-5p, miR-143-3p, and miR-223-3p) followed cortisol production over time. All miRNAs identified were tested and validated by RT-qPCR on sequenced samples. A complementary analysis on the 3'UTR sequences of the European seabass allowed to predict potential mRNA targets, some of them being particularly relevant regarding stress regulation, e.g., the glucocorticoid receptor 1 and the mineralocorticoid receptor. The present study provides new avenues and recommendations on the use of miRNAs as biomarkers of sex or stress of the European seabass, with potential application on other fish species.

microRNA regulation of skin pigmentation in golden-back mutant of crucian carp from a rice-fish integrated farming system

BACKGROUND

MicroRNAs (miRNAs) are endogenous small non-coding RNAs (21-25 nucleotides) that act as essential components of several biological processes. Golden-back crucian carp (GBCrC, Carassius auratus) is a naturally mutant species of carp that has two distinct body skin color types (golden and greenish-grey), making it an excellent model for research on the genetic basis of pigmentation. Here, we performed small RNA (sRNA) analysis on the two different skin colors via Illumina sequencing.

RESULTS

A total of 679 known miRNAs and 254 novel miRNAs were identified, of which 32 were detected as miRNAs with significant differential expression (DEMs). 23,577 genes were projected to be the targets of 32 DEMs, primarily those involved in melanogenesis, adrenergic signaling in cardiomyocytes, MAPK signaling pathway and wnt signaling pathway by functional enrichment. Furthermore, we built an interaction module of mRNAs, proteins and miRNAs based on 10 up-regulated and 13 down-regulated miRNAs in golden skin. In addition to transcriptional destabilization and translational suppression, we discovered that miRNAs and their target genes were expressed in the same trend at both the transcriptional and translational levels. Finally, we discovered that miR-196d could be indirectly implicated in regulating melanocyte synthesis and motility in the skin by targeting to myh7 (myosin-7) gene through the luciferase reporter assay, antagomir silencing in vivo and qRT-PCR techniques.

CONCLUSIONS

Our study gives a systematic examination of the miRNA profiles expressed in the skin of GBCrC, assisting in the comprehension of the intricate molecular regulation of body color polymorphism and providing insights for C. auratus breeding research.

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.

Weighted Single-Step GWAS Identifies Genes Influencing Fillet Color in Rainbow Trout

The visual appearance of the fish fillet is a significant determinant of consumers' purchase decisions. Depending on the rainbow trout diet, a uniform bright white or reddish/pink fillet color is desirable. Factors affecting fillet color are complex, ranging from the ability of live fish to accumulate carotenoids in the muscle to preharvest environmental conditions, early postmortem muscle metabolism, and storage conditions. Identifying genetic markers of fillet color is a desirable goal but a challenging task for the aquaculture industry. This study used weighted, single-step GWAS to explore the genetic basis of fillet color variation in rainbow trout. We identified several SNP windows explaining up to 3.5%, 2.5%, and 1.6% of the additive genetic variance for fillet redness, yellowness, and whiteness, respectively. SNPs are located within genes implicated in carotenoid metabolism (β,β-carotene 15,15'-dioxygenase, retinol dehydrogenase) and myoglobin homeostasis (ATP synthase subunit β, mitochondrial (ATP5F1B)). These genes are involved in processes that influence muscle pigmentation and postmortem flesh coloration. Other identified genes are involved in the maintenance of muscle structural integrity (kelch protein 41b (klh41b), collagen α-1(XXVIII) chain (COL28A1), and cathepsin K (CTSK)) and protection against lipid oxidation (peroxiredoxin, superoxide dismutase 2 (SOD2), sestrin-1, Ubiquitin carboxyl-terminal hydrolase-10 (USP10)). A-to-G single-nucleotide polymorphism in β,β-carotene 15,15'-dioxygenase, and USP10 result in isoleucine-to-valine and proline-to-leucine non-synonymous amino acid substitutions, respectively. Our observation confirms that fillet color is a complex trait regulated by many genes involved in carotenoid metabolism, myoglobin homeostasis, protection against lipid oxidation, and maintenance of muscle structural integrity. The significant SNPs identified in this study could be prioritized via genomic selection in breeding programs to improve fillet color in rainbow trout.

Physical Enrichment Triggers Brain Plasticity and Influences Blood Plasma Circulating miRNA in Rainbow Trout (Oncorhynchus mykiss)

Physical enrichment is known to improve living conditions of fish held in farming systems and has been shown to promote behavioral plasticity in captive fish. However, the brain's regulatory-mechanism systems underlying its behavioral effects remain poorly studied. The present study investigated the impact of a three-month exposure to an enriched environment (EE vs. barren environment, BE) on the modulation of brain function in rainbow trout (Oncorhynchus mykiss) juveniles. Using high-throughput RT-qPCR, we assessed mRNA genes related to brain function in several areas of the trout brain. These included markers of cerebral activity and plasticity, neurogenesis, synaptogenesis, or selected neurotransmitters pathways (dopamine, glutamate, GABA, and serotonin). Overall, the fish from EE displayed a series of differentially expressed genes (neurotrophic, neurogenesis, and synaptogenesis markers) essentially localized in the telencephalon, which could underpin the beneficial effects of complexifying the environment on fish brain plasticity. In addition, EE significantly affected blood plasma c-miRNA signatures, as revealed by the upregulation of four c-miRNAs (miR-200b/c-3p, miR-203a-3p, miR-205-1a-5p, miR-218a-5p) in fish blood plasma after 185 days of EE exposure. Overall, we concluded that complexifying the environment through the addition of physical structures that stimulate and encourage fish to explore promotes the trout's brain function in farming conditions.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Differentially Expressed miRNAs and mRNAs in Regenerated Scales of Rainbow Trout (Oncorhynchus mykiss) under Salinity Acclimation

In order to explore the potential effects of salinity acclimation on bone metabolism of rainbow trout (Oncorhynchus mykiss), transcriptional information of regenerated scales under salinity acclimation (sea water, SW) was compared to those of fish under fresh water (FW) environments. According to the high-throughput sequencing results, a total of 2620 significantly differentially expressed genes (DEGs) were identified in the data of SW vs. FW. Compared with the FW group, six significantly downregulated and 44 significantly upregulated miRNAs were identified in the SW scales (p < 0.05). Furthermore, a total of 994 significantly differentially expressed target genes (DETGs) were identified from the 50 significantly differentially expressed miRNAs (DE miRNAs). Gene ontology analysis of the aforementioned DETGs was similar to the results of the differentially expressed genes (DEGs) obtained from mRNA-seq data, these genes were mainly related to ion metabolism. KEGG enrichment analysis of the DEGs and DETGs suggested that many significantly enriched pathways were related to the energy metabolism pathway.

A cross-species comparative approach to assessing multi- and transgenerational effects of endocrine disrupting chemicals

A wide range of chemicals have been identified as endocrine disrupting chemicals (EDCs) in vertebrate species. Most studies of EDCs have focused on exposure of both male and female adults to these chemicals; however, there is clear evidence that EDCs have dramatic effects when mature or developing gametes are exposed, and consequently are associated with in multigenerational and transgenerational effects. Several publications have reviewed such actions of EDCs in subgroups of species, e.g., fish or rodents. In this review, we take a holistic approach synthesizing knowledge of the effects of EDCs across vertebrate species, including fish, anurans, birds, and mammals, and discuss the potential mechanism(s) mediating such multi- and transgenerational effects. We also propose a series of recommendations aimed at moving the field forward in a structured and coherent manner.

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.

MicroRNAs and Related Cytokine Factors Quickly Respond in the Immune Response of Channel Catfish to Lipopolysaccharides and β-Glucan Stimulation

MicroRNAs (miRNAs) are well-known as powerful regulators of gene expression, with their potential to serve for immunology widely researched in mammals and birds but rarely in fishes. To better understand fish immunology behavior, we herein investigated nine immune-related miRNAs that were reported in other animals, as well as five related cytokine factors and lysozyme (LZM) in the liver, anterior kidney, and spleen of Channel Catfish Ictalurus punctatus after being stimulated by lipopolysaccharides (LPS) and β-glucan. We also predicated the potential targets of these miRNAs via bioinformatics and further investigated nine of them via quantitative real-time PCR. Results showed that expressions of the nine miRNAs were quickly changed in varying extent after stimulation by LPS, especially for miR-122, miR-142a, miR-155, and miR-223, which were significantly changed in spleen, and the same occurred for the LZM and three cytokine factors TNF-α, IFN-γ and TLR2. Compared with LPS, although most of the miRNAs and the cytokine genes were also affected by β-glucan, the extent of the effect was weak. Bioinformatics analysis revealed many immune-related targets of the miRNAs, with some of them reported by previous studies. For the nine investigated target genes, seven targets (77.8%) were significantly upregulated after the stimulation of LPS. It therefore can be inferred that the immune-related miRNAs, LZM, and cytokine factors elicited quick immune responses of Channel Catfish to LPS stimulation as in other animals, but the regulation mechanism of miRNAs might be complex and diverse. This research will contribute to a better understanding will support further immunology research in fishes.

Circulating miRNA repertoire as a biomarker of metabolic and reproductive states in rainbow trout

Background

Circulating miRNAs (c-miRNAs) are found in most, if not all, biological fluids and are becoming well-established non-invasive biomarkers of many human pathologies. However, their features in non-pathological contexts and whether their expression profiles reflect normal life history events have received little attention, especially in non-mammalian species. The aim of the present study was to investigate the potential of c-miRNAs to serve as biomarkers of reproductive and metabolic states in fish.

Results

The blood plasma was sampled throughout the reproductive cycle of female rainbow trout subjected to two different feeding regimes that triggered contrasting metabolic states. In addition, ovarian fluid was sampled at ovulation, and all samples were subjected to small RNA-seq analysis, leading to the establishment of a comprehensive miRNA repertoire (i.e., miRNAome) and enabling subsequent comparative analyses to a panel of RNA-seq libraries from a wide variety of tissues and organs. We showed that biological fluid miRNAomes are complex and encompass a high proportion of the overall rainbow trout miRNAome. While sharing a high proportion of common miRNAs, the blood plasma and ovarian fluid miRNAomes exhibited strong fluid-specific signatures. We further revealed that the blood plasma miRNAome significantly changed depending on metabolic and reproductive states. We subsequently identified three evolutionarily conserved muscle-specific miRNAs or myomiRs (miR-1-1/2-3p, miR-133a-1/2-3p, and miR-206-3p) that accumulated in the blood plasma in response to high feeding rates, making these myomiRs strong candidate biomarkers of active myogenesis. We also identified miR-202-5p as a candidate biomarker for reproductive success that could be used to predict ovulation and/or egg quality.

Conclusions

Together, these promising results reveal the high potential of c-miRNAs, including evolutionarily conserved myomiRs, as physiologically relevant biomarker candidates and pave the way for the use of c-miRNAs for non-invasive phenotyping in various fish species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01163-5.

MicroRNA expression profiling of sea perch brain cells reveals the roles of microRNAs in autophagy induced by RGNNV infection

Nervous necrosis virus (NNV) is one of the most destructive fish viruses and affects more than 120 marine and freshwater teleost species. However, the pathogenesis of NNV has not been made clear. MicroRNAs (miRNAs) play important roles in the regulation of viral infection. To understand the roles and regulation patterns of miRNAs in NNV infection, high-throughput sequencing was carried out in Lateolabrax japonicus brain (LJB) cells with or without red-spotted grouper NNV (RGNNV) infection at 12 and 24 hr. Here, we identified 59 known and 61 novel differentially expressed miRNAs (DE miRNAs) between mock and RGNNV-infected LJB cells. KEGG pathway analysis showed that the target genes of DE miRNAs were significantly enriched in immune-related signalling pathways, such as autophagy, mitophagy and TGF-beta signalling pathways. The expression patterns of four DE miRNAs (lja-miR-145, lja-miR-182, lja-miR-183 and lja-miR-187) were verified by qRT-PCR both in vivo and in vitro. We found that lja-miR-145 promoted RGNNV proliferation, while lja-miR-183 suppressed RGNNV proliferation. Furthermore, lja-miR-145 facilitated RGNNV-induced autophagy activation, whereas lja-miR-183 repressed autophagy in LJB cells as measured by LC3B-II/I and p62 protein levels. All these results indicate the involvement of lja-miR-145 and lja-miR-183 in RGNNV-induced autophagy. In conclusion, this study provides evidence for the important roles of miRNAs in NNV infection and a basis for uncovering the molecular regulation mechanism of NNV-induced autophagy.

Comprehensive analysis of miRNA-mRNA/lncRNA during gonadal development of triploid female rainbow trout (Oncorhynchus mykiss)

Chromosomal ploidy manipulation is one of the means to create excellent germplasm. Triploid fish could provide an ideal sterile model for searching of a underlying mechanism of abnormality in meiosis. The complete understanding of the coding and noncoding RNAs regulating sterility caused by meiosis abnormality is still not well understood. By high-throughput sequencing, we compared the expression profiles of gonadal mRNA, long non-coding RNA (lncRNA), and microRNA (miRNA) at three different developmental stages between the diploid (XX) and triploid (XXX) female rainbow trout. These stages were gonads before differentiation (65 days post fertilisation, dpf), at the beginning of morphological differences (180 dpf) and showing clear difference between diploids and triploids (600 dpf), respectively. A majority of differentially expressed (DE) RNAs were identified, and 22 DE mRNAs related to oocyte meiosis and homologous recombination were characterized. The predicted miRNA-mRNA/lncRNA networks of 3 developmental stages were constructed based on the target pairs of DE lncRNA-miRNA and DE mRNA-miRNA. According to the networks, meiosis-related gene of ccne1 was targeted by dre-miR-15a-5p_R + 1, and 6 targeted DE lncRNAs were identified. Also, qRT-PCR was performed to validate the credibility of the network. Overall, this study explored the potential interplay between coding and noncoding RNAs during the gonadal development of polyploid fish. The mRNA, lncRNA and miRNA screened in this study may be helpful to identify the functional elements regulating fertility of rainbow trout, which may provide reference for character improvement in aquaculture.

Responses of microRNA and predicted mRNA and enzymatic targets in liver of two salmonids (Oncorhynchus mykiss and Salvelinus fontinalis) following air exposure

The acute stress response is well-characterized, with rainbow trout as a teleost model for physiological and molecular responses. Air exposure, which stimulates an acute stress response, modulates liver microRNAs in rainbow trout; however, these highly conserved non-coding RNAs that bind to mRNA and repress translation, have never been measured in brook trout and it is unknown how miRNA expression responds following air exposure in this less studied salmonid. Our objective was to characterize the effects of air exposure on rainbow and brook trout liver miRNA expression, as well as the mRNA expression and enzyme activity that the miRNAs are predicted to target. Brook and rainbow trout were sampled pre- and 1-, 3-, and 24-h post- a three-minute air exposure. Plasma cortisol, glucose, and lactate were measured. Relative expression of miR-21a-5p, miR-143-3p, let-7a-5p and relative expression and enzyme activities of five predicted targets (pyruvate kinase, glucokinase, citrate synthase, cytochrome c oxidase, and catalase) were measured in liver. Rainbow and brook trout both had increases in plasma cortisol and lactate, while only rainbow trout had significant post-stress increases in plasma glucose. Furthermore, both trout species had increased miR-143-3p and miR-21a-5p relative expression 24-h post-stress. Four of the five enzymes measured had altered activity following stress. Brook trout miRNAs had inverse relative expression with relative catalase mRNA expression and cytochrome c oxidase enzyme activity, but no relationship was found in rainbow trout. Therefore, we have further characterized the transcriptional and enzymatic response to air exposure in two salmonids.

PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

A comparative analysis of heart microRNAs in vertebrates brings novel insights into the evolution of genetic regulatory networks

BACKGROUND

During vertebrate evolution, the heart has undergone remarkable changes that lead to morphophysiological differences in the fully formed heart of these species, such as chamber septation, heart rate frequency, blood pressure, and cardiac output volume. Despite these differences, the heart developmental process is guided by a core gene set conserved across vertebrates. Nonetheless, the regulatory mechanisms controlling the expression of genes involved in heart development and maintenance are largely uncharted. MicroRNAs (miRNAs) have been described as important regulatory elements in several biological processes, including heart biology. These small RNA molecules are broadly conserved in sequence and genomic context in metazoans. Mutations may occur in miRNAs and/or genes that contribute to the establishment of distinct repertoires of miRNA-target interactions, thereby favoring the differential control of gene expression and, consequently, the origin of novel phenotypes. In fact, several studies showed that miRNAs are integrated into genetic regulatory networks (GRNs) governing specific developmental programs and diseases. However, studies integrating miRNAs in vertebrate heart GRNs under an evolutionary perspective are still scarce.

RESULTS

We comprehensively examined and compared the heart miRNome of 20 species representatives of the five major vertebrate groups. We found 54 miRNA families with conserved expression and a variable number of miRNA families with group-specific expression in fishes, amphibians, reptiles, birds, and mammals. We also detected that conserved miRNAs present higher expression levels and a higher number of targets, whereas the group-specific miRNAs present lower expression levels and few targets.

CONCLUSIONS

Both the conserved and group-specific miRNAs can be considered modulators orchestrating the core and peripheral genes of heart GRNs of vertebrates, which can be related to the morphophysiological differences and similarities existing in the heart of distinct vertebrate groups. We propose a hypothesis to explain evolutionary differences in the putative functional roles of miRNAs in the heart GRNs analyzed. Furthermore, we present new insights into the molecular mechanisms that could be helping modulate the diversity of morphophysiology in the heart organ of vertebrate species.

High Throughput Sequencing of MicroRNA in Rainbow Trout Plasma, Mucus, and Surrounding Water Following Acute Stress

Circulating plasma microRNAs (miRNAs) are well established as biomarkers of several diseases in humans and have recently been used as indicators of environmental exposures in fish. However, the role of plasma miRNAs in regulating acute stress responses in fish is largely unknown. Tissue and plasma miRNAs have recently been associated with excreted miRNAs; however, external miRNAs have never been measured in fish. The objective of this study was to identify the altered plasma miRNAs in response to acute stress in rainbow trout (Oncorhynchus mykiss), as well as altered miRNAs in fish epidermal mucus and the surrounding ambient water. Small RNA was extracted and sequenced from plasma, mucus, and water collected from rainbow trout pre- and 1 h-post a 3-min air stressor. Following small RNA-Seq and pathway analysis, we identified differentially expressed plasma miRNAs that targeted biosynthetic, degradation, and metabolic pathways. We successfully isolated miRNA from trout mucus and the surrounding water and detected differences in miRNA expression 1-h post air stress. The expressed miRNA profiles in mucus and water were different from the altered plasma miRNA profile, which indicated that the plasma miRNA response was not associated with or immediately reflected in external samples, which was further validated through qPCR. This research expands understanding of the role of plasma miRNA in the acute stress response of fish and is the first report of successful isolation and profiling of miRNA from fish mucus or samples of ambient water. Measurements of miRNA from plasma, mucus, or water can be further studied and have potential to be applied as non-lethal indicators of acute stress in fish.

Meta-analysis of differentially-regulated hepatic microRNAs identifies candidate post-transcriptional regulation networks of intermediary metabolism in rainbow trout

MicroRNAs (miRNAs) are small non-coding RNAs which act as post-transcriptional regulators by decreasing targeted mRNA translation and stability. Principally targeting small 3' UTR elements of protein-coding mRNAs through complementary base-pairing, miRNAs are promiscuous regulators of the transcriptome. While potent roles for hepatic miRNAs in the regulation of energy metabolism have emerged in rodent models, comparative roles in other vertebrates remain largely unexplored. Indeed, while several miRNAs are deeply conserved among vertebrates, the acquisition of lineage- and species-specific miRNAs, as well as the rewiring between miRNA-mRNA target relationships beg the question of regulatory and functional conservation and innovation of miRNAs and their targets involved in energy metabolism. Here we provide a meta-analysis of differentially expressed hepatic miRNAs in rainbow trout, a scientifically and economically important teleost species with a 'glucose-intolerant' phenotype. Following exposure to nutritional and social context-dependent metabolic challenges, we analyzed differential miRNA expression from small-RNA-sequencing datasets generated with a consistent bioinformatics pipeline in conjunction with an in silico target prediction of metabolic transcripts and pathways. We provide evidence for evolutionary conserved (let-7, miRNA-27 family) and rewired (miRNA-30 family, miRNA-152, miRNA-722) miRNA-metabolic target gene networks in the context of the salmonid genome. These findings represent important first steps in our understanding of the comparative regulation and function of hepatic miRNAs in rainbow trout energy metabolism. We propose that the identified miRNA families should be prioritized for future comparative functional investigation in the context of hepatic energy- and glucose metabolism in rainbow trout.

Non-coding RNA Expression Patterns of Two Different Teleost Gonad Maturation Stages

Non-coding RNAs (ncRNAs) are involved in several different regulatory pathways including reproduction. In teleost fish, efficacious reproduction is heavily dependent on the completion of the reproductive cycle. The presence of ncRNA, however, and their expression dynamics and putative regulatory role in mature and immature gonads have not yet been extensively explored. Therefore, the abundance of ncRNAs in mature and immature female sharpsnout seabream (Diplodus puntazzo) was investigated. The sharpsnout seabream is a rudimentary hermaphrodite which, in captivity, displays dysfunctions in the gonad maturation process. Our analyses revealed a gonad specific read length distribution with two main peaks representing miRNAs (21-26 nt) and PIWI RNA (27-34 nt). Besides, distinct expression patterns for several ncRNA biotypes including microRNAs (miRNAs), PIWI RNAs (piRNAs), and ribosomal RNAs (rRNAs) were detected. Identified miRNA accounted to 938, corresponding to ~ 13% of obtained transcripts. Among the differential expressed ncRNAs, 10 (~ 7%) were annotated as miRNA, out of which 2 were found in higher abundance in immature gonads (miR-125c and miR-24) and 8 (miR-451, miR-7a, miR-122-1, miR190a, miR129, ENSGACT00000029608, ENSGACT00000029489, and ENSGACT00000029667) were found to be higher expressed in mature gonads. Putative miRNA targets, including long non-coding RNAs (lncRNAs) and genes, are proposed. Target genes are involved in several processes of fish oocyte development, such as steroidogenesis, proteolysis, and apoptosis, and may explain hormone regulation. This study demonstrates a gonad maturation biased ncRNA profile which in turn may support the role of ncRNAs in ovarian physiology and reproductive performance of fish, stressing the specific function of each RNA biotype in oocyte development.

An estradiol-17β/miRNA-26a/cyp19a1a regulatory feedback loop in the protogynous hermaphroditic fish, Epinephelus coioides

Cyp19a1a is a key gene responsible for the production of estradiol-17β (E₂), the main functional estrogen and a major downstream regulator of reproduction in teleost fish. It is widely known that CYP19 gene expression, aromatase activity, and E₂ production can influence gonadal differentiation and sex reversal in teleost fish, but the feedback mechanisms whereby E₂ regulates cyp19a1a remain poorly understood, especially regarding the potential roles of endogenous small RNA molecules (miRNAs). Here, we identified miR-26a-5p as a regulatory factor of its predicted target gene (cyp19a1a). In vitro and in vivo studies showed that miR-26a-5p can decrease cyp19a1a expression. Furthermore, high doses of E₂ act as a repressor of miR-26a-5p. This study proposes a regulatory feedback loop whereby E₂ regulates cyp19a1a through miR-26a-5p, and suggests that this positive feedback is an important aspect of the control of E₂ production.

Potential immune regulatory role of miR-146a upon Aeromonas hydrophila and Edwardsiella piscicida infections in zebrafish

This study was aimed to understand the expression of miR-146a in zebrafish (Danio rerio) and its role in regulating immune responses during Aeromonas hydrophila and Edwardsiella piscicida infections. The miR-146a expression was observed from the 1-h post fertilization (hpf) stage and gradually increased up to the early larval stage of zebrafish. The ubiquitous expression of miR-146a was detected in all tested tissues, with the highest level in gills. The expression of miR-146a was significantly increased in larvae when exposed to E. piscicida infection at 24 and 48 h post exposure (hpe). Intraperitoneally (i.p.) injected A. hydrophila and E. piscicida into adult zebrafish showed significant upregulation of miR-146a in gills. Furthermore, immune-related genes, toll-like receptor, tlr-4, transducing signaling pathway molecules, traf-6 and myd88 (bacteria-infected larvae and adults), transcription factor relA and mcp-1b (bacteria-infected adults), pro-inflammatory, il-6 (A. hydrophila-exposed larvae) and mmp-9 (bacteria-exposed larvae) were significantly repressed. In contrast, il-1β, tnf-α, cxcl-18b, and ccl-34a.4 were induced in both bacteria-challenged larvae and adults. Based on the results, it is suggested that endogenous miR-146a could act as an infection inducible miRNA in zebrafish upon A. hydrophila and E. piscicida infections; also, it could potentially regulate the immune responses in zebrafish.

Acute air exposure modulates the microRNA abundance in stress responsive tissues and circulating extracellular vesicles in rainbow trout (Oncorhynchus mykiss)

The hypothalamic-pituitary-interrenal axis is an important regulator of stress and metabolism in teleosts. Cortisol is secreted by the head kidney where it increases gluconeogenesis in the liver to increase circulating glucose levels. MicroRNAs (miRNAs) are small, non-coding RNA molecules that bind to the 3' untranslated region of specific mRNA to regulate their expression. MicroRNAs can also be secreted into circulation by association with extracellular vesicles (EVs) where they can influence the phenotype of other tissues. In this study, adult rainbow trout were exposed to a 3-minute acute air stress and allowed to recover for 1-, 3-, or 24-h to determine how miRNAs were altered. MicroRNAs measured in this study were chosen based on their high relative abundance in tissues that drive the stress response (miR-21a-3p, let-7a-5p, miR-143-3p) or their role in regulating DNA methylation (miR-29a-3p). In general, miRNAs increased in circulating EVs during the recovery period while decreasing in head kidney and liver at the same timepoints. Predicted targets for these miRNAs were analyzed using KEGG and DAVID functional enrichment analysis. Pathways involved in metabolism and cell signaling were predicted to be upregulated. Future studies can use these results to investigate how pathways are regulated after stress. Overall, our results indicate that miRNAs are regulated during teleost stress responses and could be supporting the cortisol-mediated changes that occur.

Profiling the rainbow trout hepatic miRNAome under diet-induced hyperglycemia

Carnivorous rainbow trout exhibit prolonged postprandial hyperglycemia when fed a diet exceeding 20% carbohydrate content. This poor capacity to utilize carbohydrates has led to rainbow trout being classified as "glucose-intolerant" (GI). The metabolic phenotype has spurred research to identify the underlying cellular and molecular mechanisms of glucose intolerance, largely because carbohydrate-rich diets provide economic and ecological advantages over traditionally used fish meal, considered unsustainable for rainbow trout aquaculture operations. Evidence points to a contribution of hepatic intermediary carbohydrate and lipid metabolism, as well as upstream insulin signaling. Recently, microRNAs (miRNAs), small noncoding RNAs acting as negative posttranscriptional regulators affecting target mRNA stability and translation, have emerged as critical regulators of hepatic control of glucose-homeostasis in mammals, revealing that dysregulated hepatic miRNAs might play a role in organismal hyperglycemia in metabolic disease. To determine whether hepatic regulatory miRNA networks may contribute to GI in rainbow trout, we induced prolonged postprandial hyperglycemia in rainbow trout by using a carbohydrate-rich diet and profiled genome-wide hepatic miRNAs in hyperglycemic rainbow trout compared with fasted trout and trout fed a diet devoid of carbohydrates. Using small RNA next-generation sequencing and real-time RT-PCR validation, we identified differentially regulated hepatic miRNAs between these groups and used an in silico approach to predict bona fide mRNA targets and enriched pathways. Diet-induced hyperglycemia resulted in differential regulation of hepatic miRNAs compared with fasted fish. Some of the identified miRNAs, such as miRNA-27b-3p and miRNA-200a-3p, are known to be responsive to hyperglycemia in the liver of hyperglycemic glucose-tolerant fish and mammals, suggesting an evolutionary conserved regulation. Using Gene Ontology term-based enrichment analysis, we identify intermediate carbohydrate and lipid metabolism and insulin signaling as potential targets of posttranscriptional regulation by hyperglycemia-regulated miRNAs and provide correlative expression analysis of specific predicted miRNA-target pairs. This study identifies hepatic miRNAs in rainbow trout that exhibit differential postprandial expression in response to diets with different carbohydrate content and predicts posttranscriptionally regulated target mRNAs enriched for pathways involved in glucoregulation. Together, these results provide a framework for testable hypotheses of functional involvement of specific hepatic miRNAs in GI in rainbow trout.

Social status regulates the hepatic miRNAome in rainbow trout: Implications for posttranscriptional regulation of metabolic pathways

Juvenile rainbow trout develop social hierarchies when held in dyads, and the development of socially subordinate (SS) and social dominance (SD) phenotypes in this context has been linked to specific changes in the hepatic energy metabolism of all major macronutrients. Following our recently reported finding that transcript abundance of drosha, a key component of the microRNA (miRNA) biogenesis pathway, is increased in paired juvenile rainbow trout irrespective of social status compared to socially isolated (SI) controls, we here determined global changes of the hepatic miRNA pathway genes in detail at the transcript and protein level. Both SD and SS rainbow trout exhibited increased Ago2 protein abundance compared to SI rainbow trout, suggesting that hepatic miRNA function is increased in rainbow trout maintained in dyads. Given the well-described differences in hepatic intermediary metabolism between SD and SS rainbow trout, and the important role of miRNAs in the posttranscriptional regulation of metabolic pathways, we also identified changes in hepatic miRNA abundance between SS and SD rainbow trout using small RNA next generation sequencing. We identified a total of 24 differentially regulated miRNAs, with 15 miRNAs that exhibited increased expression, and 9 miRNAs that exhibited decreased expression in the liver of SS trout compared to SD trout. To identify potential miRNA-dependent posttranscriptional regulatory pathways important for social status-dependent regulation of hepatic metabolism in rainbow trout, we used an in silico miRNA target prediction and pathway enrichment approach. We identified enrichment for pathways related to metabolism of carbohydrates, lipids and proteins in addition to organelle-specific processes involved in energy metabolism, especially mitochondrial fusion and fission. Select predicted miRNA-mRNA target pairs within these categories were quantitatively analyzed by real-time RT-PCR to validate candidates for future studies that will probe the functional metabolic roles of specific hepatic miRNAs in the development of SD and SS metabolic phenotypes.

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.

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

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

New Insights on Intermediary Metabolism for a Better Understanding of Nutrition in Teleosts

The rapid development of aquaculture production throughout the world over the past few decades has led to the emergence of new scientific challenges to improve fish nutrition. The diet formulations used for farmed fish have been largely modified in the past few years. However, bottlenecks still exist in being able to suppress totally marine resources (fish meal and fish oil) in diets without negatively affecting growth performance and flesh quality. A better understanding of fish metabolism and its regulation by nutrients is thus mandatory. In this review, we discuss four fields of research that are highly important for improving fish nutrition in the future: ( a) fish genome complexity and subsequent consequences for metabolism, ( b) microRNAs (miRNAs) as new actors in regulation of fish metabolism, ( c) the role of autophagy in regulation of fish metabolism, and ( d) the nutritional programming of metabolism linked to the early life of fish.

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.

Circulating miRNA measurements are reflective of cholesterol-based changes in rainbow trout (Oncorhynchus mykiss)

MicroRNAs (miRNAs) are a class of small non-coding RNAs which are known to posttranscriptionally regulate the expression of most genes in both animals and plants. Meanwhile, studies have shown that numbers of miRNAs are present in body fluids including the plasma. Despite the mode of action of these circulating miRNAs still remains unknown, they have been found to be promising biomarkers for disease diagnosis, prognosis and response to treatment. In order to evaluate the potential of miRNAs as non-invasive biomarkers in aquaculture, a time-course experiment was implemented to investigate the postprandial regulation of miRNAs levels in liver and plasma as well as the hepatic expression of genes involved in cholesterol metabolism. We showed that miR-1, miR-33a, miR-122, miR-128 and miR-223 were expressed in the liver of rainbow trout and present at detectable level in the plasma. We also demonstrated that hepatic expression of miR-1, miR-122 and miR-128 were regulated by feed intake and reached their highest levels 12 hours after the meal. Interestingly, we observed that circulating levels of miR-128 and miR-223 are subjected to postprandial regulations similar to that observed in their hepatic counterparts. Statistical correlations were observed between liver and plasma for miR-128 and miR-223 and between hepatic and circulating miR-122, miR-128 and miR-223 and expression of genes related to cholesterol synthesis and efflux or glucose phosphorylation. These results demonstrated that circulating miR-122, miR-128 and miR-223 are potential biomarkers of cholesterol metabolism in rainbow trout.

Identification and analysis of differentially expressed microRNAs in rainbow trout (Oncorhynchus mykiss) responding to infectious hematopoietic necrosis virus infection

MicroRNAs (miRNAs) are a class of regulators essential for numerous biological processes. Infectious hematopoietic necrosis virus (IHNV) is one of the most important viral pathogens in salmon and trout. In this study, the miRNA expression profiles of rainbow trout upon IHNV infection were explored. In total, 392 known miRNAs and 936 novel miRNAs were identified. Twelve known and 13 novel miRNAs were differentially expressed between infected and uninfected fish. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses revealed that certain miRNA target genes were associated with biological regulation, the immune system, and signal transduction. In addition, over- and suppressed expression of miR-146a-3p, miR-155-5p, miR-216a-5p, and miR-499b-5p could respectively increase and decrease viral gene expression in cells and viral titers. MiR-146a-3p and miR-216a-5p inhibited the expression of type-I IFN and the Mx1 gene induced by IHNV. These results provide preliminary insights into the IHNV-host interactions mediated by miRNAs.

MiR-202 controls female fecundity by regulating medaka oogenesis

Female gamete production relies on coordinated molecular and cellular processes that occur in the ovary throughout oogenesis. In fish, as in other vertebrates, these processes have been extensively studied both in terms of endocrine/paracrine regulation and protein expression and activity. The role of small non-coding RNAs in the regulation of animal reproduction remains however largely unknown and poorly investigated, despite a growing interest for the importance of miRNAs in a wide variety of biological processes. Here, we analyzed the role of miR-202, a miRNA predominantly expressed in male and female gonads in several vertebrate species. We studied its expression in the medaka ovary and generated a mutant line (using CRISPR/Cas9 genome editing) to determine its importance for reproductive success with special interest for egg production. Our results show that miR-202-5p is the most abundant mature form of the miRNA and that it is expressed in granulosa cells and in the unfertilized egg. The knock out (KO) of mir-202 gene resulted in a strong phenotype both in terms of number and quality of eggs produced. Mutant females exhibited either no egg production or produced a dramatically reduced number of eggs that could not be fertilized, ultimately leading to no reproductive success. We quantified the size distribution of the oocytes in the ovary of KO females and performed a large-scale transcriptomic analysis approach to identified dysregulated molecular pathways. Together, cellular and molecular analyses indicate that the lack of miR-202 impairs the early steps of oogenesis/folliculogenesis and decreases the number of large (i.e. vitellogenic) follicles, ultimately leading to dramatically reduced female fecundity. This study sheds new light on the regulatory mechanisms that control the early steps of follicular development, including possible targets of miR-202-5p, and provides the first in vivo functional evidence that a gonad-predominant microRNA may have a major role in female reproduction.

The role of small non-coding RNAs in the regulation of animal reproduction remains poorly investigated, despite a growing interest for the importance of miRNAs in a wide variety of biological processes. Here, we analyzed the role of miR-202, a miRNA predominantly expressed in gonads in vertebrate. We studied its expression in the medaka ovary and knocked out the mir-202 gene to study its importance for reproductive success. We showed that the lack of miR-202 results in the sterility of both females and males. In particular, it led to a drastic reduction of both the number and the quality of eggs produced by females. Mutant females exhibited either no egg production or produced a drastically reduced number of eggs that could not be fertilized, ultimately leading to no reproductive success. Quantitative histological and molecular analyses indicated that mir-202 KO impairs oocyte development and is also associated with the dysregulation of many genes that are critical for reproduction. This study sheds new light on the regulatory mechanisms that control oogenesis, including possible targets of miR-202-5p, and provides the first in vivo functional evidence that a gonad-predominant microRNA may have a major role in female reproduction.

Identification and differential regulation of microRNAs during thyroid hormone-dependent metamorphosis in Microhyla fissipes

BACKGROUND

Anuran metamorphosis, which is obligatorily initiated and sustained by thyroid hormone (TH), is a dramatic example of extensive morphological, biochemical and cellular changes occurring during post-embryonic development. Thus, it provides an ideal model to understand the actions of the hormone and molecular mechanisms underlying these developmental and apoptotic processes. In addition to transcriptional factors, microRNAs (miRNAs) play key roles in diverse biological processes via post-transcriptional repression of mRNAs. However, the possible role of miRNAs in anuran metamorphosis is not well understood. Screening and identification of TH-responding miRNAs are required to reveal the integrated regulatory mechanisms of TH during metamorphosis. Given the specific role of TRs during M. fissipes metamorphosis and the characteristics of M. fissipes as an ideal model, Illumina sequencing technology was employed to get a full scope of miRNA in M. fissipes metamorphosis treated by T3.

RESULTS

Morphological and histological analysis revealed that 24 h T3 treatment M. fissipes tadpoles resembled that at the climax of natural metamorphosis. Thus, small RNA libraries were constructed from control and 24 h T3 treatment groups. A total of 164 conserved miRNAs and 36 predicted novel miRNAs were characterized. Furthermore, 5' first and ninth nucleotides of miRNAs were significantly enriched in U in our study. In all, 21 miRNAs were differentially expressed between the T3 and control groups (p < 0.01). A total of 10,206 unigenes were identified as target genes of these differentially expressed miRNAs. KEGG pathway analysis indicated that the most overrepresented miRNA target genes were enriched in the "PI3k-Akt signaling pathway". In addition, a network associated with the TH signaling pathway provides an opportunity to further understand the complex biological processes that occur in metamorphosis.

CONCLUSIONS

We identified a large number of miRNAs during M. fissipes metamorphosis, and 21 of them were differentially expressed in the two groups that represented two different metamorphic stages. These miRNAs may play important roles during metamorphosis. The study gives us clues for further studies of the mechanisms of anuran metamorphosis and provides a model to study the mechanism of TH-affected biological processes in humans.

Dietary methionine restriction: Effects on glucose tolerance, lipid content and micro-RNA composition in the muscle of rainbow trout

Lean muscle mass plays an important role in overall health, as altered skeletal muscle metabolism can impact both the incidence and prevention of conditions related to metabolic health. Intriguingly, dietary methionine restriction (MR) has been shown to ameliorate this phenotype over time potentially through mechanisms related to changes in myogenic precursor cell (MPC) differentiation status. Recently the role of micro-RNAs (miRs) in regulating the expression of muscle specific transcription factors myoD and myogenin as well as signaling molecules involved in skeletal muscle differentiation has been reported in vitro. We performed an 8week feeding trial to determine if MR in vivo could alter miR abundance as well as change metabolic markers. Results show changes in muscle miR abundance for miR-133a at 4weeks with no significant difference seen in miR-210 or miR-206. After 8weeks of MR feeding fish demonstrated increased clearance of glucose, increased fat accumulation in the liver, and decreased fat accumulation in the muscle. These data demonstrate conservation of MR effects on fish metabolism, and suggest, for the first time, that miR-133a might play a role in tissue response to MR.

Crosstalk among lncRNAs, microRNAs and mRNAs in the muscle 'degradome' of rainbow trout

In fish, protein-coding and noncoding genes involved in muscle atrophy are not fully characterized. In this study, we characterized coding and noncoding genes involved in gonadogenesis-associated muscle atrophy, and investigated the potential functional interplay between these genes. Using RNA-Seq, we compared expression pattern of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs of atrophying skeletal muscle from gravid females and control skeletal muscle from age-matched sterile individuals. A total of 852 mRNAs, 1,160 lncRNAs and 28 microRNAs were differentially expressed (DE) between the two groups. Muscle atrophy appears to be mediated by many genes encoding ubiquitin-proteasome system, autophagy related proteases, lysosomal proteases and transcription factors. Transcripts encoding atrogin-1 and mir-29 showed exceptional high expression in atrophying muscle, suggesting an important role in bulk muscle proteolysis. DE genes were co-localized in the genome with strong expression correlation, and they exhibited extensive 'lncRNA-mRNA', 'lncRNA-microRNA', 'mRNA-microRNA' and 'lncRNA-protein' physical interactions. DE genes exhibiting potential functional interactions comprised the highly correlated 'lncRNA-mRNA-microRNA' gene network described as 'degradome'. This study pinpoints extensive coding and noncoding RNA interactions during muscle atrophy in fish, and provides valuable resources for future mechanistic studies.

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.

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.

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.

Identification and expression analysis of microRNAs in medaka gonads

Gonad development is a highly regulated, coordinated biological process and increasing evidences have indicated that microRNA (miRNA) may be involved in this dynamic program. Medaka (Oryzias latipes) is a good model for reproductive research as it has distinct sex determining genes, however, research in gonadal miRNAs is lacked. In this study, two small RNA libraries from the ovaries and testes were constructed and sequenced. A total of 285 conserved and 388 novel miRNAs were obtained, among which 142 mature miRNAs were significantly (> two-fold change) up or down regulated in the testis compared to the ovary. Quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR) analysis showed that miR-430c, miR-26a and miR-202-5p were expressed in a gonad-specific or sex-biased pattern. Fluorescence in situ hybridization (FISH) indicated that miR-202-5p was present throughout spermatogenesis and was only detected at the early stages of oogenesis, this sex biased expression pattern suggested that miR-202-5p might be a crucial candidate in male differentiation and development. Our study provides the repertoire, a comprehensive annotation of miRNAs from gonads and a reference for functional studies of miRNAs in medaka.

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.

Rainbow trout exposed to benzo[a]pyrene yields conserved microRNA binding sites in DNA methyltransferases across 500 million years of evolution

The objective of this study was to examine the regulation of DNA methylation following acute (24 h) and prolonged (14 d) exposure to low (1 ng/L) and high (10 ng/L) benzo[a]pyrene. However, with the recent release of the rainbow trout genome, we were able to conduct a more detailed analysis regarding the regulation of the enzymes involved in DNA methylation; DNA methyltransferases (DNMTs). Bioinformatic approaches were used to identify candidate microRNA (miRNA) that potentially bind to the DNMT1 and DNMT3a 3′UTR. Results indicated a significant decrease in global methylation in both liver and muscle, with an associated decrease in DNA methyltransferase activity and DNMT3a transcript abundance. There was a significant increase in one specific candidate miRNA (miR29a) that was predicted to bind to DNMT3a. Taking a comparative genomics approach, the binding sites of miR29a to the DNMT3a 3′UTR was compared across species, spanning fish to mammals, and revealed a highly conserved binding motif that has been maintained since the vertebrate ancestor, approximately 500 million years ago. This research establishes that miRNA act as an essential mediator between the environment and DNA methylation patterns via DNMTs, which is further confirmed by a genomic regulatory mechanism that has been deeply conserved throughout evolution.

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.

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.

Regulation of genes related to cholesterol metabolism in rainbow trout (Oncorhynchus mykiss) fed a plant-based diet

When compared with fish meal and fish oil, plant ingredients differ not only in their protein content and amino acid and fatty acid profiles but are also devoid of cholesterol, the major component of cell membrane and precursor of several bioactive compounds. Based on these nutritional characteristics, plant-based diets can affect fish physiology and cholesterol metabolism. To investigate the mechanisms underlying cholesterol homeostasis, rainbow trout were fed from 1 g body wt for 6 mo with a totally plant-based diet (V), a marine diet (M), and a marine-restricted diet (MR), with feed intake adjusted to that of the V group. The expression of genes involved in cholesterol synthesis, esterification, excretion, bile acid synthesis, and cholesterol efflux was measured in liver. Results showed that genes involved in cholesterol synthesis were upregulated in trout fed the V diet, whereas expression of genes related to bile acid synthesis ( cyp7a1) and cholesterol elimination ( abcg8) were reduced. Feeding trout the V diet also enhanced the expression of srebp-2 while reducing that of lxrα and miR-223. Overall, these data suggested that rainbow trout coped with the altered nutritional characteristics and absence of dietary cholesterol supply by increasing cholesterol synthesis and limiting cholesterol efflux through molecular mechanisms involving at least srebp-2, lxrα, and miR-223. However, plasma and body cholesterol levels in trout fed the V diet were lower than in fish fed the M diet, raising the question of the role of cholesterol in the negative effect of plant-based diet on growth.

RNA-Seq reveals MicroRNA expression signature and genetic polymorphism associated with growth and muscle quality traits in rainbow trout

The role of microRNA expression and genetic variation in microRNA-binding sites of target genes on growth and muscle quality traits is poorly characterized. We used RNA-Seq approach to investigate their importance on 5 growth and muscle quality traits: whole body weight (WBW), muscle yield, muscle crude-fat content, muscle shear force and whiteness. Phenotypic data were collected from 471 fish, representing 98 families (~5 fish/family) from a growth-selected line. Muscle microRNAs and mRNAs were sequenced from 22 families showing divergent phenotypes. Ninety microRNAs showed differential expression between families with divergent phenotypes, and their expression was strongly associated with variation in phenotypes. A total of 204 single nucleotide polymorphisms (SNPs) present in 3′ UTR of target genes either destroyed or created novel illegitimate microRNA target sites; of them, 78 SNPs explained significant variation in the aforementioned 5 muscle traits. Majority of the phenotype-associated SNPs were present in microRNA-binding sites of genes involved in energy metabolism and muscle structure. These findings suggest that variation in microRNA expression and/or sequence variation in microRNA binding sites in target genes play an important role in mediating differences in fish growth and muscle quality phenotypes.