Thermal plasticity of the miRNA transcriptome during Senegalese sole development
© Campos et al.; licensee BioMed Central Ltd. 2014
Received: 20 October 2013
Accepted: 17 June 2014
Published: 25 June 2014
Several miRNAs are known to control myogenesis in vertebrates. Some of them are specifically expressed in muscle while others have a broader tissue expression but are still involved in establishing the muscle phenotype. In teleosts, water temperature markedly affects embryonic development and larval growth. It has been previously shown that higher embryonic temperatures promoted faster development and increased size of Senegalese sole (Solea senegalensis) larvae relatively to a lower temperature. The role of miRNAs in thermal-plasticity of growth is hitherto unknown. Hence, we have used high-throughput SOLiD sequencing to determine potential changes in the miRNA transcriptome in Senegalese sole embryos that were incubated at 15°C or 21°C until hatching and then reared at a common temperature of 21°C.
We have identified 320 conserved miRNAs in Senegalese sole, of which 48 had not been previously described in teleosts. mir-17a-5p, mir-26a, mir-130c, mir-206-3p, mir-181a-5p, mir-181a-3p and mir-199a-5p expression levels were further validated by RT- qPCR. The majority of miRNAs were dynamically expressed during early development, with peaks of expression at pre-metamorphosis or metamorphosis. Also, a higher incubation temperature (21°C) was associated with expression of some miRNAs positively related with growth (e.g., miR-17a, miR-181-5p and miR-206) during segmentation and at hatching. Target prediction revealed that these miRNAs may regulate myogenesis through MAPK and mTOR pathways. Expression of miRNAs involved in lipid metabolism and energy production (e.g., miR-122) also differed between temperatures. A miRNA that can potentially target calpain (miR-181-3p), and therefore negatively regulate myogenesis, was preferentially expressed during segmentation at 15°C compared to 21°C.
Temperature has a strong influence on expression of miRNAs during embryonic and larval development in fish. Higher expression levels of miR-17a, miR-181-5p and miR-206-3p and down-regulation of miR-181a-3p at 21°C may promote myogenesis and are in agreement with previous studies in Senegalese sole, which reported enhanced growth at higher embryonic temperatures compared to 15°C. Moreover, miRNAs involved in lipid metabolism and energy production may also contribute to increased larval growth at 21°C compared to 15°C. Taken together, our data indicate that miRNAs may play a role in temperature-induced phenotypic plasticity of growth in teleosts.
MicroRNAs (miRNAs) are short non-coding regulatory RNAs that regulate gene expression post-transcriptionally. Mature miRNAs are often conserved across a wide range of species, and mutations in proteins required for miRNA function or biogenesis have shown to impair animal development . They are involved in the majority of physiological processes, including stem cell differentiation, cell lineage specification, haematopoiesis, neurogenesis, immune response and myogenesis . MiRNA-mediated gene regulation involves binding to the 3’-UTR of their mRNA target and sequestration of mRNA by the miRNA/RISC complex (RNA-induced silencing complex) in the cytoplasmic processing bodies, which results in repression of translation initiation or mRNA degradation reviewed in .
In vertebrates, somatic and skeletal muscle growth is strongly stimulated by Igf-I (Insulin-like growth factor I) , whose functions are mediated by the Igf-I receptor (Igf1R) through activation of two major intracellular signalling pathways: mitogen-activated protein kinases (MAPKs) and mTOR (mammalian target of rapamycin) through phosphatidylinositol 3 kinase (PI3K)/Akt . mTOR mediates signalling in response to nutrient availability, cellular energy, mitogenic signals and various types of stressors. Phosphorylation of mTOR increases levels of protein synthesis by regulating essential proteins controlling mRNA translation. The MAPK pathway is essential for muscle cell proliferation and differentiation [5, 6], and the closely related hepatic insulin pathway acts upstream of Akt/mTOR. It recruits IRS proteins, which in turn recruit PI3K to phosphorylate Akt. Downstream of Akt, subsequent metabolic effects are largely mediated by mTOR and also FoxO1 . mTOR then regulates S6K1 and 4-EBP1 to induce hepatic protein synthesis and to stimulate lipogenesis .
The physiological actions of these complex pathways are multiple and several miRNAs have already been identified as intervenients on their regulation. For instance, in human, miR-199a-3p and miR-101 can suppress gene expression by directly binding to the 3’-UTR of mTOR [8, 9]. Furthermore, miR-181 regulates PTEN, the principal negative regulator of PI3K . MiRNAs that influence mTOR and MAPK pathways play crucial roles in myogenesis, since they are essential for muscle development and growth. For example, miR-181 is involved in myoblast differentiation and in establishing the muscle phenotype . Additionally, some miRNAs (e.g., miR-1, miR-133 and miR-206) are specifically expressed in muscle and positively related with its growth . Interestingly, in teleost like zebrafish (Danio rerio), fast muscle miRNA populations are differentially regulated during the transition from hyperplastic to hypertrophic muscle phenotype in adult fish . Also, the dynamic expression of miR-206, miR-26a and miR-214 throughout skeletal muscle development in the common carp (Cyprinus carpio), points to an active role of these miRNAs in the myogenic process .
miRNAs can also modulate protein abundance and/or phosphorylation status of key signaling components of the insulin pathway. Expression of miR-122  and miR-33  is inversely correlated with expression and/or phosphorylation status of the α-subunit of the metabolic sensor AMPK. In fact, miR-122, which is highly abundant in mammalian  and fish liver , exerts a stimulatory role on lipogenesis and cholesterol synthesis, and inhibits fatty-acid oxidation and therefore energy production .
In fish, changes in miRNA expression have been documented during early ontogeny [19, 20] or in response to food ingestion . However, the potential effect of water temperature on the miRNA transcriptome has never been investigated in teleosts, in spite of its profound impact on larval and juvenile growth . The flatfish Senegalese sole (Solea senegalensis) is a species of commercial interest for marine aquaculture, particularly in Southern Europe, and it can be exposed to large thermal variations in natural and aquaculture conditions . In a previous study it was shown that higher embryonic temperatures (18 and 21°C) promoted a faster development and increased the size of sole larvae by 30 days post hatch (dph) relatively to a lower temperature (15°C) . Furthermore, muscle cellularity at equivalent developmental stages was affected by temperature and a transient differential gene expression associated with incubation temperature was also observed at several stages. For example, myf5, fst, myHC and mrf4 transcript levels were higher at 21°C compared to 15°C . It is likely that this differential expression of myogenic genes is regulated by miRNAs to some extent but their role in thermal plasticity of development and growth remains to be investigated. Hence, we have used the SOLiD deep-sequencing platform to determine potential changes in the miRNA transcriptome in Senegalese sole embryos and larvae exposed to two different temperatures (15°C or 21°C) during embryonic development. Moreover, we have predicted targets of miRNAs that were differentially expressed with temperature, in order to understand their potential role in the phenotypic plasticity of growth in teleosts.
SOLiD sequencing of small RNAs
Summary of the sequenced miRNA transcriptomes in Senegalese sole embryos and larvae incubated at 15°C or 21°C
# Total reads
# Trimmed (% Total)
Average length (nt)
# Annotated miRBase (% Total)
# Conserved miRNAs
6 642 089
2 802 986 (42.2)
200 224 (3.0)
8 842 550
2 169 069 (24.5)
174 027 (2.0)
6 566 636
3 320 266 (50.6)
993 981 (15.1)
6 101 835
2 564 327 (42.0)
689 954 (11.3)
7 268 950
4 312 425 (59.3)
1 467 039 (20.2)
6 200 277
3 256 433 (52.5)
1 162 239 (18.7)
7 105 697
4 965 849 (69.9)
2 391 674 (33.7)
6 325 932
4 247 722 (67.2)
804 587 (12.7)
Differential expression of conserved miRNAs during Senegalese sole development and with embryonic temperature
There were several significant differences in miRNA expression between embryonic temperatures. miR-17a-5p, miR-26a, miR-130c and miR-206-3p transcripts were significantly higher at 21°C than at 15°C during 75Ep (P < 0.05) (Figure 3A, B, C, G). At the 20S stage, miR-130c and miR-181a-3p showed a 1.5- and 4.9-fold higher expression at 15°C than 21°C (P < 0.05), respectively (Figure 3C, E). At hatching, miR-17a-5p, miR-26a, miR-130c, miR-181a-5p, miR-199a-5p and miR-206-3p were highly expressed at 21°C compared to 15°C (Figure 3A, B, C, D, F, G; P < 0.05). No significantly differences between temperatures were found after hatching.
miRNA target prediction in Senegalese sole
Predicted targets of five differentially expressed miRNAs during early development of Senegalese sole
Ensembl Gene ID
Characterisation of Senegalese sole miRNAs
Using high-throughput SOLiD sequencing, 320 conserved miRNA guide and passenger strands were identified and their expression profiles analysed in Senegalese sole embryos and larvae subjected to two different incubation temperatures (15°C and 21°C). There was some tendency to observe a higher number of conserved miRBase reads at later developmental stages and at a lower water temperature. This is unlikely to be a technical bias, since all libraries were prepared simultaneously and run on the same slide. Instead, it reflects a dynamic miRNA transcriptome throughout early ontogeny, which is influenced by environmental factors. From the 320 conserved miRNAs, 48 were not previously found in fish but were known in mammals, birds or nematodes. Up- or down-regulation of some of them is known to be involved in cancer-associated gene deregulation, namely miR-211  and miR-509 . miR-71 mediates the effects of germ cell loss on life span, and its over-expression was found to extend the life span of Caenorhabditis elegans lacking germ cells . Most of these novel miRNAs were found in more than one small RNA library and were regulated, which indicates that they may play a key role in Senegalese sole development. Nevertheless, the specific functions of these miRNAs in embryogenesis remain to be ascertained.
There is increasing evidence that miRNAs that possess both highly conserved -5p and -3p sequences, especially within seed and anchor sequences, can originate two mature functional miRNAs, such as the case of miR-18a, miR-140 or miR-17 . In the present study, several passenger strands were expressed at higher levels than the respective guide strands at one or more developmental stages (Figure 2 and 3), which suggests some involvement of these strands in Senegalese sole development.
Expression profile and thermal plasticity of conserved miRNAs during Senegalese sole development
Critical early ontogeny events like organogenesis, hatching or metamorphosis in fish involve dramatic changes in signalling, physiology and morphology [28–31]. Well defined expression windows of myogenic and growth-related genes during Senegalese sole ontogeny and larval development have been previously characterised, since different genes are activated or down-regulated at different stages of development [23, 32]. It is plausible that allometric differences in growth of specific tissues may influence the observed miRNA expression profiles. Nevertheless, substantial changes in the Senegalese sole miRNA transcriptome during development (Table 1) highlight the importance of miRNAs in tissue and organ differentiation. miRNA differential expression has been found during the development of other vertebrate species, and similarly to our work, generally an increasing number of miRNAs are observed from the earliest to the latest stages of development [19, 33]. Furthermore, we found that whilst some Senegalese sole miRNAs were expressed throughout early ontogeny, others were stage-specific. The miR-430 family is an example of the latter (highly expressed only in earlier stages) and it is well known to accelerate deadenylation and clearance of maternal transcripts during zygotic stages [34, 35]. miRNAs mostly observed after hatching or during metamorphosis included miR-125, miR-199, miR-192, miR-214 and the let-7 family (Figure 2). The latter comprises a family of miRNAs that are highly conserved from worms to humans  and are involved in various processes, such as regulating cell proliferation and differentiation during development.
Several miRNAs showed an expression peak during pre-metamorphosis or metamorphosis (Figures 2 and 3). For instance, amongst the miRNAs validated by qPCR, miR-17a-5p, miR-26a, miR-130c,miR-181a-5p and miR-181a-3p showed a peak of expression mainly during the pre-metamorphic stage (Figure 3A, B, C, D, E), whereas, miR-199a-5p and miR-206-3p had higher levels at metamorphosis or later (Figure 3F, G). Similar results were found for some miRNAs in the olive flounder (Paralichthys olivaceus) . Flatfish larvae undergo a particularly impressive metamorphosis, resulting in an asymmetrical craniofacial remodelling and lateralized behaviour [38, 39]. In olive flounder, an involvement of miR-206a has also been suggested in the metamorphic process, which is mediated by the thyroid hormone . It is known that Senegalese sole larvae have a high growth rate and accumulate energetic compounds until the onset of metamorphosis  but growth rates decrease significantly once metamorphosis starts . It is possible that up-regulation of some miRNAs at a pre-metamorphic/metamorphic stages is associated with changes in lipid metabolism and with the complex metamorphic process. Such alterations in lipid metabolism may also be correlated with the high expression of miR-122 during metamorphosis, as previously observed in other flatfish like halibut  or flounder . miR-122, which is a liver-specific miRNA, has a crucial role in lipid metabolism in teleosts . It is therefore plausible that changes in this miRNA abundance might be related with physiological changes occurring during metamorphosis.
Some of the miRNAs identified in Senegalese sole are known to play a stimulatory role in somatic growth and specifically in skeletal muscle development. In mammals, miR-26a up-regulation is described to post-transcriptionally repress the histone methyltransferase Ezh2, which is a suppressor of skeletal muscle cell differentiation . In fish, such as the common carp, this miRNA was also associated with myogenesis . In our study, the higher miR-26a expression observed in the 21°C group during hatching (and a similar trend during 75Ep and 20S) indicates a possible activation of the myogenic programme at a higher temperature. A similar conclusion can be drawn from miR-181a-5p expression in Senegalese sole embryos. miR-181 is strongly up-regulated in regenerating muscle from an in vivo mouse model of muscle injury  and endurance exercise was also found to significantly increase the expression of miR-181 in mice skeletal muscle . These results are consistent with previous findings in Senegalese sole, where a lower incubation temperature (15°C) produced lighter larvae with smaller fast fibres .
The muscle-specific miR-206, known to reinforce the myogenic program by inhibiting the expression of DNA polymerase α (responsible for cell proliferation) and indirectly down-regulating the MyoD inhibitors Id1-3 and MyoR , showed increased expression levels during embryogenesis at 21°C compared to 15°C, further supporting our observations that this temperature enhanced muscle growth [23, 47]. Interestingly, MyoD was found to regulate the expression of miR-133  and miR-206 . Transcript levels of myod2 at the 20S stage in Senegalese sole embryos incubated at 15°C or 21°C  are consistent with the expression of miR-206 observed in the present study, which is higher at 21°C than 15°C. A similar conclusion can be drawn for miR-133, based on the digital expression profiling from SOLiD sequencing. In zebrafish and medaka (Oryzias latipes), miR-17 is ubiquitously expressed during embryonic development but becomes restricted to proliferative tissues later in development . In the present study, we observed increased expression of miR-17a-5p in the 21°C group at some developmental stages (Figure 3A); whether this is related with increased cellular proliferation is however unknown. Activation of the miR-214/199 cluster coincides with the recruitment of myoD and myogenin . Furthermore, miR-199, miR-206 and miR-17 are positively associated with muscle regeneration and remodelling . This ability to stimulate myogenesis may explain the differences in transcript levels of miR-199, miR-206 and miR-17 between the two temperature groups of Senegalese sole larvae at hatching. Higher levels of these miRNAs are consistent with the observed muscle phenotype in the 21°C group compared to larvae reared at 15°C .
The observed growth differences between temperature groups may also be associated with other metabolic functions (e.g., energy production), since the abundance of miRNAs related to lipid metabolism differed across temperature groups. For instance, the guide strand of the liver-specific miR-122 was expressed at higher levels at 15°C than at 21°C. Given the role of this miRNA in fatty-acid oxidation, and therefore in energy production, it is plausible that temperature affected lipid metabolism and energy utilisation of sole larvae, thus affecting growth. Furthermore, since the passenger strand was expressed at high levels, one can hypothesise that it may also play an active role in lipid metabolism.
Prediction of potential targets of the differentially regulated Senegalese sole miRNAs
Several miRNAs were computationally predicted with high confidence to target mRNAs directly related to myogenesis (Figure 4, Table 2). However, it should be noted that miRNAs can have multiple mRNAs targets and whole larvae were used to investigate miRNA expression in this study. Hence, it is likely that expression levels of some Senegalese sole miRNAs reflect a broader gene regulation, except for those that are muscle-specific, such as miR-206.
mTOR and p38 MAPK pathways have overlapping yet distinct roles during myogenesis: both are required for myogenic differentiation and p38 MAPK, but not mTOR, can inhibit proliferation . It was also suggested that the normoxia activation of MAPK/ERK signalling pathway not only stimulates myoblast proliferation  but also suppresses myogenic differentiation . MAPKs can be activated by a wide variety of different stimuli, but in general MAPK1 and 3 (also called ERK2 and ERK1, respectively) are preferentially activated in response to growth factors whereas the JNK and p38 kinases (which include amongst others MAPK9 and MAPK13) are more responsive to stress stimuli like osmotic shock or ionizing radiation . The major direct or indirect targets of mTOR appear to be mechanisms responsible for ribosome recruitment to mRNA, such as the eukaryotic translation initiation factor 4E (eIF4E), which is activated by the phosphorylation of its repressors, the 4EBP proteins [54, 55]. Phosphorylation of 4EBP1 by mTOR stimulates translation initiation through the release of eIF4E from 4EBP1, allowing eIF4E to associate with eIF4G to enhance cap-dependent translation .
The DNA-damage-inducible transcript 4 (REDD1) suppresses mTOR activity through the TSC1/TSC2 heterodimer [56, 57]. miR-181a-5p up-regulation in Senegalese sole hatchlings at 21°C may down-regulate this mTOR inhibitor. miR-181a-5p was also predicted to target the 3’ UTR of mapk3. Thus, the up-regulation of miR-181a-5p at 21°C during early development may have led to higher myogenic differentiation of this group. Mapk3 and mapk13 are potential targets of miR-17a-5p, which are involved in MAPK/ERK and JNK/p38 MAPK signalling pathways, respectively. In human cell lines, miR-17-5p targets the 3’UTR of Mapk9 . MAPK9 is a negative regulator of cellular proliferation through the JNK/P38 MAPK pathway  and this interaction contributed to the proliferative phenotype caused by miR-17-5p . In the present study, miR-130c was predicted to target MAPK9 mRNA and may therefore have a role in regulating cell proliferation at some developmental stages.
miR-206 was predicted to target sestrin1, which codes for a negative regulator of mTOR signalling through activation of AMPK and TSC2 phosphorylation . The up-regulation of miR-206 at 21°C during early development may have down-regulated this myogenic inhibitor. Furthermore, sestrin3 was predicted to be regulated by miR-17a. Like Sestrin1 and Sestrin2, Sestrin3 promotes AMPK activity [60, 61] and it has been shown that the induced expression of Sestrin3 by FoxO1 transcription factor inhibits mTOR Complex 1 (mTORC1), which regulates cell growth mostly through increased protein synthesis . miR-17a up-regulation at hatching in larvae reared at 21°C is consistent with a possible down-regulation of sestrin3 and a concomitant increase in growth.
Calpains are known to play a role in myogenesis, namely in promoting myoblast fusion . Senegalese sole miR-181a-3p was predicted with high confidence to target the 3’ UTR of calpain 1. Its 4.9-fold up-regulation at 15°C during somitogenesis indicates a possible repression of calpain 1, which would impair the fusion of myoblasts to multinucleated myotubes at this temperature and result in the observed muscle phenotype differences between the two temperature groups .
We have identified 320 conserved miRNAs in Senegalese sole, most of them dynamically expressed during early development and 48 not previously described in teleosts. Embryonic temperature affected the expression of several miRNAs, and it seems that a higher incubation temperature (21°C) promoted expression of miRNAs positively related with growth at specific developmental stages. In particular, higher expression levels of miR-17a, miR-181-5p and miR-206 and down-regulation of miR-181a-3p at 21°C may promote myogenesis and are in agreement with previous studies in Senegalese sole, which reported enhanced growth at higher embryonic temperatures compared to 15°C. Also, miRNAs related to lipid metabolism and energy production (e.g., miR-122) may also be involved in the differential growth of fish with temperature. Taken together, our data indicate that miRNAs may play a role in temperature-induced phenotypic plasticity of growth in teleosts.
All experiments were conducted by trained scientists following FELASA category C recommendations and were carried out in accordance with the clear boundaries of EU legal frameworks, specifically those relating to the protection of animals used for scientific purposes (i.e. Directive 2010/63/EU) and under the Portuguese legislation regarding the protection of animals used for scientific purposes (Decree-Law No. 113/2013). This study was performed at the research facilities of the CCMAR (Algarve, Portugal), which has been certified for animal experimentation by the competent authority (Direcção Geral de Alimentação e Veterinária, Portugal).
Senegalese sole embryos were incubated at two temperatures (15.2 ± 0.3°C and 20.9 ± 0.3°C) in triplicate groups at the LEOA facility, CCMAR/University of Algarve, Portugal. Eggs obtained by natural spawning of captive broodstock fish were randomly distributed amongst 6 fibreglass conical tanks (100 l) at a density of 100 eggs · l-1. After hatching all larvae were reared at the same temperature (20.8 ± 0.1°C) until the metamorphosis climax. Water temperature, O2, salinity, pH and nitrogenous compounds were monitored regularly during the entire trial, and larvae were exposed to an artificial photoperiod of 12 h:12 h light:dark. They were reared in 100 L conical tanks at a density of 40 larvae · L-1. At mouth opening larvae started to feed on rotifers (Brachionus sp.). Artemia nauplii were introduced at 5 days post-hatch (dph) and their density was gradually increased from 4 to 8 ml-1 until 8 dph. Live Artemia meta-nauplii were then introduced at this time point on increasing from 12 to 20 · ml-1 until 18 dph. Frozen meta-nauplii were offered from when larvae settled at the bottom of the tanks (around 18 dph) until 30 dph.
Sampling and RNA extraction
Embryo and larvae samples were taken in triplicate according to the developmental stages described in Fernández-Díaz et al. . The sampling points were: 75% epiboly (75Ep), 20 somites (20S), hatching (0 dph), pre-metamorphosis (Pre-Met, stage 0 according to the eye-translocation stage; 8 and 9 dph for larvae from 21°C and 15°C, respectively), metamorphosis climax (Met, stage 3 according to the eye-translocation stage; 14 dph and 15 dph in larvae incubated at 21°C and 15°C, respectively) and 30 dph. The embryos were snap-frozen in liquid nitrogen, and larvae were sacrificed by terminal anaesthesia with MS-222 (400 mg · l-1) and snap-frozen in liquid nitrogen. All samples were stored at -80°C until further analysis. Total RNA was extracted using the mirVana miRNA isolation kit (Ambion) according to the manufacturer’s instructions. Assessment of RNA quality was performed by agarose gel electrophoresis on a 1.2% (w/v) agarose gel containing SYBR safe DNA gel stain. RNA samples were then quantified with a Nanodrop spectrophotometer (Nanodrop Technologies/Saven Werner). Absorbance ratios (260/280 nm) were greater than 1.9, indicating high purity RNA.
SOLiD sequencing of Senegalese sole small RNAs
Small RNA reads were generated from eight SOLiD (Life Technologies, Austin, Texas, USA) libraries prepared from pools of triplicate samples of embryos and larvae from the following stages and at both 15 and 21°C temperature groups: 75Ep, 20S, hatching and metamorphosis climax stage 3. The enriched small RNA fraction of less than 40 nucleotides (nt) was recovered by sodium acetate precipitation, assessed on Agilent 2100 Bioanalyzer with small RNA chip (Agilent Technologies), and used for the amplified small RNA library template preparation using the SOLiD™ total RNA-seq kit with the barcoded SOLiD™ 3’ primers (Life Sciences). The average length of the amplified cDNA template was assessed in Agilent 2100 Bioanalyzer with the high sensitivity DNA chip kit (Agilent Technologies) to calculate the ratio of the desired miRNA ligation products (120–130 bp). To quantify the precise cDNA yield of miRNA ligation products for template bead preparation, real-time quantitative PCR (qPCR) was carried out using KAPA SYBR FAST Light Cycler 480 qPCR kit with SOLiD™ primer mix and DNA standards (Kapa Biosystems). Each library template was clonally amplified on SOLiD™ P1 beads by emulsion PCR and then sequenced on the SOLiD4 platform (Life Sciences) on a quadrant of a slide with 50 nt read length.
Sequencing data analysis
The SOLiD raw colour space reads were pre-processed using the CLC Genomics Workbench 4.9 (CLCbio) to remove low quality reads and clip internal adaptor sequence (CTGCTGTACGGCCAAGGCC), thus leaving reads ranging between 10 and 35 nucleotides (nt) in length. Reads corresponding to small RNAs other than miRNAs (such as tRNAS, rRNA, mtRNA or snRNAs) were discarded following BLAST searches against ncRNA databases. Unannotated clusters and singletons were mapped against miRNAs sequences deposited in miRBase 18.0 (http://www.mirbase.org/) in order to identify conserved miRNAs. The remaining reads were mapped against Pleurinectiformes and Solea EST databases (NCBI) to identify possible non-conserved miRNAs. For each mappable sequence between 16 and 28 nt and with 5 counts or above, hairpin folding was evaluated by sequence analysis to identify the presence of a stem loop using the RNA secondary structure prediction tool in the CLC Genomics Workbench. The flanking EST sequence of each analysed read should be at least 60 nt, and the free energy of the fold-back structure (ΔG) should be ≤ -18 kcal · mol-1.
The number of reads in each sample was normalised by dividing the number of reads of a given miRNA in a sample by the total number of trimmed miRNA reads in the same sample. To be able to differentiate the expression of each miRNA and categorise them according to their expression pattern, a heatmap chart was drawn and hierarchical clustering with Euclidean distance and complete linkage algorithms was performed using the MeV online software (http://www.tm4.org/mev.html).
Quantitative real-time PCR assay
List of primers used in quantification of miRNAs by qPCR
Mature sequence (5’ – 3’)
ID qPCR primers (Exiqon)
cDNA was synthesised from RNA diluted to 5 ng · μl-1 using the miRCURY LNA™ Universal RT microRNA PCR Kit (Exiqon) and following the manufacturer’s instructions. Reactions were performed using SYBR Green chemistry (Exiqon) on a Light Cycler 480® (Roche Applied Science) using white 96-well plates and the following running conditions: denaturation at 95°C: 10 min followed by 45 amplification cycles at 95°C: 10 s and 60°C: 1 min. Specificity of reaction was analysed by melting curve analysis . Five-point standard curves of a 5-fold dilution series (1:1–1:625) of pooled RNA were used for PCR efficiency calculation . All samples were run in duplicate. CT values were determined using the fit-point method using the LightCycler® 480 software with a fluorescence threshold arbitrarily set to 0.3. The provided control RNA Spike-in (Exiqon) was used as exogenous reference to examine the miRNAs expression profile during development and across temperatures, using the ΔΔCT method, according to Fernandes et al. . Differences in the expression of the target miRNAs between temperatures were evaluated by Student’s t-test using the SigmaPlot 11.0 statistical software (Systat Software, San Jose, CA, USA), since the normality assumption was met.
Target prediction analysis
3’-UTRs for 3,946 annotated ESTs representing 3,397 Senegalese sole genes were obtained from NCBI-SRA and dbEST. Raw data from Senegalese sole 454 sequencing have been submitted to SRA under the experiment ID SRX246914. MiRNA target gene prediction was performed using miRanda v3.3a  and the RNAhybrid web server (bibiserv.techfak.uni-bielefeld.de/rnahybrid) . MiRanda implements a dynamic search programme for maximal local complementarity alignments between the miRNA and the reference sequence and generate scores based on the position-specific empirical models . The free energy of optimal strand-strand interaction was calculated using the ViennaRNA package incorporated within the program. Targets for Senegalese sole miRNAs were predicted using the following default parameters: Gap open penalty, -9.0; gap extend penalty, -4.0; score threshold, 140.0; energy threshold, 1.0 kcal · mol-1; scaling parameter, 4.0. Targets predicted to have an energy threshold below -20 kcal · mol-1 were selected for further analysis. The RNAhybrid web server was further used to calculate the minimum free energy of the stable duplex between miRNA and target site within the 3’-UTR.
The authors acknowledge J. Sendão, M.F. Castanheira, P.Gavaia, E. Cabrita, H. Teixeira, F. Rocha and A. Santos (CCMAR, Portugal), for their invaluable help during the experimental setup and sampling of Senegalese sole embryos and larvae. This study was funded by Project EPISOLE - PTDC/MAR/110547/2009, granted by Fundação para a Ciência e a Tecnologia and Programa Operacional Temático Factores de Competitividade (COMPETE), FEDER (Portugal). J. Fernandes acknowledges the support provided by the Research Council of Norway (GrowCod project, ref. 190350/S40). C. Campos and S. Engrola were financially supported by Fundação para a Ciência e Tecnologia, Portugal, through grants SFRH/BD/43633/2008 and SFRH/BPD/49051/2008, respectively.
- Chang TC, Mendell JT: microRNAs in vertebrate physiology and human disease. Annu Rev Genomics Hum Genet. 2007, 8: 215-239.PubMedView ArticleGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297.PubMedView ArticleGoogle Scholar
- Wood AW, Duan C, Bern HA: Insulin-like growth factor signaling in fish. Int Rev Cytol. 2005, 243: 215-285.PubMedView ArticleGoogle Scholar
- Nave BT, Ouwens M, Withers DJ, Alessi DR, Shepherd PR: Mammalian target of rapamycin is a direct target for protein kinase B: identification of a convergence point for opposing effects of insulin and amino-acid deficiency on protein translation. Biochem J. 1999, 344 (Pt 2): 427-431.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones NC, Fedorov YV, Rosenthal RS, Olwin BB: ERK1/2 is required for myoblast proliferation but is dispensable for muscle gene expression and cell fusion. J Cell Physiol. 2001, 186 (1): 104-115.PubMedView ArticleGoogle Scholar
- Ren H, Accili D, Duan C: Hypoxia converts the myogenic action of insulin-like growth factors into mitogenic action by differentially regulating multiple signaling pathways. Proc Natl Acad Sci U S A. 2010, 107 (13): 5857-5862.PubMed CentralPubMedView ArticleGoogle Scholar
- Yecies JL, Zhang HH, Menon S, Liu S, Yecies D, Lipovsky AI, Gorgun C, Kwiatkowski DJ, Hotamisligil GS, Lee CH, Manning BD: Akt stimulates hepatic SREBP1c and lipogenesis through parallel mTORC1-dependent and independent pathways. Cell Metab. 2011, 14 (1): 21-32.PubMed CentralPubMedView ArticleGoogle Scholar
- Fornari F, Milazzo M, Chieco P, Negrini M, Calin GA, Grazi GL, Pollutri D, Croce CM, Bolondi L, Gramantieri L: MiR-199a-3p regulates mTOR and c-Met to influence the doxorubicin sensitivity of human hepatocarcinoma cells. Cancer Res. 2010, 70 (12): 5184-5193.PubMedView ArticleGoogle Scholar
- Merkel O, Hamacher F, Laimer D, Sifft E, Trajanoski Z, Scheideler M, Egger G, Hassler MR, Thallinger C, Schmatz A, Turner SD, Greil R, Kenner L: Identification of differential and functionally active miRNAs in both anaplastic lymphoma kinase (ALK) + and ALK- anaplastic large-cell lymphoma. Proc Natl Acad Sci U S A. 2010, 107 (37): 16228-16233.PubMed CentralPubMedView ArticleGoogle Scholar
- Williams A, Henao-Mejia J, Harman CC, Flavell RA: miR-181 and Metabolic Regulation in the Immune System. Cold Spring Harb Symp Quant Biol. 2013, doi:10.1101/sqb.2013.78.020024Google Scholar
- Naguibneva I, Ameyar-Zazoua M, Polesskaya A, Ait-Si-Ali S, Groisman R, Souidi M, Cuvellier S, Harel-Bellan A: The microRNA miR-181 targets the homeobox protein Hox-A11 during mammalian myoblast differentiation. Nat Cell Biol. 2006, 8 (3): 278-284.PubMedView ArticleGoogle Scholar
- Chen JF, Mandel EM, Thomson JM, Wu Q, Callis TE, Hammond SM, Conlon FL, Wang DZ: The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation. Nat Genet. 2006, 38 (2): 228-233.PubMed CentralPubMedView ArticleGoogle Scholar
- Johnston IA, Lee HT, Macqueen DJ, Paranthaman K, Kawashima C, Anwar A, Kinghorn JR, Dalmay T: Embryonic temperature affects muscle fibre recruitment in adult zebrafish: genome-wide changes in gene and microRNA expression associated with the transition from hyperplastic to hypertrophic growth phenotypes. J Exp Biol. 2009, 212 (Pt 12): 1781-1793.PubMedView ArticleGoogle Scholar
- Yan X, Ding L, Li Y, Zhang X, Liang Y, Sun X, Teng CB: Identification and profiling of microRNAs from skeletal muscle of the common carp. PLoS One. 2012, 7 (1): e30925-PubMed CentralPubMedView ArticleGoogle Scholar
- Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, Watts L, Booten SL, Graham M, McKay R, Subramaniam A, Propp S, Lollo BA, Freier S, Bennett CF, Bhanot S, Monia BP: miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006, 3 (2): 87-98.PubMedView ArticleGoogle Scholar
- Rayner KJ, Suarez Y, Davalos A, Parathath S, Fitzgerald ML, Tamehiro N, Fisher EA, Moore KJ, Fernandez-Hernando C: MiR-33 contributes to the regulation of cholesterol homeostasis. Science. 2010, 328 (5985): 1570-1573.PubMed CentralPubMedView ArticleGoogle Scholar
- Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, Tuschl T: Identification of tissue-specific microRNAs from mouse. Curr Biol. 2002, 12 (9): 735-739.PubMedView ArticleGoogle Scholar
- Mennigen JA, Panserat S, Larquier M, Plagnes-Juan E, Medale F, Seiliez I, Skiba-Cassy S: Postprandial regulation of hepatic microRNAs predicted to target the insulin pathway in rainbow trout. PLoS One. 2012, 7 (6): e38604-PubMed CentralPubMedView ArticleGoogle Scholar
- Bizuayehu TT, Lanes CF, Furmanek T, Karlsen BO, Fernandes JM, Johansen SD, Babiak I: Differential expression patterns of conserved miRNAs and isomiRs during Atlantic halibut development. BMC Genomics. 2012, 13: 11-PubMed CentralPubMedView ArticleGoogle Scholar
- Mennigen JA, Skiba-Cassy S, Panserat S: Ontogenesis of expression of metabolic genes and microRNAs in rainbow trout alevins during the transition from the endogenous to the exogenous feeding period. J Exp Biol. 2013, 216: 1597-1608.PubMedView ArticleGoogle Scholar
- Johnston IA: Environment and plasticity of myogenesis in teleost fish. J Exp Biol. 2006, 209 (Pt 12): 2249-2264.PubMedView ArticleGoogle Scholar
- Imsland AK, Foss A, Conceição LEC, Dinis MT, Delbare D, Schram E, Kamstra A, Rema P, White P: A review of the culture potential of Solea solea and S. senegalensis. Rev Fish Biol Fisher. 2003, 13 (4): 379-407.View ArticleGoogle Scholar
- Campos C, Valente LM, Conceicao LE, Engrola S, Sousa V, Rocha E, Fernandes JM: Incubation temperature induces changes in muscle cellularity and gene expression in Senegalese sole (Solea senegalensis). Gene. 2013, 516 (2): 209-217.PubMedView ArticleGoogle Scholar
- Levy C, Khaled M, Iliopoulos D, Janas MM, Schubert S, Pinner S, Chen PH, Li S, Fletcher AL, Yokoyama S, Scott KL, Garraway LA, Song JS, Granter SR, Turley SJ, Fisher DE, Novina CD: Intronic miR-211 assumes the tumor suppressive function of its host gene in melanoma. Mol Cell. 2010, 40 (5): 841-849.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhai Q, Zhou L, Zhao C, Wan J, Yu Z, Guo X, Qin J, Chen J, Lu R: Identification of miR-508-3p and miR-509-3p that are associated with cell invasion and migration and involved in the apoptosis of renal cell carcinoma. Biochem Biophys Res Commun. 2012, 419 (4): 621-626.PubMedView ArticleGoogle Scholar
- Boulias K, Horvitz HR: The C. elegans microRNA mir-71 acts in neurons to promote germline-mediated longevity through regulation of DAF-16/FOXO. Cell Metab. 2012, 15 (4): 439-450.PubMed CentralPubMedView ArticleGoogle Scholar
- Guo L, Lu Z: The fate of miRNA* strand through evolutionary analysis: implication for degradation as merely carrier strand or potential regulatory molecule?. PLoS One. 2010, 5 (6): e11387-PubMed CentralPubMedView ArticleGoogle Scholar
- Heisenberg CP, Tada M, Rauch GJ, Saude L, Concha ML, Geisler R, Stemple DL, Smith JC, Wilson SW: Silberblick/Wnt11 mediates convergent extension movements during zebrafish gastrulation. Nature. 2000, 405 (6782): 76-81.PubMedView ArticleGoogle Scholar
- Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF: Stages of embryonic development of the zebrafish. Dev Dyn. 1995, 203 (3): 253-310.PubMedView ArticleGoogle Scholar
- Stickney HL, Barresi MJ, Devoto SH: Somite development in zebrafish. Dev Dyn. 2000, 219 (3): 287-303.PubMedView ArticleGoogle Scholar
- Trikic MZ, Monk P, Roehl H, Partridge LJ: Regulation of zebrafish hatching by tetraspanin cd63. PLoS One. 2011, 6 (5): e19683-PubMed CentralPubMedView ArticleGoogle Scholar
- Funes V, Asensio E, Ponce M, Infante C, Canavate JP, Manchado M: Insulin-like growth factors I and II in the sole Solea senegalensis: cDNA cloning and quantitation of gene expression in tissues and during larval development. Gen Comp Endocrinol. 2006, 149 (2): 166-172.PubMedView ArticleGoogle Scholar
- Darnell DK, Kaur S, Stanislaw S, Konieczka JH, Yatskievych TA, Antin PB: MicroRNA expression during chick embryo development. Dev Dyn. 2006, 235 (11): 3156-3165.PubMedView ArticleGoogle Scholar
- Giraldez AJ, Mishima Y, Rihel J, Grocock RJ, Van Dongen S, Inoue K, Enright AJ, Schier AF: Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006, 312 (5770): 75-79.PubMedView ArticleGoogle Scholar
- Thatcher EJ, Flynt AS, Li N, Patton JR, Patton JG: MiRNA expression analysis during normal zebrafish development and following inhibition of the Hedgehog and Notch signaling pathways. Dev Dyn. 2007, 236 (8): 2172-2180.PubMedView ArticleGoogle Scholar
- Pasquinelli AE, Reinhart BJ, Slack F, Martindale MQ, Kuroda MI, Maller B, Hayward DC, Ball EE, Degnan B, Muller P, Spring J, Srinivasan A, Fishman M, Finnerty J, Corbo J, Levine M, Leahy P, Davidson E, Ruvkun G: Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA. Nature. 2000, 408 (6808): 86-89.PubMedView ArticleGoogle Scholar
- Xie CX, Xu SL, Yang LL, Ke ZH, Xing JB, Gai JW, Gong XL, Xu LX, Bao BL: mRNA/microRNA Profile at the Metamorphic Stage of Olive Flounder (Paralichthys olivaceus). Comp Funct Genomics. 2011, 2011: 256038-PubMed CentralPubMedView ArticleGoogle Scholar
- Fernández-Díaz C, Yúfera M, Cañavate JP, Moyano FJ, Alarcón FJ, Díaz M: Growth and physiological changes during metamorphosis of Senegal sole reared in the laboratory. J Fish Biol. 2001, 58 (4): 1086-1097.View ArticleGoogle Scholar
- Schreiber AM: Asymmetric craniofacial remodeling and lateralized behavior in larval flatfish. J Exp Biol. 2006, 209 (Pt 4): 610-621.PubMedView ArticleGoogle Scholar
- Fu YS, Shi ZY, Wang GY, Li WJ, Zhang JL, Jia L: Expression and regulation of miR-1, -133a, -206a, and MRFs by thyroid hormone during larval development in Paralichthys olivaceus. Comp Biochem Physiol B. 2012, 161 (3): 226-232.PubMedView ArticleGoogle Scholar
- Yufera M, Parra G, Santiago R, Carrascosa M: Growth, carbon, nitrogen and caloric content of Solea senegalensis (Pisces: Soleidae) from egg fertilization to metamorphosis. Marine Biol. 1999, 134: 43-49.View ArticleGoogle Scholar
- Fu Y, Shi Z, Wu M, Zhang J, Jia L, Chen X: Identification and differential expression of microRNAs during metamorphosis of the Japanese flounder (Paralichthys olivaceus). PLoS One. 2011, 6 (7): e22957-PubMed CentralPubMedView ArticleGoogle Scholar
- Mennigen JA, Martyniuk CJ, Seiliez I, Panserat S, Skiba-Cassy S: Metabolic consequences of microRNA-122 inhibition in rainbow trout: Oncorhynchus mykiss. BMC Genomics. 2014, 15: 70-PubMed CentralPubMedView ArticleGoogle Scholar
- Wong CF, Tellam RL: MicroRNA-26a targets the histone methyltransferase Enhancer of Zeste homolog 2 during myogenesis. J Biol Chem. 2008, 283 (15): 9836-9843.PubMedView ArticleGoogle Scholar
- Safdar A, Abadi A, Akhtar M, Hettinga BP, Tarnopolsky MA: miRNA in the regulation of skeletal muscle adaptation to acute endurance exercise in C57Bl/6 J male mice. PLoS One. 2009, 4 (5): e5610-PubMed CentralPubMedView ArticleGoogle Scholar
- Kim HK, Lee YS, Sivaprasad U, Malhotra A, Dutta A: Muscle-specific microRNA miR-206 promotes muscle differentiation. J Cell Biol. 2006, 174 (5): 677-687.PubMed CentralPubMedView ArticleGoogle Scholar
- Campos C, Valente LMP, Conceição LEC, Engrola S, Fernandes JMO: Temperature affects methylation of the myogenin putative promoter, its expression and muscle cellularity in Senegalese sole larvae. Epigenetics. 2013, 8 (4): 389-397.PubMed CentralPubMedView ArticleGoogle Scholar
- Liu N, Williams AH, Kim Y, McAnally J, Bezprozvannaya S, Sutherland LB, Richardson JA, Bassel-Duby R, Olson EN: An intragenic MEF2-dependent enhancer directs muscle-specific expression of microRNAs 1 and 133. Proc Natl Acad Sci U S A. 2007, 104 (52): 20844-20849.PubMed CentralPubMedView ArticleGoogle Scholar
- Rosenberg MI, Georges SA, Asawachaicharn A, Analau E, Tapscott SJ: MyoD inhibits Fstl1 and Utrn expression by inducing transcription of miR-206. J Cell Biol. 2006, 175 (1): 77-85.PubMed CentralPubMedView ArticleGoogle Scholar
- Ason B, Darnell DK, Wittbrodt B, Berezikov E, Kloosterman WP, Wittbrodt J, Antin PB, Plasterk RH: Differences in vertebrate microRNA expression. Proc Natl Acad Sci U S A. 2006, 103 (39): 14385-14389.PubMed CentralPubMedView ArticleGoogle Scholar
- Juan AH, Kumar RM, Marx JG, Young RA, Sartorelli V: Mir-214-dependent regulation of the polycomb protein Ezh2 in skeletal muscle and embryonic stem cells. Mol Cell. 2009, 36 (1): 61-74.PubMed CentralPubMedView ArticleGoogle Scholar
- Maciotta S, Meregalli M, Cassinelli L, Parolini D, Farini A, Fraro GD, Gandolfi F, Forcato M, Ferrari S, Gabellini D, Bicciato S, Cossu G, Torrente Y: Hmgb3 is regulated by microrna-206 during muscle regeneration. PLoS One. 2012, 7 (8): e43464-PubMed CentralPubMedView ArticleGoogle Scholar
- Pearson G, Robinson F, Beers Gibson T, Xu BE, Karandikar M, Berman K, Cobb MH: Mitogen-activated protein (MAP) kinase pathways: regulation and physiological functions. Endocr Rev. 2001, 22 (2): 153-183.PubMedGoogle Scholar
- Gingras AC, Raught B, Gygi SP, Niedzwiecka A, Miron M, Burley SK, Polakiewicz RD, Wyslouch-Cieszynska A, Aebersold R, Sonenberg N: Hierarchical phosphorylation of the translation inhibitor 4E-BP1. Genes Dev. 2001, 15 (21): 2852-2864.PubMed CentralPubMedGoogle Scholar
- Marcotrigiano J, Gingras AC, Sonenberg N, Burley SK: Cap-dependent translation initiation in eukaryotes is regulated by a molecular mimic of eIF4G. Mol Cell. 1999, 3 (6): 707-716.PubMedView ArticleGoogle Scholar
- Brugarolas J, Lei K, Hurley RL, Manning BD, Reiling JH, Hafen E, Witters LA, Ellisen LW, Kaelin WG: Regulation of mTOR function in response to hypoxia by REDD1 and the TSC1/TSC2 tumor suppressor complex. Genes Dev. 2004, 18 (23): 2893-2904.PubMed CentralPubMedView ArticleGoogle Scholar
- DeYoung MP, Horak P, Sofer A, Sgroi D, Ellisen LW: Hypoxia regulates TSC1/2-mTOR signaling and tumor suppression through REDD1-mediated 14-3-3 shuttling. Genes Dev. 2008, 22 (2): 239-251.PubMed CentralPubMedView ArticleGoogle Scholar
- Cloonan N, Brown MK, Steptoe AL, Wani S, Chan WL, Forrest AR, Kolle G, Gabrielli B, Grimmond SM: The miR-17-5p microRNA is a key regulator of the G1/S phase cell cycle transition. Genome Biol. 2008, 9 (8): R127-PubMed CentralPubMedView ArticleGoogle Scholar
- Sabapathy K, Hochedlinger K, Nam SY, Bauer A, Karin M, Wagner EF: Distinct roles for JNK1 and JNK2 in regulating JNK activity and c-Jun-dependent cell proliferation. Mol Cell. 2004, 15 (5): 713-725.PubMedView ArticleGoogle Scholar
- Budanov AV, Karin M: p53 target genes sestrin1 and sestrin2 connect genotoxic stress and mTOR signaling. Cell. 2008, 134 (3): 451-460.PubMed CentralPubMedView ArticleGoogle Scholar
- Chen CC, Jeon SM, Bhaskar PT, Nogueira V, Sundararajan D, Tonic I, Park Y, Hay N: FoxOs inhibit mTORC1 and activate Akt by inducing the expression of Sestrin3 and Rictor. Dev Cell. 2010, 18 (4): 592-604.PubMed CentralPubMedView ArticleGoogle Scholar
- Honda M, Masui F, Kanzawa N, Tsuchiya T, Toyo-oka T: Specific knockdown of m-calpain blocks myogenesis with cDNA deduced from the corresponding RNAi. Am J Physiol Cell Physiol. 2008, 294 (4): 957-965.View ArticleGoogle Scholar
- Campos C, Valente LM, Borges P, Bizuayehu T, Fernandes JM: Dietary lipid levels have a remarkable impact on the expression of growth-related genes in Senegalese sole (Solea senegalensis Kaup). J Exp Biol. 2010, 213 (2): 200-209.PubMedView ArticleGoogle Scholar
- Fernandes JM, Mackenzie MG, Wright PA, Steele SL, Suzuki Y, Kinghorn JR, Johnston IA: Myogenin in model pufferfish species: Comparative genomic analysis and thermal plasticity of expression during early development. Comp Biochem Physiol Part D Genomics Proteomics. 2006, 1 (1): 35-45.PubMedView ArticleGoogle Scholar
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5 (1): R1-PubMed CentralPubMedView ArticleGoogle Scholar
- Kruger J, Rehmsmeier M: RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006, 34 (Web Server issue): W451-454.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.