MiR-30c: A novel regulator of salt tolerance in tilapia

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

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

Identification and functional characterization of a novel TRPA1 gene from sea cucumber Apostichopus japonicus and interaction with miR-2013 in response to salt stress

Salinity is important abiotic factor influencing sea cucumber aquaculture. This study aimed to identify and functional study of a novel transient receptor potential cation channel subfamily A member 1 (TRPA1) involved in salinity stress through interaction with miR-2013 in the sea cucumber. The full-length cDNA sequence was 1369 bp in length and encoded 138 amino acids. The TRPA1 homolog protein was a hydrophilic protein without a signal peptide and was predicted to a spatial structure of seven helices and eight random coils and two major ANK functional domains. Bioinformatic analysis and luciferase reporter assays confirmed TRPA1 as a target gene of miR-2013. Quantitative PCR revealed that miR-2013 was induced upregulation after salinity stress, while TRPA1 showed upregulated expression with maximum expression at 24 h. The expression of miR-2013 and TRPA1 was negatively regulated. Transfection experiments were conducted to validate the role of miR-2013 and TRPA1 in salinity response. The results showed that miR-2013 was upregulated and TRPA1 was downregulated after transfection with miR-2013 mimics, while miR-2013 was downregulated and TRPA1 was upregulated after transfection with miR-2013 inhibitor. Transfection with si-TRPA1 homolog resulted in upregulation of miR-2013 and downregulation of TRPA1 homolog. These findings suggest that miR-2013 can regulate the expression of TRPA1 under salt stress, and highlight the importance of miR-2013 and TRPA1 in salt stress response. miR-2013 mimics improved the survival rate, while miR-2013 inhibitor and si-TRPA1 reduced it. These findings suggest that miR-2013 and TRPA1 play important roles in sea cucumbers adaptation to salinity changes.

Integrative miRNA-mRNA network analysis to identify crucial pathways of salinity adaptation in brain transcriptome of Labeo rohita

Introduction: Brain being the master regulator of the physiology of animal, the current study focuses on the gene expression pattern of the brain tissue with special emphasis on regulation of growth, developmental process of an organism and cellular adaptation of Labeo rohita against unfavourable environmental conditions. Methods: RNA-seq study was performed on collected brain samples at 8ppt salt concentration and analyzed for differential gene expression, functional annotation and miRNA-mRNA regulatory network. Results: We found that 2450 genes were having significant differential up and down regulation. The study identified 20 hub genes based on maximal clique centrality algorithm. These hub genes were mainly involved in various signaling pathways, energy metabolism and ion transportation. Further, 326 up and 1214 down regulated genes were found to be targeted by 7 differentially expressed miRNAs i.e., oni-miR-10712, oni-miR-10736, ssa-miR-221-3p, ssa-miR-130d-1-5p, ssa-miR-144-5p and oni-miR-10628. Gene ontology analysis of these differentially expressed genes led to the finding that these genes were involved in signal transduction i.e., calcium, FOXO, PI3K-AKT, TGF-β, Wnt and p53 signalling pathways. Differentially expressed genes were also involved in regulation of immune response, environmental adaptation i.e., neuroactive ligand-receptor interaction, ECM-receptor interaction, cell adhesion molecules and circadian entrainment, osmoregulation and energy metabolism, which are critical for salinity adaptation. Discussion: The findings of whole transcriptomic study on brain deciphered the miRNA-mRNA interaction patterns and pathways associated with salinity adaptation of L. rohita.

Profiling miRNAs of Teleost Fish in Responses to Environmental Stress: A Review

miRNAs are a class of endogenous and evolutionarily conserved noncoding short RNA molecules that post-transcriptionally regulate gene expression through sequence-specific interactions with mRNAs and are capable of controlling gene expression by binding to miRNA targets and interfering with the final protein output. The miRNAs of teleost were firstly reported in zebrafish development, but there are recent studies on the characteristics and functions of miRNAs in fish, especially when compared with mammals. Environmental factors including salinity, oxygen concentration, temperature, feed, pH, environmental chemicals and seawater metal elements may affect the transcriptional and posttranscriptional regulators of miRNAs, contributing to nearly all biological processes. The survival of aquatic fish is constantly challenged by the changes in these environmental factors. Environmental factors can influence miRNA expression, the functions of miRNAs and their target mRNAs. Progress of available information is reported on the environmental effects of the identified miRNAs, miRNA targets and the use of miRNAs in fish.

The Effect of Salinity Stress on Enzyme Activities, Histology, and Transcriptome of Silver Carp (Hypophthalmichthys molitrix)

A 56-day study was performed to examine the effect of freshwater (FW) and brackish water (BW 6‱ salinity) on the antioxidant ability, Na+/K+-ATPase (NKA) activities, histology, and transcriptome of the gill and kidney tissue in juvenile silver carp (Hypophthalmichthys molitrix). The results show that when juvenile silver carp were exposed to 6‱ salinity, the activities of superoxide dismutase (SOD) and catalase (CAT) were shown to be substantially increased (p < 0.05), while glutathione peroxidase (GSH-PX) activities in gill were not significantly affected (p < 0.05). In kidney tissue, SOD, CAT, and GSH-PX, enzyme activities peaked at 24, 8, and 4 h, respectively, but were not significantly different compared with the control group (p < 0.05). In addition, significant effects of salinity were observed for the NKA level in both the gills and kidney tissues (p < 0.05). The gill filaments of juvenile silver carp under the BW group all underwent adverse changes within 72 h, such as cracks and ruptures in the main part of the gill filaments, bending of the gill lamellae and enlargement of the gaps, and an increase in the number of mucus and chloride-secreting cells. Transcriptome sequencing showed 171 and 261 genes in the gill and kidney tissues of juvenile silver carp compared to the BW group, respectively. Based on their gene ontology annotations, transcripts were sorted into four functional gene groups, each of which may play a role in salt tolerance. Systems involved in these processes include metabolism, signal transduction, immunoinflammatory response, and ion transport. The above findings indicate that the regulation processes in juvenile silver carp under brackish water conditions are complex and multifaceted. These processes and mechanisms shed light on the regulatory mechanism of silver carp osmolarity and provide a theoretical foundation for future research into silver carp growth in brackish water aquaculture area.

Modulation of miR-429 during osmotic stress in the silverside Odontesthes humensis

Silverside fish inhabit marine coastal waters, coastal lagoons, and estuarine regions in southern South America. Although silversides are not fully adapted to freshwater, they can tolerate a wide range of salinity variations. MicroRNAs (miRNAs) are a class of ∼22 nucleotide noncoding RNAs, which are crucial regulators of gene expression at post-transcriptional level. Current data indicate that miRNAs biogenesis is altered by situations of environmental stress, thereby altering the expression of target mRNAs. Foremost, the silversides were acutely exposed to 30 g.L−1 of salt to reveal in which tissue miR-429 could be differentially expressed. Thus, fish were acclimated to freshwater (0 g.L−1) and to brackish water (10 g.L−1), and then exposed to opposite salinity treatment. Here, we reveal that miR-429, a gill-enriched miRNA, emerges as a prime osmoregulator in silversides. Taken together, our findings suggest that miR-429 is an endogenous regulator of osmotic stress, which may be developed as a biomarker to assist silverside aquaculture.

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.

Variant expression signatures of microRNAs and protein related to growth in a crossbreed between two strains of Nile tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is a species of worldwide importance for aquaculture. A crossbred lineage was developed through introgressive backcross breeding techniques and combines the high growth performance of the Chitralada (CHIT) lwith attractive reddish color of the Red Stirling (REDS) strains. Since the crossbreed has an unknown genetically improved background, the objective of this work was to characterize expression signatures that portray the advantageous phenotype of the crossbreeds. We characterized the microRNA transcriptome by high throughput sequencing (RNA-seq) and the proteome through mass spectrometry (ESI-Q-TOF-MS) and applied bioinformatics for the comparative analysis of such molecular data on the three strains. Crossbreed expressed a distinct set of miRNAs and proteins compared to the parents. They comprised several microRNAs regulate traits of economic interest. Proteomic profiles revealed differences between parental and crossbreed in expression of proteins associated with glycolisis. Distinctive miRNA and protein signatures contribute to the phenotype of crossbreed.

Novel lncRNA lncRNA001074 participates in the low salinity-induced response in the sea cucumber Apostichopus japonicus by targeting the let-7/NKAα axis

Salinity fluctuations have severe impacts on sea cucumbers and therefore important consequences in sea cucumber farming. The responses of sea cucumbers to salinity changes are reflected in the expression profiles of multiple genes and non-coding RNAs (ncRNAs). The microRNA (let-7) which is a developmental regulator, the ion transporter gene sodium potassium ATPase gene (NKAα), and the long ncRNA lncRNA001074 were previously shown to be involved in responses to salinity changes in various marine species. To better understand the relationship between ncRNAs and target genes, the let-7/NKAα/lncRNA001074 predicted interaction was investigated in this study using luciferase reporter assays and gene knockdowns in the sea cucumber Apostichopus japonicus. The results showed that NKAα was the target gene of let-7 and NKAα expression levels were inversely correlated with let-7 expression based on the luciferase reporter assays and western blots. The let-7 abundance was negatively regulated by lncRNA001074 and NKAα both in vitro and in vivo. Knockdown of lncRNA001074 led to let-7 overexpression. These results demonstrated that lncRNA001074 binds to the 3'-UTR binding site of let-7 in a regulatory manner. Furthermore, the expression profiles of let-7, NKAα, and lncRNA001074 were analyzed in sea cucumbers after the knockdown of each of these genes. The results found that lncRNA001074 competitively bound let-7 to suppress NKAα expression under low salinity conditions. The downregulation of let-7, in conjunction with the upregulation of lncRNA001074 and NKAα, may be essential for the response to low salinity change in sea cucumbers. Therefore, the dynamic balance of the lncRNA001074, NKAα, and let-7 network might be a potential response mechanism to salinity change in sea cucumbers.

Identification and Characterization of Osmoregulation Related MicroRNAs in Gills of Hybrid Tilapia Under Three Types of Osmotic Stress

Researchers have increasingly suggested that microRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression and protein translation in organs and respond to abiotic and biotic stressors. To understand the function of miRNAs in osmotic stress regulation of the gills of hybrid tilapia (Oreochromis mossambicus ♀ × Oreochromis urolepis hornorum ♂), high-throughput Illumina deep sequencing technology was used to investigate the expression profiles of miRNAs under salinity stress (S, 25‰), alkalinity stress (A, 4‰) and salinity-alkalinity stress (SA, S: 15‰, A: 4‰) challenges. The results showed that 31, 41, and 27 upregulated and 33, 42, and 40 downregulated miRNAs (P < 0.05) were identified in the salt stress, alkali stress, and saline-alkali stress group, respectively, which were compared with those in the control group (C). Fourteen significantly differently expressed miRNAs were selected randomly and then validated by a quantitative polymerase chain reaction. On the basis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, genes related to osmoregulation and biosynthesis were enriched in the three types of osmotic stress. In addition, three miRNAs and three predicted target genes were chosen to conduct a quantitative polymerase chain reaction in the hybrid tilapia and its parents during 96-h osmotic stress. Differential expression patterns of miRNAs provided the basis for research data to further investigate the miRNA-modulating networks in osmoregulation of teleost.

Fluoride exposure changed the expression of microRNAs in gills of male zebrafish (Danio rerio)

Fluoride has been found to cause detrimental effects on fish gills. Despite essential roles in various metabolism activities, whether and how miRNAs participate in the course of toxicity caused by fluoride in gills is still unclear. In this study, male zebrafish were exposed to 0, 20, 40 mg/L fluoride for 60 days to study the underlying osmotic regulatory mechanism by determining the influences of fluoride on the miRNAs and regulated genes in gills. mRNAs were isolated from the gills and the expression profiles were analyzed by using Illumina Hiseq 2500 platforms. Expressions of 7 differentially miRNAs and some related-genes in gills were validated by qRT-PCR. The results showed that miRNAs expressions were notably altered by fluoride. A total of 584 and 327 miRNAs were remarkably changed after 20 and 40 mg/L fluoride exposure, of which 322 were increased and 262 were decreased in 20 mg/L fluoride group, whereas 219 were elevated and 108 were reduced in 40 mg/L fluoride group. The differentially expressive miRNAs confirmed by qRT-PCR were consistent with micro-assay data. Cluster of Orthologous Groups of proteins (COG) function classification showed that the target genes of differentially expressive miRNAs are mainly related to signal transduction mechanisms, replication, transcription, inorganic ion transport and metabolism, repair and recombination, and energy formation and transformation. In addition, fluoride disturbed the expressions of target genes involved in the osmoregulation of the gill in the fluoride-exposed zebrafish, such as the increased expressions of OSTF1 and the decreased expressions of Na⁺-K⁺-ATPase, CFTR, and AQP-3, which provides a possibility that miRNAs regulation induced by fluoride has an effects on osmotic regulation, providing new hints to the osmotic regulatory mechanism of the toxicity caused by fluoride in zebrafish, and distinguishes new biomarkers of miRNAs for fluoride toxicity.

A novel tonicity-responsive microRNA miR-23a-5p modulates renal cell survival under osmotic stress through targeting heat shock protein 70 HSPA1B

There is growing evidence that microRNAs (miRNAs) are implicated in cellular adaptation to osmotic stress, but the underlying osmosignaling pathways are still not completely understood. In this study, we found that a passenger strand miRNA, miR-23a-5p, was significantly downregulated in response to high NaCl treatment in mouse inner medullary collecting duct cells (mIMCD3) through an miRNA profiling assay. The decrease of miR-23a-5p is hypertonicity-dependent and osmotolerant cell type-specific. Knockdown of miR-23a-5p increased cellular survival and proliferation in mIMCD3. In contrast, miR-23a-5p overexpression repressed cell viability and proliferation under hypertonic stress. RNA deep-sequencing revealed that a heat shock protein 70 (HSP70) isoform, HSP70 member 1B (HSPA1B), was significantly increased under hypertonic treatment. Based on the prediction analysis by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and TargetScan, and a further validation via a dual-luciferase assay, HSPA1B was identified as a potential target of miR-23a-5p. Overexpressed miR-23a-5p suppressed HSPA1B, whereas downregulated miR-23a-5p promoted HSPA1B expression in mIMCD3. In addition, an in vivo study demonstrated that there is a reverse correlation between the levels of miR-23a-5p and HSPA1B in mouse renal inner medulla (papilla) that is exposed to extremely high osmolality. In summary, this study elucidates that passenger strand miR-23a-5p is a novel tonicity-responsive miRNA. The downregulation of miR-23a-5p facilitates cellular adaptation to hypertonic stress in mammalian renal cells through modulating HSPA1B.

Thermal stress accelerates mercury chloride toxicity in Oreochromis niloticus via up-regulation of mercury bioaccumulation and HSP70 mRNA expression

Mercury (Hg) is an environmental pollutant that threatens aquatic life. Many environmental factors, including water temperature, are reported to influence the toxicity of dissolved chemicals in the aquatic ecosystem. Therefore, we investigated the impact of thermal stress on Hg-induced subchronic toxicity in Nile tilapia (Oreochromis niloticus). Fish were randomly allocated into five groups. Group I served as the control and kept at 25 °C. Groups II, III, IV, and V were reared at 25, 28, 31, and 34 °C, respectively, and co-exposed to HgCl₂ (1/10 LC₅₀) for 42 days. Blood and tissue samples were collected after 21 and 42 days. All HgCl₂-exposed groups exhibited significant elevations in serum levels of alanine aminotransferase (ALT), aspartate aminotransferase (AST), urea, and creatinine, along with decreases in the serum total protein and albumin. In addition, marked reductions in antioxidant enzymes, including superoxide dismutase (SOD), catalase (CAT), and glutathione peroxidase (GSPx), were observed. Remarkable increases in Hg tissue concentrations were detected along with increases in heat shock protein (HSP) 70 mRNA expression. Interestingly, the patterns data that were recorded were more coincident with the water temperature than the period of exposure. In conclusion, water temperature and exposure period are two crucial factors modulating HgCl₂-induced toxicity and bioaccumulation in Nile tilapia. Our findings provide new insights concerning the impact of thermal stress as an environmental factor on Hg toxicity and bioaccumulation in Nile tilapia and, in turn, on fish and fish consumer health.

Genome-wide analysis of MicroRNA-messenger RNA interactome in ex-vivo gill filaments, Anguilla japonica

Background

Gills of euryhaline fishes possess great physiological and structural plasticity to adapt to large changes in external osmolality and to participate in ion uptake/excretion, which is essential for the re-establishment of fluid and electrolyte homeostasis. The osmoregulatory plasticity of gills provides an excellent model to study the role of microRNAs (miRs) in adaptive osmotic responses. The present study is to characterize an ex-vivo gill filament culture and using omics approach, to decipher the interaction between tonicity-responsive miRs and gene targets, in orchestrating the osmotic stress-induced responses.

Results

Ex-vivo gill filament culture was exposed to Leibovitz’s L-15 medium (300 mOsmol l^− 1) or the medium with an adjusted osmolality of 600 mOsmol l^− 1 for 4, 8 and 24 h. Hypertonic responsive genes, including osmotic stress transcriptional factor, Na⁺/Cl⁻-taurine transporter, Na⁺/H⁺ exchange regulatory cofactor, cystic fibrosis transmembrane regulator, inward rectifying K⁺ channel, Na⁺/K⁺-ATPase, and calcium-transporting ATPase were significantly upregulated, while the hypo-osmotic gene, V-type proton ATPase was downregulated. The data illustrated that the ex-vivo gill filament culture exhibited distinctive responses to hyperosmotic challenge. In the hyperosmotic treatment, four key factors (i.e. drosha RNase III endonuclease, exportin-5, dicer ribonuclease III and argonaute-2) involved in miR biogenesis were dysregulated (P < 0.05). Transcriptome and miR-sequencing of gill filament samples at 4 and 8 h were conducted and two downregulated miRs, miR-29b-3p and miR-200b-3p were identified. An inhibition of miR-29b-3p and miR-200b-3p in primary gill cell culture led to an upregulation of 100 and 93 gene transcripts, respectively. Commonly upregulated gene transcripts from the hyperosmotic experiments and miR-inhibition studies, were overlaid, in which two miR-29b-3p target-genes [Krueppel-like factor 4 (klf4), Homeobox protein Meis2] and one miR-200b-3p target-gene (slc17a5) were identified. Integrated miR-mRNA-omics analysis revealed the specific binding of miR-29b-3p on Klf4 and miR-200b-3p on slc17a5. The target-genes are known to regulate differentiation of gill ionocytes and cellular osmolality.

Conclusions

In this study, we have characterized the hypo-osmoregulatory responses and unraveled the modulation of miR-biogenesis factors/the dysregulation of miRs, using ex-vivo gill filament culture. MicroRNA-messenger RNA interactome analysis of miR-29b-3p and miR-200b-3p revealed the gene targets are essential for osmotic stress responses.

Transcriptome changes for Nile tilapia (Oreochromis niloticus) in response to alkalinity stress

Nile tilapia is an important economic fish in the world because of its fast growth, high meat yield and strong adaptability. It is more adaptable to high alkalinity than common freshwater fish and provides valuable material for developing alkaline-tolerant strains and understanding the adaptation mechanism of fish to extreme environmental stress. In this study, we employed high throughput RNA sequencing to reveal the tissues (gill, kidney and liver) transcriptome differences of O. niloticus at different carbonate alkalinities (FW, AW40 and AW60). A total of 1,369,381,790 raw reads were obtained, including 496,441,232 reads in FW group, 437,907,696 reads in AW40 and 435,032,862 reads in AW60. In addition, 484,555,626 reads in gill, 451,618,224 reads in kidney and 433,207,940 reads in liver. A large number of stress-regulated changes were detected comprehensively. We focused on 3 significantly change pathways (steroid biosynthesis, drug metabolism and protein digestion/absorption) and 17 DEGs (HMG-CoA reductase, UDP-glucuronosyltransferase, and carbonic anhydrase etc.) which were shared among compared groups (AW40 vs FW, AW60 vs FW, AW40 vs 60 AW60) in gill, kidney and liver, respectively. These pathways/genes are sensitive to alkalinity stress and crucial to the alkalinity adaptation of tilapia. Overall, we found a large number of candidate genes, which encode important regulators of stress tolerance and ultimately contribute to future alkaline-tolerant fish breeding. Among these genes, lipid metabolism (involving signal transduction), detoxification and immune related genes are more prominent to the response and adaptability of fish to alkalinity stress.

miR-10 involved in salinity-induced stress responses and targets TBC1D5 in the sea cucumber, Apostichopus japonicas

The sea cucumber is an economically important aquaculture species in China, where it encounter hypo-saline conditions caused by freshwater outflow from rivers and rainfall. MicroRNAs (miRNA) are small noncoding RNAs of about 22 nucleotides, which are crucial regulators of gene expression at the post-transcriptional level and are involved in a variety of physiological and pathophysiological processes. miR-10 is differentially expressed in salinity acclimation, and has a seed-region match with TBC1D5. The expression profiles of miR-10 and TBC1D5 indicate that miR-10 negatively regulates the expression of TBC1D5 in coelomocytes and sea cucumbers with a miR-10 agomir or antagomir. During salinity acclimation, up-regulation of miR-10 was induced after transfection in coelomocytes with a miR-10 inhibitor, while down-regulation of TBC1D5 was induced. The miR-10 expression maximum in coelomocytes appeared at 48 h post-transfection with a miR-10 inhibitor, was later than that of in sea cucumbers, which appeared 24 h after miR-10 antagomir injection. There was no longer a negative relationship between miR-10 and TBC1D5 expression in coelomocytes and sea cucumbers with miR-10 mimics or agomir during salinity acclimation. The miR-10 antagomir or agomir only affected sodium and NKA enzyme activities, and has little effect on other chloride and potassium ions. Our results demonstrate miR-10 directly regulates TBC1D5 by targeting its 3'-UTR, and that miR-10 suppression substantially increases TBC1D5 mRNA levels in vivo and in vitro. Furthermore, miR-10 and TBC1D5 fluctuating expression patterns after treatment with a miR-10 inhibitor or mimics during salinity acclimation may indicate that they are required for adaptation to salinity stress caused by environmental change. Especially, the miR-10 up-regulation in coelomocytes with miR-10 inhibitor during salinity acclimation indicated that they are required for adaptation to salinity stress caused by environmental change. We propose that miR-10 participates in a regulatory circuit that allows for rapid gene program transitions in response to osmotic stress.

Generation of microRNA-30e-producing recombinant viral hemorrhagic septicemia virus (VHSV) and its effect on in vitro immune responses

MicroRNAs (miRNAs) are non-coding small RNAs involved in the regulation of gene expression. In the present study, we firstly reported the use of a fish RNA virus, viral hemorrhagic septicemia virus (VHSV), as a delivery vehicle of a miRNA-30e, and the effect of miR-30e produced by the recombinant VHSV on the immune responses of Epithelioma papulosum cyprini (EPC) cells was investigated. The expression of functional miR-30e using a CMV promoter-driven vector was verified by the significantly lower eGFP expression in cells transfected with a vector containing miR-30e sponge sequence than that in cells transfected with a control vector that had mutated miR-30e sponge sequence. Furthermore, the down-regulation of reporter gene containing 3'-UTR of NF-κb inhibitor α-like protein B (NFκbiαb) by miR-30e was demonstrated, suggesting that miR-30e overexpression can increase immune responses related to NF-κB activation through inhibition of IκB. A miR-30e-expressing recombinant VHSV (rVHSV-A-miR30e) that had primary microRNA-30e sequence between N and P genes was rescued using the reverse genetic method, and the successful expression of miR-30e in the cells infected with rVHSV-A-miR30e was demonstrated using Northern blot and qRT-PCR. Cells infected with rVHSV-A-miR30e showed the increase of NF-κB activation and type I interferon induced genes expression, suggesting that rVHSV-A-miR30e can produce functional miR-30e in fish cells, and VHSV can be used as a vehicle to deliver functional microRNAs in fish.

Salinity stress-induced differentially expressed miRNAs and target genes in sea cucumbers Apostichopus japonicus

Environmental salinity is an important abiotic factor influencing normal physiological functions and productive performance in the sea cucumber Apostichopus japonicus. It is therefore important to understand how changes in salinity affect sea cucumbers in the face of global climate change. In this study, we investigated the responses to salinity stress in sea cucumbers using mRNA and miRNA sequencing. The regulatory network of mRNAs and miRNAs involved in salinity stress was examined, and the metabolic pathways enriched for differentially expressed miRNAs and target mRNAs were identified. The top 20 pathways were involved in carbohydrate metabolism, fatty acid metabolism, degradation, and elongation, amino acid metabolism, genetic information processing, metabolism of cofactors and vitamins, transport and catabolism, and environmental information processing. A total of 22 miRNAs showed differential expression during salinity acclimation. The predicted 134 target genes were enriched in functions consistent with the results of gene enrichment based on transcriptome analysis. These results suggested that sea cucumbers deal with salinity stress via changes in amino acid metabolism, ion channels, transporters, and aquaporins, under stimulation by environmental signals, and that this process requires energy from carbohydrate and fatty acid metabolism. Salinity challenge also induced miRNA expression. These results provide a valuable genomic resource that extends our understanding of the unique biological characteristics of this economically important species under conditions of salinity stress.

Identifying a Long QTL Cluster Across chrLG18 Associated with Salt Tolerance in Tilapia Using GWAS and QTL-seq

Understanding the genetic mechanism of osmoregulation is important for the improvement of salt tolerance in tilapia. In our previous study, we have identified a major quantitative trait locus (QTL) region located at 23.0 Mb of chrLG18 in a Nile tilapia line by QTL-seq. However, the conservation of these QTLs in other tilapia populations or species is not clear. In this study, we successfully investigated the QTLs associated with salt tolerance in a mass cross population from the GIFT line of Nile tilapia (Oreochromis niloticus) using a ddRAD-seq-based genome-wide association study (GWAS) and in a full-sib family from the Malaysia red tilapia strain (Oreochromis spp) using QTL-seq. Our study confirmed the major QTL interval that is located at nearly 23.0 Mb of chrLG18 in Nile tilapia and revealed a long QTL cluster across chrLG18 controlling for the salt-tolerant trait in both red tilapia and Nile tilapia. This is the first GWAS analysis on salt tolerance in tilapia. Our finding provides important insights into the genetic architecture of salinity tolerance in tilapia and supplies a basis for fine mapping QTLs, marker-assisted selection, and further detailed functional analysis of the underlying genes for salt tolerance in tilapia.

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

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

Genome-wide identification and differentially expression analysis of lncRNAs in tilapia

Background

Long noncoding RNAs (LncRNAs) play important roles in fundamental biological processes. However, knowledge about the genome-wide distribution and stress-related expression of lncRNAs in tilapia is still limited.

Results

Genome-wide identification of lncRNAs in the tilapia genome was carried out in this study using bioinformatics tools. 103 RNAseq datasets that generated in our laboratory or collected from NCBI database were analyzed. In total, 72,276 high-confidence lncRNAs were identified. The averaged positive correlation coefficient (r_mean = 0.286) between overlapped lncRNA and mRNA pairs showed significant differences with the values for all lncRNA-mRNA pairs (r_mean = 0.176, z statistics = − 2.45, p value = 0.00071) and mRNA-mRNA pairs (r_mean = 0.186, z statistics = − 2.23, p value = 0.0129). Weighted correlation network analysis of the lncRNA and mRNA datasets from 12 tissues identified 21 modules and many interesting mRNA genes that clustered with lncRNAs. Overrepresentation test indicated that these mRNAs enriched in many biological processes, such as meiosis (p = 0.00164), DNA replication (p = 0.00246), metabolic process (p = 0.000838) and in molecular function, e.g., helicase activity (p = 0.000102) and catalytic activity (p = 0.0000612). Differential expression (DE) analysis identified 99 stress-related lncRNA genes and 1955 tissue-specific DE lncRNA genes. MiRNA-lncRNA interaction analysis detected 72,267 lncRNAs containing motifs with sequence complementary to 458 miRNAs.

Conclusions

This study provides an invaluable resource for further studies on molecular bases of lncRNAs in tilapia genomes. Further function analysis of the lncRNAs will help to elucidate their roles in regulating stress-related adaptation in tilapia.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5115-x) contains supplementary material, which is available to authorized users.

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.

Identifying a Major QTL Associated with Salinity Tolerance in Nile Tilapia Using QTL-Seq

Selection of new lines with high salinity tolerance allows for economically feasible production of tilapias in brackish water areas. Mapping QTLs and identifying the markers linked to salinity-tolerant traits are the first steps in the improvement of the tolerance in tilapia through marker-assisted selection techniques. By using QTL-seq strategy and linkage-based analysis, two significant QTL intervals (chrLG4 and chrLG18) on salinity-tolerant traits were firstly identified in the Nile tilapia. Fine mapping with microsatellite and SNP markers suggested a major QTL region that located at 23.0 Mb of chrLG18 and explained 79% of phenotypic variation with a LOD value of 95. Expression analysis indicated that at least 10 genes (e.g., LACTB2, KINH, NCOA2, DIP2C, LARP4B, PEX5R, and KCNJ9) near or within the QTL interval were significantly differentially expressed in intestines, brains, or gills under 10, 15, or 20 ppt challenges. Our findings suggest that QTL-seq can be effectively utilized in QTL mapping of salinity-tolerant traits in fish. The identified major QTL is a promising locus to improve our knowledge on the genetic mechanism of salinity tolerance in tilapia.

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.

Differential miRNA expression in the three-spined stickleback, response to environmental changes

miRNAs play important role in the various physiological and evolutionary processes, however, there is no data allowing comparison of evolutionary differences between various ecotypes adapted to different environmental conditions and specimen demonstrating immediate physiological response to the environmental changes. We compared miRNA expression profiles between marine and freshwater stickleback populations of the three-spined stickleback to identify the evolutionary differences. To study the immediate physiological response to foreign environment, we explored the changes induced by transfer of marine sticklebacks into freshwater environment and vice versa. Comparative analysis of changes in miRNA expression suggested that they are driven by three independent factors: (1) non-specific changes in miRNA expression under different environmental conditions; (2) specific response to freshwater conditions in the marine stickleback ecotype; (3) specific response to extreme osmotic conditions for both marine and freshwater ecotypes during the contact with non-native environment. Gene Ontology enrichment analysis of differential expressed miRNA targets supports our current hypothesis.

Interaction of heat shock protein 60 (HSP60) with microRNA in Chinese mitten crab during Spiroplasma eriocheiris infection

Heat shock protein 60 from the Chinese mitten crab Eriocheir sinensis (EsHSP60) was previously identified in relation to Spiroplasma eriocheiris infection by isobaric tags for relative and absolute quantitation labelling followed by liquid chromatography-tandem mass spectrometry. In the present study, to validate the immune function of this protein, the cDNA of the EsHSP60 gene was cloned. Various crab tissues were assessed using real-time PCR, which showed that EsHSP60 transcription occurred in all tissues examined. The expression profiles of EsHSP60 in haemolymph at transcription and protein levels when infected with S. eriocheiris were investigated by real-time PCR and Western blot analysis, respectively. A significant increase of EsHSP60 transcription and protein expression appeared post-injection in response to S. eriocheiris infection when compared to the control group. The double-luciferase reporter gene assay showed that the microRNA PC-533-3p interacted with the 3'-untranslated region of EsHSP60 and inhibited the translation of EsHSP60. The expression profiles of PC-533-3p during S. eriocheiris infection were also investigated by real-time PCR. However, the change tendency of PC-533-3p was opposite to that of the EsHSP60 after S. eriocheiris challenge. These data indicate that the EsHSP60 proteins may play an important role in mediating the immune responses of E. sinensis to an S. eriocheiris challenge.

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

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

Identification and characterization of microRNAs and their target genes from Nile tilapia (Oreochromis niloticus)

MicroRNAs (miRNAs) are a class of small single-stranded, endogenous 21–22 nt non-coding RNAs that regulate their target mRNA levels by causing either inactivation or degradation of the mRNAs. In recent years, miRNA genes have been identified from mammals, insects, worms, plants, and viruses. In this research, bioinformatics approaches were used to predict potential miRNAs and their targets in Nile tilapia from the expressed sequence tag (EST) and genomic survey sequence (GSS) database, respectively, based on the conservation of miRNAs in many animal species. A total of 19 potential miRNAs were detected following a range of strict filtering criteria. To test the validity of the bioinformatics method, seven predicted Nile tilapia miRNA genes were selected for further biological validation, and their mature miRNA transcripts were successfully detected by stem–loop RT-PCR experiments. Using these potential miRNAs, we found 56 potential targets in this species. Most of the target mRNAs appear to be involved in development, metabolism, signal transduction, transcription regulation and stress responses. Overall, our findings will provide an important foundation for further research on miRNAs function in the Nile tilapia.

Role of miR-21 in alkalinity stress tolerance in tilapia

MicroRNAs (miRNAs) are a class of short, evolutionary conserved non-coding RNA molecules, which are shown as the key regulators of many biological functions. External stress can alter miRNA expression levels, thereby changing the expression of mRNA target genes. Here, we show that miR-21 is involved in the regulation of alkalinity tolerance in Nile tilapia. Alkalinity stress results in a marked reduction in miR-21 levels. miR-21 loss of function could affect ion balance regulation, ROS production, and antioxidant enzyme activity in vivo. Moreover, miR-21 knockdown protects cell against alkalinity stress-induced injury in vitro. miR-21 directly regulates VEGFB and VEGFC expression by targeting the 3'-untranslated regions (UTRs) of their mRNAs, and inhibition of miR-21 significantly increases the levels of VEGFB and VEGFC expression in vivo. Taken together, our study reveals that miR-21 knockdown plays a protective role in alkalinity tolerance in tilapia.

MicroRNA-Sequence Profiling Reveals Novel Osmoregulatory MicroRNA Expression Patterns in Catadromous Eel Anguilla marmorata

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that regulate gene expression by post-transcriptional repression of mRNAs. Recently, several miRNAs have been confirmed to execute directly or indirectly osmoregulatory functions in fish via translational control. In order to clarify whether miRNAs play relevant roles in the osmoregulation of Anguilla marmorata, three sRNA libraries of A. marmorata during adjusting to three various salinities were sequenced by Illumina sRNA deep sequencing methods. Totally 11,339,168, 11,958,406 and 12,568,964 clear reads were obtained from 3 different libraries, respectively. Meanwhile, 34 conserved miRNAs and 613 novel miRNAs were identified using the sequence data. MiR-10b-5p, miR-181a, miR-26a-5p, miR-30d and miR-99a-5p were dominantly expressed in eels at three salinities. Totally 29 mature miRNAs were significantly up-regulated, while 72 mature miRNAs were significantly down-regulated in brackish water (10‰ salinity) compared with fresh water (0‰ salinity); 24 mature miRNAs were significantly up-regulated, while 54 mature miRNAs were significantly down-regulated in sea water (25‰ salinity) compared with fresh water. Similarly, 24 mature miRNAs were significantly up-regulated, while 45 mature miRNAs were significantly down-regulated in sea water compared with brackish water. The expression patterns of 12 dominantly expressed miRNAs were analyzed at different time points when the eels transferred from fresh water to brackish water or to sea water. These miRNAs showed differential expression patterns in eels at distinct salinities. Interestingly, miR-122, miR-140-3p and miR-10b-5p demonstrated osmoregulatory effects in certain salinities. In addition, the identification and characterization of differentially expressed miRNAs at different salinities can clarify the osmoregulatory roles of miRNAs, which will shed lights for future studies on osmoregulation in fish.

MicroRNA in teleost fish

MicroRNAs (miRNAs) are transcriptional and posttranscriptional regulators involved in nearly all known biological processes in distant eukaryotic clades. Their discovery and functional characterization have broadened our understanding of biological regulatory mechanisms in animals and plants. They show both evolutionary conserved and unique features across Metazoa. Here, we present the current status of the knowledge about the role of miRNA in development, growth, and physiology of teleost fishes, in comparison to other vertebrates. Infraclass Teleostei is the most abundant group among vertebrate lineage. Fish are an important component of aquatic ecosystems and human life, being the prolific source of animal proteins worldwide and a vertebrate model for biomedical research. We review miRNA biogenesis, regulation, modifications, and mechanisms of action. Specific sections are devoted to the role of miRNA in teleost development, organogenesis, tissue differentiation, growth, regeneration, reproduction, endocrine system, and responses to environmental stimuli. Each section discusses gaps in the current knowledge and pinpoints the future directions of research on miRNA in teleosts.

MicroRNA repertoire for functional genome research in tilapia identified by deep sequencing

The Nile tilapia (Oreochromis niloticus; Cichlidae) is an economically important species in aquaculture and occupies a prominent position in the aquaculture industry. MicroRNAs (miRNAs) are a class of noncoding RNAs that post-transcriptionally regulate gene expression involved in diverse biological and metabolic processes. To increase the repertoire of miRNAs characterized in tilapia, we used the Illumina/Solexa sequencing technology to sequence a small RNA library using pooled RNA sample isolated from the different developmental stages of tilapia. Bioinformatic analyses suggest that 197 conserved and 27 novel miRNAs are expressed in tilapia. Sequence alignments indicate that all tested miRNAs and miRNAs* are highly conserved across many species. In addition, we characterized the tissue expression patterns of five miRNAs using real-time quantitative PCR. We found that miR-1/206, miR-7/9, and miR-122 is abundantly expressed in muscle, brain, and liver, respectively, implying a potential role in the regulation of tissue differentiation or the maintenance of tissue identity. Overall, our results expand the number of tilapia miRNAs, and the discovery of miRNAs in tilapia genome contributes to a better understanding the role of miRNAs in regulating diverse biological processes.

miR-206 regulates the growth of the teleost tilapia (Oreochromis niloticus) through the modulation of IGF-1 gene expression

MicroRNAs (miRNAs) are ~22-nucleotide noncoding RNAs that play a crucial role in regulating muscle development. Our previous study shows that miR-206 is specifically expressed in tilapia skeletal muscle, and exhibits a dynamic expression pattern at different developmental stages. Here, we reveal that miR-206 emerges as a crucial regulator of tilapia growth. miR-206 loss of function leads to the acceleration of tilapia growth. IGF-1 is identified as the target gene of miR-206. miR-206 directly changes IGF-1 expression by targeting its 3' UTR, and inhibition of miR-206 substantially increases the IGF-1 mRNA level in vivo. Thus, miR-206 could be developed as a molecular marker to assist fish breeding.

Mechanism of osmoregulatory adaptation in tilapia

The shortage of freshwater resource in many countries leads to a shift to develop aquaculture in brackish water and sea water. Tilapias are euryhaline that can thrive from freshwater to full sea water. They and their hybrids are the best candidate species for cultivation in brackish habitats. Thus, understanding their osmoregulatory mechanisms will help to breed or genetically engineer salt tolerant species. In this paper, we review recent progress in understanding the mechanisms of osmoregulatory adaptations in tilapia.