This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Unable to use rtracklayer import of bigwig files in older chipseq pipeline - work around/alternatives?

Hello and thank you for anyone willing to spare a moment to help.

I am really struggling trying to execute a pipeline that worked two months ago but has since stopped.

https://github.com/BoevaLab/LILY

It fails when trying to import a BigWigFile using rtracklayer. Trying to troubleshoot with the good'ol chatbot, it says that a bioconductor update has broken an internal methods call.

Below is the function call at which it is failing and the main function itself...

It's the import.bw that is failing. I can load these bigwigs into IGV and visualize them, so I don't think it's the file itself...

> total_densities = apply(files_info,1,get_peaks_density_commonPeaks)
Error: unable to find an inherited method for function ‘seqinfo’ for signature ‘x = "BigWigFile"’

get_peaks_density_commonPeaks <- function(info){
  peaks = import.bed3(info["peak"])
  rangeToConsider=c(100:5000)
  peaks = peaks[rangeToConsider]
  peaks= peaks[which(countOverlaps(peaks,commonPeaks)>0)]
  densities = import.bw(info["density"])
  peak_density = subsetByOverlaps(densities,peaks)  
  total_density = median(peak_density$score)
  gc()
  return(total_density)
}





sessionInfo( )

R version 4.5.1 (2025-06-13)
Platform: x86_64-pc-linux-gnu
Running under: Ubuntu 24.04.3 LTS

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.12.0 
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.12.0  LAPACK version 3.12.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8     LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                  LC_ADDRESS=C               LC_TELEPHONE=C             LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

time zone: America/Toronto
tzcode source: system (glibc)

attached base packages:
[1] stats4    stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] rtracklayer_1.68.0   GenomicRanges_1.62.0 Seqinfo_1.0.0        IRanges_2.44.0       S4Vectors_0.48.0     BiocGenerics_0.56.0  generics_0.1.4      

loaded via a namespace (and not attached):
 [1] Matrix_1.7-3                jsonlite_2.0.0              BiocManager_1.30.27         compiler_4.5.1              rjson_0.2.23               
 [6] crayon_1.5.3                SummarizedExperiment_1.40.0 Biobase_2.70.0              Rsamtools_2.24.1            bitops_1.0-9               
[11] Biostrings_2.78.0           GenomicAlignments_1.44.0    parallel_4.5.1              BiocParallel_1.42.2         yaml_2.3.12                
[16] lattice_0.22-5              R6_2.6.1                    XVector_0.50.0              S4Arrays_1.10.0             curl_7.0.0                 
[21] GenomeInfoDb_1.44.3         XML_3.99-0.20               DelayedArray_0.36.0         MatrixGenerics_1.22.0       GenomeInfoDbData_1.2.14    
[26] rlang_1.1.6                 pkgload_1.4.1               cli_3.6.5                   SparseArray_1.10.3          grid_4.5.1                 
[31] rstudioapi_0.17.1           codetools_0.2-20            abind_1.4-8                 RCurl_1.98-1.17             restfulr_0.0.16            
[36] httr_1.4.7                  matrixStats_1.5.0           tools_4.5.1                 BiocIO_1.18.0               UCSC.utils_1.4.0
bioconductor rtracklayer chipseq lily

1 answer

I encountered a similar error and fixed it by updating some of the problematic packages - This was done via BiocManager::valid() and following the instructions from there. I hope this helps.

Log in to answer this question.