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.
The Pacific oyster Crassostrea gigas belongs to one of the most species-rich but genomically poorly explored phyla, the Mollusca. Here we report the sequencing and assembly of the oyster genome using ...short reads and a fosmid-pooling strategy, along with transcriptomes of development and stress response and the proteome of the shell. The oyster genome is highly polymorphic and rich in repetitive sequences, with some transposable elements still actively shaping variation. Transcriptome studies reveal an extensive set of genes responding to environmental stress. The expansion of genes coding for heat shock protein 70 and inhibitors of apoptosis is probably central to the oyster's adaptation to sessile life in the highly stressful intertidal zone. Our analyses also show that shell formation in molluscs is more complex than currently understood and involves extensive participation of cells and their exosomes. The oyster genome sequence fills a void in our understanding of the Lophotrochozoa.
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.
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.
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.
Overall, this study provides insights into the genetic correlation among complex traits, which will facilitate future oyster functional studies and breeding through molecular design.
Oysters are notably susceptible to accumulating substantial amounts of cadmium (Cd), leading to food safety concerns, yet the molecular mechanisms are not fully understood. Previous study identified ...a significant correlation between polymorphisms in the Zip1 and Zip3 of Fujian oyster ( Crassostrea angulata ) and the levels of Cd accumulation. Zip1 and Zip3 , belonging to the Zip family, are primarily responsible for the cellular transport of zinc and other metal ions. To delve deeper into their roles in Cd transport within oysters, this research cloned and analyzed the sequences of Fujian oyster Zip1 and Zip3 . The open reading frames for Zip1 and Zip3 were 960 and 978 bp, respectively, encoding proteins of 319 and 325 amino acids. Following Cd exposure, the expression levels of these genes in the mantle, gills, and other tissues initially increased significantly before decreasing, demonstrating a dose-dependent response. The lowest expression levels were observed in the adductor muscle, with Zip1 reaching its highest expression in the mantle and Zip3 in the visceral mass. In vitro RNAi experiments revealed that silencing Zip1 led to a 1.26-fold decrease in new 113 Cd accumulation in the mantle, whereas silencing Zip3 resulted in a 1.29-fold increase. Immunohistochemistry showed weak expression of Zip1 and Zip3 proteins in the connective tissue of the mantle. In summary, Zip1 and Zip3 in Fujian oyster can transport Cd as a substrate, and played different roles in the absorption and efflux of Cd.
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.
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
HSP
s 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.
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.
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.
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.
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.
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.