Differentially Methylated Region (DMR) detection identifies contiguous genomic regions with methylation differences between conditions, using methods like tile-based (methylKit), smoothing-based (bsseq BSmooth), or kernel-based (DMRcate) approaches.
if (!require('BiocManager', quietly = TRUE))
install.packages('BiocManager')
BiocManager::install(c('methylKit', 'bsseq', 'DMRcate', 'DSS'))Tell your AI agent what you want to do:
- "Find differentially methylated regions between my treatment groups"
- "Detect DMRs using BSmooth smoothing for my WGBS data"
- "Identify promoter regions with altered methylation"
"Which DMR method should I use for my WGBS data with low coverage?"
"Compare tile-based vs smoothing-based DMR detection for my samples"
"Find DMRs using 500bp tiles with methylKit"
"Run tile-based DMR analysis with stringent cutoffs"
"Detect DMRs using BSmooth smoothing for my WGBS experiment"
"Run bsseq DMR analysis with appropriate smoothing parameters"
"Annotate my DMRs with gene names and genomic features"
"Find DMRs overlapping promoters and CpG islands"
"Compare DMRs with differentially expressed genes"
"Identify DMRs in enhancer regions from my ChIP-seq data"
- Help select appropriate DMR method based on your data type
- Set tile size or smoothing parameters for your analysis
- Run DMR detection with chosen method
- Filter DMRs by q-value, methylation difference, and number of CpGs
- Annotate DMRs with genomic features
- Export results for visualization or integration
| Method | Pros | Cons | Best For |
|---|---|---|---|
| methylKit tiles | Simple, fast | Fixed windows | Quick exploration |
| BSmooth | Handles low coverage | Computationally intensive | WGBS |
| DMRcate | Array-optimized | Less flexible | 450K/EPIC arrays |
| DSS | Statistical rigor | Complex | Multi-factor designs |
| Parameter | Lenient | Moderate | Stringent |
|---|---|---|---|
| qvalue | < 0.1 | < 0.05 | < 0.01 |
| meth.diff | > 10% | > 25% | > 40% |
| CpGs | >= 3 | >= 5 | >= 10 |
- Run single-CpG analysis first to understand data quality before DMR detection
- Tile size recommendations: WGBS 500-1000bp, RRBS 100-500bp
- BSmooth works better with low coverage data due to smoothing
- For publication, use BSmooth or DSS rather than simple tiles
- If few DMRs found: relax thresholds, use smaller tiles, check sample quality
- If too many DMRs: increase meth.diff threshold, lower q-value cutoff
- Always annotate DMRs to promoters, CpG islands, and gene bodies for interpretation
- Consider integrating with expression data to prioritize functionally relevant DMRs