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

Background

Studies of DNA methylomes in a wide range of eukaryotes have revealed both conserved and divergent characteristics of DNA methylation among phylogenetic groups. However, data on invertebrates particularly molluscs are limited, which hinders our understanding of the evolution of DNA methylation in metazoa. The sequencing of the Pacific oyster Crassostrea gigas genome provides an opportunity for genome-wide profiling of DNA methylation in this model mollusc.

Results

Homologous searches against the C. gigas genome identified functional orthologs for key genes involved in DNA methylation: DNMT1, DNMT2, DNMT3, MBD2/3 and UHRF1. Whole-genome bisulfite sequencing (BS-seq) of the oyster’s mantle tissues revealed that more than 99% methylation modification was restricted to cytosines in CpG context and methylated CpGs accumulated in the bodies of genes that were moderately expressed. Young repeat elements were another major targets of CpG methylation in oysters. Comparison with other invertebrate methylomes suggested that the 5’-end bias of gene body methylation and the negative correlation between gene body methylation and gene length were the derived features probably limited to the insect lineage. Interestingly, phylostratigraphic analysis showed that CpG methylation preferentially targeted genes originating in the common ancestor of eukaryotes rather than the oldest genes originating in the common ancestor of cellular organisms.

Conclusions

Comparative analysis of the oyster DNA methylomes and that of other animal species revealed that the characteristics of DNA methylation were generally conserved during invertebrate evolution, while some unique features were derived in the insect lineage. The preference of methylation modification on genes originating in the eukaryotic ancestor rather than the oldest genes is unexpected, probably implying that the emergence of methylation regulation in these ''relatively young’ genes was critical for the origin and radiation of eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1119) contains supplementary material, which is available to authorized users.     

Methylation-Sensitive Expression of a DNA Demethylase Gene Serves As an Epigenetic Rheostat

Genomes must balance active suppression of transposable elements (TEs) with the need to maintain gene expression. In Arabidopsis, euchromatic TEs are targeted by RNA-directed DNA methylation (RdDM). Conversely, active DNA demethylation prevents accumulation of methylation at genes proximal to these TEs. It is unknown how a cellular balance between methylation and demethylation activities is achieved. Here we show that both RdDM and DNA demethylation are highly active at a TE proximal to the major DNA demethylase gene ROS1. Unexpectedly, and in contrast to most other genomic targets, expression of ROS1 is promoted by DNA methylation and antagonized by DNA demethylation. We demonstrate that inducing methylation in the ROS1 proximal region is sufficient to restore ROS1 expression in an RdDM mutant. Additionally, methylation-sensitive expression of ROS1 is conserved in other species, suggesting it is adaptive. We propose that the ROS1 locus functions as an epigenetic rheostat, tuning the level of demethylase activity in response to methylation alterations, thus ensuring epigenomic stability.     

Reduced activity of Arabidopsis chromosome-cohesion regulator gene CTF7/ECO1 alters cytosine methylation status and retrotransposon expression

Multicellular organisms such as higher plants require timely regulation of DNA replication and cell division to grow and develop. Recent work in Arabidopsis has shown that chromosome segregation during meiosis and mitosis depends on the activity of several genes that in yeast are involved in the establishment of chromosomal cohesion. In this process, proteins of the STRUCTURAL MAINTENANCE OF CHROMOSOMES (SMC) family tether chromosomes and establish inter- and intrachromosomal connections. In Arabidopsis, recruitment of SMC proteins and establishment of cohesion during key stages of the cell cycle depend on the activity of CHROMOSOME TRANSMISSION FIDELITY 7/ESTABLISHMENT OF COHESION 1 (CTF7/ECO1). Here we show that loss of CTF7/ECO1 activity alters the status of cytosine methylation in both intergenic regions and transposon loci. An increase in expression was also observed for transposon copia28, which suggests a link between CTF7/ECO1 activity, DNA methylation and gene silencing. More work is needed to determine the mechanistic relationships that intervene in this process.     

Characterization of a Tc1-like transposable element in zebrafish (Danio rerio)

We have characterized Tdr1, a family of Tc1-like transposable elements found in the genome of zebrafish (Danio rerio). The copy number and distribution of the sequence in the zebrafish genome have been determined, and by these criteria Tdr1 can be classified as a moderately repetitive, interspersed element. Examination of the sequences and structures of several copies of Tdr1 revealed that a particular deletion derivative, 1250 by long, of the transposon has been amplified to become the dominant form of Tdr1. The deletion in these elements encompasses sequences encoding the N-terminal portion of the putative Tdr1 transposase. Sequences corresponding to the deleted region were also detected, and thus allowed prediction of the nucleotide sequence of a hypothetical full-length element. Well conserved segments of Tc1-like transposons were found in the flanking regions of known fish genes, suggesting that these elements have a long evolutionary history in piscine genomes. Tdr1 elements have long, 208 by inverted repeats, with a short DNA motif repeated four times at the termini of the inverted repeats. Although different from that of the prototype C. elegans transposon Tc1, this inverted repeat structure is shared by transposable elements from salmonid fish species and two Drosophila species. We propose that these transposons form a subgroup within the Tc1-like family. Comparison of Tc1-like transposons supports the hypothesis that the transposase genes and their flanking sequences have been shaped by independent evolutionary constraints. Although Tc1-like sequences are present in the genomes of several strains of zebrafish and in salmonid fishes, these sequences are not conserved in the genus Danio, thus raising the possibility that these elements can be exploited for gene tagging and genome mapping.     

AtMBD6, a methyl CpG binding domain protein,maintains gene silencing in <Emphasis Type="Italic">Arabidopsis</Emphasis> by interacting with RNA binding proteins

DNA methylation, mediated by double-stranded RNA, is a conserved epigenetic phenomenon that protects a genome from transposons, silences unwanted genes and has a paramount function in plant or animal development. Methyl CpG binding domain proteins are members of a class of proteins that bind to methylated DNA. The Arabidopsis thaliana genome encodes 13 methyl CpG binding domain (MBD) proteins, but the molecular/biological functions of most of these proteins are still not clear. In the present study, we identified four proteins that interact with AtMBD6. Interestingly, three of them contain RNA binding domains and are co-localized with AtMBD6 in the nucleus. The interacting partners includes AtRPS2C (a 40S ribosomal protein), AtNTF2 (nuclear transport factor 2) and AtAGO4 (Argonoute 4). The fourth protein that physically interacts with AtMBD6 is a histone-modifying enzyme, histone deacetylase 6 (AtHDA6), which is a known component of the RNA-mediated gene silencing system. Analysis of genomic DNA methylation in the atmbd6, atrps2c and atntf2 mutants, using methylation-sensitive PCR detected decreased DNA methylation at miRNA/siRNA producing loci, pseudogenes and other targets of RNA-directed DNA methylation. Our results indicate that AtMBD6 is involved in RNA-mediated gene silencing and it binds to RNA binding proteins like AtRPS2C, AtAGO4 and AtNTF2. AtMBD6 also interacts with histone deacetylase AtHDA6 that might have a role in chromatin condensation at the targets of RdDM.     

Genetic Perturbation of the Maize Methylome

DNA methylation can play important roles in the regulation of transposable elements and genes. A collection of mutant alleles for 11 maize (Zea mays) genes predicted to play roles in controlling DNA methylation were isolated through forward- or reverse-genetic approaches. Low-coverage whole-genome bisulfite sequencing and high-coverage sequence-capture bisulfite sequencing were applied to mutant lines to determine context- and locus-specific effects of these mutations on DNA methylation profiles. Plants containing mutant alleles for components of the RNA-directed DNA methylation pathway exhibit loss of CHH methylation at many loci as well as CG and CHG methylation at a small number of loci. Plants containing loss-of-function alleles for chromomethylase (CMT) genes exhibit strong genome-wide reductions in CHG methylation and some locus-specific loss of CHH methylation. In an attempt to identify stocks with stronger reductions in DNA methylation levels than provided by single gene mutations, we performed crosses to create double mutants for the maize CMT3 orthologs, Zmet2 and Zmet5, and for the maize DDM1 orthologs, Chr101 and Chr106. While loss-of-function alleles are viable as single gene mutants, the double mutants were not recovered, suggesting that severe perturbations of the maize methylome may have stronger deleterious phenotypic effects than in Arabidopsis thaliana.     

Epigenetic changes and transposon reactivation in Thai rice hybrids

Inter- or intraspecific hybridization is the first step in transferring exogenous traits to the germplasm of a recipient crop. One of the complicating factors is the occurrence of epigenetic modifications of the hybrids, which in turn can change their gene expression and phenotype. In this study we present an analysis of epigenome changes in rice hybrids that were obtained by crossing rice cultivars, most of them of indica type and Thai origin. Comparing amplified fragment length polymorphism (AFLP) fingerprints of twenty-four cultivars, we calculated Nei’s indexes for measuring genetic relationships. Epigenetic changes in their hybrids were established using methylation-sensitive AFLP fingerprinting and transposon display of the rice transposable elements (TEs) Stowaway Os-1 and Mashu, leading to the question whether the relationship between parental genomes is a predictor of epigenome changes, TE reactivation and changes in TE methylation. Our study now reveals that the genetic relationship between the parents and DNA methylation changes in their hybrids is not significantly correlated. Moreover, genetic distance correlates only weakly with Mashu reactivation, whereas it does not correlate with Stowaway Os-1 reactivation. Our observations also suggest that epigenome changes in the hybrids are localized events affecting specific chromosomal regions and transposons rather than affecting the genomic methylation landscape as a whole. The weak correlation between genetic distance and Mashu methylation and reactivation points at only limited influence of genetic background on the epigenetic status of the transposon. Our study further demonstrates that hybridizations between and among specific japonica and indica cultivars induce both genomic DNA methylation and reactivation/methylation change in the Stowaway Os-1 and Mashu transposons. The observed epigenetic changes seem to affect the transposons in a clear manner, partly driven by stochastic processes, which may account for a broader phenotypic plasticity of the hybrids. A better understanding of the epigenome changes leading to such transposon activation can lead to the development of novel tools for more variability in future rice breeding.     

Advances in biotechnology and informatics to link variation in the genome to phenotypes in plants and animals

Advances in our understanding of genome structure provide consistent evidence for the existence of a core genome representing species classically defined by phenotype, as well as conditionally dispensable components of the genome that shows extensive variation between individuals of a given species. Generally, conservation of phenotypic features between species reflects conserved features of the genome; however, this is evidently not necessarily always the case as demonstrated by the analysis of the tunicate chordate Oikopleura dioica. In both plants and animals, the methylation activity of DNA and histones continues to present new variables for modifying (eventually) the phenotype of an organism and provides for structural variation that builds on the point mutations, rearrangements, indels, and amplification of retrotransposable elements traditionally considered. The translation of the advances in the structure/function analysis of the genome to industry is facilitated through the capture of research outputs in “toolboxes” that remain accessible in the public domain.     

A longitudinal study of epigenetic variation in twins

DNA methylation is a key epigenetic mechanism involved in the developmental regulation of gene expression. Alterations in DNA methylation are established contributors to inter-individual phenotypic variation and have been associated with disease susceptibility. The degree to which changes in loci-specific DNA methylation are under the influence of heritable and environmental factors is largely unknown. In this study, we quantitatively measured DNA methylation across the promoter regions of the dopamine receptor 4 gene (DRD4), the serotonin transporter gene (SLC6A4/SERT) and the X-linked monoamine oxidase A gene (MAOA) using DNA sampled at both ages 5 and 10 years in 46 MZ twinpairs and 45 DZ twin-pairs (total n = 182). Our data suggest that DNA methylation differences are apparent already in early childhood, even between genetically identical individuals, and that individual differences in methylation are not stable over time. Our longitudinal-developmental study suggests that environmental influences are important factors accounting for interindividual DNA methylation differences, and that these influences differ across the genome. The observation of dynamic changes in DNA methylation over time highlights the importance of longitudinal research designs for epigenetic research.Key words: epigenetics, DNA methylation, twin, heritability, dynamic, environment     

Extensive Natural Epigenetic Variation at a De Novo Originated Gene

Epigenetic variation, such as heritable changes of DNA methylation, can affect gene expression and thus phenotypes, but examples of natural epimutations are few and little is known about their stability and frequency in nature. Here, we report that the gene Qua-Quine Starch (QQS) of Arabidopsis thaliana, which is involved in starch metabolism and that originated de novo recently, is subject to frequent epigenetic variation in nature. Specifically, we show that expression of this gene varies considerably among natural accessions as well as within populations directly sampled from the wild, and we demonstrate that this variation correlates negatively with the DNA methylation level of repeated sequences located within the 5′end of the gene. Furthermore, we provide extensive evidence that DNA methylation and expression variants can be inherited for several generations and are not linked to DNA sequence changes. Taken together, these observations provide a first indication that de novo originated genes might be particularly prone to epigenetic variation in their initial stages of formation.     

Exploring the Roles of DNA Methylation in the Metal-Reducing Bacterium Shewanella oneidensis MR-1

We performed whole-genome analyses of DNA methylation in Shewanella oneidensis MR-1 to examine its possible role in regulating gene expression and other cellular processes. Single-molecule real-time (SMRT) sequencing revealed extensive methylation of adenine (N6mA) throughout the genome. These methylated bases were located in five sequence motifs, including three novel targets for type I restriction/modification enzymes. The sequence motifs targeted by putative methyltranferases were determined via SMRT sequencing of gene knockout mutants. In addition, we found that S. oneidensis MR-1 cultures grown under various culture conditions displayed different DNA methylation patterns. However, the small number of differentially methylated sites could not be directly linked to the much larger number of differentially expressed genes under these conditions, suggesting that DNA methylation is not a major regulator of gene expression in S. oneidensis MR-1. The enrichment of methylated GATC motifs in the origin of replication indicates that DNA methylation may regulate genome replication in a manner similar to that seen in Escherichia coli. Furthermore, comparative analyses suggest that many Gammaproteobacteria, including all members of the Shewanellaceae family, may also utilize DNA methylation to regulate genome replication.     

Assessing the Impact of Transgenerational Epigenetic Variation on Complex Traits

Loss or gain of DNA methylation can affect gene expression and is sometimes transmitted across generations. Such epigenetic alterations are thus a possible source of heritable phenotypic variation in the absence of DNA sequence change. However, attempts to assess the prevalence of stable epigenetic variation in natural and experimental populations and to quantify its impact on complex traits have been hampered by the confounding effects of DNA sequence polymorphisms. To overcome this problem as much as possible, two parents with little DNA sequence differences, but contrasting DNA methylation profiles, were used to derive a panel of epigenetic Recombinant Inbred Lines (epiRILs) in the reference plant Arabidopsis thaliana. The epiRILs showed variation and high heritability for flowering time and plant height (~30%), as well as stable inheritance of multiple parental DNA methylation variants (epialleles) over at least eight generations. These findings provide a first rationale to identify epiallelic variants that contribute to heritable variation in complex traits using linkage or association studies. More generally, the demonstration that numerous epialleles across the genome can be stable over many generations in the absence of selection or extensive DNA sequence variation highlights the need to integrate epigenetic information into population genetics studies.     

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.     

Dnmt2 is the most evolutionary conserved and enigmatic cytosine DNA methyltransferase in eukaryotes

Dnmt2 is the most strongly conserved cytosine DNA methyltransferase in eukaryotes. It has been found in all organisms possessing methyltransferases of the Dnmt1 and Dnmt3 families, whereas in many others Dnmt2 is the sole cytosine DNA methyltransferase. The Dnmt2 molecule contains all conserved motifs of cytosine DNA methyltransferases. It forms 3D complexes with DNA very similar to those of bacterial DNA methyltransferases and performs cytosine methylation by a catalytic mechanism common to all cytosine DNA methyltransferases. Catalytic activity of the purified Dnmt2 with DNA substrates is very low and could hardly be detected in direct biochemical assays. Dnmt2 is the sole cytosine DNA methyltransferase in Drosophila and other dipteran insects. Its overexpression as a transgene leads to DNA hypermethylation in all sequence contexts and to an extended life span. On the contrary, a null-mutation of the Dnmt2 gene leads to a diminished life span, though no evident anomalies in development are observed. Dnmt2 is also the sole cytosine DNA methyltransferase in several protists. Similar to Drosophila these protists have a very low level of DNA methylation. Some limited genome compartments, such as transposable sequences, are probably the methylation targets in these organisms. Dnmt2 does not participate in genome methylation in mammals, but seems to be an RNA methyltransferase modifying the 38th cytosine residue in anticodon loop of certain tRNAs. This modification enhances stability of tRNAs, especially in stressful conditions. Dnmt2 is the only enzyme known to perform RNA methylation by a catalytic mechanism characteristic of DNA methyltransferases. The Dnmt2 activity has been shown in mice to be necessary for paramutation establishment, though the precise mechanisms of its participation in this form of epigenetic heredity are unknown. It seems likely, that either of the two Dnmt2 activities could become a predominant one during the evolution of different species. The high level of the Dnmt2 evolutionary conservation proves its activity to have a significant adaptive value in natural environment.