Characterization of paramyosin protein structure and gene expression during myogenesis in Pacific oyster (Crassostrea gigas).

Paramyosin is a key component of thick filaments in invertebrate muscles. In this study, we isolated the full length cDNA of paramyosin from Pacific oyster (Crassostrea gigas), and determined its pattern of expression during myogenesis. The full length paramyosin (CgPM) cDNA contains an open reading frame (ORF) of 2586 bp encoding a 861-amino acid protein. Sequence analysis revealed an assembly competence domain (ACD) and a heptad repeat (d-e-f-g-a-b-c) with 28-residue repeat zones in the CgPM primary structure, a characteristic of coiled-coil protein. Quantitative analysis of CgPM expression revealed a sharp increase in trochophore stage, and peaked at the D-shaped stage. Strong CgPM expression was found in smooth adductor muscle, followed by striated adductor muscle and mantle tissue. Whole-mount in situ hybridization (WISH) showed a restricted pattern of CgPM expression in adductor muscle, larval velum retractor and foot muscles at the umbo and eyed larval stages. These data indicate that CgPM is strongly expressed during larval myogenesis in C. gigas, which provides the basis for further functional studies of paramyosin in oyster to better understand the molecular and cellular mechanisms of muscle formation in mollusks.

Cytoskeletal polarization and cytokinetic signaling drives polar lobe formation in spiralian embryos.

In many spiralians, asymmetry in the first two cleavages is achieved through the formation of a polar lobe (PL), which transiently constricts to sequester vegetal cytoplasm into the CD and D blastomeres. While microtubules and actin filaments are required for polar lobe formation, little else is known regarding the structural and functional similarities with the contractile ring, or how the PL constriction is able to form perpendicular to the cleavage plane. Examination of scallop embryos revealed that while activated myosin II could be detected in both the cleavage furrow and early PL constriction, astral or central spindle microtubules were not observed associated with the PL neck until the constriction was nearly complete. Further, inhibition of Aurora B had no effect on polar lobe initiation, but blocked both contractile ring ingression and PL constriction beyond phase II. The cortex destined for PL sequestration was marked by enrichment of the Arp2/3 complex, which was first detected during meiosis and remained enriched at the vegetal pole through the first two cleavages. Inhibition of Arp2/3 affected PL formation and partitioning of cytoplasm into the two daughter cells, suggesting that Arp2/3 plays a functional role in defining the zone of cortex to be sequestered into the polar lobe. Together, these data offer for the first time a mechanism by which a cytoskeletal specialization defines the polar lobe in this atypical form of asymmetric cell division.

Gene Expression Does Not Support the Developmental Hourglass Model in Three Animals with Spiralian Development.

It has been proposed that animals have a pattern of developmental evolution resembling an hourglass because the most conserved development stage-often called the phylotypic stage-is always in midembryonic development. Although the topic has been debated for decades, recent studies using molecular data such as RNA-seq gene expression data sets have largely supported the existence of periods of relative evolutionary conservation in middevelopment, consistent with the phylotypic stage and the hourglass concepts. However, so far this approach has only been applied to a limited number of taxa across the tree of life. Here, using established phylotranscriptomic approaches, we found a surprising reverse hourglass pattern in two molluscs and a polychaete annelid, representatives of the Spiralia, an understudied group that contains a large fraction of metazoan body plan diversity. These results suggest that spiralians have a divergent midembryonic stage, with more conserved early and late development, which is the inverse of the pattern seen in almost all other organisms where these phylotranscriptomic approaches have been reported. We discuss our findings in light of proposed reasons for the phylotypic stage and hourglass model in other systems.

Developmental dynamics of myogenesis in Pacific oyster Crassostrea gigas.

Pacific oyster (Crassostrea gigas) is a sessile bivalve living in the intertidal zone. It has become an attractive model for developmental studies due to its metamorphic transition from a mobile planktonic larvae to a sessile adults. To determine the effect of metamorphosis on muscle development in oyster larvae, we characterized myogenesis during larval development and metamorphosis by phalloidin staining which labels filamentous actin filaments. Our data revealed a dynamic pattern of myogenesis during embryonic and larval development. It appears that simple "U-shaped" muscle ring first developed at the trochophore stage. This was followed by a more complex musculature including an anterior adductor, velum ventral retractors at the veliger stage, and the addition of posterior adductors and foot retractors at the veliger and pediveliger stages. During metamorphosis, muscle structures in the anterior adductor, velum retractors and ventral retractors were degenerated. At the same time, mantle and gill musculature appeared and became the primary muscle system in juveniles together with the posterior adductor. In addition, indirect immunofluorescence with the monoclonal antibody against C. gigas muscle proteins (myosin heavy chains (MYHC) and α-actinin) were used to monitor changes in the developing musculature at different larval stages. The immunofluorescence staining results of muscle proteins were consistent with phalloidin staining. The expression locations of two muscle proteins were similar and mainly located in larval velum retractor and adductor muscle. The α-actinin expression positions were located in Z-lined of velum striated retractors. Data from these studies provide a comprehensive description of myogenesis in C. gigas embryos and larvae. Moreover, our data showed that metamorphosis has a significant impact on remolding the musculature after transition from a mobile planktonic larvae to a sessile mollusk, associated with certain muscle group degradation.

Nanoplastics impaired oyster free living stages, gametes and embryos.

In the marine environment, most bivalve species base their reproduction on external fertilization. Hence, gametes and young stages face many threats, including exposure to plastic wastes which represent more than 80% of the debris in the oceans. Recently, evidence has been produced on the presence of nanoplastics in oceans, thus motivating new studies of their impacts on marine life. Because no information is available about their environmental concentrations, we performed dose-response exposure experiments with polystyrene particles to assess the extent of micro/nanoplastic toxicity. Effects of polystyrene with different sizes and functionalizations (plain 2-µm, 500-nm and 50-nm; COOH-50 nm and NH2-50 nm) were assessed on three key reproductive steps (fertilization, embryogenesis and metamorphosis) of Pacific oysters (Crassostrea gigas). Nanoplastics induced a significant decrease in fertilization success and in embryo-larval development with numerous malformations up to total developmental arrest. The NH2-50 beads had the strongest toxicity to both gametes (EC50 = 4.9 µg/mL) and embryos (EC50 = 0.15 µg/mL), showing functionalization-dependent toxicity. No effects of plain microplastics were recorded. These results highlight that exposures to nanoplastics may have deleterious effects on planktonic stages of oysters, presumably interacting with biological membranes and causing cyto/genotoxicity with potentially drastic consequences for their reproductive success.

Comparative sensitivity of Crassostrea angulata and Crassostrea gigas embryo-larval development to As under varying salinity and temperature.

Oysters are a diverse group of marine bivalves that inhabit coastal systems of the world's oceans, providing a variety of ecosystem services, and represent a major socioeconomic resource. However, oyster reefs have become inevitably impacted from habitat destruction, overfishing, pollution and disease outbreaks that have pushed these structures to the break of extinction. In addition, the increased frequency of climate change related events promise to further challenge oyster species survival worldwide. Oysters' early embryonic development is likely the most vulnerable stage to climate change related stressors (e.g. salinity and temperature shifts) as well as to pollutants (e.g. arsenic), and therefore can represent the most important bottleneck that define populations' survival in a changing environment. In light of this, the present study aimed to assess two important oyster species, Crassostrea angulata and Crassostrea gigas embryo-larval development, under combinations of salinity (20, 26 and 33), temperature (20, 24 and 28 °C) and arsenic (As) exposure (0, 30, 60, 120, 240, 480, 960 and 1920 µg. As L-1), to infer on different oyster species capacity to cope with these environmental stressors under the eminent threat of climate change and increase of pollution worldwide. Results showed differences in each species range of salinity and temperature for successful embryonic development. For C angulata, embryo-larval development was successful at a narrower range of both salinity and temperature, compared to C. gigas. Overall, As induced higher toxicity to C. angulata embryos, with calculated EC50 values at least an order of magnitude lower than those calculated for C. gigas. The toxicity of As (measured as median effective concentration, EC50) showed to be influenced by both salinity and temperature in both species. Nonetheless, salinity had a greater influence on embryos' sensitivity to As. This pattern was mostly noticed for C. gigas, with lower salinity inducing higher sensitivity to As. Results were discussed considering the existing literature and suggest that C. angulata populations are likely to become more vulnerable under near future predictions for temperature rise, salinity shifts and pollution.

Oocyte maturation and origin of the germline as revealed by the expression of Nanos-like in the Pacific oyster Crassostrea gigas.

Nanos gene plays an important role in germline development in animals. However, the molecular mechanisms involved in germline development in Mollusca, the second largest animal phylum, are still poorly understood. Here we identified the Nanos orthologue from the Pacific oyster Crassostrea gigas (Cg-Nanos-like), and investigated the expression patterns of Nanos during gametogenesis and embryogenesis in C. gigas. Tissue expression analysis showed that Cg-Nanos-like was specifically expressed in female gonads. During the reproductive cycle, the expression of Cg-Nanos-like mRNA increased matching the seasonal development of the ovarian tissues in diploids, while the expression levels were significantly lower in the ovaries of sterile triploids compared to diploids. High expression of Cg-Nanos-like transcripts were detected in early embryonic stages, while the expression significantly dropped at gastrulation and was barely detectable in veliger stages. In situ hybridization showed that Cg-Nanos-like was expressed at different stages of developing oocytes, whereas positive signals were detected only in spermatogonia during the spermatogenic cycle. These findings indicated that Cg-Nanos-like was involved in the development of germ cells, and maintenance of oocyte maturation. In early embryogenesis, the transcripts were broadly expressed; following gastrulation, the expression was restricted to two cell clumps, which might be the putative primordial germ cells (PGCs) or their precursors. Based on the results, the formation of the PGCs in C. gigas was consistent with the model of transition from epigenesis to preformation.

Combined effects of arsenic, salinity and temperature on Crassostrea gigas embryotoxicity.

The combined effects of different salinity and temperature levels on the toxicity of Arsenic (As) were studied on the embryonic development of the oyster Crassostrea gigas. A standardized embryotoxicity test was performed to assess the interactive effects of these stressors, in a full factorial design experiment including a range of salinities (15, 19, 24, 28 and 32), temperatures (16, 20, 24, 28 and 32°C) and As concentrations (100, 300, 600, 1200, 2400µgL-1). The embryotoxicity endpoint was about the determination of normal larvae development rates at various conditions, and median effect concentration (EC50) determination for each As exposure condition. Results showed that toxicity induced by As was characterized by retardation of embryonic development observing toxic effects at lower concentrations than previously reported studies. The presence of As in seawater resulted in a narrower range of tolerance to both salinity and temperature. These findings bring new insights on the impacts of a common contaminant on an important shellfish species having a planktonic early life stage development, with potential implications for population survival and ecosystem functioning in a changing environment.

Expansion of TALE homeobox genes and the evolution of spiralian development.

Spiralians, including molluscs, annelids and platyhelminths, share a unique development process that includes the typical geometry of early cleavage and early segregation of cell fate in blastomeres along the animal-vegetal axis. However, the molecular mechanisms underlying this early cell fate segregation are largely unknown. Here, we report spiralian-specific expansion of the three-amino-acid loop extension (TALE) class of homeobox genes. During early development, some of these TALE genes are expressed in staggered domains along the animal-vegetal axis in the limpet Nipponacmea fuscoviridis and the polychaete Spirobranchus kraussii. Inhibition or overexpression of these genes alters the developmental fate of blastomeres, as predicted by the gene expression patterns. These results suggest that the expansion of novel TALE genes plays a critical role in the establishment of a novel cell fate segregation mechanism in spiralians.

Individual and Combined Toxicities of Benzo[a]pyrene and 2,2',4,4'-Tetrabromodiphenyl Ether on Early Life Stages of the Pacific Oyster, Crassostrea gigas.

Persistent organic pollutants (POPs) are ubiquitous and coexisted in the aquatic environment. Individual and combined toxic effects of benzo[a]pyrene (BaP) and 2,2',4,4'-tetrabromodiphenyl ether (BDE-47) on embryogenesis, and larval survival of the Pacific oyster were investigated. The EC50 values of BaP, BDE-47 and their mixture on embryogenesis were 18.4, 203.3 and 72.0 µg/L respectively, while the LC50 values for 96 h larval mortality were 26.8, 244.5 and 108.9 µg/L respectively. The Marking-Dawson additive toxicity indices were -0.02 and -0.19, indicating an additive effect with a trend to antagonism. In addition, DNA strand breaks were also observed in oyster embryos after exposure. Our study suggests that BaP and BDE-47 exposure can cause developmental abnormalities, DNA damage and larval mortality. Furthermore, the toxicity of the mixture is slightly lower than individual pollutant. These data will be helpful to predict the toxicity of organic pollutants, and provide criteria for marine water quality standards.

Effects of phenanthrene on early development of the Pacific oyster Crassostrea gigas (Thunberg, 1789).

Phenanthnere (PHE) is a polycyclic aromatic hydrocarbon continuously discarded in the marine environment and bioavailable to many aquatic species. Although studies about PHE toxicity have been documented for adult oysters, the effects on early developmental stages are poorly characterized in bivalves. In this study, the effects of PHE (0.02 and 2.0µg.L-1) were evaluated on the embryogenesis and larval development of Crassostrea gigas. Toxicity bioassays, growth and deformities assessment, analysis of shell calcium abundance and transcript levels of genes related to xenobiotic biotransformation (CYP2AU2, CYP30C1), immune system (Cg-Tal) and tissue growth and shell formation (Ferritin, Insulin-like, Cg-Try, Calmodulin and Nacrein) were assayed in D-shape larvae after 24h of PHE exposure. At the highest concentration (2.0µg.L-1), PHE decreased the frequency of normal development (19.7±2.9%) and shell size (53.5±2.8mm). Developmental deformities were mostly related to abnormal mantle and shell formation. Lower calcium levels in oyster shells exposed to PHE 2.0µg.L-1 were observed, suggesting effects on shell structure. At this same PHE concentration, CYP30C1, Cg-Tal, Cg-Tyr, Calmodulin were upregulated and CYP2AU2, Ferritin, Nacrein, and Insulin-Like were downregulated compared to control larvae. At the lowest PHE concentration (0.02µg.L-1), it was observed a minor decrease in normal larval development (89,6±6%) and the remaining parameters were not affected. This is the first study to provide evidences that exposure to PHE can affect early oyster development at the molecular and morphological levels, possibly threatening this bivalve species.

Molecular characterization and expression profiles of myosin essential light chain gene in the Pacific oyster Crassostrea gigas.

In molluscs, muscle contraction is triggered by a direct binding of Ca2+ to myosin. The myosin essential light chain (MELC) provides the Ca2+ binding site and is of importance for motor function and regulation. In this study, the complete cDNA sequence of MELC gene of the Pacific oyster, Crassostrea gigas, was obtained, and the expression profiles were performed in different tissues and different embryo-larval development stages. The results showed that the full length of C. gigas MELC (CgMELC) cDNA was 659bp, containing a 5'-untranslated region of 73 bp, a 3'-untranslated region of 112bp. The open reading frame encoded a 157 amino acid peptide. The protein sequence of CgMELC contained a conserved EF-hand calcium binding motif, and showed a high sequence identity (68.4-100%) with other bivalves. Quantitative analysis of CgMELC mRNA in tissues indicated that CgMELC was expressed at the highest level in the striated adductor muscle, followed by the smooth adductor muscle and mantle. During the development of embryos and larvae, quantitative analysis and whole-mount in situ hybridization revealed that CgMELC was expressed starting from the blastula stage and abundantly expressed in trochophore and D-shaped larvae, indicating that CgMELC might also be involved in regulation of larval muscle system development. Our data provided valuable information for further researches on the function of MELC and regulation of muscle contraction in oysters.

Expression patterns indicate that BMP2/4 and Chordin, not BMP5-8 and Gremlin, mediate dorsal-ventral patterning in the mollusk Crassostrea gigas.

Though several bilaterian animals use a conserved BMP2/4-Chordin antagonism to pattern the dorsal-ventral (DV) axis, the only lophotrochozoan species in which early DV patterning has been studied to date, the leech Helobdella robusta, appears to employ BMP5-8 and Gremlin. These findings call into question the conservation of a common DV patterning mechanism among bilaterian animals. To explore whether the unusual DV patterning mechanism in H. robusta is also used in other lophotrochozoan species, we investigated the expression of orthologous genes in the early embryo of a bivalve mollusk, Crassostrea gigas. Searching of the genome and phylogenetic analysis revealed that C. gigas possesses single orthologs of BMP2/4, Chordin, and BMP5-8 and no Gremlin homolog. Whole mount in situ hybridization revealed mRNA localization of BMP2/4 and Chordin on the opposite sides of embryos, suggesting the potential involvement of a BMP2/4-Chordin antagonism in DV patterning in this species. Furthermore, universal BMP5-8 expression and the absence of a Gremlin homolog in the C. gigas genome called into question any major contribution by BMP5-8 and Gremlin to early DV patterning in this species. Additionally, we identified seven genes showing asymmetric expression along the DV axis, providing further insight into DV patterning in C. gigas. We present the first report of a Chordin gene in a lophotrochozoan species and of the opposite expression of BMP2/4 (dorsal) and Chordin (ventral) along the D/V axis of a lophotrochozoan embryo. The findings of this study further the knowledge of axis formation in lophotrochozoan species and provide insight into the evolution of the animal DV patterning mechanism.

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.

Screening for contaminant hotspots in the marine environment of Kuwait using ecotoxicological and chemical screening techniques.

Kuwait is a country with low rainfall and highly concentrated industrial and domestic effluents entering its coastal waters. These can be both treated and untreated. In this study we sampled a series of coastal and open-sea sites and used a variety of analyses to identify those sites requiring the most attention. We used a high throughput GC-MS screen to look for over 1000 chemicals in the samples. Estrogen and androgen screens assessed the potential to disrupt endocrine activity. An oyster embryo development screen was used to assess biological effect potential. The chemical screen identified sites which had high numbers of identified industrial and domestic chemicals. The oyster screen showed that these sites had also caused high levels of developmental abnormalities with 100% of embryos affected at some sites. The yeast screen showed that estrogenic chemicals were present in outfalls at 2-3 ng/l E2 equivalent, and detectable even in some open water sites.

Impacts of Deepwater Horizon oil and associated dispersant on early development of the Eastern oyster Crassostrea virginica.

The explosion of the Deepwater Horizon (DWH) oil platform resulted in large amounts of crude oil and dispersant Corexit 9500A® released into the Gulf of Mexico and coincided with the spawning season of the oyster, Crassostrea virginica. The effects of exposing gametes and embryos of C. virginica to dispersant alone (Corexit), mechanically (HEWAF) and chemically dispersed (CEWAF) DWH oil were evaluated. Fertilization success and the morphological development, growth, and survival of larvae were assessed. Gamete exposure reduced fertilization (HEWAF: EC201h=1650µg tPAH50L(-1); CEWAF: EC201h=19.4µg tPAH50L(-1); Corexit: EC201h=6.9mgL(-1)). CEWAF and Corexit showed a similar toxicity on early life stages at equivalent nominal concentrations. Oysters exposed from gametes to CEWAF and Corexit experienced more deleterious effects than oysters exposed from embryos. Results suggest the presence of oil and dispersant during oyster spawning season may interfere with larval development and subsequent recruitment.

Lightsticks content toxicity: effects of the water soluble fraction on the oyster embryonic development.

Lightsticks are artifacts used as attractors in a type of commercial fishery, known as surface longline gear. Despite the excessive use, the contamination risks of these devices have not yet been properly investigated. This research aimed to fill up this gap by determining the chemical composition and the toxicity of lightsticks recently activated, compared to those one year after activation and to the ones collected on the beaches. The analyzes were carried out by Gas Chromatography coupled with Mass Spectrometry (GC-MS). Additionally, the variations in composition and the toxicity of their sea Water Soluble Fractions (WSF) were evaluated based on the WSF-effects of Crassostrea rhizophorae embryonic development. The GC-MS analysis made possible the identification of nineteen substances in the water soluble fraction of the lightsticks, such as dibutyl phthalate (DBP) and dimethyl phthalate (DMP). The value of the WSF-effective concentration (EC50) was in an average of 0.35%. After one year of the lightsticks activation, the toxicity was even higher (0.65%). Furthermore, other substances, also present in the lightsticks-WSF caused persistent toxicity even more dangerous to the environment than DBP and DMP. This essay discusses their toxicity effects and possible environment damages.

Regulation of Hox orthologues in the oyster Crassostrea gigas evidences a functional role for promoter DNA methylation in an invertebrate.

DNA methylation within promoter regions (PRDM) controls vertebrate early gene transcription and thereby development, but is neglected outside this group. However, epigenetic features in the oyster Crassostrea gigas suggest functional significance of PDRM in invertebrates. To investigate this, reporter constructs containing in vitro methylated oyster Hox gene promoters were transfected into oyster embryos. The influence of in vivo methylation was studied using bisulfite sequencing and DNA methyltransferase inhibition during development. Our results demonstrate that methylation controls the transcriptional activity of the promoters investigated, unraveling a functional role for PRDM in a lophotrochozoan, an important finding regarding the evolution of epigenetic regulation.

The genome of the Pacific oyster Crassostrea gigas brings new insights on the massive expansion of the C1q gene family in Bivalvia.

C1q domain-containing (C1qDC) proteins are regarded as important players in the innate immunity of bivalve mollusks and other invertebrates and their highly adaptive binding properties indicate them as efficient pathogen recognition molecules. Although experimental studies support this view, the molecular data available at the present time are not sufficient to fully explain the great molecular diversification of this family, present in bivalves with hundreds of C1q coding genes. Taking advantage of the fully sequenced genome of the Pacific oyster Crassostrea gigas and more than 100 transcriptomic datasets, we: (i) re-annotated the oyster C1qDC loci, thus identifying the correct genomic organization of 337 C1qDC genes, (ii) explored the expression pattern of oyster C1qDC genes in diverse developmental stages and adult tissues of unchallenged and experimentally treated animals; (iii) investigated the expansion of the C1qDC gene family in all major bivalve subclasses. Overall, we provide a broad description of the functionally relevant features of oyster C1qDC genes, their comparative expression levels and new evidence confirming that a gene family expansion event has occurred during the course of Bivalve evolution, leading to the diversification of hundreds of different C1qDC genes in both the Pteriomorphia and Heterodonta subclasses.

Temperature influences histone methylation and mRNA expression of the Jmj-C histone-demethylase orthologues during the early development of the oyster Crassostrea gigas.

In many groups, epigenetic mechanisms influence developmental gene regulation under environmental inputs. The Pacific oyster Crassostrea gigas belongs to lophotrochozoans and its larval development is highly dependent on temperature, but the role of epigenetic mechanisms in this context is unknown despite high levels of the recently characterized Jumonji histone demethylase (JHDM) orthologues (Cg_Jumonji) suggesting a physiological relevance of histone methylation in the oyster development. Because in other species alterations of the histone methylation pattern have deleterious outcomes, we investigated the influence of temperature during the oyster larval life on histone methylation and JHDM expression. To shed light on this point, oyster embryonic and early larval development experiments were carried out at different temperatures (18 °C, 25 °C and 32 °C). Histone methylation levels were investigated using fluorescent ELISA at 6 and 24h post-fertilization. When compared to the 25 °C group, at 18 °C H3K4, H3K9 and H3K27 residues were hypomethylated at 6h post fertilization (hpf) and hypermethylated at 24 hpf. In contrast, at 32 °C, 6hpf animals present a dramatic hypermethylation (ca. 4-fold) of all examined residues, which is minored but sustained at 24 hpf. RT-qPCR investigations of the mRNA expression of the nine oyster JHDMs, showed gene- and stage-specific temperature sensitivities throughout the early life of oysters. This study provides evidence of the biological significance of histone methylation during development in a lophotrochozoan species. Our results also indicate that temperature influences histone methylation, possibly through the expression level of putative actors of its regulation, which might participate in developmental control. To our knowledge, this is the first report indicating a direct relationship between an epigenetic mark and an environmental parameter in marine molluscs. Such investigations could help better understand the molecular mechanisms of development and adaptation in lophotrochozoans.