- Research article
- Open Access
Characterisation and expression analysis of cathepsins and ubiquitin-proteasome genes in gilthead sea bream (Sparus aurata) skeletal muscle
BMC Research Notesvolume 8, Article number: 149 (2015)
The proteolytic enzymes involved in normal protein turnover in fish muscle are also responsible for post-mortem softening of the flesh and are therefore potential determinants of product quality. The main enzyme systems involved are calpains, cathepsins, and the ubiquitin-proteasome (UbP). In this study on Sparus aurata (Sa), the coding sequences of cathepsins (SaCTSB and SaCTSDb) and UbP family members (SaN3 and SaUb) were cloned from fast skeletal muscle, and their expression patterns were examined during ontogeny and in a fasting/re-feeding experiment.
The amino acid sequences identified shared 66-100% overall identity with their orthologues in other vertebrates, with well conserved characteristic functional domains and catalytic residues. SaCTSDb showed phylogenetic, sequence and tissue distribution differences with respect to its paralogue SaCTSDa, previously identified in the ovary. Expression of gilthead sea bream cathepsins (B, L, Da, Db) and UbP members (N3, Ub, MuRF1 and MAFbx) in fast skeletal muscle was determined at three different life-history stages and in response to fasting and re-feeding in juveniles. Most of the proteolytic genes analysed were significantly up-regulated during fasting, and down-regulated with re-feeding and, between the fingerling (15 g) and juvenile/adult stages (~50/500 g), consistent with a decrease in muscle proteolysis in both later contexts. In contrast, SaCTSDa and SaMuRF1 expression was relatively stable with ontogeny and SaUb had higher expression in fingerlings and adults than juveniles.
The data obtained in the present study suggest that cathepsins and UbP genes in gilthead sea bream are co-ordinately regulated during ontogeny to control muscle growth, and indicate that feeding regimes can modulate their expression, providing a potential dietary method of influencing post-mortem fillet tenderisation, and hence, product quality.
Gilthead sea bream (Sparus aurata L.) is widely farmed around the Mediterranean area  and because of its commercial importance, muscle growth regulation in this species has been the subject of recent research [2-9]. In common with most teleosts, gilthead sea bream exhibits indeterminate growth, with muscle mass increasing by hyperplasia (production of new fibres) until 40-50% of the maximum length, and by hypertrophy (increase of fibres size) until mortality or senescence occur [10,11].
Muscle growth represents the balance between protein synthesis and degradation, and sarcomeric components have a range of half-lives . In vertebrates, four catabolic systems are known to be involved in muscle proteolysis: a) the Ca2+-dependent proteinases (calpains), b) the autophagy-lysosome system (cathepsins), c) the ATP-dependent ubiquitin-proteasome (UbP) pathway and, d) the apoptosis protease system (caspases) [13-16]. Calpains are considered a system of primary protein degradation, with a regulatory or signalling function, since they do not cleave proteins to amino acids or small peptides . Cathepsins and the UbP pathway are required for the complete degradation of protein substrates .
The UbP pathway operates through a multi-subunit proteolytic complex, the proteasome, and it targets specific proteins for destruction through a three-step enzymatic process that covalently links a poly-ubiquitin (Ub) chain to the protein substrate to be degraded. Two E3 Ub ligases, MuRF1 and MAFbx (also known as Atrogin-1), are considered transcriptional markers involved in muscle wasting in vertebrates; however, they seem to contribute differently to muscle loss. MuRF1 is involved in the breakdown of myofibrillar proteins such as myosin heavy chain, whereas MAFbx participates in the control of protein synthesis, regulating transcription factors such as MyoD [19,20]. The UbP pathway in vertebrates is particularly important during muscle atrophy, either caused by starvation or wasting diseases [21-25], and it is also involved in the age-related loss of muscle mass (sarcopenia) in mammals [26,27]. In fish, prolonged fasting has been shown to induce muscle atrophy [28-30], whereas signs of sarcopenia as the fish age have been reported only for species with determinate growth such as the zebrafish (Danio rerio) . Previous studies in rainbow trout (Oncorhynchus mykiss) have shown that fasting increased 20S proteasome activity  whereas the mRNA expression of several UbP pathway members decreased with re-feeding [33,34]. In the same species, UbP gene expression was unchanged following muscle wasting associated with spawning [16,35].
Cathepsins are lysosomal proteases optimally active in a slightly acidic environment, and are classified as aspartic, serine or cysteine proteases according to the amino acid found in their active site . The role of cathepsins in muscle proteolysis and post-mortem degradation has been investigated in several fish species. Calpains and cathepsins are thought to be involved in enzymatic degradation of key structural and extracellular matrix proteins during post-mortem tenderisation . In sockeye salmon (Oncorhynchus nerka), the muscle protein degradation that occurs during the final stages of maturation following migration and fasting involves lysosomal cathepsin D, and to a lesser extent, cathepsin L . In rainbow trout, fasting and re-feeding respectively increased and decreased the mRNA expression of cathepsins B, D and S , and spawning-induced muscle wasting was associated with increased cathepsin L and D mRNA levels [16,35]. In gilthead sea bream, distinct changes in the expression and enzyme activities of cathepsins B, D and L have been observed in the maturing ovary, indicating a specific function for these enzymes during the follicle maturation process [39-41].
Following slaughter the endogenous activity of proteolytic enzymes contributes to post-mortem softening and a loss of product quality [37,42,43]. We have recently characterised several members of the calpain system in gilthead sea bream muscle, and demonstrated that their expression is modulated by nutritional status and diet composition . In the present study, complete and partial coding sequences (CDS) for the cathepsins (SaCTSB and SaCTSDb) and UbP family members (SaN3 and SaUb) were cloned from fast skeletal muscle, and their expression patterns examined during ontogeny and with fasting and re-feeding.
All animal handling procedures were approved by the Ethics and Animal Care Committee of the University of Barcelona (CEEA 239/09) and the “Departament de Medi Ambient i Habitatge” (DMAH permit number 5420, Generalitat de Catalunya, Spain) following the European Union, Spanish and Catalan Government established norms and procedures.
Fish and experimental trials
Fish used for the cloning and expression analysis during ontogeny were obtained from a fish farm in Northern Spain. Gilthead sea bream for the fasting/re-feeding experiment were obtained from the “Institut de Recerca i Tecnologia Agroalimentàries” (IRTA) facilities (Sant Carles de la Ràpita, Spain). All fish were acclimatized to the facilities at the University of Barcelona (Barcelona, Spain) for a minimum of two weeks prior to sampling or experimental manipulations, fed ad libitum twice daily with commercial pellets (Excel, Skretting, Burgos, Spain) and held at 21 ± 1°C (range) and pH 7.5-8 in 200 or 400 L recirculating seawater tanks with 12 h dark: 12 h light photoperiod. For cloning 10 juvenile gilthead sea bream of 43 ± 3 g and for the tissue screening 6 juveniles of 151 ± 12 g (mean ± SEM) were used. For the life-history stages study, groups of 5 fish each of 15 ± 1 g (fingerlings, FL), 47 ± 5 g (juveniles, JV) and 503 ± 37 g (adults, AD) (mean ± SEM) were used. For the fasting/re-feeding experiment 120 juvenile gilthead sea bream (50 ± 6 g) were used, and the trial performed as previously described  (mean ± SEM). First, acclimatized fish were sampled for time 0 (D0C) and then, divided into two conditions: control fed group at 3% (m/m d−1) body weight (C) and a fasted group (F). Samples were collected at days 15 and 30 (D15C/F and D30C/F) and subsequently, the remaining fasted animals (D0R) were re-fed at 2% (m/m d−1) body weight (lower ration than the control to facilitate correct adaptation of the digestive system) and sampled at days 7 and 14 (D7R and D14R). Before sampling all fish were fasted 24 h, anesthetized with tricaine methane sulphonate (MS-222 0.1 g/L, Sigma, Tres Cantos, Spain) and sacrificed with a blow on the head and medullar section. Samples of fast skeletal muscle (FM), slow skeletal muscle (SM), heart (HE), brain (BR), liver (LI), spleen (SP) and immature gonad (GO) were taken by sterile dissection and immediately snap-frozen in liquid nitrogen and stored at −80°C for cloning or gene expression analyses.
RNA extraction and cDNA synthesis
Total RNA was extracted from 20–130 mg of tissue following the guanidinium thiocyanate-phenol-chloroform method  using TRIreagent (Applied Biosystems, Alcobendas, Spain), quantified using a NanoDrop2000 spectrophotometer (Thermo Scientific, Alcobendas, Spain) and, RNA quality was analysed by 1% (m/v) agarose gel electrophoresis. One μg of total RNA per sample was used to synthesise first-strand cDNA using the Transcriptor First Strand cDNA Synthesis Kit (Roche, Sant Cugat del Vallès, Spain) or the AffinityScript QPCR cDNA Synthesis Kit (Agilent Technologies, Las Rozas, Spain) in the case of the tissue screening, following the manufacturers’ instructions. cDNA samples were diluted 1:5 in milliQ H2O for conventional polymerase chain reaction (PCR) and diluted 1:40 to 1:100 in milliQ H2O for real-time quantitative PCR (qPCR).
Cloning and sequencing
To obtain the complete sequences of gilthead sea bream (Sa) cathepsin B (SaCTSB), the new paralogue of cathepsin D (SaCTSDb), the proteasome beta type-4 subunit (SaN3, also known as PSMB4) and the ubiquitin (SaUb) from fast skeletal muscle, specific primers for PCR were designed using the DNAman software package (Lynnon, Quebec, Canada) and Net primer (Premier BioSoft) (www.premierbiosoft.com/netprimer/) using the gilthead sea bream ESTs (Expressed Sequence Tags) from the NCBI database (SaCTSB: [GenBank: HS985610] and [GenBank: FG26781]; SaCTSDb: [GenBank: FM146030] and [GenBank: FG26194]; SaN3: [GenBank: HS988518]; SaUb: [GenBank: AM955423]) as templates (Table 1). The cloning was performed as previously described . Briefly, PCR products were separated by gel electrophoresis and purified from the agarose gel using the PureLink Quick Gel Extraction Kit (Invitrogen, Alcobendas, Spain). The purified PCR product was ligated into T/A pCR4-TOPO vector and transformed into chemically competent TOP10 Escherichia coli cells (all from Invitrogen, Alcobendas, Spain). One to three clones of each PCR product were sequenced using BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Alcobendas, Spain) and analysed at the “Serveis Cientificotècnics” of the University of Barcelona (Barcelona, Spain). Sequenced products were joined in silico using the sequence alignment editor and sequence analysis program BioEdit  to produce contigs with a single open reading frame (ORF). Sequences generated were analysed for similarity with other known sequences using the BLAST programs (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
Sequence and phylogenetic analyses
The putative protein architecture (domains, active sites and other important motifs) from the sequences generated was determined according to the literature and the conserved domain search program of NCBI  and the simple modular architecture research tool (SMART) version 4.0 (http://smart.embl-heidelberg.de) . Compute pI/Mw tool (ExPASy, Switzerland, http://www.expasy.org/tools/pi_tool.html) was used to estimate the molecular weight (Mw) of the predicted proteins and NetNGlyc 1.0 Server (http://www.cbs.dtu.dk/services/NetNGlyc/) was used to predict the N-Glycosylation sites. All the alignments were created with MAFFT version 7.058b (http://mafft.cbrc.jp/alignment/software/) and G-INS-i (recommended for <200 sequences with global homology) strategy. Sequences used in the present study other than those cloned from gilthead sea bream fast skeletal muscle were obtained from NCBI. Phylogenetic analyses, including multiple sequence alignment and Maximum Parsimony (MP) tree prediction of 27 vertebrate cathepsin orthologs of CTSB, CTSD and CTSL were conducted using MEGA version 5.0 . MP trees were obtained using the Close-Neighbour-Interchange algorithm and bootstrap values were inferred from 1000 replicates. The gilthead sea bream calpain 1 [GenBank: KF444899], a cytosolic cysteine protease, was used to root the phylogenetic tree.
The mRNA levels of both SaCTSD paralogues under normal physiological conditions were measured in gilthead sea bream tissues using qualitative PCR and elongation factor 1-α (SaEF1α) was used as a loading control. Reactions were performed in a final volume of 50 μL, containing 1 μL of first-strand cDNA (equivalent to 4 ng of reverse transcribed total RNA), 1.5 U of Taq polymerase (Sigma, Tres Cantos, Spain) and 200 nM (final concentration) of sense and antisense primers (Table 2). Reactions proceeded in a MyiQ Thermal Cycler (Bio-Rad, El Prat de Llobregat, Spain) with the following protocol: 1 cycle at 95°C for 5 min, 35 cycles at 95°C for 30 s, 56-60°C (primer dependent, see Table 2) for 30 s, 72°C for 30 s and 1 cycle at 72°C for 7 min. Each reaction product was separated by agarose gel electrophoresis and visualised using SYBR Safe DNA gel stain (Life Technologies, Alcobendas, Spain) in a LAS-3000 (Fujifilm, Madrid, Spain) to confirm that a single product was amplified, and then sequenced to confirm the specificity of each assay as explained above in 2.4. Semi-quantification of SaCTSD paralogues relative expression normalized to SaEF1α was performed determining band intensity using the ImageJ software version 1.47 (National Institutes of Health, Bethesda, USA).
Quantitative real-time PCR (qPCR)
To characterise the two cathepsin D paralogues in gilthead sea bream, mRNA abundance in fast skeletal muscle and immature gonad from the tissue screening samples was evaluated using qPCR. mRNA levels of the different cathepsins and UbP genes in fast skeletal muscle from different life-history stages and the fasting/re-feeding experiment were also determined. In addition, Saβ-actin, the 18S ribosomal RNA (Sa18S) and SaEF1α were tested as reference genes. The specific gilthead sea bream qPCR primers not previously reported (see Table 2) were designed using the DNAman software package and Net primer (Premier BioSoft) (www.premierbiosoft.com/netprimer/). The qPCR assay was conducted as previously described . Briefly, reactions contained first-strand cDNA (equivalent to 2.5 ng of reverse transcribed total RNA), iQ SYBR Green Supermix (Bio-Rad, El Prat de Llobregat, Spain) and 250 nM (final concentration) of sense and antisense primers (Table 2), and were performed in triplicate using a MyiQ or a CFX384 thermocycler (Bio-Rad, El Prat de Llobregat, Spain). The protocol consisted on 1 cycle of 3 min at 95°C and 40 cycles of 10 s at 95°C and 30 s at 54-60°C (primer dependent, see Table 2), followed by an amplicon dissociation analysis from 55 to 95°C at 0.5°C increase each 30 s, where a single peak was observed confirming the specificity of the reaction. SYBR Green fluorescence was recorded during the annealing-extending phase of cycling. Expression results were normalized to Saβ-actin, the most stable of the three reference genes tested, and analysed by the delta-delta method .
Statistical analyses of all parameters were performed in PASW Statistics 17.0 (IBM, Chicago, USA). Normality was analysed according to the Shapiro-Wilk test and homogeneity in variance according to Levene’s test. Statistical differences among groups in the ontogeny as well as the fasting and re-feeding experiments were assessed by one-way ANOVA, followed by Tukey’s post hoc test, or between the two SaCTSD paralogues in the tissue expression analysis by Student’s t-test. A significance of p < 0.05 was applied to all statistical tests performed. Data are presented as mean ± standard error of the mean (SEM).
Sequence analysis of SaCTSB, SaCTSDb, SaN3 and SaUb
The CDS of the first cathepsin identified, SaCTSB, consisted of 993 base-pairs (bp) that encode a protein of 330 amino acids (aa) with a putative molecular mass (Mw) of 36.41 kDa [GenBank: KJ524457]. BLAST analysis showed that SaCTSB shares 70-82% overall sequence identity with the cathepsin B proteins of a number of fish, amphibian and mammalian species (Table 3). In silico analysis identified in SaCTSB an N-terminal signal peptide (I29/Propeptide C1), a cathepsin B propeptide region, and a papain family cysteine protease domain (Additional file 1: Figure S1). The protease domain contains the four essential residues for catalysis, i.e. Q101, C107, H277 and N297, which are highly conserved among vertebrates plus a predicted N-glycosylation site located at position 190. The CDS of the second cathepsin cloned, SaCTSDb, was 1191 bp encoding a 396 aa protein with a putative Mw of 42.98 kDa [GenBank: KJ524456]. SaCTSDb consists of a putative N-terminal signal peptide, a cathepsin D propeptide region (A1_Propeptide), and an aspartyl protease domain (Additional file 1: Figure S2), with the two aspartyl, D94 and D281 catalytic residues, conserved. A comparison of the SaCTSDb with other vertebrate cathepsin D protein sequences revealed a high degree of sequence similarity (66-88%; Table 3). Moreover, the sequence analysis of SaCTSDb predicted two possible N-glycosylation sites, located at positions 131 and 249 (Figure 1). Compared to the previously reported cathepsin D in gilthead sea bream ovary (SaCTSDa), the newly identified SaCTSDb has 57% aa identity, but lacks the third residue of N-glycosylation at position 337 (Figure 1). The partial CDS of SaN3 (97% of the putative molecule) was 756 bp and encoded 252 aa [GenBank: KJ524458]. SaN3 showed high levels of identity (74-93%) with other vertebrate proteasome N3/PSMB4 sequences (Table 3). The partial SaN3 protein sequenced contains the proteasome domain with threonine endopeptidase activity and shows the characteristic aa in the active site and interaction site typical of proteasome beta type-4 subunits (Additional file 1: Figure S3). Finally, we identified SaUb with a CDS of 918 bp, encoding a 305 aa protein with a putative Mw of 34.32 kDa [GenBank: KJ524459]. Four identical aa repeat units, termed R1 to R4, were present in the SaUb sequence each consisting of 216 bp coding for a 72 aa Ub monomer (Additional file 1: Figure S4). Each Ub monomer showed an interaction site with the Ub-conjugating enzyme (E2) and also an interaction site with the C-terminal hydrolase (UCH), as well as binding sites to the CUE domain of the Cue2 protein (Additional file 1: Figure S4). The protein sequence of SaUb was 100% identical to all vertebrate Ub protein sequences analysed (Table 3).
Phylogenetic analysis of SaCTSs
Phylogenetic analysis of cathepsins B, D and L from representative mammals, amphibians and fish, including SaCTSB, SaCTSDa, SaCTSDb and SaCTSL, produced an MP-phylogenetic tree that contained two distinct branches. The cysteine proteases, cathepsins B and L, clustered together separated from the aspartic proteases, the cathepsin Ds (Figure 2). The SaCTSB, SaCTSL, SaCTSDa and SaCTSDb were clustered together with their orthologous vertebrate cathepsins (Figure 2). The relationships revealed in the phylogenetic tree were in agreement with the concept of traditional taxonomy. In addition, the new cathepsin D paralogue cloned from gilthead sea bream fast skeletal muscle, the SaCTSDb, appeared phylogenetically related to its CTSD orthologues from other teleost and tetrapod species, whereas the other paralogue previously cloned in ovary, SaCTSDa, formed a clade with only CTSD sequences of teleost species (Figure 2).
Tissue expression of SaCTSDa and SaCTSDb
To further characterise the two cathepsin D paralogues identified, and to see if they have other differences (e.g. tissue distribution) that can help postulate differences in function and explain why they have both been retained throughout evolution, gene expression analyses of a tissue panel was performed.
Qualitative PCR analysis showed that under normal physiological conditions, both SaCTSDa and SaCTSDb were ubiquitously expressed in all tissues analysed, although SaCTSDb mRNA was relatively more abundant in fast and slow skeletal muscle (Figure 3A). Moreover, both paralogues showed lower levels of expression in the liver compared to the other tissues.
Quantitative analysis of the expression of both paralogues in gilthead sea bream fast skeletal muscle and immature gonad confirmed that SaCTSDb has 15-fold higher expression in muscle than SaCTSDa, whereas a similar level of expression for both paralogues was observed in the gonad (Figure 3B).
Expression of SaCTSs and SaUbP genes during life-history stages
To study the transcriptional regulation of the different proteolytic members present in the fast skeletal muscle of gilthead sea bream, we first analysed the expression of the SaCTSs and SaUbP members at three different life-history stages. qPCR analysis showed that the mRNA expression of SaCTSB, SaCTSL, SaCTSDb and SaMAFbx decreased significantly as fish grow, being greater in the muscle of fingerlings than in juveniles or adult fish (Figure 4A, B, D and H). SaN3 expression was also down-regulated with ontogeny, being significantly higher in fingerlings than in adult fish (Figure 4E). In contrast, SaUb expression was significantly higher in the muscles of both fingerlings and adult fish relative to juveniles (Figure 4F), whereas the expression of SaCTSDa and SaMuRF1 remained unchanged (Figure 4C and G).
Expression of SaCTSs and SaUbP genes in response to fasting and re-feeding
The transcriptional regulation of SaCTSs and SaUbP members was studied in fish skeletal muscle during fasting (15 and 30 days) and subsequent re-feeding for 7 and 14 days. The mRNA levels of all genes analysed were increased during the two time points of fasting except for SaCTSB that remained unchanged (Figure 5A). The gene expression of SaCTSL (Figure 5B), SaCTSDb (Figure 5D), SaN3 (Figure 5E), SaMuRF1 (Figure 5G) and SaMAFbx (Figure 5H) was significantly higher during the whole fasting period compared to the corresponding control condition; whereas SaCTSDa (Figure 5C) and SaUb (Figure 5F) increased at day 15 but decreased again after 30 days of fasting to values similar to those of the control. All the genes studied were down-regulated in response to re-feeding. SaCTSL (Figure 6B), SaCTSDb (Figure 6D), SaN3 (Figure 6E), SaUb (Figure 6F), SaMuRF1 (Figure 6G) and SaMAFbx (Figure 6H) decreased significantly after 7 days of re-feeding and remained low at day 14, while the expression of SaCTSB (Figure 6A) and SaCTSDa (Figure 6C) decreased significantly only after 14 days of re-feeding.
In the present study, the complete CDS of cathepsin B (SaCTSB), a new paralogue of cathepsin D (SaCTSDb), ubiquitin (SaUb), and partially the CDS of the proteasome subunit beta type-4 known as N3 or PSMB4 (SaN3), were cloned from gilthead sea bream fast skeletal muscle. The deduced amino acid sequences obtained, shared high levels of overall identity with orthologues of CTSB, CTSD, Ub and N3/PSMB4 in teleosts and other vertebrate species. For example, SaCTSB and SaCTSDb, possess an N-terminal signal peptide, a propeptide region, and a cysteine or aspartyl catalytic protease domain, respectively, with complete conservation at the catalytic residues [36,51]. SaN3 has a histidine residue replacing the mammalian arginine before the key threonine essential for proteolytic function  as previously reported in the N3 sequence of rainbow trout . In order to maintain the subunit in an inactive form until the proteasome is assembled the beta type subunits are synthesized as inactive precursors and activated after autocatalytic internal cleavage . Analysis of the SaN3 sequence also showed that the active site and the beta subunit interaction site were conserved as in other vertebrate N3 sequences . The SaUb cDNA identified had four 228 bp repeats each encoding a 72 amino acids Ub monomer, which has 100% identity with other vertebrate Ub sequences. The final repeat has an additional asparagine, rather than a tyrosine residue as in mammals , as has been shown to be the case also in Ub from rainbow trout . The structural features and the high levels of identity with the respective vertebrate sequences indicate that SaN3 and SaUb are the corresponding structural orthologues of N3 and Ub, respectively.
Phylogenetic analysis of cathepsin sequences produced a tree with two main branches; one with the cysteine CTSB and CTSL proteases and the other one with the aspartic CTSDs. The cysteine proteases arose early during evolution and the high percent identity among residues across taxa suggests that both cathepsin subfamilies evolved from gene duplication events . Interestingly, SaCTSDb seems to be the more common and less derived SaCTSD paralogue from the teleost whole genome duplication . Moreover, this new paralogue SaCTSDb clustered with other fish and mammalian cathepsin Ds, whereas SaCTSDa clustered on a separate branch with only fish members.
SaCTSDb lacks the third putative N-glycosylation site present in SaCTSDa whereas other vertebrate sequences possess only one or two, as previously reported . A previous study in human cells with mono-glycosylated and non-glycosylated cathepsin D mutants has shown that glycosylation is not necessary for folding or enzyme activity but it is required for targeting the enzyme to lysosomes . SaCTSDa and b were ubiquitously expressed in all tissues analysed as previously reported in other fish, avian and mammalian species [60-65]. The SaCTSDa gene was expressed in gonad as previously observed, in agreement with its role during the cleavage of vitellogenin to yield yolk proteins [39,62,66-70]. The SaCTSDb gene was more highly expressed in skeletal muscle tissue than SaCTSDa, consistent with it having probably some specific role in the degradation of muscle proteins, which should be addressed in future studies.
In the present study, the expression of most cathepsins and UbP members’ decreased with ontogeny. Dietary protein requirements decrease as the fish grow, e.g. the required protein and energy content in diets are higher in smaller (<100 g) than in larger fish . Therefore, it appears that mRNA levels reflect a decrease in proteolytic activity and protein turnover, where fingerlings have higher proteolysis rates than juvenile and adult fish, as well as higher rates of protein synthesis, as previously shown for example in rainbow trout [72,73]. Interestingly, expression of SaCTSDb and SaMAFbx decreased during ontogeny whereas SaCTSDa and SaMuRF1 expression was unchanged, indicating different transcriptional regulation and perhaps function between the cathepsin D paralogues and the E3 Ub ligases. MAFbx mRNA levels decreased in gilthead sea bream muscle with ontogeny, as previously observed in rat gastrocnemius muscle during the age-related loss of muscle mass . On the other hand, lack of transcriptional regulation of MuRF1 has been previously reported also in the human muscle of the elderly , in agreement with our results. Furthermore, SaUb mRNA expression was higher in fingerling and adult fish than juveniles. In mammals, a poor capacity to regenerate muscle with advanced age is thought to be due to impaired signalling, exhaustion of the satellite cell pool and/or changes in the extracellular matrix . Consistent with this observation, a previous study found that Ub protein expression was also up-regulated in rat and human fast muscle fibres during aging, and it was shown using C2C12 muscle cells that Ub suppresses proliferation, which may be associated with the poor healing potential in older individuals . The increase in SaUb expression between juveniles and adults may reflect the early stages of a similar age-related impairment of regenerative capacities in gilthead sea bream. However, Ub has been shown to have numerous proteolytic functions, including the regulation of proteasomal and lysosomal degradation, and also non-proteolytic functions, such as the regulation of protein interactions, activity and localization within the cell . Thus, the significance of increased SaUb expression in adult fish may well be related to other factors and needs to be further explored.
qPCR analyses revealed a general significant increase and decrease in the expression of all the studied SaCTSs and SaUbP genes with fasting and re-feeding, respectively, with the only exception of SaCTSB that was unaffected in response to fasting. In a recent study on gilthead sea bream using a similar experimental design, fasting had little effect on several catalytic and regulatory members of the calpain system whereas these genes were down-regulated with re-feeding . Thus there are likely generalised decreases in proteolytic enzyme expression in this species under anabolic conditions, assuming changes in transcript abundance are reflected in protein levels. In agreement with the present data, in a previous study in rainbow trout, cathepsin B, D and S and genes from the UbP pathway were up-regulated during fasting, and were also significantly down-regulated following re-feeding , and the same was observed with regards to Atlantic salmon (Salmo salar) cathepsin L , and the E3 ligases, MAFbx and MuRF1 in zebrafish , Atlantic salmon [22,77-79], pacu (Piaractus mesopotamicus)  and gilthead sea bream . Similarly in another study in rainbow trout, the mRNA expression of the Ub-ligase MAFbx and the level of poly-ubiquitinated proteins in the muscle was significantly increased with fasting and decreased after re-feeding, although major changes were not observed in the activity of the main proteasomal peptidases (trypsin and chymotrypsin-like) . Moreover in the same species, another study has demonstrated that the expression of MuRF genes was up-regulated with fasting as well as during spawning, suggesting an important role of the UbP system in fish during situations of induced muscle atrophy .
In summary, the present study has shown that the different proteolytic systems are transcriptionally regulated during ontogeny and according to the physiological status in gilthead sea bream muscle, and adds to the evidence that feeding regimes and/or diet can be used to alter proteolysis in the muscle of farmed fish. This may provide a practical means of manipulating the extent of protein breakdown during post-mortem storage of fish so as to reduce the problems of soft flesh and gaping, which reduce economic value.
APROMAR-ESACUA: La Acuicultura en España. Ministerio de Agricultura. Alimentación y Medio Ambiente: Madrid; 2013.
García de la serrana D, Vieira VL, Andree KB, Darias M, Estévez A, Gisbert E, et al. Development temperature has persistent effects on muscle growth responses in gilthead sea bream. PLoS One. 2012;7:e51884.
García de la serrana D, Codina M, Capilla E, Jiménez-Amilburu V, Navarro I, Du S-J, et al. Characterisation and expression of myogenesis regulatory factors during in vitro myoblast development and in vivo fasting in the gilthead sea bream (Sparus aurata). Comp Biochem Physiol A Mol Integr Physiol. 2014;167:90–9.
Ibarz A, Felip O, Fernández-Borràs J, Martín-Pérez M, Blasco J, Torrella JR. Sustained swimming improves muscle growth and cellularity in gilthead sea bream. J Comp Physiol B. 2011;181:209–17.
Jiménez-Amilburu V, Salmerón C, Codina M, Navarro I, Capilla E, Gutiérrez J. Insulin-like growth factors effects on the expression of myogenic regulatory factors in gilthead sea bream muscle cells. Gen Comp Endocrinol. 2013;188:151–8.
Rius-Francino M, Acerete L, Jiménez-Amilburu V, Capilla E, Navarro I, Gutiérrez J. Differential effects on proliferation of GH and IGFs in sea bream (Sparus aurata) cultured myocytes. Gen Comp Endocrinol. 2011;172:44–9.
Sánchez-Gurmaches J, Cruz-Garcia L, Ibarz A, Fernández-Borrás J, Blasco J, Gutiérrez J, et al. Insulin, IGF-I, and muscle MAPK pathway responses after sustained exercise and their contribution to growth and lipid metabolism regulation in gilthead sea bream. Domest Anim Endocrinol. 2013;45:145–53.
Valente LMP, Cornet J, Donnay-Moreno C, Gouygou JP, Bergé JP, Bacelar M, et al. Quality differences of gilthead sea bream from distinct production systems in Southern Europe: Intensive, integrated, semi-intensive or extensive systems. Food Control. 2011;22:708–17.
Vélez EJ, Lutfi E, Jiménez-Amilburu V, Riera-Codina M, Capilla E, Navarro I, et al. IGF-I and amino acids effects through TOR signaling on proliferation and differentiation of gilthead sea bream cultured myocytes. Gen Comp Endocrinol. 2014;205:296–304.
Johnston IA. Environment and plasticity of myogenesis in teleost fish. J Exp Biol. 2006;209:2249–64.
Rowlerson A, Mascarello A, Radaelli G, Veggetti A. Differentiation and growth of muscle in the fish Sparus aurata (L.) II. Hyperplastic and hypertrofic growth of lateral muscle from hatching to adult. J Muscle Res Cell Motil. 1995;16:223–36.
Mommsen TP. Paradigms of growth in fish. Comp Biochem Physiol B Biochem Mol Biol. 2001;129:207–19.
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.
Jackman RW, Kandarian SC. The molecular basis of skeletal muscle atrophy. Am J Physiol Cell Physiol. 2004;287:C834–43.
Kachaeva EV, Shenkman BS. Various jobs of proteolytic enzymes in skeletal muscle during unloading: facts and speculations. J Biomed Biotechnol. 2012;2012:493618.
Salem M, Kenney PB, Rexroad CE, Yao J. Molecular characterization of muscle atrophy and proteolysis associated with spawning in rainbow trout. Comp Biochem Physiol Part D Genomics Proteomics. 2006;1:227–37.
Goll DE, Thompson VF, Li H, Wei W, Cong J. The calpain system. Physiol Rev. 2003;83:731–801.
Attaix D, Baracos VE. MAFbx/Atrogin-1 expression is a poor index of muscle proteolysis. Curr Opin Clin Nutr Metab Care. 2010;13:223–4.
Ono Y, Sorimachi H. Calpains: an elaborate proteolytic system. Biochim Biophys Acta. 1824;2012:224–36.
Glass DJ. Signaling pathways perturbing muscle mass. Curr Opin Clin Nutr Metab Care. 2010;13:225–9.
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 NI, Taylor RG, Johnston IA. Phasing of muscle gene expression with fasting-induced recovery growth in Atlantic salmon. Front Zool. 2009;6:18.
Lecker SH, Goldberg AL. Slowing muscle atrophy: putting the brakes on protein breakdown. J Physiol. 2002;545:729.
Lecker SH, Jagoe RT, Gilbert A, Gomes M, Baracos V, Bailey J, et al. Multiple types of skeletal muscle atrophy involve a common program of changes in gene expression. FASEB J. 2004;18:39–51.
Medina R, Wing SS, Goldberg AL. Increase in levels of polyubiquitin and proteasome mRNA in skeletal muscle during starvation and denervation atrophy. Biochem J. 1995;307:631–7.
Altun M, Grönholdt-Klein M, Wang L, Ulfhake B. Cellular degradation machineries in age-related loss of muscle mass (Sarcopenia). In: Nagata T, editor. Senescence. Rijeka, Croatia: InTech; 2012.
Cai D, Lee KKH, Li M, Tang MK, Chan KM. Ubiquitin expression is up-regulated in human and rat skeletal muscles during aging. Arch Biochem Biophys. 2004;425:42–50.
Beaulieu MA, Guderley H. Changes in qualitative composition of white muscle with nutritional status of Atlantic cod, Gadus morhua. Comp Biochem Physiol A Mol Integr Physiol. 1998;121:135–41.
Guderley H, Lapointe D, Bédard M, Dutil J-D. Metabolic priorities during starvation: enzyme sparing in liver and white muscle of Atlantic cod, Gadus morhua L. Comp Biochem Physiol A Mol Integr Physiol. 2003;135:347–56.
Johnston IA. Quantitative analysis of muscle breakdown during starvation in the marine flatfish (Pleuronectes platessa). Cell Tissue Res. 1981;214:369–86.
Froehlich JM, Fowler ZG, Galt NJ, Smith DL, Biga PR. Sarcopenia and piscines: the case for indeterminate-growing fish as unique genetic model organisms in aging and longevity research. Front Genet. 2013;4:159.
Salem M, Silverstein J, Rexroad CE, Yao J. Effect of starvation on global gene expression and proteolysis in rainbow trout (Oncorhynchus mykiss). BMC Genomics. 2007;8:328.
Rescan P-Y, Montfort J, Rallière C, Le Cam A, Esquerré D, Hugot K. Dynamic gene expression in fish muscle during recovery growth induced by a fasting-refeeding schedule. BMC Genomics. 2007;8:438.
Seiliez I, Panserat S, Skiba-Cassy S, Fricot A, Vachot C, Kaushik S, et al. Feeding status regulates the polyubiquitination step of the ubiquitin-proteasome-dependent proteolysis in rainbow trout (Oncorhynchus mykiss) muscle. J Nutr. 2008;138:487–91.
Salem M, Kenney PB, Rexroad CE, Yao J. Microarray gene expression analysis in atrophying rainbow trout muscle: a unique non-mammalian muscle degradation model. Physiol Genomics. 2006;28:33–45.
Turk V, Stoka V, Vasiljeva O, Renko M, Sun T, Turk B, et al. Cysteine cathepsins: from structure, function and regulation to new frontiers. Biochim Biophys Acta. 1824;2012:68–88.
Delbarre-Ladrat C, Chéret R, Taylor R, Verrez-Bagnis V. Trends in postmortem aging in fish: understanding of proteolysis and disorganization of the myofibrillar structure. Crit Rev Food Sci Nutr. 2006;46:409–21.
Mommsen TP. Salmon spawning migration and muscle protein metabolism: the August Krogh principle at work. Comp Biochem Physiol B Biochem Mol Biol. 2004;139:383–400.
Carnevali O, Carletta R, Cambi A, Vita A, Bromage N. Yolk formation and degradation during oocyte maturation in seabream Sparus aurata: involvement of two lysosomal proteinases. Biol Reprod. 1999;60:140–6.
Carnevali O, Centonze F, Brooks S, Marota I, Sumpter JP. Molecular cloning and expression of ovarian cathepsin D in seabream, Sparus aurata. Biol Reprod. 1999;61:785–91.
Carnevali O, Cionna C, Tosti L, Cerdà J, Gioacchini G. Changes in cathepsin gene expression and relative enzymatic activity during gilthead sea bream oogenesis. Mol Reprod Dev. 2008;75:97–104.
Bahuaud D, Gaarder M, Veiseth-Kent E, Thomassen M. Fillet texture and protease activities in different families of farmed Atlantic salmon (Salmo salar L.). Aquaculture. 2010;310:213–20.
Caballero MJ, Betancor M, Escrig JC, Montero D, Espinosa De Los Monteros A, Castro P, et al. Post mortem changes produced in the muscle of sea bream (Sparus aurata) during ice storage. Aquaculture. 2009;291:210–6.
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.
Chomczynski P, Sacchi N. The single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction: twenty-something years on. Nat Protoc. 2006;1:581–5.
Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Marchler-Bauer A, Zheng C, Chitsaz F, Derbyshire MK, Geer LY, Geer RC, et al. CDD: conserved domains and protein three-dimensional structure. Nucleic Acids Res. 2013;41:D348–52.
Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, et al. SMART 4.0: towards genomic data integration. Nucleic Acids Res. 2004;32:D142–4.
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.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 − ΔΔCT method. Methods. 2001;25:402–8.
Zaidi N, Maurer A, Nieke S, Kalbacher H. Cathepsin D: a cellular roadmap. Biochem Biophys Res Commun. 2008;376:5–9.
Nothwang HG, Tamura T, Tanaka K, Ichihara A. Sequence analyses and inter-species comparisons of three novel human proteasomal subunits, HsN3, HsC7-I and HsC10-II, confine potential proteolytic active-site residues. Biochim Biophys Acta. 1994;1219:361–8.
Martin SAM, Blaney S, Bowman AS, Houlihan DF. Ubiquitin-proteasome-dependent proteolysis in rainbow trout (Oncorhynchus mykiss): effect of food deprivation. Pflugers Arch. 2002;445:257–66.
Heinemeyer W, Ramos PC, Dohmen RJ. The ultimate nanoscale mincer: assembly, structure and active sites of the 20S proteasome core. Cell Mol Life Sci. 2004;61:1562–78.
Unno M, Mizushima T, Morimoto Y, Tomisugi Y, Tanaka K, Yasuoka N, et al. The structure of the mammalian 20S proteasome at 2.75 Å resolution. Structure. 2002;10:609–18.
Hayashi T, Noga M, Matsuda M. Nucleotide sequence and expression of the rat polyubiquitin mRNA. Biochim Biophys Acta. 1994;1218:232–4.
Berti PJ, Storer AC. Alignment/phylogeny of the papain superfamily of cysteine proteases. J Mol Biol. 1995;246:273–83.
Taylor JS, Braasch I, Frickey T, Meyer A, Van de Peer Y. Genome duplication, a trait shared by 22,000 species of ray-finned fish. Genome Res. 2003;13:382–90.
Fortenberry SC, Schorey JS, Chirgwin JM. Role of glycosylation in the expression of human procathepsin D. J Cell Sci. 1995;108:2001–6.
Bourin M, Gautron J, Berges M, Nys Y, Réhault-Godbert S. Sex- and tissue-specific expression of “similar to nothepsin” and cathepsin D in relation to egg yolk formation in Gallus gallus. Poult Sci. 2012;91:2288–93.
Brooks S, Tyler CR, Carnevali O, Coward K, Sumpter JP. Molecular characterisation of ovarian cathepsin D in the rainbow trout, Oncorhynchus mykiss. Gene. 1997;201:45–54.
Kurokawa T, Uji S, Suzuki T. Identification of pepsinogen gene in the genome of stomachless fish, Takifugu rubripes. Comp Biochem Physiol B Biochem Mol Biol. 2005;140:133–40.
Mei Y, Chen Y, Li J, Gao P, Wang C, Zhang H, et al. Sequence identification, tissue distribution and polymorphism of the porcine cathepsin D (CTSD) gene. Anim Biotechnol. 2008;19:144–58.
Riggio M, Scudiero R, Filosa S, Parisi E. Sex- and tissue-specific expression of aspartic proteinases in Danio rerio (zebrafish). Gene. 2000;260:67–75.
Sakai H, Saku T, Kato Y, Yamamoto K. Quantitation and immunohistochemical localization of cathepsins E and D in rat tissues and blood cells. Biochim Biophys Acta. 1989;991:367–75.
Hiramatsu N, Ichikawa N, Fukada H, Fujita T, Sullivan CV, Hara A. Identification and characterization of proteases involved in specific proteolysis of vitellogenin and yolk proteins in salmonids. J Exp Zool. 2002;292:11–25.
Kwon JY, Prat F, Randall C, Tyler CR. Molecular characterization of putative yolk processing enzymes and their expression during oogenesis and embryogenesis in rainbow trout (Oncorhynchus mykiss). Biol Reprod. 2001;65:1701–9.
Park E-M, Kim YO, Nam B-H, Kong HJ, Kim W-Y, Lee S-J, et al. Cloning and expression analysis of cathepsin D in the olive flounder Paralichthys olivaceus. Biosci Biotechnol Biochem. 2009;73:1856–9.
Sire M-F, Babin PJ, Vernier J-M. Involvement of the lysosomal system in yolk protein deposit and degradation during vitellogenesis and embryonic development in trout. J Exp Zool. 1994;269:69–83.
Yamamura J, Adachi T, Aoki N, Nakajima H, Nakamura R, Matsuda T. Precursor-product relationship between chicken vitellogenin and the yolk proteins: the 40 kDa yolk plasma glycoprotein is derived from the C-terminal cysteine-rich domain of vitellogenin II. Biochim Biophys Acta. 1995;1244:384–94.
Lupatsch I. Protein and energy requirements in Mediterranean species. Cah Options Méditerranéennes. 2005;63:9–18.
Peragón J, Barroso JB, de la Higuera M, Lupiáñez JA. Relationship between growth and protein turnover rates and nucleic acids in the liver of rainbow trout (Oncorhynchus mykiss) during development. Can J Fish Aquat Sci. 1998;55:649–57.
Peragón J, Barroso JB, García-Salguero L, de la Higuera M, Lupiáñez JA. Growth, protein-turnover rates and nucleic-acid concentrations in the white muscle of rainbow trout during development. Int J Biochem Cell Biol. 2001;33:1227–38.
Edström E, Altun M, Hägglund M, Ulfhake B. Atrogin-1/MAFbx and MuRF1 are downregulated in aging-related loss of skeletal muscle. J Gerontol A Biol Sci Med Sci. 2006;61:663–74.
Foletta VC, White LJ, Larsen AE, Léger B, Russell AP. The role and regulation of MAFbx/atrogin-1 and MuRF1 in skeletal muscle atrophy. Pflugers Arch. 2011;461:325–35.
Komander D, Rape M. The ubiquitin code. Annu Rev Biochem. 2012;81:203–29.
Bower NI, Johnston IA. Discovery and characterization of nutritionally regulated genes associated with muscle growth in Atlantic salmon. Physiol Genomics. 2010;42A:114–30.
Bower NI, García de la serrana D, Johnston IA. Characterisation and differential regulation of MAFbx/Atrogin-1 α and β transcripts in skeletal muscle of Atlantic salmon (Salmo salar). Biochem Biophys Res Commun. 2010;396:265–71.
Valente LMP, Bower NI, Johnston IA. Postprandial expression of growth-related genes in Atlantic salmon (Salmo salar L.) juveniles fasted for 1 week and fed a single meal to satiation. British J Nutr. 2012;108:2148–57.
Mareco EA. García de la serrana D, Johnston IA, Dal-Pai-Silva M. Characterization of the transcriptome of fast and slow muscle myotomal fibres in the pacu (Piaractus mesopotamicus). BMC Genomics. 2015;16:182.
Wang J, Salem M, Qi N, Kenney PB, Rexroad CE, Yao J. Molecular characterization of the MuRF genes in rainbow trout: Potential role in muscle degradation. Comp Biochem Physiol B Biochem Mol Biol. 2011;158:208–15.
Kaitetzidou E, Crespo D, Vraskou Y, Antonopoulou E, Planas JV. Transcriptomic response of skeletal muscle to lipopolysaccharide in the gilthead seabream (Sparus aurata). Mar Biotechnol NY. 2012;14:605–19.
Pinto PI, Matsumura H, Thorne MA, Power DM, Terauchi R, Reinhardt R, et al. Gill transcriptome response to changes in environmental calcium in the green spotted puffer fish. BMC Genomics. 2010;11:476.
The authors would like to thank Dr. Tom Ashton for his assistance with the cloning and Natàlia Riera and José Rives for their help with the tissue screening. We thank Carlos Mazorra from Tinamenor S.L. (Pesués, Spain) and the “Institut de Recerca i Tecnologia Agroalimentàries” (IRTA, Sant Carles de la Ràpita, Spain) for the gilthead sea bream used in the study and the personnel from the animal facilities at the Faculty of Biology for its maintenance. EC is a Ramón y Cajal researcher fellow from the “Ministerio de Ciencia e Innovación” (MICINN). This work was supported by funds from the MICINN (AGL2009-12427 and AGL2010-17324), the Catalonian Government (2009SGR-00402) and the “Xarxa de Referència d’R + D + I en Aqüicultura” and, the European Union through the project LIFECYCLE (FP7-222719).
The authors declare that they have no competing interests.
IAJ, JG and EC conceived the study. CS and EC performed the in vivo experiments and all the laboratory analyses. All authors analysed the data, wrote, read and approved the final manuscript.
Nucleotide and deduced amino acid sequences of gilthead sea bream cathepsin B (SaCTSB, [GenBank: KJ524457]). Figure S2. Nucleotide and deduced amino acid sequences of a new paralogue of gilthead sea bream cathepsin D (SaCTSDb, [GenBank: KJ524456]). Figure S3. Partial nucleotide and deduced amino acid sequences of gilthead sea bream proteasome subunit beta type-4 (SaN3, also known as PSMB4, [GenBank: KJ524458]). Figure S4. Nucleotide and deduced amino acid sequences of gilthead sea bream ubiquitin (SaUb, [GenBank: KJ524459]).
About this article
- White muscle
- Life-history stages
- Lysosomal proteolysis
- Ubiquitin-proteasome pathway