The hard-shell mussel (Mytilus coruscus) is widespread in the temperate coastal areas of the northwest Pacific and holds a significant position in the shellfish aquaculture market in China. However, the natural resources of this species have been declining, and population genetic studies of M. coruscus are also lacking. In this study, we conducted whole-genome resequencing (WGR) of M. coruscus from eight different latitudes along the Chinese coast and identified a total of 25,859,986 single nucleotide polymorphism (SNP) markers. Our findings indicated that the genetic diversity of M. coruscus from the Zhoushan region was lower compared with populations from other regions. Furthermore, we observed that the evolutionary tree clustered into two primary branches, and the Zhangzhou (ZZ) population was in a separate branch. The ZZ population was partly isolated from populations in other regions, but the distribution of branches was not geographically homogeneous, and a nested pattern emerged, consistent with the population differentiation index (F_ST) results. To investigate the selection characteristics, we utilized the northern M. coruscus populations (Dalian and Qingdao) and the central populations (Zhoushan and Xiangshan) as reference populations and the southern ZZ population as the target population. Our selection scan analysis identified several genes associated with thermal responses, including Hsp70 and CYP450. These genes may play important roles in the adaptation of M. coruscus to different living environments. Overall, our study provides a comprehensive understanding of the genomic diversity of coastal M. coruscus in China and is a valuable resource for future studies on genetic breeding and the evolutionary adaptation of this species.

The homologous E6-AP carboxy-terminal structural domain (HECT) contained in E3 ubiquitin ligases (E3s) is a key factor in protein degradation and maintenance of cellular homeostasis in animals. However, the functional roles and evolutionary aspects of the HECT gene family in bivalve mussels remain unclear and warrant further investigation. In this study, we identified 22 HECT genes within the genome of Mytilus coruscus Gould, all containing a conserved HECT structural domain derived from dispersed repeats, distributed unevenly across 11 chromosomes. Phylogenetic analysis classified M. coruscus HECT genes into six major classes, with amino acid sequences within the same evolutionary clade displaying similar conserved motifs. Homology analysis with HECT genes of four bivalve species revealed that M. coruscus and Mytilus galloprovincialis possessed the largest number of homologous gene pairs, showing a significant correlation between the two in the evolution of the HECT gene family. Homology analysis with HECT genes of four bivalve species revealed that M. coruscus and M. galloprovincialis possessed the largest number of homologous gene pairs, showing a significant correlation between the two in the evolution of the HECT gene family. M. coruscus exhibited pronounced and specific expression in gills and blood tissues. Notably, Mco_UPL3 gene expression was significantly upregulated after 12 h of acute heat stress (33 °C) and 24 h of Vibrio injection (0.4 OD). Gene ontology analysis of the HECT genes in M. coruscus revealed that it is primarily enriched in protein modification and degradation functions. This suggests that HECT genes may play a key role in protein degradation and immunomodulation in M. coruscus. These findings offer valuable insights for the breeding of stress-tolerant traits in M. coruscus. In summary, our data shed light on the potential functions of HECT E3 ligases in response to heat stress and Vibrio infection, providing practical guidance for enhancing resilience through breeding in M. coruscus.

Spontaneous mutations are the primary source of genetic variation and play a central role in shaping evolutionary processes. To investigate mutational dynamics in Daphnia obtusa, we generated a chromosome-level genome assembly spanning 129.4 Mb across 12 chromosomes, encompassing 15,321 predicted protein-coding genes. Leveraging whole-genome sequencing of eight mutation-accumulation (MA) lines propagated for an average of 482 generations (spanning over 20 years), we estimated a spontaneous single-nucleotide mutation (SNM) rate of 2.45 × 10−9 and an indel mutation rate of 3.34 × 10−10 per site per generation. The SNM spectrum was strongly biased toward C:G > T:A transitions. Despite the design of MA experiments to minimize selection, nonsynonymous mutations were strongly underrepresented, providing rare evidence that purifying selection can act detectably even during mutation accumulation. Comparative analyses with natural population data revealed that exonic mutations observed in the MA lines were significantly less likely to be present in standing variation than intronic or intergenic mutations, suggesting that purifying selection in natural populations acts to remove deleterious alleles. We also identified 48 loss-of-heterozygosity (LOH) events, comprising 8 heterozygous deletions and 40 gene-conversion events, yet found no evidence of GC-biased gene conversion. Instead, while mutation predicts a substantially lower equilibrium GC content, the observed GC level is maintained at higher values, implicating natural selection as the primary force stabilizing base composition. Together, these results provide one of the most comprehensive assessments of the interplay among mutation, selection, and genome stability in an ecologically important species.