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.     

Genome‐wide DNA methylation analysis: no evidence for stable hemimethylation in the sheep muscle genome

The importance of maintaining DNA methylation patterns and faithful transmission of these patterns during cell division to ensure appropriate gene expression has been known for many decades now. It has largely been assumed that the symmetrical nature of CpG motifs, the most common site for DNA methylation in mammals, together with the presence of maintenance methylases able to methylate newly synthesised DNA, ensures that there is concordance of methylation on both strands. However, although this assumption is compelling in theory, little experimental evidence exists that either supports or refutes this assumption. Here, we have undertaken a genome‐wide single‐nucleotide resolution analysis to determine the frequency with which hemimethylated CpG sites exist in sheep muscle tissue. Analysis of multiple independent samples provides strong evidence that stably maintained hemimethylation is a very rare occurrence, at least in this tissue. Given the rarity of stably maintained hemimethylation, next‐generation sequencing data from both DNA strands may be carefully combined to increase the accuracy with which DNA methylation can be measured at single‐nucleotide resolution.     

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.     

Both hypomethylation and hypermethylation in a 0.2-kb region of a DNA repeat in cancer

NBL2 is a tandem 1.4-kb DNA repeat, whose hypomethylation in hepatocellular carcinomas was shown previously to be an independent predictor of disease progression. Here, we examined methylation of all cytosine residues in a 0.2-kb subregion of NBL2 in ovarian carcinomas, Wilms' tumors, and diverse control tissues by hairpin-bisulfite PCR. This new genomic sequencing method detects 5-methylcytosine on covalently linked complementary strands of a DNA fragment. All DNA clones from normal somatic tissues displayed symmetrical methylation at seven CpG positions and no methylation or only hemimethylation at two others. Unexpectedly, 56% of cancer DNA clones had decreased methylation at some normally methylated CpG sites as well as increased methylation at one or both of the normally unmethylated sites. All 146 DNA clones from 10 cancers could be distinguished from all 91 somatic control clones by assessing methylation changes at three of these CpG sites. The special involvement of DNA methyltransferase 3B in NBL2 methylation was indicated by analysis of cells from immunodeficiency, centromeric region instability, and facial anomalies syndrome patients who have mutations in the gene encoding DNA methyltransferase 3B. Blot hybridization of 33 cancer DNAs digested with CpG methylation-sensitive enzymes confirmed that NBL2 arrays are unusually susceptible to cancer-linked hypermethylation and hypomethylation, consistent with our novel genomic sequencing findings. The combined Southern blot and genomic sequencing data indicate that some of the cancer-linked alterations in CpG methylation are occurring with considerable sequence specificity. NBL2 is an attractive candidate for an epigenetic cancer marker and for elucidating the nature of epigenetic changes in cancer.     

Sensitive quantitative analysis of murine LINE1 DNA methylation using high resolution melt analysis

We present here the first high resolution melt (HRM) assay to quantitatively analyze differences in murine DNA methylation levels utilizing CpG methylation of Long Interspersed Elements-1 (LINE1 or L1). By calculating the integral difference in melt temperature between samples and a methylated control, and biasing PCR primers for unmethylated CpGs, the assay demonstrates enhanced sensitivity to detect changes in methylation in a cell line treated with low doses of 5-aza-2’-deoxycytidine (5-aza). The L1 assay was confirmed to be a good marker of changes in DNA methylation of L1 elements at multiple regions across the genome when compared with total 5-methyl-cytosine content, measured by Liquid Chromatography-Mass Spectrometry (LC-MS). The assay design was also used to detect changes in methylation at other murine repeat elements (B1 and Intracisternal-A-particle Long-terminal Repeat elements). Pyrosequencing analysis revealed that L1 methylation changes were non-uniform across the CpGs within the L1-HRM target region, demonstrating that the L1 assay can detect small changes in CpG methylation among a large pool of heterogeneously methylated DNA templates. Application of the assay to various tissues from Balb/c and CBA mice, including previously unreported peripheral blood (PB), revealed a tissue hierarchy (from hypermethylated to hypomethylated) of PB > kidney > liver > prostate > spleen. CBA mice demonstrated overall greater methylation than Balb/c mice, and male mice demonstrated higher tissue methylation compared with female mice in both strains. Changes in DNA methylation have been reported to be an early and fundamental event in the pathogenesis of many human diseases, including cancer. Mouse studies designed to identify modulators of DNA methylation, the critical doses, relevant time points and the tissues affected are limited by the low throughput nature and exorbitant cost of many DNA methylation assays. The L1 assay provides a high throughput, inexpensive and sensitive screening tool for identifying and characterizing DNA methylation changes to L1 elements at multiple regions across the genome.     

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.     

L1 repeat elements in the human epsilon-G gamma-globin gene intergenic region: sequence analysis and concerted evolution within this family

We have deduced the sequence of a composite long interspersed repeated DNA in primates and herein describe its relationship to a complex repeat element (L1Heg) located in the interval linking the human epsilon- and G gamma-globin genes. The main element of L1Heg is 3' truncated and interrupted by the insertion of the 3' end of a second L1 element. Transposition of L1Heg into this intergenic locus generated a 62-bp duplication of flanking sequences. In contrast, insertion of the second repeat may have been mediated by homology between donor and target sequences. The main repeat represents a novel class of abundant elements whose sequences have diverged from other rodent and primate LINES approximately 1.3 kb downstream from the 5' terminus of L1Heg. Comparison of L1Heg with the sequences of two other related L1 members revealed a complex set of rearrangements confined within a region that resembles the long terminal repeats of other types of retroposons. The boundaries of conversion-like events were defined on the basis of the clustering of nucleotide sequence variants common to two or more nonallelic 3' L1H elements. Several of these events are apparently initiated or resolved within a common 150-bp region that coincides with the 3' terminus of a pan-mammalian open reading frame. This analysis showed that concerted genetic interactions and random drift both contribute appreciably to sequence variation within this set of L1H members.     

A comparison of existing global DNA methylation assays to low-coverage whole-genome bisulfite sequencing for epidemiological studies

DNA methylation is an epigenetic mark at the interface of genetic and environmental factors relevant to human disease. Quantitative assessments of global DNA methylation levels have therefore become important tools in epidemiology research, particularly for understanding effects of environmental exposures in complex diseases. Among the available methods of quantitative DNA methylation measurements, bisulfite sequencing is considered the gold standard, but whole-genome bisulfite sequencing (WGBS) has previously been considered too costly for epidemiology studies with high sample numbers. Pyrosequencing of repetitive sequences within bisulfite-treated DNA has been routinely used as a surrogate for global DNA methylation, but a comparison of pyrosequencing to WGBS for accuracy and reproducibility of methylation levels has not been performed. This study compared the global methylation levels measured from uniquely mappable (non-repetitive) WGBS sequences to pyrosequencing assays of several repeat sequences and repeat assay-matched WGBS data and determined uniquely mappable WGBS data to be the most reproducible and accurate measurement of global DNA methylation levels. We determined sources of variation in repetitive pyrosequencing assays to be PCR amplification bias, PCR primer selection bias in methylation levels of targeted sequences, and inherent variability in methylation levels of repeat sequences. Low-coverage, uniquely mappable WGBS showed the strongest correlation between replicates of all assays. By using multiplexing by indexed bar codes, the cost of WGBS can be lowered significantly to improve the accuracy of global DNA methylation assessments for human studies.     

DNA methylation polymorphisms precede any histological sign of atherosclerosis in mice lacking apolipoprotein E

The present work investigates the occurrence and significance of aberrant DNA methylation patterns during early stages of atherosclerosis. To this end, we asked whether the genetically atherosclerosis-prone APOE-null mice show any changes in DNA methylation patterns before the appearance of histologically detectable vascular lesion. We exploited a combination of various techniques: DNA fingerprinting, in vitro methyl-accepting assay, 5-methylcytosine quantitation, histone post-translational modification analysis, Southern blotting, and PCR. Our results show that alterations in DNA methylation profiles, including both hyper- and hypomethylation, were present in aortas and PBMC of 4-week-old mutant mice with no detectable atherosclerotic lesion. Sequencing and expression analysis of 60 leukocytic polymorphisms revealed that epigenetic changes involve transcribed genic sequences, as well as repeated interspersed elements. Furthermore, we showed for the first time that atherogenic lipoproteins promote global DNA hypermethylation in a human monocyte cell line. Taken together, our results unequivocally show that alterations in DNA methylation profiles are early markers of atherosclerosis in a mouse model and may play a causative role in atherogenesis.     

Cooperativity between DNA methyltransferases in the maintenance methylation of repetitive elements.

We used mouse embryonic stem (ES) cells with systematic gene knockouts for DNA methyltransferases to delineate the roles of DNA methyltransferase 1 (Dnmt1) and Dnmt3a and -3b in maintaining methylation patterns in the mouse genome. Dnmt1 alone was able to maintain methylation of most CpG-poor regions analyzed. In contrast, both Dnmt1 and Dnmt3a and/or Dnmt3b were required for methylation of a select class of sequences which included abundant murine LINE-1 promoters. We used a novel hemimethylation assay to show that even in wild-type cells these sequences contain high levels of hemimethylated DNA, suggestive of poor maintenance methylation. We showed that Dnmt3a and/or -3b could restore methylation of these sequences to pretreatment levels following transient exposure of cells to 5-aza-CdR, whereas Dnmt1 by itself could not. We conclude that ongoing de novo methylation by Dnmt3a and/or Dnmt3b compensates for inefficient maintenance methylation by Dnmt1 of these endogenous repetitive sequences. Our results reveal a previously unrecognized degree of cooperativity among mammalian DNA methyltransferases in ES cells.     

Deep Enzymology Studies on DNA Methyltransferases Reveal Novel Connections between Flanking Sequences and Enzyme Activity

DNA interacting enzymes recognize their target sequences embedded in variable flanking sequence context. The influence of flanking sequences on enzymatic activities of DNA methyltransferases (DNMTs) can be systematically studied with “deep enzymology” approaches using pools of double-stranded DNA substrates, which contain target sites in random flanking sequence context. After incubation with DNMTs and bisulfite conversion, the methylation states and flanking sequences of individual DNA molecules are determined by NGS. Deep enzymology studies with different human and mouse DNMTs revealed strong influences of flanking sequences on their CpG and non-CpG methylation activity and the structures of DNMT-DNA complexes. Differences in flanking sequence preferences of DNMT3A and DNMT3B were shown to be related to the prominent role of DNMT3B in the methylation of human SATII repeat elements. Mutational studies in DNMT3B discovered alternative interaction networks between the enzyme and the DNA leading to a partial equalization of the effects of different flanking sequences. Structural studies in DNMT1 revealed striking correlations between enzymatic activities and flanking sequence dependent conformational changes upon DNA binding. Correlation of the biochemical data with cellular methylation patterns demonstrated that flanking sequence preferences are an important parameter that influences genomic DNA methylation patterns together with other mechanisms targeting DNMTs to genomic sites.     

Advances in mechanisms of genetic instability related to hereditary neurological diseases

Substantial progress has been realized in the past several years in our understanding of the molecular mechanisms responsible for the expansions and deletions (genetic instabilities) of repeating tri-, tetra- and pentanucleotide repeating sequences associated with a number of hereditary neurological diseases. These instabilities occur by replication, recombination and repair processes, probably acting in concert, due to slippage of the DNA complementary strands relative to each other. The biophysical properties of the folded-back repeating sequence strands play a critical role in these instabilities. Non-B DNA structural elements (hairpins and slipped structures, DNA unwinding elements, tetraplexes, triplexes and sticky DNA) are described. The replication mechanisms are influenced by pausing of the replication fork, orientation of the repeat strands, location of the repeat sequences relative to replication origins and the flap endonuclease. Methyl-directed mismatch repair, nucleotide excision repair, and repair of damage caused by mutagens are discussed. Genetic recombination and double-strand break repair advances in Escherichia coli, yeast and mammalian models are reviewed. Furthermore, the newly discovered capacities of certain triplet repeat sequences to cause gross chromosomal rearrangements are discussed.     

Methylation pattern of Activator transposase binding sites in maize endosperm.

The maize transposable element Activator (Ac) transposes after replication from only one of the two daughter chromatids. It has been suggested that DNA methylation in conjunction with methylation-sensitive transposase binding to DNA may control the association of Ac transposition and replication. We present here a detailed genomic sequencing analysis of the cytosine methylation patterns of the transposase binding sites within both Ac ends in the wx-m9::Ac allele, where Ac is inserted into the tenth exon of the Waxy gene. The Ac elements in wx-m9::Ac kernels exhibit intriguing methylation patterns and fall into two distinct groups. Approximately 50% of the elements are fully unmethylated at cytosine residues through the 256 nucleotides at the 5' end (the promoter end). The other half is partially methylated between Ac residues 27 and 92. In contrast, at the 3' end, all Ac molecules are heavily methylated between residues 4372 and 4554. The more internally located Ac sequences and the flanking Waxy DNA are unmethylated. Although most methylated cytosines in Ac are in the symmetrical CpG and CpNpG arrangements, nonsymmetrical cytosine methylation is also common in the hypermethylated regions of Ac. These results suggest a model in which differential activation of transposon ends by hemimethylation controls the chromatid selectivity of transposition and the association with replication.     

Characterizing the strand-specific distribution of non-CpG methylation in human pluripotent cells

DNA methylation is an important defense and regulatory mechanism. In mammals, most DNA methylation occurs at CpG sites, and asymmetric non-CpG methylation has only been detected at appreciable levels in a few cell types. We are the first to systematically study the strand-specific distribution of non-CpG methylation. With the divide-and-compare strategy, we show that CHG and CHH methylation are not intrinsically different in human embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs). We also find that non-CpG methylation is skewed between the two strands in introns, especially at intron boundaries and in highly expressed genes. Controlling for the proximal sequences of non-CpG sites, we show that the skew of non-CpG methylation in introns is mainly guided by sequence skew. By studying subgroups of transposable elements, we also found that non-CpG methylation is distributed in a strand-specific manner in both short interspersed nuclear elements (SINE) and long interspersed nuclear elements (LINE), but not in long terminal repeats (LTR). Finally, we show that on the antisense strand of Alus, a non-CpG site just downstream of the A-box is highly methylated. Together, the divide-and-compare strategy leads us to identify regions with strand-specific distributions of non-CpG methylation in humans.     

Distribution of Alu and L1 repeats in human YAC recombinants

Evidence is accumulating that the two major families of interspersed repeated human DNA sequences, Alu and L1, are not randomly distributed. However, only limited information is available on their relative long-range distribution. We have analyzed a set of randomly selected, human Chromosome (Chr) 11-specific YAC recombinants constituting a total length of about 2 Mbp for the local and global distribution of Alu and L1 repeats: the data show a strong asymmetry in the distribution of these two repeat classes and give weight, at the long-range molecular level, to previous studies indicating their partition in the human genome; they also suggest a strong tendency for L1 repeats to cluster, with a higher proportion of full-length elements than expected.     

Human DNA mismatch repair in vitro operates independently of methylation status at CpG sites

Whereas in Escherichia coli DNA mismatch repair is directed to the newly synthesized strand due to its transient lack of adenine methylation, the molecular determinants of strand discrimination in eukaryotes are presently unknown. In mammalian cells, cytosine methylation within CpG sites may represent an analogous and mechanistically plausible means of targeting mismatch correction. Using HeLa nuclear extracts, we conducted a systematic analysis in vitro to determine whether cytosine methylation participates in human DNA mismatch repair. We prepared a set of A·C heteroduplex molecules that were either unmethylated, hemimethylated or fully methylated at CpG sequences and found that the methylation status persisted under the assay conditions. However, no effect on either the time course or the magnitude of mismatch repair events was evident; only strand discontinuities contributed to strand bias. By western analysis we demonstrated that the HeLa extract contained MED1 protein, which interacts with MLH1 and binds to CpG-methylated DNA; supplementation with purified MED1 protein was without effect. In summary, human DNA mismatch repair operates independently of CpG methylation status, and we found no evidence supporting a role for CpG hemimethylation as a strand discrimination signal.     

Human repeat element-mediated PCR: cloning and mapping of chromosome 10 DNA markers.

Repeat element-mediated PCR can facilitate rapid cloning and mapping of human chromosomal region-specific DNA markers from somatic cell hybrid DNA. PCR primers directed to human repeat elements result in human-specific DNA synthesis; template DNA derived from a somatic cell hybrid containing the human chromosomal region of interest provides region specificity. We have generated a series of repeat element-mediated PCR clones from a reduced complexity somatic cell hybrid containing a portion of human chromosome 10. The cloning source retains the centromere and tightly linked flanking markers, plus additional chromosome 10 sequences. Twelve new inter-Alu, two inter-L1, and four inter-Alu/L1 repeat element-mediated PCR clones were mapped by hybridization to Southern blots of repeat element-mediated PCR products amplified from somatic cell hybrid DNA templates. Two inter-Alu clones mapped to the pericentromeric region. We propose that a scarcity of Alu elements in the pericentromeric region of chromosome 10 contributed to the low number of clones obtained from this region. One inter-Alu clone, pC11/A1S-6-c23, defines the D10S94 locus, which is tightly linked to MEN2A and D10Z1.     

Methylation dynamics of repetitive DNA elements in the mouse germ cell lineage

Repetitive DNA elements account for a substantial fraction of the mammalian genome. Many are subject to DNA methylation, which is known to undergo dynamic change during mouse germ cell development. We found that repeat sequences of three different classes retain high levels of methylation at E12.5, when methylation is erased from many single-copy genes. Maximal demethylation of repeats was seen later in development and at different times in male and female germ cells. At none of the time points examined (E12.5, E15.5, and E17.5) did we see complete demethylation, suggesting that methylation patterns on repeats may be passed on from one generation to the next. In male germ cells, we observed a de novo methylation event resulting in complete methylation of all the repeats in the interval between E15.5 and E17.5, which was not seen in females. These results suggest that repeat sequences undergo coordinate changes in methylation during germ cell development and give further insights into germ cell reprogramming in mice.