Skip to main content
Download PDF

Differences between fast and slow muscles in scallops revealed through proteomics and transcriptomics

BMC Genomics volume 19, Article number: 377 (2018) Cite this article

Abstract

Background

Scallops possess striated and catch adductor muscles, which have different structure and contractile properties. The striated muscle contracts very quickly for swimming, whereas the smooth catch muscle can keep the shells closed for long periods with little expenditure of energy. In this study, we performed proteomic and transcriptomic analyses of differences between the striated (fast) and catch (slow) adductor muscles in Yesso scallop Patinopecten yessoensis.

Results

Transcriptomic analysis reveals 1316 upregulated and 8239 downregulated genes in slow compared to fast adductor muscle. For the same comparison, iTRAQ-based proteomics reveals 474 differentially expressed proteins (DEPs), 198 up- and 276 downregulated. These DEPs mainly comprise muscle-specific proteins of the sarcoplasmic reticulum, extracellular matrix, and metabolic pathways. A group of conventional muscle proteins—myosin heavy chain, myosin regulatory light chain, myosin essential light chain, and troponin—are enriched in fast muscle. In contrast, paramyosin, twitchin, and catchin are preferentially expressed in slow muscle. The association analysis of proteomic and transcriptomic data provides the evidences of regulatory events at the transcriptional and posttranscriptional levels in fast and slow muscles. Among 1236 differentially expressed unigenes, 22.7% show a similar regulation of mRNA levels and protein abundances. In contrast, more unigenes (53.2%) exhibit striking differences between gene expression and protein abundances in the two muscles, which indicates the existence of fiber-type specific, posttranscriptional regulatory events in most of myofibrillar proteins, such as myosin heavy chain, titin, troponin, and twitchin.

Conclusions

This first, global view of protein and mRNA expression levels in scallop fast and slow muscles reveal that regulatory mechanisms at the transcriptional and posttranscriptional levels are essential in the maintenance of muscle structure and function. The existence of fiber-type specific, posttranscriptional regulatory mechanisms in myofibrillar proteins will greatly improve our understanding of the molecular basis of muscle contraction and its regulation in non-model invertebrates.

Background

Adult skeletal muscles are composed of two main classes of fiber types, slow, type I, and fast, type II, which are classified according to myofibrillar ATP staining and immunohistochemistry [1, 2]. Fast and slow muscles, composed of fibers with distinct physiological properties, play indispensable roles in body motion and maintenance of metabolism [3, 4]. Slow and fast muscles occur not only in vertebrates but also in many invertebrates, such as molluscs and crustaceans [5,6,7,8]. Over the past decades, extensive work has been done on invertebrate muscles, primarily in general structure of muscles, regulation of muscle contraction, and the mechanisms of contraction and motor function [9,10,11,12,13]. More recently, our knowledge of the unique physical properties of invertebrate muscles, has been reformulated with the advent of structural knowledge of thick and thin filaments [14,15,16,17,18,19], in vitro motility assay [20, 21], and new evidences on the catch mechanism [22,23,24,25]. Recent work on fast and slow skeletal muscles reveals the complex regulatory mechanisms in key components of muscle structure, including the transcriptional and posttranscriptional events [3, 4]. In contrast, invertebrate thick filaments are very different from vertebrate striated thick filaments and show great variation within invertebrates [26, 27]. Muscle-specific genes and proteins in invertebrates (e.g. paramyosin, twitchin, and catchin) have been identified in the past decades, but the regulatory mechanisms for differential expression of these muscle-specific proteins remain largely unknown [26].

Scallops possess fast (striated) and slow (catch) adductor muscles, which lie closely apposed to one another but are divided by a connective tissue sheet. The striated adductor muscle contracts very quickly for swimming, whereas smooth catch adductor muscle lacks striations, contracts for long periods, keeping shells closed with little expenditure of energy [9, 10]. For this reason, scallops are a key model for studies on muscle structure and function.

For scallops, myosin and actin filaments are most abundant in adductor muscles, where they form the sarcomere in fast adductor muscle and the less ordered contractile units in catch muscle [9]. For fast adductor muscle, calcium binds not only to regulatory components of the thin filament but also directly to specific regulatory sites in order to activate the interaction of myosin with actin during phasic contraction (See the review from Chantler [10]). In contrast, twitchin phosphorylation mechanisms in scallop catch muscle may occur in place of direct activation by calcium binding during catch contraction [28]. Although we know the structural and functional distinctions between the two muscles, little information is available on their molecular components and regulatory mechanisms.

Our current knowledge on muscle specific proteins in scallops has been established primarily by traditional biochemical and immunohistochemical techniques, which require laborious work and yield a limited number of enzymes or proteins [10]. Recently, the advances in iTRAQ (isobaric tags for relative and absolute quantitation) proteomics and RNA-Seq transcriptomics enable us to identify and quantify thousands of genes and proteins in a single experiment [3, 29,30,31]. Here, we performed the first, integrated proteomic and transcriptomic analyses of differences between the fast (striated) and slow (catch) adductor muscles in Yesso scallop (Patinopecten yessoensis; Fig. 1), with the goal to uncover the molecular components, as well as their regulatory mechanisms in fast and slow adductor muscles of scallops. The comparison of protein and mRNA data in this study provides insights into the regulatory events at transcriptional and posttranscriptional levels that shape the identity of fast and slow adductor muscles. The present results will greatly extend our knowledge on the molecular basis of muscle contraction and its regulation in scallops.

Fig. 1
figure1

The striated (fast) and catch (slow) adductor muscles of Yesso scallop Patinopecten yessoensis. a, Photograph of adductor muscles in a live scallop; b, Photograph of adductor muscles after removal of all organs and tissues

Full size image

Results

Overall statistics for transcriptome and proteome sequencing

Illumina sequencing yielded 93,868,386 raw reads (NCBI accession number: SRR4254476). Filtering low quality reads resulted in clean reads totaling 10.73 Gb, with a GC content of 45.80%. For the four sequencing samples, similar Q20 were observed ranging from 96.83 to 97.16%. After assembly, the transcripts were subsequently assembled into 56,422 unigenes, with the mean length of 584 bp and N50 of 704 bp. The maximum number of unigenes were annotated in Nr (NCBI non-redundant protein) database (18,571), accounting for 32.91% of all unigenes, followed by 25.15% in Swissprot and 21.93% in KOG (euKaryotic Ortholog Groups) database, and 14.74% in KEGG (Kyoto Encyclopedia of Genes and Genomes) database.

Mass spectrometry, iTRAQ proteomics yielded 43,280 unique spectra among the 319,156 spectra detected (ProteomeXchange Consortium dataset identifier PXD005166). Finally, 1591 proteins were identified and annotated in the protein databases. According to Gene Ontology (GO), there were 958 unigenes assigned to three main functional categories, including biological process (BP), cellular component (CC), and molecular function (MF). A total of 3277 unigenes were classified into 26 ortholog groups in KOG database. The pathways determined by KEGG database were grouped into six specific pathways, including metabolism, human diseases, genetic information processing, cellular processes, environmental information processing, and organismal systems.

Identification of DEGs between the striated and catch muscles

Counts of expressed tags were used to estimate levels of gene expression for replicate samples of catch and striated muscles. Expression levels are tightly correlated between replicates within tissue (Pearson correlation coefficients, 0.9985 and 0.9907, for catch and striated muscle tissues, respectively). Comparison of expression levels between catch and striated muscles produced 9555 differentially expressed genes (DEGs), with the filter criteria of false discovery rate (FDR) < 0.05 and |log2FoldChange| > 1. Among the DEGs, there were 1316 upregulated and 8239 downregulated genes identified between catch and striated muscles. Volcano plot illustrates the asymmetry in differentially expressed genes (Fig. 2) between downregulated (green) and upregulated (red) genes.

Fig. 2
figure2

Volcano plots for differentially expressed genes (DEGs) between slow and fast adductor muscles. Points of the plots represent transcripts that are significantly differentially expressed. Green points represent transcripts with significantly lower expression level (slow versus fast muscle), while red circles represent transcripts with significantly higher expression level (P < 0.05)

Full size image

DEGs were assigned to 41 sub-categories of GO terms (Fig. 3). Most of the unigenes were assigned to cellular process, localization, metabolic process, single-organism process, cell, cell part, binding, and catalytic activity. The DEGs were further subjected to KEGG pathway analysis, which were classified into 222 signaling pathways. The most enriched pathways include calcium signaling pathway, protein digestion and absorption, vascular smooth muscle contraction, ECM (extracellular matrix)-receptor interaction, cardiac muscle contraction, glycolysis/gluconeogenesis, regulation of actin cytoskeleton, MAPK signaling pathway and hedgehog signaling pathway.

Fig. 3
figure3

Gene Ontology (GO) enrichment results for DEGs between slow and fast muscle. HS, catch (slow) muscle; PL, striated (fast) muscle. Red bars represent upregulated genes for slow versus fast muscles, while green bars represent downregulated genes in the same comparison

Full size image

Eight DEGs related to muscle structure and function, including twitchin-like isoform X7, titin, calponin, calcium-ATPase, paramyosin, catchin, myosin, and myosin essential light chain, were selected to verify the RNA-Seq results by the quantitative RT-PCR (qPCR). Compared with the striated muscle, the catch muscle had a significantly higher expression of twitchin-like isoform X7, titin, paramyosin, and catchin, whereas the other four genes had a significantly lower expression (p < 0.01; Fig. 4a). According to the transcriptomic data, the RPKM (Reads per kb per million reads) values of the selected genes were summarized in Fig. 4b, which showed that the expression pattern of these genes was all in accordance with the qPCR results.

Fig. 4
figure4

Quantitative RT-PCR (qPCR) results of muscle-related genes between slow and fast muscles of P. yessoensis. a, the relative expression obtained from qPCR analysis; b, the RPKM (reads per kb per million reads) values of the selected genes according to transcriptomic data. HS, catch (slow) muscle; PL, striated (fast) muscle

Full size image

Identification of DEPs between the striated and catch muscles

The iTRAQ-based quantitative proteomics identified 474 differentially expressed proteins (DEPs) between the slow and fast muscles, including 198 upregulated and 276 downregulated DEPs (Additional file 1: Table S1). These DEPs include dozens of muscle-specific proteins, such as calmodulin, myosin essential light chain, myosin heavy chain (MHC) II, titin, tropomodulin, troponin C, myosin essential light chain, paramyosin, gelsolin-like protein, myopalladin, galponin, myophilin, and twitchin. Furthermore, quantitative proteomics revealed many proteins related to calcium signaling, sarcomere and cytoskeleton, including sarcoplasmic calcium-binding protein, smoothelin-like protein, actinin, filamin, actin-interacting protein, PDZ and LIM domain protein.

Similarly, the most enriched KEGG pathways include ECM-receptor interaction, glycolysis/gluconeogenesis, cardiac muscle contraction, protein digestion and absorption, regulation of actin cytoskeleton, calcium signaling pathway, vascular smooth muscle contraction, hedgehog signaling pathway, etc. In the pathway of ECM-receptor interaction, ECM related proteins, such as collagen, laminin, and perlecan, were expressed significantly higher in catch muscle as compared with striated muscle.

Association analysis of transcriptome and proteome

Association analysis of transcriptome and proteome data for slow versus fast adductor muscles revealed a nonlinear relationship between mRNA and protein expression, with the Pearson’s correlation coefficient of 0.5249 (Fig. 5). The colored dots representing 1236 differentially expressed unigenes (slow vs. fast muscle) were unevenly distributed in nine groups. Red dots in groups 3 and 7 represent the elevated expression of mRNA and protein levels in fast (group 3) and slow (group 7) muscles; Red dots in groups 1 and 9 show the opposite changes of mRNA levels and protein abundances; Green dots (group 2 and 8) indicate changes of mRNA expression levels only; Blue dots (group 4 and 6) denote changes of expression in protein levels only; Grey dots (group 5) display no significant change of expression in either mRNA or proteins.

Fig. 5
figure5

Association analysis of transcriptome and proteome differences between slow and fast muscles in P. yessoensis (protein fold changes > 1.2 and mRNA fold change > 2). Red dots (group 1, 3, 7, and 9) represent significant changes of expression in both mRNA and protein; Green dots (group 2 and 8) show significant changes in mRNA expression levels only; Blue dots (group 4 and 6) denote significant changes of expression in protein levels only; Grey dots (group 5) display no significant change of expression in either mRNA or proteins

Full size image

In groups 3 and 7, the similar changes in mRNA and protein levels suggest that expression of these proteins in the two muscles is regulated at the transcriptional level. 281 unigenes (22.7%) showed a positive relationship between protein abundances and mRNA enrichment in the striated and catch muscles, including myofibrillar proteins, sarcoplasmic reticulum (SR) proteins, metabolism related enzymes, and membrane and extracellular proteins (Additional file 2: Table S2). In group 3, the identified unigenes mainly represent proteins enriched in striated muscle, such as myosin, actinin, calmodulin, enolase, titin, tropomodulin, and troponin. These muscle proteins were predicted to locate in the sarcomere, as seen in vertebrate striated muscles (Fig. 6). Significantly higher expression of the main components of myosin, including myosin heavy chain (MHC) and myosin essential light chain (E-LC), was detected in striated muscle compared to catch muscle. The conventional (class II) and unconventional myosin classes were found to be expressed in a tissue-specific manner. The MHC class II consistently expressed at a much higher level in striated muscle than in catch muscle, whereas unconventional myosin classes (e.g. unconventional myosin-XVI-like isoform X3 and non-muscle myosin) showed significantly higher expression in catch muscle than in striated muscle. Three troponin subunits, troponin C, troponin I and troponin T, also showed higher expression, at mRNA and protein levels, in striated muscle compared to catch muscle. Moreover, three important sarcoplasmic reticulum (SR) proteins, sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA), calcium-transporting ATPase, and sarcoplasmic calcium-binding protein (SCP), showed significantly higher mRNA and protein expression in striated muscle than in catch muscle. In addition, many proteins involved in calcium signaling pathway and metabolic pathways were also detected in group 3, such as calcium binding protein and arginine kinase. Moreover, unigenes encoding key enzymes involved in glycolysis were detected in group 3, such as glycogen phosphorylase, glycogen debranching enzyme and pyruvate kinase. The levels of all of these enzymes were positively correlated with mRNA levels in striated muscle. For arginine kinase, mRNA and protein expression showed approximately 4.5-fold and 3-fold higher expression in striated muscle than in catch muscle. Identified unigenes in group 7, in contrast to those in group 3, mainly represent proteins enriched in catch muscle, including paramyosin, twitchin, catchin, calponin, gelsolin-like protein, myophilin, smoothelin-like protein, and titin. Overall, expression patterns of mRNA and proteins were similar in most of the myofibrillar proteins, except for actin, tropomyosin, AMP deaminase, filamin, crystallins, and CapZ (Fig. 6).

Fig. 6
figure6

Schematic representation of fiber structure and composition for striated (fast) adductor muscle in P. yessoensis, based on its vertebrate counterparts. Localization of proteins identified was predicted from transcriptomic and proteomic analyses. Red fonts represent those proteins with significantly higher expression in fast muscle; Green fonts represent proteins with significantly lower expression in fast muscle; Black fonts represent no quantitative difference in protein levels between slow and fast muscles. Numbers in parentheses are the group distributions, according to the association analysis of transcriptome and proteome data

Full size image

In contrast, 658 unigenes (53.2%), in groups 1, 2, 4, 6, 8 and 9, display the striking differences between protein and mRNA expression levels, which suggest the existence of regulatory events at the posttranscriptional level. For instance, 89 unigenes in groups 1 and 9 had a negative correlation of mRNA expression and protein abundances, such as filamin, 26S proteasome, ATP-dependent RNA helicase, calpain, E3 ubiquitin-protein ligase, collagen, laminin, etc. Additionally, 473 genes (blue dots) in groups 4 and 6 showed elevated protein levels in either the catch or striated muscles, but with no difference in mRNA expression. They were mainly associated with collagen, ribosomal proteins, glycogen synthase and mitochondrial genes. Finally, 96 unigenes in groups 2 and 8 (green dots), such as insulin-like growth factor, glycogen synthase, AMP deaminase, 26S protease regulatory subunit, and lamin, showed elevated mRNA levels in either the catch or striated muscles but had similar protein expression in the two muscles.

Discussion

In this study, we performed quantitative transcriptomic and iTRAQ-based proteomic analyses of differences between the striated and catch adductor muscles of Yesso scallop P. yessoensis. Proteomic analysis quantified more than 1500 proteins, comprising the largest quantitative data set to date for proteins in scallop adductor muscles. The discovery of ~ 500 DEPs and 9555 DEGs greatly extends our knowledge on the molecular components and complex regulatory events in molluscan muscle physiology. We summarize our findings below about some of the major molecular components of the striated and catch adductor muscles, and highlight that regulatory mechanisms at the transcriptional and posttranscriptional levels are essential in the maintenance of muscle structure and function in scallops.

Myosin heavy chain

Scallop adductor myosin is a regulatory myosin, which possesses the necessary machinery within its own structure to control its interaction with actin by the cross-bridge cycle [10]. Similar to other conventional myosins, the scallop myosin is composed of two heads with myosin heavy chains (MHCs), followed by a neck domain to which the light chains (LCs) bind, and a long coiled-coil tail. A recent study revealed that the myosin heads in Ca2+-regulated myosin filaments of scallop striated muscles interact in a similar way to those in phosphorylation-regulated filaments of vertebrate smooth muscles [17]. In this study, relatively higher expression levels of E-LC and R-LC are detected in the striated muscle as compared with the catch muscle, which may correspond to high actin-dependent MgATPase activity in the striated muscle.

Scallop muscles have a similar myosin-linked regulation mechanism as vertebrates, because they share a common structural basis for switching off thick-filament activity in relaxed muscles [10, 17]. However, unlike their vertebrate counterparts, the striated and catch muscle MHCs in scallop myosin heads are expressed from a single gene as alternatively spliced products [32]. Both conventional (Class II) and unconventional myosin classes were detected in the mantle tissue of scallop P. yessoensis in previous studies [33,34,35], as well as in the present study. However, the mantle MHC II is similar to vertebrate non-muscle MHC II rather than scallop muscle MHC II, having only 22% sequence identity with scallop muscle MHC II at the amino acid level. Similarly, we also found the extremely low sequence identity (< 20%) between the unconventional MHC and conventional MHC sequences. In the present study, four MHC mRNA variants produced from a single gene, in groups 2, 3 and 6, display the differential expression between mRNA levels and protein abundances, which indicate that regulation events may occur at the posttranscriptional level, such as alternative RNA splicing.

Troponin

Troponin is the tropomyosin-binding protein, a central element in the thin-filament-linked Ca2+ regulatory system of vertebrate striated muscles [36,37,38]. It is often composed of three subunits: troponin C, troponin I and troponin T. As in vertebrate striated muscle, we detect an equal ratio of C:I:T troponin subunits in both striated and catch muscles. The high expression of three troponin subunits in the striated muscle supports the prevailing view being that scallop striated muscles possess some form of thin filament-linked regulation [10]. However, we find that sequences of scallop troponin subunits exhibit low (< 30%) homology with their vertebrate counterparts, as reported in previous studies [39, 40]. Thus, the different structures and functional adaptations of the three molluscan troponin subunits suggest a very different mechanism of action to that seen in vertebrates. As indicated in molluscs, troponin may regulate contraction through activating mechanisms involving the structural troponin C binding site as well as the inhibitory region of troponin I [10, 41]. In the present study, troponin C and T show the increased mRNA and protein expression in striated muscle compared to catch muscle, while troponin C exhibit significant changes of expression in protein levels only. The distinct protein and isoform expression of troponin in striated and catch muscles suggests the existence of fiber-type specific, posttranscriptional regulatory mechanisms.

Titin

Titin is an exceptionally large protein in vertebrate muscles, which can link filamin and α-actinin together in the Z-line periphery of the striated muscle, or in dense bodies of the smooth muscle [42, 43]. In the present study, multiple isoforms of titin are enriched in either the striated or catch muscle (seven dots in group 3 and three dots in group 7; Fig. 5), which are probably generated by alternative RNA splicing from the same gene, as suggested in vertebrate muscles [42, 44]. Therefore, the differential protein and isoform expression of titin in striated and catch muscles may be regulated by complex mechanisms, including the transcriptional and fiber-type specific, posttranscriptional regulatory events.

Other striated and catch muscle components

The similar regulations of mRNA and protein abundances are found in many muscle-specific components, which indicate that the regulatory events occur at the transcriptional level. Calponin is responsible for cross-linking actin and myosin filaments and can organize the contractile filaments into a three-dimensional network for the proper orientation and spatial distribution during force development in vertebrate smooth muscles [45]. In contrast to vertebrate, calponin in molluscan catch muscle is thought to be involved in catch regulation, serving as a competitive inhibitor of actomyosin ATPase [24, 46]. It is evidenced that the abundance of calponin in catch muscle is predominately regulated by the transcription activation of mRNA according to the present results. Moreover, smoothelin-like 1 protein (SMTNL1) may act as a physiological regulator of muscle contraction through cAMP-activated kinases in vertebrate smooth muscles [47]. In this study, SMTNL1 is enriched in catch muscle by the transcriptional regulation and may play some indirect roles in catch muscle contraction. In addition, actinin is localized at dense bodies in the smooth muscle and Z-disk in striated muscle of vertebrates, where it forms a lattice-like structure and stabilizes the muscle contractile apparatus [48,49,50]. It is suggested that the enrichment of actinin protein in striated muscle of P. yessoensis by transcriptional regulatory events may be devoted to cross-link actin filaments from adjacent sarcomeres during Z-disk assembly. Another example for the transcriptional regulation in muscle components is sarcoplasmic reticulum (SR) proteins. As indicated in vertebrates, SR proteins in striated muscle is specialized for releasing Ca2+, following sarcolemma depolarization, in order to activate muscle contraction [51]. For scallops, the striated adductor muscle has twice as much of its surface covered with sarcoplasmic reticulum as catch muscle, which is likely associated with a higher contraction rate and reduced relaxation time in striated cells [52]. The present results reveal that the elevated mRNA and protein levels of SR-related genes (e.g. SERCA, calcium-transporting ATPase, and SCP) in striated muscle are mainly controlled by the transcriptional regulatory events. The abundance of SR proteins in striated adductor muscle of scallops may be responsible for the increasing uptake of calcium and maintenance of calcium homeostasis during fast contraction and relaxation cycles [52].

Besides the transcriptional regulatory events, the striking differences between protein abundances and mRNA expression provide evidences for regulation events at the posttranscriptional level in striated and catch adductor muscles. For instance, filamin, an actin-binding protein, may be involved in reorganizing the cytoskeleton in response to signaling events and fulfill structural functions at the Z-disk in vertebrate muscles [53, 54]. In the present study, protein abundances of filamin in catch muscle (in groups 1 and 4; Fig. 5) may be involved in the slow contraction and relaxation peculiar to catch muscle, which are probably achieved by the regulatory events at the posttranscriptional level.

Implications for catch mechanism

Catch is a mechanism found in molluscan catch muscles, in which tension is maintained at a relatively low energy cost [10, 55]. In this study, a variety of proteins preferentially expressed in the catch muscle, such as paramyosin, twitchin, catchin, and other muscle specific proteins, shed potential lights on the molecular components of catch regulation in scallops. As summarized by a recent review [10], three main theories have been put forward as possible explanations for the catch mechanism. First, the unique structure of paramyosin-rich thick filaments is thought to be responsible for direct interactions between adjacent thick filaments [26, 56, 57]. Although some new evidence shows that paramyosin may not play a direct role in the maintenance of catch, an indirect role, such as the formation of a rigid network of inter-myofilament connections [22], is more difficult to exclude. A second theory holds that catch involves the formation of a long-lived actomyosin state [58].

More recently, new theories have arisen to take into account the unique roles for catchin and twitchin in the catch mechanism [59,60,61]. Catchin has been determined as an alternatively spliced product of the MHC gene with a unique non-helical N-terminal sequence, which is identical with the C-terminal 830 residues of the MHC [62]. It has been hypothesized that its unique globular N-terminus could be involved in tethering actin while the MHC-derived C-terminal coiled-coil interacts with myosin or paramyosin at the core [10]. Similarly, we reveal that the enrichment of catchin mRNA and protein levels in catch muscle compared to striated muscle in P. yessoensis indicates a potential role of catchin in catch regulation at the transcriptional level.

Twitchins are giant kinase molecules found in both the striated and catch muscles, where they are located within the A-band and at the A-I junction. In the unphosphorylated state, twitchins could interact directly with myosin, paramyosin, and catchin during the catch state [60, 63]. A previous study shows that twitchin mRNA is expressed at higher levels in striated muscle than in catch muscle [64]. However, we show that expression of twitchin proteins is significantly higher in catch muscle than in striated muscle, which supports the central role of twitchin in the catch mechanism. More interestingly, we find that three mRNA variants encoding twitchin, in groups 1, 4 and 7 (Fig. 5), are resulted from alternative splicing of the same gene in catch muscle. The striking differences of twitchin between protein and mRNA expression levels reflect the existence of regulatory events at the transcriptional and posttranscriptional levels. The present study reveals, for the first time, the fiber-type specific, posttranscriptional regulation may be involved in the maintenance of catch muscle structure and catch regulation.

Implications for energy metabolism of striated and catch muscles

The regulatory events at the transcriptional and posttranscriptional levels are not only found in the muscle-specific genes and proteins, but also in key enzymes involved in glycogen synthesis and glycolysis. A previous study indicates that most ATP equivalents are derived from arginine phosphate, in the striated muscle during the snap response, while the energy is mainly supplied from glycolysis, in the catch muscle during valve closure [65]. Arginine phosphate is the major fuel powering phasic contractions by the striated muscles of scallops, which can generate approximately 70% of the ATP used for phasic contractions [66,67,68]. As indicated, arginine kinase catalyses the reversible conversion of ADP and arginine phosphate into arginine and ATP, the activities of which reflect reliance on rapid initial bursts of phasic contraction [65, 66]. In this study, the central role of arginine phosphate in supporting phasic contraction is therefore supported by the transcriptional regulation of arginine kinase in striated muscle. Furthermore, the similar regulation of mRNA expression and protein abundances for glycolytic enzymes suggests that catabolic processes in relation to valve snap and closure responses are predominantly regulated by the transcriptional activation of the metabolic genes. However, the striking differences of glycogen synthase between protein abundances and mRNA expression levels (in groups 2 and 6; Fig. 5) indicate the existence of the posttranscriptional regulatory event in glycogen synthesis of striated and catch muscles. Overall, it is therefore suggested that energy metabolism in scallop adductor muscles is controlled by complex mechanisms, including the transcriptional and posttranscriptional regulatory events.

Conclusions

The integrated proteomics and transcriptomics study reveals a number of muscle-specific genes and proteins in striated and catch adductor muscles of Yesso scallop P. yessoensis, including muscle contractible proteins, membrane and extracellular matrix, and enzymes of metabolic pathways. The comparison of protein and mRNA data in this study provides the insights into complex regulatory events at the transcriptional and posttranscriptional levels that shape the identity of fast and slow adductor muscles. The striking differences between protein abundances and mRNA expression levels support the existence of fiber-type specific, posttranscriptional regulatory mechanisms in striated and catch adductor muscles. These findings will greatly improve our understanding of the molecular basis of muscle contraction and its regulation in non-model invertebrates.

Methods

Biological material

Two-year-old live individuals of P. yessoensis (n = 6) were obtained from a commercial hatchery in Yantai, China. To obtain high quality of gene expression data, all of these scallops were held under the same conditions in seawater at 16 ± 2 °C for three days at YSFRI (Yellow Sea Fisheries Research Institute). Scallops were fed with Isochrysis galbana and two-thirds of the culture water was exchanged every day. The striated and catch adductor muscles were separately dissected from each scallop and frozen and stored in liquid nitrogen individually. Two biological replicates for transcriptomic and proteomic analyses were prepared by dividing the same muscle sample into two parts for total RNA and proteins extraction, respectively.

RNA-Seq library construction, sequencing and quality control

Total RNA was isolated from the collected adductor muscles with Trizol Reagent (Invitrogen), following manufacturer’s instruction. RNA purity and quality were checked using a NanoPhotometer™ spectrophotometer (Implen, CA, USA) and 1% agarose electrophoresis. Extracted mRNA was enriched by Oligo(dT) beads, with removal of rRNA using a Ribo-ZeroTM Magnetic Kit (Epicentre). Then the enriched mRNA was fragmented into short fragments and reverse transcribed into cDNA with random primers. After synthesis of second-strand cDNA, fragments were then purified, end repaired, polyadenylated, and ligated to Illumina sequencing adapters. Ligation products were size-selected by agarose gel electrophoresis and sequenced by Gene Denovo Biotechnology Co. (Guangzhou, China), using Illumina HiSeqTM 2500. Low quality reads were removed, according to the previous study [31].

De novo assembly and annotations

De novo assembly the transcriptome was carried out with the short-read assembling program Trinity [69]. The unigene expression was calculated and normalized to RPKM (Reads Per kb per Million reads). To annotate the unigenes, we used the BLASTx program [70] with an E-value threshold of e− 5 to NCBI non-redundant protein (Nr) database (http://www.ncbi.nlm.nih.gov) [71], the Swiss-Prot protein database (http://www.expasy.ch/sprot) [72], the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.genome.jp/kegg) [73], and the COG/KOG database (http://www.ncbi.nlm.nih.gov/COG) [74]. GO annotation of unigenes was analyzed by Blast2GO software.

Analysis of differentially expressed genes (DEGs) and qPCR validation

Clean reads of expressed data were mapped to the reference transcriptome using the short-read alignment tool, Bowtie2, with default parameters [75]. We used the edgeR package [76] to identify differentially expressed genes between the two muscle types, as those with a significant fold-change of ≥2 in expression level and a false discovery rate (FDR) < 0.05. All DEGs were mapped to GO terms in the Gene Ontology database (http://www.geneontology.org/) [77]. The enriched GO terms were selected with FDR (False discovery rate) < 0.05. KEGG pathway enrichment analysis was also performed to identify significantly enriched metabolic pathways or signal transduction pathways in DEGs at the same condition.

To verify data accuracy of high-throughput sequencing, quantitative real-time PCR (qPCR) was used to compare the relative mRNA expression of significantly expressed genes between the striated and catch adductor muscles. Total RNAs were extracted from the striated and catch muscles using Trizol Reagent (TAKARA). cDNA was synthesized, using the Bestar qPCR RT Kit, and Real time PCR was performed using DBI SybrGreen qPCRmasterMix according manufacture’s instructions on a Mx3000P (Agilent Stratagene). In this study, β-actin was selected as the reference gene because its RPKM values were stable among samples according to the transcriptomic data. The qPCR analysis was performed with three biological and three technical replicates. The comparative Ct method (2-△△Ct method) was used to calculate the relative gene expression of the samples, which was normalized to β-actin mRNA level. The expression data were subsequently subjected to independent t-test in SPSS 17.0 to determine whether there was any difference significant at the P < 0.05 level.

Protein extraction, digestion, and iTRAQ analysis

Total protein was extracted from adductor muscles by the cold acetone method, following previous studies [78, 79]. Briefly, muscle samples were subjected to a mixture of EDTA (Ethylene Diamine Tetraacetic Acid), 1 mM PMSF (Phenylmethanesulfonyl fluoride), and 10 mM DTT (DL-Dithiothreitol) and were the ground to disrupt the cells. After centrifugation at 25,000×g for 20 min at 4 °C, 10 mM DTT was added and incubated at − 20 °C for overnight. The remnant pellet was treated with 10 mM DDT in 1.5 ml cold acetone, followed by centrifugation and air drying. The obtained pellet was then suspended and centrifuged for 20 min at 25,000×g. Next, 10 mM DDT and 55 mM IAM was added to the sample, which was incubated in a dark room for 45 min, followed by another step of suspension. The protein was measured using the 2-D Quant Kit (General Electric Company, USA) and confirmed by SDS-PAGE method according to Qin et al. [78].

Proteins were labeled, using the iTRAQ labeling kit (Applied Biosystems) according to manufacture’s instructions with minor modification. Digested protein products were then dried, dissolved and labeled. The labeled striated and catch adductor samples were fractionated using an SCX column on an HPLC system (LC-20AB, Shimadzu, Japan). The fractionation procedure for each sample was performed as in the previous study [79]. Briefly, peptides were dissolved and then flowed into columns (1 mL/min), followed by elution using buffer A (25 mM NaH2PO4 and 25% ACN) and Buffer B (a mixture of Buffer A and 1 M KCl).

The peptides were injected and separated in a Nano-LC system, using a C18 analytical reverse-phase column (300 nL/min). The column was equilibrated with a gradient schedule of different concentrations (5, 45, 80, and 5%) of Solution B (a mixture of 95% acetonitrile and 0.1% formic acid [78]). A Triple TOF 5600 instrument was used for the mass spectrometry. All spectrometry data were collected and summarized using Bruker Daltonics micrOTOFcontrol. The raw data from the LC-MS/MS were transformed into MGF files and analyzed using the Data Analysis Software. The Mascot search engine was set to the following conditions, such as trypsin (digestion enzyme), cysteine carbamidomethylation (fixed modification), glutamine (pyroglutamic acid), and iTRAQ 8Plex on tyrosine (variable modification). The scanning range for MS was from 50 to 2000 m/z, with nitrogen as the collision gas (voltage, 1250 V; interface temperature, 150 °C). These experimental procedures were performed by Guangzhou Gene denovo Biotechnology Co., Ltd. (Guangzhou, China).

Functional analysis of differentially expressed proteins (DEPs) and its association analysis with transcriptome

To annotate the differentially expressed proteins, WEGO and Blast2GO were both used to search against the GO database according to the previous studies [78, 79]. KEGG pathway analysis was also performed to identify significantly enriched metabolic pathways in the striated or catch adductor muscles of scallops. A significance threshold of 0.05 was selected for the functional annotation of DEPs in the two muscles.

We used Pearson correlation to assess the relationships between mRNA and protein expression levels in the striated and catch muscles. The log2 transformation of the average fold changes was obtained for the transcriptome and proteome data. The transformed data was calculated and displayed as the scatter plots using R program, with the screening criterions (for mRNA, fold change > 2; for proteins, fold change > 1.2).

Abbreviations

BP:

Biological process

CC:

Cellular component

COG:

Clusters of Orthologous Groups of proteins

DEGs:

Differentially expressed genes

DEPs:

Differentially expressed proteins

DTT:

DL-Dithiothreitol

ECM:

Extracellular matrix

E-LC:

Myosin essential light chain

FDR:

False discovery rate

GO:

Gene Ontology

HPLC:

High performance liquid chromatography

iTRAQ:

Isobaric tags for relative and absolute quantitation

KEGG:

Kyoto Encyclopedia of Genes and Genomes

KOG:

euKaryotic Ortholog Groups

MF:

Molecular function

MHC:

Myosin heavy chain

Nr:

NCBI non-redundant protein

PMSF:

Phenylmethanesulfonyl fluoride

R-LC:

Myosin regulatory light chain

SCP:

Sarcoplasmic calcium-binding protein

SCX:

Strong cation exchange

SERCA:

Sarcoplasmic/endoplasmic reticulum calcium ATPase

SMTNL1:

Moothelin-like 1 protein

SR:

Sarcoplasmic reticulum

WEGO:

Web Gene Ontology Annotation Plot

References

  1. 1.

    Kelly AM, Rubinstein NA. Why are fetal muscles slow? Nature. 1980;288:266–9.

    Article  PubMed  CAS  Google Scholar 

  2. 2.

    Lefaucheur L, Ecolan P, Plantard L, Gueguen N. New insights into muscle fiber types in the pig. J Histochem Cytochem. 2002;50:719–30.

    Article  PubMed  CAS  Google Scholar 

  3. 3.

    Drexler HC, Ruhs A, Konzer A, Mendler L, Bruckskotten M, Looso M, et al. On marathons and sprints: an integrated quantitative proteomics and transcriptomics analysis of differences between slow and fast muscle fibers. Mol Cell Proteomics. 2012;50:M111–010801.

    Google Scholar 

  4. 4.

    Rakus D, Gizak A, Deshmukh A, Wiśniewski JR. Absolute quantitative profiling of the key metabolic pathways in slow and fast skeletal muscle. J Proteome Res. 2015;14:1400–11.

    Article  PubMed  CAS  Google Scholar 

  5. 5.

    Mykles DL. Heterogeneity of myofibrillar proteins in lobster fast and slow muscles: variants of troponin, paramyosin, and myosin light chains comprise four distinct protein assemblages. J Exp Zool. 1985;234:23–32.

    Article  PubMed  CAS  Google Scholar 

  6. 6.

    Mykles DL. Histochemical and biochemical characterization of two slow fiber types in decapod crustacean muscles. J Exp Zool. 1988;245:232–43.

    Article  PubMed  CAS  Google Scholar 

  7. 7.

    Cotton JL, Mykles DL. Cloning of a crustacean myosin heavy chain isoform: exclusive expression in fast muscle. J Exp Zool. 1993;267:578–86.

    Article  PubMed  CAS  Google Scholar 

  8. 8.

    Grimaldi A, Tettamanti G, Brivio MF, Valvassori R, Eguileor MD. Differentiation of slow and fast fibers in tentacles of sepia officinalis, (mollusca). Develop Growth Differ. 2004;46:181–93.

    Article  Google Scholar 

  9. 9.

    Chantler PD. Scallop adductor muscles: structure and function. Dev Aquac Fish Sci. 2006;35:229–316.

    Article  Google Scholar 

  10. 10.

    Chantler PD. Scallop adductor muscles: structure and function. In: Shumway SE, Parsons GJ, editors. Scallops: Biology, Ecology, Aquaculture, and Fisheries (Vol. 40). Amsterdam: Elsevier; 2016. p. 161–207.

    Google Scholar 

  11. 11.

    Szent-Györgyi AG, Chantler PD. Control of contraction by calcium binding to myosin. Myology. 1994;1:506–28.

    Google Scholar 

  12. 12.

    Geeves MA, Holmes KC. Structural mechanism of muscle contraction. Annu Rev Biochem. 1999;68:687–728.

    Article  PubMed  CAS  Google Scholar 

  13. 13.

    Hooper SL, Thuma JB. Invertebrate muscles: muscle specific genes and proteins. Physiol Rev. 2005;85:1001–60.

    Article  PubMed  CAS  Google Scholar 

  14. 14.

    Yang Y, Gourinath S, Kovacs M, Nyitray L, Reutzel R, Himmel DM, et al. Rigor-like structures from muscle myosins reveal key mechanical elements in the transduction pathways of this allosteric motor. Structure. 2007;15:553–64.

    Article  PubMed  CAS  Google Scholar 

  15. 15.

    Zhao FQ, Craig R. Millisecond time-resolved changes occurring in Ca2+-regulated myosin filaments upon relaxation. J Mol Biol. 2008;381:256–60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. 16.

    Zhao FQ, Craig R, Woodhead JL. Head-head interaction characterises the relaxed state of Limulus muscle myosin filaments. J Mol Biol. 2009;385:423–31.

    Article  PubMed  CAS  Google Scholar 

  17. 17.

    Woodhead JL, Zhao FQ, Craig R. Structural basis of the relaxed state of a Ca2+-regulated myosin filament and its evolutionary implications. Proc Natl Acad Sci U S A. 2013;110:8561–6.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Hu Z, Taylor DW, Reedy MK, Edwards RJ, Taylor KA. Structure of myosin filaments from relaxed Lethocerus flight muscle by cryo-EM at 6 Å resolution. Sci Adv. 2016;2:e1600058.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. 19.

    Mattisson AGM, Beechey RB. Some studies on cellular fractions of the adductor muscle of Pecten maximus. Exp Cell Res. 1966;41:227–43.

    Article  PubMed  CAS  Google Scholar 

  20. 20.

    Yamada A, Oiwa K. Myosin Mg-ATPase of molluscan muscles is slightly activated by F-actin under catch state in vitro. J Muscle Res Cell M. 2013;34:115–23.

    Article  CAS  Google Scholar 

  21. 21.

    Sulbarán G, Alamo L, Pinto A, Márquez G, Méndez F, Padrón R, et al. An invertebrate smooth muscle with striated muscle myosin filaments. Proc Natl Acad Sci U S A. 2015;112:5660–8.

    Article  CAS  Google Scholar 

  22. 22.

    Galler S. Molecular basis of the catch state in molluscan smooth muscles: a catchy challenge. J Muscle Res Cell M. 2008;29:7399.

    Google Scholar 

  23. 23.

    Galler S, Litzlbauer J, Kröss M, Grassberger H. The highly efficient holding function of the mollusc ‘catch’ muscle is not based on decelerated myosin head cross-bridge cycles. Proc R Soc B. 2010;277:803–8.

    Article  PubMed  Google Scholar 

  24. 24.

    Funabara D, Osawa R, Ueda M, Kanoh S, Hartshorne DJ, Watabe S. Myosin loop 2 is involved in the formation of a trimeric complex of twitchin, actin, and myosin. J Biol Chem. 2009;284:18015–20.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. 25.

    Vyatchin IG, Shevchenko UV, Lazarev SS, Matusovsky OS, Shelud'ko NS. Troponin-like regulation in muscle thin filaments of the mussel Crenomytilus grayanus (Bivalvia: Mytiloida). BBA-Proteins Proteom. 1854;2015:1444–50.

    Google Scholar 

  26. 26.

    Hooper SL, Hobbs KH, Thuma JB. Invertebrate muscles: thin and thick filament structure; molecular basis of contraction and its regulation, catch and asynchronous muscle. Prog Neurobiol. 2008;86:72–127.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. 27.

    Guderley HE, Tremblay I. Swimming in scallops. In: Shumway SE, Parsons GJ, editors. Scallops: Biology, Ecology, Aquaculture, and Fisheries (Vol. 40). Amsterdam: Elsevier; 2016. p. 535–66.

    Google Scholar 

  28. 28.

    Funabara D, Watabe S, Kanoh S. Phosphorylation properties of twitchin from yesso scallop catch and striated muscles. Fisheries Sci. 2015;81:1–10.

    Article  CAS  Google Scholar 

  29. 29.

    Sun X, Liu Z, Zhou L, Wu B, Dong Y, Yang A. Integration of next generation sequencing and EPR analysis to uncover molecular mechanism underlying shell color variation in scallops. PLoS One. 2016;11:e0161876.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. 30.

    Zhang X, Chen Y, Pan J, Liu X, Chen H, Zhou X, et al. Itraq-based quantitative proteomic analysis reveals the distinct early embryo myofiber type characteristics involved in landrace and miniature pig. BMC Genomics. 2016;17:1–10.

    Article  CAS  Google Scholar 

  31. 31.

    Sun X, Yang A, Wu B, Zhou L, Liu Z. Characterization of the mantle transcriptome of Yesso scallop (Patinopecten yessoensis): identification of genes potentially involved in biomineralization and pigmentation. PLoS One. 2015;10:e0122967.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. 32.

    Nyitray L, Jancso A, Ochiai Y, Graf L, Szent-Gyorgyi AG. Scallop striated and smooth muscle myosin heavy-chain isoforms are produced by alternative RNA splicing from a single gene. Proc Natl Acad Sci U S A. 1994;91:12686–90.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. 33.

    Hasegawa Y. Isolation of a cDNA encoding the motor domain of nonmuscle myosin which is specifically expressed in the mantle pallial cell layer of scallop (Patinopecten yessoensis). J Biochem. 2000;128:983–8.

    Article  PubMed  CAS  Google Scholar 

  34. 34.

    Hasegawa Y, Araki T. Identification of a novel unconventional myosin from scallop mantle tissue. J Biochem. 2002;131:113–9.

    Article  PubMed  CAS  Google Scholar 

  35. 35.

    Hasegawa Y, Ikeda Y. Cloning of a cDNA encoding the tail region of non-muscle myosin II from the mantle tissue of scallop Patinopecten yessoensis. Fisheries Sci. 2008;74:1201–3.

    Article  CAS  Google Scholar 

  36. 36.

    Ogut O, Granzier H, Jin JP. Acidic and basic troponin T isoforms in mature fast-twitch skeletal muscle and effect on contractility. Am J Phys. 1999;276:1162–70.

    Article  Google Scholar 

  37. 37.

    Perry SV. Activation of the contractile mechanism by calcium. In: Engel AG, Franzini-Armstrong C, editors. Myology, vol. 1. 3rd ed. New York: McGraw-Hill; 2004. p. 281–306.

    Google Scholar 

  38. 38.

    Yumoto F, Tanokura M. Structural and functional analysis of troponins from scallop striated and human cardiac muscles. Adv Exp Med Biol. 2007;592:163–73.

    Article  PubMed  Google Scholar 

  39. 39.

    Ojima T, Tanaka H, Nishita K. Cloning and sequence of a cDNA encoding Akazara scallop troponin C. Arch Biochem Biophys. 1994;311:272–6.

    Article  PubMed  CAS  Google Scholar 

  40. 40.

    Ojima T, Koizumi N, Ueyama K, Inoue A, Nishita K. Functional role of Ca2+-binding site IV of scallop troponin C. J Biochem. 2000;128:803–9.

    Article  PubMed  CAS  Google Scholar 

  41. 41.

    Tanaka H, Takeya Y, Doi T, Yumoto F, Tanokura M, Ohtsuki I, et al. Comparative studies on the functional roles of N- and C-terminal regions of molluskan and vertebrate troponin-I. FEBS J. 2005;272:4475–86.

    Article  PubMed  CAS  Google Scholar 

  42. 42.

    Labeit S, Lahmers S, Burkart C, Chi F, Mcnabb M, Witt S, et al. Expression of distinct classes of titin isoforms in striated and smooth muscles by alternative splicing, and their conserved interaction with filamins. J Mol Biol. 2006;362:664–81.

    Article  PubMed  CAS  Google Scholar 

  43. 43.

    Chi RJH. Smooth muscle titin interactions with alpha-actinin. The Florida State University, ProQuest Dissertations Publishing, 2007.

  44. 44.

    Granzier H, Labeit S. Structure-function relations of the giant elastic protein titin in striated and smooth muscle cells. Muscle Nerve. 2007;36:740–55.

    Article  PubMed  CAS  Google Scholar 

  45. 45.

    Winder SJ, Walsh MP. Calponin: thin filament-linked regulation of smooth muscle contraction. Cell Signal. 1993;5:677–86.

    Article  PubMed  CAS  Google Scholar 

  46. 46.

    Sirenko VV, Dobrzhanskaya AV, Shelud’Ko NS, Borovikov YS. Calponin-like protein from mussel smooth muscle is a competitive inhibitor of actomyosin ATPase. Biochemistry Biokhimiia. 2016;81:28–33.

    Article  PubMed  CAS  Google Scholar 

  47. 47.

    Borman MA, Freed TA, Haystead TAJ, Macdonald JA. The role of the calponin homology domain of smoothelin-like 1 (smtnl1) in myosin phosphatase inhibition and smooth muscle contraction. Mol Cell Biochem. 2009;327:93–100.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. 48.

    Takeshi ENDO, Masaki T. Molecular properties and functions in vitro of chicken smooth-muscle α-actinin in comparison with those of striated-muscle α-actinins. J Biochem. 1982;92:1457–68.

    Article  Google Scholar 

  49. 49.

    Sjöblom B, Salmazo A, Djinović-Carugo K. α-Actinin structure and regulation. Cell Mol Life Sci. 2008;65:2688–701.

    Article  PubMed  CAS  Google Scholar 

  50. 50.

    Atkinson RA, Joseph C, Dal PF, Birolo L, Stier G, Pucci P, et al. Binding of alpha-actinin to titin: implications for z-disk assembly. Biochemistry. 2000;39:5255–64.

    Article  PubMed  CAS  Google Scholar 

  51. 51.

    Barone V, Randazzo D, Re VD, Sorrentino V, Rossi D. Organization of junctional sarcoplasmic reticulum proteins in skeletal muscle fibers. J Muscle Res Cell M. 2015;36:1–15.

    Article  Google Scholar 

  52. 52.

    Sanger JW, Sanger JM. Sarcoplasmic reticulum in the adductor muscles of a Bermuda scallop: comparison of smooth versus cross-striated portions. Biol Bull. 1985;168:447–60.

    Article  Google Scholar 

  53. 53.

    Bönnemann CG, Thompson TG, Pf VDV, Goebel HH, Warlo I, Vollmers B, et al. Filamin c accumulation is a strong but nonspecific immunohistochemical marker of core formation in muscle. J Neurol Sci. 2003;206:71–8.

    Article  PubMed  Google Scholar 

  54. 54.

    Small JV, Fürst DO, De MJ. Localization of filamin in smooth muscle. J Cell Biol. 1986;102:210–20.

    Article  PubMed  CAS  Google Scholar 

  55. 55.

    Castellani L, Cohen C. Myosin rod phosphorylation and the catch state of molluscan muscles. Science. 1987;235:334–7.

    Article  PubMed  CAS  Google Scholar 

  56. 56.

    Johnson WH, Kahn JS, Szent-Gyorgyi AG. Paramyosin and contraction of catch muscles. Science. 1959;130:160–1.

    Article  PubMed  CAS  Google Scholar 

  57. 57.

    Chen GX, Tan RY, Gong ZX, Huang YP, Wang SZ, Cao TG. Paramyosin and the catch mechanism. Biophys Chem. 1988;29:147–53.

    Article  PubMed  CAS  Google Scholar 

  58. 58.

    Baguet F, Gillis JM. Energy cost of tonic contraction in lamellibranch catch muscle. J Physiol. 1968;198:127–43.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  59. 59.

    Shelud’ko NS, Matusovsky OS, Permyakova TV, Matusovskaya GG. ‘Twitchin-actin linkage hypothesis’ for the catch mechanism in molluscan muscles: evidence that twitchin interacts with myosin, myorod, and the paramyosin core and affects properties of actomyosin. Archiv Biochem Biophys. 2007;466:125–35.

    Article  CAS  Google Scholar 

  60. 60.

    Butler TM, Mooers SU, Siegman MJ. Catch force links and the low to high force transition of myosin. Biophys J. 2006;90:3193–202.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  61. 61.

    Butler TM, Mooers SU, Narayan SR, Siegman MJ. The N-terminal region of twichin binds thick and thin contractile filaments: redundant mechanisms of catch force maintenance. J Biol Chem. 2010;285:40654–65.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  62. 62.

    Yamada A, Yoshio M, Oiwa K, Nyitray L. Catchin, a novel protein in molluscan catch muscles, is produced by alternative splicing from the myosin heavy chain gene. J Mol Biol. 2000;295:169–78.

    Article  PubMed  CAS  Google Scholar 

  63. 63.

    Vibert P, Edelstein SM, Castellani L, Elliott BW. Mini-titins in striated and smooth molluscan muscles: structure, location and immunological crossreactivity. J Muscle Res Cell M. 1993;14:598–607.

    Article  CAS  Google Scholar 

  64. 64.

    Andersen O, Torgersen JS, Pagander HH, Magnesen T, Johnston IA. Gene expression analyses of essential catch factors in the smooth and striated adductor muscles of larval, juvenile and adult great scallop (Pecten maximus). J Muscle Res Cell M. 2009;30:233–42.

    Article  CAS  Google Scholar 

  65. 65.

    Zwaan AD, Thompson RJ, Livingstone DR. Physiological and biochemical aspects of the valve snap and valve closure responses in the giant scallop Placopecten magellanicus. J Comp Physiol B. 1980;137:97–104.

    Article  Google Scholar 

  66. 66.

    Tremblay I, Guderley HE. Scallops show that muscle metabolic capacities reflect locomotor style and morphology. Physiol Biochem Zool. 2014;87:231–44.

    Article  PubMed  CAS  Google Scholar 

  67. 67.

    Storey KB. Effects of arginine phosphate and octopine on glycolytic enzyme activities from Sepia officinalis mantle muscle. J Comp Physiol B. 1981;142:501–7.

    Article  CAS  Google Scholar 

  68. 68.

    Tremblay I, Guderley HE, Himmelman JH. Swimming away or clamming up: the use of phasic and tonic adductor muscles during escape responses varies with shell morphology in scallops. J Exp Biol. 2012;215:4131–43.

    Article  PubMed  Google Scholar 

  69. 69.

    Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29:644.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. 70.

    Lobo. Basic local alignment search tool (blast). J Mol Biol. 2008;215:403–10.

    Google Scholar 

  71. 71.

    Pruitt KD, Tatusova T, Maglott DR. NCBI reference sequence (refseq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005;33:D501–4.

    Article  PubMed  CAS  Google Scholar 

  72. 72.

    The UniProt Consortium. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45:D158–69.

    Article  CAS  Google Scholar 

  73. 73.

    Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:d480–4.

    Article  PubMed  CAS  Google Scholar 

  74. 74.

    Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, et al. The COG database: an updated version includes eukaryotes. BMC Bioinformatics. 2003;4:41.

    Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  76. 76.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    Article  PubMed  CAS  Google Scholar 

  77. 77.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  78. 78.

    Qin J, Gu F, Liu D, Yin C, Zhao S, Chen H, et al. Proteomic analysis of elite soybean jidou17 and its parents using itraq-based quantitative approaches. Proteome Sci. 2013;11:334–45.

    Article  CAS  Google Scholar 

  79. 79.

    Song H, Wang HY, Zhang T. Comprehensive and quantitative proteomic analysis of metamorphosis-related proteins in the veined rapa whelk, Rapana venosa. Int J Mol Sci. 2016;17:924.

    Article  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank three anonymous reviewers for their extensive and constructive comments on the manuscript. We also would like to thank Dennis Hedgecock and Xiaoshen Yin from University of Southern California for the assistance in preparing the revision of this manuscript. We are grateful to the Guangzhou Gene denovo Biotechnology Co., Ltd. for assisting in the iTRAQ and RNA sequencing.

Funding

This study was supported by research grants from the National Natural Science Foundation of China (31602153), Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, China P. R. China (2016LMFS-B02), Natural Science Foundation of Shandong Province (ZR2016CQ32) and the Key Laboratory of Mariculture, Ministry of Education, Ocean University of China (KLM2018004). The funding bodies did not contribute to the design of the study, sample collection, analysis, interpretation of data, or in writing the manuscript.

Availability of data and materials

For transcriptome sequencing, the raw reads have been deposited in the NCBI SRA database (SRR4254476). For iTRAQ proteomics, the mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium with the dataset identifier PXD005166. Other data sets and results generated in this project are included within the article and its additional files.

Author information

Affiliations

  1. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao, 266071, China

    Xiujun Sun, Zhihong Liu, Biao Wu, Liqing Zhou, Wei Wu & Aiguo Yang

  2. Function Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, Qingdao, 266200, People’s Republic of China

    Xiujun Sun, Zhihong Liu, Biao Wu, Liqing Zhou, Wei Wu & Aiguo Yang

  3. College of Fisheries, Ocean University of China, Qingdao, 266003, China

    Qi Wang

Authors
  1. Xiujun Sun
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zhihong Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Biao Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Liqing Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Qi Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Wei Wu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Aiguo Yang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

XS, ZL and AY conceived the project and designed the experiment. XS and BW performed the experiment. LZ, QW, and WW analyzed the data. AY, ZL, BW, LZ, QW and WW made substantial contributions to the interpretation of data. XS wrote the draft manuscript. All authors contributed to the production of the revised manuscript. All authors have read and approved the final manuscript for publication, and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Corresponding author

Correspondence to Aiguo Yang.

Ethics declarations

Ethics approval and consent to participate

The scallops used in the current study were marine-cultured animals. All the experiments were conducted following institutional and national guidelines. Since no endangered or protected species was involved in the experiments of this study, no specific permission was required for the animal material.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. The identification and annotation results of differentially expressed proteins (DEPs) between the striated and catch adductor muscles by the iTRAQ-based quantitative proteomics. There are a total of 474 DEPs identified between the two muscles, which include 198 upregulated and 276 downregulated DEPs. The functional annotation for these DEPs is associated with significantly enriched GO terms and KEGG pathways. (XLSX 135 kb)

Additional file 2:

Table S2. The most enriched genes and proteins in the striated and catch adductor muscles of Yesso scallop Patinopecten yessoensis. List of the selected unigenes is divided into four categories, including muscle proteins, metabolism related enzymes, calcium signaling, membrane and extracellular proteins. The related information on quadrant (see Fig. 5), protein annotation, and false discovery rate (FDR) for these selected unigenes is summarized in this table. (DOC 60 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sun, X., Liu, Z., Wu, B. et al. Differences between fast and slow muscles in scallops revealed through proteomics and transcriptomics. BMC Genomics 19, 377 (2018). https://doi.org/10.1186/s12864-018-4770-2

Download citation

Keywords

  • Scallops
  • Striated adductor muscle
  • Catch adductor muscle
  • Proteomics
  • Transcriptomics

BMC Genomics

ISSN: 1471-2164

Contact us