Epigenetic features in the oyster Crassostrea gigas suggestive of functionally relevant promoter DNA methylation in invertebrates

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.

Cadmium nitrate and DNA methylation in gastropods: comparison between ovotestis and hepatopancreas

Dietary ingestion is the main route of exposure to hazardous contaminants in land animals. Cadmium, a high-profile toxic metal, affects living systems at different organismal levels, including major storage organs (liver, kidneys), key organs for species survival (gonads), and epigenetic networks regulating gene expression. 5-methylcytosine (5mC) is the most common and best-characterized epigenetic mark among different modified nucleosides in DNA. This important player in methylation-driven gene expression is impacted by cadmium in sentinel terrestrial vertebrates. However, limited information exists regarding its impact on macroinvertebrates, especially land snails commonly used as (eco)toxicological models. We first investigate the methylomic effects of dietary cadmium given as cadmium nitrate on terrestrial mollusks. Mature specimens of the common brown garden snail, Cornu aspersum, were continuously exposed for four weeks to environmentally-relevant cadmium levels. We determined global genomic DNA methylation in hepatopancreas and ovotestis, as well as changes in the methylation status of CG pairs at the 5' region close to the transcription site of gene encoding the Cd-selective metallothionein (Cd-MT). Weight gain/loss, hypometabolism tendency, and survival rates were also assessed. Although this exposure event did not adversely affect survival, gastropods exposed to the highest Cd dose revealed a significant reduction in body weight and a significant increase in hypometabolic behavior. The hepatopancreas, but not the ovotestis, displayed significant hypermethylation, but only for the aforementioned specimens. We also found that the 5' end of the Cd-MT gene was unmethylated in both organs and its methylation status was insensitive to cadmium exposure. Our results are important since they provide scientists, for the first time, with quantitative data on DNA methylation in gastropod ovotestis and refine our understanding of Cd epigenetic effects on terrestrial mollusks.

Genome-Wide Differential DNA Methylomes Provide Insights into the Infertility of Triploid Oysters

Chromosomal incompatibility and gene expression changes would affect the development of polyploid gonad and gamete formation. The role of epigenetics like DNA methylation in reproductive development is fully demonstrated in diploid animals. The lack of polyploid species and the infertility of polyploid animals, especially the odd ploidy, limit the study of epigenetic regulation mechanism of polyploid reproduction. Fertile and infertile individuals exist in triploid Pacific oyster Crassostrea gigas, which provide an interesting model for studies on the effect of epigenetic regulation on gonadal development. The whole genome single base resolution DNA methylomes in gonads of triploid females α (F-3nα), triploid females β (F-3nβ), triploid males α (M-3nα), triploid hermaphrodite predominantly males (HPM-3n), diploid females (F-2n), and diploid males (M-2n) were generated by using bisulfite-sequencing. The overall DNA methylation profiles in gene regions and transposable regions of fertile and infertile triploid oysters were consistent with those of diploid oysters. The DNA methylation level of CG context decreased in infertile triploid oysters, with more hypomethylated than hypermethylated regions, and the opposite is true in fertile triploid oysters. Genes harbored with differentially methylated regions (DMRs) in infertile triploids were mainly related to the metabolism pathways and the signal pathways. Correlation analysis indicated that the expression of gene transcriptions was generally positively associated with DNA methylation in gene body regions, and DMRs in infertile triploid oysters played significant roles in gonadal development as a possible critical epigenetic regulator of gonadal development gene transcriptional activity. These findings indicate a potential relationship between DNA methylation variability and gene expression plasticity in newly formed polyploidy. As far as we know, this is the first study revealing the epigenetic regulation of gonadal development in invertebrates based on fertile and infertile models, meanwhile providing a new mentality to explore the regulatory mechanisms of infertility in triploids.

The methylome of Biomphalaria glabrata and other mollusks: enduring modification of epigenetic landscape and phenotypic traits by a new DNA methylation inhibitor

BACKGROUND

5-Methylcytosine (5mC) is an important epigenetic mark in eukaryotes. Little information about its role exists for invertebrates. To investigate the contribution of 5mC to phenotypic variation in invertebrates, alteration of methylation patterns needs to be produced. Here, we apply new non-nucleoside DNA methyltransferase inhibitors (DNMTi) to introduce aleatory changes into the methylome of mollusk species.

RESULTS

Flavanone inhibitor Flv1 was efficient in reducing 5mC in the freshwater snails Biomphalaria glabrata and Physa acuta, and to a lesser degree, probably due to lower stability in sea water, in the oyster Crassostrea gigas. Flv1 has no toxic effects and significantly decreased the 5mC level in the treated B. glabrata and in its offspring. Drug treatment triggers significant variation in the shell height in both generations. A reduced representation bisulfite-sequencing method called epiGBS corroborates hypomethylation effect of Flv1 in both B. glabrata generations and identifies seven Differential Methylated Regions (DMR) out of 32 found both in Flv1-exposed snails and its progeny, from which 5 were hypomethylated, demonstrating a multigenerational effect. By targeted bisulfite sequencing, we confirmed hypomethylation in a locus and show that it is associated with reduced gene expression.

CONCLUSIONS

Flv1 is a new and efficient DNMTi that can be used to induce transient and heritable modifications of the epigenetic landscape and phenotypic traits in mollusks, a phylum of the invertebrates in which epigenetics is understudied.

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.

Sequence determinants, function, and evolution of CpG islands

In vertebrates, cytosine-guanine (CpG) dinucleotides are predominantly methylated, with ∼80% of all CpG sites containing 5-methylcytosine (5mC), a repressive mark associated with long-term gene silencing. The exceptions to such a globally hypermethylated state are CpG-rich DNA sequences called CpG islands (CGIs), which are mostly hypomethylated relative to the bulk genome. CGIs overlap promoters from the earliest vertebrates to humans, indicating a concerted evolutionary drive compatible with CGI retention. CGIs are characterised by DNA sequence features that include DNA hypomethylation, elevated CpG and GC content and the presence of transcription factor binding sites. These sequence characteristics are congruous with the recruitment of transcription factors and chromatin modifying enzymes, and transcriptional activation in general. CGIs colocalize with sites of transcriptional initiation in hypermethylated vertebrate genomes, however, a growing body of evidence indicates that CGIs might exert their gene regulatory function in other genomic contexts. In this review, we discuss the diverse regulatory features of CGIs, their functional readout, and the evolutionary implications associated with CGI retention in vertebrates and possibly in invertebrates.

Comparative Epigenomics Reveals Host Diversity of the Trichinella Epigenomes and Their Effects on Differential Parasitism

Comparative epigenomics provides new insights on evolutionary biology in relation with complex interactions between species and their environments. In the present study, we focus on deciphering the conservation and divergence of DNA methylomes during Trichinella evolution. Whole-genome bisulfite sequencing and RNA-seq were performed on the two clades of Trichinella species, in addition to whole-genome sequencing. We demonstrate that methylation patterns of sing-copy orthologous genes (SCOs) of the 12 Trichinella species are host-related and can mirror known phylogenetic relationships. Among these SCOs, we identify a panel of genes exhibiting hyper-/hypo-methylated features in gene-bodies or respective promoters that play pivotal roles in transcriptome regulation. These hyper-/hypo-methylated SCOs are also of functional significance across developmental stages, as they are highly enriched species-specific and stage-specific expressed genes both in Ad and ML stages. We further identify a set of parasitism-related functional genes that exhibit host-related differential methylation and expression among those SCOs, including p53-like transcription factor and Cdc37 that are of functional significance for elucidating differential parasitology between the two clades of Trichinella. This comparative epigenome study can help to decipher the environmental effects on differential adaptation and parasitism of the genus Trichinella.

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.

Histone Methylation Participates in Gene Expression Control during the Early Development of the Pacific Oyster Crassostrea gigas

Histone methylation patterns are important epigenetic regulators of mammalian development, notably through stem cell identity maintenance by chromatin remodeling and transcriptional control of pluripotency genes. But, the implications of histone marks are poorly understood in distant groups outside vertebrates and ecdysozoan models. However, the development of the Pacific oyster Crassostrea gigas is under the strong epigenetic influence of DNA methylation, and Jumonji histone-demethylase orthologues are highly expressed during C. gigas early life. This suggests a physiological relevance of histone methylation regulation in oyster development, raising the question of functional conservation of this epigenetic pathway in lophotrochozoan. Quantification of histone methylation using fluorescent ELISAs during oyster early life indicated significant variations in monomethyl histone H3 lysine 4 (H3K4me), an overall decrease in H3K9 mono- and tri-methylations, and in H3K36 methylations, respectively, whereas no significant modification could be detected in H3K27 methylation. Early in vivo treatment with the JmjC-specific inhibitor Methylstat induced hypermethylation of all the examined histone H3 lysines and developmental alterations as revealed by scanning electronic microscopy. Using microarrays, we identified 376 genes that were differentially expressed under methylstat treatment, which expression patterns could discriminate between samples as indicated by principal component analysis. Furthermore, Gene Ontology revealed that these genes were related to processes potentially important for embryonic stages such as binding, cell differentiation and development. These results suggest an important physiological significance of histone methylation in the oyster embryonic and larval life, providing, to our knowledge, the first insights into epigenetic regulation by histone methylation in lophotrochozoan development.

Elevated temperature attenuates ovarian functions and induces apoptosis and oxidative stress in the American oyster, Crassostrea virginica: potential mechanisms and signaling pathways

Global climate change is predicted to intensify thermal stress in marine and coastal organisms, affecting their development, growth, and reproductive functions. In this study, we performed histological observations on ovarian development, immunohistochemical analyses of ovarian heat shock protein-70 (HSP70), nitrotyrosine protein (NTP, an indicator of reactive nitrogen species (RNS)), and dinitrophenyl protein (DNP, an indicator of protein oxidation) expressions, in situ TUNEL assay for cellular apoptosis, biochemical analyses of ovarian caspase-3/7 activity and protein carbonyl (PC, a measure of reactive oxygen species (ROS)) contents, nitrate/nitrite (NOx) levels, and extrapallial fluid (EPF, an important body fluid) pH in the American oyster, Crassostrea virginica. Oysters were exposed to medium (28 °C) and high (32 °C) temperatures under controlled laboratory conditions for 1 week. Oysters exposed to higher temperatures significantly decreased the number and diameter of eggs, and EPF protein concentrations compared with controls (24 °C). In contrast, EPF pH, ovarian HSP70 mRNA levels, and protein expression were increased after heat exposure, consistent with increased ovarian apoptosis. The enhanced apoptosis in ovaries was associated with increased ovarian caspase-3/7 activity, PC contents, NOx levels, and NTP and DNP expressions in heat-exposed oysters. Collectively, these results suggest that higher temperatures drastically increase RNS and ROS levels, increasing incidence of apoptosis and subsequently reducing ovarian functions in oysters.

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.

DNA methylation in adults and during development of the self-fertilizing mangrove rivulus, Kryptolebias marmoratus

In addition to genetic variation, epigenetic mechanisms such as DNA methylation might make important contributions to heritable phenotypic diversity in populations. However, it is often difficult to disentangle the contributions of genetic and epigenetic variation to phenotypic diversity. Here, we investigated global DNA methylation and mRNA expression of the methylation-associated enzymes during embryonic development and in adult tissues of one natural isogenic lineage of mangrove rivulus fish, Kryptolebias marmoratus. Being the best-known self-fertilizing hermaphroditic vertebrate affords the opportunity to work with genetically identical individuals to examine, explicitly, the phenotypic effects of epigenetic variance. Using the LUminometric Methylation Assay (LUMA), we described variable global DNA methylation at CpG sites in adult tissues, which differed significantly between hermaphrodite ovotestes and male testes (79.6% and 87.2%, respectively). After fertilization, an immediate decrease in DNA methylation occurred to 15.8% in gastrula followed by re-establishment to 70.0% by stage 26 (liver formation). Compared to zebrafish, at the same embryonic stages, this reprogramming event seems later, deeper, and longer. Furthermore, genes putatively encoding DNA methyltransferases (DNMTs), Ten-Eleven Translocation (TET), and MeCP2 proteins showed specific regulation in adult gonad and brain, and also during early embryogenesis. Their conserved domains and expression profiles suggest that these proteins play important roles during reproduction and development. This study raises questions about mangrove rivulus' peculiar reprogramming period in terms of epigenetic transmission and physiological adaptation of individuals to highly variable environments. In accordance with the general-purpose genotype model, epigenetic mechanisms might allow for the expression of diverse phenotypes among genetically identical individuals. Such phenotypes might help to overcome environmental challenges, making the mangrove rivulus a valuable vertebrate model for ecological epigenetic studies. The mangrove rivulus, Kryptolebias marmoratus, is the best-known self-fertilizing hermaphroditic vertebrate that allows to work with genetically identical individuals to examine, explicitly, the phenotypic effects of epigenetic variance. The reprogramming event is later, more dramatic and longer than in other described vertebrates. High evolutionary conservation and expression patterns of DNMT, TET, and MeCP2 proteins in K. marmoratus suggest biological roles for each member in gametogenesis and development.

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.

Parental diuron-exposure alters offspring transcriptome and fitness in Pacific oyster Crassostrea gigas

One of the primary challenges in ecotoxicology is to contribute to the assessment of the ecological status of ecosystems. In this study, we used Pacific oyster Crassostrea gigas to explore the effects of a parental exposure to diuron, a herbicide frequently detected in marine coastal environments. The present toxicogenomic study provides evidence that exposure of oyster genitors to diuron during gametogenesis results in changes in offspring, namely, transcriptomic profile alterations, increased global DNA methylation levels and reduced growth and survival within the first year of life. Importantly, we highlighted the limitations to identify particular genes or gene expression signatures that could serve as biomarkers for parental herbicide-exposure and further for multigenerational and transgenerational effects of specific chemical stressors. By analyzing samples from two independent experiments, we demonstrated that, due to complex confounding effects with both tested solvent vehicles, diuron non-specifically affected the offspring transcriptome. These original results question the potential development of predictive genomic tools for detecting specific indirect impacts of contaminants in environmental risk assessments. However, our results indicate that chronic environmental exposure to diuron over several generations may have significant long term impacts on oyster populations with adverse health outcomes.

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.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Bisulfite Sequencing with Daphnia Highlights a Role for Epigenetics in Regulating Stress Response to Microcystis through Preferential Differential Methylation of Serine and Threonine Amino Acids

Little is known about the influence that environmental stressors may have on genome-wide methylation patterns, and to what extent epigenetics may be involved in environmental stress response. Yet, studies of methylation patterns under stress could provide crucial insights on stress response and toxicity pathways. Here, we focus on genome-wide methylation patterns in the microcrustacean Daphnia magna, a model organism in ecotoxicology and risk assessment, exposed to the toxic cyanobacterium Microcystis aeruginosa. Bisulfite sequencing of exposed and control animals highlighted differential methylation patterns in Daphnia upon exposure to Microcystis primarily in exonic regions. These patterns are enriched for serine/threonine amino acid codons and genes related to protein synthesis, transport and degradation. Furthermore, we observed that genes with differential methylation corresponded well with genes susceptible to alternative splicing in response to Microcystis stress. Overall, our results suggest a complex mechanistic response in Daphnia characterized by interactions between DNA methylation and gene regulation mechanisms. These results underscore that DNA methylation is modulated by environmental stress and can also be an integral part of the toxicity response in our study species.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Distribution of CpG Motifs in Upstream Gene Domains in a Reef Coral and Sea Anemone: Implications for Epigenetics in Cnidarians

Coral reefs are under assault from stressors including global warming, ocean acidification, and urbanization. Knowing how these factors impact the future fate of reefs requires delineating stress responses across ecological, organismal and cellular scales. Recent advances in coral reef biology have integrated molecular processes with ecological fitness and have identified putative suites of temperature acclimation genes in a Scleractinian coral Acropora hyacinthus. We wondered what unique characteristics of these genes determined their coordinate expression in response to temperature acclimation, and whether or not other corals and cnidarians would likewise possess these features. Here, we focus on cytosine methylation as an epigenetic DNA modification that is responsive to environmental stressors. We identify common conserved patterns of cytosine-guanosine dinucleotide (CpG) motif frequencies in upstream promoter domains of different functional gene groups in two cnidarian genomes: a coral (Acropora digitifera) and an anemone (Nematostella vectensis). Our analyses show that CpG motif frequencies are prominent in the promoter domains of functional genes associated with environmental adaptation, particularly those identified in A. hyacinthus. Densities of CpG sites in upstream promoter domains near the transcriptional start site (TSS) are 1.38x higher than genomic background levels upstream of -2000 bp from the TSS. The increase in CpG usage suggests selection to allow for DNA methylation events to occur more frequently within 1 kb of the TSS. In addition, observed shifts in CpG densities among functional groups of genes suggests a potential role for epigenetic DNA methylation within promoter domains to impact functional gene expression responses in A. digitifera and N. vectensis. Identifying promoter epigenetic sequence motifs among genes within specific functional groups establishes an approach to describe integrated cellular responses to environmental stress in reef corals and potential roles of epigenetics on survival and fitness in the face of global climate change.

Inheritance and Variation of Genomic DNA Methylation in Diploid and Triploid Pacific Oyster (Crassostrea gigas)

DNA methylation is an important epigenetic mechanism that could be responsive to environmental changes indicating a potential role in natural selection and adaption. In order to evaluate an evolutionary role of DNA methylation, it is essential to first gain a better insight into inheritability. To address this question, this study investigated DNA methylation variation from parents to offspring in the Pacific oyster Crassostrea gigas using fluorescent-labeled methylation-sensitive amplified polymorphism (F-MSAP) analysis. Most of parental methylated loci were stably transmitted to offspring segregating following Medelian expectation. However, methylated loci deviated more often than non-methylated loci and offspring showed a few de novo methylated loci indicating DNA methylation changes from parents to offspring. Interestingly, some male-specific methylated loci were found in this study which might help to explore sex determination in oyster. Despite environmental stimuli, genomic stresses such as polyploidization also can induce methylation changes. This study also compared global DNA methylation level and individual methylated loci between diploid and triploid oysters. Results showed no difference in global methylation state but a few ploidy-specific loci were detected. DNA methylation variation during polyploidization was less than autonomous methylation variation from parents to offspring.

GigaTON: an extensive publicly searchable database providing a new reference transcriptome in the pacific oyster Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is one of the most important aquaculture shellfish resources worldwide. Important efforts have been undertaken towards a better knowledge of its genome and transcriptome, which makes now C. gigas becoming a model organism among lophotrochozoans, the under-described sister clade of ecdysozoans within protostomes. These massive sequencing efforts offer the opportunity to assemble gene expression data and make such resource accessible and exploitable for the scientific community. Therefore, we undertook this assembly into an up-to-date publicly available transcriptome database: the GigaTON (Gigas TranscriptOme pipeliNe) database.

DESCRIPTION

We assembled 2204 million sequences obtained from 114 publicly available RNA-seq libraries that were realized using all embryo-larval development stages, adult organs, different environmental stressors including heavy metals, temperature, salinity and exposure to air, which were mostly performed as part of the Crassostrea gigas genome project. This data was analyzed in silico and resulted into 56621 newly assembled contigs that were deposited into a publicly available database, the GigaTON database. This database also provides powerful and user-friendly request tools to browse and retrieve information about annotation, expression level, UTRs, splice and polymorphism, and gene ontology associated to all the contigs into each, and between all libraries.

CONCLUSIONS

The GigaTON database provides a convenient, potent and versatile interface to browse, retrieve, confront and compare massive transcriptomic information in an extensive range of conditions, tissues and developmental stages in Crassostrea gigas. To our knowledge, the GigaTON database constitutes the most extensive transcriptomic database to date in marine invertebrates, thereby a new reference transcriptome in the oyster, a highly valuable resource to physiologists and evolutionary biologists.

Persistence of Positive Carryover Effects in the Oyster, Saccostrea glomerata, following Transgenerational Exposure to Ocean Acidification

Ocean acidification (OA) is predicted to have widespread implications for marine organisms, yet the capacity for species to acclimate or adapt over this century remains unknown. Recent transgenerational studies have shown that for some marine species, exposure of adults to OA can facilitate positive carryover effects to their larval and juvenile offspring that help them to survive in acidifying oceanic conditions. But whether these positive carryover effects can persist into adulthood or the next generation is unknown. Here we tested whether positive carryover effects found in larvae of the oyster, Saccostrea glomerata following transgenerational exposure to elevated CO2, could persist into adulthood and whether subsequent transgenerational exposure of adults to elevated CO2 would facilitate similar adaptive responses in the next generation of larvae and juveniles. Following our previous transgenerational exposure of parental adults and first generation (F1) larvae to ambient (385 μatm) and elevated (856 μatm) CO2, newly settled F1 juveniles were transferred to the field at ambient CO2 for 14 months, until they reached reproductive maturity. At this time, the F1 adults were returned to the laboratory and the previous transgenerational CO2 exposure was repeated to produce F2 offspring. We found that the capacity of adults to regulate extracellular pH at elevated CO2 was improved if they had a prior history of transgenerational exposure to elevated CO2. In addition, subsequent transgenerational exposure of these adults led to an increase in the resilience of their larval and juvenile offspring. Offspring with a history of transgenerational exposure to elevated CO2 had a lower percentage abnormality, faster development rate, faster shell growth and increased heart rate at elevated CO2 compared with F2 offspring with no prior history of exposure to elevated CO2. Our results suggest that positive carryover effects originating during parental and larval exposure will be important in mediating some of the impacts of OA for later life-history stages and generations.