SUMMARY: DNA methylation (DNAm) of neighbouring CpG sites is highly correlated, making DNAm function in terms of blocks. DNAm patterns and functionality are linked to both chromatin structure of DNA and gene regulation. Defining biologically meaningful DNA methylation blocks from whole-genome bisulfite sequencing (WGBS) data remains challenging, as most existing methods rely on fixed genomic windows rather than the observed methylation pattern. We present FUSE, a data-driven segmentation method that captures intrinsic methylation segments directly from WGBS data by jointly analyzing multiple samples. FUSE identifies spatially homogeneous methylation blocks shared across the input cohort while allowing different methylation states across samples. Applied to 61 WGBS samples from the ENCODE database, FUSE identified segments which overlap significantly with promoters, enhancers, and repetitive elements. FUSE was able to recover the true segment breakpoints in synthetic data with high sensitivity under increased levels of noise. As such, FUSE facilitates post hoc methylation analyses by aggregating coherent CpG sites into candidate segments for downstream differential methylation testing or other comparative studies. AVAILABILITY AND IMPLEMENTATION: FUSE is implemented as an R-package methFuse, available at https://github.com/holmsusa/methFuse and https://cran.r-project.org/package=methFuse. A GenomeSpy visualization of the data is available at https://csbi.ltdk.helsinki.fi/p/fuse_encode_gs/.

MOTIVATION: RNA sequencing and other high-throughput technologies are essential in understanding complex diseases, such as cancers, but are susceptible to technical factors manifesting as patterns in the measurements. These batch patterns hinder the discovery of biologically relevant patterns. Unbiased batch effect correction in heterogeneous populations currently requires special experimental designs or phenotypic labels, which are not readily available for patient samples in existing datasets. RESULTS: We present POIBM, an RNA-seq batch correction method, which learns virtual reference samples directly from the data. We use a breast cancer cell line dataset to show that POIBM exceeds or matches the performance of previous methods, while being blind to the phenotypes. Further, we analyze The Cancer Genome Atlas RNA-seq data to show that batch effects plague many cancer types; POIBM effectively discovers the true replicates in stomach adenocarcinoma; and integrating the corrected data in endometrial carcinoma improves cancer subtyping. AVAILABILITY AND IMPLEMENTATION: https://bitbucket.org/anthakki/poibm/ (archived at https://doi.org/10.5281/zenodo.6122436). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.