Developmental Transcriptome Analysis and Identification of Genes Involved in Larval Metamorphosis of the Razor Clam, Sinonovacula constricta

Transcriptome Profiling Based at Different Time Points after Hatching Deepened Our Understanding on Larval Growth and Development of Amphioctopus fangsiao

As the quality of life improves, there is an increasing demand for nutrition-rich marine organisms like fish, shellfish, and cephalopods. To address this, artificial cultivation of these organisms is being explored along with ongoing research on their growth and development. A case in point is Amphioctopus fangsiao, a highly valued cephalopod known for its tasty meat, nutrient richness, and rapid growth rate. Despite its significance, there is a dearth of studies on the A. fangsiao growth mechanism, particularly of its larvae. In this study, we collected A. fangsiao larvae at 0, 4, 12, and 24 h post-hatching and conducted transcriptome profiling. Our analysis identified 4467, 5099, and 4181 differentially expressed genes (DEGs) at respective intervals, compared to the 0 h sample. We further analyzed the expression trends of these DEGs, noting a predominant trend of continuous upregulation. Functional exploration of this trend entailed GO and KEGG functional enrichment along with protein-protein interaction network analyses. We identified GLDC, DUSP14, DPF2, GNAI1, and ZNF271 as core genes, based on their high upregulation rate, implicated in larval growth and development. Similarly, CLTC, MEF2A, PPP1CB, PPP1R12A, and TJP1, marked by high protein interaction numbers, were identified as hub genes and the gene expression levels identified via RNA-seq analysis were validated through qRT-PCR. By analyzing the functions of key and core genes, we found that the ability of A. fangsiao larvae to metabolize carbohydrates, lipids, and other energy substances during early growth may significantly improve with the growth of the larvae. At the same time, muscle related cells in A. fangsiao larvae may develop rapidly, promoting the growth and development of larvae. Our findings provide preliminary insights into the growth and developmental mechanism of A. fangsiao, setting the stage for more comprehensive understanding and broader research into cephalopod growth and development mechanisms.

Stage-Specific Transcriptomes of the Mussel Mytilus coruscus Reveals the Developmental Program for the Planktonic to Benthic Transition

Many marine invertebrate larvae undergo complex morphological and physiological changes during the planktonic-benthic transition (a.k.a. metamorphosis). In this study, transcriptome analysis of different developmental stages was used to uncover the molecular mechanisms underpinning larval settlement and metamorphosis of the mussel, Mytilus coruscus. Analysis of highly upregulated differentially expressed genes (DEGs) at the pediveliger stage revealed enrichment of immune-related genes. The results may indicate that larvae co-opt molecules of the immune system to sense and respond to external chemical cues and neuroendocrine signaling pathways forecast and trigger the response. The upregulation of adhesive protein genes linked to byssal thread secretion indicates the anchoring capacity required for larval settlement arises prior to metamorphosis. The results of gene expression support a role for the immune and neuroendocrine systems in mussel metamorphosis and provide the basis for future studies to disentangle gene networks and the biology of this important lifecycle transformation.

Expression of p38MAPK and its regulation of apoptosis under high temperature stress in the razor clam Sinonovacula constricta

p38MAPK is a key branch of the MAPK (mitogen-activated protein kinase) pathway that plays an important role in physiological processes such as apoptosis, cell proliferation and growth. In this experiment, we screened and identified one p38MAPK gene in the razor clam Sinonovacula constricta, which encoded 359 amino acids and was widely expressed in various adult tissues. After 24 h of high temperature stress at 34 °C, the transcript expression of p38MAPK showed significant changes in all tested tissues. In particular in the gill and hepatopancreas tissues, where the expression increased 1.81 and 7.83 times compared with the control group, respectively (P < 0.01). Furthermore, we examined the expression of the apoptosis suppressor gene Bcl-2 and pro-apoptosis gene Bax by knock-down of p38MAPK with dsRNA interference in the gill and hepatopancreas tissues. The obvious up-regulation expression of Bcl-2 and significant suppression of Bax were observed, respectively (P < 0.01). Moreover, the TUNEL staining technique was used to detect apoptosis before and after interference. The degree of apoptosis in the gill and hepatopancreas tissues was reduced after interference with p38MAPK, and the ROS content was significantly reduced (P < 0.01). The results suggested that p38MAPK had a regulatory role in the heat tolerance of razor clams.

Temporally balanced selection during development of larval Pacific oysters (Crassostrea gigas) inherently preserves genetic diversity within offspring

Balancing selection is one of the mechanisms which has been proposed to explain the maintenance of genetic diversity in species across generations. For species with large populations and complex life histories, however, heterogeneous selection pressures may create a scenario in which the net effects of selection are balanced across developmental stages. With replicated cultures and a pooled sequencing approach, we show that genotype-dependent mortality in larvae of the Pacific oyster (Crassostrea gigas) is largely temporally dynamic and inconsistently in favour of a single genotype or allelic variant at each locus. Overall, the patterns of genetic change we observe to be taking place are more complex than what would be expected under classical examples of additive or dominant genetic interactions. They are also not easily explained by our current understanding of the effects of genetic load. Collectively, temporally heterogeneous selection pressures across different larval developmental stages may act to maintain genetic diversity, while also inherently sheltering genetic load within oyster populations.

Transcriptome Analysis of Hormone-and Cuticle-Related Genes in the Development Process of Deutonymph in Tetranychus urticae

Tetranychus urticae is an important agricultural pest that feeds on more than 1100 plant species. To investigate gene expression network in development process of deutonymph, a comprehensive transcriptome analysis of different developmental time points of deutonymph in T. urticae was performed. Comparing with expression profile of 7 h, 309, 876, 2736, and 3432 differential expression genes were detected at time points 14 h, 21 h, 28 h, and 35 h, respectively. The expression dynamic analysis indicated that genes in hormone- (ecdysteroid and juvenile hormone) and cuticle- (chitin and cuticle proteins) related pathways were indispensable for development process in deutonymph. Among hormone related pathway genes, the ecdysteroid biosynthesis pathway genes were highly expressed at the growth period of development process, which is opposite to the expression patterns of juvenile hormone biosynthesis pathway genes. For cuticle related pathway genes, 13 chitinase genes were identified in the genome of T. urticae, and 8 chitinase genes were highly expressed in different time points of developmental process in the deutonymph of T. urticae. Additionally, 59 cuticle protein genes were identified from genome, and most of the cuticle protein genes were expressed in the molting period of developmental process in deutonymph. This study reveals critical genes involved in the development process of deutonymph and also provides comprehensive development transcriptome information for finding more molecular targets to control this pest.

Constructing a de novo transcriptome and a reference proteome for the bivalve Scrobicularia plana: Comparative analysis of different assembly strategies and proteomic analysis

Scrobicularia plana is a coastal and estuarine bivalve widely used in ecotoxicological studies. However, the underlying molecular mechanisms for S. plana pollutant responses are hardly known due to the lack of molecular databases. Thus, in this study we present a holistic approach to assess a robust reference transcriptome and proteome of this clam. A mixture of control and metal-exposed individuals was used for mRNA isolation. Four sets of high quality filtered preprocessed reads were generated (two quality scores and two sequenced lengths) and assembled with Mira, Ray and Trinity algorithms. The sixty-four generated assemblies were refined, filtered and evaluated for their proteomic quality. Eight assemblies presented top Detonate scores but one was selected due to its compactness and biological representation, which was generated: (i) from the highest quality dataset (Q20L100), (ii) using Trinity algorithm with all k-mers (AtKa), (iii) removing redundancy by CD-HIT (RR80), and (iv) filtering out poor contigs (F), that was subsequently named Q20L100AtKaRR80F. S. plana proteomic analysis revealed 10,017 peptide groups that corresponded to 2066 proteins with a wide coverage of molecular functions and biological processes, confirming the strength of the database generated.

Integrated mRNA and Small RNA Sequencing Reveals Regulatory Expression of Larval Metamorphosis of the Razor Clam

The razor clam, Sinonovacula constricta, is an important economic marine shellfish, and its larval development involves obvious morphological and physiological changes. MicroRNA plays a key role in the physiological changes of the organism through regulating targeted mRNA. This study performed miRNA-mRNA sequencing for eight different developmental stages of S. constricta using Illumina sequencing. A total of 2156 miRNAs were obtained, including 2069 known miRNAs and 87 novel miRNAs. In addition, target genes were predicted for key miRNAs differentially expressed between adjacent development samples by integrating the mRNA transcriptome. Further analysis revealed that the differentially expressed genes were enriched in complement activation, alternative pathways, translation, and negative regulation of monocyte molecular protein-1 production. KEGG pathway annotation showed significant enrichment in the regulation of the ribosome, phagosome, tuberculosis and fluid shear stress, and atherosclerosis. Ten mRNAs and ten miRNAs that are related to larval metamorphosis were identified using real-time PCR. Furthermore, the double luciferase experiment validated the negative regulatory relationship between miR-133 and peroxisome proliferator-activated receptor-γ (PPAR-γ). These results indicated that the target genes regulated by these differentially expressed miRNAs may play an important regulatory role in the metamorphosis development of S. constricta.

Chromosome-level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818)

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high-quality chromosome-level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single-molecule real-time sequencing, Illumina paired-end sequencing, 10X Genomics linked-reads and Hi-C reads. The genome size was 1,220.85 Mb, containing scaffold N50 of 65.93 Mb and contig N50 of 976.94 Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single-Copy Orthologs were determined in this genome assembly. And Hi-C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein-coding genes were predicted in the S. constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S. constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S. constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.

Identification of a dopamine receptor in Sinonovacula constricta and its antioxidant responses

The dopamine (DA) D2 receptor is a member of the G protein-coupled receptors of organisms and plays an important role in immune system regulation. The presence of DA receptors has been widely reported in vertebrates, but few studies have been conducted in shellfish. Here, we identified a novel DA-D2 receptor gene, ScDopR2-1, in the razor clam Sinonovacula constricta. ScDopR2-1 belongs to the family of G protein-coupled receptors, containing seven hydrophobic transmembrane domains, along with 16 predicted N-glycosylation sites and 69 phosphorylation sites. A longer third intracellular loop and a shorter C-terminus in ScDopR2-1 are characteristic features of D2 receptors. ScDopR2-1 is widely expressed in tissues from adult clams, showing high expression in siphon and foot tissues. Furthermore, in response to Vibrio anguillarum challenge, ScDopR2-1 expression levels are significantly increased in liver tissue. Moreover, changes in the activities of catalase (CAT) and superoxide dismutase (SOD) also indicate that the organism causes an immune response. In summary, the results indicate that ScDopR2-1 plays a pivotal role in antioxidant responses in S. constricta.

Dynamic response in the larval geoduck (Panopea generosa) proteome to elevated pCO2

Pacific geoducks (Panopea generosa) are clams found along the northeast Pacific coast where they are important components of coastal and estuarine ecosystems and a major aquaculture product. The Pacific coastline, however, is also experiencing rapidly changing ocean habitat, including significant reductions in pH. To better understand the physiological impact of ocean acidification on geoduck clams, we characterized for the first time the proteomic profile of this bivalve during larval development and compared it to that of larvae exposed to low pH conditions. Geoduck larvae were reared at pH 7.5 (ambient) or pH 7.1 in a commercial shellfish hatchery from day 6 to day 19 postfertilization and sampled at six time points for an in-depth proteomics analysis using high-resolution data-dependent analysis. Larvae reared at low pH were smaller than those reared at ambient pH, especially in the prodissoconch II phase of development, and displayed a delay in their competency for settlement. Proteomic profiles revealed that metabolic, cell cycle, and protein turnover pathways differed between the two pH and suggested that differing phenotypic outcomes between pH 7.5 and 7.1 are likely due to environmental disruptions to the timing of physiological events. In summary, ocean acidification results in elevated energetic demand on geoduck larvae, resulting in delayed development and disruptions to normal molecular developmental pathways, such as carbohydrate metabolism, cell growth, and protein synthesis.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

Inhibition of growth in juvenile manila clam Ruditapes philippinarum: Potential adverse outcome pathway of TBBPA

Tetrabromobisphenol A (TBBPA) is ubiquitous and its contents showing an increasing trend in the coastal environment. In order to investigate the effects of TBBPA on marine bivalves, juvenile manila clams Ruditapes phillipinarum were exposed to TBBPA for 28 days. The results showed that shell growth rate of juvenile clams after exposure to 62.5-1000 μg L^-1 TBBPA for 28 d were significantly inhibited (p < 0.05). Then in order to link the changes in filtration rate, mRNA expression of insulin-like growth factor homologue (IGF) and tissue thyroid hormone (TH) contents to growth, juvenile clams were exposed to 62.5 and 500 μg L^-1 TBBPA for 14 days. The transcriptional levels of neuroendocrine signals (NPF and insulin homologue) associated with filter feeding regulation, and genes of TH synthesis-related enzymes were also examined. The results showed that filtration rates was significantly reduced to 44.1% and 14% of controls after 14 d of exposure. In parallel, exposure to TBBPA significantly increased the expression levels of insulin which may elicit the filter feeding inhibition. TBBPA exposure caused alterations in tissue content of THs and mRNA expression of TH synthesis-related enzymes. However, the data showed increased T₃ content, T₃/T₄ ratio and mRNA expression of IGF. These data demonstrated that the most important key event of TBBPA could be linked to growth impairment in juveniles was the reduction of filtration rate. These results provide a robust framework towards revealing the underlying mechanism of the growth inhibition caused by TBBPA on bivalves and understanding the adverse outcome pathway across taxonomic phyla.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

Effects of Alkalinity and pH on Survival, Growth, and Enzyme Activities in Juveniles of the Razor Clam, Sinonovacula constricta

In order to clarify the possibility of rearing razor clams (Sinonovacula constricta) in inland saline water (ISW) and to facilitate their breeding under these stressful conditions, we performed semi-static acute and chronic toxicity tests to determine the effects of carbonate alkalinity (CA) and pH on the survival and growth rate, and critical metabolic enzyme activity in juvenile of S. constricta (JSC). (1) Acute toxicity test. As the water CA increased from 1.22 to 45.00 mmol L-1, the survival rate decreased significantly, which was exacerbated by the increase in the pH. When the water CA was set at 2.5 mmol L-1, the 48 h lethal concentration 50% (LC50) for JSCs with respect to pH was 9.86. When the water pH was 9.0, 9.5, and 10.0, the 48 h LC50 values for JSCs with respect to CA were 10.38, 8.79, and 3.11 mmol L-1, respectively. (2) Chronic toxicity test. Four experimental groups comprising the control, CAS, pHS, and CA-pHS were designated according to the target ISW data. After 3 months of stress, the JSC survival rate in each group exceeded 85%, but survival was significantly lower in the CA-pHS group than the control group (p < 0.05) in the first month. For the JSCs in various groups, the shell length growth rate (SGR) and weight gain (WG) rate were significantly lower in the CA-pHS group than the other groups (p < 0.05 for SGR; p < 0.001 for WG) in the first month. However, the difference in the growth rate among groups decreased in the next 2 months. For the JSCs in the CA-pHS group, the oxygen consumption, ammonia-N excretion, Na+/K+-ATPase, aspartate aminotransferase, and superoxide dismutase levels were significantly higher than those in the other groups during the first month, but there were no significant differences between the groups subsequently. The acetylcholinesterase and lysozyme levels did not differ significantly among groups during stress for 3 months. The integrated biomarker response index showed that stressors comprising high pH and CA could be tolerated well by JSCs over long periods of stress. These results indicate that water CA and pH together affect the survival, growth, and physiological activity of JSCs. S. constricta is suitable for culture in ISW.

Comparative Transcriptome Analysis Reveals Growth-Related Genes in Juvenile Chinese Sea Cucumber, Russian Sea Cucumber, and Their Hybrids

Heterosis is important for sea cucumber breeding, but its molecular mechanism remains largely unexplored. In this study, parental lines of Apostichopus japonicus from Russia (R) and China (C) were used to construct hybrids (CR and RC) by reciprocal crossing. We examined the transcriptional profiles of the hybrids (CR and RC) and the purebreds (CC and RR) at different developmental times. A total of 60.27 Gb of clean data was obtained, and 176,649 unigenes were identified, of which 50,312 unigenes were annotated. A total of 414,536 SNPs were identified. A total of 7011 differentially expressed genes (DEGs) were obtained between the purebreds and hybrids at 45 days after fertilization (DAF), and a total of 8218 DEGs were obtained between the purebreds and hybrids at 75 DAF. In addition, a total of 7652 DEGs were obtained between 45 DAF and 75 DAF. The significant DEGs were mainly involved in the MAPK and FOXO signaling pathways, especially in the Ras-Raf-MEK1/2-ERK module, which may be a key regulator of development and growth in juvenile A. japonicus. In addition, we also identified key growth-related genes, such as fgfs, igfs, megfs and hgfs, which were upregulated in the hybrids (RC and CR); these genes may play important roles in heterosis in A. japonicus. Our study provides fundamental information on the molecular mechanisms underlying heterosis in sea cucumber and might suggest strategies for the selection of rapidly growing strains of sea cucumber in aquaculture.

Two I84 family protease inhibitors from Chinese razor clams Sinonovacula constricta expressed in response to environmental challenges

Protease inhibitors play critical roles in numerous biological processes including host defense in all multicellular organisms. Eighty three evolutionary families of protease inhibitors are currently accommodated in the MEROPS database and the I84 family currently consists of 3 novel serine protease inhibitors from the eastern oyster Crassostrea virginica. In this study, we identified 2 new I84 family members from the Chinese razor clam Sinonovacula constricta, scSI-1 and scSI-2, using cDNA cloning and sequencing. The scSI-1 cDNA consisted of 494 bp with a 282 bp ORF encoding a 93-amino acid polypeptide that was predicted to have a 19-amino acid signal peptide and a 74-residue mature protein with a calculated molecular mass of 8248.5 Da. The scSI-2 cDNA was 490 bp long with a 273 bp ORF encoding a 90-amino acid polypeptide that was predicted to have an 18-amino acid signal peptide and a 72-residue nature protein with a calculated molecular mass of 7528.4 Da. ScSI-1 and scSI-2 shared high sequence similarity with the 3 known members of I84 family and both expressed primarily in the clam digestive glands. Protease inhibitory activity in the clam plasma also exhibited the signature kinetic characteristics of the I84 members from the oyster. In addition, levels of scSI-1 and scSI-2 gene expression in digestive glands and the protease inhibitory activity in plasma elevated significantly in clams challenged by bacterial injections and Vibrio harveyi was more effective than Staphylococcus epidermidis in inducing the gene expression and plasma protease inhibitory activity. Moreover, drastic changes of salinity and temperature also caused significant changes in the gene expression and plasma activity. These results indicated that scSI-1 and scSI-2 represented 2 new members of the I84 family and they likely play a role in clam host defense against infections and in reactions against physiochemical stressors.

Impact of High Hydrostatic Pressure on the Shelling Efficacy, Physicochemical Properties, and Microstructure of Fresh Razor Clam (Sinonovacula constricta)

The effects of high hydrostatic pressure (HHP) treatments (200, 300, and 400 MPa for 1, 3, 5 and 10 min) on the shelling efficacy (the rate of shelling, the rate of integrity and yield of razor clam meat) and the physicochemical (drip loss, water-holding capacity, pH, conductivity, lipid oxidation, Ca²⁺ -ATPase activity, myofibrillar protein content), microbiological (total viable counts) and microstructural properties of fresh razor clam (Sinonovacula constricta) were investigated. HHP treatments significantly (P < 0.05) increased shelling efficiency, water-holding capacity, pH, conductivity, and lipid oxidation, and HHP-treated razor clam showed lower levels of microorganisms and drip loss than untreated razor clam. Levels of thiobarbituric acid reacting substances (TBA) in HHP-treated razor clam were greatly increased (up to 0.93 ± 0.09 mg MDA/kg at 400 MPa for 10 min) which was caused by the formation of hydroperoxides during HHP treatment. All HHP treatments were found to have adverse effects on the activity of Ca²⁺ -ATPase and the content of myofibrillar protein (MP), which might be due to the substantial damage to the tertiary structure of proteins at high pressure. Moreover, scanning electron microscopy (SEM) revealed the compaction of the muscle fibers and a decrease in the extracellular space with increasing pressure and holding time. This phenomenon was mainly correlated with the compaction of muscle fibers and denaturation, aggregation, and gelation of muscle protein triggered by high pressure. In general, HHP could be applied as a safe and effective nonthermal technology to produce high-quality shelled razor clam.

PRACTICAL APPLICATION

High hydrostatic pressure (HHP) is now well known as a nonthermal processing technology and becoming increasingly acknowledged. However, it has not been widely applied to shell seafood due to its uncertain influence on its quality and shelling property. This study could provide valuable information regarding the shelling efficacy, physicochemical properties, and microstructure of razor clam treated by HHP. And it demonstrated that HHP showed a positive impact on quality of razor clam treated by HHP.

Construction of the First High-Density Genetic Linkage Map and Analysis of Quantitative Trait Loci for Growth-Related Traits in Sinonovacula constricta

The razor clam (Sinonovacula constricta) is an important aquaculture species, for which a high-density genetic linkage map would play an important role in marker-assisted selection (MAS). In this study, we constructed a high-density genetic map and detected quantitative trait loci (QTLs) for Sinonovacula constricta with an F₁ cross population by using the specific locus amplified fragment sequencing (SLAF-seq) method. A total of 315,553 SLAF markers out of 467.71 Mreads were developed. The final linkage map was composed of 7516 SLAFs (156.60-fold in the parents and 20.80-fold in each F₁ population on average). The total distance of the linkage map was 2383.85 cM, covering 19 linkage groups with an average inter-marker distance of 0.32 cM. The proportion of gaps less than 5.0 cM was on average 96.90%. A total of 16 suggestive QTLs for five growth-related traits (five QTLs for shell height, six QTLs for shell length, three QTLs for shell width, one QTL for total body weight, and one QTL for soft body weight) were identified. These QTLs were distributed on five linkage groups, and the regions showed overlapping on LG9 and LG13. In conclusion, the high-density genetic map and QTLs for S. constricta provide a valuable genetic resource and a basis for MAS.

RNA-seq analysis revealed ROS-mediated related genes involved in cadmium detoxification in the razor clam Sinonovacula constricta

The razor clam Sinonovacula constricta is a eukaryotic benthic intertidal bivalve species that is tolerant to different heavy metals, such as cadmium ion (Cd(2+)). However, the mechanism by which S. constricta responds to Cd(2+)-induced stress remains unclear. In this study, eight transcriptome libraries were constructed and sequenced using razor clams exposed to Cd(2+) for 12 and 48 h. A total of 18,330 unigenes with an average length of 500 bp were annotated. Among these 18,330 unigenes, 582 and 649 displayed differential expression profiles at 12 and 48 h in the gill, respectively. The corresponding differential unigenes in the hepatopancreas were 1056 and 382. Gene Ontology annotation revealed that these unigenes were highly pronounced in metabolic process, cellular process, binding, and catalytic activity. Notably, ROS production-related genes, such as heat shock proteins 32, metallothionein, and glutathione, were synchronously enriched in all experimental samples with induced expression profiles, which was also validated by qPCR. Our results highlighted the relation between immune regulation and Cd(2+)-induced stress in razor clam and provided new insights into the molecular mechanisms of heavy toxicology.