This vignette exemplifies how to perform unsupervised footprint detection and quantification using FootprintCharter as per Baderna & Barzaghi et al., 2024 and Barzaghi et al., 2024.
FootprintCharter partitions molecules by their methylation patterns without relying on orthogonal genomic annotations such as TF motifs.
Methylation = qs::qread(system.file("extdata", "Methylation_4.qs", package="SingleMoleculeFootprinting"))
RegionOfInterest = GRanges("chr6", IRanges(88106000, 88106500))
TFBSs = qs::qread(system.file("extdata", "TFBSs_1.qs", package="SingleMoleculeFootprinting"))
PlotAvgSMF(MethGR = Methylation[[1]], RegionOfInterest = RegionOfInterest, TFBSs = TFBSs)
## No sorted reads passed...plotting counts from all reads
## R version 4.4.2 (2024-10-31)
## Platform: x86_64-pc-linux-gnu
## Running under: Ubuntu 24.04.1 LTS
##
## Matrix products: default
## BLAS: /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so.3
## LAPACK: /usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.26.so; LAPACK version 3.12.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_US.UTF-8 LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: Etc/UTC
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] BSgenome.Mmusculus.UCSC.mm10_1.4.3 BSgenome_1.75.0
## [3] rtracklayer_1.67.0 BiocIO_1.17.1
## [5] Biostrings_2.75.1 XVector_0.47.0
## [7] GenomicRanges_1.59.1 GenomeInfoDb_1.43.2
## [9] IRanges_2.41.1 S4Vectors_0.45.2
## [11] BiocGenerics_0.53.3 generics_0.1.3
## [13] SingleMoleculeFootprinting_2.1.0
##
## loaded via a namespace (and not attached):
## [1] RColorBrewer_1.1-3 sys_3.4.3
## [3] jsonlite_1.8.9 magrittr_2.0.3
## [5] GenomicFeatures_1.59.1 farver_2.1.2
## [7] rmarkdown_2.29 zlibbioc_1.52.0
## [9] vctrs_0.6.5 memoise_2.0.1
## [11] Rsamtools_2.23.1 RCurl_1.98-1.16
## [13] QuasR_1.47.0 ggpointdensity_0.1.0
## [15] htmltools_0.5.8.1 S4Arrays_1.7.1
## [17] progress_1.2.3 curl_6.0.1
## [19] SparseArray_1.7.2 sass_0.4.9
## [21] bslib_0.8.0 plyr_1.8.9
## [23] httr2_1.0.7 cachem_1.1.0
## [25] buildtools_1.0.0 GenomicAlignments_1.43.0
## [27] lifecycle_1.0.4 pkgconfig_2.0.3
## [29] Matrix_1.7-1 R6_2.5.1
## [31] fastmap_1.2.0 GenomeInfoDbData_1.2.13
## [33] MatrixGenerics_1.19.0 digest_0.6.37
## [35] colorspace_2.1-1 ShortRead_1.65.0
## [37] patchwork_1.3.0 AnnotationDbi_1.69.0
## [39] RSQLite_2.3.8 hwriter_1.3.2.1
## [41] labeling_0.4.3 filelock_1.0.3
## [43] fansi_1.0.6 httr_1.4.7
## [45] abind_1.4-8 compiler_4.4.2
## [47] Rbowtie_1.47.0 withr_3.0.2
## [49] bit64_4.5.2 BiocParallel_1.41.0
## [51] viridis_0.6.5 DBI_1.2.3
## [53] qs_0.27.2 biomaRt_2.63.0
## [55] rappdirs_0.3.3 DelayedArray_0.33.2
## [57] rjson_0.2.23 tools_4.4.2
## [59] glue_1.8.0 restfulr_0.0.15
## [61] grid_4.4.2 gtable_0.3.6
## [63] tidyr_1.3.1 RApiSerialize_0.1.4
## [65] hms_1.1.3 stringfish_0.16.0
## [67] xml2_1.3.6 utf8_1.2.4
## [69] ggrepel_0.9.6 pillar_1.9.0
## [71] stringr_1.5.1 dplyr_1.1.4
## [73] BiocFileCache_2.15.0 lattice_0.22-6
## [75] bit_4.5.0 deldir_2.0-4
## [77] tidyselect_1.2.1 maketools_1.3.1
## [79] knitr_1.49 gridExtra_2.3
## [81] SummarizedExperiment_1.37.0 xfun_0.49
## [83] Biobase_2.67.0 matrixStats_1.4.1
## [85] stringi_1.8.4 UCSC.utils_1.3.0
## [87] yaml_2.3.10 evaluate_1.0.1
## [89] codetools_0.2-20 interp_1.1-6
## [91] GenomicFiles_1.43.0 tibble_3.2.1
## [93] cli_3.6.3 RcppParallel_5.1.9
## [95] munsell_0.5.1 jquerylib_0.1.4
## [97] Rcpp_1.0.13-1 dbplyr_2.5.0
## [99] tidyverse_2.0.0 png_0.1-8
## [101] XML_3.99-0.17 parallel_4.4.2
## [103] ggplot2_3.5.1 blob_1.2.4
## [105] prettyunits_1.2.0 latticeExtra_0.6-30
## [107] jpeg_0.1-10 plyranges_1.27.0
## [109] bitops_1.0-9 pwalign_1.3.0
## [111] txdbmaker_1.3.1 viridisLite_0.4.2
## [113] VariantAnnotation_1.53.0 scales_1.3.0
## [115] purrr_1.0.2 crayon_1.5.3
## [117] rlang_1.1.4 KEGGREST_1.47.0