Thermal plasticity of the miRNA transcriptome during Senegalese sole development

Embryonic temperature has long-term effects on muscle circRNA expression and somatic growth in Nile tilapia

Embryonic temperature has a lasting impact on muscle phenotype in vertebrates, involving complex molecular mechanisms that encompass both protein-coding and non-coding genes. Circular RNAs (circRNAs) are a class of regulatory RNAs that play important roles in various biological processes, but the effect of variable thermal conditions on the circRNA transcriptome and its long-term impact on muscle growth plasticity remains largely unexplored. To fill this knowledge gap, we performed a transcriptomic analysis of circRNAs in fast muscle of Nile tilapia (Oreochromis niloticus) subjected to different embryonic temperatures (24°C, 28°C and 32°C) and then reared at a common temperature (28°C) for 4 months. Nile tilapia embryos exhibited faster development and subsequently higher long-term growth at 32°C compared to those reared at 28°C and 24°C. Next-generation sequencing data revealed a total of 5,141 unique circRNAs across all temperature groups, of which 1,604, 1,531, and 1,169 circRNAs were exclusively found in the 24°C, 28°C and 32°C groups, respectively. Among them, circNexn exhibited a 1.7-fold (log2) upregulation in the 24°C group and a 1.3-fold (log2) upregulation in the 32°C group when compared to the 28°C group. Conversely, circTTN and circTTN_b were downregulated in the 24°C groups compared to their 28°C and 32°C counterparts. Furthermore, these differentially expressed circRNAs were found to have multiple interactions with myomiRs, highlighting their potential as promising candidates for further investigation in the context of muscle growth plasticity. Taken together, our findings provide new insights into the molecular mechanisms that may underlie muscle growth plasticity in response to thermal variation in fish, with important implications in the context of climate change, fisheries and aquaculture.

Plasma cortisol and production of miRNAs in red drum (Sciaenops ocellatus) exposed to three distinct challenges

The red drum Sciaenops ocellatus is a marine fish species of high commercial interest. Despite improvements in current aquaculture practices, there are still concerns about the impact of daily manipulations regarding fish welfare. To investigate how does fish respond to various challenges, S. ocellatus juveniles were submitted to two acute challenges, namely a confinement stress and a cold-temperature shock, as well as a chronic stress challenge consisting of 18 days of repetitive challenge events. The level of cortisol produced by individuals was used as a measure of activation hypothalamic-pituitary-interrenal (HPI) axis. A significant increase in cortisol levels was detected only after the confinement stress. Interestingly, the fish exposed to a chronic stress for 18 days exhibited cortisol levels significantly lower than those of non-challenged fish. The small RNA-sequencing conducted for the chronic stress experiment only allowed us to identify two plasmatic microRNAs more abundant in non-challenged fish (miR-205-1-5p and let-7b-5p) compared to challenged fish. The miR-205-1-5p was of particular interest since it was already detected in previous studies on other fish species. In silico analysis allowed to predict potentially highly conserved mRNA targets of this specific miRNA, among which is tnfrsfa that plays a key role in the secondary stress response.

Identification of serum microRNA alterations associated with long-term exercise-induced motor improvements in patients with Parkinson disease

BACKGROUND

Long-term physical exercise has been shown to benefit patients with Parkinson disease (PD), but there is a lack of evidence regarding the underlying mechanism. A better understanding of how such benefits are induced by exercise might contribute to the development of therapeutic targets for improving the motor function in individuals with PD. The purpose of this study was therefore to investigate the possible association between exercise-induced motor improvements and the changes in serum microRNA (miRNA) levels of PD patients through small RNA sequencing for the first time.

METHODS

Thirteen PD patients completed our 3-month home-and-community-based exercise program, while 6 patients were assigned to the control group. Motor functions were measured, and small RNA sequencing with data analysis was performed on serum miRNAs both before and after the program. The results were further validated by quantitative real-time polymerase chain reaction. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were then conducted to determine the role of differentially expressed miRNAs.

RESULTS

The 3-month home-and-community-based exercise program induced significant motor improvements in PD patients in terms of Unified Parkinson's Disease Rating Scale activities of daily living and Motor Subscale (P < .05), comfortable walking speed (P = .003), fast walking speed (P = .028), Six-Minute Walk Test (P = .004), Berg Balance Scale (P = .039), and Timed Up and Go (P = .002). A total of 11 miRNAs (10 upregulated and one downregulated) were identified to be remarkably differentially expressed after intervention in the exercise group, but not in the control group. The results of miRNA sequencing were further validated by quantitative real-time polymerase chain reaction. It was found that the targets of altered miRNAs were mostly enriched in the mitogen-activated protein kinase, Wnt, and Hippo signaling pathways and the GO annotations mainly included binding, catalytic activity, and transcription regulator activity.

CONCLUSION

The exercise-induced motor improvements were possibly associated with changes in circulating miRNA levels in PD patients. These miRNAs, as well as the most enriched pathways and GO terms, may play a critical role in the mechanism of exercise-induced benefits in PD and serve as novel treatment targets for the disease, although further investigations are needed.

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

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

Are fish immunocompetent enough to face climate change?

Ongoing climate change has already been associated with increased disease outbreaks in wild and farmed fish. Here, we evaluate the current knowledge of climate change-related ecoimmunology in teleosts with a focus on temperature, hypoxia, salinity and acidification before exploring interactive effects of multiple stressors. Our literature review reveals that acute and chronic changes in temperature and dissolved oxygen can compromise fish immunity which can lead to increased disease susceptibility. Moreover, temperature and hypoxia have already been shown to enhance the infectivity of certain pathogens/parasites and to accelerate disease progression. Too few studies exist that have focussed on acidification, but direct immune effects seem to be limited while salinity studies have led to contrasting results. Likewise, multi-stressor experiments essential for unravelling the interactions of simultaneously changing environmental factors are still scarce. This ultimately impedes our ability to estimate to what extent climate change will hamper fish immunity. Our review about epigenetic regulation mechanisms highlights the acclimation potential of the fish immune response to changing environments. However, due to the limited number of epigenetic studies, overarching conclusions cannot be drawn. Finally, we provide an outlook on how to better estimate the effects of realistic climate change scenarios in future immune studies in fish.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

The novel circular RNA CircMef2c is positively associated with muscle growth in Nile tilapia

Muscle growth in teleosts is a complex biological process orchestrated by numerous protein-coding genes and non-coding RNAs. A few recent studies suggest that circRNAs are involved in teleost myogenesis, but the molecular networks involved remain poorly understood. In this study, an integrative omics approach was used to determine myogenic circRNAs in Nile tilapia by quantifying and comparing the expression profile of mRNAs, miRNAs, and circRNAs in fast muscle from full-sib fish with distinct growth rates. There were 1947 mRNAs, 9 miRNAs, and 4 circRNAs differentially expressed between fast- and slow-growing individuals. These miRNAs can regulate myogenic genes and have binding sites for the novel circRNA circMef2c. Our data indicate that circMef2c may interact with three miRNAs and 65 differentially expressed mRNAs to form multiple competing endogenous RNA networks that regulate growth, thus providing novel insights into the role of circRNAs in the regulation of muscle growth in teleosts.

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

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

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

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

miRNA and mRNA expression analysis reveals the effects of continuous heat stress on antibacterial responses to Aeromonas hydrophila lipopolysaccharide (LPS) in grass carp (Ctenopharyngodon idella)

Grass carp (Ctenopharyngodon idella) is the largest economic fish in freshwater culture in China, which is predisposed to infectious diseases under high temperature. Under the background of global warming, the industrialization of the Pearl River Delta region has led to aggravated thermal pollution, which has increasingly serious impacts on the aquatic ecological environment. This will result in more frequent exposure of grass carp to overheated water temperatures. Previous studies have only identified the regulatory genes of fish that respond to pathogens or temperature stress, but the transcriptional response to both is unknown. In this study, the histopathological analysis showed heat stress exacerbated spleen damage induced by Aeromonas hydrophila. The transcriptional responses of the spleens from A. hydrophila lipopolysaccharide (LPS) -injected grass carp undergoing heat stress and at normal temperatures for 6, 24, and 72 h were investigated by mRNA and microRNA sequencing. We identified 28, 20, and 141 differentially expressed (DE) miRNAs and 126, 383, and 4841 DE mRNAs between the two groups after 6, 24, and 72 h, respectively. There were 67 DE genes mainly involved in the cytochrome P450 pathway, antioxidant defense, inflammatory response, pathogen recognition pathway, antigen processing and presentation, and the ubiquitin-proteasome system. There were 5 DE miRNAs involved in regulating apoptosis and inflammation. We further verified 17 DE mRNAs and 5 DE miRNAs using quantitative real-time PCR. Based on miRNAs and mRNAs analysis, continuous heat stress will affect the antibacterial responses of grass carp spleens, resulting in aggravation of spleen injury. Together, these results provide data for further understanding of the decreased tolerance of fish to pathogen infection in persistent high-temperature environments.

Temperature-Biased miRNA Expression Patterns during European Sea Bass (Dicentrarchus labrax) Development

Environmental effects and, particularly, temperature changes have been demonstrated to influence the activity, function, and well-being of teleosts. Temperature may change seasonally in the wild, and in captivity under aquaculture operations. Moreover, climate change is expected to shift temperature profiles worldwide. MicroRNAs (miRNA) are important temperature-sensitive gene-expression regulators acting at the post-transcriptional level. They are known to be key regulators in development, reproduction, and immune responses. Therefore, early larval development of the European sea bass (Dicentrarchus labrax), one of the most extensively cultured species in Mediterranean aquaculture, was investigated at early rearing temperatures, i.e., 15, 17.5, and 20 °C, in regard to the impact of temperatures on miRNAs through sncRNA high-throughput sequencing but also at the phenotypic level in terms of growth, sex, vision, and skeletal deformities. Expression profiling revealed stage- and temperature-specific miRNA expression targeting genes with roles in reproduction and immune response mainly at the flexion and all-fins stages. Similar stage- and temperature-specific results were also observed concerning the number of rod cells and lower jaw elongation. The present work presents for the first time highly promising results on the influence of early rearing temperature at the post-transcriptional level during European sea bass development, with a putative impact on reproduction and immune response, as well as regarding teleost vision and larval development.

Biochemical indices, gene expression, and SNPs associated with salinity adaptation in juvenile chum salmon (Oncorhynchus keta) as determined by comparative transcriptome analysis

Chum salmon (Oncorhynchus keta) migrate from freshwater to saltwater, and incur developmental, physiological and molecular adaptations as the salinity changes. The molecular regulation for salinity adaptation in chum salmon is currently not well defined. In this study, 1-g salmon were cultured under 0 (control group, D0), 8‰ (D8), 16‰ (D16), and 24‰ (D24) salinity conditions for 42 days. Na⁺/K⁺-ATPase and Ca²⁺/Mg²⁺-ATPase activities in the gill first increased and then decreased in response to higher salinity environments where D8 exhibited the highest Na⁺/K⁺ATPase and Ca²⁺/Mg²⁺-ATPase activity and D24 exhibited the lowest. Alkaline phosphatase (AKP) activity was elevated in all salinity treatment groups relative to controls, while no significant difference in acid phosphatase (ACP) activity was observed across treatment groups. De novo transcriptome sequencing in the D0 and D24 groups using RNA-Seq analysis identified 187,836 unigenes, of which 2,143 were differentially expressed in response to environmental salinity (71 up-regulated and 2,072 down-regulated). A total of 56,020 putative single nucleotide polymorphisms (SNPs) were also identified. The growth, development, osmoregulation and maturation factors of N-methyl-D-aspartate receptors (nmdas) expressed in memory formation, as well as insulin-like growth factor 1 (igf-1) and igf-binding proteins (igfbps) were further investigated using targeted qRT-PCR. The lowest expression of all these genes occurred in the low salinity environments (D8 or D16), while their highest expression occurred in the high salinity environments (D24). These results provide preliminary insight into salinity adaptation in chum salmon and a foundation for the development of marker-assisted breeding for this species.

Insights Into miRNA-mRNA Regulatory Mechanisms of Cold Adaptation in Gymnocypris eckloni: Ubiquitin-Mediated Proteolysis Is Pivotal for Adaptive Energy Metabolism

This study aimed to understand cold stress adaptations mechanism in fish. Thus, the transcriptional response to cold conditions in Gymnocypris eckloni was evaluated using RNA-seq and microRNA (miRNA)-seq analyses. Low-temperature (LT) group G. eckloni was cultivated outdoors in waters cooled to 2–4°C for 3 weeks, while individuals in the control temperature (CT) group were exposed to 14–16°C. Significantly different responses were observed in both mRNA and miRNA expression profiles, with more mRNAs (1,833 and 1,869 mRNAs were up- and downregulated, respectively) and fewer miRNAs (15 and 6 were up- and downregulated, respectively) observed in the LT group individuals relative to the CT group individuals. A miRNA-mRNA network involved in the regulation of G. eckloni responses to cold stress was constructed; this network included ubiquitin-mediated proteolysis, protein processing, and oxidative phosphorylation. These results provided new insights into mechanisms of cold tolerance by fish, including decreased metabolic activity in addition to proteolysis.

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

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

Epigenetic mapping of the somatotropic axis in Nile tilapia reveals differential DNA hydroxymethylation marks associated with growth

In vertebrates, the somatotropic axis comprising the pituitary gland, liver and muscle plays a major role in myogenesis. Its output in terms of muscle growth is highly affected by nutritional and environmental cues, and thus likely epigenetically regulated. Hydroxymethylation is emerging as a DNA modification that modulates gene expression but a holistic characterization of the hydroxymethylome of the somatotropic axis has not been investigated to date. Using reduced representation 5-hydroxymethylcytosine profiling we demonstrate tissue-specific localization of 5-hydroxymethylcytosines at single nucleotide resolution. Their abundance within gene bodies and promoters of several growth-related genes supports their pertinent role in gene regulation. We propose that cytosine hydroxymethylation may contribute to the phenotypic plasticity of growth through epigenetic regulation of the somatotropic axis.

Molecular Mechanisms Regulating Muscle Plasticity in Fish

Growth rates in fish are largely dependent on genetic and environmental factors, of which the latter can be highly variable throughout development. For this reason, muscle growth in fish is particularly dynamic as muscle structure and function can be altered by environmental conditions, a concept referred to as muscle plasticity. Myogenic regulatory factors (MRFs) like Myogenin, MyoD, and Pax7 control the myogenic mechanisms regulating quiescent muscle cell maintenance, proliferation, and differentiation, critical processes central for muscle plasticity. This review focuses on recent advancements in molecular mechanisms involving microRNAs (miRNAs) and DNA methylation that regulate the expression and activity of MRFs in fish. Findings provide overwhelming support that these mechanisms are significant regulators of muscle plasticity, particularly in response to environmental factors like temperature and nutritional challenges. Genetic variation in DNA methylation and miRNA expression also correlate with variation in body weight and growth, suggesting that genetic markers related to these mechanisms may be useful for genomic selection strategies. Collectively, this knowledge improves the understanding of mechanisms regulating muscle plasticity and can contribute to the development of husbandry and breeding strategies that improve growth performance and the ability of the fish to respond to environmental challenges.

Effect of increased embryonic temperature during developmental windows on survival, morphology and oxygen consumption of rainbow trout (Oncorhynchus mykiss)

Temperature is a crucial environmental factor that influences physiological functions in fishes, and increased temperature during development can shape an organism's phenotype. An active line of inquiry in comparative developmental physiology is whether short-term exposure to thermal changes have lasting phenotypic effects. This is the first study to apply a developmental 3-dimensional critical window experimental design for a vertebrate, using time, temperature, and phenotypic response (i.e., variables measured). Rainbow trout (Oncorhynchus mykiss) are an anadromous species for which resident populations occupy freshwater environments that are likely impacted by variable and rising temperatures, particularly during embryonic development. To assess thermal effects on fish development, we examined trout hatchling phenotypes following rearing in constant temperatures (5, 10, 15 and 17.5 °C) and following exposure to increased temperature above 5 °C during specific developmental windows. Time to 50% hatch, hatchling mass and body length showed general trends of decreasing with increasing constant temperature, and survival was highest in constant 10 °C incubation. Thermally shifting embryos into 17.5 °C during gastrulation and organogenesis reduced survival at hatch compared to 10 °C, and exposure to 15 and 17.5 °C only late in development produced lighter and shorter hatchlings. Oxygen consumption rate (V̇o₂) at organogenesis differed between embryos incubated constantly in increased temperature or exposed only during organogenesis, but generally we found limited temperature effects on V̇o₂ that may be due to high data variability. Collectively, these results suggest that survival of rainbow trout hatchlings is most sensitive to 17.5 °C exposure during gastrulation and organogenesis, while warm water exposure later in development has greater impacts on morphology. Thus, trait-specific critical windows of thermal sensitivity exist for rainbow trout embryos that alter the hatching phenotype.

The genomic structure of the highly-conserved dmrt1 gene in Solea senegalensis (Kaup, 1868) shows an unexpected intragenic duplication

Knowing the factors responsible for sex determination in a species has significant theoretical and practical implications; the dmrt1 gene (Doublesex and Mab-3 (DM)-related Transcription factor 1) plays this role in diverse animal species. Solea senegalensis is a commercially important flat fish in which females grow 30% faster than males. It has 2n = 42 chromosomes and an XX / XY chromosome system for sex determination, without heteromorph chromosomes but with sex proto-chromosome. In the present study, we are providing the genomic structure and nucleotide sequence of dmrt1 gene obtained from cDNA from male and female adult gonads. A cDNA of 2027 containing an open-reading frame (ORF) of 1206 bp and encoding a 402 aa protein it is described for dmrt1 gene of S. senegalensis. Multiple mRNA isoforms indicating a high variable system of alternative splicing in the expression of dmrt1 of the sole in gonads were studied. None isoforms could be related to sex of individuals. The genomic structure of the dmrt1 of S. senegalensis showed a gene of 31400 bp composed of 7 exons and 6 introns. It contains an unexpected duplication of more than 10399 bp, involving part of the exon I, exons II and III and a SINE element found in the sequence that it is proposed as responsible for the duplication. A mature miRNA of 21 bp in length was localized at 336 bp from exon V. Protein-protein interacting networks of the dmrt1 gene showed matches with dmrt1 protein from Cynoglossus semilaevis and a protein interaction network with 11 nodes (dmrt1 plus 10 other proteins). The phylogenetic relationship of the dmrt1 gene in S. senegalensis is consistent with the evolutionary position of its species. The molecular characterization of this gene will enhance its functional analysis and the understanding of sex differentiation in Solea senegalensis and other flatfish.

Embryonic canalization and its limits-A view from temperature

Many animals are able to produce similar offspring over a range of environmental conditions. This property of the developmental process has been termed canalization-the channeling of developmental pathways to generate a stable outcome despite varying conditions. Temperature is one environmental parameter that has fundamental effects on cell physiology and biochemistry, yet developmental programs generally result in a stable phenotype under a range of temperatures. On the other hand, there are typically upper and lower temperature limits beyond which the developmental program is unable to produce normal offspring. This review summarizes data on how development is affected by temperature, particularly high temperature, in various animal species. It also brings together information on potential cell biological and developmental genetic factors that may be responsible for developmental stability in varying temperatures, and likely critical mechanisms that break down at high temperature. Also reviewed are possible means for studying temperature effects on embryogenesis and how to determine which factors are most critical at the high-temperature limits for normal development. Increased knowledge of these critical factors will point to the targets of selection under climate change, and more generally, how developmental robustness in varying environments is maintained.

Analyses of mRNA-seq and miRNA-seq of the brain reveal the sex differences of gene expression and regulation before and during gonadal differentiation in 17β-estradiol or 17α-methyltestosterone-induced olive flounder (Paralichthys olivaceus)

Olive flounder (Paralichthys olivaceus) is a commercially important flatfish species cultured in East Asia. Female flounders generally grow more rapidly than males, therefore control of the sex ratio seems to be a proposed way to increase production. However, the sex determination gene and sex determination mechanism have yet been elucidated. The brain is an important organ that is involved in gonadal development. To explore the sex differences of gene expression in the brain before and during the flounder gonadal differentiation, we used messenger RNA (mRNA)-seq technology to investigate transcriptomes of male and female brains. Between female and male brains, 103 genes were differentially expressed before ovarian differentiation, 16 genes were differentially expressed before testicular differentiation, and 64 genes were differentially expressed during gonadal differentiation. According to annotation and Kyoto Encyclopedia of Genes and Genomes information, the differentially expressed genes (DEGs) were involved in circadian rhythm, circadian rhythm-fly, circadian entrainment, dopaminergic synapse, calcium signaling, glutamatergic synapse, taste transduction, herpes simplex infection, long-term depression, retrograde endocannabinoid signaling, and the synaptic vesicle cycle pathways. MicroRNA (miRNA)-seq was performed during the gonadal differentiation and the target genes of miRNAs were predicted. Integrated analysis of mRNA-seq and miRNA-seq showed that 29 of the 64 DEGs were regulated by the differentially expressed miRNAs during the gonadal differentiation. Our study provides a basis for further studies of brain sex differentiation and the molecular mechanism of sex determination in olive flounder.

Chronic exposure to venlafaxine and increased water temperature reversibly alters microRNA in zebrafish gonads (Danio rerio)

MicroRNA (miRNA) are short, non-coding RNA that act by downregulating targeted mRNA transcripts. Only recently have they been used as endpoints in studies of aquatic toxicology. The purpose of this study was to determine the effect of an antidepressant contaminant, venlafaxine (VFX), and increased temperature on specific microRNA levels in zebrafish (Danio rerio) reproductive tissue. Adult zebrafish were exposed to one of four conditions; control, 1 μg/L VFX (VFX), 32 °C (Temp), or 1 μg/L VFX + 32 °C (VFX & Temp) for 21 days. Half of the fish were returned to control conditions for a 21-day recovery period. RT-qPCR was performed to measure relative abundances of several miRNAs known to respond to antidepressant exposure: dre-miR-22b-3p, dre-miR-301a, dre-miR-140-5p, dre-let-7d-5p, dre-miR-210-5p, and dre-miR-457b-5p. After the exposure period, dre-miR-22b-3p and dre-miR-301a showed a significant downregulation in response to all treatments. In contrast, after the recovery period, there were no significant differences in microRNA abundance. These altered microRNA are predicted to target several genes, including phosphofructokinase, and are associated with ovarian pathologies. Combined, we have shown that VFX and increased water temperature alter miRNA abundances in zebrafish reproductive tissue, an effect correlated with a functional stress response and cell cycle dysregulation.

Atrazine induces necroptosis by miR-181-5p targeting inflammation and glycometabolism in carp lymphocytes

Atrazine (ATR) causes environmental problems and damages the health of fish and aquatic animals. MicroRNAs (miRNAs) play important roles in immune regulation. However, the immunotoxicity mechanism of ATR in fish lymphocytes and the role of miRNA in this process remain unclear. To further study these mechanisms, spleen lymphocytes were exposed to 20, 40 and 60 μg/ml ATR for 18 h. Fluorescence staining and flow cytometry showed that the number of necrotic lymphocytes increased after ATR exposure. Compared with the control group, the mRNA expression of miR-181-5p was inhibited and the mRNA levels of TNF-α and HK2 were increased after ATR exposure. Additionally, the NF-κB inflammatory pathway and the levels of glycometabolism-related genes were upregulated. These results suggest that ATR induces inflammation and elevates glycometabolism in lymphocytes. We further found that the mRNA levels of receptor-interacting serine-threonine kinase 1 (RIP1), receptor-interacting serine-threonine kinase 3 (RIP3), mixed lineage kinase domain-like pseudokinase (MLKL), cylindromatosis (CYLD) and Fas-Associated protein with Death Domain (FADD) and the protein levels of RIP3 and MLKL in the treatment groups were significantly increased compared to those in control group, suggesting that ATR causes lymphocyte necroptosis. We conclude that miR-181-5p plays a key role in necroptosis in carp lymphocytes exposed to ATR by downregulating the expression of HK and TNF-α, which increases the level of glycometabolism and induces the inflammatory response, respectively.

Inventory of European Sea Bass (Dicentrarchus labrax) sncRNAs Vital During Early Teleost Development

During early animal ontogenesis, a plethora of small non-coding RNAs (sncRNAs) are greatly expressed and have been shown to be involved in several regulatory pathways vital to proper development. The rapid advancements in sequencing and computing methodologies in the last decade have paved the way for the production of sequencing data in a broad range of organisms, including teleost species. Consequently, this has led to the discovery of sncRNAs as well as the potentially novel roles of sncRNA in gene regulation. Among the several classes of sncRNAs, microRNAs (miRNAs) have, in particular, been shown to play a key role in development. The present work aims to identify the miRNAs that play important roles during early European sea bass (Dicentrarchus labrax) development. The European sea bass is a species of high commercial impact in European and especially Mediterranean aquaculture. This study reports, for the first time, the identification and characterization of small RNAs that play a part in the 10 developmental stages (from morula to all fins) of the European sea bass. From 10 developmental stages, more than 135 million reads, generated by next-generation sequencing, were retrieved from publicly available databases as well as newly generated. The analysis resulted in about 2,000 sample grouped reads, and their subsequently annotation revealed that the majority of transcripts belonged to the class of miRNAs followed by small nuclear RNAs and small nucleolar RNAs. The analysis of small RNA expression among the developmental stages under study revealed that miRNAs are active throughout development, with the main activity occurring after the earlier stages (morula and 50% epiboly) and at the later stages (first feeding, flexion, and all fins). Furthermore, investigating miRNAs exclusively expressed in one of the stages unraveled five miRNAs with a higher abundance only in the morula stage (miR-155, miR-430a, d1, d2, and miR-458), indicating possible important key roles of those miRNAs in further embryonic development. An additional target search showed putative miRNA-mRNA interactions with possible direct and indirect regulatory functions of the identified miRNAs.

High-throughput sequencing reveals microRNAs in response to heat stress in the head kidney of rainbow trout (Oncorhynchus mykiss)

Recently, the research of animal microRNAs (miRNAs) has attracted wide attention for its regulatory effect in the development process and the response to abiotic stresses. Rainbow trout is a commercially and cold water fish species, and usually encounters heat stress, which affects its growth and leads to a huge economic loss. But there were few investigations about the roles of miRNAs in heat stress in rainbow trout. In this study, miRNAs of rainbow trout which were involved in heat stress were identified by high-throughput sequencing of six small RNA libraries from head kidney tissues under control (18 °C) and heat-treated (24 °C) conditions. A total of 392 conserved miRNAs and 989 novel miRNAs were identified, of which 78 miRNAs were expressed in different response to heat stress. Ten of these miRNAs were further validated by quantitative real-time PCR. In addition to, including 393 negative correlation miRNA-target gene pairs, several important regulatory pathways were involved in heat stress of the potential target genes, including protein processing in endoplasmic reticulum, NOD-like receptor signaling pathway, and phagosome. Our data significantly advance understanding of heat stress regulatory mechanism of miRNA in the head kidney of rainbow trout, which provide a useful resource for the cultivation of rainbow trout.

Network of microRNA-transcriptional factor-mRNA in cold response of turbot Scophthalmus maximus

The aim of this study is to understand fish cold-tolerant mechanism. We analyzed the transcriptional reactions to the cold condition in turbot Scophthalmus maximus by using RNA-seq and microRNA (miRNA)-seq. Meio-gynogenetic diploid turbots were treated at 0 °C to distinguish the cold-tolerant (CT) and cold-sensitive (CS) groups. The results showed that there were quite different responses at both mRNA and miRNA levels, with more up-regulated mRNAs (1069 vs. 194) and less down-regulated miRNAs (4 vs. 1) in CT versus CS relative to the control group. The network of miRNA-transcription factor-mRNA, regulating turbot different response to cold stress, was constructed, which involved in cell cycle, component of cell membrane, signal transduction, and circadian rhythm pathways. The above information demonstrates mechanisms by which cold tolerance is increased in fish.

Characterization and expression profiling of microRNAs in response to plant feeding in two host-plant strains of the lepidopteran pest Spodoptera frugiperda

Background

A change in the environment may impair development or survival of living organisms leading them to adapt to the change. The resulting adaptation trait may reverse, or become fixed in the population leading to evolution of species. Deciphering the molecular basis of adaptive traits can thus give evolutionary clues. In phytophagous insects, a change in host-plant range can lead to emergence of new species. Among them, Spodoptera frugiperda is a major agricultural lepidopteran pest consisting of two host-plant strains having diverged 3 MA, based on mitochondrial markers. In this paper, we address the role of microRNAs, important gene expression regulators, in response to host-plant change and in adaptive evolution.

Results

Using small RNA sequencing, we characterized miRNA repertoires of the corn (C) and rice (R) strains of S. frugiperda, expressed during larval development on two different host-plants, corn and rice, in the frame of reciprocal transplant experiments. We provide evidence for 76 and 68 known miRNAs in C and R strains and 139 and 171 novel miRNAs. Based on read counts analysis, 34 of the microRNAs were differentially expressed in the C strain larvae fed on rice as compared to the C strain larvae fed on corn. Twenty one were differentially expressed on rice compared to corn in R strain. Nine were differentially expressed in the R strain compared to C strain when reared on corn. A similar ratio of microRNAs was differentially expressed between strains on rice. We could validate experimentally by QPCR, variation in expression of the most differentially expressed candidates. We used bioinformatics methods to determine the target mRNAs of known microRNAs. Comparison with the mRNA expression profile during similar reciprocal transplant experiment revealed potential mRNA targets of these host-plant regulated miRNAs.

Conclusions

In the current study, we performed the first systematic analysis of miRNAs in Lepidopteran pests feeding on host-plants. We identified a set of the differentially expressed miRNAs that respond to the plant diet, or differ constitutively between the two host plant strains. Among the latter, the ones that are also deregulated in response to host-plant are molecular candidates underlying a complex adaptive trait.

Electronic supplementary material

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

miRNAs regulate acute transcriptional changes in broiler embryos in response to modification of incubation temperature

MicroRNAs are post-transcriptional regulators that play critical roles in diverse biological processes. We hypothesize that miRNAs may be involved in regulating transcriptome responses to changes in embryonic incubation temperature in chickens affecting differentiation and proliferation processes during tissue development. Therefore, we conducted comparative transcriptome profiling of miRNAs to examine altered expression in breast and hind muscle of embryos and day 35 chickens experiencing high (38.8 °C), control (37.8 °C), or low (36.8 °C) embryonic incubation temperature during embryonic day (ED) 7–10 or ED10–13. The results revealed differential expression of miRNAs due to modification of embryonic incubation temperature in a muscle type-specific and a developmental stage-specific manner. The immediate effects of thermal change observed in embryos were substantial compared to the subtle long-term effects in chickens at day 35 post-hatch. Upregulation of miR-133 in breast muscle and downregulation of miR-199a-5p, miR-1915, and miR-638 in hind muscle post ED7–10 high-temperature treatment are functionally associated with myogenesis and body size. ED10–13 low-temperature treatment led to downregulation of let-7, miR-93, and miR-130c that are related to proliferation and differentiation. The results provide insight into the dynamics of miRNA expression at variable embryonic incubation temperatures during developmental processes and indicate a major regulatory role of miRNAs in acute responses to modified environmental conditions that affect remodelling of cells and tissues.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.

Small ocean temperature increases elicit stage-dependent changes in DNA methylation and gene expression in a fish, the European sea bass

In natural fish populations, temperature increases can result in shifts in important phenotypic traits. DNA methylation is an epigenetic mechanism mediating phenotypic changes. However, whether temperature increases of the magnitude predicted by the latest global warming models can affect DNA methylation is unknown. Here, we exposed European sea bass to moderate temperature increases in different periods within the first two months of age. We show that increases of even 2 °C in larvae significantly changed global DNA methylation and the expression of ecologically-relevant genes related to DNA methylation, stress response, muscle and organ formation, while 4 °C had no effect on juveniles. Furthermore, DNA methylation changes were more marked in larvae previously acclimated to a different temperature. The expression of most genes was also affected by temperature in the larvae but not in juveniles. In conclusion, this work constitutes the first study of DNA methylation in fish showing that temperature increases of the magnitude predicted by the latest global warming models result in stage-dependent alterations in global DNA methylation and gene expression levels. This study, therefore, provides insights on the possible consequences of climate change in fish mediated by genome-wide epigenetic modifications.

Phylogeny, expression patterns and regulation of DNA Methyltransferases in early development of the flatfish, Solea senegalensis

Background

The identification of DNA methyltransferases (Dnmt) expression patterns during development and their regulation is important to understand the epigenetic mechanisms that modulate larval plasticity in marine fish. In this study, dnmt1 and dnmt3 paralogs were identified in the flatfish Solea senegalensis and expression patterns in early developmental stages and juveniles were determined. Additionally, the regulation of Dnmt transcription by a specific inhibitor (5-aza-2′-deoxycytidine) and temperature was evaluated.

Results

Five paralog genes of dnmt3, namely dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1 and dnmt3bb.2 and one gene for dnmt1 were identified. Phylogenetic analysis revealed that the dnmt gene family was highly conserved in teleosts and three fish-specific genes, dnmt3aa, dnmt3ba and dnmt3bb.2 have evolved. The spatio-temporal expression patterns of four dnmts (dnmt1, dnmt3aa, dnmt3ab and dnmt3bb.1) were different in early larval stages although all of them reduced expression with the age and were detected in neural organs and dnmt3aa appeared specific to somites. In juveniles, the four dnmt genes were expressed in brain and hematopoietic tissues such as kidney, spleen and gills. Treatment of sole embryos with 5-aza-2′-deoxycytidine down-regulated dntm1 and up-regulated dntm3aa. Moreover, in lecithotrophic larval stages, dnmt3aa and dnmt3ab were temperature sensitive and their expression was higher in larvae incubated at 16 °C relative to 20 °C.

Conclusion

Five dnmt3 and one dnmt1 paralog were identified in sole and their distinct developmental and tissue-specific expression patterns indicate that they may have different roles during development. The inhibitor 5-aza-2′-deoxycytidine modified the transcript abundance of dntm1 and dntm3aa in embryos, which suggests that a regulatory feedback mechanism exists for these genes. The impact of thermal regime on expression levels of dnmt3aa and dnmt3ab in lecithotrophic larval stages suggests that these paralogs might be involved in thermal programing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-017-0154-0) contains supplementary material, which is available to authorized users.

Thermal imprinting modifies bone homeostasis in cold challenged sea bream (Sparus aurata, L.)

Fish are ectotherms and temperature plays a determinant role in their physiology, biology and ecology and is a driver of seasonal responses. The present study assessed how thermal imprinting during embryonic and larval stages modified the response of adult fish to low water temperature. We targeted the gilthead sea bream that develops a condition known as winter syndrome when it is exposed to low water temperatures. Eggs and larvae of sea bream were exposed to four different thermal regimes and then the response of the resulting adults to a low temperature challenge was assessed. Sea bream exposed to a high-low thermal regime as eggs and larvae (HLT, 22°C until hatch and then 18°C until larvae-juvenile transition) had increased plasma cortisol and lower sodium and potassium in response to a cold challenge compared to the other thermal history groups. Plasma glucose and osmolality were increased in cold challenge HLT fish relative to the unchallenged HLT fish. Cold challenge modified bone homeostasis/responsiveness in the low-high thermal regime group (LHT) relative to other groups and ocn, ogn1/2, igf1, gr and trα/β transcripts were all down-regulated. In the low temperature group (LT) and HLT group challenged with a low temperature, ALP/TRAP activities were decreased relative to unchallenged groups and bone calcium content also decreased in the LT group. Overall, the results indicate that thermal imprinting during early development of sea bream causes a change in the physiological response of adults to a cold challenge.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Resuscitation of a dead cardiomyocyte

Although cardiac resuscitation can revive the whole body, the mechanisms are unclear. To this end, we propose that reviving a dead/dysfunctional cardiomyocyte will shed light on resuscitation mechanisms and pave the way to treat cardiac myopathies. The degradation of the myocyte cytoskeleton by the proteasome system which involves calpains, ubiquitin, caspases and matrix metalloproteases is the main focus of this review. The activation of calpains beyond the calpastatin-mediated inhibition due to extensive calcium harbor can lead to titin degradation, damage to the sarcomere and contractile dysfunction. The ubiquitin proteasome system can disturb the protein homeostasis within the cell and generate a dysfunctional myocyte. The matrix metalloproteases disrupt the collagen/elastin ratio and connexins to generate arrhythmias. The concept of cardiac resuscitation stems from protecting the myocyte cytoskeleton and keeping the protein homeostasis intact through management of the degradation machinery. In this regard, proteasome inhibitors for the degradation machinery have an elegant space. Recently exosomes have been identified potentially, as carriers of microRNAs or proteins that can modify the target cells. Exosomes loaded with the inhibitor "cargo" which comprises microRNAs, siRNAs or proteins to inhibit the degradation machinery can be a method of choice for cardiac resuscitation-a process difficult to execute.

Micromanaging metabolism-a role for miRNAs in teleost energy metabolism

MicroRNAs (miRNAs) are small, non-protein coding RNA sequences, which are found in most eukaryotes. Since their initial discovery, miRNAs have emerged as important regulators of many biological processes. One of the most important processes profoundly regulated by miRNAs is energy metabolism. Traditionally, metabolic functions of miRNAs have been studied in genome-sequenced mammalian organisms, especially the mouse model. However, partially driven by commercial interest in aquaculture, increasingly feasible large-scale molecular techniques have resulted in the characterization of miRNA repertoires, and importantly, several genome sequences of several (commercially important) teleost species, which also hold important roles as research models in the comparative physiology of energy metabolism. This review aims to introduce the recent advances in miRNA research in teleost fish and to describe the current knowledge of miRNA function in teleost energy metabolism. The most pressing research needs and questions to determine metabolic roles of miRNAs in teleost models are presented, as well as applicable technical approaches and current bottlenecks. Rainbow trout, which possess the advantages of newly available molecular tools and a long history as comparative research model in teleost energy metabolism, are discussed as a promising research model to address these questions.

Temperature during early development has long-term effects on microRNA expression in Atlantic cod

Background

Environmental temperature has serious implications in life cycle of aquatic ectotherms. Understanding the molecular mechanisms of temperature acclimation and adaptation of marine organisms is of the uttermost importance for ecology, fisheries, and aquaculture, as it allows modeling the effects of global warming on population dynamics. Regulatory molecules are major modulators of acclimation and adaptation; among them, microRNAs (miRNAs) are versatile and substantial contributors to regulatory networks of development and adaptive plasticity. However, their role in thermal plasticity is poorly known. We have asked whether the temperature and its shift during the early ontogeny (embryonic and larval development) affect the miRNA repertoire of Atlantic cod (Gadus morhua), and if thermal experience has long-term consequences in the miRNA profile.

Results

We characterized miRNA during different developmental stages and in juvenile tissues using next generation sequencing. We identified 389 putative miRNA precursor loci, 120 novel precursor miRNAs, and 281 mature miRNAs. Some miRNAs showed stage- or tissue-enriched expression and miRNAs, such as the miR-17 ~ 92 cluster, myomiRs (miR-206), neuromiRs (miR-9, miR-124), miR-130b, and miR-430 showed differential expression in different temperature regimes. Long-term effect of embryonic incubation temperature was revealed on expression of some miRNAs in juvenile pituitary (miR-449), gonad (miR-27c, miR-30c, and miR-200a), and liver (let-7 h, miR-7a, miR-22, miR-34c, miR-132a, miR-192, miR-221, miR-451, miR-2188, and miR-7550), but not in brain. Some of differentially expressed miRNAs in the liver were confirmed using LNA-based rt-qPCR. The effect of temperature on methylation status of selected miRNA promoter regions was mostly inconclusive.

Conclusions

Temperature elevation by several degrees during embryonic and larval developmental stages significantly alters the miRNA profile, both short-term and long-term. Our results suggest that a further rise in seas temperature might affect life history of Atlantic cod.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1503-7) contains supplementary material, which is available to authorized users.