- Research article
- Open Access
Genome-wide analysis of MicroRNA-messenger RNA interactome in ex-vivo gill filaments, Anguilla japonica
BMC Genomics volume 21, Article number: 208 (2020)
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.
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.
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.
The catadromous fish Japanese eels have a complex life cycle in freshwater and seawater environments. The fish is euryhaline and actively engage physiological responses to oppose osmotic perturbations to stabilize body osmolality. Fish gills therefore possess great structural and functional plasticity in response to variations of water osmolality, to support the process of fluid and electrolyte homeostasis [1, 2]. In response to osmotic stress, a distinct suite of modifications in gill epithelia is activated for functional adaptation, which involves cell proliferation and differentiation, changes in the activities, expressions and trafficking of different ion transporters/channels [3, 4]. In the past decades, considerable numbers of reports identified the underlying molecular and physiological factors associated with ion-osmoregulatory functions of gills. In recent, advances in the understanding of small non-coding RNAs such as microRNAs (miRs) on regulatory circuit in physiological and pathophysiological functions shed light on roles of miRs in osmotic stress responses.
In mammals, the involvement of miRs in fluid and electrolyte transport has become apparent. A number of studies have illustrated the implications of miRs in ion transport functions. For instances, a direct regulatory role of miR-9 in the expression of large-conductance calcium and voltage-activated potassium channels in neuronal ion transport in alcohol tolerance phenomenon in an adult rat model was identified . The involvement of miR-133a and miR-133b in primary culture of canine cardiomyocytes  and miR-1  on potassium channels for cardiac arrhythmias in human samples and rat models were demonstrated. Moreover, indirect regulatory roles of miRs (i.e. miR-155, miR124 & miR-135a) on ion transport were reported via their inhibitory actions on receptors of electrolyte-regulated hormones, angiotensin-I receptor or mineralocorticoid receptor, in screening samples of patients with hypertension [8, 9]. The expression level of the with-no-lysine kinase (wnk1, a regulator of electrolyte homeostasis) in mouse nephrons, was found to be regulated by miR-192. Its expression was modulated by physiological stimuli (i.e. aldosterone or salts loading) . In renal medullary epithelial mIMCD3 cells, tonicity was also identified as a stimulus to regulate miR-200b and miR-717 expression, which inhibit the expression of the transcriptional factor, osmotic response element binding protein (orebp) . In comparison, studies of miRs in osmotic stress responses in fish are limited. Using zebrafish embryos, the expression of the miR-8 family (miR-200) in ionocytes was found to inhibit the expression of Na+/K+-exchanger regulatory cofactor (nherf1) to impair cellular responses to osmotic stress . Using tilapia, decreased expression levels of miR-429 in gills under hyperosmotic stress was observed. An in-vivo inhibitory effect of miR-429 on the expression of osmoregulatory transcription factor (ostf1) in gills  and an inhibitory action of miR-30c on the expression of renal hsp70 under hyperosmotic stress  were reported. Moreover, differential expression patterns of miRs in gills of marbled eels, Anguilla marmorata adapted at different salinities [i.e. freshwater (FW), brackish-water or seawater (SW)] were described . Japanese eels are euryhaline fish. The osmoregulatory tissue - gills provides an excellent model to study role of miRs in the regulation on plasticity of adaptive osmotic responses in vertebrates. The present study aimed to identify and characterize the involvement of miRs and messenger RNAs (mRNA) under osmotic perturbations. In this report, we used both ex-vivo gill filament and primary gill cell culture models, accompanied with miR-, transcriptome-sequencing and miR inhibition, to identify tonicity-sensitive miRs and to characterize their expressions in gill filaments.
Hyperosmotic treatment induced differential expression of transcriptional Factor-1/regulators, ion channels/transporters and miR biogenesis factors in ex-vivo gill filament culture
Gill filaments of eels were ex-vivo cultured in isotonic (300 mOsmol l− 1) or hypertonic (600 mOsmol l− 1) L-15 medium (penicillin-streptomycin (PS), gentamycin, FBS) at 23 °C for 4, 8 and 24 h. Transcripts levels of ostf1, Na+/Cl−-taurine transporter (taut), Na+/H+ exchange regulatory cofactor (nherf1), cystic fibrosis transmembrane regulator (cftr), and inward rectifying K+ channel (kir), were significantly upregulated at hypertonic medium, in a time-dependent manner (Fig. 1). Moreover, the adenosine triphosphatases (ATPases), including the subunits of Na+/K+-ATPase [ATPase Na+/K+ transporting subunit α1 (at1α1), ATPase Na+/K+ transporting subunit α3 (at1α3)], and calcium-transporting ATPase 2 (at2b2) were increased (p < 0.05). The mRNA expression level of V-type proton ATPase 1 (vpp1) was significantly reduced (p < 0.05) while the levels of aquaporin-3 (aqp-3) and Na+/K+/2Cl− cotransporter (nkcc, s12a2) mRNA showed no noticeable changes (Fig. 2).
For the key factors involved in miR biogenesis, the mRNA expression levels of drosha RNase III endonuclease (drosha), exportin-5 (xpo5), dicer ribonuclease III (dicer-1) and argonaute-2 (ago2) in ex-vivo gill filaments were measured at 4, 8 and 24 h of post-hypertonic treatment (Fig. 3). In general, the expression levels of drosha, and xpo5 transcripts were significantly reduced under hyperosmotic stress. The expression levels of dicer-1, however was significantly increased.
Hyperosmotic treatment induced differential expression of miRs in ex-vivo gill filament model
Gill filaments of eels were ex-vivo cultured in isotonic (300 mOsmol l− 1) or hypertonic (600 mOsmol l− 1) L-15 medium (PS, gentamycin, FBS) at 23 °C for 4 and 8 h. Total RNA was extracted and isolated using mirVana (Invitrogen). The A260/A280 value of the isolated RNA was > 1.8, and the RIN was over 8. Library construction was prepared for miR sequencing. A total of 568 million quality-trimmed raw reads were obtained from the small RNA sequencing (Additional file 1: ST1). De novo analysis identified 658 and 662 miR precursors from the sequencing samples at 4 and 8 h respectively (Table 1). In consideration of Randford p-value, there were 82 and 84 miR precursors found to be not significant at 4 and 8 h treatment respectively. The corrected miR precursors were 576 and 578 at 4 and 8 h treatments. Figures 4a, b show the volcanic plots of deregulated miRs of the samples. In the 4 h treatment, the expression levels of 81 miR precursors were found to be significantly different (DESeq2 adjusted p-value < 0.05). There were 42 miRs upregulated and 39 miRs downregulated. In 8 h treatment, there were 55 differentially expressed miR precursors in which 24 miR precursors were upregulated and 31 were downregulated (DESeq2 adjusted p < 0.05). To prioritize potential differentially regulated miRs, the significance of p-value was set to p < 0.001. This led to a reduction in the numbers of differentially expressed miRs from 81 to 18 at 4 h, and 79 to 10 at 8 h treatments. The Venn diagram (Fig. 4c) shows the common differentially expressed miRs at 4 and 8 h. The two miRs, miR-29b-3p and miR-200b-3p were commonly downregulated at both time points of the hyperosmotic treatments (Fig. 4d). The differentially expressed miRs were validated using real-time PCR of the samples from isotonic and hypertonic treated gill filaments (Fig. 4e).
Since the suppression of gene expression is one of the major functions of miR, a transcriptome analysis of differential gene expression in ex-vivo gill filament-culture was conducted. There were 2085 differentially expressed genes (DEGs) at 4 h of the hypertonic treatment (Fig. 5a), including 890 upregulated genes and 1195 downregulated genes (Additional file 1: ST3). In 8 h of hyperosmotic treatment, 1670 DEGs including 841 upregulated genes and 829 downregulated genes were identified (Fig. 5b, Additional file 1: ST4). In the comparison of DEGs from the 4 h and 8 h treatments, 577 commonly upregulated and 711 downregulated genes were observed (Fig. 5c).
Integrated miR-inhibition and Transcriptome analysis in primary gill culture model
To underpin the gene targets, specific miR-inhibitors were used to block the activities of miR-29b-3p and miR-200b-3p. Figure 6a showed the effects of the inhibition on the individual miRs in primary gill cell culture. The inhibition of miR-29b-3p did not affect the expression level of miR-200b-3p, and vice versa. Supplementary Tables 5 and 6 showed the list of upregulated gene targets in cells after the treatment with miR-29b-3p and miR-200b-3p inhibitors, respectively. We overlaid the transcriptomic data of 4- and 8 h-upregulated genes (ST3 & ST4) and of miR-inhibition, to select two target-genes [i.e., krueppel-like factor 4 (klf4) and homeobox protein meis2] for miR-29b-3p (Fig. 6b), and one target-gene (slc17a5) for miR-200b-3p (Fig. 6c). miRanda algorithm was then used to predict the binding of the individual miRs to the target genes. TransDecoder (version 5.0.2) was used to determine the coding regions including 5′ UTRs and 3′ UTRs of transcripts . The analysis showed that miR-29b-3p and miR-200b-3p could bind to the 3′ untranslated region (3’UTR) of klf4 and slc17a5 (Fig. 6d). The binding site of klf4 is at 3’UTR between 447 to 468, while the binding site of s17a5 is at the 3’UTR between 53 to 73.
miRs are a class of endogenous and conserved small RNA molecules (~ 22 nt) that play gene regulatory roles in targeting mRNAs for cleavage or translational repression . Based on the finding of the recent studies, miRs are found to have a pivotal role in the regulation of fluid and electrolyte balance [18, 19]. It warrants further investigation to identify tonicity-inducible miRs and to determine their potential roles in osmoregulatory responses. In the present study, using ex-vivo gill filament culture model, we characterized the hyperosmotic responses with respective to the expression levels of some well-defined hyper- and hypo-osmotic genes (i.e. regulators and transporters). The involvement of miR-biogenesis genes and tonicity-responsive miRs were characterized. Possible mRNA candidates were predicted using genome-wide analysis of microRNA-mRNA interactome.
In the first part of the study, we characterized the ex-vivo gill filament culture model with regard to its responses to hyper-osmotic stress. The purpose is to establish a culture model, which retains the three-dimensional organization (primary and secondary lamellae) and takes into account of all the cell types [pavement cells (PVCs), chloride cells (CCs), mucous cells and undifferentiated cells] of gill tissues in experiments. In this study, numerous well-characterized hyper-osmotic and hypo-osmotic inducible genes were chosen to evaluate the functional responses and the validity of the culture. The hyperosmotic inducible genes, like ostf1 [20, 21], nherf , wnk1  are regulatory proteins known to modulate gene transcription, apical trafficking of transmembrane G-protein coupled receptors/ion-transporters and epithelial chloride (Cl−) transport respectively. The hyperosmotic ostf1 induction was demonstrated in intact fish and primary gill epithelial cell culture . Its expression could also be stimulated by the seawater-adapting hormone, cortisol . The orchestrating role of ostf1 to integrate environmental and hormonal signals for osmosensory function was documented. In our study, ostf1 mRNA was stimulated in our ex-vivo culture model, under hyperosmotic stress. Our data also showed the upregulation of nherf1, suggesting its function in hyperosmotic adaptation. Intriguingly, it was reported that nherf1 negatively regulated the apical localization of Na+/H+ exchanger (nhe) in renal brush border cells . Nhe facilitates Na+ retention in hypotonic solution, in exchange of H+, the process is important for freshwater adaptation. Nhef1 was found to be expressed in ionocytes in zebrafish embryos . The downregulation of nherf1 blocked Na+ accumulation in ionocytes. Presumably, in our study the upregulation of nherf1 in hyperosmotic gill filaments might downregulate nhe. Hyperosmotic acclimation is known to stimulate Na+ and Cl− secretion in ionocytes via the upregulation of kir, Na+/K+-ATPase (at1a1 and at1a3), and cftr. The sodium pump generated an extracellular Na+ gradient to facilitate the transcellular co-transport of Cl− to cell cytoplasm, followed by Cl− secretion via apical cftr. All these transporters were upregulated in the hyperosmotic treated gill filaments. The calcium-transporting ATPase 2 (at2b2) is known to reduce intracellular Ca2+ to maintain calcium homeostasis, which could be perturbed by osmotic stress and seawater adaptation in eels [27, 28]. The Na+/K+/2Cl−-cotransporter (s12a2) aids in the secondary transport of Na+, K+ and Cl−, driven by Na+-gradient established by Na+/K+-ATPase . However, there was no noticeable induction of sl2a2 in this culture model. In addition to ion secretion, hyperosmotic acclimation involves the accumulation of organic osmolytes, like taurine to increase cellular osmolality. The expression levels of the taurine transporter, taut was upregulated in the hyperosmotic treated gill filaments. The aquaporin-3 (aqp-3) that mediates the translocation of water, glycerol, urea and other small solutes, was highly expressed in freshwater eel gills . No significant reduction in its expression levels was measured in this study. On the other hand, the transporter (vpp1) involved in hypo-osmotic acclimation, was downregulated in the hyperosmotic treated gill filaments. As compared with the primary gill cell model, it is a two-dimensional culture and comprises mostly PVCs [31, 32]. In the past studies, using primary gill PVC culture under hyperosmotic treatment, upregulation of ostf1, Na+/K+-ATPase and taut mRNA expression were detected [31, 33]. The stimulation of cftr, kir, and nherf1 expression in the hyperosmotic-treated culture, has not yet been reported. Indeed, the expression of cftr and kir are known to be cell-specific, localized in seawater CCs. The low percentage of CCs in the primary cell culture made the model not representative as compared with the ex-vivo gill filament culture for studying gill physiology. Collectively, our data demonstrated that the ex-vivo gill filament culture exhibited distinctive responses to hyperosmotic acclimation. It would be a useful model to investigate underlying mechanisms of osmotic responses in fish gills.
With hindsight, the involvement of miRs in fluid and electrolyte transport has become apparent. A number of studies have illustrated the implications of miRs in ion transport in renal and non-renal tissues, showing a pivotal role of miRs in the regulation of fluid and electrolyte balance. In the present study, we studied the expression of four major factors involved in miR biogenesis, using the ex-vivo gill filament culture. Briefly, RNA polymerases II and III are involved in primary-miR (pri-miR) transcription [34, 35]. Upon nuclear cleavage of pri-miR by the Drosha RNase III endonuclease, a ~ 60-70 nt miRNA precursor (pre-miRNA, a stem loop intermediate) is produced and transported from cell nucleus to cytoplasm in a Ran-GTPase dependent manner by the export receptor Exportin-5 [36, 37]. The pre-miRNA is then processed in the cytoplasm by the action of another RNase III endonuclease, Dicer to produce mature miR, followed by Argonaut proteins to target through sequence complementarity in RNA-induced silencing complex . In this study, the mRNA expression levels of drosha, xpo5 and dicer1 were dysregulated under hyperosmotic treatment, indicating that the treatment modulated miR-biogenesis in gill cells. Since miRs take part in various aspects of physiological functions, including tissue remodeling and cell survival, it is anticipated an alternation of miR expression in the hyperosmotic treated gill filaments. Thus, miR sequencing was conducted using the gill filaments at two consecutive time-points. Bioinformatics analysis of the dysregulated miRs from the two time-points at high stringency of filtering (p < 0.001) shortlisted two miR candidates, miR-29b-3p and miR-200b-3p. In the literatures, the family of miR-29 was found to respond to environmental stress and cellular repairing processes. In zebrafish model, the involvement of miR-29 in acute environmental stress (cold stress) was studied , in which the miR targeted on a core clock gene per2 to enhance cold tolerance of the fish larvae. Moreover, miR-29b was reported to play a role in cell regeneration, associated with the processes of cell survival and cytoskeleton reorganization in zebrafish . Furthermore, miR-29 was recognized to modulate iron transport and oxidative stress in neurons of killifish brain . Likewise, miR-200 family comprises of the -a, −b and -c members, those were reported to regulate somatic growth and neurogenesis in zebrafish [42, 43]. A study showed the involvement of miR-200 in growth hormone (GH)/insulin-like growth factor (IGF)-signaling pathway and a high expression level of miR-200 was associated with p53 to induce apoptosis . The similar functional association was found to play roles in sperm motility of zebrafish . Intriguingly, miR-200 was involved in epithelial ion transport to inhibit the expression of nherf1 in ionocytes of zebrafish embryos . Our data agreed with the study to show that the downregulation of miR-200b-3p was accompanied with an upregulation of nherf1 in hyperosmotic gill filaments. Nonetheless, immunological function of miR-200 on the regulation of Toll-like receptor pathways in sea cucumber was described .
In the last part of this study, we integrated the transcriptome data of time-course hyperosmotic treatment and of RNA inhibition, followed by miRanda algorithm analysis to determine the miR-mRNA binding. Klf4 and slc17a5 were identified to be the target-genes of miR-29b-3p and miR-200b-3p, respectively. Klf4 is a zinc-finger transcriptional factor essential for terminal differentiation of intestinal epithelium and skin epidermis [46, 47]. A recent study demonstrated the role of Klf4 in the maintenance of ionocyte progenitor’s population in zebrafish embryos . Chloride cells in fish gills are ionocytes responsible for ion transport. During the course of seawater acclimation, there are significant increase of chloride cells in fish gills for hypo-osmotic adaptation . Presumably, in this study, the hyperosmotic treatment of gill filament culture induced the expression of klf4 to increase chloride cell densities. The miR-200b-3p target-gene, slc17a5 belongs to the solute carrier family of transporters to transport amino acids , which can serve as organic osmolytes to regulate cellular osmolality, against hyperosmotic stress . The implication of the data are consistent with the physiological responses of fish gills.
This is the first study to characterize the hypo-osmotic responses, and the dysregulation of miR-biogenesis factors/miRs in ex-vivo gill filament culture under hyperosmotic challenge. The data showed the culture retained the in-vivo molecular responses to upregulate hyper-osmotic factors/transporters and to downregulate hypo-osmotic transporter. Using multiple experimental approaches, two miRs and two target genes were identified to be involved in the hyperosmotic challenge of gill filaments. The two miRs were predicted to target a cluster of genes involved in osmoregulation. MiR inhibition analysis further identified two gene targets, those are essential for differentiation of gill chloride cells and the regulation of cellular osmolality. There are many other tonicity-responsive miRs and genes were identified in this study. The data support further analysis of other miRs and transcripts to elucidate their regulatory roles in osmotic responses.
Fish maintenance, in-vitro and ex-vivo gill culture model
Japanese eels were obtained from fish vendor Lok Fu market. Japanese eels (600-800 g) was reared in fiberglass tanks supplied with charcoal-filtering aerated freshwater water at 18–20 °C under a 12 L:12D photoperiod for at least 2 weeks to facilitate acclimation before experiments. The eel was handled in accordance to the guidelines and regulation of Hong Kong Baptist University. The fish were anesthetized with 0.1% MS-222 (Sigma) and perfused with phosphate buffered saline (PBS, pH 7.7). Gill tissue was excised and washed for the procedures of ex-vivo gill filament culture. Briefly, gill arches were cut every 5 mm interval with intact primary filaments, washed with PBS, pH 7.7 and maintained in Leibovitz’s L-15 medium (Gibco, Invitrogen, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS, HyClone®, Perbio Sciences, Logan, UT, USA), 1% penicillin/streptomycin, 0.5% fungizone (Gibco, Invitrogen). The culture was incubated at 23 °C.
For primary gill cell culture, the gill filaments were then subject to hyperosmotic-treatment. For primary gill cell culture, gill filaments were cut into small fragments and subjected to two cycles of tryptic digestion (0.5% trypsin + 5.3 mmol l− 1 EDTA) (Sigma), each for 20 min at room temperature in a rotator (300 r.p.m.). Cell suspension was filtered through stainless steel mesh (104 and 73.7 μm, Sigma), washed with PBS, pH 7.7 and finally resuspended in Leibovitz’s L-15 medium (Gibco, Invitrogen, Grand Island, NY, USA) supplemented with 10% fetal bovine serum (FBS, HyClone®, Perbio Sciences, Logan, UT, USA), 1% penicillin/streptomycin, 0.5% fungizone (Gibco, Invitrogen) and seeded at a density of 2 × 106 cells cm− 2 onto collagen coated culture plates (Iwaki, Chiba, Japan). The culture was incubated at 23 °C in a growth chamber with humidified air. One day after seeding, each culture well was rinsed with PBS (pH 7.7) to remove mucous and unattached cells. The culture was subject to miR-knockdown analysis.
Hyperosmotic challenge and RNA preparation
Hypertonic medium was prepared by an addition of 3% of 5 M NaCl to the Leibovitz’s L-15 medium (300 mOsmol l− 1). The osmolality of freshly prepared hypertonic media was measured by a vapor pressure osmometer (Wescor, 5500XR, Logan, UT, USA). For the gill filament culture, tissue samples were exposed to the 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. After the hypertonic treatment, total RNA of gill filaments were extracted using mirVana miRNA isolation kit (Invitrogen) according to the manufacturer’s instruction. RNA quality was assessed using the Agilent 2100 Bioanalyzer system and samples with a RNA Integrity Number (RIN) greater than 8 was used for small RNA and RNA library constructions.
Small RNA sequencing
Twenty small RNA libraries were constructed [five replicates per individual time points (4 and 8 h) per individual osmotic conditions (300 and 600 mOsmol l− 1)] using 5 μg of total RNA from each sample. Short RNA transcripts (18–30 nucleotides long) were resolved and isolated using Polyacrylamide gel electrophoresis (PAGE) gels. The selected small RNA molecules were ligated with 3′ and 5′ adapters, followed by the synthesis of single stranded complementary DNA (cDNA) using SuperScript II Reverse Transcriptase (Invitrogen). The cDNA was then amplified using indexing primers. After quantitative real-time PCR (qRT-PCR) amplification, the library size was determined using the Bioanalyzer (Agilent Technologies 2100). Concentrations of the libraries were tested using qRT-PCR (EvaGreen). The libraries were then processed by the Beijing Genomics Institute (BGI). Single-end read (50 bases pair (bp) read-length) was sequenced using the BGISEQ-500RS sequencer to produce at least 25 M clean reads per sample. Sequencing statistic was shown in Additional file 1: ST1. The adaptor sequences were trimmed and Phred quality score less than 20 were removed. The sequencing reads were mapped to the Anguilla japonica reference genome from NCBI (AVPY00000000.1, https://www.ncbi.nlm.nih.gov/Traces/wgs/?val=AVPY01) using mapper.pl module of the mirDeep2 package . MicroRNAs were identified using the miRDeep2 procedure and final results were annotated for sequence similarity to known miRNAs in miRBase databases . The quantification results produced by the quantifier.pl module of the mirDeep2 package were then used to compare the expression of miRNA in the different time points with the DESeq2 .
Transcriptome sequencing and miRNA target genes prediction
Total RNA extracted from the same set of gill filament samples were subjected for cDNA library construction. The libraries were processed by BGI. Single-end reads, 50 bp read-length, were sequenced on the BGISEQ-500RS sequencer produce at least 25 M clean reads per sample. Sequence-reads were in turn trimmed according to BWA’s – q algorithm  and the sequence reads were quantified using the Anguilla japonica transcriptome assembly (Kallisto version 0.43.0) [56, 57]. Read-count data were then subjected to differential expression analysis using edgeR package . Genes with |log2 (fold change: treatment/control)| > 0.5 and q-value < 0.01 were considered as differentially expressed genes (DEGs). Furthermore, Ingenuity Pathway Analysis (IPA®, QIAGEN Redwood City, www.qiagen.com/ingenuity) was used to decipher the effect of hypertonic treatment in gill filament. In the bioinformatics analysis, all the pathways, diseases or biofunctions with P < 0.05 were considered as statistically significant.
The miRanda algorithm  was used to predict genes targeted the miRNAs as previously described . Briefly, the whole length of the miRNA was first aligned against the mRNA sequences. Thermodynamic stability of the RNA duplexes based on these alignments was calculated. Strict alignment in the seed region at position 2–7 bp of the miRNA to the mRNA sequence is required, of which the region is the key determinant in miRNA–mRNA recognition . The target genes were then subjected to Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.7  to understand the functional roles of hypertonic treatment-deregulated miRNAs.
Real-time PCR analysis of identified miRNAs and differential gene expression under hyperosmotic stress
Gill filaments were exposed to hyperosmotic treatment (600 mOsmol l− 1). After 4 h or 8 h incubation, total RNA of gill filaments were extracted using mirVana miRNA isolation kit (Invitrogen). For mature miRNA detection and quantification, 250 ng of total RNA was reverse transcribed in HiSpec buffer with miScript II RT kit (Qiagen). The cDNA generated is used for RT-PCR with miScript SYBR green PCR kit (Qiagen). Human small nuclear RNA 6 (RNU 6) (Qiagen: MS00033740) was used for normalizing RT-PCR of miRNA. Sequences for primers targeting miR-29b-3p and miR-200b-3p are miR-29b-3p: 5′- UAGCACCAUUUGAAAUCAGUGU-3′; miR-200b-3p: 5′-UAAUACUGCCUGGUAAU GAUG-3′. For differential gene expression analysis, cDNA generated with SuperScript VILO master mix (Invitrogen) was used for qRT-PCR analysis on hyper-osmotic and hypo-osmotic inducible genes. Sequences for specific primers targeting Japanese eel genes are shown in Additional file 1: ST2. Relative quantification of miRNA and mRNA expression was calculated based on delta delta Ct (ΔΔCt) values .
Inhibition of miR activity
For experiments of miR-inhibition (miRi), primary gill cells were transfected with either mirVana™ miRNA inhibitor negative control #1 (Thermo Fisher, Cat #4464076) or mirVana™ target miRNA inhibitors (Thermo Fisher) at 100 nM in L-15 culture medium using Lipofectamine 2000 according to the manufacturer’s instructions (Invitrogen). After 24 h treatment, total RNA of miR-transfected cells was harvested for transcriptome sequencing. Quantitative RT-PCR was performed to assess the efficiency of inhibition. Sequences for miRNA inhibitors targeting miR-29b-3p, and miR-200-3p are 5′-UAGCACC AUUUGAAAUCAGUGU-3′ (Thermo Fisher, Assay ID: MH10432) and 5′-UAAUACUGC CUGGUAAUGAUG-3′ (Thermo Fisher, Assay ID: MH11741), respectively.
Statistical analyses of DEGs and differentially expressed miRNA profiles were performed using GraphPad Prism 6 (GraphPad Software, San Diego, CA). Student’s t-test was used to assess for potential differences between control (isotonic) and treatment (hypertonic) groups. We reported statistical significance at *P < 0.05, **P < 0.01 and ***P < 0.001. Values were expressed as mean ± S.E.M.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article and its additional files.
ATPase Na+/K+ transporting subunit a1
ATPase Na+/K+ transporting subunit a3
Calcium-transporting ATPase 2
Beijing Genomics Institute
Cystic fibrosis transmembrane regulator
Differentially expressed genes
Dicer ribonuclease III
Drosha RNase III endonuclease
Fetal bovine serum
Insulin-like growth factor
Inward rectifying K+ channel
Kruppel-like factor 4
Na+/K+-exchanger regulatory cofactor
- Nkcc, s12a2:
Osmotic response element binding protein
Osmoregulatory transcription factor
Polyacrylamide gel electrophoresis
RNA Integrity Number
- S17a5, Slc17a5:
Solute carrier family 17 member 5
V-type proton ATPase 1
Evans DH. Teleost fish osmoregulation: what have we learned since August Krogh, Homer Smith, and Ancel Keys. Am J Physiol Regul Integr Comp Physiol. 2008;295(2):R704–13.
Evans DH, Piermarini PM, Choe KP. The multifunctional fish gill: dominant site of gas exchange, osmoregulation, acid-base regulation, and excretion of nitrogenous waste. Physiol Rev. 2005;85(1):97–177.
Hwang PP, Lee TH, Lin LY. Ion regulation in fish gills: recent progress in the cellular and molecular mechanisms. Am J Physiol Regul Integr Comp Physiol. 2011;301(1):R28–47.
Takei Y, Hiroi J, Takahashi H, Sakamoto T. Diverse mechanisms for body fluid regulation in teleost fishes. Am J Physiol Regul Integr Comp Physiol. 2014;307(7):R778–92.
Pietrzykowski AZ, Friesen RM, Martin GE, Puig SI, Nowak CL, Wynne PM, et al. Posttranscriptional regulation of BK channel splice variant stability by miR-9 underlies neuroadaptation to alcohol. Neuron. 2008;59(2):274–87.
Xiao L, Xiao J, Luo X, Lin H, Wang Z, Nattel S. Feedback remodeling of cardiac potassium current expression: a novel potential mechanism for control of repolarization reserve. Circulation. 2008;118(10):983–92.
Yang B, Lin H, Xiao J, Lu Y, Luo X, Li B, et al. The muscle-specific microRNA miR-1 regulates cardiac arrhythmogenic potential by targeting GJA1 and KCNJ2. Nat Med. 2007;13(4):486–91.
Ceolotto G, Papparella I, Bortoluzzi A, Strapazzon G, Ragazzo F, Bratti P, et al. Interplay between miR-155, AT1R A1166C polymorphism, and AT1R expression in young untreated hypertensives. Am J Hypertens. 2011;24(2):241–6.
Sober S, Laan M, Annilo T. MicroRNAs miR-124 and miR-135a are potential regulators of the mineralocorticoid receptor gene (NR3C2) expression. Biochem Biophys Res Commun. 2010;391(1):727–32.
Elvira-Matelot E, Zhou XO, Farman N, Beaurain G, Henrion-Caude A, Hadchouel J, et al. Regulation of WNK1 expression by miR-192 and aldosterone. J Am Soc Nephrol. 2010;21(10):1724–31.
Huang W, Liu H, Wang T, Zhang T, Kuang J, Luo Y, et al. Tonicity-responsive microRNAs contribute to the maximal induction of osmoregulatory transcription factor OREBP in response to high-NaCl hypertonicity. Nucleic Acids Res. 2011;39(2):475–85.
Flynt AS, Thatcher EJ, Burkewitz K, Li N, Liu Y, Patton JG. miR-8 microRNAs regulate the response to osmotic stress in zebrafish embryos. J Cell Biol. 2009;185(1):115–27.
Yan B, Zhao LH, Guo JT, Zhao JL. miR-429 regulation of osmotic stress transcription factor 1 (OSTF1) in tilapia during osmotic stress. Biochem Biophys Res Commun. 2012;426(3):294–8.
Yan B, Guo JT, Zhao LH, Zhao JL. MiR-30c: a novel regulator of salt tolerance in tilapia. Biochem Biophys Res Commun. 2012;425(2):315–20.
Wang X, Yin D, Li P, Yin S, Wang L, Jia Y, et al. MicroRNA-sequence profiling reveals novel Osmoregulatory MicroRNA expression patterns in Catadromous eel Anguilla marmorata. PLoS One. 2015;10(8):e0136383.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136(2):215–33.
Chandrasekaran K, Karolina DS, Sepramaniam S, Armugam A, Wintour EM, Bertram JF, et al. Role of microRNAs in kidney homeostasis and disease. Kidney Int. 2012;81(7):617–27.
Elvira-Matelot E, Jeunemaitre X, Hadchouel J. Regulation of ion transport by microRNAs. Curr Opin Nephrol Hypertens. 2011;20(5):541–6.
Chow SC, Wong CK. Regulatory function of hyperosmotic stress-induced signaling cascades in the expression of transcription factors and osmolyte transporters in freshwater Japanese eel primary gill cell culture. J Exp Biol. 2011;214(Pt 8):1264–70.
Fiol DF, Kultz D. Rapid hyperosmotic coinduction of two tilapia (Oreochromis mossambicus) transcription factors in gill cells. Proc Natl Acad Sci U S A. 2005;102(3):927–32.
Li J, Callaway DJ, Bu Z. Ezrin induces long-range interdomain allostery in the scaffolding protein NHERF1. J Mol Biol. 2009;392(1):166–80.
Shekarabi M, Zhang J, Khanna AR, Ellison DH, Delpire E, Kahle KT. WNK kinase signaling in ion homeostasis and human disease. Cell Metab. 2017;25(2):285–99.
Fiol DF, Chan SY, Kultz D. Regulation of osmotic stress transcription factor 1 (Ostf1) in tilapia (Oreochromis mossambicus) gill epithelium during salinity stress. J Exp Biol. 2006;209(Pt 16):3257–65.
Chow SC, Tse WK, Wong CK. Dexamethasone (DEX) induces osmotic stress transcription factor 1 (Ostf1) through the Akt-GSK3beta pathway in freshwater Japanese eel gill cell cultures. Biol Open. 2013;2(5):487–91.
Weinman EJ, Steplock D, Donowitz M, Shenolikar S. NHERF associations with sodium-hydrogen exchanger isoform 3 (NHE3) and ezrin are essential for cAMP-mediated phosphorylation and inhibition of NHE3. Biochemistry. 2000;39(20):6123–9.
Gu J, Law AY, Yeung BH, Wong CK. Activation of gill Ca2+−sensing receptor as a protective pathway to reduce Ca2+−induced cytotoxicity. J Mol Endocrinol. 2014;53(2):155–64.
Gu J, Law AY, Yeung BH, Wong CK. Characterization of stanniocalcin 1 binding and signaling in gill cells of Japanese eels. J Mol Endocrinol. 2015;54(3):305–14.
Tse WK, Au DW, Wong CK. Characterization of ion channel and transporter mRNA expressions in isolated gill chloride and pavement cells of seawater acclimating eels. Biochem Biophys Res Commun. 2006;346(4):1181–90.
Cutler CP, Cramb G. Branchial expression of an aquaporin 3 (AQP-3) homologue is downregulated in the European eel Anguilla anguilla following seawater acclimation. J Exp Biol. 2002;205(Pt 17):2643–51.
Tse WK, Au DW, Wong CK. Effect of osmotic shrinkage and hormones on the expression of Na+/H+ exchanger-1, Na+/K+/2Cl- cotransporter and Na+/K+ −ATPase in gill pavement cells of freshwater adapted Japanese eel, Anguilla japonica. J Exp Biol. 2007;210(Pt 12):2113–20.
Tse WK, Chow SC, Wong CK. The cloning of eel osmotic stress transcription factor and the regulation of its expression in primary gill cell culture. J Exp Biol. 2008;211(Pt 12):1964–8.
Chow SC, Ching LY, Wong AM, Wong CK. Cloning and regulation of expression of the Na+−cl--taurine transporter in gill cells of freshwater Japanese eels. J Exp Biol. 2009;212(Pt 20):3205–10.
Dieci G, Fiorino G, Castelnuovo M, Teichmann M, Pagano A. The expanding RNA polymerase III transcriptome. Trends Genet. 2007;23(12):614–22.
Lin SL, Kim H, Ying SY. Intron-mediated RNA interference and microRNA (miRNA). Front Biosci. 2008;13:2216–30.
Lee Y, Jeon K, Lee JT, Kim S, Kim VN. MicroRNA maturation: stepwise processing and subcellular localization. EMBO J. 2002;21(17):4663–70.
Lund E, Guttinger S, Calado A, Dahlberg JE, Kutay U. Nuclear export of microRNA precursors. Science. 2004;303(5654):95–8.
Lee Y, Ahn C, Han J, Choi H, Kim J, Yim J, et al. The nuclear RNase III Drosha initiates microRNA processing. Nature. 2003;425(6956):415–9.
Hung IC, Hsiao YC, Sun HS, Chen TM, Lee SJ. MicroRNAs regulate gene plasticity during cold shock in zebrafish larvae. BMC Genomics. 2016;17(1):922.
Fuller-Carter PI, Carter KW, Anderson D, Harvey AR, Giles KM, Rodger J. Integrated analyses of zebrafish miRNA and mRNA expression profiles identify miR-29b and miR-223 as potential regulators of optic nerve regeneration. BMC Genomics. 2015;16:591.
Ripa R, Dolfi L, Terrigno M, Pandolfini L, Savino A, Arcucci V, et al. MicroRNA miR-29 controls a compensatory response to limit neuronal iron accumulation during adult life and aging. BMC Biol. 2017;15(1):9.
Choi PS, Zakhary L, Choi WY, Caron S, Alvarez-Saavedra E, Miska EA, et al. Members of the miRNA-200 family regulate olfactory neurogenesis. Neuron. 2008;57(1):41–55.
Jing J, Xiong S, Li Z, Wu J, Zhou L, Gui JF, et al. A feedback regulatory loop involving p53/miR-200 and growth hormone endocrine axis controls embryo size of zebrafish. Sci Rep. 2015;5:15906.
Xiong S, Ma W, Jing J, Zhang J, Dan C, Gui JF, et al. An miR-200 cluster on chromosome 23 regulates sperm motility in Zebrafish. Endocrinology. 2018;159(5):1982–91.
Lv Z, Li C, Zhang P, Wang Z, Zhang W, Jin CH. MiR-200 modulates coelomocytes antibacterial activities and LPS priming via targeting Tollip in Apostichopus japonicus. Fish Shellfish Immunol. 2015;45(2):431–6.
Chen X, Whitney EM, Gao SY, Yang VW. Transcriptional profiling of Kruppel-like factor 4 reveals a function in cell cycle regulation and epithelial differentiation. J Mol Biol. 2003;326(3):665–77.
Segre JA, Bauer C, Fuchs E. Klf4 is a transcription factor required for establishing the barrier function of the skin. Nat Genet. 1999;22(4):356–60.
Chen YC, Liao BK, Lu YF, Liu YH, Hsieh FC, Hwang PP, et al. Zebrafish Klf4 maintains the ionocyte progenitor population by regulating epidermal stem cell proliferation and lateral inhibition. PLoS Genet. 2019;15(4):e1008058.
Wong CK, Chan DK. Isolation of viable cell types from the gill epithelium of Japanese eel Anguilla japonica. Am J Phys. 1999;276(2):R363–72.
Sreedharan S, Shaik JH, Olszewski PK, Levine AS, Schioth HB, Fredriksson R. Glutamate, aspartate and nucleotide transporters in the SLC17 family form four main phylogenetic clusters: evolution and tissue expression. BMC Genomics. 2010;11:17.
Wiesenthal AA, Muller C, Harder K, Hildebrandt JP. Alanine, proline and urea are major organic osmolytes in the snail Theodoxus fluviatilis under hyperosmotic stress. J Exp Biol. 2019;222(Pt 3):jeb193557.
Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ. miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006;34(Database issue):D140–4.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Wang SY, Lau K, Lai KP, Zhang JW, Tse AC, Li JW, et al. Hypoxia causes transgenerational impairments in reproduction of fish. Nat Commun. 2016;7:12114.
Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34(5):525–7.
Lai KP, Li JW, Gu J, Chan TF, Tse WK, Wong CK. Transcriptomic analysis reveals specific osmoregulatory adaptive responses in gill mitochondria-rich cells and pavement cells of the Japanese eel. BMC Genomics. 2015;16:1072.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in drosophila. Genome Biol. 2003;5(1):R1.
Lai KP, Li JW, Tse AC, Chan TF, Wu RS. Hypoxia alters steroidogenesis in female marine medaka through miRNAs regulation. Aquat Toxicol. 2016;172:1–8.
Lambert NJ, Gu SG, Zahler AM. The conformation of microRNA seed regions in native microRNPs is prearranged for presentation to mRNA targets. Nucleic Acids Res. 2011;39(11):4827–35.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Keng Po Lai is supported by Hong Kong SAR, Macao SAR and Taiwan Province Talent Young Scientist Program of Guangxi.
This whole project including the sample preparation, sequencing and data analysis was supported by General Research Fund (HKBU12162016), Research Grant Council Hong Kong, to Dr. Chris KC Wong (Hong Kong Baptist University).
Ethics approval and consent to participate
All fish experiments were conducted following ethical principles and procedures approved by the guidelines and regulation at Hong Kong Baptist University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Clean read stats: Small RNA sequencing statistic; Table S2. Primer sequence for RT-PCR: List of primer used for Real-time PCR analysis. Table S3. mRNA seq-4 h hyper: List of DEGs identified in eel gill filament culture under 4 h hypertonic stress treatment. Table S4. mRNAseq-8 h hyper: List of DEGs in eel gill filament culture under 8 h hypertonic stress treatment. Table S5. mRNAseq-miR-29b inhibition. List of upregulated genes in primary gill cell culture, after miR-29b-3p inhibition. Table S6. mRNAseq-miR-200b inhibition: List of upregulated genes in primary gill cell culture, after miR-200b-3p inhibition.
About this article
Cite this article
Ng, H.M., Ho, J.C.H., Nong, W. et al. Genome-wide analysis of MicroRNA-messenger RNA interactome in ex-vivo gill filaments, Anguilla japonica. BMC Genomics 21, 208 (2020). https://doi.org/10.1186/s12864-020-6630-0
- Fish gill osmoregulation
- Hyperosmotic stress
- Japanese eel
- microRNA inhibitors