Constitutively methylated CpG dinucleotides as mutation hot spots in the retinoblastoma gene (RB1).

A wide spectrum of mutations, ranging from point mutations to large deletions, have been described in the retinoblastoma gene (RB1). Mutations have been found throughout the gene; however, these genetic alterations do not appear to be homogeneously distributed. In particular, a significant proportion of disease-causing mutations results in the premature termination of protein synthesis, and the majority of these mutations occur as C-->T transitions at CpG dinucleotides (CpGs). Such recurrent CpG mutations, including those found in RB1, are likely the result of the deamination of 5-methylcytosine within these CpGs. In the present study, we used the sodiumbisulfite conversion method to detect cytosine methylation in representative exons of RB1. We analyzed DNA from a variety of tissues and specifically targeted CGA codons in RB1, where recurrent premature termination mutations have been reported. We found that DNA methylation within RB1 exons 8, 14, 25, and 27 appeared to be restricted to CpGs, including six CGA codons. Other codons containing methylated cytosines have not been reported to be mutated. Therefore, disease-causing mutations at CpGs in RB1 appear to be determined by several factors, including the constitutive presence of DNA methylation at cytosines within CpGs, the specific codon within which the methylated cytosine is located, and the particular region of the gene within which that codon resides.     

An evaluation of analysis pipelines for DNA methylation profiling using the Illumina HumanMethylation450 BeadChip platform

The proper identification of differentially methylated CpGs is central in most epigenetic studies. The Illumina HumanMethylation450 BeadChip is widely used to quantify DNA methylation; nevertheless, the design of an appropriate analysis pipeline faces severe challenges due to the convolution of biological and technical variability and the presence of a signal bias between Infinium I and II probe design types. Despite recent attempts to investigate how to analyze DNA methylation data with such an array design, it has not been possible to perform a comprehensive comparison between different bioinformatics pipelines due to the lack of appropriate data sets having both large sample size and sufficient number of technical replicates. Here we perform such a comparative analysis, targeting the problems of reducing the technical variability, eliminating the probe design bias and reducing the batch effect by exploiting two unpublished data sets, which included technical replicates and were profiled for DNA methylation either on peripheral blood, monocytes or muscle biopsies. We evaluated the performance of different analysis pipelines and demonstrated that: (1) it is critical to correct for the probe design type, since the amplitude of the measured methylation change depends on the underlying chemistry; (2) the effect of different normalization schemes is mixed, and the most effective method in our hands were quantile normalization and Beta Mixture Quantile dilation (BMIQ); (3) it is beneficial to correct for batch effects. In conclusion, our comparative analysis using a comprehensive data set suggests an efficient pipeline for proper identification of differentially methylated CpGs using the Illumina 450K arrays.     

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.     

Developmental methylation of the coding region of c-fos occurs perinatally, stepwise and sequentially in the liver of laboratory mouse

We have studied the dynamics of de novo DNA methylation of 16 contiguous CpGs in the non-CpG island-coding region of the proto-oncogene c-fos during mouse development by Na-bisulfite sequencing. Methylation commences from 16.5 dpc and occurs in stepwise-manner. In liver 7 sites are methylated between 16.5 dpc and day 5 after birth, but all the sites are completely methylated on 20 dpp and remain so in the adult liver. The present study provides evidence that (1) pattern of methylation of c-fos is distinct from those DNA sequences which methylate pre- and post-implantation, both in terms of the timing and spreading, and (2) spacing of CpGs is an important factor in determining the course of methylation. We suggest that there could be other isoforms of Dnmtases for the c-fos like embryonic genes, not only because they methylate later in development but also because of the difference in kinetics of the reaction, and that the nucleation of certain methylated sites facilitate methylation of neighbouring sites and their maintenance in subsequent cell generations.     

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.     

Methylated site display (MSD)-AFLP,a sensitive and affordable method for analysis of CpG methylation profiles

Background

It has been pointed out that environmental factors or chemicals can cause diseases that are developmental in origin. To detect abnormal epigenetic alterations in DNA methylation, convenient and cost-effective methods are required for such research, in which multiple samples are processed simultaneously. We here present methylated site display (MSD), a unique technique for the preparation of DNA libraries. By combining it with amplified fragment length polymorphism (AFLP) analysis, we developed a new method, MSD-AFLP.

Results

Methylated site display libraries consist of only DNAs derived from DNA fragments that are CpG methylated at the 5′ end in the original genomic DNA sample. To test the effectiveness of this method, CpG methylation levels in liver, kidney, and hippocampal tissues of mice were compared to examine if MSD-AFLP can detect subtle differences in the levels of tissue-specific differentially methylated CpGs. As a result, many CpG sites suspected to be tissue-specific differentially methylated were detected. Nucleotide sequences adjacent to these methyl-CpG sites were identified and we determined the methylation level by methylation-sensitive restriction endonuclease (MSRE)-PCR analysis to confirm the accuracy of AFLP analysis. The differences of the methylation level among tissues were almost identical among these methods. By MSD-AFLP analysis, we detected many CpGs showing less than 5% statistically significant tissue-specific difference and less than 10% degree of variability. Additionally, MSD-AFLP analysis could be used to identify CpG methylation sites in other organisms including humans.

Conclusion

MSD-AFLP analysis can potentially be used to measure slight changes in CpG methylation level. Regarding the remarkable precision, sensitivity, and throughput of MSD-AFLP analysis studies, this method will be advantageous in a variety of epigenetics-based research.

     

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.     

Maternal cigarette smoking during pregnancy and offspring DNA methylation in midlife

Maternal smoking in pregnancy (MSP) has been associated with DNA methylation in specific CpG sites (CpGs) in infants and children. We investigated whether MSP, independent of own personal active smoking, was associated with midlife DNA methylation in CpGs that were previously identified in studies of MSP-DNA methylation in children. We used data on MSP collected from pregnant mothers of 89 adult women born in 1959–1964 and measured DNA methylation in blood (granulocytes) collected in 2001–2007 (mean age: 43 years). Seventeen CpGs were differentially methylated by MSP, with multiple CpGs mapping to CYP1A1, MYO1G, AHRR, and GFI1. These associations were consistent in direction with prior studies (e.g., MSP associated with more and less methylation in AHRR and CYP1A1, respectively) and, with the exception of AHRR CpGs, were not substantially altered by adjustment for active smoking. These preliminary results confirm prior prospective reports that MSP influences the offspring DNA methylation, and extends the timeframe to midlife, and suggest that these effects may persist into adulthood, independently of active smoking.     

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.     

DNA methylation at mammalian replication origins.

In Escherichia coli, DNA methylation regulates both origin usage and the time required to reassemble prereplication complexes at replication origins. In mammals, at least three replication origins are associated with a high density cluster of methylated CpG dinucleotides, and others whose methylation status has not yet been characterized have the potential to exhibit a similar DNA methylation pattern. One of these origins is found within the approximately 2-kilobase pair region upstream of the human c-myc gene that contains 86 CpGs. Application of the bisulfite method for detecting 5-methylcytosines at specific DNA sequences revealed that this region was not methylated in either total genomic DNA or newly synthesized DNA. Therefore, DNA methylation is not a universal component of mammalian replication origins. To determine whether or not DNA methylation plays a role in regulating the activity of origins that are methylated, the rate of remethylation and the effect of hypomethylation were determined at origin beta (ori-beta), downstream of the hamster DHFR gene. Remethylation at ori-beta did not begin until approximately 500 base pairs of DNA was synthesized, but it was then completed by the time that 4 kilobase pairs of DNA was synthesized (<3 min after release into S phase). Thus, DNA methylation cannot play a significant role in regulating reassembly of prereplication complexes in mammalian cells, as it does in E. coli. To determine whether or not DNA methylation plays any role in origin activity, hypomethylated hamster cells were examined for ori-beta activity. Cells that were >50% reduced in methylation at ori-beta no longer selectively activated ori-beta. Therefore, at some loci, DNA methylation either directly or indirectly determines where replication begins.     

Epigenome-wide profiling identifies significant differences in DNA methylation between matched-pairs of T- and B-lymphocytes from healthy individuals

Multiple reports now describe changes to the DNA methylome in rheumatoid arthritis and in many cases have analyzed methylation in mixed cell populations from whole blood. However, these approaches may preclude the identification of cell type-specific methylation, which may subsequently bias identification of disease-specific changes. To address this possibility, we conducted genome-wide DNA methylation profiling using HumanMethylation450 BeadChips to identify differences within matched pairs of T-lymphocytes and B-lymphocytes isolated from the peripheral blood of 10 healthy females. Array data were processed and differential methylation identified using NIMBL software. Validation of array data was performed by bisulfite pyrosequencing. Genome-wide DNA methylation was initially determined by analysis of LINE-1 sequences and was higher in B-lymphocytes than matched T-lymphocytes (69.8% vs. 65.2%, P ≤ 0.01). Pairwise analysis identified 679 CpGs, representing 250 genes, which were differentially methylated between T-lymphocytes and B-lymphocytes. The majority of sites (76.6%) were hypermethylated in B-lymphocytes. Pyrosequencing of selected candidates confirmed the array data in all cases. Hierarchical clustering revealed perfect segregation of samples into two distinct clusters based on cell type. Differentially methylated genes showed enrichment for biological functions/pathways associated with leukocytes and T-lymphocytes. Our work for the first time shows that T-lymphocytes and B-lymphocytes possess intrinsic differences in DNA methylation within a restricted set of functionally related genes. These data provide a foundation for investigating DNA methylation in diseases in which these cell types play important and distinct roles.     

Methylation, mutation and cancer.

The fifth base in human DNA, 5-methylcytosine, is inherently mutagenic. This has led to marked changes in the distribution of the CpG methyl acceptor site and an 80% depletion in its frequency of occurrence in vertebrate DNA. The coding regions of many genes contain CpGs which are methylated in sperm and serve as hot spots for mutation in human genetic diseases. Fully 30-40% of all human germline point mutations are thought to be methylation induced even though the CpG dinucleotide is under-represented and efficient cellular repair systems exist. Importantly, tumor suppressor genes such as p53 also contain methylated CpGs and these serve as hot spots for mutations in some, but not all, human cancers. Comparison of the spectrum of mutations present in this gene in different human cancers allows for predictions to be made on the molecular mechanisms of tumorigenesis.