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.

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

ASMdb: a comprehensive database for allele-specific DNA methylation in diverse organisms

Genome-scale sequence data processing and epigenetic analysis of DNA methylation]

ASMS: finding allele specific methylation in human genomes without phasing

Large-scale Comparative Epigenomics Reveals Hierarchical Regulation of Non-Cg Methylation in Arabidopsis.

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

Wemics: A Single‐Base Resolution Methylation Quantification Method for Enhanced Prediction of Epigenetic Regulation

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

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

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

The DNA Methylome of Human Peripheral Blood Mononuclear Cells.

Genome-wide profiling of DNA methylation in human cancer cells.

Investigating the potential of single-cell DNA methylation data to detect allele-specific methylation and imprinting

A Bayesian Approach for Analysis of Whole-Genome Bisulfite Sequencing Data Identifies Disease-Associated Changes in DNA Methylation

Base-resolution Analyses of Sequence and Parent-of-origin Dependent DNA Methylation in the Mouse Genome.

Identification of Differentially Methylated Sites with Weak Methylation Effects.

Guide Positioning Sequencing Identifies Aberrant DNA Methylation Patterns That Alter Cell Identity and Tumor-Immune Surveillance Networks

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

Identification of Methylation States of DNA Regions for Illumina Methylation BeadChip

Surface ligation-based resonance light scattering analysis of methylated genomic DNA on a microarray platform.

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