- Research article
- Open Access
Characterization of the transcriptome of fast and slow muscle myotomal fibres in the pacu (Piaractus mesopotamicus)
BMC Genomicsvolume 16, Article number: 182 (2015)
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.
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.
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.
Fish myotomes are composed of anatomically segregated muscle fibre types each with distinct contractile and metabolic phenotypes . Based on their contractile speed skeletal muscle fibres are classified as slow, intermediate or fast [2,3]. Fast twitch (white) muscle fibres comprise the bulk of the myotome and are recruited for energetic movements associated with prey capture and escape behaviour . Fast fibres have elevated densities of myofibrils, reduced myoglobin content and higher a capillary density than slow fibres and utilise phosphagen breakdown and anaerobic glycolysis to power contraction . Sustained swimming activity is supported by superficial layers of slow (red) and intermediate (pink) twitch fibres which are recruited at slow and high cruising speeds respectively . Slow fibres have extensive lipid and glycogen stores, abundant mitochondria and high capillary densities reflecting their reliance on aerobic metabolism . Intermediate (pink) muscle fibre types are found between the slow and fast fibre layers, but express distinct isoforms of myosin heavy chains  and have intermediate contraction speeds and metabolic characteristics [3,7].
Muscle growth reflects the balance between protein synthesis and degradation. These two processes are influenced by numerous biotic and abiotic factors including food availability, growth factors, age, sex, diet composition, swimming activity, oxygen saturation, light and temperature [8,9]. The Insulin-like growth factor (Igf) network, composed of Igfs, binding proteins (Igfbp) and receptors (Igf1r and Igf2r), plays a pivotal role in integrating internal and external inputs to regulate muscle mass . Igf1 regulates several signalling pathways including the Pi3k/Akt/mTor network that controls protein synthesis [8,9]. Typically, fibre production continues until 45% of the maximum body length of the fish, and subsequent growth is entirely by fibre elongation and hypertrophy [10-12]. Myogenesis involves the activation, proliferation and fusion of a resident myoblast population involving hundreds of structural and regulatory genes [10,11].
Muscle protein degradation occurs through three major pathways [13,14] namely: membrane-bound lysosomal enzymes, calpain proteases [14,15], and the Ubiquitin Proteasome (Ub) Pathway (UPP) . UPP represents the most important system for degradation of unnecessary or damaged proteins. Targeted proteins are linked to ubiquitin, which acts as a recognition tag for the proteasome. Ubiquitin tagging of proteins requires the coordinated activity of three classes of enzymes known as E1, E2 and E3. It is the E3 enzymes, also known as E3-ubiquitin ligases, which conjugate ubiquitin to the target protein . Ubiquitin mediated degradation can be reversed through action of deubiquitinating enzymes (DUB), a large group of proteases that cleave ubiquitin-protein conjugates removing the UPP signal, and play an essential role in the regulation of protein degradation . DUBs are subdivided into four families: ubiquitin C-terminal hydrolases (UCHs), ovarian tumour proteases (OTUs), JAB1/MPN/MOV34 metalloenzymes and ubiquitin-specific proteases (USPs) . USP is the largest family of DUB and regulates a wide variety of cellular processes. Although the essential role of USPs in protein degradation is well established, less is known about the function and regulation of specific family members: for example Usp7 has been associated with p53 and Akt turnover, Usp8 with receptor tyrosine endocytosis, Usp33 with the Von Hippel–Lindau disease (VHL) pathway and Usp19 is thought to have a role in muscle development . There have been several studies of the expression of protein degradation related genes in fish, mostly in response to varying nutrition and focused on E3-ubiquitin ligases (particularly Fbox32 and Murf1) [20,21]. In contrast, nothing is known about the transcriptional regulation of USP family members.
The pacu (Piaractus mesopotamicus) is a member of the characiform family (SuperOder Ostariophysi) native to the Prata Basin (South America) and is a target for finfish aquaculture in Brazil. One of the main limitations for the development of a selective breeding program for this species is a lack of genetic information as well as limited knowledge about its physiology. Next Generation Sequencing (NGS) technologies have dramatically increased the amount of sequence data for teleosts and helped to overcome the lack of annotated genomes (so far only 12 fish annotated genomes are publicly available) .
The primary objectives of the present study were 1) to increase the genetic resources available for pacu 2) to exploit the anatomical separation of fibre types to characterise the expression signatures of fast and slow muscles and 3) to systematically investigate the expression of USPs in fast and slow muscle in response to fasting and refeeding. Teleost fish underwent a teleost-specific genome duplication (TSGD) event around 450 million years ago (Mya) which was followed by diploidisation and gene loss [23,24]. It is thought that around 15 to 20% of TSGD paralogues have been retained in the diploid genome of extant species [23-26]. Several studies have demonstrated that paralogues from the TSGD or the salmonid-specific whole-genome duplication (WGD), which occurred 88Mya , can display different patterns of expression during myogenesis and muscle growth [27-29]. A secondary objective was therefore to search for evidence of differential expression between teleost specific paralogues within and between muscle fibre types.
Results and discussion
De novo assembly
Individual barcoded libraries from fast and slow skeletal muscle tissues were generated from 5 adult pacu (1.5 ± 0.61 kg; mean ± SD) and sequenced using HiSeq 2000 platform yielding a total of 374,952,267 and 390,706,230 pair-end reads respectively (~0.6 Tb) (raw reads were deposited in the European Nucleotide Archive accession number PRJEB6656). 665,042,497 trimmed paired-end reads (86% of the total) were successfully assembled into 504,065 contigs with an average length of 1,334 bp and N50 = 2,772 bp (Table 1). The present study significantly improves upon previous teleost transcriptomes . A total of 237,637 contigs (47%) were successfully annotated, representing over 15,000 unique genes with an average coverage of 53 and over 8,000 genes with >90% of the coding sequence (CDS) represented (Table 1; annotation results can be found in Additional file 1). The number of genes found represent between 56-65% of all protein-coding genes currently identified in Danio rerio (Cypriniformes) and Astyanax mexicanus (Characiformes) the only two other Ostariophysi genomes available . Pacu contigs were also annotated and classified into 319 different vertebrate signalling and metabolic pathways (Additional file 2). In some cases, such as for the Pi3k/Akt/mTor network we found that over 90% of the genes were represented.
The current transcriptome dramatically increases the genetic resources available for the future development of genetic-based breeding programmes in pacu. Simple sequence repeats (SSRs) or microsatellites are widely used for parentage identification and stock management in family selection programs . Screening for SSRs was focused on those contigs for which over 90% of the CDS was present to assure that repeated sequences were correctly identified in the UTR regions. A total of 380 SSRs were identified in the transcriptome, and the great majority of them (54%) were repeats of a dinucleotide motif (Table 2; Additional file 3) increasing the number of potential microsatellites for this species . We also screened the annotated transcriptome to identify TSGD-paralogues. It has been estimated that around 15 to 20% of teleost-specific paralogues have been retained in the diploid genomes of extant species [23-28] with a slightly higher proportion retained in Ostariphysi genomes . We identified a total of 442 paralogues pairs in the pacu when compare with zebrafish (884 genes; Figure 1A, Additional file 4). Although this is a significant improvement when compared with previous studies  the number of paralogues identified was still lower than what we could expect from a 16,000 gene transcriptome.
Digital gene expression analysis
Individual fast and slow muscle libraries were mapped to the complete transcriptome. Mapped reads were normalized by contig length, library size and only those with an average of more than 10 mapped were used to estimate digital gene expression (DGE). Whereas there are several transcriptomic studies in teleosts [29,32,33] to our knowledge this is the first one to compare global patterns of DGE in muscle different fibre types. Global DGE analysis revealed 956 and 604 genes differentially expressed between slow and fast skeletal muscle respectively (Additional file 5). Gene ontology (GO) enrichment analysis of the differently expressed genes (Table 3; Additional file 6) revealed a significant enrichment in genes involved in oxidative phosphorylation and lipid metabolism (GO:0005739; 0006629; 0009055; 0005811; 0018925; Table 3) in slow skeletal muscle and carbohydrate metabolism and kinase activity (GO:0005975; 0019752; 0016301) in fast skeletal muscle. These results are in agreement with the main differences described between slow and fast muscle with respect to their predominant means of energy supply . To further validate GO analysis results, reads mapped were individually explored in all components related to general metabolic pathways present in our transcriptome: glycolysis, ß-oxidation, Krebs cycle and the electron transport chain (Additional file 7). The majority of glycolysis-related enzymes found, such as glucose-6-phosphate isomerase (gpi), 6-phosphofructokinase (pfkm) or enolase (eno), were significantly more highly expressed in fast skeletal muscle (FDR < 0.05) (Additional file 7). As expected slow muscle had a significantly higher DGE in almost all components of the Krebs cycle, ß-oxidation pathway and the electron transport chain (Additional file 8). Therefore DGE results captured fundamental physiological and metabolic differences between tissues. We also studied DGE in relation to the Pi3k/mTor pathway, myogenesis-related genes and TSGD-paralogues between fibre types. Previous studies have reported that TSGD-paralogues can display different expression profiles in muscle during myogenesis and in response to varying nutrition, however those studies were always restricted to fast skeletal muscle and to a limited number of paralogues [34,35]. By using DGE we found that 301 and 245 TSGD-paralogues pairs were differentially expressed for slow and fast skeletal muscle respectively (FDR < 0.05; example in Figure 1B; Additional file 9). Similarly, we found 124 paralogues with significant differences in DGE between muscle types (FDR < 0.05; Figure 1B). Various mechanisms have been postulated to explain TSGD paralogue retention during evolution including mutations leading to differential regulation in expression (sub-functionalization) and/or the evolution of some novel function (neofunctionalization) . Whilst we cannot distinguish between these possibilities in the present study it is apparent that TSGD paralogue retention has significantly contributed to phenotypic patterning of gene expression between fast and slow muscle fibre types.
The Pi3k/Akt/mTor network is involved in several cellular processes: muscle growth, cell cycle, muscle differentiation and myoblast proliferation [8,37,38]. The majority of the pathway’s components were identified in the transcriptome, including several paralogues (Figure 2; Additional file 10). We found differences in DGE between muscle types for 28 components (FDR < 0.05; Figure 2). Only 6 genes were more abundant in fast skeletal muscle including insulin-like growth factor receptor 1a (igfr1a) and tuberous sclerosis 2 (tsc2) (Figure 2, empty circles). A total of 22 components of the pathway were found to be more highly expressed in slow skeletal muscle including insulin-like growth factor 2b (igf2b), several insulin-like growth factors binding proteins (igfbp1a, igfbp2a, igfbp5a and igfbp7), protein kinase B gamma (akt3), ribosomal protein S6 kinase beta 1a (rps6kb1a), eukaryotic translation initiation factor 4E-binding protein 1 (eif4ebp1) and mitogen-activated protein kinases (mapk14a and mapk14b) (Figure 2, red circles). Those genes significantly up regulated in slow skeletal muscle are directly involved in the stimulation of protein synthesis, suggesting a higher protein synthesis potential in this fibre type [39,40]. Some of the genes significantly up regulated in fast skeletal muscle are associated with the regulation of protein synthesis during fasting. For instance, lkb1 is expressed when energy levels are low, mlst8 binds mTor to stabilize the complex when amino acids are not optimal and tsc2 is related with protein synthesis inhibition .
DGE was also analysed for the genes associated with myoblast activation, proliferation, differentiation and migration (Figure 3A). Stac3 is a multifunctional signal adaptor protein which shows a strong correlation with myogenin expression both in cell culture and in vivo and is required for normal myotube formation and sarcomere assembly . Stac3 localises to T-tubules and also functions in the mediation of voltage-induced Ca2+ release and contractility . The higher expression of stac3, in correlation with myogenin, in fast than slow muscle (Figure 3E) may therefore reflect the more extensive development of T-tubules and sarcoplasmic reticulum in this fibre type which are required to achieve shorter contraction cycles during high speed swimming . In contrast, the other myogenic regulatory factors (Figure 3B) and members of the myocyte enhancer gene family (except mef2ca) (Figure 3C) had similar expression patterns in fast and slow muscles. Differences in expression for myostatin paralogues has been previously reported in other fish species including rainbow trout (Oncorhynchus mykiss) and tilapia (Oreochromis niloticus) with mstna expression (also known as mstn2) confined to brain, testes and spleen and mstnb (also known as mstn1) more expressed in heart and muscle [44,45].
Transcriptional regulation of USP and E3-ubiquitin ligases
The quality of the present transcriptome allowed us to study the regulation of the Ubiquitin Specific Proteases (USP) family, a group of DUBs for which very little is known in skeletal muscle . USPs are a highly diverse family with a common Ubiquitin carboxyl-terminal hydrolase (UCH 2_3) functional domain. Several USPs members contain a variety of other functional domains including dual specific phosphatases (DUSP), zinc-finger in ubiquitin carboxyl-terminal hydrolase (ZF-UBP), ubiquitin associated domains (UBA) or even an Hsp90 binding motif (MEEVD) (Additional file 11). A total of 45 USPs, representing almost the complete repertoire identified zebrafish, were found in the pacu transcriptome. Again, the normalized number of reads mapped was used to estimate their abundance in fast and slow skeletal muscle (Figure 4). We found that 4 USPs were significantly more abundant in slow (usp2b, usp10, usp43a usp48) (Figure 4A) and 9 USPs were more highly expressed in fast skeletal muscle (usp5b, usp9, usp14, usp19, usp21, usp24, usp28, usp45, usp47) (Figure 4B). The lack of studies concerning USP function in fish makes it difficult to explain these differences in a physiological context. With the objective of gaining knowledge about the role of USPs in muscle growth we performed a nutritional challenge experiment to analyse their expression during the transition from a catabolic to an anabolic state. To this end pacu juveniles were fasted for 4 days followed by a short period of satiation feeding (24 h) to stimulate protein synthesis and a subset of muscle expressed USP genes were measured by real-time PCR (GE) (Figure 4). The sub-set was chosen to reflect their diversity with respect to functional domains present (usp2a, usp2b, usp4, usp5a, usp5b, usp8, usp9, usp11, usp12a, usp12b, usp14, usp16, usp19, usp21, usp28, usp30, usp36 and usp46). The expression of several Pi3k/Akt pathway components (igf1, igf2a, igf2b, igf1ra, igf1rb and igf3) and E3-ubiquiting ligases (mafbx, murf1a, murf1b, fbox-25, huwe, ufd2, trip12, syyna and syah) were also measured to gain a better understanding of the metabolic context in which USP expression was occurring (Figure 5A). Transcript abundance from the transcriptome was highly correlated with that determined by qPCR (R2 = 0.76; n = 21; P < 0.0001; data not show).
The E3-ubiquitin ligases are an essential part of the proteasome system, directly involved in protein degradation [45,46]. All E3-ubiquitin ligases increased their relative abundance during fasting, for example, fbox-32 a muscle specific E3-ubiquitin ligase increased 100-fold (Figure 5B) and fbox25 increased 10-fold, in line with results from other similar studies [21,47]. The majority of the USPs increased their expression around 2-fold with fasting and recovered pre-fasting levels 24 h after re-feeding (Figure 5A). In contrast, Usp12a, showed a transient increase 6 h after re-feeding, and usp46 and usp5b did not change in expression (Figure 5D and E). The overlapping expression profiles of USPs and E3-ubiquitin ligases suggests that USPs play an important role during muscle atrophy. USPs can cleave ubiquitin from proteins, effectively removing the proteasome signal . It is possible that the increase in their abundance in fasted fish is related to fine tuning of the regulation of protein degradation. For example, many USPs targets are essential to maintain cell homeostasis including mdm4, p53, h2a, h2b, fbw7, fancd2 or brca2 [17,48-50], and it is possible that these proteins may be relatively spared during fasting.
We have produced an in depth transcriptome for fast and slow myotomal muscle for the pacu (Piaractus mesopotamicus), an important South American aquaculture species. This resource allowed us to characterise the expression signatures of the main myotomal muscle fibre types and identify candidate microsatellite sequences that could be used in breeding programs. The availability of the transcriptome allowed us to identify and analyse the expression of teleost-specific paralogues retained in the Ostariophysi lineage. The transcriptome also enabled a comprehensive study of E3 ubiquitin ligase and USP gene expression in the context of the transition between the fasting (catabolic state) and satiation feeding (anabolic state). We identified differences in expression within gene family members thereby identifying candidates for further investigation.
For the generation of de novo transcriptome fast and slow skeletal muscles were dissected from 5 adult pacu (Piaractus mesopotamicus) (1.50 ± 0.61 kg; mean ± SD body mass). Fish were maintained in 1000 litres fibreglass tanks at the Aquaculture Centre of the University of West Paulista (Unoeste) Presidente Prudente, São Paulo, Brazil, under natural photoperiod (12 L: 12D) and temperature (28°C ± 1°C, range). Fast skeletal muscle was dissected from the dorsal epaxial region at 0.5 fork length (FL) (FL, tip of snout to fork in the tail) whereas slow skeletal muscle was dissected from the lateral line and any remains of fast skeletal muscle carefully removed under a dissection microscope to obtain pure slow muscle. Tissues were preserve in RNAlater (Ambion/Applied Biosystems, Oslo, Norway) and frozen at −20°C until further analysis.
For the fasting-re-feeding experiment, 15 g pacu (n = 80) were maintained in duplicate fibreglass 500 litre tanks as described above and fed with a commercial diet until the start the experiment. Fish were fasted for 4 days followed by a period of satiation feeding for 24 h. Fast skeletal muscle was sampled before fasting (−4d), daily during food deprivation (−3d, −2d, −1d, 0d; n = 8) and 6, 12 and 24 h (n = 8) after re-feeding. Dissected fast skeletal muscle was preserved in RNAlater at −20°C until further analysis. All fish were sacrificed according to the Ethical Principles In Animal Research adopted by Brazilian College of Animal Experimentation (COBEA) and was approved by the Ethics Committee on Use of Animals/ Bioscience Institute/Unesp (CEUA = 506).
Samples sequencing and de novo assembly
Total RNA from adult pacu fast and slow skeletal muscle was used to prepare 10 individual Illumina libraries. Libraries preparation and sequencing was performed at the Centre for Applied Genomics of the Hospital for Sick Children (SickKids), Toronto, Canada. The resulting libraries were paired end sequenced using in an Illumina HiSeq2000.
Raw paired end reads generated were processed by the Department of Informatics of the Centre for Applied Genomics of SickKids hospital. After removing low quality reads, 86% of the paired end reads were de novo assembled using Trinity software . RSEM application was used to identify transcript abundance by estimating the number of reads mapped per contig. The DEseq algorithm from the Bioconductor/R packages was used to identify differentially expressed transcripts .
Contigs were annotated using Blast2GO software . Sequences were blasted against the NCBI non-redundant (nr) database using BLASTx with an e-value cut-off of 10−3 followed by functional annotation using software default parameters . Contigs were mapped against the known vertebrate metabolic and molecular pathways using the online KEGG Automatic Annotation Server (KAAS) . KAAS annotation was performed using the single-directional Best Hit (SBH) method against Homo sapiens, Pan troglodites, Mus musculus, Rattus norvegicus, Sus scrofa, Gallus gallus, Meleagris gallopavo, Danio rerio and Xenopus laevis.
Identification of complete coding sequence
Annotated contigs were blasted against the Zebrafish complete proteome  using tBLASTn algorithm in BioEdit software . BLAST alignments were explored to evaluate the percentage of conding sequence cover by the contig compared with its zebrafish orthologue. Sequences with more than 90% of coverage were considered as complete coding sequence (CDS). The CDS amino acids sequence was predicted using the Virtual Ribosome server .
Sequences successfully annotated covering >90% of the CDS were investigated for SSR using msatcommander-1.0.2-alpha .
RNA Extraction and cDNA synthesis
Total RNA was extracted using 1 ml TRIsure (Bioline, London, UK) following the manufacturer’s recommendations. Integrity was confirmed by ethidium bromide gel electrophoresis of 1 μg of total RNA. Concentration, 260/280 and 260/230 ratios were estimated using a NanoDrop 1000 spectrophotometer (Thermo Fischer Scientific, Waltman, MA). All RNA samples used had 260/280 nm and 260/230 ratios above 1.9 and 2.2 respectively. 1 μg of total RNA was reverse transcribed into cDNA for 30 min at 42°C using a Quantitect (QIAGEN, Manchester, UK) reverse transcription kit following manufacture’s recommendations including a genomic DNA wipe-out step. To ensure that no genomic DNA was present in the samples a RT- control without the reverse transcriptase enzyme was performed simultaneously.
Quantitative Real-time PCR
The following procedures were compliant with the minimal information requirements for publication of quantitative PCR guidelines . Primers were designed to have a Tm of 60°C using Net primer online server (Premier BioSoft). Where possible primers were designed to cross exon-exon junctions. Exon-exon junctions were predicted by aligning the pacu contig against their zebrafish orthologue complete gene sequence retrieved from Ensembl  using Spidey online server . Primers pairs, amplicon size and efficiency are listed in Additional file 12.
Quantitative PCR (qPCR) was performed using a MX3005P qPCR machine (Agilent, La Jolla, CA, USA). Each qPCR reaction contained 7.5 μl of SensiFast (Bioline) Master Mix, 6 μl cDNA (80-fold dilution and 40-fold dilution for igf genes) and 0.75 μl of each primer at 500 nM to a final volume of 15 μl. Duplicate reactions were performed in 96-well plates (Agilent) with the following protocol: initial activation 95°C for 2 min followed by 40 cycles of 5 s at 95°C, 20s at 65°C. The qPCR was followed by a dissociation-melting curve from 60 to 95°C to confirm that a single product was amplified. Control reactions included no-template and RT- were simultaneously amplified.
Ribosomal protein 13 and 19 (rpl13, rpl19), Peptidylprolyl isomerase Aa (ppiaa), elongation factor 1 alpha, Glyceraldehyde-3-phosphate dehydrogenase (gapdh) and hypoxanthine phosphoribosyltransferase 1 (hprt1) were tested as reference genes. BestKeeper  analysis showed that hprt1 was the most stable reference gene and was used for data normalisation. Relative expression was calculated using the Pfalff method .
Teleost specific paralogues identification
Contigs were blasted (BLASTx) against the Zebrafish proteome using BioEdit software with an e-value threshold of e−40 .
To confirm that contigs found were truly paralogues amino acids sequences from potential pacu paralogues were blasted (BLASTp) against the zebrafish (Danio rerio), stickleback (Gasterosteus aculeatus), takifugu (Takifugu rubripes), medaka (Oryzias latipes), green pufferfish (Tetraodon nigroviridis), chicken (Gallus gallus), frog (Xenopus laevis) and human (Homo sapiens) proteomes using Ensembl BLAST server . Best hits amino acids sequences from each proteome were retrieved. Sequences were aligned using the MAFFT online server . Phylogenetic trees were constructed using a Maximum Likelihood analysis using PhyML online server combined with the G-Blocks option to cure unreliable aligned sections . For each case the best evolutionary model was estimated using MEGA5 .
Global DGE statistic analysis was performed using DEseq package from R-Bioconductor [52,65]. For testing specific hypothesis involving differential mapping of specific pathways gene expression significance was tested using t-test, or Mann–Whitney U test when parametric were not fulfilled, followed by a Benjamin-Hochberg correction (False Discovery Rate, FDR). Differences between time-points in qPCR expression during fasting re-feeding were tested using Kruskal-Wallis test. Differences were considered significant when FDR, for differential mapping, or p-value, for qPCR expression, were <0.05. Gene expression data was clustered using an unsupervised hierarchical clustering algorithm using PermutMatrix .
Sänger A, Stoiber W. Muscle fiber diversity and plasticity. Fish Physiol. 2001;18:187–250.
Altringham JD, Johnston IA. Modelling muscle power output in swimming fish. J Exp Biol. 1990;148:395–402.
Coughlin DJ, Rome LC. The Roles of Pink and Red Muscle in Powering Steady Swimming in Scup, stenotomus chrysops. Am Zool. 1996;36:666–77.
Johnston IA, Calvo J, Guderley H, Fernandez D, Palmer L. Latitudinal variation in the abundance and oxidative capacities of muscle mitochondria in perciform fishes. J Exp Biol. 1998;201:1–12.
Wakeling JM, Johnston IA. Muscle power output limits fast-start performance in fish. J Exp Biol. 1998;201(Pt 10):1505–26.
Mascarello F, Rowlerson A, Radaelli G, Scapolo PA, Veggetti A. Differentiation and growth of muscle in the fish Sparus aurata (L): I. Myosin expression and organization of fibre types in lateral muscle from hatching to adult. J Muscle Res Cell Motil. 1995;16:213–22.
Johnston IA, Dvison W, Goldspink G. Energy metabolism of carp swimming muscle. J Comp Physiol. 1977;144:203–16.
Johnston IA, Bower NI, Macqueen DJ. Growth and the regulation of myotomal muscle mass in teleost fish. J Exp Biol. 2011;214:1617–28.
Cleveland BM, Weber GM, Blemings KP, Silverstein JT. Insulin-like growth factor-I and genetic effects on indexes of protein degradation in response to feed deprivation in rainbow trout (Oncorhynchus mykiss). Am J Physiol Regul Integr Comp Physiol. 2009;297:R1332–42.
Rescan PY, Montfort J, Fautrel A, Rallière C, Lebret V. Gene expression profiling of the hyperplastic growth zones of the late trout embryo myotome using laser capture microdissection and microarray analysis. BMC Genomics. 2013;14:173.
Johnston IA. Environment and plasticity of myogenesis in teleost fish. J Exp Biol. 2006;209:2249–64.
Johnston IA. Muscle development and growth: potential implications for flesh quality in fish. Aquaculture. 1999;177:99–115.
Tacchi L, Bickerdike R, Secombes CJ, Martin SA. Muscle-specific RING finger (MuRF) cDNAs in Atlantic salmon (Salmo salar) and their role as regulators of muscle protein degradation. Mar Biotechnol. 2012;14:35–45.
Salmerón C, García de la Serrana D, Jiménez-Amilburu V, Fontanillas R, Navarro I, Johnston IA, et al. Characterisation and Expression of Calpain Family Members in Relation to Nutritional Status, Diet Composition and Flesh Texture in Gilthead Sea Bream (Sparus aurata). PLoS One. 2013;8:e75349.
Argilés JM, López-Soriano FJ, Busquets S. Apoptosis signalling is essential and precedes protein degradation in wasting skeletal muscle during catabolic conditions. Int J Biochem Cell Biol. 2008;40:1674–8.
Lecker SH, Goldberg AL, Mitch WE. Protein degradation by the ubiquitin-proteasome pathway in normal and disease states. J Am Soc Nephrol. 2006;17:1807–19.
David Y, Ternette N, Edelmann MJ, Ziv T, Gayer B, Sertchook R, et al. E3 ligases determine ubiquitination site and conjugate type by enforcing specificity on E2 enzymes. J Biol Chem. 2011;286:44104–15.
Komander D, Clague MJ, Urbé S. Breaking the chains: structure and function of the deubiquitinases. Nat Rev Mol Cell Biol. 2009;10:550–63.
Daviet L, Colland F. Targeting ubiquitin specific proteases for drug discovery. Biochimie. 2008;90:270–83.
Amaral IPG, Johnston IA. Insulin-like growth factor (IGF) signalling and genome-wide transcriptional regulation in fast muscle of zebrafish following a single-satiating meal. J Exp Biol. 2011;214:2125–39.
Bower N, Johnston IA. Characterisation and differential regulation of MAFbx/Atrogin-1 a and b transcripts in skeletal muscle of Atlantic salmon (Salmo salar). Physiol Genomics. 2010;396:114–30.
Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, et al. Ensembl 2015. Nucleic Acids Res. 2014;43:D662–9. October 2014.
Garcia de la Serrana D, Mareco E, Johnston IA. Systematic variation in the pattern of gene paralogue retention between the teleost super-orders Ostariophysi and Acanthopterygii. Genome Biol Evol. 2014;6:981–7.
Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, et al. Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004;431:946–57.
Sato Y, Nishida M. Teleost fish with specific genome duplication as unique models of vertebrate evolution. Environ Biol Fishes. 2010;88:169–88.
Braasch I, Postlethwait J. Polyploidy in fish and the teleost genome duplication. In: Soltis PS, Soltis DE, editors. Polyploidy genome Evol. Berlin Heidelberg: Springer; 2012. p. 341–83.
Macqueen DJ, Johnston IA. A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification. Proc Biol Sci. 2014;281:20132881.
Amaral IPG, Johnston IA. Experimental selection for body size at age modifies early life-history traits and muscle gene expression in adult zebrafish. J Exp Biol. 2012;215:3895–904.
Magnanou E, Klopp C, Noirot C, Besseau L, Falcón J. Generation and characterization of the sea bass Dicentrarchus labrax brain and liver transcriptomes. Gene. 2014;544:56–66.
Garcia de la serrana Castillo D, Estévez A, Andree K, Johnston I. Fast skeletal muscle transcriptome of the Gilthead sea bream (Sparus aurata) determined by next generation sequencing. BMC Genomics. 2012;13:181.
Mojekwu TO, Anumudu CI. Microsatellite markers in Aquaculture: Application in Fish population genetics. IOSR J Environ Sci Toxicol Food Technol. 2013;5:43–8.
Calduch-Giner JA, Bermejo-Nogales A, Benedito-Palos L, Estensoro I, Ballester-Lozano G, Sitjà-Bobadilla A, et al. Deep sequencing for de novo construction of a marine fish (Sparus aurata) transcriptome database with a large coverage of protein-coding transcripts. BMC Genomics. 2013;14:178.
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.
Ganassi M, Badodi S, Polacchini A, Baruffaldi F, Battini R, Hughes SM, et al. Distinct functions of alternatively spliced isoforms encoded by zebrafish mef2ca and mef2cb. Biochim Biophys Acta. 1839;2014:559–70.
Garcia De La Serrana D, Johnston IA. Expression of heat shock protein (Hsp90) paralogues is regulated by amino acids in skeletal muscle of Atlantic salmon. PLoS One. 2013;8:e74295.
Maere S, Peer YV. Duplicate Retention After Small and Large-Scale Duplications. In: Dittmar K, Liberles D, editors. Evol after Gene Duplic. 1st ed. Hoboken, NJ, USA: John Wiley & Sons, Inc; 2010. p. 31–56.
Rommel C, Bodine SC, Clarke BA, Rossman R, Nunez L, Stitt TN, et al. Mediation of IGF-1-induced skeletal myotube hypertrophy by PI(3)K/Akt/mTOR and PI(3)K/Akt/GSK3 pathways. Nat Cell Biol. 2001;3:1009–13.
Glass DJ. Molecular mechanisms modulating muscle mass. Trends Mol Med. 2003;9:344–50.
Lang CH, Vary TC, Frost RA. Acute in vivo elevation of insulin-like growth factor (IGF) binding protein-1 decreases plasma free IGF-I and muscle protein synthesis. Endocrinology. 2003;144:3922–33.
Huang K, Fingar DC. Growing knowledge of the mTOR signaling network. Semin Cell Dev Biol. 2014;36C:79–90.
Bower NI, Garcia de la Serrana D, Cole NJ, Hollway GE, Lee HT, Assinder S, et al. STAC3 is required for myotube formation and myogenic differentiation in vertebrate skeletal muscle. J Biol Chem. 2012;287:43936–49.
Nelson BR, Wu F, Liu Y, Anderson DM, McAnally J, Lin W, et al. Skeletal muscle-specific T-tubule protein STAC3 mediates voltage–induced Ca2+ release and contractility. Proc Natl Acad Sci U S A. 2013;110:11881–6.
Fleming JR, Crockford T, Altringham JD, Johnston IA. Effects of termperature acclimatation on muscle relaxation in the carp: A mechanical, biochemical and ultrastructural study. J Exp Zool. 1990;255:286–95.
Gabillard JC, Biga PR, Rescan PY, Seiliez I. Revisiting the paradigm of myostatin in vertebrates: insights from fishes. Gen Comp Endocrinol. 2013;194:45–54.
Deshaies RJ, Joazeiro CA. RING domain E3 ubiquitin ligases. Annu Rev Biochem. 2009;78:399–434.
Berndsen CE, Wolberger C. New insights into ubiquitin E3 ligase mechanism. Nat Struct Mol Biol. 2014;21:301–7.
Bower NI, Taylor RG, Johnston IA. Phasing of muscle gene expression with fasting-induced recovery growth in Atlantic salmon. Front Zool. 2009;6:18.
Gomes LC, Di Benedetto G, Scorrano L. During autophagy mitochondria elongate, are spared from degradation and sustain cell viability. Nat Cell Biol. 2011;13:589–98.
Cataldo F, Peche LY, Klaric E, Brancolini C, Myers MP, Demarchi F, et al. CAPNS1 regulates USP1 stability and maintenance of genome integrity. Mol Cell Biol. 2013;33:2485–96.
Holowaty MN, Sheng Y, Nguyen T, Arrowsmith C, Frappier L. Protein interaction domains of the ubiquitin-specific protease, USP7/HAUSP. J Biol Chem. 2003;278:47753–61.
Haas B, Papanicolaou A, Yassour M. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.
Gentleman RC, Gentleman RC, Carey VJ, Carey VJ, Bates DM, Bates DM, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008;36:3420–35.
Moriya Y, Itoh M, Okuda S. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35:W182–5. Web Server issue.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999:95–98.
Wernersson R. Virtual Ribosome: a comprehensive DNA translation tool with support for integration of sequence feature annotation. Nucleic Acids Res. 2006;34:W385–8. Web Server issue.
Faircloth BC. Msatcommander: Detection of Microsatellite Repeat Arrays and Automated, Locus-Specific Primer Design. Mol Ecol Resour. 2008;8:92–4.
Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611–22.
Wheelan SJ, Church DM, Ostell JM. Spidey : A Tool for mRNA-to-Genomic Alignments. Genome Res. 2001;11:1952–7.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26:509–15.
Katoh K, Standley D. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Dereeper A, Guignon V, Blanc G. Phylogeny. fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008;36:W465–9. Web Server issue.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Caraux G, Pinloche S. PermutMatrix: a graphical environment to arrange gene expression profiles in optimal linear order. Bioinformatics. 2005;21:1280–1.
MDPS and EAM were supported by FAPESP, Proc. n° 12/02489-4 and 2011/09346-1. This work also received funding from CAPES, Proc. n° 2524/12. IAJ and DGDLS were supported by the Marine Alliance for Science and Technology for Scotland pooling initiative and Scottish Funding Council grant number HR09011.
The authors declare that they have no competing interests.
MDPS and IAJ conceived the study, DGDLS was responsible for RNA extraction, dsDNA synthesis, and bioinformatics analysis, EAM was responsible for fish husbandry, bioinformatics analysis, and gene expression analysis. All authors were involved in writing the manuscript and approved the final version.
Annotated transcriptome. Annotated nucleotide using blast2go from the pacu slow and fast muscle transcriptome.
KEEG Automatic Annotation Sever results. Hierarchical information of the sequences from annotated transcriptome mapped to the KEGG pathways maps using the automatic annotation tool KAAS .
Simple Sequence Repeats Identification. Simple Sequence Repeats identified in contigs containing over 90% of the CDS using Msatcommander software .
Teleost-specific paralogues sequences and phylogenetic trees. Supplementary file contains the amino acids sequences for the pacu paralogues and the corresponding orthologues for zebrafish, stickleback, tilapia, green puffer fish and human, their MAFFT alignments and individual phylogenetic trees.
Fast and slow skeletal muscle global digital gene expression analysis. Global digital gene expression analysis between fast and slow skeletal muscle was performed using the number of reads mapped normalized by length and library size was investigated using R-Bioconductor. Normalized counts for individual animals per contig are shown. Differences between tissues were analysed using the DESEQ package. Significant differences were considered when FDR < 0.05.
Fast and slow skeletal muscle Gene Ontology (GO) analysis. GO analysis was performed for differently expressed contigs originating from the fast and slow skeletal muscle comparison. Significantly enriched categories were investigated using Fisher Exact Test with the whole annotated transcriptome as reference and considered significant when FDR < 0.05.
Digital gene expression of metabolic pathways in fast and slow skeletal muscle. Components from the different metabolic pathways were extracted from KAAS annotation results. Analysis between fast and slow skeletal muscle was performed using the number of reads mapped normalized by length and library size. Normalized counts for individual animals per contig are shown. Differences between tissues were analysed by t-test followed by a False Discovery Rate correction. Significant differences were considered when FDR < 0.05.
Schematic representation of fast and slow skeletal muscle main metabolic pathways. Metabolic genes represented in the transcriptome mapped into a reconstruction of the main metabolic pathways. “Red circles” and “Empty circles” indicate components significantly higher in slow and fast skeletal muscle respectively muscle (FDR < 0.05). “Green circles” indicate components with no significant differences between fibre types.
Teleost-specific paralogues digital gene expression. Teleost-specific paralogues digital gene expression in slow and fast skeletal muscle. Comparisons were made between paralogue pairs and fibres types. Analysis between fast and slow skeletal muscle was performed using the number of reads mapped normalized by length and library size. Normalized counts for individual animals per contig are shown. Differences between tissues were analysed by t-test followed by a False Discovery Rate correction. Significant differences were considered when FDR < 0.05.
Pi3k/mTor pathway components digital gene expression. Pi3k/mTor components digital gene expression for slow and fast skeletal muscle. Analysis between fast and slow skeletal muscle was performed using the number of reads mapped normalized by length and library size. Normalized counts for individual animals per contig are shown. Differences between tissues were analysed by t-test followed by a False Discovery Rate correction. Significant differences were considered when FDR < 0.05.
Pacu ubiquitin specific proteases functional domains and digital gene expression. Ubiquitin specific proteases (USP) digital gene expression analysis between slow and fast skeletal muscle. Analysis between fast and slow skeletal muscle was performed using the number of reads mapped normalized by length and library size. Normalized counts for individual animals per contig are shown. Differences between tissues were analysed by t-test followed by a False Discovery Rate correction. Significant differences were considered when FDR < 0.05. Pacu, zebrafish and human protein USP functional domains were determined using InterProScan webserver. Relative gene expression was analysed for fbox25, huwe, mafbx, murf1a, murf1b, syah1, syyna, trip12, ufd2, igf1, igf2a, igf2b, igf3, igf1ra, igf1rb, usp2a, usp2b, usp4, usp5a, usp5b, usp8, usp9, usp11, usp12b, usp12b, usp14, usp16, usp19, usp21, usp24, usp28, usp30, usp36 and usp46, Values represents mean ± SE (n = 8 fish).
qPCR primers sequence of target genes: fbox25, huwe, mafbx, murf1a, murf1b, syah1, syyna, trip12, ufd2, igf1, igf2a, igf2b, igf3, igf1ra, igf1rb, usp2a, usp2b, usp4, usp5a, usp5b, usp8, usp9, usp11, usp12b, usp12b, usp14, usp16, usp19, usp21, usp24, usp28, usp30, usp36 and usp46 and housekeeping genes: Ef1a, Polr2a, Gapdh, Hprt1, ppiaa, sdha, Rpl13, Rpl19.