Estimating the total mouse DNA methylation according to the B1 repetitive elements

Measuring the degree of methylation of the B1 element in mouse may represent the global DNA methylation status because about 30,000 copies of the B1 element are randomly dispersed in the total mouse genome. Six CpG dinucleotides are located within each 163 bp size of B1 element, and each CpG dinucleotide was partially methylated. We quantitated the DNA methylation of the B1 repetitive elements by performing PCR for the methylation specific PCR (MSP) and also by the pyrosequencing. Each CpG dinucleotide was methylated at an average of 9% in the mouse genome by the pyrosequencing and MSP. Especially, we checked whether CpG methylation of the B1 element could respond to a treatment of the DNA methylation inhibitor, 5-azacytidine (5-AzaC). Consequently, the calibration graph resulting from measuring the relative CpG methylation percentage of the B1 element is linearly decreased with the increasing amount of 5-AzaC (up to 50 ng/ml concentration) in the NIH3T3 cells with a standard deviation of only 1.73% between three independent assays. Our methods can be applied to the routine analysis of the global DNA methylation changes in mouse in vivo and in vitro in pharmaceuticals and basic epigenetic research with efforts being less labor-intensive.     

Distribution of DNA methylation,CpGs, and CpG islands in human isochores

DNA methylation is a major epigenetic modification of the genome that affects basic biological functions, such as gene expression and cell development. We used the human genome sequences and the DNA methylation data that are available in order to establish a map of the levels of GC and methylation in isochores. We also looked for the correlations that hold between GC levels and the distribution of the (1) dinucleotide CpG, (2) ratio 5mC/CpG, and (3) CpG islands. Our results show that methylation levels, CpG frequencies, and the density of CpG islands are positively correlated with the GC level of isochores. In contrast, the correlation between the 5mC/CpG ratio and GC is a negative one because the increase in methylation lags behind that of CpG, to reach a plateau in the GC-richest, gene-richest isochore families H2 and H3. In conclusion, there are more CpG targets that remain unmethylated in the GC-richest, gene-richest isochores in comparison with the other isochores. This conclusion supports the idea that the widespread methylation under consideration here has a general inhibitory effect on gene expression.     

Active mammalian replication origins are associated with a high-density cluster of mCpG dinucleotides.

ori-beta is a well-characterized origin of bidirectional replication (OBR) located approximately 17 kb downstream of the dihydrofolate reductase gene in hamster cell chromosomes. The approximately 2-kb region of ori-beta that exhibits greatest replication initiation activity also contains 12 potential methylation sites in the form of CpG dinucleotides. To ascertain whether DNA methylation might play a role at mammalian replication origins, the methylation status of these sites was examined with bisulfite to chemically distinguish cytosine (C) from 5-methylcytosine (mC). All of the CpGs were methylated, and nine of them were located within 356 bp flanking the minimal OBR, creating a high-density cluster of mCpGs that was approximately 10 times greater than average for human DNA. However, the previously reported densely methylated island in which all cytosines were methylated regardless of their dinucleotide composition was not detected and appeared to be an experimental artifact. A second OBR, located at the 5' end of the RPS14 gene, exhibited a strikingly similar methylation pattern, and the organization of CpG dinucleotides at other mammalian origins revealed the potential for high-density CpG methylation. Moreover, analysis of bromodeoxyuridine-labeled nascent DNA confirmed that active replication origins were methylated. These results suggest that a high-density cluster of mCpG dinucleotides may play a role in either the establishment or the regulation of mammalian replication origins.     

Age-related epigenome-wide DNA methylation and hydroxymethylation in longitudinal mouse blood

DNA methylation at cytosine-phosphate-guanine (CpG) dinucleotides changes as a function of age in humans and animal models, a process that may contribute to chronic disease development. Recent studies have investigated the role of an oxidized form of DNA methylation – 5-hydroxymethylcytosine (5hmC) – in the epigenome, but its contribution to age-related DNA methylation remains unclear. We tested the hypothesis that 5hmC changes with age, but in a direction opposite to 5-methylcytosine (5mC), potentially playing a distinct role in aging. To characterize epigenetic aging, genome-wide 5mC and 5hmC were measured in longitudinal blood samples (2, 4, and 10 months of age) from isogenic mice using two sequencing methods – enhanced reduced representation bisulfite sequencing and hydroxymethylated DNA immunoprecipitation sequencing. Examining the epigenome by age, we identified 28,196 unique differentially methylated CpGs (DMCs) and 8,613 differentially hydroxymethylated regions (DHMRs). Mouse blood showed a general pattern of epigenome-wide hypermethylation and hypo-hydroxymethylation with age. Comparing age-related DMCs and DHMRs, 1,854 annotated genes showed both differential 5mC and 5hmC, including one gene – Nfic – at five CpGs in the same 250 bp chromosomal region. At this region, 5mC and 5hmC levels both decreased with age. Reflecting these age-related epigenetic changes, Nfic RNA expression in blood decreased with age, suggesting that age-related regulation of this gene may be driven by 5hmC, not canonical DNA methylation. Combined, our genome-wide results show age-related differential 5mC and 5hmC, as well as some evidence that changes in 5hmC may drive age-related DNA methylation and gene expression.     

Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine

DNA cytosine-5 methylation is a well-studied epigenetic pathway implicated in gene expression control and disease pathogenesis. Different technologies have been developed to examine the distribution of 5-methylcytosine (5mC) in specific sequences of the genome. Recently, substantial amounts of 5-hydroxymethylcytosine (5hmC), most likely derived from enzymatic oxidation of 5mC by TET1, have been detected in certain mammalian tissues. Here, we have examined the ability of several commonly used DNA methylation profiling methods to distinguish between 5mC and 5hmC. We show that techniques based on sodium bisulfite treatment of DNA are incapable of distinguishing between the two modified bases. In contrast, techniques based on immunoprecipitation with anti-5mC antibody (methylated DNA immunoprecipitation, MeDIP) or those based on proteins that bind to methylated CpG sequences (e.g. methylated-CpG island recovery assay, MIRA) do not detect 5hmC and are specific for 5mC unless both modified bases occur in the same DNA fragment. We also report that several methyl-CpG binding proteins including MBD1, MBD2 and MBD4 do not bind to sequences containing 5hmC. Selective mapping of 5hmC will require the development of unique tools for the detection of this modified base.     

CpG mutation rates in the human genome are highly dependent on local GC content

CpG dinucleotides mutate at a high rate because cytosine is vulnerable to deamination, cytosines in CpG dinucleotides are often methylated, and deamination of 5-methylcytosine (5mC) produces thymidine. Previous experiments have shown that DNA melting is the rate-limiting step in cytosine deamination. Here we show, through the analysis of human single-nucleotide polymorphisms (SNPs), that the mutation rate produced by 5mC deamination is highly dependent on local GC content. In fact, linear regression analysis showed that the log(10) of the 5mC mutation rates (inferred from SNP frequencies) had slopes of -3 when graphed with respect to the GC content of neighboring sequences. This is the ideal slope that would be expected if the correlation between CpG underrepresentation and GC content had been solely caused by DNA melting. Moreover, this same result was obtained regardless of the SNP locations (all SNPs versus only SNPs in noncoding intergenic regions, excluding CpG islands) and regardless of the lengths over which GC content was calculated (SNP sequences with a modal length of 564 bp versus genomic contigs with a modal length of 163 kb). Several alternative interpretations are discussed.     

Role of base excision repair in maintaining the genetic and epigenetic integrity of CpG sites

Cytosine methylation at CpG dinucleotides is a central component of epigenetic regulation in vertebrates, and the base excision repair (BER) pathway is important for maintaining both the genetic stability and the methylation status of CpG sites. This perspective focuses on two enzymes that are of particular importance for the genetic and epigenetic integrity of CpG sites, methyl binding domain 4 (MBD4) and thymine DNA glycosylase (TDG). We discuss their capacity for countering C to T mutations at CpG sites, by initiating base excision repair of G·T mismatches generated by deamination of 5-methylcytosine (5mC). We also consider their role in active DNA demethylation, including pathways that are initiated by oxidation and/or deamination of 5mC.     

The majority of methylated deoxycytidines in human DNA are not in the CpG dinucleotide

The only natural postsynthetic modification known to occur in mammalian DNA is the methylation in the 5 position of deoxycytidines. Of the four 5'-CpN-3' dinucleotides (ie. CpG, CpC, CpA, and CpT), the dinucleotide which contains the highest proportion of deoxycytidines methylated is CpG, with 40 to 80% methylation in different mammalian genomes. It has also been shown that CpA, CpT, and CpC are methylated as well but to a much lower extent. Here we report the result of a full nearest neighbour analysis (together with quantitation of methylation levels in the 4 CpN dinucleotides) for DNA from human spleen. Using the values we have calculated the overall frequencies for all the methylated dinucleotides in the human genome. Because of the relative underrepresentation (by 7 to 10 fold) of the CpG dinucleotide, only 45.5% of total mC was present in mCpG, with 54.5% in mCpA, mCpT plus mCpC. These calculations have implications for studies into the function and significance of DNA methylation in mammalian cells.     

MIRA-SNuPE,a quantitative,multiplex method for measuring allele-specific DNA methylation

5-methyl-C (5mC) and 5-hydroxymethyl-C (5hmC) are epigenetic marks with well-known and putative roles in gene regulation, respectively. These two DNA covalent modifications cannot be distinguished by bisulfite sequencing or restriction digestion, the standard methods of 5mC detection. The methylated CpG island recovery assay (MIRA), however, specifically detects 5mC but not 5hmC. We further developed MIRA for the analysis of allele-specific CpG methylation at differentially methylated regions (DMRs) of imprinted genes. MIRA specifically distinguished between the parental alleles by capturing the paternally methylated H19/Igf2 DMR and maternally methylated KvDMR1 in mouse embryo fibroblasts (MEFs) carrying paternal and maternal duplication of mouse distal Chr7, respectively. MIRA in combination with multiplex single nucleotide primer extension (SNuPE) assays specifically captured the methylated parental allele from normal cells at a set of maternally and paternally methylated DMRs. The assay correctly recognized aberrant biallelic methylation in a case of loss of imprinting. The MIRA-SNuPE assays revealed that placenta exhibited less DNA methylation bias at DMRs compared to yolk sac, amnion, brain, heart, kidney, liver and muscle. This method should be useful for the analysis of allele-specific methylation events related to genomic imprinting, X chromosome inactivation and for verifying and screening haplotype-associated methylation differences in the human population.Key words: epigenetics, imprinting, DMR, MIRA, MBD, DNA methylation, SNuPE     

Genome-wide regulation of 5hmC, 5mC, and gene expression by Tet1 hydroxylase in mouse embryonic stem cells

DNA methylation at the 5 position of cytosine (5mC) in the mammalian genome is a key epigenetic event critical for various cellular processes. The ten-eleven translocation (Tet) family of 5mC-hydroxylases, which convert 5mC to 5-hydroxymethylcytosine (5hmC), offers a way for dynamic regulation of DNA methylation. Here we report that Tet1 binds to unmodified C or 5mC- or 5hmC-modified CpG-rich DNA through its CXXC domain. Genome-wide mapping of Tet1 and 5hmC reveals mechanisms by which Tet1 controls 5hmC and 5mC levels in mouse embryonic stem cells (mESCs). We also uncover a comprehensive gene network influenced by Tet1. Collectively, our data suggest that Tet1 controls DNA methylation both by binding to CpG-rich regions to prevent unwanted DNA methyltransferase activity, and by converting 5mC to 5hmC through hydroxylase activity. This Tet1-mediated antagonism of CpG methylation imparts differential maintenance of DNA methylation status at Tet1 targets, ultimately contributing to mESC differentiation and the onset of embryonic development.     

Different measures of “genome-wide” DNA methylation exhibit unique properties in placental and somatic tissues

DNA methylation of CpGs located in two types of repetitive elements—LINE1 (L1) and Alu—is used to assess “global” changes in DNA methylation in studies of human disease and environmental exposure. L1 and Alu contribute close to 30% of all base pairs in the human genome and transposition of repetitive elements is repressed through DNA methylation. Few studies have investigated whether repetitive element DNA methylation is associated with DNA methylation at other genomic regions, or the biological and technical factors that influence potential associations. Here, we assess L1 and Alu DNA methylation by Pyrosequencing of consensus sequences and using subsets of probes included in the Illumina Infinium HumanMethylation27 BeadChip array. We show that evolutionary age and assay method affect the assessment of repetitive element DNA methylation. Additionally, we compare Pyrosequencing results for repetitive elements to average DNA methylation of CpG islands, as assessed by array probes classified into strong, weak and non-islands. We demonstrate that each of these dispersed sequences exhibits different patterns of tissue-specific DNA methylation. Correlation of DNA methylation suggests an association between L1 and weak CpG island DNA methylation in some of the tissues examined. We caution, however, that L1, Alu and CpG island DNA methylation are distinct measures of dispersed DNA methylation and one should not be used in lieu of another. Analysis of DNA methylation data is complex and assays may be influenced by environment and pathology in different or complementary ways.     

Different measures of “genome-wide” DNA methylation exhibit unique properties in placental and somatic tissues

DNA methylation of CpGs located in two types of repetitive elements—LINE1 (L1) and Alu—is used to assess “global” changes in DNA methylation in studies of human disease and environmental exposure. L1 and Alu contribute close to 30% of all base pairs in the human genome and transposition of repetitive elements is repressed through DNA methylation. Few studies have investigated whether repetitive element DNA methylation is associated with DNA methylation at other genomic regions, or the biological and technical factors that influence potential associations. Here, we assess L1 and Alu DNA methylation by Pyrosequencing of consensus sequences and using subsets of probes included in the Illumina Infinium HumanMethylation27 BeadChip array. We show that evolutionary age and assay method affect the assessment of repetitive element DNA methylation. Additionally, we compare Pyrosequencing results for repetitive elements to average DNA methylation of CpG islands, as assessed by array probes classified into strong, weak and non-islands. We demonstrate that each of these dispersed sequences exhibits different patterns of tissue-specific DNA methylation. Correlation of DNA methylation suggests an association between L1 and weak CpG island DNA methylation in some of the tissues examined. We caution, however, that L1, Alu and CpG island DNA methylation are distinct measures of dispersed DNA methylation and one should not be used in lieu of another. Analysis of DNA methylation data is complex and assays may be influenced by environment and pathology in different or complementary ways.     

奶牛产奶性状与乳房炎相关基因CpG含量及分布特征的比较分析

DNA甲基化是表观遗传学的重要组成部分,基因启动子区及第一外显子区的CpG甲基化通常抑制该基因的表达,而去甲基化则促进基因表达。已有的研究发现荷斯坦牛的乳房炎指标SCC(Somatic cell count)与产奶量呈较强负相关。文章分析并比较了这两类性状的相关基因的启动子区、第一外显子、下游2 000 bp序列中CpG含量及分布特征。结果表明,乳房炎相关基因的启动子、第一外显子中CpG含量显著低于产奶性状相关基因,而两类性状基因下游2 000 bp序列中CpG含量无显著性差异。另外,文中提出了两个量化基因序列中CpG特征的指标,一个是CpG平均距离,用来衡量序列中的CpG分布;另一个是条件概率p(G|C),用以量化序列中二核苷酸CpG随碱基C出现的可能性,并对两类基因的启动子和第一外显子区域的这两个指标做了统计检验。研究结果对产奶性状与乳房炎相关基因的DNA甲基化调控研究奠定了基础。     

In vivo control of CpG and non-CpG DNA methylation by DNA methyltransferases

The enzymatic control of the setting and maintenance of symmetric and non-symmetric DNA methylation patterns in a particular genome context is not well understood. Here, we describe a comprehensive analysis of DNA methylation patterns generated by high resolution sequencing of hairpin-bisulfite amplicons of selected single copy genes and repetitive elements (LINE1, B1, IAP-LTR-retrotransposons, and major satellites). The analysis unambiguously identifies a substantial amount of regional incomplete methylation maintenance, i.e. hemimethylated CpG positions, with variant degrees among cell types. Moreover, non-CpG cytosine methylation is confined to ESCs and exclusively catalysed by Dnmt3a and Dnmt3b. This sequence position-, cell type-, and region-dependent non-CpG methylation is strongly linked to neighboring CpG methylation and requires the presence of Dnmt3L. The generation of a comprehensive data set of 146,000 CpG dyads was used to apply and develop parameter estimated hidden Markov models (HMM) to calculate the relative contribution of DNA methyltransferases (Dnmts) for de novo and maintenance DNA methylation. The comparative modelling included wild-type ESCs and mutant ESCs deficient for Dnmt1, Dnmt3a, Dnmt3b, or Dnmt3a/3b, respectively. The HMM analysis identifies a considerable de novo methylation activity for Dnmt1 at certain repetitive elements and single copy sequences. Dnmt3a and Dnmt3b contribute de novo function. However, both enzymes are also essential to maintain symmetrical CpG methylation at distinct repetitive and single copy sequences in ESCs.     

Germline methylation patterns inferred from local nucleotide frequency of repetitive sequences in the human genome

Given the genomic abundance and susceptibility to DNA methylation, interspersed repetitive sequences in the human genome can be exploited as valuable resources in genome-wide methylation studies. To learn about the relationships between DNA methylation and repeat sequences, we performed a global measurement of CpG dinucleotide frequencies for interspersed repetitive sequences and inferred germline methylation patterns in the human genome. Although extensive CpG depletion was observed for most repeat sequences, those in the proximity to CpG islands have been relatively removed from germline methylation being the potential source of germline activation. We also investigated the CpG depletion patterns of Alu pairs to see whether they might play an active role in germline methylation. Two kinds of Alu pairs, direct or inverted pairs classified according to the orientation, showed contrast CpG depletion patterns with respect to separating distance of Alus, i.e., as two Alu elements are more closely spaced in a pair, a higher extent of CpG depletion was observed in inverted orientation and vice versa for directly repetitive Alu pairs. This suggests that specific organization of repetitive sequences, such as inverted Alu pairs, might play a role in triggering DNA methylation consistent with a homology-dependent methylation hypothesis.     

Methyl CpG Binding Domain Ultra‐Sequencing: a novel method for identifying inter‐individual and cell‐type‐specific variation in DNA methylation

Experience‐dependent changes in DNA methylation can exert profound effects on neuronal function and behaviour. A single learning event can induce a variety of DNA modifications within the neuronal genome, some of which may be common to all individuals experiencing the event, whereas others may occur in a subset of individuals. Variations in experience‐induced DNA methylation may subsequently confer increased vulnerability or resilience to the development of neuropsychiatric disorders. However, the detection of experience‐dependent changes in DNA methylation in the brain has been hindered by the interrogation of heterogeneous cell populations, regional differences in epigenetic states and the use of pooled tissue obtained from multiple individuals. Methyl CpG Binding Domain Ultra‐Sequencing (MBD Ultra‐Seq) overcomes current limitations on genome‐wide epigenetic profiling by incorporating fluorescence‐activated cell sorting and sample‐specific barcoding to examine cell‐type‐specific CpG methylation in discrete brain regions of individuals. We demonstrate the value of this method by characterizing differences in 5‐methylcytosine (5mC) in neurons and non‐neurons of the ventromedial prefrontal cortex of individual adult C57BL/6 mice, using as little as 50 ng of genomic DNA per sample. We find that the neuronal methylome is characterized by greater CpG methylation as well as the enrichment of 5mC within intergenic loci. In conclusion, MBD Ultra‐Seq is a robust method for detecting DNA methylation in neurons derived from discrete brain regions of individual animals. This protocol will facilitate the detection of experience‐dependent changes in DNA methylation in a variety of behavioural paradigms and help identify aberrant experience‐induced DNA methylation that may underlie risk and resiliency to neuropsychiatric disease.     

Methylated DNA Immunoprecipitation

The identification of DNA methylation patterns is a common procedure in the study of epigenetics, as methylation is known to have significant effects on gene expression, and is involved with normal development as well as disease ^1-4. Thus, the ability to discriminate between methylated DNA and non-methylated DNA is essential for generating methylation profiles for such studies. Methylated DNA immunoprecipitation (MeDIP) is an efficient technique for the extraction of methylated DNA from a sample of interest ^5-7. A sample of as little as 200 ng of DNA is sufficient for the antibody, or immunoprecipitation (IP), reaction. DNA is sonicated into fragments ranging in size from 300-1000 bp, and is divided into immunoprecipitated (IP) and input (IN) portions. IP DNA is subsequently heat denatured and then incubated with anti-5''mC, allowing the monoclonal antibody to bind methylated DNA. After this, magnetic beads containing a secondary antibody with affinity for the primary antibody are added, and incubated. These bead-linked antibodies will bind the monoclonal antibody used in the first step. DNA bound to the antibody complex (methylated DNA) is separated from the rest of the DNA by using a magnet to pull the complexes out of solution. Several washes using IP buffer are then performed to remove the unbound, non-methylated DNA. The methylated DNA/antibody complexes are then digested with Proteinase K to digest the antibodies leaving only the methylated DNA intact. The enriched DNA is purified by phenol:chloroform extraction to remove the protein matter and then precipitated and resuspended in water for later use. PCR techniques can be used to validate the efficiency of the MeDIP procedure by analyzing the amplification products of IP and IN DNA for regions known to lack and known to contain methylated sequences. The purified methylated DNA can then be used for locus-specific (PCR) or genome-wide (microarray and sequencing) methylation studies, and is particularly useful when applied in conjunction with other research tools such as gene expression profiling and array comparative genome hybridization (CGH) ⁸. Further investigation into DNA methylation will lead to the discovery of new epigenetic targets, which in turn, may be useful in developing new therapeutic or prognostic research tools for diseases such as cancer that are characterized by aberrantly methylated DNA ^{2, 4, 9-11}.