Isolation and enrichment of human genomic CpG islands by methylation-sensitive mirror orientation selection

CpG islands (CGIs) in human genomic DNA are GC-rich fragments whose aberrant methylation is associated with human disease development. In the current study, methylation-sensitive mirror orientation selection (MS-MOS) was developed to efficiently isolate and enrich unmethylated CGIs from human genomic DNA. The unmethylated CGIs prepared by the MS-MOS procedure subsequently were used to construct a CGI library. Then the sequence characteristics of cloned inserts of the library were analyzed by bioinformatics tools, and the methylation status of CGI clones was analyzed by HpaII PCR. The results showed that the MS-MOS method could be used to isolate up to 0.001% of differentially existed unmethylated DNA fragments in two complex genomic DNA. In the CGI library, 34.1% of clones had insert sequences satisfying the minimal criteria for CGIs. Excluding duplicates, 22.0% of the 80,000 clones were unique CGI clones, representing 60% of all the predicted CGIs (about 29,000) in human genomic DNA, and most or all of the CGI clones were unmethylated in human normal cell DNA based on the HpaII PCR analysis results of randomly selected CGI clones. In conclusion, MS-MOS was an efficient way to isolate and enrich human genomic CGIs. The method has powerful potential application in the comprehensive identification of aberrantly methylated CGIs associated with human tumorigenesis to improve understanding of the epigenetic mechanisms involved.     

Identification of differential genomic DNA Methylation in the hypothalamus of pubertal rat using reduced representation Bisulfite sequencing

Background

There are many variables affecting the onset of puberty in animals, including genetic, nutritional, and environmental factors. Recent studies suggest that epigenetic regulation, especially DNA methylation, plays a majorrole in the regulation of puberty. However, there have been no reports on DNA methylation of the pubertal genome.

Methods

We investigated DNA methylation in the female rat hypothalamus at prepuberty and puberty using reduced representation bisulfite sequencing technology. The identified genes and signaling pathways exhibiting changes to DNA methylation in pubertal rats were determined by Gene Ontogeny and Kyoto Encyclopedia of Genes and Genomes analysis.

Results

The distribution of the three types of methylated C bases in promoter and CpG island (CGI) regions in the hypothalamus was as follows: 87.79% CG, 3.05% CHG, 9.16% CHH for promoters, and 88.35% CG, 3.21% CHG, 88.35% CHH for CGI in prepubertal rats; and 90.78% CG, 2.13% CHG, 7.09% CHH for promoters, and 88.59% CG, 88.59% CHG, 8.35% CHH for CGI in pubertal animals. CG showed the highest percentage of methylation, and was the highest methylation state in CGI. Compared to prepubertal hyoyhalamus samples, we identified ten genes with altered methylation in promoter regions in the pubertal hypothalamus samples, and 43 genes with altered methylation in the CGI. Changes in DNA methylation were found in gonadotropin-releasing hormone signaling pathways, and the oocyte meiosis pathway.

Conclusion

Our results demonstrate changes in DNA methylation occur in female rats from prepuberty to puberty suggestng DNA methylation may play a crucial role in the regulation of puberty onset. This study provides essential information for future studies on the role of epigenetics in the regulation of puberty.

     

Identifying differentially methylated genes using mixed effect and generalized least square models

Background  

DNA methylation plays an important role in the process of tumorigenesis. Identifying differentially methylated genes or CpG islands (CGIs) associated with genes between two tumor subtypes is thus an important biological question. The methylation status of all CGIs in the whole genome can be assayed with differential methylation hybridization (DMH) microarrays. However, patient samples or cell lines are heterogeneous, so their methylation pattern may be very different. In addition, neighboring probes at each CGI are correlated. How these factors affect the analysis of DMH data is unknown.     

CpG methylation recruits sequence specific transcription factors essential for tissue specific gene expression

CG methylation is an epigenetically inherited chemical modification of DNA found in plants and animals. In mammals it is essential for accurate regulation of gene expression and normal development. Mammalian genomes are depleted for the CG dinucleotide, a result of the chemical deamination of methyl-cytosine in CG resulting in TpG. Most CG dinucleotides are methylated, but ~ 15% are unmethylated. Five percent of CGs cluster into ~ 20,000 regions termed CG islands (CGI) which are generally unmethylated. About half of CGIs are associated with housekeeping genes. In contrast, the gene body, repeats and transposable elements in which CGs are generally methylated. Unraveling the epigenetic machinery operating in normal cells is important for understanding the epigenetic aberrations that are involved in human diseases including cancer. With the advent of high-throughput sequencing technologies, it is possible to identify the CG methylation status of all 30 million unique CGs in the human genome, and monitor differences in distinct cell types during differentiation and development. Here we summarize the present understanding of DNA methylation in normal cells and discuss recent observations that CG methylation can have an effect on tissue specific gene expression. We also discuss how aberrant CG methylation can lead to cancer. This article is part of a Special Issue entitled: Chromatin in time and space.     

Development of a Novel Output Value for Quantitative Assessment in Methylated DNA Immunoprecipitation-CpG Island Microarray Analysis

In DNA methylation microarray analysis, quantitative assessment of intermediate methylation levels in samples with various global methylation levels is still difficult. Here, specifically for methylated DNA immunoprecipitation-CpG island (CGI) microarray analysis, we developed a new output value. The signal log ratio reflected the global methylation levels, but had only moderate linear correlation (r = 0.72) with the fraction of DNA molecules immunoprecipitated. By multiplying the signal log ratio using a coefficient obtained from the probability value that took account of signals in neighbouring probes, its linearity was markedly improved (r = 0.94). The new output value, Me value, reflected the global methylation level, had a strong correlation also with the fraction of methylated CpG sites obtained by bisulphite sequencing (r = 0.88), and had an accuracy of 71.8 and 83.8% in detecting completely methylated and unmethylated CGIs. Analysis of gastric cancer cell lines using the Me value showed that methylation of CGIs in promoters and gene bodies was associated with low and high, respectively, gene expression. The degree of demethylation of promoter CGIs after 5-aza-2''-deoxycytidine treatment had no association with that of induction of gene expression. The Me value was considered to be useful for analysis of intermediate methylation levels of CGIs.     

Exome localization of complex disease association signals

Background

Recent progress in high-throughput technologies has greatly contributed to the development of DNA methylation profiling. Although there are several reports that describe methylome detection of whole genome bisulfite sequencing, the high cost and heavy demand on bioinformatics analysis prevents its extensive application. Thus, current strategies for the study of mammalian DNA methylomes is still based primarily on genome-wide methylated DNA enrichment combined with DNA microarray detection or sequencing. Methylated DNA enrichment is a key step in a microarray based genome-wide methylation profiling study, and even for future high-throughput sequencing based methylome analysis.

Results

In order to evaluate the sensitivity and accuracy of methylated DNA enrichment, we investigated and optimized a number of important parameters to improve the performance of several enrichment assays, including differential methylation hybridization (DMH), microarray-based methylation assessment of single samples (MMASS), and methylated DNA immunoprecipitation (MeDIP). With advantages and disadvantages unique to each approach, we found that assays based on methylation-sensitive enzyme digestion and those based on immunoprecipitation detected different methylated DNA fragments, indicating that they are complementary in their relative ability to detect methylation differences.

Conclusions

Our study provides the first comprehensive evaluation for widely used methodologies for methylated DNA enrichment, and could be helpful for developing a cost effective approach for DNA methylation profiling.     

Directional telomeric silencing and lack of canonical B1 elements in two silencer Autonomously Replicating Sequences in S. cerevisiae

Background

Detection of cell-free methylated DNA in plasma is a promising tool for tumour diagnosis and monitoring. Due to the very low amounts of cell-free DNA in plasma, analytical sensitivity is of utmost importance. The vast majority of currently available methods for analysing DNA methylation are based on bisulfite-mediated deamination of cytosine. Cytosine is rapidly converted to uracil during bisulfite treatment, whereas 5-methylcytosine is only slowly converted. Hence, bisulfite treatment converts an epigenetic modification into a difference in sequence, amenable to analysis either by sequencing or PCR based methods. However, the recovery of bisulfite-converted DNA is very poor.

Results

Here we introduce an alternative method for the crucial steps of bisulfite treatment with high recovery. The method is based on an accelerated deamination step and alkaline desulfonation in combination with magnetic silica purification of DNA, allowing preparation of deaminated DNA from patient samples in less than 2 hours.

Conclusions

The method presented here allows low levels of DNA to be easily and reliably analysed, a prerequisite for the clinical usefulness of cell-free methylated DNA detection in plasma.     

Genomic landscape of CpG rich elements in human

Background

The studies on CpG islands (CGI) and Alu elements functions, evolution, and distribution in the genome started since the discovery in nineteen eighties (1981, 1986, correspondingly). Their highly skewed genome wide distribution implies the non-random retrotransposition pattern. Besides CGIs in gene promoters, CGIs clusters were observed in the homeobox gene regions and in the macrosatellites, but the whole picture of their distribution specifics was not grasped. Attempts to identify any causative features upon their (genome wide) distribution, such as the DNA context mediated preferred insertion sites of Alu repeats, have been made to ascribe their clusters location.

Methods

Recent emergence of high resolution 3D map of human genome allowed segregating the genome into the large scale chromatin domains of naturally observable nuclear subcompartments, or Topologically Associated Domains (TADs), designated by spatial chromatin distribution. We utilized the chromatin map to elucidate relations between large scale chromatin state and CpG rich elements landscape.In the course of analysis it was confirmed that genes, Alu and CGI clusters maintain obvious, albeit different in strength, preference for open chromatin. For the first time it was clearly shown that the clusters density of the Alu and CGIs monotonically depend on the chromatin accessibility rate. In particular, the highest density of these elements is found in A1 euchromatin regions characterized by a high density of small length genes replicating in the early S-phase. It implies that these elements mediate (CGIs) or are a side element (Alus) of chromatin accessibility.

Results

We elucidated that both methylated and non-methylated CGIs display the affinity to chromatin accessibility. As a part of comparative genomics section, we elucidated that the dog’s genome non-canonical structure, outstanding in mammals for its high CGIs abundance compared to gene number, is explained by the presence of dense tandem CGI extended hotspots (500 kb on average) in subtelomeric and pericentromeric regions with highly skewed CG content, and not by CGIs global distribution pattern shift.

Conclusions

The study underlines the close association of CG-rich elements distribution with the newly introduced large scale chromatin state map, proposing a refined standpoint on interrelation of aforementioned genome elements and the chromatin state. To our expertise, the TAD-associated partition model employed in the study is likely the most substantial one regarding CpG rich clusters distribution among the whole genome chromatin/isochores maps available.

     

Factors to preserve CpG-rich sequences in methylated CpG islands

Background

Mammalian CpG islands (CGIs) normally escape DNA methylation in all adult tissues and developmental stages. However, in our previous study we unexpectedly identified many methylated CGIs in human peripheral blood leukocytes. Methylated CpG dinucleotides convert to TpG dinucleotides through deaminization of their cytosine bases more frequently than hypomethylated CpG dinucleotides. Therefore, we wondered how methylated CGIs in germline or non-germline cells maintain their CpG-rich sequences. It is known that events such as germline hypomethylation, CpG selection, biased gene conversion (BGC), and frequent CpG fixation can contribute to the maintenance of CpG-rich sequences in methylated CGIs in germline or non-germline cells. However, it has not been investigated which of the processes maintain CpG-rich sequences of methylated CGIs in each genomic position.

Results

In this study, we comprehensively examined the contribution of the processes described above to the maintenance of CpG-rich sequences in methylated CGIs in germline and non-germline cells which were classified by genomic positions. Approximately 60–80% of CGIs with high methylation in H1 cell line (H1-HM) in all the genomic positions showed a low average CpG → TpG/CpA substitution rate. In contrast, fewer than half the numbers of CGIs with H1-HM in all the genomic positions showed a low average CpG → TpG/CpA substitution rate and low levels of methylation in sperm cells (SPM-LM). Furthermore, a small fraction of CGIs with a low average CpG → TpG/CpA substitution rate and high levels of methylation in sperm cells (SPM-HM) showed CpG selection.On the other hand, independent of the positions in genes, most CGIs with SPM-HM showed a slightly higher average TpG/CpA → CpG substitution rate compared with those with SPM-LM.

Conclusions

Relatively high numbers (approximately 60–80%) of CGIs with H1-HM in all the genomic positions preserve their CpG-rich sequences by a low CpG → TpG/CpA substitution rate caused mainly by their SPM-LM, and for those with SPM-HM partly by CpG selection and TpG/CpA → CpG fixation. BGC has little contribution to the maintenance of CpG-rich sequences of CGIs with SPM-HM which were classified by genomic positions.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1286-x) contains supplementary material, which is available to authorized users.