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 scallopPatinopecten 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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4770-2) contains supplementary material, which is available to authorized users.

Keywords: Scallops, Striated adductor muscle, Catch adductor muscle, Proteomics, Transcriptomics

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.

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.

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.

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.

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).

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.

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 Ca²⁺-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 scallopP. yessoensis in previous studies [33–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 Ca²⁺ regulatory system of vertebrate striated muscles [36–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–50]. It is suggested that the enrichment of actinin protein in striated muscle ofP. 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 Ca²⁺, 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–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 inP. 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–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.

Biological material

Two-year-old live individuals ofP. 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 withIsochrysis 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 theP < 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 NaH₂PO4 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).

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.

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.

• 1.Kelly AM, Rubinstein NA. Why are fetal muscles slow? Nature. 1980;288:266–269. doi: 10.1038/288266a0. [DOI] [PubMed] [Google Scholar]
• 2.Lefaucheur L, Ecolan P, Plantard L, Gueguen N. New insights into muscle fiber types in the pig. J Histochem Cytochem. 2002;50:719–730. doi: 10.1177/002215540205000513. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1074/mcp.M111.010801. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–1411. doi: 10.1021/pr5010357. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1002/jez.1402340105. [DOI] [PubMed] [Google Scholar]
• 6.Mykles DL. Histochemical and biochemical characterization of two slow fiber types in decapod crustacean muscles. J Exp Zool. 1988;245:232–243. doi: 10.1002/jez.1402450303. [DOI] [PubMed] [Google Scholar]
• 7.Cotton JL, Mykles DL. Cloning of a crustacean myosin heavy chain isoform: exclusive expression in fast muscle. J Exp Zool. 1993;267:578–586. doi: 10.1002/jez.1402670605. [DOI] [PubMed] [Google Scholar]
• 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–193. doi: 10.1111/j.1440-169X.2004.00738.x. [DOI] [PubMed] [Google Scholar]
• 9.Chantler PD. Scallop adductor muscles: structure and function. Dev Aquac Fish Sci. 2006;35:229–316. doi: 10.1016/S0167-9309(06)80031-1. [DOI] [Google Scholar]
• 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. pp. 161–207. [Google Scholar]
• 11.Szent-Györgyi AG, Chantler PD. Control of contraction by calcium binding to myosin. Myology. 1994;1:506–528. [Google Scholar]
• 12.Geeves MA, Holmes KC. Structural mechanism of muscle contraction. Annu Rev Biochem. 1999;68:687–728. doi: 10.1146/annurev.biochem.68.1.687. [DOI] [PubMed] [Google Scholar]
• 13.Hooper SL, Thuma JB. Invertebrate muscles: muscle specific genes and proteins. Physiol Rev. 2005;85:1001–1060. doi: 10.1152/physrev.00019.2004. [DOI] [PubMed] [Google Scholar]
• 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–564. doi: 10.1016/j.str.2007.03.010. [DOI] [PubMed] [Google Scholar]
• 15.Zhao FQ, Craig R. Millisecond time-resolved changes occurring in Ca2+-regulated myosin filaments upon relaxation. J Mol Biol. 2008;381:256–260. doi: 10.1016/j.jmb.2008.06.032. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–431. doi: 10.1016/j.jmb.2008.10.038. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–8566. doi: 10.1073/pnas.1218462110. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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. doi: 10.1126/sciadv.1600058. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 19.Mattisson AGM, Beechey RB. Some studies on cellular fractions of the adductor muscle of Pecten maximus. Exp Cell Res. 1966;41:227–243. doi: 10.1016/S0014-4827(66)80130-1. [DOI] [PubMed] [Google Scholar]
• 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–123. doi: 10.1007/s10974-013-9339-8. [DOI] [PubMed] [Google Scholar]
• 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–5668. doi: 10.1073/pnas.1513439112. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 22.Galler S. Molecular basis of the catch state in molluscan smooth muscles: a catchy challenge. J Muscle Res Cell M. 2008;29:7399. doi: 10.1007/s10974-008-9149-6. [DOI] [PubMed] [Google Scholar]
• 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–808. doi: 10.1098/rspb.2009.1618. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–18020. doi: 10.1074/jbc.M109.016485. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–1450. doi: 10.1016/j.bbapap.2015.07.010. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1016/j.pneurobio.2008.06.004. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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. pp. 535–566. [Google Scholar]
• 28.Funabara D, Watabe S, Kanoh S. Phosphorylation properties of twitchin from yesso scallop catch and striated muscles. Fisheries Sci. 2015;81:1–10. doi: 10.1007/s12562-014-0815-x. [DOI] [Google Scholar]
• 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. doi: 10.1371/journal.pone.0161876. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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. doi: 10.1186/s12864-015-2294-6. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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. doi: 10.1371/journal.pone.0122967. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–12690. doi: 10.1073/pnas.91.26.12686. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–988. doi: 10.1093/oxfordjournals.jbchem.a022850. [DOI] [PubMed] [Google Scholar]
• 34.Hasegawa Y, Araki T. Identification of a novel unconventional myosin from scallop mantle tissue. J Biochem. 2002;131:113–119. doi: 10.1093/oxfordjournals.jbchem.a003065. [DOI] [PubMed] [Google Scholar]
• 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–1203. doi: 10.1111/j.1444-2906.2008.01642.x. [DOI] [Google Scholar]
• 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–1170. doi: 10.1152/ajpcell.1999.276.5.C1162. [DOI] [PubMed] [Google Scholar]
• 37.Perry SV. Activation of the contractile mechanism by calcium. In: Engel AG, Franzini-Armstrong C, editors. Myology, vol. 1. 3. New York: McGraw-Hill; 2004. pp. 281–306. [Google Scholar]
• 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–173. doi: 10.1007/978-4-431-38453-3_15. [DOI] [PubMed] [Google Scholar]
• 39.Ojima T, Tanaka H, Nishita K. Cloning and sequence of a cDNA encoding Akazara scallop troponin C. Arch Biochem Biophys. 1994;311:272–276. doi: 10.1006/abbi.1994.1237. [DOI] [PubMed] [Google Scholar]
• 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–809. doi: 10.1093/oxfordjournals.jbchem.a022818. [DOI] [PubMed] [Google Scholar]
• 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–4486. doi: 10.1111/j.1742-4658.2005.04866.x. [DOI] [PubMed] [Google Scholar]
• 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–681. doi: 10.1016/j.jmb.2006.07.077. [DOI] [PubMed] [Google Scholar]
• 43.Chi RJH. Smooth muscle titin interactions with alpha-actinin. The Florida State University, ProQuest Dissertations Publishing, 2007.
• 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–755. doi: 10.1002/mus.20886. [DOI] [PubMed] [Google Scholar]
• 45.Winder SJ, Walsh MP. Calponin: thin filament-linked regulation of smooth muscle contraction. Cell Signal. 1993;5:677–686. doi: 10.1016/0898-6568(93)90029-L. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1134/S000629791601003X. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1007/s11010-009-0047-z. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–1468. doi: 10.1093/oxfordjournals.jbchem.a134070. [DOI] [PubMed] [Google Scholar]
• 49.Sjöblom B, Salmazo A, Djinović-Carugo K. α-Actinin structure and regulation. Cell Mol Life Sci. 2008;65:2688–2701. doi: 10.1007/s00018-008-8080-8. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–5264. doi: 10.1021/bi991891u. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1007/s10974-014-9401-1. [DOI] [PubMed] [Google Scholar]
• 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–460. doi: 10.2307/1541524. [DOI] [Google Scholar]
• 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–78. doi: 10.1016/S0022-510X(02)00341-6. [DOI] [PubMed] [Google Scholar]
• 54.Small JV, Fürst DO, De MJ. Localization of filamin in smooth muscle. J Cell Biol. 1986;102:210–220. doi: 10.1083/jcb.102.1.210. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 55.Castellani L, Cohen C. Myosin rod phosphorylation and the catch state of molluscan muscles. Science. 1987;235:334–337. doi: 10.1126/science.3026049. [DOI] [PubMed] [Google Scholar]
• 56.Johnson WH, Kahn JS, Szent-Gyorgyi AG. Paramyosin and contraction of catch muscles. Science. 1959;130:160–161. doi: 10.1126/science.130.3368.160. [DOI] [PubMed] [Google Scholar]
• 57.Chen GX, Tan RY, Gong ZX, Huang YP, Wang SZ, Cao TG. Paramyosin and the catch mechanism. Biophys Chem. 1988;29:147–153. doi: 10.1016/0301-4622(88)87034-0. [DOI] [PubMed] [Google Scholar]
• 58.Baguet F, Gillis JM. Energy cost of tonic contraction in lamellibranch catch muscle. J Physiol. 1968;198:127–143. doi: 10.1113/jphysiol.1968.sp008597. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–135. doi: 10.1016/j.abb.2007.07.014. [DOI] [PubMed] [Google Scholar]
• 60.Butler TM, Mooers SU, Siegman MJ. Catch force links and the low to high force transition of myosin. Biophys J. 2006;90:3193–3202. doi: 10.1529/biophysj.105.077453. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–40665. doi: 10.1074/jbc.M110.166041. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–178. doi: 10.1006/jmbi.1999.3349. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1007/BF00141557. [DOI] [PubMed] [Google Scholar]
• 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–242. doi: 10.1007/s10974-009-9192-y. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1007/BF00689208. [DOI] [Google Scholar]
• 66.Tremblay I, Guderley HE. Scallops show that muscle metabolic capacities reflect locomotor style and morphology. Physiol Biochem Zool. 2014;87:231–244. doi: 10.1086/674107. [DOI] [PubMed] [Google Scholar]
• 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–507. doi: 10.1007/BF00688982. [DOI] [Google Scholar]
• 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–4143. doi: 10.1242/jeb.075986. [DOI] [PubMed] [Google Scholar]
• 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. doi: 10.1038/nbt.1883. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 70.Lobo Basic local alignment search tool (blast) J Mol Biol. 2008;215:403–410. doi: 10.1016/S0022-2836(05)80360-2. [DOI] [PubMed] [Google Scholar]
• 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–D504. doi: 10.1093/nar/gki025. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 72.The UniProt Consortium UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2017;45:D158–D169. doi: 10.1093/nar/gkw1099. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–d484. doi: 10.1093/nar/gkm882. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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. doi: 10.1186/1471-2105-4-41. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 75.Langmead B, Salzberg S. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–359. doi: 10.1038/nmeth.1923. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 76.Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–140. doi: 10.1093/bioinformatics/btp616. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–29. doi: 10.1038/75556. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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–345. doi: 10.1186/1477-5956-11-12. [DOI] [PMC free article] [PubMed] [Google Scholar]
• 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. doi: 10.3390/ijms17060924. [DOI] [PMC free article] [PubMed] [Google Scholar]

Supplementary Materials

Data Availability Statement

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.