Abstract:5-methylcytosine (5mC) is the most common chemical modification occurring on the CpG sites across the human genome. Bisulfite conversion combined with short-read whole genome sequencing can capture and quantify the modification at single nucleotide resolution. However, the PCR amplification process could lead to duplicative methylation patterns and introduce 5mC detection bias. Additionally, the limited read length also restricts co-methylation analysis between distant CpG sites. The bisulfite conversion process presents a significant challenge for detecting variant-specific methylation due to the destruction of allele information in the sequencing reads. To address these issues, we sought to characterize the human methylation profiling with the nanopore long-read sequencing, aiming to demonstrate its potential for long-range co-methylation analysis with native modification call and intact allele information retained. In this regard, we first analyzed the nanopore demo data in the adaptive sampling sequencing run targeting all human CpG islands. We applied the linkage disequilibrium (LD) R2 to calculate the co-methylation in nanopore data, and further identified 27,875, 50,481, 26,542 and 51,189 methylation haplotype blocks (MHB) in COLO829, COLO829BL, HCC1395 and HCC1395BL cell lines, respectively. Interestingly, while we found that majority of the co-methylation were in a short range (≤200bp), a small portion (1~3%) showed long distance (≥1,000bp), suggesting potential remote regulatory mechanisms across the genome. To further characterize the epigenetic changes related to transcription factor binding, we profiled the 5mC percentage changes surrounding various motif sites in JASPAR collection and found that CTCF and KLF5 binding sites showed reduced methylation, while FOXE1 and ZNF354A sites showed increased methylation. To further investigate the allele-specific 5mCG in the prostate genome, we designed a target region covering methylation quantitative trait loci (mQTL) and genome-wide association study (GWAS) risk germline variants and generated long reads with adaptive sampling run in the 22Rv1 cell line. To identify the allele-specific methylation in the 22Rv1 cell line, we performed long-read based phasing and compared the 5mCG signals between the two haplotypes. As a result, we identified 6,390 haplotype-specific methylated regions in the 22Rv1 cell line (p-MWU ≤ 1e-5 and delta ≥ 50%). By examining haplotype-specific methylated regions near the phasing variants, we identified examples of allele-specific methylated regions that showed allele-specific accessibility in the ATAC-seq data. By further integrating the ATAC-seq data of 22Rv1, we found that methylation levels were negatively correlated with chromatin accessibility at the genome-wide scale. Our study has revealed native methylome profiling while preserving haplotype information, offering a novel approach to uncovering the regulatory mechanisms of the human prostate genome.

MethHaplo: combining allele-specific DNA methylation and SNPs for haplotype region identification

Simultaneous profiling of DNA methylation and chromatin architecture in mixed populations and in single cells

CpG_MPs: Identification of CpG Methylation Patterns of Genomic Regions from High-Throughput Bisulfite Sequencing Data

MethylGenotyper: Accurate Estimation of SNP Genotypes and Genetic Relatedness from DNA Methylation Data

Identification of Methylation States of DNA Regions for Illumina Methylation BeadChip

Accurate estimation of SNP genotypes and genetic relatedness from DNA methylation data

Long-range Polony Haplotyping of Individual Human Chromosome Molecules

High accuracy methylation identification tools on single molecular level for PacBio HiFi data

Fine Mapping Regulatory Variants by Characterizing Native CpG Methylation with Nanopore Long Read Sequencing

Functional DNA methylation differences between tissues, cell types, and across individuals discovered using the M&M algorithm

Identification of methylation and hydroxymethylation haplotype blocks as distal regulatory elements aids in deconvolution of heterogeneous brain tissues

Identification of Methylated Regions with Peak Search Based on Poisson Model from Massively Parallel Methylated Dna Immunoprecipitation-Sequencing Data

Identification of methylation haplotype blocks aids in deconvolution of heterogeneous tissue samples and tumor tissue-of-origin mapping from plasma DNA

DNA Methylome Profiling at Single-Base Resolution Through Bisulfite Sequencing of 5Mc-Immunoprecipitated DNA.

Uncovering the roles of DNA hemi-methylation in transcriptional regulation using MspJI-assisted hemi-methylation sequencing

A New Approach to Decode DNA Methylome and Genomic Variants Simultaneously from Double Strand Bisulfite Sequencing.

High resolution profiling of human exon methylation by liquid hybridization capture-based bisulfite sequencing

A DNA Methylation Haplotype Block Landscape in Human Tissues and Preimplantation Embryos Reveals Regulatory Elements Defined by Comethylation Patterns.

Development of a Method for Identifying and Functionally Analyzing Allele-Specific DNA Methylation Based on BS-seq Data

Methods for discovering genomic loci exhibiting complex patterns of differential methylation

CGmapTools Improves the Precision of Heterozygous SNV Calls and Supports Allele-Specific Methylation Detection and Visualization in Bisulfite-Sequencing Data.