Differential DNA Methylation Analysis without a Reference Genome
Ontology highlight
ABSTRACT: Genome-wide DNA methylation mapping uncovers epigenetic changes associated with animal development, environmental adaptation, and species evolution. To address the lack of high-throughput methods for studying DNA methylation in non-model organisms, we developed an integrated approach for studying DNA methylation differences without a reference genome. Experimentally, our method relies on an optimized 96-well protocol for reduced representation bisulfite sequencing (RRBS), which we have validated in nine species (human, mouse, rat, cow, dog, chicken, zebrafish, carp, and sea bass). Bioinformatically, we developed the RefFreeDMA software (http://RefFreeDMA.computational-epigenetics.org) to deduce ad hoc genomes directly from RRBS reads and to pinpoint differentially methylated regions. These regions are interpreted using motif enrichment analysis and/or cross-mapping to annotated genomes. We validated our method by reference-free analysis of cell type-specific DNA methylation in the blood of human, cow, and carp. In summary, we present a cost-effective method for epigenome analysis in ecology and evolution, which enables epigenome-wide association studies in natural populations and species without a reference genome.
ORGANISM(S): Bos taurus Homo sapiens Cyprinus carpio
Project description:Sperm cells are characterized by a unique epigenome, and an incorrect establishment of DNA methylation patterns during the differentiation of male germ cells into spermatozoa has been associated with male infertility in several species. While bull semen is widely used in artificial insemination, the literature describing DNA methylation in bovine sperm is still scarce. The purpose of this study is to compare the methylomes of sperm and somatic cell types in cattle using the RRBS technology.
Project description:Reduced representation bisulfite sequencing (RRBS) has been proven a powerful method in DNA methylome profiling. Since the initial development of this method, the RRBS protocol has been modified in order to optimize it for genomic coverage, starting material, and library-construction throughput, which has resulted in new methods such as enhanced RRBS (ERRBS), double-enzyme RRBS (dRRBS), gel-free and multiplexed RRBS (mRRBS), and single-cell RRBS (scRRBS). However, each of these methods has failed to address PCR-derived duplication artifacts, which can bias the results of DNA methylation analyses. To overcome the aforementioned complication, we developed quantitative RRBS (Q-RRBS), a method in which unique molecular identifiers (UMIs) are used to eliminate PCR-induced duplication. By performing Q-RRBS on varying amounts of starting material, we determined that duplication-induced artifacts were more severe when small quantities of the starting material were used. However, through using the UMIs, we successfully eliminated these artifacts. Our results demonstrate that Q-RRBS is an optimal strategy for DNA methylation profiling of single cells or samples containing ultra-trace amounts of cells.
Project description:We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale, base-resolution DNA methylation profiles of primary tissue samples from various organs. Reference-genome independent analysis of this comprehensive dataset defined a “genomic code” of DNA methylation, which allowed us to predict global and locus-specific DNA methylation from the DNA sequence within and across species. This code appears broadly conserved throughout vertebrate evolution, with two major transitions – once in the first vertebrates and again with the emergence of reptiles. Beyond the central role of species-specific DNA sequence composition, our dataset identified the tissue type and the individual as two main sources of DNA methylation variability within species. Tissue type was the dominant factor in fish, birds, and mammals, while in invertebrates, reptiles, and amphibians both factors were similarly strong. Cross-species comparisons focusing on heart and liver tissues supported a highly conserved role of DNA methylation for tissue type and identity and cross-mapping based promoter methylation analysis revealed divergence at specific genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.
Project description:Background: Researching the murine epigenome in disease models has been hampered by the lack of an appropriate and cost-effective DNA methylation array. Until recently, investigators have been limited to the relatively expensive and analysis intensive bisulphite sequencing methods. Here, we performed a comprehensive, comparative analysis between the new Mouse Methylation BeadChip (MMB) and reduced representation bisulphite sequencing (RRBS) in two murine models of colorectal carcinogenesis, providing insight into the utility to each platforms in a real world environment. Results: We captured 1.47x106 CpGs by RRBS and 2.64x105 CpGs by MMB, mapping to 13,778 and 13,365 CpG islands, respectively. RRBS captured significantly more CpGs per island (median 41 for RRBS versus 2 for MMB). We found that 64.4% of intra-island CpG methylation variability can be captured by measuring approximately one quarter of CpG island (CGI) CpGs. MMB was more precise in measuring DNA methylation, especially at sites that had low RRBS coverage. This impacted differential methylation analysis, with more statistically significantly differentially methylated CpG sites identified by MMB in all experimental conditions, however the difference was minute when appropriate thresholding for the magnitude of methylation change (0.2 beta value difference) was applied, providing confidence that both techniques can identify similar differential DNA methylation. Gene ontology enrichment analysis of differentially hypermethylated gene promoters identified similar biological processes and pathways by both RRBS and MMB across two murine model systems. Conclusion: MMB is an effective tool for profiling the murine methylome that performs comparably to RRBS, identifying similar differentially methylated pathways. Although MMB captures a similar proportion of CpG islands, it does so with fewer CpGs per island. We show that subsampling informative CpGs from CpG islands is an appropriate strategy to capture whole island variation. Choice of technology is experiment dependent and will be predicated on the underlying biology being probed.
Project description:Background: Researching the murine epigenome in disease models has been hampered by the lack of an appropriate and cost-effective DNA methylation array. Until recently, investigators have been limited to the relatively expensive and analysis intensive bisulphite sequencing methods. Here, we performed a comprehensive, comparative analysis between the new Mouse Methylation BeadChip (MMB) and reduced representation bisulphite sequencing (RRBS) in two murine models of colorectal carcinogenesis, providing insight into the utility to each platforms in a real world environment. Results: We captured 1.47x106 CpGs by RRBS and 2.64x105 CpGs by MMB, mapping to 13,778 and 13,365 CpG islands, respectively. RRBS captured significantly more CpGs per island (median 41 for RRBS versus 2 for MMB). We found that 64.4% of intra-island CpG methylation variability can be captured by measuring approximately one quarter of CpG island (CGI) CpGs. MMB was more precise in measuring DNA methylation, especially at sites that had low RRBS coverage. This impacted differential methylation analysis, with more statistically significantly differentially methylated CpG sites identified by MMB in all experimental conditions, however the difference was minute when appropriate thresholding for the magnitude of methylation change (0.2 beta value difference) was applied, providing confidence that both techniques can identify similar differential DNA methylation. Gene ontology enrichment analysis of differentially hypermethylated gene promoters identified similar biological processes and pathways by both RRBS and MMB across two murine model systems. Conclusion: MMB is an effective tool for profiling the murine methylome that performs comparably to RRBS, identifying similar differentially methylated pathways. Although MMB captures a similar proportion of CpG islands, it does so with fewer CpGs per island. We show that subsampling informative CpGs from CpG islands is an appropriate strategy to capture whole island variation. Choice of technology is experiment dependent and will be predicated on the underlying biology being probed.
Project description:We describe XmaI-RRBS method for rapid and affordable genome-wide DNA methylation analysis, with library preparation taking only four days and sequencing possible within four hours. Small sizes of the XmaI-RRBS libraries allow their multiplexing and sequencing on the benchtop high-throughput machines. Described here is the first RRBS protocol validated for the Ion Torrent Personal Genome Machine.
Project description:Reduced representation bisulfite sequencing (RRBS) provides an efficient method for measuring DNA methylation at single base resolution in regions of high CpG density. This technique has been extensively tested on the HiSeq2500, which uses a 4-colour detection method, however it is unclear if the method will also work on the NextSeq500 platform, which employs a 2-colour detection system. We created an RRBS library and sequenced it on both the HiSeq2500 and NextSeq500, and found no significant difference in the base composition of reads derived from either machine. Moreover, the methylation calls made from the data of each instrument were highly concordant, with methylation patterns across the genome appearing as expected. Therefore, RRBS can be sequenced on the Nextseq500 with comparable quality to that of the HiSeq2500.
Project description:We describe XmaI-RRBS method for rapid and affordable genome-wide DNA methylation analysis, with library preparation taking only four days and sequencing possible within four hours. Small sizes of the XmaI-RRBS libraries allow their multiplexing and sequencing on the benchtop high-throughput machines. Described here is the first RRBS protocol validated for the Ion Torrent Personal Genome Machine. DNA from MCF7 cell line and 6 normal breast samples (total 7 samples) were subjected to reduced representation bisulfite sequencing analysis (XmaI-RRBS) by using Ion Torrent platform.
Project description:Dietary folate is a major source of methyl groups required for DNA methylation, an epigenetic modification that is actively maintained and remodelled during spermatogenesis. While high dose folic acid supplementation (up to ten times the daily recommended dose) has been shown to improve sperm parameters in infertile men, the effects of supplementation on the sperm epigenome are unknown. To assess the impact of six months of high dose folic acid supplementation on the sperm epigenome, we studied 30 men with idiopathic infertility. Blood folate concentrations increased significantly after supplementation with no significant improvements in sperm parameters. Methylation levels of the differentially methylated regions of several imprinted loci (H19, DLK1/GTL2, MEST, SNRPN, PLAGL1, KCNQ1OT1) were normal both before and after supplementation. Reduced representation bisulfite sequencing (RRBS) revealed a significant global loss of methylation across different regions of the sperm genome. The most marked loss of DNA methylation was found in sperm from patients homozygous for the methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism, a common polymorphism in a key enzyme required for folate metabolism. RRBS analysis also showed that most of the differentially methylated tiles were located in DNA repeats, low CpG density and intergenic regions. Ingenuity Pathway Analysis revealed that methylation of promoter regions was altered in several genes involved in cancer and neurobehavioral disorders including CBFA2T3, PTPN6, COL18A1, ALDH2, UBE4B, ERBB2, GABRB3, CNTNAP4 and NIPA1. Our data reveal alterations of the human sperm epigenome associated with high dose folic acid supplementation, effects that were exacerbated by a common polymorphism in MTHFR. Reduced representation bisulfite sequencing of 28 human sperm samples before and after 6 month of high dose folic acid supplementation.
Project description:The placental epigenome plays a vital role in regulating mammalian growth and development. Aberrations in placental DNA methylation are linked to several disease states, including intrauterine growth restriction and preeclampsia. Studying the evolution and development of the placental epigenome is critical to understanding the origin and progression of such diseases. Although high resolution studies have found substantial variation between placental methylomes of different species, the nature of methylome variation has yet to be characterized within any individual species. We conducted a study of placental DNA methylation at high resolution in multiple strains and closely related species of house mice (Mus musculus musculus, Mus m. domesticus, and M. spretus), across developmental timepoints (embryonic days 15 to 18), and between the labyrinthine transport and junctional endocrine layers. We observed substantial genome-wide methylation heterogeneity in mouse placenta compared to other differentiated tissues. Species-specific methylation profiles were concentrated in retrotransposon subfamilies, specifically RLTR10 and RLTR20 subfamilies. Regulatory regions such as gene promoters and CpG islands displayed cross-species conservation, but showed strong differences between layers and developmental timepoints. Partially methylated domains exist in the mouse placenta and widen during development. Taken together, our results characterize the mouse placental methylome as a highly heterogeneous and deregulated landscape globally, intermixed with actively regulated promoter and retrotransposon sequences.