Mapping the zebrafish brain methylome using reduced representation bisulfite sequencing

Reduced representation bisulfite sequencing (RRBS) has been used to profile DNA methylation patterns in mammalian genomes such as human, mouse and rat. The methylome of the zebrafish, an important animal model, has not yet been characterized at base-pair resolution using RRBS. Therefore, we evaluated the technique of RRBS in this model organism by generating four single-nucleotide resolution DNA methylomes of adult zebrafish brain. We performed several simulations to show the distribution of fragments and enrichment of CpGs in different in silico reduced representation genomes of zebrafish. Four RRBS brain libraries generated 98 million sequenced reads and had higher frequencies of multiple mapping than equivalent human RRBS libraries. The zebrafish methylome indicates there is higher global DNA methylation in the zebrafish genome compared with its equivalent human methylome. This observation was confirmed by RRBS of zebrafish liver. High coverage CpG dinucleotides are enriched in CpG island shores more than in the CpG island core. We found that 45% of the mapped CpGs reside in gene bodies, and 7% in gene promoters. This analysis provides a roadmap for generating reproducible base-pair level methylomes for zebrafish using RRBS and our results provide the first evidence that RRBS is a suitable technique for global methylation analysis in zebrafish.     

DNA methylation analysis by pyrosequencing

Pyrosequencing is a sequencing-by-synthesis method that quantitatively monitors the real-time incorporation of nucleotides through the enzymatic conversion of released pyrophosphate into a proportional light signal. Quantitative measures are of special importance for DNA methylation analysis in various developmental and pathological situations. Analysis of DNA methylation patterns by pyrosequencing combines a simple reaction protocol with reproducible and accurate measures of the degree of methylation at several CpGs in close proximity with high quantitative resolution. After bisulfite treatment and PCR, the degree of each methylation at each CpG position in a sequence is determined from the ratio of T and C. The process of purification and sequencing can be repeated for the same template to analyze other CpGs in the same amplification product. Quantitative epigenotypes are obtained using this protocol in approximately 4 h for up to 96 DNA samples when bisulfite-treated DNA is already available as the starting material.     

A rapid, quantitative, non-radioactive bisulfite-SNuPE- IP RP HPLC assay for methylation analysis at specific CpG sites

The precise mapping and quantification of DNA methylation as an epigenetic parameter during development and in diseased tissues is of great importance for functional genomics. Here we describe a rapid, quantitative method to assess methylation levels at specific CpG sites using PCR products of bisulfite-treated genomic DNA. Using single nucleotide primer extension (SNuPE) assays in combination with ion pair reverse phase high performance liquid chromatography (IP RP HPLC) separation techniques, methylated and unmethylated CpGs can be discriminated and quantified based on the different masses and hydrophobicities of the extended primer products. The assay is linear, highly reproducible and several sites can be measured simultaneously in one reaction. It can be semi-automated and eliminates the need for cloning and sequencing of individual bisulfite PCR products.     

Collaboration between CpG sites is needed for stable somatic inheritance of DNA methylation states

Inheritance of 5-methyl cytosine modification of CpG (CG/CG) DNA sequences is needed to maintain early developmental decisions in vertebrates. The standard inheritance model treats CpGs as independent, with methylated CpGs maintained by efficient methylation of hemimethylated CpGs produced after DNA replication, and unmethylated CpGs maintained by an absence of de novo methylation. By stochastic simulations of CpG islands over multiple cell cycles and systematic sampling of reaction parameters, we show that the standard model is inconsistent with many experimental observations. In contrast, dynamic collaboration between CpGs can provide strong error-tolerant somatic inheritance of both hypermethylated and hypomethylated states of a cluster of CpGs, reproducing observed stable bimodal methylation patterns. Known recruitment of methylating enzymes by methylated CpGs could provide the necessary collaboration, but we predict that recruitment of demethylating enzymes by unmethylated CpGs strengthens inheritance and allows CpG islands to remain hypomethylated within a sea of hypermethylation.     

DNA Methylation Analysis of Chromosome 21 Gene Promoters at Single Base Pair and Single Allele Resolution

Differential DNA methylation is an essential epigenetic signal for gene regulation, development, and disease processes. We mapped DNA methylation patterns of 190 gene promoter regions on chromosome 21 using bisulfite conversion and subclone sequencing in five human cell types. A total of 28,626 subclones were sequenced at high accuracy using (long-read) Sanger sequencing resulting in the measurement of the DNA methylation state of 580427 CpG sites. Our results show that average DNA methylation levels are distributed bimodally with enrichment of highly methylated and unmethylated sequences, both for amplicons and individual subclones, which represent single alleles from individual cells. Within CpG-rich sequences, DNA methylation was found to be anti-correlated with CpG dinucleotide density and GC content, and methylated CpGs are more likely to be flanked by AT-rich sequences. We observed over-representation of CpG sites in distances of 9, 18, and 27 bps in highly methylated amplicons. However, DNA sequence alone is not sufficient to predict an amplicon's DNA methylation status, since 43% of all amplicons are differentially methylated between the cell types studied here. DNA methylation in promoter regions is strongly correlated with the absence of gene expression and low levels of activating epigenetic marks like H3K4 methylation and H3K9 and K14 acetylation. Utilizing the single base pair and single allele resolution of our data, we found that i) amplicons from different parts of a CpG island frequently differ in their DNA methylation level, ii) methylation levels of individual cells in one tissue are very similar, and iii) methylation patterns follow a relaxed site-specific distribution. Furthermore, iv) we identified three cases of allele-specific DNA methylation on chromosome 21. Our data shed new light on the nature of methylation patterns in human cells, the sequence dependence of DNA methylation, and its function as epigenetic signal in gene regulation. Further, we illustrate genotype–epigenotype interactions by showing novel examples of allele-specific methylation.     

DNA Methylation: Bisulphite Modification and Analysis

Epigenetics describes the heritable changes in gene function that occur independently to the DNA sequence. The molecular basis of epigenetic gene regulation is complex, but essentially involves modifications to the DNA itself or the proteins with which DNA associates. The predominant epigenetic modification of DNA in mammalian genomes is methylation of cytosine nucleotides (5-MeC). DNA methylation provides instruction to gene expression machinery as to where and when the gene should be expressed. The primary target sequence for DNA methylation in mammals is 5''-CpG-3'' dinucleotides (Figure 1). CpG dinucleotides are not uniformly distributed throughout the genome, but are concentrated in regions of repetitive genomic sequences and CpG "islands" commonly associated with gene promoters (Figure 1). DNA methylation patterns are established early in development, modulated during tissue specific differentiation and disrupted in many disease states including cancer. To understand the biological role of DNA methylation and its role in human disease, precise, efficient and reproducible methods are required to detect and quantify individual 5-MeCs.This protocol for bisulphite conversion is the "gold standard" for DNA methylation analysis and facilitates identification and quantification of DNA methylation at single nucleotide resolution. The chemistry of cytosine deamination by sodium bisulphite involves three steps (Figure 2). (1) Sulphonation: The addition of bisulphite to the 5-6 double bond of cytosine (2) Hydrolic Deamination: hydrolytic deamination of the resulting cytosine-bisulphite derivative to give a uracil-bisulphite derivative (3) Alkali Desulphonation: Removal of the sulphonate group by an alkali treatment, to give uracil. Bisulphite preferentially deaminates cytosine to uracil in single stranded DNA, whereas 5-MeC, is refractory to bisulphite-mediated deamination. Upon PCR amplification, uracil is amplified as thymine while 5-MeC residues remain as cytosines, allowing methylated CpGs to be distinguished from unmethylated CpGs by presence of a cytosine "C" versus thymine "T" residue during sequencing.DNA modification by bisulphite conversion is a well-established protocol that can be exploited for many methods of DNA methylation analysis. Since the detection of 5-MeC by bisulphite conversion was first demonstrated by Frommer et al.¹ and Clark et al.², methods based around bisulphite conversion of genomic DNA account for the majority of new data on DNA methylation. Different methods of post PCR analysis may be utilized, depending on the degree of specificity and resolution of methylation required. Cloning and sequencing is still the most readily available method that can give single nucleotide resolution for methylation across the DNA molecule.     

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.     

Deep Bisulfite Sequencing of Aberrantly Methylated Loci in a Patient with Multiple Methylation Defects

NLRP7 is a maternal effect gene as maternal mutations in this gene cause recurrent hydatidiform moles, spontaneous abortions and stillbirths, whereas live births are very rare. We have studied a patient with multiple anomalies born to a mother with a heterozygous NLRP7 mutation. By array-based CpG methylation analysis of blood DNA from the patient, his parents and 18 normal controls on Illumina Infinium HumanMethylation27 BeadChips we found that the patient had methylation changes (delta ß ≥ 0.3) at many imprinted loci as well as at 87 CpGs associated with 85 genes of unknown imprinting status. Using a pseudoproband (permutation) approach, we found methylation changes at only 7-24 CpGs (mean 15; standard deviation 4.84) in the controls. Thus, the number of abberantly methylated CpGs in the patient is more than 14 standard deviations higher. In order to identify novel imprinted genes among the 85 conspicuous genes in the patient, we selected 19 (mainly hypomethylated) genes for deep bisulfite amplicon sequencing on the ROCHE/454 Genome Sequencer in the patient and at least two additional controls. These controls had not been included in the array analysis and were heterozygous for a single nucleotide polymorphism at the test locus, so that allele-specific DNA methylation patterns could be determined. Apart from FAM50B, which we proved to be imprinted in blood, we did not observe allele-specific DNA methylation at the other 18 loci. We conclude that the patient does not only have methylation defects at imprinted loci but (at least in blood) also an excess of methylation changes at apparently non-imprinted loci.     

Regionally Specific and Genome-Wide Analyses Conclusively Demonstrate the Absence of CpG Methylation in Human Mitochondrial DNA

Although CpG methylation clearly distributes genome-wide in vertebrate nuclear DNA, the state of methylation in the vertebrate mitochondrial genome has been unclear. Several recent reports using immunoprecipitation, mass spectrometry, and enzyme-linked immunosorbent assay methods concluded that human mitochondrial DNA (mtDNA) has much more than the 2 to 5% CpG methylation previously estimated. However, these methods do not provide information as to the sites or frequency of methylation at each CpG site. Here, we have used the more definitive bisulfite genomic sequencing method to examine CpG methylation in HCT116 human cells and primary human cells to independently answer these two questions. We found no evidence of CpG methylation at a biologically significant level in these regions of the human mitochondrial genome. Furthermore, unbiased next-generation sequencing of sodium bisulfite treated total DNA from HCT116 cells and analysis of genome-wide sodium bisulfite sequencing data sets from several other DNA sources confirmed this absence of CpG methylation in mtDNA. Based on our findings using regionally specific and genome-wide approaches with multiple human cell sources, we can definitively conclude that CpG methylation is absent in mtDNA. It is highly unlikely that CpG methylation plays any role in direct control of mitochondrial function.     

Epigenetic profile of developmentally important genes in bovine oocytes

Assisted reproductive technologies are associated with an increased incidence of epigenetic aberrations, specifically in imprinted genes. Here, we used the bovine oocyte as a model to determine putative epigenetic mutations at three imprinted gene loci caused by the type of maturation, either in vitro maturation (IVM) in Tissue Culture Medium 199 (TCM) or modified synthetic oviduct fluid (mSOF) medium, or in vivo maturation. We applied a limiting dilution approach and direct bisulfite sequencing to analyze the methylation profiles of individual alleles (DNA molecules) for H19/IGF2, PEG3, and SNRPN, which are each associated with imprinting defects in humans and/or the mouse model, and are known to be differentially methylated in bovine embryos. Altogether, we obtained the methylation patterns of 203 alleles containing 4,512 CpG sites from immature oocytes, 213 alleles with 4,779 CpG sites from TCM-matured oocytes, 215 alleles/4,725 CpGs in mSOF-matured oocytes, and 78 alleles/1,672 CpGs from in vivo-matured oocytes. The total rate of individual CpGs and entire allele methylation errors did not differ significantly between the two IVM and the in vivo group, indicating that current IVM protocols have no or only marginal effects on these critical epigenetic marks. Furthermore, the mRNA expression profiles of the three imprinted genes and a panel of eight other genes indicative of oocyte competence were determined by quantitative real-time PCR. We found different mRNA expression profiles between in vivo-matured oocytes versus their in vitro-matured counterparts, suggesting an influence on regulatory mechanisms other than DNA methylation.