PoreMeth2: decoding the evolution of methylome alterations with Nanopore sequencing

In epigenetic analysis, identifying differentially methylated regions (DMRs) typically involves detecting groups of consecutive CpGs that show significant changes in their average methylation levels. However, the methylation state of a genomic region can also be characterized by a mixture of patterns (epialleles) with variable frequencies, and the relative proportions of such patterns can provide insights into its mechanisms of formation. Traditional methods based on bisulfite conversion and NGS, due to the read size (150 bp), allow epiallele frequency analysis only in high-CpG-density regions, limiting differential methylation studies to just 50% of the human methylome. Nanopore sequencing, with its long reads, enables the analysis of epiallele frequency across both high- and low-CpG-density regions. We introduce a novel computational approach, PoreMeth2, an R library that integrates epiallelic diversity and methylation frequency changes from Nanopore data to identify DMRs, assess their formation mechanisms, and annotate them to genic and regulatory elements. We applied PoreMeth2 to cancer and glial cell datasets, demonstrating its ability to distinguish epigenomic changes with a strong effect on gene expression from those with a weaker impact on transcriptional activity. PoreMeth2 is publicly available at https://github.com/Lab-CoMBINE/PoreMeth2.