Genome-wide investigation of the TGF-ß superfamily in scallops.

BACKGROUND: Transforming growth factor ß (TGF-ß) superfamily genes can regulate various processes, especially in embryogenesis, adult development, and homeostasis. To understand the evolution and divergence patterns of the TGF-ß superfamily in scallops, genome-wide data from the Bay scallop (Argopecten irradians), the Zhikong scallop (Chlamys farreri) and the Yesso scallop (Mizuhopecten yessoensis) were systematically analysed using bioinformatics methods. RESULTS: Twelve members of the TGF-ß superfamily were identified for each scallop. The phylogenetic tree showed that these genes were grouped into 11 clusters, including BMPs, ADMP, NODAL, GDF, activin/inhibin and AMH. The number of exons and the conserved motif showed some differences between different clusters, while genes in the same cluster exhibited high similarity. Selective pressure analysis revealed that the TGF-ß superfamily in scallops was evolutionarily conserved. The spatiotemporal expression profiles suggested that different TGF-ß members have distinct functions. Several BMP-like and NODAL-like genes were highly expressed in early developmental stages, patterning the embryonic body plan. GDF8/11-like genes showed high expression in striated muscle and smooth muscle, suggesting that these genes may play a critical role in regulating muscle growth. Further analysis revealed a possible duplication of AMH, which played a key role in gonadal growth/maturation in scallops. In addition, this study found that several genes were involved in heat and hypoxia stress in scallops, providing new insights into the function of the TGF-ß superfamily. CONCLUSION: Characteristics of the TGF-ß superfamily in scallops were identified, including sequence structure, phylogenetic relationships, and selection pressure. The expression profiles of these genes in different tissues, at different developmental stages and under different stresses were investigated. Generally, the current study lays a foundation for further study of their pleiotropic biological functions in scallops.

Functional analyses of TRAF6 gene in Argopecten scallops.

The tumor necrosis factor (TNF) receptor-associated factor (TRAF) family has been reported to be involved in many immune pathways. In a previous study, we identified 5 TRAF genes, including TRAF2, 3, 4, 6, and 7, in the bay scallop (Argopecten irradians, Air) and the Peruvian scallop (Argopecten purpuratus, Apu). Since TRAF6 is a key molecular link in the TNF superfamily, we conducted a series of studies targeting the TRAF6 gene in the Air and Apu scallops as well as their hybrid progeny, Aip (Air ♀ × Apu ♂) and Api (Apu ♀ × Air ♂). Subcellular localization assay showed that the Air-, Aip-, and Api-TRAF6 were widely distributed in the cytoplasm of the human embryonic kidney cell line (HEK293T). Additionally, dual-luciferase reporter assay revealed that among TRAF3, TRAF4, and TRAF6, only the overexpression of TRAF6 significantly activated NF-κB activity in the HEK293T cells in a dose-dependent manner. These results suggest a crucial role of TRAF6 in the immune response in Argopecten scallops. To investigate the specific immune mechanism of TRAF6 in Argopecten scallops, we conducted TRAF6 knockdown using RNA interference. Transcriptomic analyses of the TRAF6 RNAi and control groups identified 1194, 2403, and 1099 differentially expressed genes (DEGs) in the Air, Aip, and Api scallops, respectively. KEGG enrichment analyses revealed that these DEGs were primarily enriched in transport and catabolism, amino acid metabolism, peroxisome, lysosome, and phagosome pathways. Expression profiles of 28 key DEGs were confirmed by qRT-PCR assays. The results of this study may provide insights into the immune mechanisms of TRAF in Argopecten scallops and ultimately benefit scallop breeding.

Chromosome-level genome assembly of sea scallop Placopecten magellanicus provides insights into the genetic characteristics and adaptive evolution of large scallops.

Placopecten magellanicus (Gmelin, 1791), a deep-sea Atlantic scallop, holds significant commercial value as a benthic marine bivalve along the northwest Atlantic coast. Recognizing its economic importance, the need to reconstruct its genome assembly becomes apparent, fostering insights into natural resources and generic breeding potential. This study reports a high-quality chromosome-level genome of P. magellanicus, achieved through the integration of Illumina short read sequencing, PacBio HiFi sequencing, and Hi-C sequencing techniques. The resulting assembly spans 1778 Mb with a scaffold N50 of 86.71 Mb. An intriguing observation arises - the genome size of P. magellanicus surpasses that of its Pectinidae family peers by 1.80 to 2.46 times. Within this genome, 28,111 protein-coding genes were identified. Comparative genomic analysis involving five scallop species unveils the critical determinant of this expanded genome: the proliferation of repetitive sequences recently inserted, contributing to its enlarged size. The landscape of whole genome collinearity sheds light on the relationships among scallop species, enhancing our broader understanding of their genomic framework. This genome provides genomic resources for future molecular biology research on scallops and serves as a guide for the exploration of longevity-related genes in scallops.

Functional Analysis of ß-Carotene Oxygenase 2 (BCO2) Gene in Yesso Scallop (Patinopecten yessoensis).

Carotenoids are essential nutrients for humans and animals, and carotenoid coloration represents an important meat quality parameter for many farmed animals. Increasingly, studies have demonstrated that vertebrate carotenoid cleavage oxygenases (CCOs) are essential enzymes in carotenoid metabolism and are therefore potential candidate genes for improving carotenoid deposition. However, our understanding of carotenoid bioavailability and CCOs functions in invertebrates, particularly marine species, is currently quite limited. We previously identified that a CCO homolog, PyBCO-like 1, was the causal gene for carotenoid coloration in the 'Haida golden scallop', a variety of Yesso scallop (Patinopecten yessoensis) characterized by carotenoid enrichment. Here, we found that another CCO-encoding gene named PyBCO2 (ß-carotene oxygenase 2) was widely expressed in P. yessoensis organs/tissues, with the highest expression in striated muscle. Inhibiting BCO2 expression in P. yessoensis through RNA interference led to increased carotenoid (pectenolone and pectenoxanthin) deposition in the striated muscle, and the color of the striated muscle changed from white to light orange. Our results indicate that PyBCO2 might be a candidate gene used for improving carotenoid content in normal Yesso scallops, and also in 'Haida golden scallops'.

Somatostatin Receptor Gene Functions in Growth Regulation in Bivalve Scallop and Clam.

Bivalves hold an important role in marine aquaculture and the identification of growth-related genes in bivalves could contribute to a better understanding of the mechanism governing their growth, which may benefit high-yielding bivalve breeding. Somatostatin receptor (SSTR) is a conserved negative regulator of growth in vertebrates. Although SSTR genes have been identified in invertebrates, their involvement in growth regulation remains unclear. Here, we identified seven SSTRs (PySSTRs) in the Yesso scallop, Patinopecten yessoensis, which is an economically important bivalve cultured in East Asia. Among the three PySSTRs (PySSTR-1, -2, and -3) expressed in adult tissues, PySSTR-1 showed significantly lower expression in fast-growing scallops than in slow-growing scallops. Then, the function of this gene in growth regulation was evaluated in dwarf surf clams (Mulinia lateralis), a potential model bivalve cultured in the lab, via RNA interference (RNAi) through feeding the clams Escherichia coli containing plasmids expressing double-stranded RNAs (dsRNAs) targeting MlSSTR-1. Suppressing the expression of MlSSTR-1, the homolog of PySSTR-1 in M. lateralis, resulted in a significant increase in shell length, shell width, shell height, soft tissue weight, and muscle weight by 20%, 22%, 20%, 79%, and 92%, respectively. A transcriptome analysis indicated that the up-regulated genes after MlSSTR-1 expression inhibition were significantly enriched in the fat digestion and absorption pathway and the insulin pathway. In summary, we systemically identified the SSTR genes in P. yessoensis and revealed the growth-inhibitory role of SSTR-1 in bivalves. This study indicates the conserved function of somatostatin signaling in growth regulation, and ingesting dsRNA-expressing bacteria is a useful way to verify gene function in bivalves. SSTR-1 is a candidate target for gene editing in bivalves to promote growth and could be used in the breeding of fast-growing bivalves.

The Glucose-Succinate Pathway: A Crucial Anaerobic Metabolic Pathway in the Scallop Chlamys farreri Experiencing Heat Stress.

Recently, the increase in marine temperatures has become an important global marine environmental issue. The ability of energy supply in marine animals plays a crucial role in avoiding the stress of elevated temperatures. The investigation into anaerobic metabolism, an essential mechanism for regulating energy provision under heat stress, is limited in mollusks. In this study, key enzymes of four anaerobic metabolic pathways were identified in the genome of scallop Chlamys farreri, respectively including five opine dehydrogenases (CfOpDHs), two aspartate aminotransferases (CfASTs) divided into cytoplasmic (CfAST1) and mitochondrial subtype (CfAST2), and two phosphoenolpyruvate carboxykinases (CfPEPCKs) divided into a primitive type (CfPEPCK2) and a cytoplasmic subtype (CfPEPCK1). It was surprising that lactate dehydrogenase (LDH), a key enzyme in the anaerobic metabolism of the glucose-lactate pathway in vertebrates, was absent in the genome of scallops. Phylogenetic analysis verified that CfOpDHs clustered according to the phylogenetic relationships of the organisms rather than substrate specificity. Furthermore, CfOpDHs, CfASTs, and CfPEPCKs displayed distinct expression patterns throughout the developmental process and showed a prominent expression in muscle, foot, kidney, male gonad, and ganglia tissues. Notably, CfASTs displayed the highest level of expression among these genes during the developmental process and in adult tissues. Under heat stress, the expression of CfASTs exhibited a general downregulation trend in the six tissues examined. The expression of CfOpDHs also displayed a downregulation trend in most tissues, except CfOpDH1/3 in striated muscle showing significant up-regulation at some time points. Remarkably, CfPEPCK1 was significantly upregulated in all six tested tissues at almost all time points. Therefore, we speculated that the glucose-succinate pathway, catalyzed by CfPEPCK1, serves as the primary anaerobic metabolic pathway in mollusks experiencing heat stress, with CfOpDH3 catalyzing the glucose-opine pathway in striated muscle as supplementary. Additionally, the high and stable expression level of CfASTs is crucial for the maintenance of the essential functions of aspartate aminotransferase (AST). This study provides a comprehensive and systematic analysis of the key enzymes involved in anaerobic metabolism pathways, which holds significant importance in understanding the mechanism of energy supply in mollusks.

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora.

BACKGROUND: Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation. RESULTS: In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR. CONCLUSIONS: This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

The intestinal bacterial community over seasons and its relationship with physiological status of Yesso scallop Patinopecten yessoensis.

Emerging evidence indicates that the intestinal bacterial communities associated with eukaryotes play critical roles in the physiological activities and health of their hosts. Yesso scallop Patinopecten yessoensis, one of the cold-water aquaculture species in the North Yellow Sea of China, has suffered from massive mortality in recent years. In the present study, P. yessoensis were collected from Zhangzi Island, Dalian from March 2021 to January 2022 to investigate the intestinal bacterial community and physiological indices. 16S rRNA gene sequencing data revealed that the diversity of intestinal bacteria changed significantly over seasons, with the highest Chao1 (237.42) and Shannon (6.13) indices detected in January and the lowest Chao1 (115.44) and Shannon (2.73) indices detected in July. Tenericutes, Proteobacteria and Firmicutes were dominant phyla in the intestinal bacteria of P. yessoensis, among which Firmicutes and Proteobacteria significantly enriched in August and January, respectively. Mycoplasma was the most abundant genus during the sampling period, which exhibited the highest abundance in October (75.26%) and lowest abundance in August (13.15%). The functional profiles of intestinal bacteria also exhibited seasonal variation, with the pathways related to pentose phosphate and deoxyribonucleotides biosynthesis enriched in August while the glycogen biosynthesis pathway enriched in October. Redundancy analysis showed that seawater pH, dissolved inorganic nitrogen and silicate were major environmental factors driving the temporal succession of scallop intestinal bacteria. Correlation clustering analysis suggested that the relative abundances of Endozoicomonas and Vibrio in the intestine were positively correlated with superoxide dismutase activity in hepatopancreas while negatively correlated with malondialdehyde content in hepatopancreas and glycogen content in adductor muscle. All the results revealed that the intestine harbored a lower bacterial diversity and a higher abundance of Vibrio in August, compared to January, which were closely related to the oxidative stress status of scallop in summer. These findings will advance our understanding of the relationship between seasonal alteration in the intestinal bacteria and the physiological status of scallops.

Effect of high temperature stress on glycogen metabolism in gills of Yesso scallop Patinopecten yessoensis.

Glycogen is the main energy storage material in mollusc, and the regulation of its metabolism is essential for the response against high temperature stress. In the present study, the alternation of lactic acid (LD) content, glycogen reserves, mRNA expression level of genes encoding glycogen metabolism enzymes and activities of glycogen metabolism enzymes in gills of Yesso scallop Patinopecten yessoensis after an acute high temperature treatment at 25 °C were examined to understand the effect of high temperature on glycogen metabolism. The activity of T-ATPase in gills of scallops presented a gradual increase trend especially at 6 h after the acute high temperature treatment (p < 0.05). The glycogen reserves did not change significantly even there was a downward trend at 24 h after the acute high temperature treatment (p > 0.05). The mRNA transcripts of glycogen synthase (PyGCS) in gills of scallops decreased significantly at 1, 3, 6 and 12 h (p < 0.05), and recovered to normal level at 24 h (p > 0.05) after the acute high temperature treatment, while that of glycogen phosphorylase a (PyGPa) and phosphoenol pyruvate carboxy kinase (PyPEPCK) were both significantly down-regulated from 1 h to 24 h (p < 0.05) after the acute high temperature treatment. The activity of PyGPa at 1, 12 and 24 h and the content of LD at 3 and 24 h in gills of scallops after the acute high temperature treatment both increased significantly (p < 0.05). Furthermore, the mRNA transcripts of hexokinase (PyHK) and pyruvate kinase (PyPK) in gills of scallops increased significantly (p < 0.05) after the acute high temperature treatment, and the response of PyHK was stronger. However, there was no significant difference on the activity of PyPK in gills of scallops between the experimental samples and the blank samples (p > 0.05). In addition, the mRNA transcripts of citrate synthase (PyCS) in gills of scallops were significantly down-regulated at 6 h and 12 h (p < 0.05), and finally returned to normal level at 24 h (p > 0.05) after the acute high temperature treatment. These results collectively indicated acute high temperature stress leaded the alternation of glycogen metabolism in the gills of Yesso scallop, glycogenesis, gluconeogenesis and TCA cycle were inhibited, and the glycolysis pathway of glycogen was enhanced to produce more energy for coping with environmental pressure.

Characterization of TRAF genes and their responses to Vibrio anguillarum challenge in Argopecten scallops.

The tumor necrosis factor receptor-related factor (TRAF) family has been reported to be involved in many immune pathways, such as TNFR, TLR, NLR, and RLR in animals. However, little is known about the roles of TRAF genes in the innate immune of Argopecten scallops. In this study, we first identified five TRAF genes, including TRAF2, TRAF3, TRAF4, TRAF6 and TRAF7, but not TRAF1 and TRAF5, from both the bay scallop A. irradians (Air) and the Peruvian scallop A. purpuratus (Apu). The phylogenetic analysis showed that the TRAF genes in Argopecten scallops (AiTRAF) belong to the branch of molluscan TRAF family, which lacks TRAF1 and TRAF5. Since TRAF6 is a key bridge factor in the tumor necrosis factor superfamily and plays an important role in innate and adaptive immunity, we cloned the ORFs of the TRAF6 gene in both A. irradians and A. purpuratus, as well as in two reciprocal hybrids (Aip for the hybrid Air × Apu and Api for the hybrid Apu × Air). Differences in conformational and post-translational modification resulted from the variation in amino acid sequences may cause differences in activity among them. Analysis of conserved motifs and protein structural domains revealed that AiTRAF contains typical structural domains similar to those of other mollusks and has the same conserved motifs. Tissue expression of TRAF in Argopecten scallops challenged by Vibrio anguillarum was examined by qRT-PCR. The results showed that AiTRAF were higher in gill and hepatopancreas. When challenged by Vibrio anguillarum, the expression of AiTRAF was significantly increased compared with the control group, indicating that AiTRAF may play an important role in the immunity of scallops. In addition, the expression of TRAF was higher in Api and Aip than in Air when challenged by Vibrio anguillarum, suggesting that TRAF may have contributed to the high resistance of Api and Aip to Vibrio anguillarum. The results of this study may provide new insights into the evolution and function of TRAF genes in bivalves and ultimately benefit scallop breeding.

The involvement of AMP-activated protein kinase α in regulating glycolysis in Yesso scallop Patinopecten yessoensis under high temperature stress.

AMP-activated protein kinase α subunit (AMPKα), the central regulatory molecule of energy metabolism, plays an important role in maintaining energy homeostasis and helping cells to resist the influence of various adverse factors. In the present study, an AMPKα was identified from Yesso scallop Patinopecten yessoensis (PyAMPKα). The open reading frame (ORF) of PyAMPKα was of 1599 bp encoding a putative polypeptide of 533 amino acid residues with a typical KD domain, a α-AID domain and a α-CTD domain. The deduced amino acid sequence of PyAMPKα shared 59.89-74.78% identities with AMPKαs from other species. The mRNA transcripts of PyAMPKα were found to be expressed in haemocytes and all the examined tissues, including gill, mantle, gonad, adductor muscle and hepatopancreas, with the highest expression level in adductor muscle. PyAMPKα was mainly located in cytoplasm of scallop haemocytes. At 3 h after high temperature stress treatment (25 °C), the mRNA transcripts of PyAMPKα, the phosphorylation level of PyAMPKα at Thr170 and the lactic acid (LD) content in adductor muscle all increased significantly, while the glycogen content decreased significantly. The activity of pyruvate kinase (PyPK) and the relative mRNA expression level of phosphofructokinase (PyPFK) were significantly up-regulated at 3 h after high temperature stress treatment (25 °C). Furthermore, the PyAMPKα activator AICAR could effectively upregulate the phosphorylation level of PyAMPKα, and increase activities of PyPFK and pyruvate kinase (PyPK). Meanwhile the glycogen content also declined under AICAR treatment. These results collectively suggested that PyAMPKα was involved in the high temperature stress response of scallops by enhancing glycolysis pathway of glycogen. These results would be helpful for understanding the functions of PyAMPKα in maintaining energy homeostasis under high temperature stress in scallops.

MicroRNA expression signature in the striated and smooth adductor muscles of Yesso scallop Patinopecten yessoensis.

Increasing evidences point to the potential role of microRNAs (miRNAs) in muscle growth and development in animals. However, knowledge on the identity of miRNAs and their targets in molluscs remains largely unknown. Scallops have one large adductor muscle, composed of fast (striated) and slow (smooth) muscle types, which display great differences in muscle fibers, meat quality, cell types and molecular components. In the present study, we performed a comprehensive investigation of miRNA transcriptomes in fast and slow adductor muscles of Yesso scallop Patinopecten yessoensis. As a result, 47 differentially expressed miRNAs representing ten miRNA families were identified between the striated and smooth adductor muscles. The KEGG enrichment analysis of their target genes were mainly associated with amino acid metabolism, energy metabolism and glycan biosynthesis. The target genes of miR-133 and miR-71 were validated by the dual-luciferase reporter assays and miRNA antagomir treatment in vivo. The identification and functional validation of these different miRNAs in scallops will greatly help our understanding of miRNA regulatory mechanism that achieves the unique muscle phenotypes in scallops. The present findings provide the direct evidences for muscle-specific miRNAs involved in muscle growth and differentiation in molluscs.

Additive and epistatic effects influence spectral tuning in molluscan retinochrome opsin.

The relationship between genotype and phenotype is non-trivial because of the often complex molecular pathways that make it difficult to unambiguously relate phenotypes to specific genotypes. Photopigments, comprising an opsin apoprotein bound to a light-absorbing chromophore, present an opportunity to directly relate the amino acid sequence to an absorbance peak phenotype (λmax). We examined this relationship by conducting a series of site-directed mutagenesis experiments of retinochrome, a non-visual opsin, from two closely related species: the common bay scallop, Argopecten irradians, and the king scallop, Pecten maximus. Using protein folding models, we identified three amino acid sites of likely functional importance and expressed mutated retinochrome proteins in vitro. Our results show that the mutation of amino acids lining the opsin binding pocket is responsible for fine spectral tuning, or small changes in the λmax of these light-sensitive proteins. Mutations resulted in a blue or red shift as predicted, but with dissimilar magnitudes. Shifts ranged from a 16 nm blue shift to a 12 nm red shift from the wild-type λmax. These mutations do not show an additive effect, but rather suggest the presence of epistatic interactions. This work highlights the importance of binding pocket shape in the evolution of spectral tuning and builds on our ability to relate genotypic changes to phenotypes in an emerging model for opsin functional analysis.

Transcriptional expression analysis reveals multiple effects of nonylphenol exposure on scallop immune system.

Nonylphenol (NP) is an endocrine disruptor and environmental hormone representing alkylphenol compounds. Marine mollusks are an important source of protein for people worldwide. Many researchers have begun to study the effect of NP on marine mollusks immune system in view of its toxicity; however, the underlying molecular mechanisms require in-depth analysis. In this study, we focused on the transcriptional expression change of immune-related genes and antioxidant enzymes activities variation after NP exposure in a marine bivalve mollusk, Chlamys farreri, to explore the immunomodulatory capacity of NP in marine mollusks. We identified MAVS (Mitochondrial antiviral signaling protein), a key adaptor molecule in the RLR (RIG-I like receptor) pathway, and studied the expression of multiple immune-related genes in response to different concentrations of NP. The key genes involved in RLR/TLR (Toll like receptor) innate immune pathway, apoptosis, and cellular antioxidation mechanism were investigated. Changes in the enzymatic activities of scallop antioxidant enzymes after NP exposure were also examined. The results revealed that the genes expression and the antioxidant enzymes activities show significant changes, thus proving that NP stimulation affects the scallop immune system. Our research results demonstrate the immunomodulatory capacity of NP in marine bivalve mollusks and lay the foundation for further in-depth analysis of the molecular mechanism of NP toxicity.

Expression Plasticity of Peroxisomal Acyl-Coenzyme A Oxidase Genes Implies Their Involvement in Redox Regulation in Scallops Exposed to PST-Producing Alexandrium.

Filter-feeding bivalves can accumulate paralytic shellfish toxins (PST) produced by toxic microalgae, which may induce oxidative stress and lipid peroxidation. Peroxisomal acyl-coenzyme A oxidases (ACOXs) are key enzymes functioning in maintaining redox and lipid homeostasis, but their roles in PST response in bivalves are less understood. Herein, a total of six and six ACOXs were identified in the Chlamys farreri and Patinopecten yessoensis genome, respectively, and the expansion of ACOX1s was observed. Gene expression analysis revealed an organ/tissue-specific expression pattern in both scallops, with all ACOXs being predominantly expressed in the two most toxic organs, digestive glands and kidneys. The regulation patterns of scallop ACOXs after exposure to different PST-producing algaes Alexandrium catenella (ACDH) and A. minutum (AM-1) were revealed. After ACDH exposure, more differentially expressed genes (DEGs) were identified in C. farreri digestive glands (three) and kidneys (five) than that in P. yessoensis (two), but the up-regulated DEGs showed similar expression patterns in both species. In C. farreri, three DEGs were found in both digestive glands and kidneys after AM-1 exposure, with two same CfACOX1s being acutely and chronically induced, respectively. Notably, these two CfACOX1s also showed different expression patterns in kidneys between ACDH (acute response) and AM-1 (chronic response) exposure. Moreover, inductive expression of CfACOXs after AM-1 exposure was observed in gills and mantles, and all DEGs in both tissues were up-regulated and their common DEGs exhibited both acute and chronic induction. These results indicate the involvement of scallop ACOXs in PST response, and their plasticity expression patterns between scallop species, among tissues, and between the exposure of different PST analogs.

Cell type diversity in scallop adductor muscles revealed by single-cell RNA-Seq.

Studies on cell atlas in marine invertebrates provide a better understanding of cell types, stem cell maintenance, and lineages of cell differentiation. To investigate the molecular features of various cell types in molluscan muscles, we performed single-cell RNA sequencing (scRNA-seq) to map cell types in scallop adductor muscles. We uncovered the cell type-specific features of 20 cell clusters defined by the expression of multiple specific molecular markers. These cell clusters are mainly classified into four broad classes, including mesenchymal stem cells, muscle cells, neurons, and haemolymph cells. In particular, we identified a diverse repertoire of neurons in the striated adductor muscle, but not in the smooth muscle. We further reconstructed the cell differentiation events using all the cell clusters by single-cell pseudotemporal trajectories. By integrating dual BrdU-PCNA immunodetection, neuron-specific staining and electron microscopy observation, we showed the spatial distribution of mesenchymal stem cells and neurons in striated adductor muscle of scallops. The present findings will not only be useful to address the cell type-specific gene expression profiles in scallop muscles, but also provide valuable resources for cross-species comparison of marine organisms.

Transcriptomic and metabonomic analyses reveal roles of VPS 29 in carotenoid accumulation in adductor muscles of QN Orange scallops.

BACKGROUND: In our previous studies, we demonstrated that the accumulation of carotenoids in QN Orange scallops might be regulated by the vacuolar protein sorting 29 (VPS29) gene. VPS genes are involved in pigments accumulation (including carotenoids) in some species and VPS29 is known as the core component of the membrane transport complex Retromer. However, the possible mechanism of carotenoids accumulation underlying the VPS29 remains unexplored. This study aimed to further elucidate the roles of VPS29 in the carotenoid deposition. RESULTS: Transcriptomic analyses revealed four differentially expressed genes related to carotenoid accumulation, including three down-regulated genes, low-density lipoprotein receptor domain class, scavenger receptor, Niemann Pick C1-like 1, and one up-regulated gene, ATP binding cassette transporter in RNAi group. Results from metabonomic analyses indicated increased profiles of retinol and decreased fatty acids between the RNAi and the control group. CONCLUSIONS: It thus speculated that VPS may be related to the accumulation of carotenoids as RNAi of VPS 29 seemed to result in a reduction in pectenolone through the blockage in the absorption of carotenoids and an accelerated cleavage of carotenoids into retinol.

Integrated miRNA-mRNA analysis provides potential biomarkers for selective breeding in bay scallop (Argopecten irradians).

Growth traits were compared between selected Argopecten irradians (BA) and non-selected A. irradians (NA; as a control). The results indicated that 1) the BA line exhibited greater average body weight and adductor muscle wet weight increase compared with the NA line at the same age of 10 months. 2) Comparative and integrated microRNA (miRNA) and mRNA transcriptome analyses identified 3373 differentially expressed genes (DEGs), 33 differentially expressed miRNAs (DEMs), and 39 "DEM-DEG" pairs in the BA line compared with the control. DEGs, DEMs, and "DEM-DEG" pairs involved in insulin signaling, immune related pathways, and actin cytoskeleton regulation were identified as candidates correlated with growth improvement in the BA line. A total of 259 positively selected genes were also identified. Collectively, our observations in this study will enrich the molecular information for A. irradians and provide potential biomarkers for future selective breeding and new seed creation in scallops.

Transcriptome and Network Analyses Reveal the Gene Set Involved in PST Accumulation and Responses to Toxic Alexandrium minutum Exposure in the Gills of Chlamys farreri.

Bivalve molluscs are filter-feeding organisms that can accumulate paralytic shellfish toxins (PST) through ingesting toxic marine dinoflagellates. While the effects of PST accumulation upon the physiology of bivalves have been documented, the underlying molecular mechanism remains poorly understood. In this study, transcriptomic analysis was performed in the gills of Zhikong scallop (Chlamys farreri) after 1, 3, 5, 10, and 15 day(s) exposure of PST-producing dinoflagellate Alexandrium minutum. Higher numbers of differentially expressed genes (DEGs) were detected at day 1 (1538) and day 15 (989) than that at day 3 (77), day 5 (82), and day 10 (80) after exposure, and most of the DEGs were only regulated at day 1 or day 15, highlighting different response mechanisms of scallop to PST-producing dinoflagellate at different stages of exposure. Functional enrichment results suggested that PST exposure induced the alterations of nervous system development processes and the activation of xenobiotic metabolism and substance transport processes at the acute and chronic stages of exposure, respectively, while the immune functions were inhibited by PST and might ultimately cause the activation of apoptosis. Furthermore, a weighted gene co-expression network was constructed, and ten responsive modules for toxic algae exposure were identified, among which the yellow module was found to be significantly correlated with PST content. Most of the hub genes in the yellow module were annotated as solute carriers (SLCs) with eight being OCTN1s, implying their dominant roles in regulating PST accumulation in scallop gills. Overall, our results reveal the gene set responding to and involved in PST accumulation in scallop gills, which will deepen our understanding of the molecular mechanism of bivalve resistance to PST.

Genome-wide DNA methylation profile changes associated with shell colouration in the Yesso scallop (Patinopecten yessoensis) as measured by whole-genome bisulfite sequencing.

BACKGROUND: Mollusca, a phylum of highly rich species, possess vivid shell colours, but the underlying molecular mechanism remains to be elucidated. DNA methylation, one of the most common epigenetic modifications in eukaryotes, is believed to play a vital role in various biological processes. However, analysis of the effects of DNA methylation on shell colouration has rarely been performed in molluscs, limiting the current knowledge of the molecular mechanism of shell colour formation. RESULTS: In the present study, to reveal the role of epigenetic regulation in shell colouration, WGBS, the "gold standard" of DNA methylation analysis, was first performed on the mantle tissues of Yesso scallops (Patinopecten yessoensis) with different shell colours (brown and white), and DNA methylomes at single-base resolution were generated. About 3% of cytosines were methylated in the genome of the Yesso scallop. A slight increase in mCG percentage and methylation level was found in brown scallops. Sequence preference of nearby methylated cytosines differed between high and low methylation level sites and between the brown- and white-shelled scallops. DNA methylation levels varied among the different genomic regions; all the detected regions in the brown group exhibited higher methylation levels than the white group. A total of 41,175 DMRs (differentially methylated regions) were detected between brown and white scallops. GO functions and pathways associated with the biosynthesis of melanin and porphyrins were significantly enriched for DMRs, among which several key shell colour-related genes were identified. Further, different correlations between mRNA expression levels and DNA methylation status were found in these genes, suggesting that DNA methylation regulates shell colouration in the Yesso scallop. CONCLUSIONS: This study provides genome-wide DNA methylation landscapes of Yesso scallops with different shell colours, offering new insights into the epigenetic regulatory mechanism underlying shell colour.