Investigating the Impact of Humic Acid on Copper Accumulation in Sinonovacula constricta Using a Toxicokinetic-Toxicodynamic Model.

In aquatic ecosystems, the interaction between heavy metals and dissolved organic carbon (DOC) plays a pivotal role in modifying the bioavailability of these metals. This study, employing a toxicokinetic-toxicodynamic model, delves into the interactive effects of humic acid (HA), a significant component of DOC, on the bioaccumulation and toxicity of copper (Cu) in the estuarine economic bivalve Sinonovacula constricta. Utilizing the stable isotope 65Cu as a tracer, we evaluated Cu uptake in S. constricta under varied DOC concentrations in a controlled laboratory setting. Our findings reveal that at DOC concentrations below 3.05 mg L-1, the bioavailability of Cu is reduced due to shifts in the speciation distribution of Cu, resulting in decreased bioaccumulation within S. constricta. Conversely, at DOC levels exceeding 3.05 mg L-1, the formation of colloidal Cu-HA complexes allows its entry into the bivalves' digestive system. Moreover, toxicity assays demonstrate an increase in S. constricta survival rates with higher DOC concentrations, suggesting a protective effect of DOC against Cu toxicity. The integration of accumulation and toxicity data infers that Cu-HA complexes, when ingested via the digestive tract, exhibit lower toxicity compared to Cu directly assimilated from the water phase. These findings emphasize the need to consider environmental DOC levels in assessing Cu pollution risks and provide insights for managing heavy metal toxicity in estuarine aquaculture.

DNA methylation mediates differentiation in thermal responses of Pacific oyster (Crassostrea gigas) derived from different tidal levels.

Epigenetic mechanisms such as DNA methylation have the potential to affect organism acclimatization and adaptation to environmental changes by influencing their phenotypic plasticity; however, little is known about the role of methylation in the adaptive phenotypic divergence of marine invertebrates. Therefore, in this study, a typical intertidal species, the Pacific oyster (Crassostrea gigas), was selected to investigate the epigenetic mechanism of phenotypic plasticity in marine invertebrates. Intertidal and subtidal oysters subjected to one-generation common garden experiments and exhibited phenotypic divergence were used. The methylation landscape of both groups of oysters was investigated under temperate and high temperature. The two tidal oysters exhibited divergent methylation patterns, regardless of the temperature, which was mainly original environment-induced. Intertidal samples exhibited significant hypomethylation and more plasticity of methylation in response to heat shock, while subtidal samples showed hypermethylation and less plasticity. Combined with RNA-seq data, a positive relationship between methylation and expression in gene bodies was detected on a genome-wide scale. In addition, approximately 11% and 7% of differentially expressed genes showed significant methylation variation under high temperatures in intertidal and subtidal samples, respectively. Genes related to apoptosis and organism development may be regulated by methylation in response to high temperature in intertidal oysters, whereas oxidation-reduction and ion homeostasis-related genes were involved in subtidal oysters. The results also suggest that DNA methylation mediates phenotypic divergence in oysters adapting to different environments. This study provides new insight into the epigenetic mechanisms underlying phenotypic plasticity in adaptation to rapid climate change in marine organisms.

Regulation Between HSF1 Isoforms and HSPs Contributes to the Variation in Thermal Tolerance Between Two Oyster Congeners.

Heat shock transcription factor 1 (HSF1) plays an important role in regulating heat shock, which can activate heat shock proteins (HSPs). HSPs can protect organisms from thermal stress. Oysters in the intertidal zone can tolerate thermal stress. The Pacific oyster (Crassostrea gigas gigas) and Fujian oyster (C. gigas angulata)-allopatric subspecies with distinct thermal tolerances-make good study specimens for analyzing and comparing thermal stress regulation. We cloned and compared HSF1 isoforms, which is highly expressed under heat shock conditions in the two subspecies. The results revealed that two isoforms (HSF1a and HSF1d) respond to heat shock in both Pacific and Fujian oysters, and different heat shock conditions led to various combinations of isoforms. Subcellular localization showed that isoforms gathered in the nucleus when exposed to heat shock. The co-immunoprecipitation revealed that HSF1d can be a dimer. In addition, we selected HSPs that are expressed under the heat shock response, according to the RNA-seq and proteomic analyses. For the HSPs, we analyzed the coding part and the promoter sequences. The result showed that the domains of HSPs are conserved in two subspecies, but the promoters are significantly different. The Dual-Luciferase assay showed that the induced expression isoform HSF1d had the highest activity in C. gigas gigas, while the constitutively-expressed HSF1a was most active in C. gigas angulata. In addition, variation in the level of HSP promoters appeared to be correlated with gene expression. We argue that this gene is regulated based on the different expression levels between the two subspecies' responses to heat shock. In summary, various stress conditions can yield different HSF1 isoforms and respond to heat shock in both oyster subspecies. Differences in how the isoforms and promoter are activated may contribute to their differential expressions. Overall, the results comparing C. gigas gigas and C. gigas angulata suggest that these isoforms have a regulatory relationship under heat shock, providing valuable information on the thermal tolerance mechanism in these commercially important oyster species.

Functional characterization of retinoid X receptor with an emphasis on the mediation of organotin poisoning in the Pacific oyster (Crassostrea gigas).

Marine mollusks suffer harmful effects due to environmental organotin compounds such as tributyltin (TBT) and triphenyltin (TPT). It is known that gastropod imposex caused by organotins is mediated by a key nuclear receptor, retinoid X receptor (RXR). The organotin-mediated toxic effects on oysters grown in seawater include a thicker shell, incomplete growth, disrupted development and a high rate of mortality. However, few studies have been conducted to determine the role of RXR in the toxic effects of organotins on bivalves. Here, we cloned an RXR homolog (CgRXR) from the Pacific oyster (Crassostrea gigas) and characterized its molecular function. Expression of the CgRXR RNA transcripts was assessed in whole developmental stages and tissues, with the highest expression detected in the blastula and mantle, respectively. The subcellular localization experiment confirmed that CgRXR protein was expressed in the nucleus exclusively as a nuclear receptor. Electrophoretic mobility shift assay indicated that CgRXR could bind to the DNA motifs DR0-DR5. The dual-luciferase reporter assay demonstrated that the transcriptional activity of CgRXR was activated by conserved ligands (9-cis retinoic acid and cis-4,7,10,13,16,19-docosahexanoic acid) and endocrine-disrupting chemicals (TBT and TPT). These results revealed the conserved gene function involved in protein localization, ligand binding and heterodimer formation with thyroid hormone receptor. However, the DNA binding properties of CgRXR differed from those of other invertebrate and vertebrate RXRs. CgRXR had the highest expression level in the blastula and mantle, and the disrupted development or shell malformation induced by organotins suggested a possible correlation of CgRXR with shell formation in bivalves. The results indicated the potential involvement of CgRXR in the toxic effects of organotins (TBT and TPT) through signaling pathway disruption. Functional characterization of CgRXR will help us better understand the endocrinology of bivalves.

Identification of SNPs involved in Zn and Cu accumulation in the Pacific oyster (Crassostrea gigas) by genome-wide association analysis.

Oysters accumulate high concentrations of zinc (Zn) and copper (Cu), which can be transferred to human due to sea food consumption. Breeding new oyster varieties with low Zn and Cu accumulations is one important way to improve food safety. However, the genetic basis for metal accumulation in mollusks is not well understood. To address this issue, oysters collected in the field were used for genome-wide association study (GWAS) and then the identified genes were used for mRNA expressions analysis in laboratory. First, GWAS were conducted for Zn and Cu accumulation in 288 wild Pacific oysters (Crassostrea gigas) farmed in the same ocean environment. The oysters did not show obvious population structure or kinship but exhibited 8.43- and 10.0- fold changes of Zn and Cu contents respectively. GWAS have identified 11 and 12 single nucleotide polymorphisms (SNPs) associated with Zn and Cu, respectively, as well as 16 genes, which were Zn-containing proteins or participated in caveolae-dependent endocytosis. Second, the mRNA expressions of these 16 genes were observed under Zn and Cu exposure. After 9 days of Zn exposure, Zn contents increased 3.1-fold, while the mRNA expression of cell number regulator 3 increased 1.65-fold. Under 9 days of Cu exposure, Cu contents increased 1.97-fold, while the mRNA expression of caveolin-1 decreased 0.61-fold. These provide the evidence for their roles in regulating physiological levels of these two metals. The findings advance our understanding of the genetic basis of Zn and Cu accumulation in mollusks, which can be useful for breeding new, less toxic varieties of oysters.

Genome-wide association analysis of nutrient traits in the oyster Crassostrea gigas: genetic effect and interaction network.

BACKGROUND: Oyster is rich in glycogen and free amino acids and is called "the milk of sea". To understand the main genetic effects of these traits and the genetic networks underlying their correlation, we have conducted the whole genome resequencing with 427 oysters collected from the world-wide scale. RESULTS: After association analysis, 168 clustered significant single nucleotide polymorphism (SNP) loci were identified for glycogen content and 17 SNPs were verified with 288 oyster individuals in another wide populations. These were the most important candidate loci for oyster breeding. Among 24 genes in the 100-kb regions of the leading SNP loci, cytochrome P450 17A1 (CYP17A1) contained a non-synonymous SNP and displayed higher expressions in high glycogen content individuals. This might enhance the gluconeogenesis process by the transcriptionally regulating the expression of phosphoenolpyruvate carboxykinase (PEPCK) and glucose 6-phosphatase (G6Pase). Also, for amino acids content, 417 clustered significant SNPs were identified. After genetic network analysis, three node SNP regions were identified to be associated with glycogen, protein, and Asp content, which might explain their significant correlation. CONCLUSION: Overall, this study provides insights into the genetic correlation among complex traits, which will facilitate future oyster functional studies and breeding through molecular design.

Divergence and plasticity shape adaptive potential of the Pacific oyster.

The interplay between divergence and phenotypic plasticity is critical to our understanding of a species' adaptive potential under rapid climate changes. We investigated divergence and plasticity in natural populations of the Pacific oyster Crassostrea gigas with a congeneric oyster Crassostrea angulata from southern China used as an outgroup. Genome re-sequencing of 371 oysters revealed unexpected genetic divergence in a small area that coincided with phenotypic divergence in growth, physiology, heat tolerance and gene expression across environmental gradients. These findings suggest that selection and local adaptation are pervasive and, together with limited gene flow, influence population structure. Genes showing sequence differentiation between populations also diverged in transcriptional response to heat stress. Plasticity in gene expression is positively correlated with evolved divergence, indicating that plasticity is adaptive and favoured by organisms under dynamic environments. Divergence in heat tolerance-partly through acetylation-mediated energy depression-implies differentiation in adaptive potential. Trade-offs between growth and survival may play an important role in local adaptation of oysters and other marine invertebrates.

Population resequencing reveals candidate genes associated with salinity adaptation of the Pacific oyster Crassostrea gigas.

The Pacific oyster Crassostrea gigas is an important cultivated shellfish. As a euryhaline species, it has evolved adaptive mechanisms responding to the complex and changeable intertidal environment that it inhabits. To investigate the genetic basis of this salinity adaptation mechanism, we conducted a genome-wide association study using phenotypically differentiated populations (hyposalinity and hypersalinity adaptation populations, and control population), and confirmed our results using an independent population, high-resolution melting, and mRNA expression analysis. For the hyposalinity adaptation, we determined 24 genes, including Cg_CLCN7 (chloride channel protein 7) and Cg_AP1 (apoptosis 1 inhibitor), involved in the ion/water channel and transporter mechanisms, free amino acid and reactive oxygen species metabolism, immune responses, and chemical defence. Three SNPs located on these two genes were significantly differentiated between groups, as was Cg_CLCN7. For the hypersalinity adaptation, the biological process for positive regulating the developmental process was enriched. Enriched gene functions were focused on transcriptional regulation, signal transduction, and cell growth and differentiation, including calmodulin (Cg_CaM) and ficolin-2 (Cg_FCN2). These genes and polymorphisms possibly play an important role in oyster hyposalinity and hypersalinity adaptation. They not only further our understanding of salinity adaptation mechanisms but also provide markers for highly adaptable oyster strains suitable for breeding.

Characterization and functional analysis of two inhibitor of apoptosis genes in Zhikong scallop Chlamys farreri.

The proteins of inhibitor of apoptosis (IAP) family play important roles in regulation of apoptosis, immunological response and cell proliferation. Here we reported two IAP genes (named CfIAP1 and CfIAP2) in Zhikong scallop Chlamys farreri. The full-length CfIAP1 cDNA contained 1552 nucleotides, encoding a predicted protein of 251 amino acids with two BIR domains. The full-length CfIAP2 cDNA contained 1243 nt, encoding a 356-aa protein with one BIR domain and one RING domain. The two genes are ubiquitously expressed in six types of tissue of C. farreri. The expression levels of CfIAP1 and CfIAP2 were significantly up-regulated after challenged with acute viral necrobiotic disease virus, lipopolysaccharide and exposure to air. Subcellular localization assay showed that CfIAP1 was mainly distributed in cytoplasm and CfIAP2 was in cytoplasm and nucleus. As assessed using a kit designed to test Caspase3 function in mammalian cells, the activity of CfCaspase3 was enhanced as a result of the down-regulation of CfIAP2 expression by dsRNA-mediated gene silencing. Our study indicated that CfIAP1 and CfIAP2 may participate in the innate immunity and stress responses and that CfIAP2 might block apoptosis via inhibiting CfCaspase3 indirectly through an unexplored mechanism in C. farreri.

Expression Characterization of Stress Genes Under High and Low Temperature Stresses in the Pacific Oyster, Crassostrea gigas.

As a characteristic sessile inhabitant of the intertidal zone, the Pacific oyster Crassostrea gigas occupies one of the most physically stressful environments on earth. With high exposure to terrestrial conditions, oysters must tolerate broad fluctuations in temperature range. However, oysters' cellular and molecular responses to temperature stresses have not been fully characterized. Here, we analyzed oyster transcriptome data under high and low temperatures. We also identified over 30 key temperature stress-responsive candidate genes, which encoded stress proteins such as heat shock proteins and apoptosis-associated proteins. The expression characterization of these genes under short-term cold and hot environments (5 and 35 °C) and long-term cold environments (5 °C) was detected by quantitative real-time PCR. Most of these genes reached expression peaks during the recovery stage after 24 h of heat stress, and these genes were greatly induced around day 3 in long-term cold stress while responded little to short-term cold stress. In addition, in the second heat stress after 2 days of recovery, oysters showed milder expression in these genes and a lower mortality rate, which indicated the existence of plasticity in the oyster's response to heat stress. We confirmed that homeostatic flexibility and anti-apoptosis might be crucial centers of temperature stress responses in oysters. Furthermore, we analyzed stress gene families in 11 different species and found that the linage-specific expansion of stress genes might be implicated in adaptive evolution. These results indicated that both plasticity and evolution played an important role in the stress response adaptation of oysters.

Identification of Thyroid Hormones and Functional Characterization of Thyroid Hormone Receptor in the Pacific Oyster Crassostrea gigas Provide Insight into Evolution of the Thyroid Hormone System.

Thyroid hormones (THs) play important roles in development, metamorphosis, and metabolism in vertebrates. During the past century, TH functions were regarded as a synapomorphy of vertebrates. More recently, accumulating evidence has gradually convinced us that TH functions also occur in invertebrate chordates. To date, however, TH-related studies in non-chordate invertebrates have been limited. In this study, THs were qualitatively detected by two reliable methods (HPLC and LC/MS) in a well-studied molluscan species, the Pacific oyster Crassostrea gigas. Quantitative measurement of THs during the development of C. gigas showed high TH contents during embryogenesis and that oyster embryos may synthesize THs endogenously. As a first step in elucidating the TH signaling cascade, an ortholog of vertebrate TH receptor (TR), the most critical gene mediating TH effects, was cloned in C. gigas. The sequence of CgTR has conserved DNA-binding and ligand-binding domains that normally characterize these receptors. Experimental results demonstrated that CgTR can repress gene expression through binding to promoters of target genes and can interact with oyster retinoid X receptor. Moreover, CgTR mRNA expression was activated by T4 and the transcriptional activity of CgTR promoter was repressed by unliganded CgTR protein. An atypical thyroid hormone response element (CgDR5) was found in the promoter of CgTR, which was verified by electrophoretic mobility shift assay (EMSA). These results indicated that some of the CgTR function is conserved. However, the EMSA assay showed that DNA binding specificity of CgTR was different from that of the vertebrate TR and experiments with two dual-luciferase reporter systems indicated that l-thyroxine, 3,3',5-triiodothyronine, and triiodothyroacetic acid failed to activate the transcriptional activity of CgTR. This is the first study to functionally characterize TR in mollusks. The presence of THs and the functions of CgTR in mollusks contribute to better understanding of the evolution of the TH system.

Cloning and characterization of a novel caspase-8-like gene in Crassostrea gigas.

Cysteine-dependent aspartate-directed proteases, or caspases, play key roles in apoptosis and immune defense. In this study, we cloned the first caspase-8-like gene (CgCaspase8-2) identified in the pacific oyster, Crassostrea gigas. The 2572-bp cDNA encodes a putative protein of 714 amino acids that contains two tandem death effector domains (DEDs) at the N-terminal, and P20 and P10 domains at the C-terminal. The conserved pentapeptide motif QACQG was also identified in the deduced CgCaspase8-2 protein. Phylogenetic analysis indicated that CgCaspase8-2 was clustered with initiator caspases in the invertebrate subgroup, but the similarity between CgCaspase8-2 and other invertebrate caspase-8s was low. CgCaspase8-2 protein was localized in the cytoplasm, and over-expression of CgCaspase8-2 in HEK293T cells induced cell death, suggesting a role in apoptosis. Quantitative real-time PCR results demonstrated that CgCaspase8-2 was widely expressed in various tissues and developmental stages, with the highest CgCaspase8-2 expression levels detected in hemolymph and the blastula stage. Furthermore, CgCaspase8-2 transcripts showed no change in response to a bacterial challenge but exhibited notable up-regulation post-poly (I:C) challenge, suggesting that CgCaspase8-2 is specifically involved in immune responses against viruses. In summary, CgCaspase8-2 is involved in both apoptotic and immune function.

Evolution of a novel nuclear receptor subfamily with emphasis on the member from the Pacific oyster Crassostrea gigas.

Nuclear receptors (NRs) belong to the transcription factor superfamily that regulates development, homeostasis, differentiation, and reproduction in metazoans via control of gene expression. Recently, rapid advances in genome projects on various metazoans have provided new opportunities for studying the evolution and function of NRs. Typically structured NRs are divided into six subfamilies. Here, the gene for a typically structured NR (CgNR8A1) was cloned from the Pacific oyster Crassostrea gigas. However, this novel receptor could not be assigned to a known NR subfamily. By data mining, nine other CgNR8A1 gene homologs were identified in metazoans such as cnidarians, mollusks, annelids, echinoderms, hemichordates, and cephalochordates. Phylogenetic analysis showed that these receptors belonged to a novel NR subfamily, hereafter designated as NR8. Evolutionary analysis revealed that the NR8 subfamily was phylogenetically the third-oldest NR subfamily, and it originated from a common ancestor of Eumetazoa; several gene loss events occurred independently in ancestors of vertebrates, ecdysozoans, and platyhelminths, which do not have NR8 members. Furthermore, the function of CgNR8A1 was investigated to provide an insight into the functions of this novel NR subfamily. A nuclear localization signal peptide, GKHRNKKPRLD, was identified in CgNR8A1, and a recombinant full-length protein of CgNR8A1 was localized in the nuclei of HeLa cells. The mRNA expression profile of CgNR8A1 suggested that it might be involved in the embryogenesis of C. gigas. The electrophoretic mobility shift assay showed that CgNR8A1 binds strongly to conserved DNA core motifs DR0, DR2, and DR4 and weakly to DR1, DR3, DR5, Half, and Pal0. In summary, the novel NR8 subfamily identified in this study improves our understanding of NR evolution, and the functional analysis of CgNR8A1 provided further insights into the functions of NR8A1s.

Candidate Gene Polymorphisms and their Association with Glycogen Content in the Pacific Oyster Crassostrea gigas.

BACKGROUND: The Pacific oyster Crassostrea gigas is an important cultivated shellfish that is rich in nutrients. It contains high levels of glycogen, which is of high nutritional value. To investigate the genetic basis of this high glycogen content and its variation, we conducted a candidate gene association analysis using a wild population, and confirmed our results using an independent population, via targeted gene resequencing and mRNA expression analysis. RESULTS: We validated 295 SNPs in the 90 candidate genes surveyed for association with glycogen content, 86 of were ultimately genotyped in all 144 experimental individuals from Jiaonan (JN). In addition, 732 SNPs were genotyped via targeted gene resequencing. Two SNPs (Cg_SNP_TY202 and Cg_SNP_3021) in Cg_GD1 (glycogen debranching enzyme) and one SNP (Cg_SNP_4) in Cg_GP1 (glycogen phosphorylase) were identified as being associated with glycogen content. The glycogen content of individuals with genotypes TT and TC in Cg_SNP_TY202 was higher than that of individuals with genotype CC. The transcript abundance of both glycogen-associated genes was differentially expressed in high glycogen content and low glycogen content individuals. CONCLUSIONS: This study identified three polymorphisms in two genes associated with oyster glycogen content, via candidate gene association analysis. The transcript abundance differences in Cg_GD1 and Cg_GP1 between low- and the high-glycogen content individuals suggests that it is possible that transcript regulation is mediated by variations of Cg_SNP_TY202, Cg_SNP_3021, and Cg_SNP_4. These findings will not only provide insights into the genetic basis of oyster quality, but also promote research into the molecular breeding of oysters.

Identification and functional characterization of two Bcl-2 family protein genes in Zhikong scallop Chlamys farreri.

Apoptosis plays significant roles in maintenance of homeostasis, immune defense and development. The Bcl-2 family proteins are important regulators of the intrinsic apoptosis. In the study, we have characterized a Bcl-2-like gene (named CfBcl-2) and a Bax-like gene (named CfBax) from the Zhikong scallop Chlamys farreri. The full-length of the CfBcl-2 cDNA is 944 nucleotides (nt) encoding a putative protein of 225 amino acid residues (aa) that contains four Bcl-2 homology (BH) domains, and the CfBax cDNA is 505 nt encoding a putative protein of 115 aa that contains three Bcl-2 BH domains. Sequence and phylogenetic analysis demonstrate that CfBcl-2 and CfBax present typical domain organization of the corresponding Bcl-2 related proteins and are more similar and clustered with their homologues of other molluscs. The two genes are ubiquitously expressed in six tissues of C. farreri, with the highest expression level of CfBcl-2 in adductor muscle and highest expression level of CfBax in gill. The expressions of CfBcl-2 and CfBax in hemocytes were both significantly up-regulated after an in vivo exposure of scallops to air, injection with lipopolysaccharide and infection with acute viral necrobiotic disease virus, and the expression patterns of the two genes after the three treatments vary in different change magnitude and up-regulation timespan. Yeast two-hybrid assay reveals a direct interaction between the CfBcl-2 and CfBax proteins. These results indicate that the CfBcl-2 and CfBax may participate in the apoptosis-based stress and immune responses against noxious stimulation.

Discovery and validation of genic single nucleotide polymorphisms in the Pacific oyster Crassostrea gigas.

The economic and ecological importance of the oyster necessitates further research on the molecular mechanisms, which both regulate the commercially important traits of the oyster and help it to survive in the variable marine environment. Single nucleotide polymorphisms (SNPs) have been widely used to assess genetic variation and identify genes underlying target traits. In addition, high-resolution melting (HRM) analysis is a potentially powerful method for validating candidate SNPs. In this study, we adopted a rapid and efficient pipeline for the screening and validation of SNPs in the genic region of Crassostrea gigas based on transcriptome sequencing and HRM analysis. Transcriptomes of three wild oyster populations were sequenced using Illumina sequencing technology. In total, 50-60 million short reads, corresponding to 4.5-5.4 Gbp, from each population were aligned to the oyster genome, and 5.8 × 10(5) SNPs were putatively identified, resulting in a predicted SNP every 47 nucleotides on average. The putative SNPs were unevenly distributed in the genome and high-density (≥2%), nonsynonymous coding SNPs were enriched in genes related to apoptosis and responses to biotic stimuli. Subsequently, 1,671 loci were detected by HRM analysis, accounting for 64.7% of the total selected candidate primers, and finally, 1,301 polymorphic SNP markers were developed based on HRM analysis. All of the validated SNPs were distributed into 897 genes and located in 672 scaffolds, and 275 of these genes were stress inducible under unfavourable salinity, temperature, and exposure to air and heavy metals. The validated SNPs in this study provide valuable molecular markers for genetic mapping and characterization of important traits in oysters.

Genome-wide and single-base resolution DNA methylomes of the Pacific oyster Crassostrea gigas provide insight into the evolution of invertebrate CpG methylation.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Identification of conserved and novel microRNAs in the Pacific oyster Crassostrea gigas by deep sequencing.

MicroRNAs (miRNAs) play important roles in regulatory processes in various organisms. To date many studies have been performed in the investigation of miRNAs of numerous bilaterians, but limited numbers of miRNAs have been identified in the few species belonging to the clade Lophotrochozoa. In the current study, deep sequencing was conducted to identify the miRNAs of Crassostrea gigas (Lophotrochozoa) at a genomic scale, using 21 libraries that included different developmental stages and adult organs. A total of 100 hairpin precursor loci were predicted to encode miRNAs. Of these, 19 precursors (pre-miRNA) were novel in the oyster. As many as 53 (53%) miRNAs were distributed in clusters and 49 (49%) precursors were intragenic, which suggests two important biogenetic sources of miRNAs. Different developmental stages were characterized with specific miRNA expression patterns that highlighted regulatory variation along a temporal axis. Conserved miRNAs were expressed universally throughout different stages and organs, whereas novel miRNAs tended to be more specific and may be related to the determination of the novel body plan. Furthermore, we developed an index named the miRNA profile age index (miRPAI) to integrate the evolutionary age and expression levels of miRNAs during a particular developmental stage. We found that the swimming stages were characterized by the youngest miRPAIs. Indeed, the large-scale expression of novel miRNAs indicated the importance of these stages during development, particularly from organogenetic and evolutionary perspectives. Some potentially important miRNAs were identified for further study through significant changes between expression patterns in different developmental events, such as metamorphosis. This study broadened the knowledge of miRNAs in animals and indicated the presence of sophisticated miRNA regulatory networks related to the biological processes in lophotrochozoans.

Aragonite shells are more ancient than calcite ones in bivalves: new evidence based on omics.

Two calcium carbonate crystal polymorphs, aragonite and calcite, are the main inorganic components of mollusk shells. Some fossil evidences suggest that aragonite shell is more ancient than calcite shell for the Bivalvia. But, the molecular biology evidence for the above deduction is absent. In this study, we searched for homologs of bivalve aragonite-related and calcite-related shell proteins in the oyster genome, and found that no homologs of calcite-related shell protein but some homologs of aragonite-related shell proteins in the oyster genome. We explained the results as the new evidence to support that aragonite shells are more ancient than calcite shells in bivalves combined the published biogeological and seawater chemistry data.

Oyster Shell Proteins Originate from Multiple Organs and Their Probable Transport Pathway to the Shell Formation Front.

Mollusk shell is one kind of potential biomaterial, but its vague mineralization mechanism hinders its further application. Mollusk shell matrix proteins are important functional components that are embedded in the shell, which play important roles in shell formation. The proteome of the oyster shell had been determined based on the oyster genome sequence by our group and gives the chance for further deep study in this area. The classical model of shell formation posits that the shell proteins are mantle-secreted. But, in this study, we further analyzed the shell proteome data in combination with organ transcriptome data and we found that the shell proteins may be produced by multiple organs though the mantle is still the most important organ for shell formation. To identify the transport pathways of these shell proteins not in classical model of shell formation, we conducted a shell damage experiment and we determined the shell-related gene set to identify the possible transport pathways from multiple organs to the shell formation front. We also found that there may exist a remodeling mechanism in the process of shell formation. Based on these results along with some published results, we proposed a new immature model, which will help us think about the mechanism of shell formation in a different way.