Deep RNA sequencing of the skeletal muscle transcriptome in swimming fish

Short-term induced swimming activity enhanced innate immune parameters and antioxidant status of European eel (Anguilla anguilla)

The swimming activity, although an essential trait in the life cycle of fish, is still poorly understood in farmed fish. The current study aimed to investigate the impact of short-term induced swimming on the immune and antioxidant defence systems in European eel (Anguilla anguilla). Sixteen male yellow European eels (total length: 39.9 ± 0.7 cm; body weight: 108.8 ± 6.1 g) were individually placed in swimming flumes and divided into two groups: i) no swimming (n = 8); and ii) induced-swimming (n = 8) at 0.3 body lengths (BL)·s-1 for 7 h. Swimming resulted in a 2-fold lower cortisol concentration in plasma, whereas plasma glucose, lactate, and several immune-related parameters did not present variations between groups. Interestingly, swimming led to higher lysozyme, peroxidase, and protease activities in skin mucus, whereas bactericidal activity did not show differences among groups. Additionally, the gene expression of interleukin 1 beta showed an up-regulation in the skin of fish with induced swimming, while no differences were observed in the head-kidney or gills. Furthermore, modulation of the antioxidant status was observed in the liver and posterior skeletal muscle after induced swimming. Fish subjected to swimming showed lower lipid peroxidation and higher reduced glutathione levels, increasing the reduced/oxidized glutathione ratio. However, no variations in the antioxidant status were observed between groups in the anterior skeletal muscle. This study showed modulation of immune and oxidative stress markers in European eels upon short-term induced swimming compared to non-swimming fish.

Comparative transcriptomic analyses of brain-liver-muscle in channel catfish (Ictalurus punctatus) with differential growth rate

Growth is an important economically trait for aquatic animals. The popularity of farmed channel catfish (Ictalurus punctatus) in China has recently surged, prompting a need for research into the genetic mechanisms that drive growth and development to expedite the selection of fast-growing variants. In this study, the brain, liver and muscle transcriptomes of channel catfish between fast-growing and slow-growing groups were analyzed using RNA-Seq. Totally, 63, 110 and 86 differentially expressed genes (DEGs) were from brain, liver and muscle tissues. DEGs are primarily involved in growth, development, metabolism and immunity, which are related to the growth regulation of channel catfish, such as growth hormone receptor b (ghrb), fibroblast growth factor receptor 4 (fgfr4), bone morphogenetic protein 1a (bmp1a), insulin-like growth factor 2a (igf2a), collagen, type I, alpha 1a (col1a1a), acyl-CoA synthetase long chain family member 2 (acsl2) and caveolin 1 (cav1). This study advances our knowledge of the genetic mechanisms accounting for differences in growth rate and offers crucial gene resources for future growth-related molecular breeding programs in channel catfish.

Transcriptome analysis and gene expression analysis related to salinity-alkalinity and low temperature adaptation of Triplophysa yarkandensis

T. yarkandensis is a common species of Triplophysa, and it is distributed in Shule river of Hexi Corridor, of Gansu province in China. In order to enrich gene database resources and explore the environment adaptation of T. yarkandensis, fifteen tissues were collected from three adult T. yarkandensis for transcriptome sequencing and de novo assembly. Nine major international gene annotation databases (NR, COG, egg_NOG, TrEMBL, Pfam, KOG, Swiss prot, KEGG and Gene Ontology) were utilized to annotate unigenes. A detailed study was conducted to explore the gene expression and the differentially expressed genes among five tissues (brain, heart, kidney, liver and spleen). In addition, the current study showed that candidate genes involved in salinity-alkalinity and low temperature adaptation were differentially expressed in tissues of T. yarkandensis. Precisely, mapk1, abcc1, gpx1, gpx4, cat and aqp1 genes participated in the regulation process of salinity-alkalinity adaptation, and elovl4, acaca, fasn, acaa2, acox1 and acox3 genes were involved in fatty acid metabolism and closely associated with low temperature adaptation. On the one hand, it was found that the expression of these genes varied among different tissues, and the important pathways involved in these genes were mapped. Furthermore, we analyzed mapk1 and acox1 genes in depth to obtain the predicted gene structure and important amino acid sites. The transcriptome information in this study will be conducive to provide further understanding for the molecular level research and exploration of the environmental adaptation of T. yarkandensis.

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.

Comparative analyses of liver transcriptomes reveal the effect of exercise on growth-, glucose metabolism-, and oxygen transport-related genes and signaling pathways in grass carp (Ctenopharyngodon idella)

Grass carp is one of the most common farmed fish and its growth rate has been the focus of various studies. However, the impact of long-term exercise on growth rate of juvenile grass carp has not been clearly established. In this study, a four-month exercise trial and liver transcriptome analysis were performed to investigate changes in growth, liver molecular regulatory network and key genes in grass carp. When compared to the non-exercised grass carp (N-EXF), the exercised grass carp (EXF) showed a significant improvement in growth. Liver transcriptome analysis revealed 1714 significantly up-regulated and 1672 significantly down-regulated genes. These genes were enriched in various signaling pathways. These pathways included: those associated with growth, such as the PI3K-Akt and mTOR signaling pathways; those associated with glucose metabolism, such as glycolysis/gluconeogenesis, insulin and AMPK signaling pathways as well as those associated with oxygen transport, such as HIF-1, PI3K-Akt, PPAR and MAPK signaling pathways. In addition, growth-associated genes, such as ghr, igf1 and igf1r; glucose metabolism-associated genes, such as ins and insr as well as oxygen transport-associated genes, such as vhl, pdha and epo were identified. In conclusion, long-term moderate exercise improved the growth rate of grass carp. Our findings elucidate on changes in the liver molecular regulatory network and functional genes that occur during moderate exercise in fish.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

Targeting the Mild-Hypoxia Driving Force for Metabolic and Muscle Transcriptional Reprogramming of Gilthead Sea Bream (Sparus aurata) Juveniles

Simple Summary

Reduced oxygen availability generates a number of adaptive features across all the animal kingdom, and the goal of this study was targeting the mild-hypoxia driving force for metabolic and muscle transcriptional reprogramming of gilthead sea bream juveniles. Attention was focused on blood metabolic and muscle transcriptomic landmarks before and after exhaustive exercise. Our results after mild-hypoxia conditioning highlighted an increased contribution of lipid metabolism to whole energy supply to preserve the aerobic energy production, a better swimming performance regardless of changes in feed intake, as well as reduced protein turnover and improved anaerobic fitness with the restoration of normoxia.

Abstract

On-growing juveniles of gilthead sea bream were acclimated for 45 days to mild-hypoxia (M-HYP, 40–60% O₂ saturation), whereas normoxic fish (85–90% O₂ saturation) constituted two different groups, depending on if they were fed to visual satiety (control fish) or pair-fed to M-HYP fish. Following the hypoxia conditioning period, all fish were maintained in normoxia and continued to be fed until visual satiation for 3 weeks. The time course of hypoxia-induced changes was assessed by changes in blood metabolic landmarks and muscle transcriptomics before and after exhaustive exercise in a swim tunnel respirometer. In M-HYP fish, our results highlighted a higher contribution of aerobic metabolism to whole energy supply, shifting towards a higher anaerobic fitness following normoxia restoration. Despite these changes in substrate preference, M-HYP fish shared a persistent improvement in swimming performance with a higher critical speed at exercise exhaustion. The machinery of muscle contraction and protein synthesis and breakdown was also largely altered by mild-hypoxia conditioning, contributing this metabolic re-adjustment to the positive regulation of locomotion and to the catch-up growth response during the normoxia recovery period. Altogether, these results reinforce the presence of large phenotypic plasticity in gilthead sea bream, and highlights mild-hypoxia as a promising prophylactic measure to prepare these fish for predictable stressful events.

Physiological Effects of Water Flow Induced Swimming Exercise in Seabream Sparus aurata

A longer on-land rearing period of Gilthead seabream Sparus aurata before transfer to sea-cages would allow the farmer to benefit from exercise-enhanced growth, resilience, and robustness as induced by increasing water flow in the tanks. In this study, the physiological effects of flow-conditioning were investigated by subjecting large groups of experimental fish to minimal flow or to flow regimes inducing swimming exercise at 1 or 2 body length (BL) s⁻¹ for a period of 8 months (February–October) in 1,500 L tanks. Fish representing the three treatment groups were then used for: (1) a stress challenge netting test and plasma cortisol measurement (baseline, peaking, and recovery levels), (2) blood plasma measurements of glucose, triglycerides, lactate, cholesterol, growth hormone (GH), and insulin-like growth factor 1 (IGF1), and (3) heart and muscle gene expression of the GH and IGF1 receptors and the muscle transcriptome by deep RNA sequencing (RNAseq). Fish size after 8 months of flow conditioning was 92 ± 27 g body weight (BW) for fish under minimal flow, 106 ± 24 g BW (+15%) at 1 BL s⁻¹, and 125 ± 27 g BW (+36%) at 2 BL s⁻¹. Flow conditioning at 1 BL s⁻¹ provided optimal conditions for growth and uniformity, but also stress (lowest baseline plasma cortisol), robustness (higher condition factor and larger hearts), and energy mobilization (increased plasma glucose). Although flow enhanced growth linearly with swimming speed, also the percentage of lordotic fish increased with exercise, particularly high for swimming at 2 BL s⁻¹. The absence of important differences in plasma GH and IGF1, and expression levels of their receptors in heart and white skeletal muscle, indicated that other factors may be involved in growth enhancement. RNAseq of the white skeletal muscle showed upregulated expression of genes involved in muscle contraction, muscle development and its molecular regulation, and immune genes that may play a role in the muscle repair mechanism. An exercise regime of swimming at 1 BL s⁻¹ can be considered as optimal for farming robust seabream although the increase of skeletal deformities should be avoided.

Identification of mRNA Degradome Variation Dependent on Divergent Muscle Mass in Different Pig Breeds

The search is still on for the molecular processes associated with the development and metabolism of skeletal muscles. Selection conducted in farm animals is focused on high muscle mass because it delivers higher economic profit. The present study aimed to shed light on mRNA degradome signals that could be characteristic for molecular processes associated with an abundance of muscle mass and to identify miRNA regulatory networks controlling these processes in pigs applying next-generation-sequencing (NGS). In the study, over 10,000 degraded transcripts were identified per sample, with the highest abundance for genes encoding mitochondrial proteins (COXs, NDs, CYTB, ATP6 and ATP8). Moreover, only 26% of the miRNA targets were found within this degraded transcript pool, which suggested for miRNAs other molecular mechanism at different level of gene expression than mRNA degradation. On the other hand, a small share of the identified degraded transcripts associated with miRNA regulation suggests a different mechanism of mRNA degradation for identified degraded transcropts. Subsequently, most of the miRNA gene degraded targets, such as ENO3, CKM, CRYAB and ADAM19 encode proteins involved in the muscle mass control. The present study showed an interesting dependence between miRNAs and their targets. Nevertheless, the complete view of the miRNA regulatory network could be a subject of further advanced research, which would employ a miRNA transfection procedure in skeletal muscle cell cultures.

Comparative transcriptome analyses and identification of candidate genes involved in vertebral abnormality of bighead carp Hypophthalmichthys nobilis

Body deformity occurs both in wild and farmed fishes, which is one of the most challenging problems for aquaculture industry. In most cases, such body deformities are linked to skeletal deformities. Currently, very limited information is available on skeletal deformities of farmed fish species which may be caused by genetic factor. In this study, we performed muscle and vertebra transcriptome analyses in body deformity and normality of bighead carp Hypophthalmichthys nobilis (from one meiotic gynogenesis family) using RNA-Seq. A total of 43,923 and 44,416 unigenes were predicted in muscles and vertebrae, respectively. Based on these data, we further explored the gene expression profiles in gynogenetic normal and abnormal bighead carp. No differentially expressed gene (DEG) was found in transcriptome data of muscles. Totally, 20 key DEGs were identified in transcriptome data of vertebrae, such as low density lipoprotein-related protein 2 (lrp2), bone morphogenetic protein 2B (bmp2b) and collagen alpha-1(IV) (col4a1). 12 potential pathways were also identified in vertebra transcriptome data, which were mainly involved in development, growth, cytoskeleton and energy metabolism, such as MAPK signaling pathway, regulation of actin cytoskeleton and TGF-beta signaling pathway. Results of this study will be informative for the understanding of genetic mechanisms for body shape formation and also provide potential candidate genes for selection program involved in body shape and skeletal development in H. nobilis.

Transcriptome analysis of the growth performance of hybrid mandarin fish after food conversion

During recent years, China has become a hotspot for the domestication of mandarin fish, and this is of great commercial value. Although the food preference of domesticated mandarin fish has been studied, little is known about genes regulating their growth. We raised hybrid mandarin fish on artificial feed for 3 months, the results showed that the survival rate of hybrid mandarin fish was 60.00%. Their total length and body weight were 18.34 ±0.43 cm and 100.44 ±4.87 g. The absolute length and weight gain rates were 0.14 cm/d and 1.08 g/d, respectively. Finally, RNA sequencing (RNA-Seq) was performed to identify potential genes and pathways activated in response to growth performance. The transcriptome analysis generated 68, 197 transcripts and 45,871 unigenes. Among them, 1025 genes were up-regulated and 593 genes were down-regulated between the fast- and slow-growth fish. Finally, we obtained 32 differentially expressed genes, which were mainly related to fatty acid biosynthesis (e.g. FASN and ACACB), collecting duct acid secretion (e.g. ATP6E and KCC4), cell cycle (e.g. CDC20 and CCNB), and the insulin-like growth factor (IGF) system (IGFBP1). These pathways might be related to the growth of hybrid mandarin fish. In addition, more potential single nucleotide polymorphisms (SNPs) were detected in the fast-growth fish than in the slow-growth fish. The results suggest that the interaction of metabolism and abundant alleles might determine the growth of hybrid mandarin fish after food conversion.

Identifying the Related Genes of Muscle Growth and Exploring the Functions by Compensatory Growth in Mandarin Fish (Siniperca chuatsi)

How organisms display many different biochemical, physiological processes through genes expression and regulatory mechanisms affecting muscle growth is a central issue in growth and development. In Siniperca chuatsi, the growth-related genes and underlying relevant mechanisms are poorly understood, especially for difference of body sizes and compensatory growth performance. Muscle from 3-month old individuals of different sizes was used for transcriptome analysis. Results showed that 8,942 different expression genes (DEGs) were identified after calculating the RPKM. The DEGs involved in GH-IGF pathways, protein synthesis, ribosome synthesis and energy metabolisms, which were expressed significantly higher in small individuals (S) than large fish (L). In repletion feeding and compensatory growth experiments, eight more significant DEGs were used for further research (GHR2, IGFR1, 4ebp, Mhc, Mlc, Myf6, MyoD, troponin). When food was plentiful, eight genes participated in and promoted growth and muscle synthesis, respectively. Starvation can be shown to inhibit the expression of Mhc, Mlc and troponin, and high expression of GHR2, IGFR1, and 4ebp inhibited growth. Fasting promoted the metabolic actions of GHR2, IGFR1, and 4ebp rather than the growth-promoting actions. MyoD can sense and regulate the hunger, which also worked with Mhc and Mlc to accelerate the compensatory growth of S. chuatsi. This study is helpful to understand the regulation mechanisms of muscle growth-related genes. The elected genes will contribute to the selective breeding in future as candidate genes.

The Estrogen-Responsive Transcriptome of Female Secondary Sexual Traits in the Gulf Pipefish

Sexual dimorphism often results from hormonally regulated trait differences between the sexes. In sex-role-reversed vertebrates, females often have ornaments used in mating competition that are expected to be under hormonal control. Males of the sex-role-reversed Gulf pipefish (Syngnathus scovelli) develop female-typical traits when they are exposed to estrogens. We aimed to identify genes whose expression levels changed during the development and maintenance of female-specific ornaments. We performed RNA-sequencing on skin and muscle tissue in male Gulf pipefish with and without exposure to estrogen to investigate the transcriptome of the sexually dimorphic ornament of vertical iridescent bands found in females and estrogen-exposed males. We further compared differential gene expression patterns between males and females to generate a list of genes putatively involved in the female secondary sex traits of bands and body depth. A detailed analysis of estrogen-receptor binding sites demonstrates that estrogen-regulated genes tend to have nearby cis-regulatory elements. Our results identified a number of genes that differed between the sexes and confirmed that many of these were estrogen-responsive. These estrogen-regulated genes may be involved in the arrangement of chromatophores for color patterning, as well as in the growth of muscles to achieve the greater body depth typical of females in this species. In addition, anaerobic respiration and adipose tissue could be involved in the rigors of female courtship and mating competition. Overall, this study generates a number of interesting hypotheses regarding the genetic basis of a female ornament in a sex-role-reversed pipefish.

Training protocols differently affect AMPK-PGC-1α signaling pathway and redox state in trout muscle

Beneficial effects of physical exercise training are in part related to enhancement of muscle mitochondrial performance. The effects of two different trainings were investigated on transcripts and proteins of the AMPK-PGC-1α signaling pathway, the mitochondrial functioning (citrate synthase (CS), oxidative phosphorylation complexes, uncoupling proteins (UCP)) and the antioxidant defenses (superoxide dismutase (SOD), glutathione peroxidase (GPx), catalase) in rainbow trout red and white skeletal muscles. One group of trouts swam for 10 days at a moderate intensity (approximately 57% U_crit or 2.0 body lengths/s, 23.5 h/day) and another group at a high intensity (approximately 90% U_crit or 3.2 body lengths/s, 2 h/day). In the red muscle, the increase of Cs mRNA levels was significantly correlated with the transcripts of Ampkα1, Ampkα2, Pgc-1α, the oxidative phosphorylation complexes, Ucp2α, Ucp2β, Sod1, Sod2 and Gpx1. After 10 days of training, high intensity training (HIT) stimulates more the transcription of genes involved in this aerobic pathway than moderate intensity training (MIT) in the skeletal muscles, and mainly in the red oxidative muscle. However, no changes in CS, cytochrome c oxidase (COX) and antioxidant defenses activities and in oxidative stress marker (isoprostane plasmatic levels) were observed. The transcriptomic responses are fiber- and training-type dependent when proteins were not yet expressed after 10 days of training. As in mammals, our results suggest that HIT could promote benefit effects in fish.

Early transcriptomic responses associated with the membrane-initiated action of cortisol in the skeletal muscle of rainbow trout (Oncorhynchus mykiss)

Cortisol is a critical neuroendocrine regulator of the stress response in fish. Cortisol practically affects all tissues by interacting with an intracellular receptor and modulating target gene expression. However, cortisol also interacts with components of the plasma membrane in a nongenomic process that activates rapid signaling. Until now, the implication of this novel cortisol signaling for the global transcriptional response has not been explored. In the present work, we evaluated the effects of the membrane-initiated actions of cortisol on the in vivo transcriptome of rainbow trout ( Oncorhynchus mykiss) skeletal muscle. RNA-Seq analyses were performed to examine the transcriptomic changes in rainbow trout stimulated by physiological concentrations of cortisol and cortisol coupled with bovine serum albumin (cortisol-BSA), a membrane-impermeable analog of cortisol. A total of 660 million paired-ends reads were generated. Reads mapped onto the reference genome revealed that 1,737; 897; and 1,012 transcripts were differentially expressed after 1, 3, and 9 h of cortisol-BSA treatment, respectively. Gene Ontology analysis showed that this novel action of cortisol modulates several biological processes, such as mRNA processing, ubiquitin-dependent protein catabolic processes, and transcription regulation. In addition, a KEGG analysis revealed that focal adhesion was the main signaling pathway that was upregulated at all the times tested. Taking these results together, we propose that the membrane-initiated cortisol action contributes significantly in the regulation of stress-mediated gene expression.

Comparative analysis of the transcriptome of the Amazonian fish species Colossoma macropomum (tambaqui) and hybrid tambacu by next generation sequencing

Background

The C. macropomum is a characiform fish from the Amazon basin that has been hybridized with other pacu species to produce commercial hybrids, such as the tambacu. However, little is known of the functional genomics of the parental species or these hybrid forms. The transcriptome of C. macropomum and tambacu were sequenced using 454 Roche platform (pyrosequencing) techniques to characterize the domains of Gene Ontology (GO) and to evaluate the levels of gene expression in the two organisms.

Results

The 8,188,945 reads were assembled into 400,845 contigs. A total of 58,322 contigs were annotated with a predominance of biological processes for both organisms, as determined by Gene Ontology (GO). Similar numbers of metabolic pathways were identified in both the C. macropomum and the tambacu, with the metabolism category presenting the largest number of transcripts. The BUSCO analysis indicated that our assembly was more than 40% complete. We identified 21,986 genes for the two fishes. The P and Log2FC values indicated significant differences in the levels of gene expression, with a total of 600 up-regulated genes.

Conclusion

In spite of the lack of a reference genome, the functional annotation was successful, and confirmed a considerable difference in the specificity and levels of gene expression between the two organisms. This report provides a comprehensive baseline for the genetic management of these commercially important fishes, in particular for the identification of specific genes that may represent markers involved in the immunity, growth, and fertility of these organisms, with potential practical applications in aquaculture management.

De novo transcriptome based on next-generation sequencing reveals candidate genes with sex-specific expression in Arapaima gigas (Schinz, 1822), an ancient Amazonian freshwater fish

Background

The Arapaima (Arapaima gigas) is one of the world's largest freshwater bony fish, and is found in the rivers of the Amazon basin. This species is a potential aquaculture resource, although reproductive management in captivity is limited in particular due to the lack of external sexual dimorphism. In this study, using the 454 Roche platform (pyrosequencing) techniques, we evaluated a major portion of the transcriptome of this important Amazonian species.

Results

Four libraries obtained from the liver and skin tissue of juvenile specimens (representing males and females separately) were sequenced, yielding 5,453,919 high-quality reads. The de novo transcriptome assembly resulted in 175,792 contigs, with 51,057 significant blast hits. A total of 38,586 transcripts were mapped by Gene Ontology using Blast2GO. We identified 20,219 genes in the total transcriptome (9,551 in the liver and 16,818 in the skin). The gene expression analyses indicated 105 genes in the liver and 204 in the skin with differentiated expression profiles, with 95 being over-expressed in the females and 214 in the males. The log2 Fold Change and heatmap based on Reads Per Kilobase per Million mapped reads (RPKM) revealed that the gene expression in the skin is highly differentiated between male and female arapaima, while the levels of expression in the liver are similar between the sexes.

Conclusion

Transcriptome analysis based on pyrosequencing proved to be a reliable tool for the identification of genes with differentiated expression profiles between male and female arapaima. These results provide useful insights into the molecular pathways of sexual dimorphism in this important Amazonian species, and for comparative analyses with other teleosts.

Immunomodulatory Effects of Dietary Seaweeds in LPS Challenged Atlantic Salmon Salmo salar as Determined by Deep RNA Sequencing of the Head Kidney Transcriptome

Seaweeds may represent immuno-stimulants that could be used as health-promoting fish feed components. This study was performed to gain insights into the immunomodulatory effects of dietary seaweeds in Atlantic salmon. Specifically tested were 10% inclusion levels of Laminaria digitata (SW1) and a commercial blend of seaweeds (Oceanfeed^®) (SW2) against a fishmeal based control diet (FMC). Differences between groups were assessed in growth, feed conversion ratio and blood parameters hematocrit and hemoglobin. After a LPS challenge of fish representing each of the three groups, RNAseq was performed on the head kidney as major immune organ to determine transcriptomic differences in response to the immune activation. Atlantic salmon fed with dietary seaweeds did not show major differences in performance in comparison with fishmeal fed fish. RNAseq resulted in ∼154 million reads which were mapped against a NCBI Salmo salar reference and against a de novo assembled S. salar reference for analyses of expression of immune genes and ontology of immune processes among the 87,600 cDNA contigs. The dietary seaweeds provoked a more efficient immune response which involved more efficient identification of the infection site, and processing and presentation of antigens. More specifically, chemotaxis and the chemokine-mediated signaling were improved and therewith the defense response to Gram-positive bacterium reduced. Specific Laminaria digitata effects included reduction of the interferon-gamma-mediated signaling. Highly upregulated and specific for this diet was the expression of major histocompatibility complex class I-related gene protein. The commercial blend of seaweeds caused more differential expression than Laminaria digitata and improved immune processes such as receptor-mediated endocytosis and cell adhesion, and increased the expression of genes involved in response to lipopolysaccharide and inflammatory response. Particularly, expression of many important immune receptors was up-regulated illustrating increased responsiveness. NF-kappa-B inhibitor alpha is an important gene that marked the difference between both seaweed diets as Laminaria digitata inhibits the expression for this cytokine while the blend of seaweeds stimulates it. It can be concluded that the inclusion of seaweeds such as Laminaria digitata can have important modulatory effects on the immune capacity of Atlantic salmon resulting in a more efficient immune response.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Trans-omics approaches used to characterise fish nutritional biorhythms in leopard coral grouper (Plectropomus leopardus)

Aquaculture is now a major supplier of fish, and has the potential to be a major source of protein in the future. Leopard coral groupers are traded in Asian markets as superior fish, and production via aquaculture has commenced. As feeding efficiency is of great concern in aquaculture, we sought to examine the metabolism of leopard coral groupers using trans-omics approaches. Metabolic mechanisms were comprehensively analysed using transcriptomic and metabolomic techniques. This study focused on the dynamics of muscular metabolites and gene expression. The omics data were discussed in light of circadian rhythms and fasting/feeding. The obtained data suggest that branched-chain amino acids played a role in energy generation in the fish muscle tissues during fasting. Moreover, glycolysis, TCA cycles, and purine metabolic substances exhibited circadian patterns, and gene expression also varied. This study is the first step to understanding the metabolic mechanisms of the leopard coral grouper.

Genetic architecture of gene transcription in two Atlantic salmon (Salmo salar) populations

Gene expression regulation has an important role in short-term acclimation and long-term adaptation to changing environments. However, the genetic architecture of gene expression has received much less attention than that of traditional phenotypic traits. In this study, we used a 5 × 5 full-factorial breeding design within each of two Atlantic salmon (Salmo salar) populations to characterize the genetic architecture of gene transcription. The two populations (LaHave and Sebago) are being used for reintroduction efforts into Lake Ontario, Canada. We used high-throughput quantitative real-time PCR to measure gene transcription levels for 22 genes in muscle tissue of Atlantic salmon fry. We tested for population differences in gene transcription and partitioned the transcription variance into additive genetic, non-additive genetic and maternal effects within each population. Interestingly, average additive genetic effects for gene transcription were smaller than those reported for traditional phenotypic traits in salmonids, suggesting that the evolutionary potential of gene transcription is lower than that of traditional traits. Contrary to expectations for early life stage traits, maternal effects were small. In general, the LaHave population had higher additive genetic effects for gene transcription than the Sebago population had, indicating that the LaHave fish have a higher adaptive potential to respond to the novel selection pressures associated with reintroduction into a novel environment. This study highlights not only the profound variation in gene transcription possible among salmonid populations but also the among-population variation in the underlying genetic architecture of such traits.

Analysis of Gene Expression in an Inbred Line of Soft-Shell Clams (Mya arenaria) Displaying Growth Heterosis: Regulation of Structural Genes and the NOD2 Pathway

Mya arenaria is a bivalve mollusk of commercial and economic importance, currently impacted by ocean warming, acidification, and invasive species. In order to inform studies on the growth of M. arenaria, we selected and inbred a population of soft-shell clams for a fast-growth phenotype. This population displayed significantly faster growth (p < 0.0001), as measured by 35.4% greater shell size. To assess the biological basis of this growth heterosis, we characterized the complete transcriptomes of six individuals and identified differentially expressed genes by RNAseq. Pathways differentially expressed included structural gene pathways. Also differentially expressed was the nucleotide-binding oligomerization domain 2 (NOD2) receptor pathway that contributes to determination of growth, immunity, apoptosis, and proliferation. NOD2 pathway members that were upregulated included a subset of isoforms of RIPK2 (mean 3.3-fold increase in expression), ERK/MAPK14 (3.8-fold), JNK/MAPK8 (4.1-fold), and NFκB (4.08-fold). These transcriptomes will be useful resources for both the aquaculture community and researchers with an interest in mollusks and growth heterosis.

Effects of high temperatures on threatened estuarine fishes during periods of extreme drought

Climate change and associated increases in water temperatures may impact physiological performance in ectotherms and exacerbate endangered species declines. We used an integrative approach to assess the impact of elevated water temperature on two fishes of immediate conservation concern in a large estuary system, the threatened longfin smelt (Spirinchus thaleichthys) and endangered delta smelt (Hypomesus transpacificus). Abundances have reached record lows in California, USA, and these populations are at imminent risk of extirpation. California is currently impacted by a severe drought, resulting in high water temperatures, conditions that will become more common as a result of climate change. We exposed fish to environmentally relevant temperatures (14°C and 20°C) and used RNA sequencing to examine the transcriptome-wide responses to elevated water temperature in both species. Consistent with having a lower temperature tolerance, longfin smelt exhibited a pronounced cellular stress response, with an upregulation of heat shock proteins, after exposure to 20°C that was not observed in delta smelt. We detected an increase in metabolic rate in delta smelt at 20°C and increased expression of genes involved in metabolic processes and protein synthesis, patterns not observed in longfin smelt. Through examination of responses across multiple levels of biological organization, and by linking these responses to habitat distributions in the wild, we demonstrate that longfin smelt may be more susceptible than delta smelt to increases in temperatures, and they have little room to tolerate future warming in California. Understanding the species-specific physiological responses of sensitive species to environmental stressors is crucial for conservation efforts and managing aquatic systems globally.

Substantial Downregulation of Myogenic Transcripts in Skeletal Muscle of Atlantic Cod during the Spawning Period

Gonadal maturation is an extremely energy consuming process for batch spawners and it is associated with a significant decrease in growth and seasonal deterioration in flesh quality. Our knowledge about the molecular mechanisms linking sexual maturation and muscle growth is still limited. In the present study, we performed RNA-Seq using 454 GS-FLX pyrosequencing in fast skeletal muscle sampled from two-year-old Atlantic cod (Gadus morhua) at representative time points throughout the reproductive cycle (August, March and May). In total, 126,937 good quality reads were obtained, with 546 nucleotide length and 52% GC content on average. RNA-Seq analysis using the CLC Genomics Workbench with the Atlantic cod reference UniGene cDNA data revealed 59,581 (46.9%) uniquely annotated reads. Pairwise comparison for expression levels identified 153 differentially expressed UniGenes between time points. Notably, we found a significant suppression of myh13 and myofibrillar gene isoforms in fast skeletal muscle during the spawning season. This study uncovered a large number of differentially expressed genes that may be influenced by gonadal maturation, thus representing a significant contribution to our limited understanding of the molecular mechanisms regulating muscle wasting and regeneration in batch spawners during their reproductive cycle.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

mRNA-seq reveals skeletal muscle atrophy in response to handling stress in a marine teleost, the red cusk-eel (Genypterus chilensis)

Background

Fish reared under intensive conditions are repeatedly exposed to stress, which negatively impacts growth. Although most fish follow a conserved pattern of stress response, with increased concentrations of cortisol, each species presents specificities in the cell response and stress tolerance. Therefore, culturing new species requires a detailed knowledge of these specific responses. The red cusk-eel (Genypterus chilensis) is a new economically important marine species for the Chilean aquaculture industry. However, there is no information on the stress- and cortisol-induced mechanisms that decrease skeletal muscle growth in this teleost.

Results

Using Illumina RNA-seq technology, skeletal muscle sequence reads for G. chilensis were generated under control and handling stress conditions. Reads were mapped onto a reference transcriptome, resulting in the in silico identification of 785 up-regulated and 167 down-regulated transcripts. Gene ontology enrichment analysis revealed a significant up-regulation of catabolic genes associated with skeletal muscle atrophy. These results were validated by RT-qPCR analysis for ten candidates genes involved in ubiquitin-mediated proteolysis, autophagy and skeletal muscle growth. Additionally, using a primary culture of fish skeletal muscle cells, the effect of cortisol was evaluated in relation to red cusk-eel skeletal muscle atrophy.

Conclusions

The present data demonstrated that handling stress promotes skeletal muscle atrophy in the marine teleost G. chilensis through the expression of components of the ubiquitin-proteasome and autophagy-lysosome systems. Furthermore, cortisol was a powerful inductor of skeletal muscle atrophy in fish myotubes. This study is an important step towards understanding the atrophy system in non-model teleost species and provides novel insights on the cellular and molecular mechanisms that control skeletal muscle growth in early vertebrates.

Electronic supplementary material

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

Transcriptional dynamics of immune, growth and stress related genes in skeletal muscle of the fine flounder (Paralichthys adpersus) during different nutritional statuses

The effects of stress on immune activity and growth in early vertebrates have not been studied in detail. The present study used fine flounder (Paralichthys adspersus) skeletal muscle as a model to evaluate molecules involved in the stress response, including the glucocorticoid receptors, foxo1/3, and the target genes of these. Additionally, immune markers (il-1β and tnfα) and effector molecules of atrophy (bnip3, caspase-3, and lc3) were assessed. These molecules were analyzed during periods of long-term fasting and refeeding. During fasting, gene expression related to the stress response and atrophy increased; whereas immune markers were down-regulated. During refeeding, atrophy- and stress-related gene expression significantly decreased. In contrast, immune markers were up-regulated. These results provide novel insight on the control of growth in the skeletal muscle of a non-mammalian species under a stressful condition, suggesting that growth, stress, and immune activity in muscle are closely related and coordinated by orchestrated transcriptional dynamics.

Transcriptome and Molecular Pathway Analysis of the Hepatopancreas in the Pacific White Shrimp Litopenaeus vannamei under Chronic Low-Salinity Stress

The Pacific white shrimp Litopenaeus vannamei is a euryhaline penaeid species that shows ontogenetic adaptations to salinity, with its larvae inhabiting oceanic environments and postlarvae and juveniles inhabiting estuaries and lagoons. Ontogenetic adaptations to salinity manifest in L. vannamei through strong hyper-osmoregulatory and hypo-osmoregulatory patterns and an ability to tolerate extremely low salinity levels. To understand this adaptive mechanism to salinity stress, RNA-seq was used to compare the transcriptomic response of L. vannamei to changes in salinity from 30 (control) to 3 practical salinity units (psu) for 8 weeks. In total, 26,034 genes were obtained from the hepatopancreas tissue of L. vannamei using the Illumina HiSeq 2000 system, and 855 genes showed significant changes in expression under salinity stress. Eighteen top Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly involved in physiological responses, particularly in lipid metabolism, including fatty-acid biosynthesis, arachidonic acid metabolism and glycosphingolipid and glycosaminoglycan metabolism. Lipids or fatty acids can reduce osmotic stress in L. vannamei by providing additional energy or changing the membrane structure to allow osmoregulation in relevant organs, such as the gills. Steroid hormone biosynthesis and the phosphonate and phosphinate metabolism pathways were also involved in the adaptation of L. vannamei to low salinity, and the differential expression patterns of 20 randomly selected genes were validated by quantitative real-time PCR (qPCR). This study is the first report on the long-term adaptive transcriptomic response of L. vannamei to low salinity, and the results will further our understanding of the mechanisms underlying osmoregulation in euryhaline crustaceans.

How Many Genes Are Expressed in a Transcriptome? Estimation and Results for RNA-Seq

RNA-seq experiments estimate the number of genes expressed in a transcriptome as well as their relative frequencies. However, an undetermined number of genes can remain undetected due to their low expression relative to the sample size (sequence depth). Estimation of the true number of genes expressed in a transcriptome is essential in order to determine which genes are exclusively expressed in specific tissues or under particular conditions. A reliable estimate of the true number of expressed genes is also required to accurately measure transcriptome changes and to predict the sequencing depth needed to increase the proportion of detected genes. This problem is analogous to ecological sampling problems such as estimating the number of species at a given site. Here we present a non-parametric estimator for the number of undetected genes as well as for the extra sample size needed to detect a given proportion of the undetected genes. Our estimators are superior to ones already published by having smaller standard errors and biases. We applied our method to a set of 32 publicly available RNA-seq experiments, including the evaluation of 311 individually sequenced libraries. We found that in the majority of the cases more than one thousand genes are undetected, and that on average approximately 6% of the expressed genes per accession remain undetected. This figure increases to approximately 10% if individual sequencing libraries are analyzed. Our method is also applicable to metagenomic experiments. Using our method, the number of undetected genes as well as the sample size needed to detect them can be calculated, leading to more accurate and complete gene expression studies.

Global gene expression patterns of grass carp following compensatory growth

BACKGROUND

Compensatory growth is accelerated compared with normal growth and occurs when growth-limiting conditions are overcome. Most animals, especially fish, are capable of compensatory growth, but the mechanisms remain unclear. Further investigation of the mechanism of compensatory growth in fish is needed to improve feeding efficiency, reduce cost, and explore growth-related genes.

RESULTS

In the study, grass carp, an important farmed fish in China, were subjected to a compensatory growth experiment followed by transcriptome analysis by RNA-sequencing. Samples of fish from starved and re-feeding conditions were compared with the control. Under starved conditions, 4061 and 1988 differentially expressed genes (DEGs) were detected in muscle and liver tissue when compared the experimental group with control group, respectively. After re-feeding, 349 and 247 DEGs were identified in muscle and liver when the two groups were compared. Moreover, when samples from experimental group in starved and re-feeding conditions were compared, 4903 and 2444 DEGs were found in muscle and liver. Most of these DEGs were involved in metabolic processes, or encoded enzymes or proteins with catalytic activity or binding functions, or involved in metabolic and biosynthetic pathways. A number of the more significant DEGs were subjected to further analysis. Under fasting conditions, many up-regulated genes were associated with protein ubiquitination or degradation, whereas many down-regulated genes were involved in the metabolism of glucose and fatty acids. Under re-feeding conditions, genes participating in muscle synthesis and fatty acid metabolism were up-regulated significantly, and genes related to protein ubiquitination or degradation were down-regulated. Moreover, Several DEGs were random selected for confirmation by real-time quantitative PCR.

CONCLUSIONS

Global gene expression patterns of grass carp during compensatory growth were determined. To our knowledge, this is a first reported for a teleost fish. The results will enhance our understanding of the mechanism of compensatory growth in teleost fish.

Characterization of the transcriptome of fast and slow muscle myotomal fibres in the pacu (Piaractus mesopotamicus)

BACKGROUND

The Pacu (Piaractus mesopotamicus) is a member of the Characiform family native to the Prata Basin (South America) and a target for the aquaculture industry. A limitation for the development of a selective breeding program for this species is a lack of available genetic information. The primary objectives of the present study were 1) to increase the genetic resources available for the species, 2) to exploit the anatomical separation of myotomal fibres types to compare the transcriptomes of slow and fast muscle phenotypes and 3) to systematically investigate the expression of Ubiquitin Specific Protease (USP) family members in fast and slow muscle in response to fasting and refeeding.

RESULTS

We generated 0.6 Tb of pair-end reads from slow and fast skeletal muscle libraries. Over 665 million reads were assembled into 504,065 contigs with an average length of 1,334 bp and N50 = 2,772 bp. We successfully annotated nearly 47% of the transcriptome and identified around 15,000 unique genes and over 8000 complete coding sequences. 319 KEGG metabolic pathways were also annotated and 380 putative microsatellites were identified. 956 and 604 genes were differentially expressed between slow and fast skeletal muscle, respectively. 442 paralogues pairs arising from the teleost-specific whole genome duplication were identified, with the majority showing different expression patterns between fibres types (301 in slow and 245 in fast skeletal muscle). 45 members of the USP family were identified in the transcriptome. Transcript levels were quantified by qPCR in a separate fasting and refeeding experiment. USP genes in fast muscle showed a similar transient increase in expression with fasting as the better characterized E3 ubiquitin ligases.

CONCLUSION

We have generated a 53-fold coverage transcriptome for fast and slow myotomal muscle in the pacu (Piaractus mesopotamicus) significantly increasing the genetic resources available for this important aquaculture species. We describe significant differences in gene expression between muscle fibre types for fundamental components of general metabolism, the Pi3k/Akt/mTor network and myogenesis, including detailed analysis of paralogue expression. We also provide a comprehensive description of USP family member expression between muscle fibre types and with changing nutritional status.

Transcriptome analysis revealed changes of multiple genes involved in immunity in Cynoglossus semilaevis during Vibrio anguillarum infection

Half-smooth tongue sole (Cynoglossus semilaevis) is one of the most valuable marine aquatic species in Northern China. Given to the rapid development of aquaculture industry, the C. semilaevis was subjected to disease-causing bacteria Vibrio anguillarum. It therefore is indispensable and urgent to understand the mechanism of C. semilaevis host defense against V. anguillarum infection. In the present study, the extensively analysis at the transcriptome level for V. Anguillarum disease in tongue sole was carried out. In total, 94,716 high quality contigs were generated from 75,884,572 clean reads in three libraries (HOSG, NOSG, and CG). 22,746 unigenes were identified when compared with SwissProt, an NR protein database and NT nucleotide database. 954 genes exhibiting the differentially expression at least one pair of comparison in all three libraries were identified. GO enrichment for these genes revealed gene response to biotic stimulus, immune system regulation, and immune response and cytokine production. Further, the pathways such as complement and coagulation cascades and Vibrio cholerae infection pathways were enriched in defensing of pathogen. Besides, 13,428 SSRs and 118,239 SNPs were detected in tongue sole, providing further support for genetic variation and marker-assisted selection in future. In summary, this study identifies several putative immune pathways and candidate genes deserving further investigation in the context of development of therapeutic regimens and lays the foundation for selecting resistant lines of C. semilaevis against V. anguillarum.

Forced sustained swimming exercise at optimal speed enhances growth of juvenile yellowtail kingfish (Seriola lalandi)

Swimming exercise at optimal speed may optimize growth performance of yellowtail kingfish in a recirculating aquaculture system. Therefore, optimal swimming speeds (U_opt in m s⁻¹ or body lengths s⁻¹, BL s⁻¹) were assessed and then applied to determine the effects of long-term forced and sustained swimming at U_opt on growth performance of juvenile yellowtail kingfish. U_opt was quantified in Blazka-type swim-tunnels for 145, 206, and 311 mm juveniles resulting in values of: (1) 0.70 m s⁻¹ or 4.83 BL s⁻¹, (2) 0.82 m s⁻¹ or 3.25 BL s⁻¹, and (3) 0.85 m s⁻¹ or 2.73 BL s⁻¹. Combined with literature data from larger fish, a relation of U_opt (BL s⁻¹) = 234.07(BL)^−0.779 (R² = 0.9909) was established for this species. Yellowtail kingfish, either forced to perform sustained swimming exercise at an optimal speed of 2.46 BL s⁻¹ (“swimmers”) or allowed to perform spontaneous activity at low water flow (“resters”) in a newly designed 3600 L oval flume (with flow created by an impeller driven by an electric motor), were then compared. At the start of the experiment, ten fish were sampled representing the initial condition. After 18 days, swimmers (n = 23) showed a 92% greater increase in BL and 46% greater increase in BW as compared to resters (n = 23). As both groups were fed equal rations, feed conversion ratio (FCR) for swimmers was 1.21 vs. 1.74 for resters. Doppler ultrasound imaging showed a statistically significant higher blood flow (31%) in the ventral aorta of swimmers vs. resters (44 ± 3 vs. 34 ± 3 mL min⁻¹, respectively, under anesthesia). Thus, growth performance can be rapidly improved by optimal swimming, without larger feed investments.

Muscle metabolic remodeling in response to endurance exercise in salmonids

Phenotypic plasticity of skeletal muscle is relevant to swimming performance and metabolism in fishes, especially those that undergo extreme locomotory feats, such as seasonal migration. However, the influence of endurance exercise and the molecular mechanisms coordinating this remodeling are not well understood. The present study examines muscle metabolic remodeling associated with endurance exercise in fed rainbow trout as compared to migrating salmon. Trout were swum for 4 weeks at 1.5 BL/s, a speed similar to that of migrating salmon and red and white muscles were sampled after each week. We quantified changes in key enzymes in aerobic and carbohydrate metabolism [citrate synthase (CS), β-hydroxyacyl-CoA dehydrogenase (HOAD), hexokinase (HK)] and changes in mRNA expression of major regulators of metabolic phenotype (AMPK, PPARs) and lipid (carnitine palmitoyltransferase, CPT I), protein (aspartate aminotransferase, AST) and carbohydrate (HK) oxidation pathways. After 1 week of swimming substantial increases were seen in AMPK and PPARα mRNA expression and of their downstream target genes, CPTI and HK in red muscle. However, significant changes in CS and HK activity occurred only after 4 weeks. In contrast, there were few changes in mRNA expression and enzyme activities in white muscle over the 4-weeks. Red muscle results mimic those found in migrating salmon suggesting a strong influence of exercise on red muscle phenotype. In white muscle, only changes in AMPK and PPAR expression were similar to that seen with migrating salmon. However, in contrast to exercise alone, in natural migration HK decreased while AST increased suggesting that white muscle plays a role in supplying fuel and intermediates possibly through tissue breakdown during prolonged fasting. Dissecting individual and potentially synergistic effects of multiple stressors will enable us to determine major drivers of the metabolic phenotype and their impacts on whole animal performance.

Effects of changes in food supply at the time of sex differentiation on the gonadal transcriptome of juvenile fish. Implications for natural and farmed populations

Background

Food supply is a major factor influencing growth rates in animals. This has important implications for both natural and farmed fish populations, since food restriction may difficult reproduction. However, a study on the effects of food supply on the development of juvenile gonads has never been transcriptionally described in fish.

Methods and Findings

This study investigated the consequences of growth on gonadal transcriptome of European sea bass in: 1) 4-month-old sexually undifferentiated fish, comparing the gonads of fish with the highest vs. the lowest growth, to explore a possible link between transcriptome and future sex, and 2) testis from 11-month-old juveniles where growth had been manipulated through changes in food supply. The four groups used were: i) sustained fast growth, ii) sustained slow growth, iii) accelerated growth, iv) decelerated growth. The transcriptome of undifferentiated gonads was not drastically affected by initial natural differences in growth. Further, changes in the expression of genes associated with protein turnover were seen, favoring catabolism in slow-growing fish and anabolism in fast-growing fish. Moreover, while fast-growing fish took energy from glucose, as deduced from the pathways affected and the analysis of protein-protein interactions examined, in slow-growing fish lipid metabolism and gluconeogenesis was favored. Interestingly, the highest transcriptomic differences were found when forcing initially fast-growing fish to decelerate their growth, while accelerating growth of initially slow-growing fish resulted in full transcriptomic convergence with sustained fast-growing fish.

Conclusions

Food availability during sex differentiation shapes the juvenile testis transcriptome, as evidenced by adaptations to different energy balances. Remarkably, this occurs in absence of major histological changes in the testis. Thus, fish are able to recover transcriptionally their testes if they are provided with enough food supply during sex differentiation; however, an initial fast growth does not represent any advantage in terms of transcriptional fitness if later food becomes scarce.

Transcriptional assessment by microarray analysis and large-scale meta-analysis of the metabolic capacity of cardiac and skeletal muscle tissues to cope with reduced nutrient availability in Gilthead Sea Bream (Sparus aurata L.)

The effects of nutrient availability on the transcriptome of cardiac and skeletal muscle tissues were assessed in juvenile gilthead sea bream fed with a standard diet at two feeding levels: (1) full ration size and (2) 70 % satiation followed by a finishing phase at the maintenance ration. Microarray analysis evidenced a characteristic transcriptomic profile for each muscle tissue following changes in oxidative capacity (heart > red skeletal muscle > white skeletal muscle). The transcriptome of heart and secondly that of red skeletal muscle were highly responsive to nutritional changes, whereas that of glycolytic white skeletal muscle showed less ability to respond. The highly expressed and nutritionally regulated genes of heart were mainly related to signal transduction and transcriptional regulation. In contrast, those of white muscle were enriched in gene ontology (GO) terms related to proteolysis and protein ubiquitination. Microarray meta-analysis using the bioinformatic tool Fish and Chips ( http://fishandchips.genouest.org/index.php ) showed the close association of a representative cluster of white skeletal muscle with some of cardiac and red skeletal muscle, and many GO terms related to mitochondrial function appeared to be common links between them. A second round of cluster comparisons revealed that mitochondria-related GOs also linked differentially expressed genes of heart with those of liver from cortisol-treated gilthead sea bream. These results show that mitochondria are among the first responders to environmental and nutritional stress stimuli in gilthead sea bream, and functional phenotyping of this cellular organelle is highly promising to obtain reliable markers of growth performance and well-being in this fish species.

Transcriptome sequencing of Atlantic salmon (Salmo salar L.) notochord prior to development of the vertebrae provides clues to regulation of positional fate, chordoblast lineage and mineralisation

Background

In teleosts such as Atlantic salmon (Salmo salar L.), segmentation and subsequent mineralisation of the notochord during embryonic stages are essential for normal vertebrae formation. However, the molecular mechanisms leading to segmentation and mineralisation of the notochord are poorly understood. The aim of this study was to identify genes/pathways acting in gradients over time and along the anterior-posterior axis during notochord segmentation and immediately prior to mineralisation of the vertebral bodies in Atlantic salmon.

Results

Notochord samples were collected from unsegmented, pre-segmented and segmented developmental stages. In each stage, the cellular core of the notochord was cut into three pieces along the longitudinal axis (anterior, mid, posterior). RNA was sequenced (22 million pair-end 100 bp/ library) and mapped to the salmon genome. 66569 transcripts were predicted and 55775 were annotated. In order to identify possible gradients leading to segmentation of the notochord, all 71 notochord-expressed hox genes were investigated, most of them displaying a typical anterior-posterior expression pattern along the notochord axis. The clustering of hox genes revealed a pattern that could be related to notochord segmentation. We further investigated how mineralisation is initiated in the notochord, and several factors related to chondrogenic lineage were identified (sox9, sox5, sox6, tgfb3, ihhb and col2a1), suggesting a cartilage-like character of the notochord. KEGG analysis of differentially expressed genes between stages revealed down-regulation of pathways associated with ECM, cell division, metabolism and development at onset of notochord segmentation. This implies that inhibitory signals produce segmentation of the notochord. One such potential inhibitory signal was identified, col11a2, which was detected in segments of non-mineralising notochord.

Conclusions

An incomplete salmon genome was successfully used to analyse RNA-seq data from the cellular core of the Atlantic salmon notochord. In transcriptome we found; hox gene patterns possibly linked to segmentation; down-regulation of pathways in the notochord at onset of segmentation; segmented expression of col11a2 in non-mineralised segments of the notochord; and a chondroblast-like footprint in the notochord.

RNA-seq analysis of early hepatic response to handling and confinement stress in rainbow trout

Fish under intensive rearing conditions experience various stressors which have negative impacts on survival, growth, reproduction and fillet quality. Identifying and characterizing the molecular mechanisms underlying stress responses will facilitate the development of strategies that aim to improve animal welfare and aquaculture production efficiency. In this study, we used RNA-seq to identify transcripts which are differentially expressed in the rainbow trout liver in response to handling and confinement stress. These stressors were selected due to their relevance in aquaculture production. Total RNA was extracted from the livers of individual fish in five tanks having eight fish each, including three tanks of fish subjected to a 3 hour handling and confinement stress and two control tanks. Equal amount of total RNA of six individual fish was pooled by tank to create five RNA-seq libraries which were sequenced in one lane of Illumina HiSeq 2000. Three sequencing runs were conducted to obtain a total of 491,570,566 reads which were mapped onto the previously generated stress reference transcriptome to identify 316 differentially expressed transcripts (DETs). Twenty one DETs were selected for qPCR to validate the RNA-seq approach. The fold changes in gene expression identified by RNA-seq and qPCR were highly correlated (R² = 0.88). Several gene ontology terms including transcription factor activity and biological process such as glucose metabolic process were enriched among these DETs. Pathways involved in response to handling and confinement stress were implicated by mapping the DETs to reference pathways in the KEGG database.

Accession Numbers

Raw RNA-seq reads have been submitted to the NCBI Short Read Archive under accession number SRP022881.

Customized Perl Scripts

All customized scripts described in this paper are available from Dr. Guangtu Gao or the corresponding author.

Fueling the engine: induction of AMP-activated protein kinase in trout skeletal muscle by swimming

AMP-activated protein kinase (AMPK) is well known to be induced by exercise and to mediate important metabolic changes in the skeletal muscle of mammals. Despite the physiological importance of exercise as a modulator of energy use by locomotory muscle, the regulation of this enzyme by swimming has not been investigated in fish. We found that sustained swimming (40 days at 0.75 body lengths s(-1)) increased AMPK activity in red and white trout skeletal muscle (3.9- and 2.2-fold, respectively) as well as the expression of AMPK target genes involved in energy use: lipoprotein lipase and citrate synthase in red and white muscle and CPT1β1b and PGC-1α in red muscle. Furthermore, electrical pulse stimulation of cultured trout myotubes increased AMPK activity and glucose uptake (1.9- and 1.2-fold, respectively) in an AMPK-dependent manner. These results suggest that AMPK may play an important mediatory role in the metabolic adaptation to swimming in fish skeletal muscle.

RNA-Seq technology and its application in fish transcriptomics

High-throughput sequencing technologies, also known as next-generation sequencing (NGS) technologies, have revolutionized the way that genomic research is advancing. In addition to the static genome, these state-of-art technologies have been recently exploited to analyze the dynamic transcriptome, and the resulting technology is termed RNA sequencing (RNA-seq). RNA-seq is free from many limitations of other transcriptomic approaches, such as microarray and tag-based sequencing method. Although RNA-seq has only been available for a short time, studies using this method have completely changed our perspective of the breadth and depth of eukaryotic transcriptomes. In terms of the transcriptomics of teleost fishes, both model and non-model species have benefited from the RNA-seq approach and have undergone tremendous advances in the past several years. RNA-seq has helped not only in mapping and annotating fish transcriptome but also in our understanding of many biological processes in fish, such as development, adaptive evolution, host immune response, and stress response. In this review, we first provide an overview of each step of RNA-seq from library construction to the bioinformatic analysis of the data. We then summarize and discuss the recent biological insights obtained from the RNA-seq studies in a variety of fish species.

Advances in genomics of bony fish

In this review, we present an overview of the recent advances of genomic technologies applied to studies of fish species belonging to the superclass of Osteichthyes (bony fish) with a major emphasis on the infraclass of Teleostei, also called teleosts. This superclass that represents more than 50% of all known vertebrate species has gained considerable attention from genome researchers in the last decade. We discuss many examples that demonstrate that this highly deserved attention is currently leading to new opportunities for answering important biological questions on gene function and evolutionary processes. In addition to giving an overview of the technologies that have been applied for studying various fish species we put the recent advances in genome research on the model species zebrafish and medaka in the context of its impact for studies of all fish of the superclass of Osteichthyes. We thereby want to illustrate how the combined value of research on model species together with a broad angle perspective on all bony fish species will have a huge impact on research in all fields of fundamental science and will speed up applications in many societally important areas such as the development of new medicines, toxicology test systems, environmental sensing systems and sustainable aquaculture strategies.