Whole-genome Shotgun Bisulfite Sequencing Research Articles

Repetitive DNA methylome analysis by small-scale and single-cell shotgun bisulfite sequencing.

Whole-genome shotgun bisulfite sequencing (WG-SBS) is currently the most powerful tool available for understanding genomewide cytosine methylation with single-base resolution; however, the high sequencing cost limits its widespread application, particularly for mammalian genomes. We mapped high- to low-coverage SBS short reads of mouse and human female developing germ cells to consensus sequences of repetitive elements that were multiplied in the respective host genome. This mapping strategy effectively identified active and evolutionarily young retrotransposon subfamilies and centromeric satellite repeats that were resistant to DNA demethylation during the investigated progressive stages of germ cell development. Notably, quantities of only tens of thousands of uniquely mapped reads provided sufficient sensitivity to allow for methylation analyses of multiple retrotransposons and satellite repeats in mice. Furthermore, we produced SBS results from single female murine germ cells by an improved multiplexing and amplification-free SBS method (scPBAT). The scPBAT results quantitatively provided ≥5× sequencing coverage for at least 30 repeats, and the individual methylation patterns detected were similar to the bulk cell-based results. Our single-cell methylome sequencing technique will allow researchers to investigate intergenic methylation characteristics from limited amounts of mammalian cells as well as cells from other organisms with genomic annotations.

DNA methylation map of human atherosclerosis.

Epigenetic alterations may contribute to the development of atherosclerosis. In particular, DNA methylation, a reversible and highly regulated DNA modification, could influence disease onset and progression because it functions as an effector for environmental influences, including diet and lifestyle, both of which are risk factors for cardiovascular diseases. To address the role of DNA methylation changes in atherosclerosis, we compared a donor-matched healthy and atherosclerotic human aorta sample using whole-genome shotgun bisulfite sequencing. We observed that the atherosclerotic portion of the aorta was hypermethylated across many genomic loci in comparison with the matched healthy counterpart. Furthermore, we defined specific loci of differential DNA methylation using a set of donor-matched aortic samples and a high-density (>450 000 CpG sites) DNA methylation microarray. The functional importance in the disease was corroborated by crossing the DNA methylation signature with the corresponding expression data of the same samples. Among the differentially methylated CpGs associated with atherosclerosis onset, we identified genes participating in endothelial and smooth muscle functions. These findings provide new clues toward a better understanding of the molecular mechanisms of atherosclerosis. Our data identify an atherosclerosis-specific DNA methylation profile that highlights the contribution of different genes and pathways to the disorder. Interestingly, the observed gain of DNA methylation in the atherosclerotic lesions justifies efforts to develop DNA demethylating agents for therapeutic benefit.

Intronic Non-CG DNA hydroxymethylation and alternative mRNA splicing in honey bees

BackgroundPrevious whole-genome shotgun bisulfite sequencing experiments showed that DNA cytosine methylation in the honey bee (Apis mellifera) is almost exclusively at CG dinucleotides in exons. However, the most commonly used method, bisulfite sequencing, cannot distinguish 5-methylcytosine from 5-hydroxymethylcytosine, an oxidized form of 5-methylcytosine that is catalyzed by the TET family of dioxygenases. Furthermore, some analysis software programs under-represent non-CG DNA methylation and hydryoxymethylation for a variety of reasons. Therefore, we used an unbiased analysis of bisulfite sequencing data combined with molecular and bioinformatics approaches to distinguish 5-methylcytosine from 5-hydroxymethylcytosine. By doing this, we have performed the first whole genome analyses of DNA modifications at non-CG sites in honey bees and correlated the effects of these DNA modifications on gene expression and alternative mRNA splicing.ResultsWe confirmed, using unbiased analyses of whole-genome shotgun bisulfite sequencing (BS-seq) data, with both new data and published data, the previous finding that CG DNA methylation is enriched in exons in honey bees. However, we also found evidence that cytosine methylation and hydroxymethylation at non-CG sites is enriched in introns. Using antibodies against 5-hydroxmethylcytosine, we confirmed that DNA hydroxymethylation at non-CG sites is enriched in introns. Additionally, using a new technique, Pvu-seq (which employs the enzyme PvuRts1l to digest DNA at 5-hydroxymethylcytosine sites followed by next-generation DNA sequencing), we further confirmed that hydroxymethylation is enriched in introns at non-CG sites.ConclusionsCytosine hydroxymethylation at non-CG sites might have more functional significance than previously appreciated, and in honey bees these modifications might be related to the regulation of alternative mRNA splicing by defining the locations of the introns.

Amplification-free whole-genome shotgun bisulfite sequencing of mouse primordial germ cells

Background Dynamic epigenetic reprogramming occurs during mammalian germ cell development, whereas the targets of this process including DNA demethylation and de novo methylation remain poorly understood. Here, we examined genome-wide methylation profiles in developing primordial germ cells (PGCs) of mice using high-throughput shotgun sequencing of bisulfite-treated DNA (whole-genome shotgun bisulfite sequencing; WGSBS), which accurately quantifies whole-genome methylation levels at single-base resolution.

High-resolution DNA methylome analysis of primordial germ cells identifies gender-specific reprogramming in mice

Dynamic epigenetic reprogramming occurs during mammalian germ cell development, although the targets of this process, including DNA demethylation and de novo methylation, remain poorly understood. We performed genome-wide DNA methylation analysis in male and female mouse primordial germ cells at embryonic days 10.5, 13.5, and 16.5 by whole-genome shotgun bisulfite sequencing. Our high-resolution DNA methylome maps demonstrated gender-specific differences in CpG methylation at genome-wide and gene-specific levels during fetal germline progression. There was extensive intra- and intergenic hypomethylation with erasure of methylation marks at imprinted, X-linked, or germline-specific genes during gonadal sex determination and partial methylation at particular retrotransposons. Following global demethylation and sex determination, CpG sites switched to de novo methylation in males, but the X-linked genes appeared resistant to the wave of de novo methylation. Significant differential methylation at a subset of imprinted loci was identified in both genders, and non-CpG methylation occurred only in male gonocytes. Our data establish the basis for future studies on the role of epigenetic modifications in germline development and other biological processes.

Single Nucleotide Analysis of Cytosine Methylation by Whole‐Genome Shotgun Bisulfite Sequencing

Unlike other methods to assess methylation, whole-genome shotgun bisulfite sequencing is used to generate quantitative genome-wide methylation profiles at single-nucleotide resolution. As described in this unit, this method allows for the quantitative measurement of methylation at each cytosine base by treatment of genomic DNA with sodium bisulfite followed by sequencing and alignment of the reads to a reference genome.

Contribution of Intragenic DNA Methylation in Mouse Gametic DNA Methylomes to Establish Oocyte-Specific Heritable Marks

Genome-wide dynamic changes in DNA methylation are indispensable for germline development and genomic imprinting in mammals. Here, we report single-base resolution DNA methylome and transcriptome maps of mouse germ cells, generated using whole-genome shotgun bisulfite sequencing and cDNA sequencing (mRNA-seq). Oocyte genomes showed a significant positive correlation between mRNA transcript levels and methylation of the transcribed region. Sperm genomes had nearly complete coverage of methylation, except in the CpG-rich regions, and showed a significant negative correlation between gene expression and promoter methylation. Thus, these methylome maps revealed that oocytes and sperms are widely different in the extent and distribution of DNA methylation. Furthermore, a comparison of oocyte and sperm methylomes identified more than 1,600 CpG islands differentially methylated in oocytes and sperm (germline differentially methylated regions, gDMRs), in addition to the known imprinting control regions (ICRs). About half of these differentially methylated DNA sequences appear to be at least partially resistant to the global DNA demethylation that occurs during preimplantation development. In the absence of Dnmt3L, neither methylation of most oocyte-methylated gDMRs nor intragenic methylation was observed. There was also genome-wide hypomethylation, and partial methylation at particular retrotransposons, while maintaining global gene expression, in oocytes. Along with the identification of the many Dnmt3L-dependent gDMRs at intragenic regions, the present results suggest that oocyte methylation can be divided into 2 types: Dnmt3L-dependent methylation, which is required for maternal methylation imprinting, and Dnmt3L-independent methylation, which might be essential for endogenous retroviral DNA silencing. The present data provide entirely new perspectives on the evaluation of epigenetic markers in germline cells.

Pash 3.0: A versatile software package for read mapping and integrative analysis of genomic and epigenomic variation using massively parallel DNA sequencing

BackgroundMassively parallel sequencing readouts of epigenomic assays are enabling integrative genome-wide analyses of genomic and epigenomic variation. Pash 3.0 performs sequence comparison and read mapping and can be employed as a module within diverse configurable analysis pipelines, including ChIP-Seq and methylome mapping by whole-genome bisulfite sequencing.ResultsPash 3.0 generally matches the accuracy and speed of niche programs for fast mapping of short reads, and exceeds their performance on longer reads generated by a new generation of massively parallel sequencing technologies. By exploiting longer read lengths, Pash 3.0 maps reads onto the large fraction of genomic DNA that contains repetitive elements and polymorphic sites, including indel polymorphisms.ConclusionsWe demonstrate the versatility of Pash 3.0 by analyzing the interaction between CpG methylation, CpG SNPs, and imprinting based on publicly available whole-genome shotgun bisulfite sequencing data. Pash 3.0 makes use of gapped k-mer alignment, a non-seed based comparison method, which is implemented using multi-positional hash tables. This allows Pash 3.0 to run on diverse hardware platforms, including individual computers with standard RAM capacity, multi-core hardware architectures and large clusters.