Discovery and functional characterization of novel miRNAs in the marine medaka Oryzias melastigma

Comparative analysis of miRNAs and mRNAs in large yellow croaker head kidney cells (LYCK) provided novel insights into the redox regulation of fish

Reactive oxygen species (ROS) over-production and oxidative stress resulted from climate change and environmental pollution seriously endangered global fish populations and healthy development of marine aquaculture. Peroxiredoxins (Prxs), a highly conserved family of thiol-specific antioxidants, can mitigate ROS and protect cells from oxidative stress. We previously demonstrated that large yellow croaker PrxIV (LcPrxIV) could not only regulate the pro-inflammatory responses, but also scavenge ROS. However, the underlying mechanism how LcPrxIV regulated immune response and redox homeostasis remains unknown. MicroRNAs (miRNAs) are non-coding RNAs that play important roles in the regulation of various biological processes. In this study, mRNA and miRNA expression profiles from LYCK-pcDNA3.1 and LYCK-PrxIV cells, with or without oxidative stress stimulated by H2O2 were evaluated using high-throughput sequencing. A series of differentially expressed miRNAs (DEMs) and differentially expressed genes (DEGs), as well as DEM-DEG pairs were identified from each two-group comparison, respectively. GO and KEGG functional analyses indicated that most significant DEGs were associated with signaling pathways related to oxidative stress and immune response. Subsequent DEM-DEG interaction analysis revealed that miR-731 and miR-1388 may be involved in both redox regulation and immune response via synergistic effect with LcPrxIV. Interestingly, miR-731 could regulate the expression of different down-stream DEGs under different stimulations of LcPrxIV over-expression, H2O2, or both. Moreover, miR-731 could cause the DEG, γ-glutamyl hydrolase (GGH), to be expressed in opposite ways under different stimulations. On the other hand, the expression of miR-1388 could be negatively or positively regulated under the stimulation of LcPrxIV over-expression with or without oxidative stress, thus regulating gene expression of different mRNAs. Based on these results, we speculate that LcPrxIV may participate in immune response or redox regulation by regulating the expression of different down-stream genes through controlling the expression level of a certain miRNA or by regulating the varieties of expressed miRNAs.

Transgenerational bone toxicity in F3 medaka (Oryzias latipes) induced by ancestral benzo[a]pyrene exposure: Cellular and transcriptomic insights

Benzo[a]pyrene (BaP), a ubiquitous pollutant, raises environmental health concerns due to induction of bone toxicity in the unexposed offspring. Exposure of F0 ancestor medaka (Oryzias latipes) to 1 µg/L BaP for 21 days causes reduced vertebral bone thickness in the unexposed F3 male offspring. To reveal the inherited modifications, osteoblast (OB) abundance and molecular signaling pathways of transgenerational BaP-induced bone thinning were assessed. Histomorphometric analysis showed a reduction in OB abundance. Analyses of the miRNA and mRNA transcriptomes revealed the dysregulation of Wnt signaling (frzb/ola-miR-1-3p, sfrp5/ola-miR-96-5p/miR-455-5p) and bone morphogenetic protein (Bmp) signaling (bmp3/ola-miR-96-5p/miR-181b-5p/miR-199a-5p/miR-205-5p/miR-455-5p). Both pathways are major indicators of impaired bone formation, while the altered Rank signaling in osteoclasts (c-fos/miR-205-5p) suggests a potentially augmented bone resorption. Interestingly, a typical BaP-responsive pathway, the Nrf2-mediated oxidative stress response (gst/ola-miR-181b-5p/miR-199a-5p/miR-205), was also affected. Moreover, mRNA levels of epigenetic modification enzymes (e.g., hdac6, hdac7, kdm5b) were found dysregulated. The findings indicated that epigenetic factors (e.g., miRNAs, histone modifications) may directly regulate the expression of genes associated with transgenerational BaP bone toxicity and warrants further studies. The identified candidate genes and miRNAs may serve as potential biomarkers for BaP-induced bone disease and as indicators of historic exposures in wild fish for conservation purposes.

Hypoxia-induced epigenetic transgenerational miRNAs dysregulation involved in reproductive impairment of ovary

Hypoxia is a potent endocrine disruptor that is posing serious problems to the fish reproductive systems. Our previous studies reported that hypoxia could cause a transgenerational impairment of ovarian development and interfere hatching success in F2 offspring of marine medaka fish (Oryzias melastigma) through epigenetic regulation. As part of the epigenetic regulation, we investigated the involvement of microRNAs (miRNAs) in hypoxia-induced transgenerational reproductive impairments. In the present study, we used comparative small RNA sequencing to reveal that hypoxia caused miRNA dysregulation in ovaries of F0 hypoxia group and F2 transgenerational group. We found 4 common dysregulated miRNA in the F0 and F2 generations. Furthermore, integrated miRNA-mRNA analysis, followed by gene ontology enrichment analysis on the hypoxia-dysregulated miRNA-target genes further highlighted the importance of these dysregulated miRNAs in biological processes related to reproduction. More importantly, we identified 3 miRNA-mRNA pairs (novel miRNA-525-DIAPH2, novel miRNA-525-MYOCD, and novel miRNA-525-RAI14) that might play epigenetic roles in hypoxia-induced reproductive impairment. For the first time, our findings suggested the involvement of miRNA in hypoxia-induced reproductive impairments may be inherited via a transgenerational manner.

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.

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.

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.

MicroRNA expression profile analysis of skin immune response in crucian carp (Carassius auratus) infected by Aeromonas hydrophila

MicroRNAs (miRNAs) are non-coding RNA molecules that regulate gene expression in fish, but its regulatory mechanism of the skin mucosal immune response remains poorly understood. In order to investigate the immunological role of miRNAs, three sRNA libraries (mSC, mST1, mST2) from skin samples of crucian carp (Carassiusauratus) infected with Aeromonas hydrophila at three time points (0, 6 and 12 hpi) were constructed and examined using Illumina Hiseq 2000 platform. All of the identified miRNA, rRNA and tRNA were 69444 (13.39%), 29550 (5.70%) and 10704 (2.06%) in skin, respectively. At 6 and 12 hpi, 829 and 856 miRNAs were differentially expressed, respectively. Among these DEMs, 53 known and 10 novel miRNAs were all significantly differentially expressed during early infection (p < 0.01). GO and KEGG enrichment analyses revealed that 118111 target-genes were primarily involved in cellular process, metabolic process, biological regulation and stress response, such as antigen processing and presentation, complement and coagulation cascades, phagosome, MAPK, TLR, NF-κB and JAK-STAT signaling pathways. These results will help to elucidate the mechanism of miRNAs involved in the skin mucosal immune response of crucian carp against Aeromonas hydrophila infection.

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.

Effects of fasting and refeeding on protein and glucose metabolism in Arctic charr

Refeeding, following a period of food deprivation will often lead to compensatory growth. Although many studies have focused on molecular mechanisms behind this accelerated growth response in fish, little is known on the roles of protein and metabolism. We also assessed, for the first time, the potential roles of miRNAs in regulating compensatory growth. Artcic charr, Salvelinus alpinus, a northern freshwater species, was deprived of food for 101 days and then fed to satiety for 126 days. The refeeding period resulted in compensatory growth, with a partial compensation of body mass. The feed deprivation period lead to a decrease in hepatosomatic index (HSI) and intestinal somatic index (ISI). HSI and ISI were then gradually replenished during early refeeding, following a lag phase prior to the compensatory growth response. mRNA transcripts regulating protein degradation via the autophagy pathway (Cathepsin D and Cathepsin L) in muscle were upregulated during feed restriction and downregulated after refeeding, which could allow for greater protein accretion in muscle, facilitating compensatory growth. Transcript levels from the ubiquitin proteasome pathway (Mafbx and Murf1) and the calpain system (Calpain 7 and Calpastatin) suggested that these pathways were not involved in regulating compensatory growth. Furthermore, we've shown that miRNAs (miR-29a and miR-223) could be involved in fish glycogen homeostasis during the early stages of refeeding. These findings provide a deeper understanding of the molecular mechanisms regulating growth in fish.

The destiny of the resistance/susceptibility against GCRV is controlled by epigenetic mechanisms in CIK cells

Hemorrhagic disease caused by grass carp reovirus (GCRV) has severely threatened the grass carp (Ctenopharyngodon idella) cultivation industry. It is noteworthy that the resistance against GCRV infection was reported to be inheritable, and identified at both individual and cellular levels. Therefore, this work was inspired and dedicated to unravel the molecular mechanisms of fate decision post GCRV infection in related immune cells. Foremost, the resistant and susceptible CIK (C. idella kidney) monoclonal cells were established by single cell sorting, subculturing and infection screening successively. RNA-Seq, MeDIP-Seq and small RNA-Seq were carried out with C1 (CIK cells), R2 (resistant cells) and S3 (susceptible cells) groups. It was demonstrated that genome-wide DNA methylation, mRNA and microRNA expression levels in S3 were the highest among three groups. Transcriptome analysis elucidated that pathways associated with antioxidant activity, cell proliferation regulation, apoptosis activity and energy consuming might contribute to the decision of cell fates post infection. And a series of immune-related genes were identified differentially expressed across resistant and susceptible groups, which were negatively modulated by DNA methylation or microRNAs. To conclude, this study systematically uncovered the regulatory mechanism on the resistance from epigenetic perspective and provided potential biomarkers for future studies on resistance breeding.