- Research article
- Open Access
mRNA-seq reveals skeletal muscle atrophy in response to handling stress in a marine teleost, the red cusk-eel (Genypterus chilensis)
BMC Genomics volume 16, Article number: 1024 (2015)
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.
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.
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.
Over the past decades, the aquaculture industry has significantly advanced in terms of technology and rearing methods due to the growing demand of fish for human consumption . While many different species are cultured worldwide, the Chilean finfish industry is highly concentrated in salmonid farming [2, 3]. However, there is a recent local and global trend towards diversifying breeding species to maintain the sustainability of the aquaculture industry [4, 5]. In Chile, one such cultivated marine species is the red cusk-eel (Genypterus chilensis, Guichenot, 1881), a teleost of the Ophidiidae family . This fish is highly valued in national and international markets due to exceptional flesh quality and high nutritional value [7, 8]. Nevertheless, the culturing of this species is seriously hindered by scarce biological knowledge and, primarily, limited information on the negative effects of stress associated with intensive farming in marine species .
The stress response in fish occurs when a stimulus is perceived as a threat. This response increases cortisol secretion, as mediated through the hypothalamic-pituitary-interrenal axis . If the stress is chronic, affected individuals could have permanently increased circulating cortisol levels, a situation that involves short-term metabolic changes and a long-term response associated with reduced growth . An important tissue for growth regulation is skeletal muscle, which is fundamental for an organism’s metabolism and physiology [12, 13]. In mammals, it is well documented that under stressful pathological conditions increased cortisol levels are associated to skeletal muscle atrophy results from increased protein breakdown and decreased protein synthesis , however relatively little is known about this condition in teleost skeletal muscle. The two major route that increases overall rates of protein degradation during muscle atrophy are the ubiquitin-proteasome and the autophagy-lysosome systems . The stimulation by cortisol of these two protein degradation pathways are mediated through the increased expression of several atrogenes, such as Atrogin-1 , Foxo  and MURF-1 , as well as other intracellular mediators related to autophagy such as KLF15 , and REDD1 . Although a few studies have gone further in the understanding of skeletal muscle response in relation to different stressors in teleosts by using PCR arrays [21, 22], microarray [23–26], and RNA-seq approaches [27–30], there are no studies focused to understand the relevance of cortisol in the skeletal muscle response to stress.
It was only recently that the first annotated transcriptome of G. chilensis was published through the use of Illumina technology, thus providing valuable transcriptomic information for this species and for members of the Ophidiidae family . The aim of the present study was to characterize the effect of the stress caused by handling on the skeletal muscle of G. chilensis, an economically important marine fish in the Chilean aquaculture industry. For this, Illumina reads obtained from control and stressed fish were mapped onto the reference transcriptome of the red cusk-eel to identify differentially expressed transcripts (DETs). Gene ontology enrichment analysis revealed a significant up-regulation of catabolic genes associated with skeletal muscle atrophy. Using in vitro approximations, it was also determined that cortisol-mediated genomic actions are a powerful inductor of skeletal muscle atrophy in fish. In the short term, these results will aid in modifying rearing protocols, the aim of which is to improve animal welfare. Importantly, these results also contribute towards the overall understanding of the atrophy system in a non-model teleost, providing novel insights on the cellular and molecular mechanisms that control skeletal muscle growth in early vertebrates.
Physiological and transcriptomic responses of red cusk-eel skeletal muscle to handling stress
Juvenile red cusk-eels were stressed daily with a netting and chasing protocol, and samples were obtained after five days of this procedure. Plasma cortisol and glucose levels significantly increased (5.3-fold and 1.7-fold, respectively) as compared to the control fish, revealing the magnitude of the stress induced by handling. However, no significant differences in plasma lactate levels were observed (Fig. 1). Two separate cDNA libraries were constructed from the skeletal muscle of pooled juvenile red cusk-eels from control and stressed conditions. To examine sequencing variations, cDNA library replicates of each condition were constructed. In total, Illumina MiSeq sequencing generated 23,731,898 paired-end reads (Table 1). Raw data were deposited in the NCBI Sequence Read Archive under Accession Number [GenBank SRS614525: SRR2060847, SRR2063782, SRR2064146, SRR2064148]. After trimming adapters, low quality base pairs, and short reads, the two sequence sets were reduced to 22,992,184 high-quality reads (Table 1). These reads were analyzed with the CLC Genomic Workbench software v.7.0.3 using the previously reported reference transcriptome for G. chilensis , resulting in ~98.8 % of the reads mapped. The expression level of each transcript was represented as RPKM, with 785 transcripts up-regulated under stressed conditions and 167 transcripts down-regulated under stressed conditions. A complete list of the differentially expressed transcripts is included in Additional file 1: Table S1.
Gene ontology (GO) enrichment, KEGG pathway analysis and qPCR validation
The DAVID gene functional classification tool was used to identify groups of transcripts sharing common GO terms. The transcripts up-regulated under stressed conditions were significantly enriched in biological processes, such as modification-dependent macromolecule catabolic process (GO:0043632), modification-dependent protein catabolic process (GO:0019941), and cellular protein catabolic process (GO:0044257), among others (Table 2). The down-regulated transcripts under stressed conditions were significantly enriched in biological processes including striated muscle contraction (GO:0006941), muscle organ development (GO:0007517), and muscle contraction (GO:0006936), among others (Table 2). The GO terms for up-regulated transcripts were significantly enriched in molecular functions, such as ubiquitin-protein ligase activity (GO:0004842), acid-amino acid ligase activity (GO:0016881), and small conjugating protein ligase activity (GO:0019787), among others (Additional file 2: Table S2). The GO terms for down-regulated transcripts were significantly enriched in molecular functions, such as cytoskeletal protein binding (GO:0008092), structural constituent of muscle (GO:0008307), actin binding (GO:0003779), among others (Additional file 2: Table S2). The enrichment of up-regulated transcripts in cellular component distribution were in membrane-enclosed lumen (GO:0031974), organelle lumen (GO:0043233), and intracellular organelle lumen (GO:0070013), among others (Additional file 3: Table S3). The enrichment of down-regulated transcripts in cellular component distribution were in contractile fiber part (GO:0044449), sarcomere (GO:0030017), and contractile fiber (GO:0043292), among others (Additional file 3: Table S3). Finally, pathway analysis through KEGG revealed up-regulated transcripts to be enriched in ubiquitin-mediated proteolysis, regulation of autophagy, and Proteasome (Additional file 4: Table S4). The down-regulated transcripts were enriched in KEGG pathways such as muscle contraction, hypertrophic cardiomyopathy, and dilated cardiomyopathy (Additional file 4: Table S4).
PathVisio 3.0 was used to visualize and integrate the transcriptomic data obtained in KEGG pathway analysis, related to ubiquitin-mediated proteolysis, regulation of autophagy, and muscle contraction (Fig. 2). 21 transcripts associated to the ubiquitin-mediated proteolysis were up-regulated such as components of the E1- ubiquitin activating enzyme (uba1), E2-ubiquitin conjugating enzyme (ube2b, ube2g1, ube2q, ube2r), E3-ubiquitin ligase (ube3a, ube3c, nedd4, herc2, ube4b, mdm2, pirh2, fbx032), as well as components of the 26S proteasome (psmd1, psmd4, psmd5, psmd7, psmd8, psmd11, psmd14, psme1). Similarly, 10 transcripts associated to the regulation of autopaghy were up-regulated such as components of mTOR signaling (ampk, redd1, ulk1, pi3kciii) and components that control autophagosome formation (atg9, atg5, atg16l, atg4, atg7, lc3i). Oppositely, 17 transcripts associated to muscle contraction were down-regulated such as components of the Z-disk (actn2, actn3), myosin heavy chain (myh2, myh3, myh4, myh7, myh10), myosin binding (mybpc1, mybpc3), tropomyosin (tpm1), troponin (tnnt2, tnnt3, tnni1, tnnc1), actin chain (acta1), and myosin light chain (myl3, myl4). Additionally, transcription factors involved in skeletal muscle growth (myod1, myod2) and transcription factors involved in skeletal muscle atrophy (smad2, foxo1) were incorporated.
Eight up-regulated transcripts related to skeletal muscle atrophy (foxo1, ddit4, psmd1, smad2, fbx032, eif4ebp3, atg5, and atg16l1) and two down-regulated transcripts related to skeletal muscle growth (myod1 and myod2) were selected for RT-qPCR analysis so as to validate the results of RNA-seq analysis. The transcript expression fold-changes measured by these two methods (Fig. 3) were highly correlated, with a significant R2 value of 0.85 (p-value = 1E−4). Taken together, these results indicate that stress induces red cusk-eel skeletal muscle atrophy and suggest that cortisol modulates the expression of components from the ubiquitin-proteasome and autophagy systems.
Cellular and molecular response of red cusk-eel myotubes to cortisol
To determine if cortisol was the main modulator of skeletal muscle atrophy in red cusk-eel, the mRNA expression of foxo1, ddit4, psmd1, smad2, fbx032, eif4ebp3, atg5, atg16l, myod1 and myod2 were monitored in myotube lysates 12, 24, and 36 h following cortisol treatment (250 ng/mL) (Fig. 4a-i). At 12 h post-treatment, maximum increases in ddit4 and eif4ebp3 mRNA expressions were observed. At 24 h post-treatment, there was a maximum increase in foxo1 and smad2 mRNA expressions, which was followed by maximum increases of psmd1, fbx032, atg5, and atg16l mRNA expressions at 36 h post-treatment with cortisol. Pretreatment of skeletal myotubes with the glucocorticoid receptor antagonist RU486 (1 μM) significantly inhibited the cortisol-induced up-regulation of foxo1, ddit4, psmd1, smad2, fbx032, eif4ebp3, atg5 and atg16l. No changes in gene expression were observed for myod1 or myod2 following cortisol treatment (data not shown). To determine whether cortisol treatment effectively induced atrophy in fish myotubes, protein ubiquitination levels were recorded 36 h post-treatment. Increased ubiquitination occurred after treatment, and this increase was inhibited by pretreatment with RU486 (Fig. 4j). Additionally, cortisol treatment triggered a significant decrease in myotube diameters, in association with muscular atrophy, at five days post-treatment (Fig. 4k). These results indicate that cortisol, through their genomic action, is a powerful inductor of skeletal muscle atrophy in fish myotubes.
Genypterus chilensis is an economically important marine species for Chilean aquaculture industry. The first reference transcriptome for this species was previously sequenced and annotated through Illumina HiSeq paired-end sequencing . In this study, RNA-seq analysis was used to determine the effects of handling stress on the skeletal muscle of the red cusk-eel, revealing 952 DETs associated with an up-regulation of catabolic biological processes or a down-regulation of skeletal muscle contraction biological processes. The results obtained by this in silico approach were validated by a RT-qPCR analysis of genes representative of the most relevant identified processes, as well as by an in vitro approach that revealed interesting details regarding the regulation dynamics and mechanisms of cortisol-mediated skeletal muscle atrophy in fish. In general terms, the present study found the primary stress response associated to increased levels cortisol and glucose in plasma was very similar to other fish species, such as the rainbow trout (Oncorhynchus mykiss)  or Senegalese sole (Solea senegalensis) .
In teleosts, cortisol is the primary glucocorticoid and regulator of the physiological response to stress . Cortisol binds and activates specific corticoid receptors in target tissues. In fish there are two types of corticoid steroid receptors: mineralocorticoid receptors (MR) and glucocorticoid receptors (GR) . Most fish possess one MR and two GR isoforms, GR1 and GR2 . Once inside the cell, cortisol binds to cytoplasmic GRs to induce conformational changes that cause the receptor to dissociate from respective chaperone molecules. The hormone-receptor complex then translocates to the nucleus where it dimerises and binds to the glucocorticoid response elements of target genes . Depending on the co-factors recruited, this leads to the transactivation or transrepression of these genes. This mode of cortisol action is referred to as the classical or genomic pathway . Cortisol-mediated stress is a key controller in aerobic and anaerobic metabolism , increasing gluconeogenesis, and inhibiting glycogen synthesis . An increased metabolic rate contributes to reduced growth. Cortisol may also control growth by acting through elevated plasma glucose levels to reduce appetite and food ingestion . Due to this, elevated cortisol levels could be inferred through reduced food conversion efficiency, and growth .
For years it was assumed that the low growth rates experienced by fish under stressful conditions were due to metabolic changes impacting anabolic processes . However, the present work is the first to evidence that stress mediated by cortisol induces the expression of components associated with skeletal muscle catabolism. The transcripts up-regulated under stressed conditions were associated with catabolic pathways such as ubiquitin-mediated proteolysis and autophagy, while down-regulated transcripts were associated with muscle contraction or muscle hypertrophy pathways. Interestingly in humans, similar signatures have been observed in skeletal muscle atrophy induced by corticosteroids used to treat several pathological conditions or by cortisol released in response to stress conditions such as sepsis, cachexia, and starvation .
In mammals, in vitro models have been very useful in defining the glucocorticoid pathway by reducing complexity, as compared to an in vivo model [43, 44]. Similarly, fish myotubes primary culture has been used as a valid model for understanding the signaling pathways involved in muscle development and growth in teleost, including for glucose transport , TOR signaling , IGF-1 signaling [47–49], Myostatin signaling [50, 51], GH signaling , proinflammatory cytokine signaling , PAMPs signaling , and cortisol signaling . Similarly, in the present study using an in vitro model, it was determined that cortisol through the expression of components of the ubiquitin-proteasome and autophagy-lysosome systems was an inductor of skeletal muscle atrophy in fish,.
The maximum foxo1expression in myotubes occurred 24 h post-treatment with cortisol. The role of this transcription factor in skeletal muscle atrophy induced by glucocorticoids has been demonstrated in mammals through in vitro and in vivo approximations [56–58]. Moreover, a recent study demonstrated the presence of functional glucocorticoid response elements in the mammalian FoxO1 promoter . In fish, the expression and activity of FoxO1 have been observed under nutritional stress in fine flounder (Paralichthys adspersus) . Interestingly, IGF-1, an inductor of skeletal muscle hypertrophy, induced the phosphorylation of FoxO1 and FoxO4 in trout myotubes in association with decreased transcriptional activity . The described FoxO1 target genes included atrogin-1 (fbx032), a protein involved in muscle proteolysis via the ubiquitin-proteasome system . Atrogin-1 expression occurs in relation to skeletal muscle atrophy under fasting conditions in rainbow trout  and fine flounder . Recently, it was determined that FOXO1/Atrogin-1 signaling pathway is involved in the skeletal muscle atrophy induced by LPS in rainbow trout myotubes .
Other important components of protein catabolism include the 26S proteasome and autophagosome. The up-regulation of the 26S proteasome non-ATPase regulatory subunit 1 (psmd1) has been observed in atrophying fast-twitch muscles from rats treated with dexamethasone ; however, there are no reports on the differential expression of this gene in teleosts. A similar phenomenon has been reported for ATG5 and ATG16L1, where dexamethasone treatment in L6 myotubes induces the expression of several autophagy genes 24 h post-treatment, including of atg5 and atg16l1 . In teleost species, ATG5 activity in zebrafish has been linked to neurogenesis and organogenesis . Nevertheless, there are no reports that relate ATG5 gene expression with fish muscle atrophy.
Another gene with differential expression identified in this study was REDD1 (ddit4), a cytoplasmic protein with a crucial role in repressing protein synthesis mediated by TOR . Additionally, the same signaling pathway showed an over expression of 4EBP-1 (eif4ebp3). This protein represses protein synthesis by directly interacting with the eukaryotic translation initiation factor 4E . In both cases, the present study detected maximum expression 12 h post-treatment with cortisol, suggesting that the inhibition of protein synthesis precedes atrogene expression. In mammals, acute dexamethasone treatment induces REDD1 mRNA expression in rat skeletal muscle in vivo and in L6 myoblasts, as well as down-regulating mTOR signaling through the activation of 4EBP-1 . In zebrafish, REDD1 regulates dorsoventral patterning through the negative modulation of Wnt/β-catenin activity .
Among the genes down-regulated under handling stress, the present study found proteins associated with contractile functions and development, including specific muscle transcription factors such as MyoD1 and MyoD2. Both proteins belong to the myogenic regulatory factors (MRFs) protein family and are involved in the regulation of teleost myogenesis . While both transcripts had decreased expression in RNA-seq analysis, only myod1 showed a significant decrease in its expression in RT-qPCR analysis. Moreover, changes in expression were not detected for myod1 or myod2 in fish myotubes treated with cortisol. This result differs from that observed in mammals, where treatment of the C2C12 myoblast with dexamethasone reduces myod expression . This observation can be explained by the differences between mammalian and teleost myogenesis .
The present work used RNA-seq analysis to determine the effects of handling stress on the skeletal muscle transcriptome of G. chilensis, an important marine fish for the Chilean aquaculture industry. Handling stress induced physiological changes, associated with increases in circulating levels of cortisol, and major changes in global skeletal muscle gene expression. Under stressed conditions, 785 transcripts were up-regulated in association with catabolic signaling pathway such as the ubiquitin-proteasome and autophagy-lysosome systems. Conversely, 167 transcripts were down-regulated under stressed conditions in association with biological processes such as muscle contraction and muscle system processes. The transcriptional expression of components from the ubiquitin-proteasome and autophagy-lysosome pathways were analyzed by RT-qPCR and all were significantly increased under handling stress. To further analyze the relevance of cortisol in teleost skeletal muscle atrophy, an in vitro approximation was performed, providing details about the dynamics and mechanisms by which cortisol induces skeletal muscle atrophy in teleosts. There was a coordinated expression of genes related to the suppression of protein synthesis, protein ubiquitination, and autophagy. All of these were modulated by the genomic actions of cortisol. This is the first study to indicate that stress and cortisol are powerful inductors of skeletal muscle atrophy in a teleost, in addition to providing valuable information for monitoring the culturing and growth of marine fish species under intensive rearing conditions.
Fish sampling and stress conditions
Juvenile red cusk-eels (Genypterus chilensis) with an average weight of 900 ± 50 g and length of 55 ± 5 cm were collected from the Centro de Investigación Marina de Quintay (CIMARQ) (33°13′S 71°38′W, V Region, Valparaíso, Chile). Fish were maintained under natural temperature and light:dark photoperiod conditions (13 °C ± 1 °C and L:D 12:12) for the spring season. The specimens were randomly distributed between a control and stressed group, placed in two separate 90 L tanks, and acclimated for 2 weeks before the experiment. The stressed group was then subjected to a standardized handling stress protocol consisting of netting and chasing the fish for 5 min daily for 5 days . Six hours after the final handling stimulation, stressed and control fish were quickly netted, and blood samples were taken via caudal puncture using heparinized tubes (n = 4 per group). Plasma was collected by centrifugation at 5000 × g for 10 min and stored at −80 °C until analysis. Following blood sampling, fish were sacrificed through an overdose of anesthetic (3-aminobenzoic acid ethyl ester, 300 mg/L), and white muscle was collected, immediately frozen in liquid nitrogen, and stored at −80 °C until analysis. Fish were reared and sampled according to protocols approved by the Bioethical Committee of Andrés Bello University.
Cortisol, glucose, and lactate measurements
Plasma cortisol levels were measured using the Enzyme Immunoassay Kit (Cayman, MI, USA) following the manufacturer’s recommendations. Glucose and lactate plasma levels were measured using Colorimetric Assay Kits (Abcam, Cambridge, UK) following the manufacturer’s recommendations.
Skeletal muscle RNA extraction and transcriptome sequencing
Total RNA was extracted using the RNeasy Mini Kit (Qiagen, TX, USA) following the manufacturer’s instructions. RNA was quantified through spectrophotometry using NanoDrop technology with the Epoch Multi-Volume Spectrophotometer System (BioTek, VT, USA). Total RNA isolated from skeletal muscle was treated with DNase I to remove genomic DNA. RNA concentration was mesuared by Qubit® 2.0 Fluorometer (Life Technology, Carlsbad, CA, USA) and RNA integrity was determined using Fragment Analyzer™ Automated CE System (Analytical Advanced Technologies, Ames, IA, USA). Library construction of 1 μg total RNA for each sample was made using Illumina® TruSeq® RNA Sample prep kit (Illumina®, USA), according to protocol indications, with mRNA fragmentation time of 2 min at 94 °C. After PCR amplification step, mRNA libraries sizes were verified on AATI Fragment Analyzer™, and quantified by qPCR-based quantification using Kapa Library Quantification kit (#KK4824). Libraries were diluted and prepared to a final concentration of 12,5pM for loading on Illumina® MiSeq desktop sequencer, according to the MiSeq System User guide. cDNA library replicates of each condition were constructed and sequenced.
Data processing, differentially expressed transcripts, and GO enrichment analysis
Raw sequencing reads were trimmed by removing Illumina adapter sequences and low quality bases. Sequences shorter than 50 base pairs were also discarded. To analyze in silico gene expression levels, the CLC Genomics Workbench v.7.0.3 software (http://www.clcbio.com/genomics/) was used with the G. chilensis reference transcriptome . RNA-seq analysis was carried out for sequence reads obtained from control and stressed conditions using the following default parameters: two mismatches, minimum fraction length of 0.9, minimum fraction similarity of 0.8, and a maximum of 10 hits per read. Gene expressions were based on reads per kilobase of exon model per million mapped reads (RPKM) values . The up- and down-regulated genes were analyzed for gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways through the DAVID database  and categorized based on GO terms for biological processes, molecular functions, and cellular components, as well as by genes up- or down-regulated in association with KEGG pathways. To establish a relationship between red cusk-eel DETs and DAVID background, a tBLASTx search was conducted against Danio rerio Ensembl proteins for significant matches with the red cusk-eel transcriptome. Danio rerio Ensembl Gene IDs were obtained from the corresponding Ensembl protein entries. DAVID analysis was performed with custom IDs set as the ‘Background’. Standard settings for gene count (2) and ease (0.1) were used. The cut-off p-value for biological processes was 1E−6 and cut-off p-value for molecular functions and cellular components was 1E−3.
Cell cultures and cortisol treatments
Primary myoblasts were isolated from juvenile red cusk-eels. Dorsal white muscle was obtained under sterile conditions and placed in an F10 medium containing 9 mM NaHCO3, 20 mM HEPES, 15 % horse serum, 100 U/ml penicillin, and 10 mg/ml streptomycin at pH 7.4. After mechanical dissociation of the muscle, the tissue was digested with a 0.1 % collagenase and 0.1 % trypsin solution in an F10 medium for 4 h at 15 °C. The suspension was centrifuged at 300 × g for 5 min at 10 °C, and the resulting pellet was resuspended in an F10 medium. The cellular suspension was filtered through 40 μm nylon filters. The cells were centrifuged at 1000 × g for 10 min at 15 °C, resuspended in 3 ml of cold PBS, and layered on top of 4 ml of Ficoll-Paque gradient (GE Healthcare) in a 15 ml tube. Samples were then centrifuged at 1400 × g for 30 min at 15 °C. Following this, the cell layer was extracted with a pipette and washed with 10 ml of F10 medium. Cells were seeded at a density of 2 × 106 per mL in plates previously treated with poly-L-lysine and laminin. Cells were incubated for seven days at 15 °C under atmospheric air in a proliferating medium containing the F10 medium, 9 mM NaHCO3, 20 mM HEPES, 10 % fetal bovine serum, 100 U/ml penicillin, and 10 mg/ml streptomycin. Then, myoblasts were cultivated an additional seven days in a differentiating medium composed by the F10 medium, 9 mM NaHCO3, 20 mM HEPES, 100 U/ml penicillin, and 10 mg/ml streptomycin to obtain differentiated myotubes. On the fourteenth day, myotubes were treated with a medium containing either ethanol (vehicle control), cortisol (250 ng/mL), or RU486 (1 μM). After 12, 24, and 36 h, the medium was removed, and cells were harvested for total RNA isolation or protein extraction. All treatments were performed with n = 4, in two independent experiments.
Total RNA was extracted from skeletal muscle tissue and skeletal myotubes using the RNeasy Mini Kit (Qiagen, TX, USA) following the manufacturer’s recommendations. RNA was quantified by spectrophotometry using NanoDrop technology with the Epoch Multi-Volume Spectrophotometer System (BioTek, VT, USA). Only RNAs with an A260/280 ratio between 1.9 and 2.1 were used for cDNA synthesis. Residual genomic DNA was removed using the genomic DNA wipeout buffer included in the Quantitect Reverse Transcription Kit (Qiagen, TX, USA). Subsequently, 1 μg of RNA was reverse transcribed into cDNA for 30 min at 42 °C using the manufacturer’s recommendations. qPCR analysis was performed with the Stratagene MX3000P qPCR System (Stratagene, CA, USA). Each qPCR reaction mixture contained 7.5 μl of 2× Brilliant II SYBR Master Mix (Stratagene, CA, USA), 6 μl of cDNA (diluted 40-fold), and 250 nM of each primer. For gene expression normalization, the reference gene 40S ribosomal protein S30 (fau) was used . The primers used in this study are listed in Additional file 5: Table S5. Control reactions included a no template control and a control without reverse transcriptase. The QGene program was used for analysis of gene expression .
Western blot analysis
After treatment, cells were solubilized at 4 °C in 30 μL of lysis buffer containing 50 mM Tris-HCl pH 7.4, 150 mM NaCl, 1 mM EDTA, 1 % NP-40, 5 mM Na3VO4, 20 mM NaF, 10 mM sodium pyrophosphate, and a protease inhibitor cocktail (Calbiochem). Proteins extracts were resolved by 10 % SDS-PAGE, transferred to polyvinylidene difluoride membranes (Millipore, MA, USA), and blocked for 1 h at room temperature in Tris-buffered saline, 0.1 % Tween 20, and 5 % fat-free milk. Incubations with primary antibodies (1:1000) were performed at 4 °C overnight. After incubation for 1 h with horseradish peroxidase (HRP)-conjugated secondary antibodies (1:2000), membranes were developed through enhanced chemiluminescence (Amersham Biosciences, Amersham, UK). The films were scanned, and the ImageJ program was employed for densitometric analysis of the bands . Antibodies against ubiquitin (cat. no. 2880), β-actin (cat. no. 4967), and HRP-conjugated secondary anti-rabbit and anti-mouse were obtained from Cell Signaling Technology (MA, USA).
Myotube diameter measurements
Myotubes were loaded with 5.4 mM Calcein AM for 45 min at 15 °C in an F10 medium. Myotubes were transferred to the recording chamber and mounted on an inverted fluorescence microscope. Fluorescence was detected using excitation at 488 nm and emission at 540 nm, and images were collected. The diameters were measured for a total of 50 myotubes from random fields using the ImageJ program . Myotubes were measured at three points along their length, and results were expressed as a percentage of the diameter in the control group (n = 8) in two independent experiments.
Pathway visualization and Heat map summary of atrophy data
The construction of pathway maps of ubiquitin-proteasome, autophagy-lysosome and striated muscle contraction and visualization of differentially expressed gene were performed using PathVisio v3 . The appropriate gene database was selected, and the expression dataset was created after importing the saved input file using the Import Expression Dataset tab.
The relationships among relative expression profiles, a heat map summary and hierarchical clustering analysis of gene expression was performed using Permutmatrix . Clustering and seriation were based on Pearson’s correlation coefficient of z-score normalized relative transcript abundance values. McQuitty’s method as hierarchical clustering was used.
Data are expressed as the mean ± SE. Differences in means between groups were determined using one-way ANOVA followed by Bonferroni’s post-test. Data were accepted as significant at a value of P < 0.05. Correlations between RNA-seq and RT-qPCR data were assessed through multiple linear regression, which obtained the coefficient of determination (R2) and p-value. All statistical analyses were performed using the program GraphPad Prism v.5.00 (GraphPad Software, CA, USA).
Fragments per kilobase of exon model per million mapped reads
Differentially expressed transcript
FAO: Fisheries and Aquaculture Department. The State of World Fisheries and Aquaculture 2012. Rome: FAO; 2012.
Ibieta P, Tapia V, Venegas C, Hausdorf M, Takle H. Chilean Salmon Farming on the Horizon of Sustainability: Review of the Development of a Highly Intensive Production, the ISA Crisis and Implemented Actions to Reconstruct a More Sustainable Aquaculture Industry. In: Sladonja B, editor. Aquaculture and the Environment - A Shared Destiny. Croatia: InTech; 2011.
Iizuka M, Katz J. Globalisation, Sustainability and the Role of Institutions: The Case of the Chilean Salmon Industry: Globalisation, Sustainability and the Role of Institutions. Tijdschr Econ Soc Geogr. 2015;106:140–53.
François NRL, Carter C. Finfish Aquaculture Diversification. New York: CABI; 2010.
Chirivella J, Mañanós EL, Muñoz JL, Soler E. Diversification of farmed species. España: Ministerio de Medio Ambiente y Medio Rural y Marino; 2011.
Nielsen JG, editor. Ophidiiform Fishes of the World (order Ophidiiformes):: An Annotated and Illustrated Catalogue of Pearlfishes, Cusk-eels, Brotulas, and Other Ophidiiform Fishes Known to Date. Rome: Food and Agriculture Organization of the United Nations; 1999 [FAO Species Catalogue, v. 18].
Vega R, Pradenas M, Estrada JM, Ramirez D, Valdebenito I, Mardones A, et al. Evaluacion y comparacion de la eficiencia de dos sistemas de incubacion de huevos de Genypterus chilensis (Guichenot, 1848). Latin Am J Aquat Res. 2012;40:187–200.
Chong J, Gonzalez P. Reproductive cycle and maturity mean size of the red cusk eel, Genypterus chilensis (Guichenot, 1881) in the coast off Talcahuano, Chile. Rev Biol Mar Oceanogr. 1881;2009(44):257–62.
Vijayan MM, Moon TW. The stress response and the plasma disappearance of corticosteroid and glucose in a marine teleost, the sea raven. Can J Zool. 1994;72:379–86.
Wendelaar Bonga SE. The stress response in fish. Physiol Rev. 1997;77:591–625.
Ellis T, Yildiz HY, López-Olmeda J, Spedicato MT, Tort L, Øverli Ø, et al. Cortisol and finfish welfare. Fish Physiol Biochem. 2012;38:163–88.
Johnston IA. Physiology of muscle in hatchery raised fish. Comp Biochem Physiol B Biochem. 1982;73:105–24.
Weatherley AH, Gill HS. Dynamics of increase in muscle fibers in fishes in relation to size and growth. Experientia. 1985;41:353–4.
Braun TP, Marks DL. The regulation of muscle mass by endogenous glucocorticoids. Front Physiol. 2015;6:12.
Sandri M. Protein breakdown in muscle wasting: Role of autophagy-lysosome and ubiquitin-proteasome. Int J Biochem Cell Biol. 2013;45:2121–9.
Gomes MD, Lecker SH, Jagoe RT, Navon A, Goldberg AL. Atrogin-1, a muscle-specific F-box protein highly expressed during muscle atrophy. Proc Natl Acad Sci U S A. 2001;98:14440–5.
Sandri M, Sandri C, Gilbert A, Skurk C, Calabria E, Picard A, et al. Foxo Transcription Factors Induce the Atrophy-Related Ubiquitin Ligase Atrogin-1 and Cause Skeletal Muscle Atrophy. Cell. 2004;117:399–412.
Waddell DS, Baehr LM, van den Brandt J, Johnsen SA, Reichardt HM, Furlow JD, et al. The glucocorticoid receptor and FOXO1 synergistically activate the skeletal muscle atrophy-associated MuRF1 gene. Am J Physiol Endocrinol Metab. 2008;295:E785–97.
Shimizu N, Yoshikawa N, Ito N, Maruyama T, Suzuki Y, Takeda S, et al. Crosstalk between glucocorticoid receptor and nutritional sensor mTOR in skeletal muscle. Cell Metab. 2011;13:170–82.
Britto FA, Begue G, Rossano B, Docquier A, Vernus B, Sar C, et al. REDD1 deletion prevents dexamethasone-induced skeletal muscle atrophy. Am J Physiol Endocrinol Metab. 2014;307:E983–93.
Bermejo-Nogales A, Calduch-Giner JA, Pérez-Sánchez J. Unraveling the molecular signatures of oxidative phosphorylation to cope with the nutritionally changing metabolic capabilities of liver and muscle tissues in farmed fish. PLoS ONE. 2015;10:e0122889.
Calduch-Giner JA, Echasseriau Y, Crespo D, Baron D, Planas JV, Prunet P, et al. 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.). Mar Biotechnol. 2014;16:423–35.
Salem M, Kenney PB, Rexroad CE, Yao J. Microarray gene expression analysis in atrophying rainbow trout muscle: a unique nonmammalian muscle degradation model. Physiol Genomics. 2006;28:33–45.
Johansen KA, Sealey WM, Overturf K. The effects of chronic immune stimulation on muscle growth in rainbow trout. Comp Biochem Physiol B Biochem Mol Biol. 2006;144:520–31.
Magnoni LJ, Roher N, Crespo D, Krasnov A, Planas JV. In Vivo Molecular Responses of Fast and Slow Muscle Fibers to Lipopolysaccharide in a Teleost Fish, the Rainbow Trout (Oncorhynchus mykiss). Biology (Basel). 2015;4:67–87.
Nesan D, Vijayan MM. The transcriptomics of glucocorticoid receptor signaling in developing zebrafish. PLoS ONE. 2013;8:e80726.
Garcia de la Serrana D, Devlin RH, Johnston IA. RNAseq analysis of fast skeletal muscle in restriction-fed transgenic coho salmon (Oncorhynchus kisutch): an experimental model uncoupling the growth hormone and nutritional signals regulating growth. BMC Genomics. 2015;16:564.
Sun S, Xuan F, Fu H, Zhu J, Ge X, Gu Z. Transciptomic and histological analysis of hepatopancreas, muscle and gill tissues of oriental river prawn (Macrobrachium nipponense) in response to chronic hypoxia. BMC Genomics. 2015;16:491.
Salem M, Manor ML, Aussanasuwannakul A, Kenney PB, Weber GM, Yao J. Effect of sexual maturation on muscle gene expression of rainbow trout: RNA-Seq approach. Physiol Rep. 2013;1:e00120.
Palstra AP, Beltran S, Burgerhout E, Brittijn SA, Magnoni LJ, Henkel CV, et al. Deep RNA Sequencing of the Skeletal Muscle Transcriptome in Swimming Fish. PLoS ONE. 2013;8:e53171.
Aedo JE, Maldonado J, Estrada JM, Fuentes EN, Silva H, Gallardo-Escarate C, et al. Sequencing and de novo assembly of the red cusk-eel (Genypterus chilensis) transcriptome. Mar Genomics. 2014;18:105–7.
Krasnov A, Koskinen H, Pehkonen P, Rexroad CE, Afanasyev S, Mölsä H. Gene expression in the brain and kidney of rainbow trout in response to handling stress. BMC Genomics. 2005;6:3.
Cordeiro OD, Silva TS, Alves RN, Costas B, Wulff T, Richard N, et al. Changes in liver proteome expression of Senegalese sole (Solea senegalensis) in response to repeated handling stress. Mar Biotechnol. 2012;14:714–29.
Pankhurst NW. The endocrinology of stress in fish: an environmental perspective. Gen Comp Endocrinol 2011;170:265–275.
Greenwood AK, Butler PC, White RB, DeMarco U, Pearce D, Fernald RD. Multiple corticosteroid receptors in a teleost fish: distinct sequences, expression patterns, and transcriptional activities. Endocrinology. 2003;144:4226–36.
Baker ME, Funder JW, Kattoula SR. Evolution of hormone selectivity in glucocorticoid and mineralocorticoid receptors. J Steroid Biochem Mol Biol. 2013;137:57–70.
Aluru N, Vijayan MM. Stress transcriptomics in fish: A role for genomic cortisol signaling. Gen Comp Endocrinol. 2009;164:142–50.
De Boeck G, Alsop D, Wood C. Cortisol effects on aerobic and anaerobic metabolism, nitrogen excretion, and whole-body composition in juvenile rainbow trout. Physiol Biochem Zool. 2001;74:858–68.
Milligan CL. A regulatory role for cortisol in muscle glycogen metabolism in rainbow trout Oncorhynchus mykiss Walbaum. J Exp Biol. 2003;206(Pt 18):3167–73.
Bernier NJ, Bedard N, Peter RE. Effects of cortisol on food intake, growth, and forebrain neuropeptide Y and corticotropin-releasing factor gene expression in goldfish. Gen Comp Endocrinol. 2004;135:230–40.
Oliva-Teles A. Nutrition and health of aquaculture fish: Nutrition and health of aquaculture fish. J Fish Dis. 2012;35:83–108.
Calura E, Cagnin S, Raffaello A, Laveder P, Lanfranchi G, Romualdi C. Meta-analysis of expression signatures of muscle atrophy: gene interaction networks in early and late stages. BMC Genomics. 2008;9:630.
Croxtall JD, van Hal PTW, Choudhury Q, Gilroy DW, Flower RJ. Different glucocorticoids vary in their genomic and non-genomic mechanism of action in A549 cells. Br J Pharmacol. 2002;135:511–9.
Sanden S, Tripmacher R, Weltrich R, Rohde W, Hiepe F, Burmester GR, et al. Glucocorticoid dose dependent downregulation of glucocorticoid receptors in patients with rheumatic diseases. J Rheumatol. 2000;27:1265–70.
Díaz M, Vraskou Y, Gutiérrez J, Planas JV. Expression of rainbow trout glucose transporters GLUT1 and GLUT4 during in vitro muscle cell differentiation and regulation by insulin and IGF-I. Am J Physiol Regul Integr Comp Physiol. 2009;296:R794–800.
Seiliez I, Gabillard J-C, Skiba-Cassy S, Garcia-Serrana D, Gutiérrez J, Kaushik S, et al. An in vivo and in vitro assessment of TOR signaling cascade in rainbow trout (Oncorhynchus mykiss). Am J Physiol Regul Integr Comp Physiol. 2008;295:R329–35.
Castillo J, Le Bail P-Y, Paboeuf G, Navarro I, Weil C, Fauconneau B, et al. IGF-I binding in primary culture of muscle cells of rainbow trout: changes during in vitro development. Am J Physiol Regul Integr Comp Physiol. 2002;283:R647–52.
Castillo J, Codina M, Martínez ML, Navarro I, Gutiérrez J. Metabolic and mitogenic effects of IGF-I and insulin on muscle cells of rainbow trout. Am J Physiol Regul Integr Comp Physiol. 2004;286:R935–41.
Castillo J, Ammendrup-Johnsen I, Codina M, Navarro I, Gutiérrez J. IGF-I and insulin receptor signal transduction in trout muscle cells. Am J Physiol Regul Integr Comp Physiol. 2006;290:R1683–90.
Garikipati DK, Rodgers BD. Myostatin stimulates myosatellite cell differentiation in a novel model system: evidence for gene subfunctionalization. Am J Physiol Regul Integr Comp Physiol. 2012;302:R1059–66.
Garikipati DK, Rodgers BD. Myostatin inhibits myosatellite cell proliferation and consequently activates differentiation: evidence for endocrine-regulated transcript processing. J Endocrinol. 2012;215:177–187.
Levesque HM, Shears MA, Fletcher GL, Moon TW. Myogenesis and muscle metabolism in juvenile Atlantic salmon (Salmo salar) made transgenic for growth hormone. J Exp Biol. 2008;211(Pt 1):128–37.
Pooley NJ, Tacchi L, Secombes CJ, Martin SA. Inflammatory responses in primary muscle cell cultures in Atlantic salmon (Salmo salar). BMC Genomics. 2013;14:747.
Aedo JE, Reyes AE, Avendaño-Herrera R, Molina A, Valdés JA. Bacterial lipopolysaccharide induces rainbow trout myotube atrophy via Akt/FoxO1/Atrogin-1 signaling pathway. Acta Biochim Biophys Sin. 2015;47:932–7.
Galt NJ, Froehlich JM, Remily EA, Romero SR, Biga PR. The effects of exogenous cortisol on myostatin transcription in rainbow trout, Oncorhynchus mykiss. Comp Biochem Physiol A Mol Integr Physiol. 2014;175:57–63.
Cho JE, Fournier M, Da X, Lewis MI. Time course expression of Foxo transcription factors in skeletal muscle following corticosteroid administration. J Appl Physiol. 2010;108:137–45.
Gilson H, Schakman O, Combaret L, Lause P, Grobet L, Attaix D, et al. Myostatin gene deletion prevents glucocorticoid-induced muscle atrophy. Endocrinology. 2007;148:452–60.
Sanchez AMJ, Candau RB, Bernardi H. FoxO transcription factors: their roles in the maintenance of skeletal muscle homeostasis. Cell Mol Life Sci. 2014;71:1657–71.
Qin W, Pan J, Qin Y, Lee DN, Bauman WA, Cardozo C. Identification of functional glucocorticoid response elements in the mouse FoxO1 promoter. Biochem Biophys Res Commun. 2014;450:979–83.
Fuentes EN, Ruiz P, Valdes JA, Molina A. Catabolic signaling pathways, atrogenes, and ubiquitinated proteins are regulated by the nutritional status in the muscle of the fine flounder. PLoS ONE. 2012;7:e44256.
Seiliez I, Sabin N, Gabillard J-C. FoxO1 is not a key transcription factor in the regulation of myostatin (mstn-1a and mstn-1b) gene expression in trout myotubes. Am J Physiol Regul Integr Comp Physiol. 2011;301:R97–104.
Cleveland BM, Evenhuis JP. Molecular characterization of atrogin-1/F-box protein-32 (FBXO32) and F-box protein-25 (FBXO25) in rainbow trout (Oncorhynchus mykiss): Expression across tissues in response to feed deprivation. Comp Biochem Physiol B Biochem Mol Biol. 2010;157:248–57.
Combaret L, Taillandier D, Dardevet D, Béchet D, Rallière C, Claustre A, et al. Glucocorticoids regulate mRNA levels for subunits of the 19 S regulatory complex of the 26 S proteasome in fast-twitch skeletal muscles. Biochem J. 2004;378(Pt 1):239–46.
Troncoso R, Paredes F, Parra V, Gatica D, Vásquez-Trincado C, Quiroga C, et al. Dexamethasone-induced autophagy mediates muscle atrophy through mitochondrial clearance. Cell Cycle. 2014;13:2281–95.
Hu Z, Zhang J, Zhang Q. Expression pattern and functions of autophagy-related gene atg5 in zebrafish organogenesis. Autophagy. 2011;7:1514–27.
Schakman O, Kalista S, Barbé C, Loumaye A, Thissen JP. Glucocorticoid-induced skeletal muscle atrophy. Int J Biochem Cell Biol. 2013;45:2163–72.
Wang H, Kubica N, Ellisen LW, Jefferson LS, Kimball SR. Dexamethasone represses signaling through the mammalian target of rapamycin in muscle cells by enhancing expression of REDD1. J Biol Chem. 2006;281:39128–34.
Feng Q, Zou X, Lu L, Li Y, Liu Y, Zhou J, et al. The stress-response gene redd1 regulates dorsoventral patterning by antagonizing Wnt/β-catenin activity in zebrafish. PLoS ONE. 2012;7:e52674.
Gabillard JC, Sabin N, Paboeuf G. In vitro characterization of proliferation and differentiation of trout satellite cells. Cell Tissue Res. 2010;342:471–7.
Te Pas MF, de Jong PR, Verburg FJ. Glucocorticoid inhibition of C2C12 proliferation rate and differentiation capacity in relation to mRNA levels of the MRF gene family. Mol Biol Rep. 2000;27:87–98.
Rossi G, Messina G. Comparative myogenesis in teleosts and mammals. Cell Mol Life Sci. 2014;71:3081–99.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4:44–57.
Fuentes EN, Safian D, Valdés JA, Molina A. Isolation and selection of suitable reference genes for real-time PCR analyses in the skeletal muscle of the fine flounder in response to nutritional status: assessment and normalization of gene expression of growth-related genes. Fish Physiol Biochem. 2013;39:765–77.
Simon P. Q-Gene: processing quantitative real-time RT-PCR data. Bioinformatics. 2003;19:1439–40.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.
Van Iersel MP, Kelder T, Pico AR, Hanspers K, Coort S, Conklin BR, et al. Presenting and exploring biological pathways with PathVisio. BMC Bioinformatics. 2008;9:399.
Caraux G, Pinloche S. PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005;21:1280–1.
This study was funded by the National Commission for Scientific and Technological Research (CONICYT), FONDAP projects 15110027 and 15090007.
The authors declare that they have no competing interests.
AM and JV conceived and designed the study. JV wrote the manuscript. CM and MB performed the Illumina Sequencing. JA and JM analyzed the RNA Seq data and performed statistical analyses. JA and VA performed the RT-qPCR and in vitro experiments. CG, HS and JME contributed analysis, text and comments to the paper. All authors read and approved the final manuscript.
Complete list of up-regulated and down-regulated transcripts in response to handling stress. (XLSX 66 kb)
Enriched Molecular functions of up-regulated and down-regulated transcripts in response to handling stress. (XLSX 10 kb)
Enriched cellular component of up-regulated and down-regulated transcripts in response to handling stress. (XLSX 10 kb)
Enriched KEGG pathways of up-regulated and down-regulated transcripts in response to handling stress. (XLSX 10 kb)
Primer sequences for qPCR assay, amplicon size and PCR. (XLSX 10 kb)
About this article
Cite this article
Aedo, J.E., Maldonado, J., Aballai, V. et al. mRNA-seq reveals skeletal muscle atrophy in response to handling stress in a marine teleost, the red cusk-eel (Genypterus chilensis). BMC Genomics 16, 1024 (2015). https://doi.org/10.1186/s12864-015-2232-7
- Genypterus chilensis
- Red cusk-eel
- Handling stress
- Skeletal muscle atrophy