Predominant intragenic methylation is associated with gene expression characteristics in a bivalve mollusc

DNA Methylation Profiling in Genetically Selected Clarias magur (Hamilton, 1822) Provides Insights into the Epigenetic Regulation of Growth and Development

DNA methylation is an epigenetic alteration that impacts gene expression without changing the DNA sequence affecting an organism's phenotype. This study utilized a reduced representation bisulfite sequencing (RRBS) approach to investigate the patterns of DNA methylation in genetically selected Clarias magur stocks. RRBS generated 249.22 million reads, with an average of 490,120 methylation sites detected in various parts of genes, including exons, introns, and intergenic regions. A total of 896 differentially methylated regions (DMRs) were identified; 356 and 540 were detected as hyper-methylated and hypo-methylated regions, respectively. The DMRs and their association with overlapping genes were explored using whole genome data of magur, which revealed 205 genes in exonic, 210 in intronic, and 480 in intergenic regions. The analysis identified the maximum number of genes enriched in biological processes such as RNA biosynthetic process, response to growth factors, nervous system development, neurogenesis, and anatomical structure morphogenesis. Differentially methylated genes (DMGs) such as myrip, mylk3, mafb, egr3, ndnf, meis2a, foxn3, bmp1a, plxna3, fgf6, sipa1l1, mcu, cnot8, trim55b, and myof were associated with growth and development. The selected DMGs were analyzed using real-time PCR, which showed altered mRNA expression levels. This work offers insights into the epigenetic mechanisms governing growth performance regulation in magur stocks. This work provides a valuable resource of epigenetic data that could be integrated into breeding programs to select high-performing individuals.

Whole-genome DNA methylation profiling revealed epigenetic regulation of NF-κB signaling pathway involved in response to Vibrio alginolyticus infection in the Pacific oyster, Crassostrea gigas

DNA methylation, an essential epigenetic alteration, is tightly linked to a variety of biological processes, such as immune response. To identify the epigenetic regulatory mechanism in Pacific oyster (Crassostrea gigas), whole-genome bisulfite sequencing (WGBS) was conducted on C. gigas at 0 h, 6 h, and 48 h after infection with Vibrio alginolyticus. At 6 h and 48 h, a total of 11,502 and 14,196 differentially methylated regions (DMRs) were identified (p<0.05, FDR<0.001) compared to 0 h, respectively. Gene ontology (GO) analysis showed that differentially methylated genes (DMGs) were significantly enriched in various biological pathways including immunity, cytoskeleton, epigenetic modification, and metabolic processes. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that transcription machinery (ko03021) is one of the most important pathways. Integrated transcriptome and methylome analyses allowed the identification of 167 and 379 DMG-related DEGs at 6 h and 48 h, respectively. These genes were significantly enriched in immune-related pathways, including nuclear factor kappa B (NF-κB) signaling pathway (ko04064) and tumor necrosis factor (TNF) signaling pathway (ko04668). Interestingly, it's observed that the NF-κB pathway could be activated jointly by TNF Receptor Associated Factor 2 (TRAF2) and Baculoviral IAP Repeat Containing 3 (BIRC3, the homolog of human BIRC2) which were regulated by DNA methylation in response to the challenge posed by V. alginolyticus infection. Through this study, we provided insightful information about the epigenetic regulation of immunity-related genes in the C. gigas, which will be valuable for the understanding of the innate immune system modulation and defense mechanism against bacterial infection in invertebrates.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Genome-wide profiling of DNA methylome and transcriptome reveals epigenetic regulation of Urechis unicinctus response to sulfide stress

Sulfide is a well-known environmental pollutant that can have detrimental effects on most organisms. However, few metazoans living in sulfide-rich environments have developed mechanisms to tolerate and adapt to sulfide stress. Epigenetic mechanisms, including DNA methylation, have been shown to play a vital role in environmental stress adaptation. Nevertheless, the precise function of DNA methylation in biological sulfide adaptation remains unclear. Urechis unicinctus, a benthic organism inhabiting sulfide-rich intertidal environments, is an ideal model organism for studying adaptation to sulfide environments. In this study, we conducted a comprehensive analysis of the DNA methylome and transcriptome of U. unicinctus after exposure to 50 μM sulfide. The results revealed dynamic changes in the DNA methylation (5-methylcytosine) landscape in response to sulfide stress, with U. unicinctus exhibiting elevated DNA methylation levels following stress exposure. Integrating differentially expressed genes (DEGs) and differentially methylated regions (DMRs), we identified a crucial role of gene body methylation in predicting gene expression. Furthermore, using a DNA methyltransferase inhibitor, we validated the involvement of DNA methylation in the sulfide stress response and the gene regulatory network influenced by DNA methylation. The results indicated that by modulating DNA methylation levels during sulfide stress, the expression of glutathione S-transferase, glutamyl aminopeptidase, and cytochrome c oxidase could be up-regulated, thereby facilitating the metabolism and detoxification of exogenous sulfides. Moreover, DNA methylation was found to regulate and enhance the oxidative phosphorylation pathway, including NADH dehydrogenase, isocitrate dehydrogenase, and ATP synthase. Additionally, DNA methylation influenced the regulation of Cytochrome P450 and macrophage migration inhibitory factor, both of which are closely associated with oxidative stress and stress resistance. Our findings not only emphasize the role of DNA methylation in sulfide adaptation but also provide novel insights into the potential mechanisms through which marine organisms adapt to environmental changes.

Unravelling epigenetic mechanisms in Cerastoderma edule genome: a comparison of healthy and neoplastic cockles

Cancer is a multifaceted genetic disease characterized by the acquisition of several essential hallmarks. Notably, certain cancers exhibit horizontal transmissibility, observed across mammalian species and diverse bivalves, the latter referred to as hemic neoplasia. Within this complex landscape, epigenetic mechanisms such as histone modifications and cytosine methylation emerge as fundamental contributors to the pathogenesis of these transmissible cancers. Our study delves into the epigenetic landscape of Cerastoderma edule, focusing on whole-genome methylation and hydroxymethylation profiles in heathy specimens and transmissible neoplasias by means of Nanopore long-read sequencing. Our results unveiled a global hypomethylation in the neoplastic specimens compared to their healthy counterparts, emphasizing the role of DNA methylation in these tumorigenic processes. Furthermore, we verified that intragenic CpG methylation positively correlated with gene expression, emphasizing its role in modulating transcription in healthy and neoplastic cockles, as also highlighted by some up-methylated oncogenic genes. Hydroxymethylation levels were significantly more elevated in the neoplastic samples, particularly within satellites and complex repeats, likely related to structural functions. Additionally, our analysis also revealed distinct methylation and activity patterns in retrotransposons, providing additional insights into bivalve neoplastic processes. Altogether, these findings contribute to understanding the epigenetic dynamics of bivalve neoplasias and shed light on the roles of DNA methylation and hydroxymethylation in tumorigenesis. Understanding these epigenetic alterations holds promise for advancing our broader understanding of cancer epigenetics.

High temperature influences DNA methylation and transcriptional profiles in sea urchins (Strongylocentrotus intermedius)

BACKGROUND

DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level.

RESULTS

The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation.

CONCLUSIONS

High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

Epigenetic and Genetic Population Structure is Coupled in a Marine Invertebrate

Delineating the relative influence of genotype and the environment on DNA methylation is critical for characterizing the spectrum of organism fitness as driven by adaptation and phenotypic plasticity. In this study, we integrated genomic and DNA methylation data for two distinct Olympia oyster (Ostrea lurida) populations while controlling for within-generation environmental influences. In addition to providing the first characterization of genome-wide DNA methylation patterns in the oyster genus Ostrea, we identified 3,963 differentially methylated loci between populations. Our results show a clear coupling between genetic and epigenetic patterns of variation, with 27% of variation in interindividual methylation differences explained by genotype. Underlying this association are both direct genetic changes in CpGs (CpG-SNPs) and genetic variation with indirect influence on methylation (mQTLs). When comparing measures of genetic and epigenetic population divergence at specific genomic regions this relationship surprisingly breaks down, which has implications for the methods commonly used to study epigenetic and genetic coupling in marine invertebrates.

Epigenetic-associated phenotypic plasticity of the ocean acidification-acclimated edible oyster in the mariculture environment

For marine invertebrates with a pelagic-benthic life cycle, larval exposure to ocean acidification (OA) can affect adult performance in response to another environmental stressor. This carry-over effect has the potential to alter phenotypic traits. However, the molecular mechanisms that mediate "OA"-triggered carry-over effects have not been explored despite such information being key to improving species fitness and management strategies for aquafarming. This study integrated the genome-wide DNA methylome and transcriptome to examine epigenetic modification-mediated carry-over OA impacts on phenotypic traits of the ecologically and commercially important oyster species Crassostrea hongkongensis under field conditions. Larvae of C. hongkongensis were exposed to control pH 8.0 and low pH 7.4 conditions, mimicking near future OA scenario in their habitat, before being outplanted as post-metamorphic juveniles at two mariculture field sites with contrasting environmental stressors for 9 months. The larval carry-over OA effect was found to have persistent impacts on the growth and survival trade-off traits on the outplanted juveniles, although the beneficial or adverse effect depended on the environmental conditions at the outplanted sites. Site-specific plasticity was demonstrated with a diverse DNA methylation-associated gene expression profile, with signal transduction and the endocrine system being the most common and highly enriched functions. Highly methylated exons prevailed in the key genes related to general metabolic and endocytic responses and these genes are evolutionarily conserved in various marine invertebrates in response to OA. These results suggest that oysters with prior larval exposure history to OA had the ability to trigger rapid local adaptive responses via epigenetic modification to cope with multiple stressors in the field.

Comparative Methylome Analysis Reveals Epigenetic Signatures Associated with Growth and Shell Color in the Pacific Oyster, Crassostrea gigas

Fast growth is one of the most important breeding goals for all economic species such as the Pacific oyster (Crassostrea gigas), an aquaculture mollusk with top global production. Although the genetic basis and molecular mechanisms of growth-related traits have been widely investigated in the oyster, the role of DNA methylation involved in growth regulation remains largely unexplored. In this study, we performed a comparative DNA methylome analysis of two selectively bred C. gigas strains with contrasted phenotypes in growth and shell color based on whole-genome bisulfite sequencing (WGBS). Genome-wide profiling of DNA methylation at the single-base resolution revealed that DNA methylations were widely spread across the genome with obvious hotspots, coinciding with the distribution of genes and repetitive elements. Higher methylation levels were observed within genic regions compared with intergenic and promoter regions. Comparative analysis of DNA methylation allowed the identification of 339,604 differentially methylated CpG sites (DMCs) clustering in 27,600 differentially methylated regions (DMRs). Functional annotation analysis identified 11,033 genes from DMRs which were enriched in biological processes including cytoskeleton system, cell cycle, signal transduction, and protein biosynthesis. Integrative analysis of methylome and transcriptome profiles revealed a positive correlation between gene expression and DNA methylation within gene-body regions. Protein-protein interaction (PPI) analysis of differentially expressed and methylated genes allowed for the identification of integrin beta-6 (homolog of human ITGB3) as a hub modulator of the PI3K/Akt signaling pathway that was involved in various growth-related processes. This work provided insights into epigenetic regulation of growth in oysters and will be valuable resources for studying DNA methylation in invertebrates.

Acclimatory gene expression of primed clams enhances robustness to elevated pCO2

Sub-lethal exposure to environmental challenges may enhance ability to cope with chronic or repeated change, a process known as priming. In a previous study, pre-exposure to seawater enriched with pCO₂ improved growth and reduced antioxidant capacity of juvenile Pacific geoduck Panopea generosa, suggesting that transcriptional shifts may drive phenotypic modifications post-priming. To this end, juvenile clams were sampled and TagSeq gene expression data analyzed after 1) a 110-day acclimation under ambient (921 μatm, naïve) and moderately-elevated pCO₂ (2870 μatm, pre-exposed); then following 2) a second 7-day exposure to three pCO₂ treatments (ambient: 754 μatm; moderately-elevated: 2750 μatm; severely-elevated: 4940 μatm), a 7-day return to ambient pCO₂ , and a third 7-day exposure to two pCO₂ treatments (ambient: 967 μatm; moderately-elevated: 3030 μatm). Pre-exposed geoducks frontloaded genes for stress and apoptosis/innate immune response, homeostatic processes, protein degradation, and transcriptional modifiers. Pre-exposed geoducks were also responsive to subsequent encounters, with gene sets enriched for mitochondrial recycling and immune defense under elevated pCO₂ and energy metabolism and biosynthesis under ambient recovery. In contrast, gene sets with higher expression in naïve clams were enriched for fatty-acid degradation and glutathione components, suggesting naïve clams could be depleting endogenous fuels, with unsustainable energetic requirements if changes in carbonate chemistry persist. Collectively, our transcriptomic data indicates pCO₂ priming during post-larval periods could, via gene expression regulation, enhance robustness in bivalves to environmental change. Such priming approaches may be beneficial for aquaculture, as seafood demand intensifies concurrent with increasing climate change in marine systems.

Differential DNA methylation in Pacific oyster reproductive tissue in response to ocean acidification

BACKGROUND

There is a need to investigate mechanisms of phenotypic plasticity in marine invertebrates as negative effects of climate change, like ocean acidification, are experienced by coastal ecosystems. Environmentally-induced changes to the methylome may regulate gene expression, but methylome responses can be species- and tissue-specific. Tissue-specificity has implications for gonad tissue, as gonad-specific methylation patterns may be inherited by offspring. We used the Pacific oyster (Crassostrea gigas) - a model for understanding pH impacts on bivalve molecular physiology due to its genomic resources and importance in global aquaculture- to assess how low pH could impact the gonad methylome. Oysters were exposed to either low pH (7.31 ± 0.02) or ambient pH (7.82 ± 0.02) conditions for 7 weeks. Whole genome bisulfite sequencing was used to identify methylated regions in female oyster gonad samples. C- > T single nucleotide polymorphisms were identified and removed to ensure accurate methylation characterization.

RESULTS

Analysis of gonad methylomes revealed a total of 1284 differentially methylated loci (DML) found primarily in genes, with several genes containing multiple DML. Gene ontologies for genes containing DML were involved in development and stress response, suggesting methylation may promote gonad growth homeostasis in low pH conditions. Additionally, several of these genes were associated with cytoskeletal structure regulation, metabolism, and protein ubiquitination - commonly-observed responses to ocean acidification. Comparison of these DML with other Crassostrea spp. exposed to ocean acidification demonstrates that similar pathways, but not identical genes, are impacted by methylation.

CONCLUSIONS

Our work suggests DNA methylation may have a regulatory role in gonad and larval development, which would shape adult and offspring responses to low pH stress. Combined with existing molluscan methylome research, our work further supports the need for tissue- and species-specific studies to understand the potential regulatory role of DNA methylation.

DNA Methylation Analyses Unveil a Regulatory Landscape in the Formation of Nacre Color in Pearl Oyster Pinctada fucata martensii

Pearl color is regulated by genetics, biological pigments, and organic matrices and an important factor that influences the pearl economic value. The epigenetic regulation mechanism underlying pearl pigmentation remains poorly understood. In this study, we collected the mantle pallial (MP) and mantle central (MC) of the golden-lipped strain, and MP of the silver-lipped strain of pearl oyster Pinctada fucata martensii. The whole-genome bisulfite sequencing (WGBS) technology was employed to investigate the possible implication of epigenetic factors regulating nacre color variation. Our results revealed approximately 2.5% of the cytosines in the genome of the P. fucata martensii were methylated, with the CG methylation type was in most abundance. Overall, we identified 12, 621 differentially methylated regions (DMRs) corresponding to 3,471 DMR-associated genes (DMGs) between the two comparison groups. These DMGs were principally enriched into KEGG metabolic pathways including ABC transporters, Terpenoid backbone biosynthesis, and fatty acid degradation. In addition, integrating information about DMGs, DEGs, and function annotation indicated eight genes LDLR, NinaB, RDH, CYP, FADS, fn3, PU-1, KRMP as the candidate genes related to pigmentation of nacre color. A further study proved that the pigment in nacre is violaxanthin. The results of our study provide the support that there is an association between nacre color formation and DNA methylation profiles and will help to reveal the epigenetic regulation of nacre pigmentation formation in pearl oyster P. fucata martensii.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Ecotoxicological epigenetics in invertebrates: Emerging tool for the evaluation of present and past pollution burden

The effect of environmental pollution on epigenetic changes and their heredity in affected organisms is of major concern as such changes can play a significant role in adaptation to changing environmental conditions. Changes of epigenetic marks including DNA methylation, histone modifications, and non-coding RNA's can induce changes in gene transcription leading to physiological long-term changes or even transgenerational inheritance. Such mechanisms have until recently been scarcely studied in invertebrate organisms, mainly focusing on model species including Caenorhabditis elegans and Daphnia magna. However, more data are becoming available, particularly focused on DNA methylation changes caused by anthropogenic pollutants in a wide range of invertebrates. This review examines the literature from field and laboratory studies utilising invertebrate species exposed to environmental pollutants and their effect on DNA methylation. Possible mechanisms of epigenetic modifications and their role on physiology and adaptation as well as the incidence of intergenerational and transgenerational inheritance are discussed. Furthermore, critical research challenges are defined and the way forward is proposed. Future studies should focus on the use of next generation sequencing tools to define invertebrate methylomes under environmental stress in higher resolution, those data should further be linked to gene expression patterns and phenotypes and detailed studies focusing on transgenerational effects are encouraged. Moreover, studies of other epigenetic mechanisms in various invertebrate species, apart from DNA methylation would provide better understanding of interconnected cross-talk between epigenetic marks. Taken together incorporating epigenetic studies in ecotoxicology context presents a promising tool for development of sensitive biomarkers for environmental stress assessment.

Invertebrate methylomes provide insight into mechanisms of environmental tolerance and reveal methodological biases

There is a growing focus on the role of DNA methylation in the ability of marine invertebrates to rapidly respond to changing environmental factors and anthropogenic impacts. However, genome-wide DNA methylation studies in nonmodel organisms are currently hampered by a limited understanding of methodological biases. Here, we compare three methods for quantifying DNA methylation at single base-pair resolution-whole genome bisulfite sequencing (WGBS), reduced representation bisulfite sequencing (RRBS), and methyl-CpG binding domain bisulfite sequencing (MBDBS)-using multiple individuals from two reef-building coral species with contrasting environmental sensitivity. All methods reveal substantially greater methylation in Montipora capitata (11.4%) than the more sensitive Pocillopora acuta (2.9%). The majority of CpG methylation in both species occurs in gene bodies and flanking regions. In both species, MBDBS has the greatest capacity for detecting CpGs in coding regions at our sequencing depth, but MBDBS may be influenced by intrasample methylation heterogeneity. RRBS yields robust information for specific loci albeit without enrichment of any particular genome feature and with significantly reduced genome coverage. Relative genome size strongly influences the number and location of CpGs detected by each method when sequencing depth is limited, illuminating nuances in cross-species comparisons. As genome-wide methylation differences, supported by data across bisulfite sequencing methods, may contribute to environmental sensitivity phenotypes in critical marine invertebrate taxa, these data provide a genomic resource for investigating the functional role of DNA methylation in environmental tolerance.

Symbiont shuffling induces differential DNA methylation responses to thermal stress in the coral Montastraea cavernosa

Algal symbiont shuffling in favour of more thermotolerant species has been shown to enhance coral resistance to heat-stress. Yet, the mechanistic underpinnings and long-term implications of these changes are poorly understood. This work studied the modifications in coral DNA methylation, an epigenetic mechanism involved in coral acclimatization, in response to symbiont manipulation and subsequent heat stress exposure. Symbiont composition was manipulated in the great star coral Montastraea cavernosa through controlled thermal bleaching and recovery, producing paired ramets of three genets dominated by either their native symbionts (genus Cladocopium) or the thermotolerant species (Durusdinium trenchi). Single-base genome-wide analyses showed significant modifications in DNA methylation concentrated in intergenic regions, introns and transposable elements. Remarkably, DNA methylation changes in response to heat stress were dependent on the dominant symbiont, with twice as many differentially methylated regions found in heat-stressed corals hosting different symbionts (Cladocopium vs. D. trenchii) compared to all other comparisons. Interestingly, while differential gene body methylation was not correlated with gene expression, an enrichment in differentially methylated regions was evident in repetitive genome regions. Overall, these results suggest that changes in algal symbionts favouring heat tolerant associations are accompanied by changes in DNA methylation in the coral host. The implications of these results for coral adaptation, along with future avenues of research based on current knowledge gaps, are discussed in the present work.

Genome-wide DNA methylation profile changes associated with shell colouration in the Yesso scallop (Patinopecten yessoensis) as measured by whole-genome bisulfite sequencing

BACKGROUND

Mollusca, a phylum of highly rich species, possess vivid shell colours, but the underlying molecular mechanism remains to be elucidated. DNA methylation, one of the most common epigenetic modifications in eukaryotes, is believed to play a vital role in various biological processes. However, analysis of the effects of DNA methylation on shell colouration has rarely been performed in molluscs, limiting the current knowledge of the molecular mechanism of shell colour formation.

RESULTS

In the present study, to reveal the role of epigenetic regulation in shell colouration, WGBS, the "gold standard" of DNA methylation analysis, was first performed on the mantle tissues of Yesso scallops (Patinopecten yessoensis) with different shell colours (brown and white), and DNA methylomes at single-base resolution were generated. About 3% of cytosines were methylated in the genome of the Yesso scallop. A slight increase in ^mCG percentage and methylation level was found in brown scallops. Sequence preference of nearby methylated cytosines differed between high and low methylation level sites and between the brown- and white-shelled scallops. DNA methylation levels varied among the different genomic regions; all the detected regions in the brown group exhibited higher methylation levels than the white group. A total of 41,175 DMRs (differentially methylated regions) were detected between brown and white scallops. GO functions and pathways associated with the biosynthesis of melanin and porphyrins were significantly enriched for DMRs, among which several key shell colour-related genes were identified. Further, different correlations between mRNA expression levels and DNA methylation status were found in these genes, suggesting that DNA methylation regulates shell colouration in the Yesso scallop.

CONCLUSIONS

This study provides genome-wide DNA methylation landscapes of Yesso scallops with different shell colours, offering new insights into the epigenetic regulatory mechanism underlying shell colour.

Physiological Variation in Response to Vibrio and Hypoxia by Aquacultured Eastern Oysters in the Southeastern United States

Eastern oysters (Crassostrea virginica) have long been recognized as model organisms of extreme environmental tolerance, showing resilience to variation in temperature, salinity, hypoxia and microbial pathogens. These phenotypic responses, however, show variability between geographic locations or habitats (e.g., tidal). Physiological, morphological and genetic differences occur in populations throughout a species' geographical range, which may have been shaped by regional abiotic and biotic variations. Few studies of C. virginica have explored the combined factors of physiological mechanisms of divergent phenotypes between locations and the genetic relationships of individuals between these locations. To characterize genetic relationships of four locations with aquacultured oysters along the North Carolina and Virginia coast, we sequenced a portion of cytochrome oxidase subunit I (COI) that revealed significant variation in haplotype distribution between locations. We then measured mitochondrial physiology and expression of the innate immunity response of hemocytes to lab acclimation and combined stress conditions to compare basal expression and stress response in oysters between these locations. For stress sensing genes, toll-like receptors had the strongest location-specific response to hypoxia and Vibrio, whereas mannose receptor and a stress-receptor were specific to hypoxia and bacteria, respectively. The expression of stress response genes also showed location-specific and stressor-specific changes in expression, particularly for big defensin and the complement gene Cq3. Our results further suggested that genetic similarity of oysters from different locations was not clearly related to physiological and molecular responses. These results are informative for understanding the range of physiological plasticity for stress responses in this commercially important oyster species. They also have implications in the oyster farming industry as well as conservation efforts to restore endangered native oyster beds.

Evolutionary History of DNA Methylation Related Genes in Bivalvia: New Insights From Mytilus galloprovincialis

DNA methylation is an essential epigenetic mechanism influencing gene expression in all organisms. In metazoans, the pattern of DNA methylation changes during embryogenesis and adult life. Consequently, differentiated cells develop a stable and unique DNA methylation pattern that finely regulates mRNA transcription during development and determines tissue-specific gene expression. Currently, DNA methylation remains poorly investigated in mollusks and completely unexplored in Mytilus galloprovincialis. To shed light on this process in this ecologically and economically important bivalve, we screened its genome, detecting sequences homologous to DNA methyltransferases (DNMTs), methyl-CpG-binding domain (MBD) proteins and Ten-eleven translocation methylcytosine dioxygenase (TET) previously described in other organisms. We characterized the gene architecture and protein domains of the mussel sequences and studied their phylogenetic relationships with the ortholog sequences from other bivalve species. We then comparatively investigated their expression levels across different adult tissues in mussel and other bivalves, using previously published transcriptome datasets. This study provides the first insights on DNA methylation regulators in M. galloprovincialis, which may provide fundamental information to better understand the complex role played by this mechanism in regulating genome activity in bivalves.

Population differences in Chinook salmon (Oncorhynchus tshawytscha) DNA methylation: Genetic drift and environmental factors

Local adaptation and phenotypic differences among populations have been reported in many species, though most studies focus on either neutral or adaptive genetic differentiation. With the discovery of DNA methylation, questions have arisen about its contribution to individual variation in and among natural populations. Previous studies have identified differences in methylation among populations of organisms, although most to date have been in plants and model animal species. Here we obtained eyed eggs from eight populations of Chinook salmon (Oncorhynchus tshawytscha) and assayed DNA methylation at 23 genes involved in development, immune function, stress response, and metabolism using a gene-targeted PCR-based assay for next-generation sequencing. Evidence for population differences in methylation was found at eight out of 23 gene loci after controlling for developmental timing in each individual. However, we found no correlation between freshwater environmental parameters and methylation variation among populations at those eight genes. A weak correlation was identified between pairwise DNA methylation dissimilarity among populations and pairwise F ST based on 15 microsatellite loci, indicating weak effects of genetic drift or geographic distance on methylation. The weak correlation was primarily driven by two genes, GTIIBS and Nkef. However, single-gene Mantel tests comparing methylation and pairwise F ST were not significant after Bonferroni correction. Thus, population differences in DNA methylation are more likely related to unmeasured oceanic environmental conditions, local adaptation, and/or genetic drift. DNA methylation is an additional mechanism that contributes to among population variation, with potential influences on organism phenotype, adaptive potential, and population resilience.

Inference of DNA methylation patterns in molluscs

Mollusca are the second largest and arguably most diverse phylum of the animal kingdom. This is in sharp contrast to our very limited knowledge concerning epigenetic mechanisms including DNA methylation in this invertebrate group. Here, we inferred DNA methylation patterns by analysing the normalized dinucleotide CG content in protein-coding sequences and identified DNA methyltransferases (DNMT1 and 3) in published transcriptomes and genomes of 140 species across all eight classes of molluscs. Given the evolutionary age and morphological diversity of molluscs, we expected to find evidence for diverse methylation patterns. Our inferences suggest that molluscs possess substantial levels of DNA methylation in gene bodies as a rule. Yet, we found deviations from this general picture with regard to (i) the CpG observed/expected distributions indicating a reduction in DNA methylation in certain groups and (ii) the completeness of the DNMT toolkit. Reductions were evident in Caudofoveata, Solenogastres, Polyplacophora, Monoplacophora, as well as Scaphopoda. Heterobranchia and Oegopsida were remarkable as they lacked DNMT3, usually responsible for de novo methylation, yet showed signs of DNA methylation. Our survey may serve as guidance for direct empirical analyses of DNA methylation in molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Environmentally Driven Color Variation in the Pearl Oyster Pinctada margaritifera var. cumingii (Linnaeus, 1758) Is Associated With Differential Methylation of CpGs in Pigment- and Biomineralization-Related Genes

Today, it is common knowledge that environmental factors can change the color of many animals. Studies have shown that the molecular mechanisms underlying such modifications could involve epigenetic factors. Since 2013, the pearl oyster Pinctada margaritifera var. cumingii has become a biological model for questions on color expression and variation in Mollusca. A previous study reported color plasticity in response to water depth variation, specifically a general darkening of the nacre color at greater depth. However, the molecular mechanisms behind this plasticity are still unknown. In this paper, we investigate the possible implication of epigenetic factors controlling shell color variation through a depth variation experiment associated with a DNA methylation study performed at the whole genome level with a constant genetic background. Our results revealed six genes presenting differentially methylated CpGs in response to the environmental change, among which four are linked to pigmentation processes or regulations (GART, ABCC1, MAPKAP1, and GRL101), especially those leading to darker phenotypes. Interestingly, the genes perlucin and MGAT1, both involved in the biomineralization process (deposition of aragonite and calcite crystals), also showed differential methylation, suggesting that a possible difference in the physical/spatial organization of the crystals could cause darkening (iridescence or transparency modification of the biomineral). These findings are of great interest for the pearl production industry, since wholly black pearls and their opposite, the palest pearls, command a higher value on several markets. They also open the route of epigenetic improvement as a new means for pearl production improvement.

Genotoxic and epigenetic effects of diuron in the Pacific oyster: in vitro evidence of interaction between DNA damage and DNA methylation

Recently, research has contributed to better knowledge on the occurrence of pesticides in coastal water by identifying frequently detected substances, their concentration range and their acute and chronic toxicity for organisms. Pesticide pollution is of particular concern in France due to important agricultural activities and presence of several exoreic catchment areas that vehicle pesticides up to coastal waters, impacting non-target marine species. Several ecotoxicology questions remain to be addressed concerning the long-term effects of chronic pesticide exposure and the mechanisms involved in adaptation to chemical stress. In the present study, we brought new insights on the genetic and epigenetic effects of the herbicide diuron in oyster genitors. During gametogenesis, we exposed Crassostrea gigas to environmentally realistic herbicide concentrations (0.2-0.3 μg L^-1 during two 7-day periods at half-course and end of gametogenesis). Diuron exposure was shown to decrease global DNA methylation and total methyltransferase activity in whole oyster tissue; this is consistent with the previous observation of a significant decrease in DNMT1 gene expression. Diuron effect seemed to be tissue-specific; hypermethylation was detected in the digestive gland, whereas diuron exposure had no effect on gill and gonad tissue. The genotoxicity of diuron was confirmed by the detection of one adduct in gonad DNA. By using in vitro approaches and human DNMT1 (DNMT1 has not been purified yet in bivalves), the presence of DNA lesions (adduct, 8-oxodGuo) was shown to interfere with DNMT1 activity, indicating a complex interaction between DNA damage and DNA methylation. Based on our results, we propose mechanisms to explain the effect of diuron exposure on DNA methylation, a widespread epigenetic mark.

Benchmarking DNA methylation assays in a reef-building coral

Interrogation of chromatin modifications, such as DNA methylation, has the potential to improve forecasting and conservation of marine ecosystems. The standard method for assaying DNA methylation (whole genome bisulphite sequencing), however, is currently too costly to apply at the scales required for ecological research. Here, we evaluate different methods for measuring DNA methylation for ecological epigenetics. We compare whole genome bisulphite sequencing (WGBS) with methylated CpG binding domain sequencing (MBD-seq), and a modified version of MethylRAD we term methylation-dependent restriction site-associated DNA sequencing (mdRAD). We evaluate these three assays in measuring variation in methylation across the genome, between genotypes, and between polyp types in the reef-building coral Acropora millepora. We find that all three assays measure absolute methylation levels similarly for gene bodies (gbM), as well as exons and 1 Kb windows with a minimum Pearson correlation 0.66. Differential gbM estimates were less correlated, but still concurrent across assays. We conclude that MBD-seq and mdRAD are reliable and cost-effective alternatives to WGBS. The considerably lower sequencing effort required for mdRAD to produce comparable methylation estimates makes it particularly useful for ecological epigenetics.

Genomic Landscape of Mutational Biases in the Pacific Oyster Crassostrea gigas

Mutation is a driving force of evolution that has been shaped by natural selection and is universally biased. Previous studies determined genome-wide mutational patterns for several species and investigated the heterogeneity of mutational patterns at fine-scale levels. However, little evidence of the heterogeneity of mutation rates over large genomic regions was shown. Hence, the mutational patterns of different large-scale genomic regions and their association with selective pressures still need to be explored. As the second most species-rich animal phylum, little is known about the mutational patterns in Mollusca, especially oysters. In this study, the mutational bias patterns are characterized by using whole-genome resequencing data in the Crassostrea gigas genome. I studied the genome-wide relative rates of the pair mutations and found that the predominant mutation is GC -> AT, irrespective of the genomic regions. This analysis reveals that mutational biases were associated with gene expression levels across the C. gigas genome. Genes with higher expression levels and breadth expression patterns, longer coding length, and more exon numbers had relatively higher GC -> AT rates. I also found that genes with larger dN/dS values had relatively higher GC -> AT rates. This work represents the first comprehensive research on the mutational biases in Mollusca species. Here, I comprehensively investigated the relationships between mutational biases with some intrinsic genetic factors and evolutionary indicators and proposed that selective pressures are important forces shaping the mutational biases across the C. gigas genome.

Insights into the DNA methylation of sea cucumber Apostichopus japonicus in response to skin ulceration syndrome infection

DNA methylation is an important epigenetic modification that regulates gene expression in many biological processes, including immune response. In this study, whole-genome bisulfite sequencing (WGBS) was carried out on healthy body wall (HB) and skin ulceration syndrome (SUS) infected body wall (SFB) to gain insights into the epigenetic regulatory mechanism in sea cucumber Apostichopus japonicus. After comparison, a total of 116,522 differentially methylated regions (DMRs) were obtained including 67,269 hyper-methylated and 49,253 hypo-methylated DMRs (p < 0.05, FDR < 0.001). GO enrichment analysis indicated that regulation of DNA-templated transcription (GO: 0006355), where DNA methylation occurred, was the most significant term in the biology process. The integration of methylome and transcriptome analysis revealed that 10,499 DMRs were negatively correlated with 496 differentially expressed genes (DEGs). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that these DEGs were enriched in the phosphoinositide 3-kinase-protein kinase B (PI3K/Akt)/mammalian target of rapamycin (mTOR) signaling pathway. Interestingly, two serine/threonine-protein kinases, nemo-like kinase (NLK) and mTOR, were highlighted after functional analysis. The variations of methylation in these two genes were associated with SUS infection and immune regulation. They regulated gene expression at different levels and showed interaction during response process. The validation of methylation sites showed high consistency between pyrosequencing and WGBS. WGBS analysis not only revealed the changes of DNA methylation, but also presented important information about the regulation of key genes after SUS infection in A. japonicus.

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.

Genome-wide DNA Methylation Analysis of Mantle Edge and Mantle Central from Pearl Oyster Pinctada fucata martensii

DNA methylation is a type of epigenetic modification that alters gene expression without changing the DNA sequence and mediates some cases of phenotypic plasticity. In this study, we identified six DNA methyltransferase (DNMT) genes and two methyl-CpG binding domain protein2 (MBD2) gene from Pinctada fucata martensii. We also analyzed the genome-wide DNA methylation levels of mantle edge (ME) and mantle central (MC) from P. f. martensii via methylated immunoprecipitation sequencing (MeDIP-Seq). Results revealed that both ME and MC had 122 million reads, and had 58,702 and 55,721 peaks, respectively. The obtained methylation patterns of gene elements and repeats showed that the methylation of the protein-coding genes, particularly intron and coding exons (CDSs), was more frequent than that of other genomic elements in the pearl oyster genome. We combined the methylation data with the RNA-seq data of the ME and MC of P. f. martensii and found that promoter, CDS, and intron methylation levels were positively correlated with gene expression levels except the highest gene expression level. We also identified 313 differential methylation genes (DMGs) and annotated 212 of them. These DMGs were significantly enriched in 30 pathways, such as amino acid and protein metabolism, energy metabolism, terpenoid synthesis, and immune-related pathways. This study comprehensively analyzed the methylomes of biomineralization-related tissues and helped enhance our understanding of the regulatory mechanism underlying shell formation.

Population epigenetic divergence exceeds genetic divergence in the Eastern oyster Crassostrea virginica in the Northern Gulf of Mexico

Populations may respond to environmental heterogeneity via evolutionary divergence or phenotypic plasticity. While evolutionary divergence occurs through DNA sequence differences among populations, plastic divergence among populations may be generated by changes in the epigenome. Here, we present the results of a genome-wide comparison of DNA methylation patterns and genetic structure among four populations of Eastern oyster (Crassostrea virginica) in the northern Gulf of Mexico. We used a combination of restriction site-associated DNA sequencing (RADseq) and reduced representation bisulfite sequencing (RRBS) to explore population structure, gene-wide averages of F ST, and DNA methylation differences between oysters inhabiting four estuaries with unique salinity profiles. This approach identified significant population structure despite a moderately low F ST (0.02) across the freshwater boundary of the Mississippi river, a finding that may reflect recent efforts to restore oyster stock populations. Divergence between populations in CpG methylation was greater than for divergence in F ST, likely reflecting environmental effects on DNA methylation patterns. Assessment of CpG methylation patterns across all populations identified that only 26% of methylated DNA was intergenic; and, only 17% of all differentially methylated regions (DMRs) were within these same regions. DMRs within gene bodies between sites were associated with genes known to be involved in DNA damage repair, ion transport, and reproductive timing. Finally, when assessing the correlation between genomic variation and DNA methylation between these populations, we observed population-specific DNA methylation profiles that were not directly associated with single nucleotide polymorphisms or broader gene-body mean F ST trends. Our results suggest that C. virginica may use DNA methylation to generate environmentally responsive plastic phenotypes and that there is more divergence in methylation than divergence in allele frequencies.

DNA methylation is maintained with high fidelity in the honey bee germline and exhibits global non-functional fluctuations during somatic development

BACKGROUND

DNA methylation of active genes, also known as gene body methylation, is found in many animal and plant genomes. Despite this, the transcriptional and developmental role of such methylation remains poorly understood. Here, we explore the dynamic range of DNA methylation in honey bee, a model organism for gene body methylation.

RESULTS

Our data show that CG methylation in gene bodies globally fluctuates during honey bee development. However, these changes cause no gene expression alterations. Intriguingly, despite the global alterations, tissue-specific CG methylation patterns of complete genes or exons are rare, implying robust maintenance of genic methylation during development. Additionally, we show that CG methylation maintenance fluctuates in somatic cells, while reaching maximum fidelity in sperm cells. Finally, unlike universally present CG methylation, we discovered non-CG methylation specifically in bee heads that resembles such methylation in mammalian brain tissue.

CONCLUSIONS

Based on these results, we propose that gene body CG methylation can oscillate during development if it is kept to a level adequate to preserve function. Additionally, our data suggest that heightened non-CG methylation is a conserved regulator of animal nervous systems.

What drives phenotypic divergence among coral clonemates of Acropora palmata?

Evolutionary rescue of populations depends on their ability to produce phenotypic variation that is heritable and adaptive. DNA mutations are the best understood mechanisms to create phenotypic variation, but other, less well-studied mechanisms exist. Marine benthic foundation species provide opportunities to study these mechanisms because many are dominated by isogenic stands produced through asexual reproduction. For example, Caribbean acroporid corals are long lived and reproduce asexually via breakage of branches. Fragmentation is often the dominant mode of local population maintenance. Thus, large genets with many ramets (colonies) are common. Here, we observed phenotypic variation in stress responses within genets following the coral bleaching events in 2014 and 2015 caused by high water temperatures. This was not due to genetic variation in their symbiotic dinoflagellates (Symbiodinium "fitti") because each genet of this coral species typically harbours a single strain of S. "fitti". Characterization of the microbiome via 16S tag sequencing correlated the abundance of only two microbiome members (Tepidiphilus, Endozoicomonas) with a bleaching response. Epigenetic changes were significantly correlated with the host's genetic background, the location of the sampled polyps within the colonies (e.g., branch vs. base of colony), and differences in the colonies' condition during the bleaching event. We conclude that long-term microenvironmental differences led to changes in the way the ramets methylated their genomes, contributing to the differential bleaching response. However, most of the variation in differential bleaching response among clonemates of Acropora palmata remains unexplained. This research provides novel data and hypotheses to help understand intragenet variability in stress phenotypes of sessile marine species.

Mollusc shells as metagenomic archives: The true treasure is the chest itself

Mollusc shells, beyond the treasure of information inherently conveyed through their morphology and chemical composition also have the capacity to preserve DNA sequences over the long term in their inner structure. This has been clearly demonstrated for the first time in the study published in this issue of Molecular Ecology Resources by Der Sarkissian et al. (). With a methodology specifically dedicated to ancient DNA and solid matrices, the authors were able to successfully extract and amplify DNA from marine shells spanning the last 7,000 years. Furthermore, using metagenomic analyses, they could identify important factors affecting DNA recovery. Using reference genomes and sequences in a targeted approach to assign high-throughput sequencing reads, the authors revealed both the presence of endogenous mollusc DNA and a potent pathogen of Manilla clam. Collectively, the results presented in this study open extremely promising research avenues, from palaeogenomics and evolutionary biology to ecological genomics at population and community levels, as well as the opportunity to fine-tune diagnostic tools for conservation and aquaculture purposes. Last but not least, this study also offers exciting perspectives in epigenomics and the evolution of regulatory processes in the context of adaptation to global change. It can be easily expected that the approach developed by Der Sarkissian et al. () will be pursued and extensively investigated in the near future by the scientific community interested in these issues.

Dynamics of DNA Methylation and DNMT Expression During Gametogenesis and Early Development of Scallop Patinopecten yessoensis

DNA methylation reprograms during gametogenesis and embryo development, which is essential for germ cell specification and genomic imprinting in mammals. Corresponding process remains poorly investigated in molluscs. Here, we examined global DNA methylation level in the gonads of scallop Patinopecten yessoensis during gametogenesis and in embryos/larvae at different stages. DNA methylation level fluctuates during gametogenesis and early development, peaking at proliferative stage of ovary, growing stage of testis, and in blastulae. To understand the mechanisms underlying these changes, we conducted genome-wide characterization of DNMT family and investigated their expression profiles based on transcriptomes and in situ hybridization. Three genes were identified, namely PyDNMT1, PyDNMT2, and PyDNMT3. Expression of PyDnmt3 agrees with DNA methylation level during oogenesis and early development, suggesting PyDNMT3 may participate in de novo DNA methylation that occurs mainly at proliferative stage of ovary and testis, and in blastulae and gastrulae. PyDnmt1 expression is positively correlated with DNA methylation level during spermatogenesis, and is higher at maturation stage of ovary and in 2-8 cell embryos than other stages, implying possible involvement of PyDNMT1 in DNA methylation maintenance during meiosis and embryonic development. This study will facilitate better understanding of the developmental epigenetic reprogramming in bivalve molluscs.

Marine Environmental Epigenetics

Marine organisms' persistence hinges on the capacity for acclimatization and adaptation to the myriad of interacting environmental stressors associated with global climate change. In this context, epigenetics-mechanisms that facilitate phenotypic variation through genotype-environment interactions-are of great interest ecologically and evolutionarily. Our comprehensive review of marine environmental epigenetics guides our recommendations of four key areas for future research: the dynamics of wash-in and wash-out of epigenetic effects, the mechanistic understanding of the interplay of different epigenetic marks and the interaction with the microbiome, the capacity for and mechanisms of transgenerational epigenetic inheritance, and the evolutionary implications of the interaction of genetic and epigenetic features. Emerging insights in marine environmental epigenetics can be applied to critical issues such as aquaculture, biomonitoring, and biological invasions, thereby improving our ability to explain and predict the responses of marine taxa to global climate change.

Coral epigenetic responses to nutrient stress: Histone H2A.X phosphorylation dynamics and DNA methylation in the staghorn coral Acropora cervicornis

Nutrient pollution and thermal stress constitute two of the main drivers of global change in the coastal oceans. While different studies have addressed the physiological effects and ecological consequences of these stressors in corals, the role of acquired modifications in the coral epigenome during acclimatory and adaptive responses remains unknown. The present work aims to address that gap by monitoring two types of epigenetic mechanisms, namely histone modifications and DNA methylation, during a 7-week-long experiment in which staghorn coral fragments (Acropora cervicornis) were exposed to nutrient stress (nitrogen, nitrogen + phosphorus) in the presence of thermal stress. The major conclusion of this experiment can be summarized by two main results: First, coral holobiont responses to the combined effects of nutrient enrichment and thermal stress involve the post-translational phosphorylation of the histone variant H2A.X (involved in responses to DNA damage), as well as nonsignificant modifications in DNA methylation trends. Second, the reduction in H2A.X phosphorylation (and the subsequent potential impairment of DNA repair mechanisms) observed after prolonged coral exposure to nitrogen enrichment and thermal stress is consistent with the symbiont-driven phosphorus limitation previously observed in corals subject to nitrogen enrichment. The alteration of this epigenetic mechanism could help to explain the synergistic effects of nutrient imbalance and thermal stress on coral fitness (i.e., increased bleaching and mortality) while supporting the positive effect of phosphorus addition to improving coral resilience to thermal stress. Overall, this work provides new insights into the role of epigenetic mechanisms during coral responses to global change, discussing future research directions and the potential benefits for improving restoration, management and conservation of coral reef ecosystems worldwide.

Natural variability and modulation by environmental stressors of global genomic cytosine methylation levels in a freshwater crustacean, Gammarus fossarum

To improve the assessment of aquatic organism responses to environmental stressors, there is an interest in studying epigenetic marks in addition to other validated biomarkers. Indeed, the epigenetic marks may be influenced by the surrounding environment. Non-model invertebrates such as gammarids are sentinel organisms representative of the diversity of natural stream communities. Despite their ecologically relevance, the epigenetic responses have been to date poorly documented in these species. The present study explores the measurement of the global cytosine methylation level in the genome of the freshwater crustacean Gammarus fossarum. In a first step, natural variability of global cytosine methylation level (basal level) was assessed by studying the effect of sex, age and sampling site of organisms. Results showed a significant effect of age and sampling site. In a second step, effects of water temperature and food starvation were studied. For both factors, a hypermethylation was observed after 1 month of exposure. In a third step, gammarids were exposed to a range of environmentally relevant cadmium concentrations (0.05-5 μg/L) in order to assess the effect of a chemical stress. Whatever the cadmium concentration used, a significant hypomethylation was observed after 14 days followed by a trend for hypermethylation after 1 month of exposure. These results are the first ones dealing with the 5C-methylation status in gammarids. The results constitute potential markers of environmental stresses in relevant sentinel species widely used in ecotoxicological studies.

Pharmacoepigenomic Interventions as Novel Potential Treatments for Alzheimer's and Parkinson's Diseases

Cerebrovascular and neurodegenerative disorders affect one billion people around the world and result from a combination of genomic, epigenomic, metabolic, and environmental factors. Diagnosis at late stages of disease progression, limited knowledge of gene biomarkers and molecular mechanisms of the pathology, and conventional compounds based on symptomatic rather than mechanistic features, determine the lack of success of current treatments, including current FDA-approved conventional drugs. The epigenetic approach opens new avenues for the detection of early presymptomatic pathological events that would allow the implementation of novel strategies in order to stop or delay the pathological process. The reversibility and potential restoring of epigenetic aberrations along with their potential use as targets for pharmacological and dietary interventions sited the use of epidrugs as potential novel candidates for successful treatments of multifactorial disorders involving neurodegeneration. This manuscript includes a description of the most relevant epigenetic mechanisms involved in the most prevalent neurodegenerative disorders worldwide, as well as the main potential epigenetic-based compounds under investigation for treatment of those disorders and their limitations.

Insights into Epigenome Evolution from Animal and Plant Methylomes

Evolutionary studies of DNA methylation offer insights into the mechanisms governing the variation of genomic DNA methylation across different species. Comparisons of gross levels of DNA methylation between distantly related species indicate that the size of the genome and the level of genomic DNA methylation are positively correlated. In plant genomes, this can be reliably explained by the genomic contents of repetitive sequences. In animal genomes, the role of repetitive sequences on genomic DNA methylation is less clear. On a shorter timescale, population-level comparisons demonstrate that genetic variation can explain the observed variability of DNA methylation to some degree. The amount of DNA methylation variation that has been attributed to genetic variation in the human population studies so far is substantially lower than that from Arabidopsis population studies, but this disparity might reflect the differences in the computational and experimental techniques used. The effect of genetic variation on DNA methylation has been directly examined in mammalian systems, revealing several causative factors that govern DNA methylation. On the other hand, studies from Arabidopsis have furthered our understanding of spontaneous mutations of DNA methylation, termed “epimutations.” Arabidopsis has an extremely high rate of spontaneous epimutations, which may play a major role in shaping the global DNA methylation landscape in this genome. Key missing information includes the frequencies of spontaneous epimutations in other lineages, in particular animal genomes, and how population-level variation of DNA methylation leads to species-level differences.

Epigenetic considerations in aquaculture

Epigenetics has attracted considerable attention with respect to its potential value in many areas of agricultural production, particularly under conditions where the environment can be manipulated or natural variation exists. Here we introduce key concepts and definitions of epigenetic mechanisms, including DNA methylation, histone modifications and non-coding RNA, review the current understanding of epigenetics in both fish and shellfish, and propose key areas of aquaculture where epigenetics could be applied. The first key area is environmental manipulation, where the intention is to induce an ‘epigenetic memory’ either within or between generations to produce a desired phenotype. The second key area is epigenetic selection, which, alone or combined with genetic selection, may increase the reliability of producing animals with desired phenotypes. Based on aspects of life history and husbandry practices in aquaculture species, the application of epigenetic knowledge could significantly affect the productivity and sustainability of aquaculture practices. Conversely, clarifying the role of epigenetic mechanisms in aquaculture species may upend traditional assumptions about selection practices. Ultimately, there are still many unanswered questions regarding how epigenetic mechanisms might be leveraged in aquaculture.

The Role of Epigenomics in Aquatic Toxicology

Over the past decade, the field of molecular biology has rapidly incorporated epigenetic studies to evaluate organism-environment interactions that can result in chronic effects. Such responses arise from early life stage stress, the utilization of genetic information over an individual's life time, and transgenerational inheritance. Knowledge of epigenetic mechanisms provides the potential for a comprehensive evaluation of multigenerational and heritable effects from environmental stressors, such as contaminants. Focused studies have provided a greater understanding of how many responses to environmental stressors are driven by epigenetic modifiers. We discuss the promise of epigenetics and suggest future research directions within the field of aquatic toxicology, with a particular focus on the potential for identifying key heritable marks with consequential impacts at the organism and population levels. Environ Toxicol Chem 2017;36:2565-2573. © 2017 SETAC.

The association between DNA methylation and exon expression in the Pacific oyster Crassostrea gigas

Background

DNA methylation is one of the most important epigenetic modifications of eukaryotic genomes and is believed to play integral roles in diverse biological processes. Although DNA methylation has been well studied in mammals, data are limited in invertebrates, particularly Mollusca. The Pacific oyster Crassostrea gigas is an emerging genetic model for functional analysis of DNA methylation in Mollusca. Recent studies have shown that there is a positive association between methylation status and gene expression in C. gigas; however, whether this association exists at the exon level remains to be determined.

Results

In this study, we characterized the genome-wide methylation pattern across two different tissues of C. gigas and found that methylated genes are expressed in more tissues and development stages than unmethylated genes. Furthermore, we found that different types of exons had different methylation levels, with the lowest methylation levels in the first exons, followed by the last exons, and the internal exons. We found that the exons included in the gene transcript contained significantly higher DNA methylation levels than skipped exons. We observed that the DNA methylation levels increased slowly after the start sites and end sites of exons seperately, and then decreased quickly towards the middle sites of exons. We also found that methylated exons were significantly longer than unmethylated exons.

Conclusion

This study constitutes the first genome-wide analysis to show an association between exon-level DNA methylation and mRNA expression in the oyster. Our findings suggest that exon-level DNA methylation may play a role in the construction of alternative splicing by positively influencing exon inclusion during transcription.

Molecular and Biochemical Methods Useful for the Epigenetic Characterization of Chromatin-Associated Proteins in Bivalve Molluscs

Bivalve molluscs constitute a ubiquitous taxonomic group playing key functions in virtually all ecosystems, and encompassing critical commercial relevance. Along with a sessile and filter-feeding lifestyle in most cases, these characteristics make bivalves model sentinel organisms routinely used for environmental monitoring studies in aquatic habitats. The study of epigenetic mechanisms linking environmental exposure and specific physiological responses (i.e., environmental epigenetics) stands out as a very innovative monitoring strategy, given the role of epigenetic modifications in acclimatization and adaptation. Furthermore, the heritable nature of many of those modifications constitutes a very promising avenue to explore the applicability of epigenetic conditioning and selection in management and restoration strategies. Chromatin provides a framework for the study of environmental epigenetic responses. Unfortunately, chromatin and epigenetic information are very limited in most non-traditional model organisms and even completely lacking in most environmentally and ecologically relevant organisms. The present work aims to provide a comprehensive and reproducible experimental workflow for the study of bivalve chromatin. First, a series of guidelines for the molecular isolation of genes encoding chromatin-associated proteins is provided, including information on primers suitable for conventional PCR, Rapid Amplification of cDNA Ends (RACE), genome walking and quantitative PCR (qPCR) experiments. This section is followed by the description of methods specifically developed for the analysis of histone and SNBP proteins in different bivalve tissues, including protein extraction, purification, separation and immunodetection. Lastly, information about available antibodies, their specificity and performance is also provided. The tools and protocols described here complement current epigenetic analyses (usually limited to DNA methylation) by incorporating the study of structural elements modulating chromatin dynamics.

Perkinsus infection is associated with alterations in the level of global DNA methylation of gills and gastrointestinal tract of the oyster Crassostrea gasar

Bivalves are filter feeders that obtain food from seawater that may contain infectious agents, such as the protozoan parasites Perkinsus marinus and P. olseni that are associated with massive mortalities responsible for losses in the aquaculture industry. Despite all physical and chemical barriers, microorganisms cross epithelia and infect host tissues to cause pathologies. Epigenetics mechanisms play important roles in a variety of human processes, from embryonic development to cell differentiation and growth. It is currently emerging as crucial mechanism involved in modulation of host-parasite interactions and pathogenesis, promoting discovery of targets for drug treatment. In bivalves, little is known about epigenetic mechanism in host parasite interactions. The objective of the present study was to evaluate the effect of Perkinsus sp. infections on DNA methylation levels in tissues of Crassostrea gasar oysters. Samples were collected in 2015 and 2016 in the Mamanguape River estuary (PB). Oyster gills were removed and used for Perkinsus sp.

DIAGNOSIS

Gills (G) and gastrointestinal tract (GT), as well as cultured P. marinus trophozoites were preserved in 95% ethanol for DNA extractions. DNA methylation levels were estimated from G and GT tissues of uninfected (n=60) and infected oysters (n=60), and from P. marinus trophozoites, by ELISA assays. Results showed that the mean prevalence of Perkinsus sp. infections was high (87.3%) in 2015 and moderate (59.6%) in 2016. DNA methylation levels of G and GT tissues were significantly lower in infected oyster than in uninfected oysters, suggesting that infections are associated with hypomethylation. Methylation level was significantly higher in G than in GT tissues, indicating a likely tissue-specific mechanism. P. marinus trophozoites showed 33% methylation. This was the first study that confirms alterations of DNA methylation in two tissues of C. gasar oysters in association with Perkinsus sp. infections.

Methylation profile of a satellite DNA constituting the intercalary G+C-rich heterochromatin of the cut trough shell Spisula subtruncata (Bivalvia, Mactridae)

Tandemly repeated DNAs usually constitute significant portions of eukaryotic genomes. In bivalves, however, repetitive DNAs are habitually not widespread. In our search for abundant repetitive DNAs in trough shells, we discovered a novel satellite DNA, SSUsat, which constitutes at least 1.3% of the genome of Spisula subtruncata. As foreseen by the satellite DNA library hypothesis, we confirmed that this satellite DNA is also present in two other Mactridae species, showing a highly conserved nucleotide sequence together with a dramatic diminution in the number of repeats. Predominantly located at the G + C-rich intercalary heterochromatin of S. subtruncata, SSUsat displays several DNA methylation peculiarities. The level of methylation of SSUsat is high (3.38%) in comparison with bivalve standards and triplicates the mean of the S. subtruncata genome (1.13%). Methylation affects not only the cytosines in CpG dinucleotides but also those in CHH and CHG trinucleotides, a feature common in plants but scarce and without any clear known relevance in animals. SSUsat segments enriched in methylated cytosines partly overlap those showing higher sequence conservation. The presence of a chromosome pair showing an accumulation of markedly under-methylated SSUsat monomers additionally indicates that the methylation processes that shape repetitive genome compartments are quite complex.

Dynamics of DNA methylomes underlie oyster development

DNA methylation is a critical epigenetic regulator of development in mammals and social insects, but its significance in development outside these groups is not understood. Here we investigated the genome-wide dynamics of DNA methylation in a mollusc model, the oyster Crassostrea gigas, from the egg to the completion of organogenesis. Large-scale methylation maps reveal that the oyster genome displays a succession of methylated and non methylated regions, which persist throughout development. Differentially methylated regions (DMRs) are strongly regulated during cleavage and metamorphosis. The distribution and levels of methylated DNA within genomic features (exons, introns, promoters, repeats and transposons) show different developmental lansdscapes marked by a strong increase in the methylation of exons against introns after metamorphosis. Kinetics of methylation in gene-bodies correlate to their transcription regulation and to distinct functional gene clusters, and DMRs at cleavage and metamorphosis bear the genes functionally related to these steps, respectively. This study shows that DNA methylome dynamics underlie development through transcription regulation in the oyster, a lophotrochozoan species. To our knowledge, this is the first demonstration of such epigenetic regulation outside vertebrates and ecdysozoan models, bringing new insights into the evolution and the epigenetic regulation of developmental processes.

Elucidating the mechanisms which govern the development of multicellular animals and their evolution is a fundamental task. Epigenetic mechanisms like DNA methylation have recently emerged as critical regulators of mammalian development through the control of genes that determine the identity of cells and the transmission of parental imprints. In invertebrates however, DNA is mostly unmethylated and does not play a role in development except in the peculiar case of social insects. Therefore the significance of DNA methylation in development is thought to be restricted to vertebrates, and thereby considered a recent evolutionary acquisition, and the situation in more distant organisms is unknown. Here we investigated the dynamics of genome-wide DNA methylation patterns in a mollusc, the oyster C. gigas, throughout its development. We found that the dynamics of DNA methylation correspond to the expression dynamics of distinct functional gene clusters that control two critical development steps, cleavage and metamorphosis, and we provide insights into the underlying molecular mechanisms in a non-vertebrate species. These findings challenge the present considerations on the evolution of developmental processes and their epigenetic regulation, and open a new area of research in molecular and developmental biology in invertebrates.

Epigenetic DNA Methylation Mediating Octopus vulgaris Early Development: Effect of Essential Fatty Acids Enriched Diet

The common octopus, Octopus vulgaris, is a good candidate for aquaculture but a sustainable production is still unviable due to an almost total mortality during the paralarvae stage. DNA methylation regulates gene expression in the eukaryotic genome, and has been shown to exhibit plasticity throughout O. vulgaris life cycle, changing profiles from paralarvae to adult stages. This pattern of methylation could be sensitive to small alterations in nutritional and environmental conditions during the species early development, thus impacting on its health, growth and survival. In this sense, a full understanding of the epigenetic mechanisms operating during O. vulgaris development would contribute to optimizing the culture conditions for this species. Paralarvae of O. vulgaris were cultured over 28 days post-hatching (dph) using two different Artemia sp. based diets: control and a long chain polyunsaturated fatty acids (LC-PUFA) enriched diet. The effect of the diets on the paralarvae DNA global methylation was analyzed by Methyl-Sensitive Amplification Polymorphism (MSAP) and global 5-methylcytosine enzyme-linked immunosorbent assay (ELISA) approaches. The analysis of different methylation states over the time revealed a global demethylation phenomena occurring along O. vulgaris early development being directly driven by the age of the paralarvae. A gradual decline in methylated loci (hemimethylated, internal cytosine methylated, and hypermethylated) parallel to a progressive gain in non-methylated (NMT) loci toward the later sampling points was verified regardless of the diet provided and demonstrate a pre-established and well-defined demethylation program during its early development, involving a 20% of the MSAP loci. In addition, a differential behavior between diets was also observed at 20 dph, with a LC-PUFA supplementation effect over the methylation profiles. The present results show significant differences on the paralarvae methylation profiles during its development and a diet effect on these changes. It is characterized by a process of demethylation of the genome at the paralarvae stage and the influence of diet to favor this methylation loss.

Effects of Florida Red Tides on histone variant expression and DNA methylation in the Eastern oyster Crassostrea virginica

Massive algal proliferations known as Harmful Algal Blooms (HABs) represent one of the most important threats to coastal areas. Among them, the so-called Florida Red Tides (FRTs, caused by blooms of the dinoflagellate Karenia brevis and associated brevetoxins) are particularly detrimental in the southeastern U.S., causing high mortality rates and annual losses in excess of $40 million. The ability of marine organisms to cope with environmental stressors (including those produced during HABs) is influenced by genetic and epigenetic mechanisms, the latter resulting in phenotypic changes caused by heritable modifications in gene expression, without involving changes in the genetic (DNA) sequence. Yet, studies examining cause-effect relationships between environmental stressors, specific epigenetic mechanisms and subsequent responses are still lacking. The present work contributes to increase this knowledge by investigating the effects of Florida Red Tides on two types of mechanisms participating in the epigenetic memory of Eastern oysters: histone variants and DNA methylation. For that purpose, a HAB simulation was conducted in laboratory conditions, exposing oysters to increasing concentrations of K. brevis. The obtained results revealed, for the first time, the existence of H2A.X, H2A.Z and macroH2A genes in this organism, encoding histone variants potentially involved in the maintenance of genome integrity during responses to the genotoxic effect of brevetoxins. Additionally, an increase in H2A.X phosphorylation (γH2A.X, a marker of DNA damage) and a decrease in global DNA methylation were observed as the HAB simulation progressed. Overall, the present work provides a basis to better understand how epigenetic mechanisms participate in responses to environmental stress in marine invertebrates, opening new avenues to incorporate environmental epigenetics approaches into management and conservation programs.

The Biomphalaria glabrata DNA methylation machinery displays spatial tissue expression, is differentially active in distinct snail populations and is modulated by interactions with Schistosoma mansoni

BACKGROUND

The debilitating human disease schistosomiasis is caused by infection with schistosome parasites that maintain a complex lifecycle alternating between definitive (human) and intermediate (snail) hosts. While much is known about how the definitive host responds to schistosome infection, there is comparably less information available describing the snail's response to infection.

METHODOLOGY/PRINCIPLE FINDINGS

Here, using information recently revealed by sequencing of the Biomphalaria glabrata intermediate host genome, we provide evidence that the predicted core snail DNA methylation machinery components are associated with both intra-species reproduction processes and inter-species interactions. Firstly, methyl-CpG binding domain protein (Bgmbd2/3) and DNA methyltransferase 1 (Bgdnmt1) genes are transcriptionally enriched in gonadal compared to somatic tissues with 5-azacytidine (5-AzaC) treatment significantly inhibiting oviposition. Secondly, elevated levels of 5-methyl cytosine (5mC), DNA methyltransferase activity and 5mC binding in pigmented hybrid- compared to inbred (NMRI)- B. glabrata populations indicate a role for the snail's DNA methylation machinery in maintaining hybrid vigour or heterosis. Thirdly, locus-specific detection of 5mC by bisulfite (BS)-PCR revealed 5mC within an exonic region of a housekeeping protein-coding gene (Bg14-3-3), supporting previous in silico predictions and whole genome BS-Seq analysis of this species' genome. Finally, we provide preliminary evidence for parasite-mediated host epigenetic reprogramming in the schistosome/snail system, as demonstrated by the increase in Bgdnmt1 and Bgmbd2/3 transcript abundance following Bge (B. glabrata embryonic cell line) exposure to parasite larval transformation products (LTP).

CONCLUSIONS/SIGNIFICANCE

The presence of a functional DNA methylation machinery in B. glabrata as well as the modulation of these gene products in response to schistosome products, suggests a vital role for DNA methylation during snail development/oviposition and parasite interactions. Further deciphering the role of this epigenetic process during Biomphalaria/Schistosoma co-evolutionary biology may reveal key factors associated with disease transmission and, moreover, enable the discovery of novel lifecycle intervention strategies.

Effects of soybean meal on digestive enzymes activity, expression of inflammation-related genes, and chromatin modifications in marine fish (Sparus aurata L.) larvae

The effects of soybean meal (SBM) in early diet of Sparus aurata larvae at two developmental windows were assessed. Prolonged (beyond 14 days post-hatch, dph) feeding with SBM decreased the activity of pancreatic enzymes of larvae. In the absence of SBM these larvae later resumed enzyme activities, but exhibited a significant delay in development. Larvae response to SBM involved up-regulation of extracellular matrix remodeling enzymes and pro-inflammatory cytokines, coupled with a drop in putative intestinal enzymes. Larvae receiving SBM at first feeding appear later to have lower expression of inflammation-related genes, especially those fed SBM until 14 dph. Multivariate analysis confirmed that the duration of the SBM early feeding period drives the physiology of larvae in different directions. Feeding larvae with SBM increased global histone H3 acetylation, whereas upon removal of SBM the process was reverted. A more in deep analysis revealed a dynamic interplay among several reversible histone modifications such as H3K14ac and H3K27m3. Finally, we showed that SBM feeding of larvae results in global hypomethylation that persist after SBM removal. This study is the first demonstrating an effect of diet on marine fish epigenetics. It is concluded that there are limitations for extending SBM feeding of S. aurata larvae beyond 14 dph even under co-feeding with live feed, affecting key physiological processes and normal growth. However, up to 14 dph, SBM does not affect normal development, and produces apparently lasting effects on some key enzymes, genes, and chromatin modifications.