Blood clams differ from their molluscan kins by exhibiting a unique red-blood (RB) phenotype; however, the genetic basis and biochemical machinery subserving this evolutionary innovation remain unclear. As a fundamental step toward resolving this mystery, we presented the first chromosome-level genome and comprehensive transcriptomes of the blood clam Tegillarca granosa for an integrated genomic, evolutionary, and functional analyses of clam RB phenotype. We identified blood clam-specific and expanded gene families, as well as gene pathways that are of RB relevant. Clam-specific RB-related hemoglobins (Hbs) showed close phylogenetic relationships with myoglobins (Mbs) of blood clam and other molluscs without the RB phenotype, indicating that clam-specific Hbs were likely evolutionarily derived from the Mb lineage. Strikingly, similar to vertebrate Hbs, blood clam Hbs were present in a form of gene cluster. Despite the convergent evolution of Hb clusters in blood clam and vertebrates, their Hb clusters may have originated from a single ancestral Mb-like gene as evidenced by gene phylogeny and synteny analysis. A full suite of enzyme-encoding genes for heme synthesis was identified in blood clam, with prominent expression in hemolymph and resembling those in vertebrates, suggesting a convergence of both RB-related Hb and heme functions in vertebrates and blood clam. RNA interference experiments confirmed the functional roles of Hbs and key enzyme of heme synthesis in the maintenance of clam RB phenotype. The high-quality genome assembly and comprehensive transcriptomes presented herein serve new genomic resources for the super-diverse phylum Mollusca, and provide deep insights into the origin and evolution of invertebrate RB.

The Chromosome-Level Genome Assembly and Comprehensive Transcriptomes of the Razor Clam (Sinonovacula constricta)

MicroRNAs (miRNAs) are endogenous non-coding small RNAs (sRNAs) that can base pair with their target mRNAs, which represses their translation or induces their degradation in various biological processes. To identify miRNAs regulated by heavy metal stress, we constructed two sRNA libraries for the blood clam Tegillarca granosa: one for organisms exposed to toxic levels of cadmium (Cd) and one for a control group. Sequencing of the two libraries and subsequent analysis revealed 215 conserved and 39 new miRNAs. Most of the new miRNAs in T. granosa were up- or down-regulated in response to Cd exposure. There were significant differences in expression between the Cd and control groups for 16 miRNAs. Of these, five miRNAs were significantly up-regulated and 11 were significantly down-regulated in the Cd stress library. Potential targets were predicted for the 16 differential miRNAs in pre-miRNAs identified according to sequence homology. Some of the predicted miRNA targets are associated with regulation of the response to stress induced by heavy metals. Five differentially expressed miRNAs (Tgr-nmiR-8, Tgr-nmiR-21, Tgr-miR-2a, Tgr-miR-10a-5p, and Tgr-miR-184b) were validated by qRT-PCR. Our study is the first large-scale identification of miRNAs in T. granosa haemocytes. Our findings suggest that some miRNAs and their target genes and pathways may play critical roles in the responses of this species to environmental heavy metal stresses.

Oysters are an important aquaculture species distributed worldwide, including in Zhejiang Province, located on the east coast of China. Because of the high diversity and complicated introduction history of oysters and their seedlings, there has been much disagreement regarding the origin of each species, and the dominant and indigenous species remain unclear. We sampled 16 batches of oysters from seven sites in three aquaculture bays and found two main oyster species, Crassostrea sikamea and Crassostrea angulata. The former occupied the higher intertidal zone and comprised more than 70% of the cultured oysters. Based on the cytochrome oxidase C subunit I (COI) and mitochondrial noncoding region (MNR), C. sikamea showed higher genetic diversity than C. angulata. The analysis of molecular variance among COI sequences of these species from the Xiangshan Bay populations were comparable to those of other populations and showed that most of the molecular variance was within groups, which was consistent with the low pairwise fixation index F_ST values. The neutrality test revealed that C. sikamea experienced population expansion events, whereas for C. angulata, the significant Fu’s Fs and non-significant Tajima’s D test results may indicate a possible population expansion event, implying that C. sikamea is likely an indigenous species. The method established based on internal transcribed spacer 1 digestion by the HindIII restriction enzyme is useful for identifying C. sikamea and C. angulata in the local region. The specific primers on the MNR sequence show potential for distinguishing C. sikamea from four other important Crassostrea oysters. These results highlight the abundance of C. sikamea on the Zhejiang coast and lay the foundation for protecting and utilizing the local oyster germplasm resources and for the sustainable development of the oyster industry.

Accumulation of excessive ammonia is a big threat to aquatic animals, which causes adverse effects on the health, production reduction, and even high mortality. The razor clam Sinonovacula constricta, a bivalve living in intertidal mudflat with a deep-burrowing lifestyle, often faces a high concentration of ambient ammonia. However, there is less available information concerning the toxic effects of ammonia on razor clam and its molecular mechanisms of adaptation to ammonia stress. The aim of this study was to investigate the effects of ammonia exposure on the gill and hepatopancreas of razor clam by transcriptome sequencing. The results showed that the median lethal concentration of ammonia was 244.55 mg/L for 96 h. A total of 1,415 and 306 differentially expressed genes (DEGs) were identified in the gill and hepatopancreas, respectively. The functional annotation showed that DEGs of the gill were mainly involved in the regulation of nitrogen compound metabolic process, nitrogen compound transport, and amide transport. The DEGs of the hepatopancreas were mostly enriched in oxidation-reduction process, response to stress, and amine metabolic process. The expression levels of NH₃/NH₄⁺ transporting channels and H⁺ excreting-related genes, including Rhesus glycoproteins (Rh), Na⁺/K⁺-ATPase (NKA), Na⁺/H⁺ exchanger, V-ATPase (VHA), and carbonic anhydrase (CA), were upregulated significantly in the gill (p < 0.05). In addition, the expression levels of glutamine and urea synthesis-related genes that played vital roles in ammonia detoxification, such as glutamine synthetase (GS), arginase (ARG), and argininosuccinate synthetase (ASS), were also increased obviously in the hepatopancreas (p < 0.05). Taken together, our results indicate that the synergistic action of ammonia excretion in the gill and ammonia metabolism in the hepatopancreas might be the mechanism through which the clams tolerate to environmental ammonia. This study provides a molecular basis for the better evaluation of the responding mechanism of ammonia tolerance.