Studies of DNA methylomes in a wide range of eukaryotes have revealed both conserved and divergent characteristics of DNA methylation among phylogenetic groups. However, data on invertebrates particularly molluscs are limited, which hinders our understanding of the evolution of DNA methylation in metazoa. The sequencing of the Pacific oyster Crassostrea gigas genome provides an opportunity for genome-wide profiling of DNA methylation in this model mollusc. Homologous searches against the C. gigas genome identified functional orthologs for key genes involved in DNA methylation: DNMT1, DNMT2, DNMT3, MBD2/3 and UHRF1. Whole-genome bisulfite sequencing (BS-seq) of the oyster’s mantle tissues revealed that more than 99% methylation modification was restricted to cytosines in CpG context and methylated CpGs accumulated in the bodies of genes that were moderately expressed. Young repeat elements were another major targets of CpG methylation in oysters. Comparison with other invertebrate methylomes suggested that the 5’-end bias of gene body methylation and the negative correlation between gene body methylation and gene length were the derived features probably limited to the insect lineage. Interestingly, phylostratigraphic analysis showed that CpG methylation preferentially targeted genes originating in the common ancestor of eukaryotes rather than the oldest genes originating in the common ancestor of cellular organisms. Comparative analysis of the oyster DNA methylomes and that of other animal species revealed that the characteristics of DNA methylation were generally conserved during invertebrate evolution, while some unique features were derived in the insect lineage. The preference of methylation modification on genes originating in the eukaryotic ancestor rather than the oldest genes is unexpected, probably implying that the emergence of methylation regulation in these 'relatively young’ genes was critical for the origin and radiation of eukaryotes.

The condition index (CI) is an economically important tool for assessing the quality of oysters, such as the Pacific oyster Crassostrea gigas. However, little is known about the mechanisms that underlie differences in CI between different C. gigas populations. In this study, we integrated transcriptomic and metabolomic profiling to investigate the mechanisms that underlie the differences between high- and low-CI groups in one- and two-year-old populations of C. gigas. The results indicate that differences in CI were associated with the regulation of growth-related genes, the FoxO signaling pathway, and the complex regulation of carbohydrate, lipid, amino acid, and energy metabolism. Moreover, the mechanisms underlying these differences differed between the populations. This study is the first to elucidate the molecular and chemical mechanisms associated with CI, and the results will be helpful for breeding higher quality oysters.

Worldwide, the bivalve aquaculture industry has realized or recognized the potential gains from selective breeding programs using phenotypic and pedigree data. Yesso scallop Patinopecten yessoensis are among the most important commercial shellfish in China. A family-based breeding program to investigate the genetic variations for growth and survival-related traits at suspended and bottom environments of sea-based culture systems was reported in this study. We proposed and conducted a novel phenotyping technique to longitudinally evaluate the shell heights in the yesso scallop. At harvest after rearing for 20 months, the individual shell heights at 6, 10, 16, and 20 months were simultaneously obtained by the growth rings in the outer shells of the animals. Meanwhile, the body weight (BW) of the survived individuals at harvest was also recorded. Variance components and genetic parameters for growth and survival-related traits were estimated using an animal and threshold model, respectively. In the suspended environment, the heritability estimates for BW and shell heights at specific ages ranged from moderate to high (0.328–0.853). The estimated correlations between shell heights at contiguous ages were consistently high, ranging from 0.890 to 0.958 but decreased with increasing intervals between ages (0.496–0.828). The estimated correlations between shell heights at contiguous ages and BW at harvest were similar, ranging from 0.535 to 0.983. The heritability estimates for individual survival at harvest were at a low level of 0.128 by the probit-threshold model. While at the bottom environment, estimates of heritability for growth and survival-related traits were similar but slightly lower than those at the suspended environment. Furthermore, the genetic correlation for BW between the two environments was very small, which probably indicates genotype-by-environment interaction effects for growth in the yesso scallop. The study can provide prior information, which might develop a new idea for selection in this species. The results are discussed concerning selection work with yesso scallop, and solutions for accurate estimation of genetic parameters and increasing genetic gain are also outlined.