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

Temperature is an important environmental factor affecting metabolism, growth and life activities of marine bivalves. To reveal the molecular and metabolic responses to heat stress, comparative transcriptomics and metabolomics were performed in the commercially important bivalve, Manila clam Ruditapes philippinarum. Comparative transcriptomics revealed a total of 3980 differential genes (DEGs) mainly involved in energy metabolism, protein processing, glycerolphospholipid metabolism, inoxidizability and anti-apoptosis. Comparative metabolomic analysis identified 293 differential metabolites (DEMs) in the heat-stressed clams, including Acetyl-CoA, Beta-D-Glucose, phosphatidylcholine and free amino acids. For the heat-stress clams, they showed the significantly increased enzyme activities of AKP and ACP, as well as antioxidant-related enzyme activities (SOD and CAT) and oxidative damage substance (MDA). The changes in gene expression of glycolysis-related genes (GLUT1, GLK, HK1 and GPI6) and metabolites (Beta-D-Glucose and Acetyl-CoA) reflect the adjusted energy metabolism in the heat-stressed clams. The protein processing in endoplasmic reticulum may play the important roles in the stability of cell membrane structure during heat adaptation. The increased contents of glycerophospholipid (PC) may be essential to maintain the stability and permeability of cell membrane in the heat-stressed clams. Molecular chaperones (HSP70 and HSP90) may participate in the degradation of misfolded proteins to maintain cell homeostasis during the heat stress. To defense the heat stress, clams also initiate the protection and defense mechanisms by activating antioxidant system, innate immune system and anti-apoptotic system. These findings will shed new lights on molecular and metabolic adaptation to heat stress in the intertidal mollusks in a warming climate at the global scale.

The Manila clam (Ruditapes philippinarum) is one of the most commercially important bivalves along the coast of China. With the continuous expansion of clam farming scale, it may lead to some serious problems, including loss of genetic variation, inbreeding depression, and reduced effective population size (N_e). In the present study, eleven microsatellite markers were used to investigate the genetic diversity and differentiation among 13 clam populations along the coast of China. As a result, 150 alleles were detected according to the genotyping results of eleven microsatellite loci. The observed heterozygosity (H_o) was estimated to be ranging from 0.437 to 0.678, while the expected heterozygosity (H_e) was calculated to be varying from 0.587 to 0.700. F_st values between populations ranged from 0.0046-0.1983. In particular, the Laizhou population had the highest genetic variability, which was significantly different from the others (all F_st values > 0.1). For all the clam populations, there was no significant linear regression between genetic and geographic distance, indicating that these populations do not follow a pattern of isolation by distance (IBD). Genetic structure was estimated according to NJ, principal coordinates (PCoA), and structure-based clustering. Estimates of effective population size range from dozens to thousands among different populations, based on linkage-disequilibrium and molecular coancestry methods. The results reveal the genetic diversity of clams and verify the hypothesis that clam population differentiation may be influenced by the mode of southern breeding and northern culture, providing guiding information for natural resource conservation and genetic breeding of clams.

China has the largest production yield of Manila clam Ruditapes philippinarum in the world. Most of the clam seeds for aquaculture are mainly derived from artificial breeding in southern China, likely resulting in the loss of genetic variation and inbreeding depression. To understand the genetic and haplotype diversity of R. philippinarum, 14 clam populations sampled from different regions of China were analyzed by three molecular markers, including COI, 16SrRNA and ITS. Based on the results of the COI and ITS genes, the 14 populations showed a moderate to high level of genetic diversity, with an average haplotype diversity of 0.9242 and nucleotide diversity of 0.05248. AMOVA showed that there was significant genetic differentiation among all populations (mean F_ST of the total population was 0.4534). Pairwise F_ST analysis showed that genetic differentiation reached significant levels between Laizhou and other populations. Two Laizhou populations showed great divergence from other populations, forming an independent branch in the phylogenetic tree. The shared haplotypes Hap_2 and Hap_4 of COI appeared most frequently in most clam populations. In contrast, 16SrRNA analysis of the clam populations revealed the dominated haplotype Hap_2, accounting for 70% of the total number of individuals. The haplotype diversity of the Laizhou population (Laizhou shell-wide (KK) and Laizhou dock (LZMT)) was relatively higher than other populations, showing multiple unique haplotypes (e.g., Hap_40, Hap_41 and Hap_42). These findings of genetic and haplotype diversity of clam populations provide guiding information for genetic resource conservation and genetic improvement of the commercially important R. philippinarum.

Glycogen, a stored form of glucose, is an important form of energy for aquatic shellfish, contributing to the flavor and quality of the oyster. The glycogen synthase (GYS) and glycogen synthase kinase 3β (GSK3β) are two major enzymes in the glycogenesis. However, the information of the two genes in the Jinjiang oyster Crassostrea ariakensis remains limited. In this study, we identified the genes of GYS and GSK3β and further explored their function in the glycogen synthesis of C. ariakensis. The GYS and GSK3β were distributed in all tested tissues, and high expression of GYS and glycogen content were detected in the gonad, labial palp, hepatopancreas, and mantle, while the high expression of GSK3β was observed in the gill and adductor muscle. The expression of GYS was positively correlated with the glycogen content, while GSK3β was negatively correlated. Additionally, knockdown of GSK3β using RNAi decreased the GYS expression, revealing the negative regulatory effect of GSK3β on GYS. These findings enrich the research data of GSK3β and GYS involved in glycogen synthesis, providing valuable information for further research on the function of GSK3β and GYS in the glycogen synthesis process of oyster.