Abstract

DNA methylation in gene or promoter or gene body could restrict/promote the gene transcription. Moreover, methylation in the gene regions along with CpG island regions could modulate the transcription to undetectable gene expression levels. Therefore, it is necessary to investigate the methylation levels within the gene, gene body, CpG island regions and their overlapped regions and then identify the gene-based differentially methylated regions (GeneDMRs). Here, R package GeneDMRs aims to facilitate computing gene based methylation rate using next generation sequencing (NGS)-based methylome data. The user-friendly R package GeneDMRs is presented to analyze the methylation levels in each gene/promoter/exon/intron/CpG island/CpG island shore or each overlapped region (e.g., gene-CpG island/promoter-CpG island/exon-CpG island/intron-CpG island/gene-CpG island shore/promoter-CpG island shore/exon-CpG island shore/intron-CpG island shore). Here, we used the public reduced representation bisulfite sequencing (RRBS) data of mouse (GSE62392) for evaluating software and found novel biologically significant results to supplement the previous research. The R package GeneDMRs can facilitate computing gene based methylation rate to interpret complex interplay between methylation levels and gene expression differences or similarities across physiological conditions or disease states.