Background: Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa. Results: Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors. Conclusions: Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.

The pen shell, Atrina pectinata, is one of the commercial bivalves in East Asia and thought to be recently affected by anthropogenic pressure (habitat destruction and/or fishing pressure). Information on its population genetic structure is crucial for the conservation of A. pectinata. Considering its long pelagic larval duration and iteroparity with high fecundity, the genetic structure for A. pectinata could be expected to be weak at a fine scale. However, the unusual oceanography in the coasts of China and Korea suggests potential for restricted dispersal of pelagic larvae and geographical differentiation. In addition, environmental changes associated with Pleistocene sea level fluctuations on the East China Sea continental shelf may also have strongly influenced historical population demography and genetic diversity of marine organisms. Here, partial sequences of the mitochondrial Cytochrome c oxidase subunit I (COI) gene and seven microsatellite loci were used to estimate population genetic structure and demographic history of seven samples from Northern China coast and one sample from North Korea coast. Despite high levels of genetic diversity within samples, there was no genetic differentiation among samples from Northern China coast and low but significant genetic differentiation between some of the Chinese samples and the North Korean sample. A late Pleistocene population expansion, probably after the Last Glacial Maximum, was also demonstrated for A. pectinata samples. No recent genetic bottleneck was detected in any of the eight samples. We concluded that both historical recolonization (through population range expansion and demographic expansion in the late Pleistocene) and current gene flow (through larval dispersal) were responsible for the weak level of genetic structure detected in A. pectinata.

• Metabolomics reveals the metabolic responses of hard clam to hypo-salinity stress. • cAMP-PKA pathway and NKA play osmoregulatory roles in hard clams. • Activated antioxidant responses and reorganized membrane lipids occurred at 5 d. • Alaine and lactate accumulation suggest the onset of anaerobic metabolism at 1 d. • Fatty acids β-oxidation is promoted to provide energy for osmoregulation.

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group’s phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella , Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

Heat shock protein 70 (HSP70) members participate in a wide range of housekeeping and stress-related activities in eukaryotic cells. In marine ecosystems, bivalves encounter abiotic stresses, including high temperatures and low dissolved oxygen. Here, 133 MmHSP70 genes were identified through combined methods including Blastp, HMM and manual filtration, based on the whole Mercenaria mercenaria genome. The MmHSP70 genes were unevenly distributed, and 41 genes (33.08%) were located on Chr 7. Phylogenetic analyses indicated that the MmHSP70 gene family mainly consisted of two clusters and the Hspa12 subfamily underwent lineage-specific expansion. A high-density collinear gene block was observed between M. mercenaria Chr 7 and Cyclina sinensis Chr 14. Tandem duplication MmHSP70 gene pairs experienced different levels of purifying selection, which could be an important source of sequence and functional constraints. MmHSP70 genes showed tissue-specific and stress-specific expression. Most tandem duplication HSP70 gene pairs had high expression under hypoxia stress. HSP70 B2 tandem duplication gene pairs showed significantly increased expression under heat plus severe hypoxia stress. This study provided a comprehensive understanding of the MmHSP70 gene family in the M. mercenaria and laid a significant foundation for further studies on the functional characteristics of MmHSP70 genes during exposure to heat and hypoxia stress.