Skip to contents

methfuse hexagon methFuse

methFuse implements FUSE: FUnctional SEgmentation of DNA methylation data through hierarchical clustering.

Features

  • Hierarchical clustering based on methylation count matrices
  • Detection of optimal number of clusters
  • Summarization of segments, including stability flag
  • Per-segment methylation estimates per sample

Installation

From GitHub (Development Version)

Either using remotes: (recommended)

# Install remotes if needed
install.packages("remotes")

# Install fuseR from GitHub
remotes::install_github("holmsusa/methFuse")

or using devtools:

# Install devtools if needed
install.packages("devtools")

# Install fuseR from GitHub
devtools::install_github("holmsusa/methFuse")

System Requirements

  • R version ≥ 4.0
  • C++ toolchain for native code compilation

You may need platform-specific tools:

  • macOS: Xcode Command Line Tools (xcode-select –install)
  • Linux: build-essential
  • Windows: Rtools

Supported input types

fuse.segment() supports the following input formats:

  • Matrix input
    • Unmethylated counts (K0)
    • Methylated counts (K1)
    • Explicit chromosome (chr) and position (pos) vectors
  • BSseq objects (Bioconductor)
    • Counts and genomic coordinates are extracted automatically
    • Requires bsseq package
  • methrix objects (Bioconductor)
    • Supports large datasets via DelayedMatrix
    • Counts and genomic coordinates are extracted automatically
    • Requires methrix and DelayedArray packages

Install needed packages with

BiocManager::install(c("bsseq", "methrix", "DelayedArray"))

Quick Example 

library(fuseR)
set.seed(1234)

# Generate sample data
# Unmethylated counts, T's
K0 <- matrix(
  rep(c(sample(0:20, 200, replace = TRUE), sample(20:40, 200, replace = TRUE)), 2),
  nrow = 100, byrow = TRUE
)
# Methylated counts, C's
K1 <- matrix(
  rep(c(sample(20:40, 200, replace = TRUE), sample(0:20, 200, replace = TRUE)), 2),
  nrow = 100, byrow = TRUE
)

# Perform segmentation
segment_result <- fuse.segment(
  K0, K1, 
  chr = sub("\\..*$", "", rownames(K0)), 
  pos = as.numeric(sub("^.*\\.", "", rownames(K0)))
)

# Access summary and per-segment betas
head(segment_result$summary)
head(segment_result$betas_per_segment)

Check out a full example workflow in the vignette.

License

This package is licensed under the MIT License. See LICENSE for details.

Author

Susanna Holmstrom