library(ggplot2)
library(purrr)
library(GenomicRanges)
## Loading required package: stats4
## Loading required package: BiocGenerics
##
## Attaching package: 'BiocGenerics'
## The following objects are masked from 'package:stats':
##
## IQR, mad, sd, var, xtabs
## The following objects are masked from 'package:base':
##
## Filter, Find, Map, Position, Reduce, anyDuplicated, aperm,
## append, as.data.frame, basename, cbind, colnames, dirname,
## do.call, duplicated, eval, evalq, get, grep, grepl, intersect,
## is.unsorted, lapply, mapply, match, mget, order, paste, pmax,
## pmax.int, pmin, pmin.int, rank, rbind, rownames, sapply,
## saveRDS, setdiff, table, tapply, union, unique, unsplit,
## which.max, which.min
## Loading required package: S4Vectors
##
## Attaching package: 'S4Vectors'
## The following object is masked from 'package:utils':
##
## findMatches
## The following objects are masked from 'package:base':
##
## I, expand.grid, unname
## Loading required package: IRanges
##
## Attaching package: 'IRanges'
## The following object is masked from 'package:purrr':
##
## reduce
## Loading required package: GenomeInfoDb
library(InteractionSet)
## Loading required package: SummarizedExperiment
## Loading required package: MatrixGenerics
## Loading required package: matrixStats
##
## Attaching package: 'MatrixGenerics'
## The following objects are masked from 'package:matrixStats':
##
## colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse,
## colCounts, colCummaxs, colCummins, colCumprods, colCumsums,
## colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs,
## colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats,
## colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds,
## colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads,
## colWeightedMeans, colWeightedMedians, colWeightedSds,
## colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet,
## rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods,
## rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps,
## rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins,
## rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks,
## rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs,
## rowVars, rowWeightedMads, rowWeightedMeans, rowWeightedMedians,
## rowWeightedSds, rowWeightedVars
## Loading required package: Biobase
## Welcome to Bioconductor
##
## Vignettes contain introductory material; view with
## 'browseVignettes()'. To cite Bioconductor, see
## 'citation("Biobase")', and for packages 'citation("pkgname")'.
##
## Attaching package: 'Biobase'
## The following object is masked from 'package:MatrixGenerics':
##
## rowMedians
## The following objects are masked from 'package:matrixStats':
##
## anyMissing, rowMedians
library(HiCExperiment)
## Consider using the `HiContacts` package to perform advanced genomic operations
## on `HiCExperiment` objects.
##
## Read "Orchestrating Hi-C analysis with Bioconductor" online book to learn more:
## https://js2264.github.io/OHCA/
##
## Attaching package: 'HiCExperiment'
## The following object is masked from 'package:SummarizedExperiment':
##
## metadata<-
## The following object is masked from 'package:S4Vectors':
##
## metadata<-
## The following object is masked from 'package:ggplot2':
##
## resolution
library(HiContactsData)
## Loading required package: ExperimentHub
## Loading required package: AnnotationHub
## Loading required package: BiocFileCache
## Loading required package: dbplyr
##
## Attaching package: 'AnnotationHub'
## The following object is masked from 'package:Biobase':
##
## cache
library(multiHiCcompare)
##
## Attaching package: 'multiHiCcompare'
## The following object is masked from 'package:HiCExperiment':
##
## resolution
## The following object is masked from 'package:ggplot2':
##
## resolutionWorkflow 3: Inter-centromere interactions in yeast
Pre-loading packages and objects 📦
Aims
This chapter illustrates how to plot the aggregate signal over pairs of genomic ranges, in this case pairs of yeast centromeres.
Datasets
We leverage two yeast datasets in this notebook.
- One from a WT yeast strain in G1 phase
- One from a WT yeast strain in G2/M phase
Importing Hi-C data and plotting contact matrices
library(HiContactsData)
library(HiContacts)
## Registered S3 methods overwritten by 'readr':
## method from
## as.data.frame.spec_tbl_df vroom
## as_tibble.spec_tbl_df vroom
## format.col_spec vroom
## print.col_spec vroom
## print.collector vroom
## print.date_names vroom
## print.locale vroom
## str.col_spec vroom
library(purrr)
library(ggplot2)
hics <- list(
'G1' = import(HiContactsData('yeast_g1', 'mcool'), format = 'cool', resolution = 4000),
'G2M' = import(HiContactsData('yeast_g2m', 'mcool'), format = 'cool', resolution = 4000)
)
## see ?HiContactsData and browseVignettes('HiContactsData') for documentation
## loading from cache
## see ?HiContactsData and browseVignettes('HiContactsData') for documentation
## loading from cache
imap(hics, ~ plotMatrix(
.x, use.scores = 'balanced', limits = c(-4, -1), caption = FALSE
) + ggtitle(.y))
## $G1
##
## $G2M
We can visually appreciate that inter-chromosomal interactions, notably between centromeres, are less prominent in G2/M.
Checking P(s) and cis/trans interactions ratio
library(dplyr)
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:dbplyr':
##
## ident, sql
## The following object is masked from 'package:Biobase':
##
## combine
## The following object is masked from 'package:matrixStats':
##
## count
## The following objects are masked from 'package:GenomicRanges':
##
## intersect, setdiff, union
## The following object is masked from 'package:GenomeInfoDb':
##
## intersect
## The following objects are masked from 'package:IRanges':
##
## collapse, desc, intersect, setdiff, slice, union
## The following objects are masked from 'package:S4Vectors':
##
## first, intersect, rename, setdiff, setequal, union
## The following objects are masked from 'package:BiocGenerics':
##
## combine, intersect, setdiff, union
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
pairs <- list(
'G1' = PairsFile(HiContactsData('yeast_g1', 'pairs')),
'G2M' = PairsFile(HiContactsData('yeast_g2m', 'pairs'))
)
## see ?HiContactsData and browseVignettes('HiContactsData') for documentation
## loading from cache
## see ?HiContactsData and browseVignettes('HiContactsData') for documentation
## loading from cache
ps <- imap_dfr(pairs, ~ distanceLaw(.x, by_chr = TRUE) |>
mutate(sample = .y)
)
## Importing pairs file /home/biocbuild/.cache/R/ExperimentHub/195d383c745512_8630 in memory. This may take a while...
## Importing pairs file /home/biocbuild/.cache/R/ExperimentHub/195d3854b3f765_8631 in memory. This may take a while...
plotPs(ps, aes(x = binned_distance, y = norm_p, group = interaction(sample, chr), color = sample)) +
scale_color_manual(values = c('black', 'red'))
## Warning: Removed 2133 rows containing missing values or values outside the scale
## range (`geom_line()`).
plotPsSlope(ps, ggplot2::aes(x = binned_distance, y = slope, group = interaction(sample, chr), color = sample)) +
scale_color_manual(values = c('black', 'red'))
## Warning: Removed 2183 rows containing missing values or values outside the scale
## range (`geom_line()`).
This confirms that interactions in cells synchronized in G2/M are enriched for 10-30kb-long interactions.
ratios <- imap_dfr(hics, ~ cisTransRatio(.x) |> mutate(sample = .y))
ggplot(ratios, aes(x = chr, y = trans_pct, fill = sample)) +
geom_col() +
labs(x = 'Chromosomes', y = "% of trans interactions") +
scale_y_continuous(labels = scales::percent) +
facet_grid(~sample)
We can also highlight that trans (inter-chromosomal) interactions are proportionally decreasing in G2/M-synchronized cells.
Centromere virtual 4C profiles
data(centros_yeast)
v4c_centro <- imap_dfr(hics, ~ virtual4C(.x, GenomicRanges::resize(centros_yeast[2], 8000)) |>
as_tibble() |>
mutate(sample = .y) |>
filter(seqnames == 'IV')
)
ggplot(v4c_centro, aes(x = start, y = score, fill = sample)) +
geom_area() +
theme_bw() +
labs(
x = "chrIV position",
y = "Contacts with chrII centromere",
title = "Interaction profile of chrII centromere"
) +
coord_cartesian(ylim = c(0, 0.015))
Aggregated 2D signal over all pairs of centromeres
We can start by computing all possible pairs of centromeres.
centros_pairs <- lapply(1:length(centros_yeast), function(i) {
lapply(1:length(centros_yeast), function(j) {
S4Vectors::Pairs(centros_yeast[i], centros_yeast[j])
})
}) |>
do.call(c, args = _) |>
do.call(c, args = _) |>
InteractionSet::makeGInteractionsFromGRangesPairs()
centros_pairs <- centros_pairs[anchors(centros_pairs, 'first') != anchors(centros_pairs, 'second')]
centros_pairs
## GInteractions object with 240 interactions and 0 metadata columns:
## seqnames1 ranges1 seqnames2 ranges2
## <Rle> <IRanges> <Rle> <IRanges>
## [1] I 151583-151641 --- II 238361-238419
## [2] I 151583-151641 --- III 114322-114380
## [3] I 151583-151641 --- IV 449879-449937
## [4] I 151583-151641 --- V 152522-152580
## [5] I 151583-151641 --- VI 147981-148039
## ... ... ... ... ... ...
## [236] XVI 556255-556313 --- XI 440229-440287
## [237] XVI 556255-556313 --- XII 151366-151424
## [238] XVI 556255-556313 --- XIII 268222-268280
## [239] XVI 556255-556313 --- XIV 628588-628646
## [240] XVI 556255-556313 --- XV 326897-326955
## -------
## regions: 16 ranges and 0 metadata columns
## seqinfo: 17 sequences (1 circular) from R64-1-1 genomeThen we can aggregate the Hi-C signal over each pair of centromeres.
aggr_maps <- purrr::imap(hics, ~ {
aggr <- aggregate(.x, centros_pairs, maxDistance = 1e999)
plotMatrix(
aggr, use.scores = 'balanced', limits = c(-5, -1),
cmap = HiContacts::rainbowColors(),
caption = FALSE
) + ggtitle(.y)
})
## Warning in valid.GenomicRanges.seqinfo(x, suggest.trim = TRUE): GRanges object contains 120 out-of-bound ranges located on sequences
## I, III, V, VI, VIII, IX, XII, and XIV. Note that ranges located on a
## sequence whose length is unknown (NA) or on a circular sequence are
## not considered out-of-bound (use seqlengths() and isCircular() to
## get the lengths and circularity flags of the underlying sequences).
## You can use trim() to trim these ranges. See
## ?`trim,GenomicRanges-method` for more information.
## Warning in valid.GenomicRanges.seqinfo(x, suggest.trim = TRUE): GRanges object contains 120 out-of-bound ranges located on sequences
## III, V, VI, VIII, IX, XII, XIV, and I. Note that ranges located on a
## sequence whose length is unknown (NA) or on a circular sequence are
## not considered out-of-bound (use seqlengths() and isCircular() to
## get the lengths and circularity flags of the underlying sequences).
## You can use trim() to trim these ranges. See
## ?`trim,GenomicRanges-method` for more information.
## Warning in valid.GenomicRanges.seqinfo(x, suggest.trim = TRUE): GRanges object contains 240 out-of-bound ranges located on sequences
## I, III, V, VI, VIII, IX, XII, and XIV. Note that ranges located on a
## sequence whose length is unknown (NA) or on a circular sequence are
## not considered out-of-bound (use seqlengths() and isCircular() to
## get the lengths and circularity flags of the underlying sequences).
## You can use trim() to trim these ranges. See
## ?`trim,GenomicRanges-method` for more information.
## Going through preflight checklist...
## Parsing the entire contact matrice as a sparse matrix...
## Modeling distance decay...
## Filtering for contacts within provided targets...
## Warning in valid.GenomicRanges.seqinfo(x, suggest.trim = TRUE): GRanges object contains 120 out-of-bound ranges located on sequences
## I, III, V, VI, VIII, IX, XII, and XIV. Note that ranges located on a
## sequence whose length is unknown (NA) or on a circular sequence are
## not considered out-of-bound (use seqlengths() and isCircular() to
## get the lengths and circularity flags of the underlying sequences).
## You can use trim() to trim these ranges. See
## ?`trim,GenomicRanges-method` for more information.
## Warning in valid.GenomicRanges.seqinfo(x, suggest.trim = TRUE): GRanges object contains 120 out-of-bound ranges located on sequences
## III, V, VI, VIII, IX, XII, XIV, and I. Note that ranges located on a
## sequence whose length is unknown (NA) or on a circular sequence are
## not considered out-of-bound (use seqlengths() and isCircular() to
## get the lengths and circularity flags of the underlying sequences).
## You can use trim() to trim these ranges. See
## ?`trim,GenomicRanges-method` for more information.
## Warning in valid.GenomicRanges.seqinfo(x, suggest.trim = TRUE): GRanges object contains 240 out-of-bound ranges located on sequences
## I, III, V, VI, VIII, IX, XII, and XIV. Note that ranges located on a
## sequence whose length is unknown (NA) or on a circular sequence are
## not considered out-of-bound (use seqlengths() and isCircular() to
## get the lengths and circularity flags of the underlying sequences).
## You can use trim() to trim these ranges. See
## ?`trim,GenomicRanges-method` for more information.
## Going through preflight checklist...
## Parsing the entire contact matrice as a sparse matrix...
## Modeling distance decay...
## Filtering for contacts within provided targets...
cowplot::plot_grid(plotlist = aggr_maps, nrow = 1)
Aggregated 1D interaction profile of centromeres
One can generalize the previous virtual 4C plot, by extracting the interaction profile between all possible pairs of centromeres in each dataset.
df <- map_dfr(1:{length(centros_yeast)-1}, function(i) {
centro1 <- GenomicRanges::resize(centros_yeast[i], fix = 'center', 8000)
map_dfr({i+1}:length(centros_yeast), function(j) {
centro2 <- GenomicRanges::resize(centros_yeast[j], fix = 'center', 80000)
gi <- InteractionSet::GInteractions(centro1, centro2)
imap_dfr(hics, ~ .x[gi] |>
interactions() |>
as_tibble() |>
mutate(
sample = .y,
center = center2 - start(GenomicRanges::resize(centro2, fix = 'center', 1))
) |>
select(sample, seqnames1, seqnames2, center, balanced)
)
})
})
ggplot(df, aes(x = center/1e3, y = balanced)) +
geom_line(aes(group = interaction(seqnames1, seqnames2)), alpha = 0.03, col = "black") +
geom_smooth(col = "red", fill = "red") +
theme_bw() +
theme(legend.position = 'none') +
labs(
x = "Distance from centromere (kb)", y = "Normalized interaction frequency",
title = "Centromere pairwise interaction profiles"
) +
facet_grid(~sample)
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: Removed 25 rows containing non-finite outside the scale range
## (`stat_smooth()`).
## Warning: Removed 3 rows containing missing values or values outside the scale range
## (`geom_line()`).
Session info
Click to expand 👇
sessioninfo::session_info(include_base = TRUE)
## ─ Session info ────────────────────────────────────────────────────────────
## setting value
## version R version 4.4.1 (2024-06-14)
## os Ubuntu 24.04.1 LTS
## system x86_64, linux-gnu
## ui X11
## language (EN)
## collate C
## ctype en_US.UTF-8
## tz America/New_York
## date 2024-10-30
## pandoc 2.7.3 @ /usr/bin/ (via rmarkdown)
##
## ─ Packages ────────────────────────────────────────────────────────────────
## package * version date (UTC) lib source
## abind 1.4-8 2024-09-12 [2] CRAN (R 4.4.1)
## aggregation 1.0.1 2018-01-25 [2] CRAN (R 4.4.1)
## AnnotationDbi 1.68.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## AnnotationHub * 3.14.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## base * 4.4.1 2024-09-25 [3] local
## beeswarm 0.4.0 2021-06-01 [2] CRAN (R 4.4.1)
## Biobase * 2.66.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## BiocFileCache * 2.14.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## BiocGenerics * 0.52.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## BiocIO 1.16.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## BiocManager 1.30.25 2024-08-28 [2] CRAN (R 4.4.1)
## BiocParallel 1.40.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## BiocVersion 3.20.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## Biostrings 2.74.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## bit 4.5.0 2024-09-20 [2] CRAN (R 4.4.1)
## bit64 4.5.2 2024-09-22 [2] CRAN (R 4.4.1)
## blob 1.2.4 2023-03-17 [2] CRAN (R 4.4.1)
## cachem 1.1.0 2024-05-16 [2] CRAN (R 4.4.1)
## Cairo 1.6-2 2023-11-28 [2] CRAN (R 4.4.1)
## calibrate 1.7.7 2020-06-19 [2] CRAN (R 4.4.1)
## cli 3.6.3 2024-06-21 [2] CRAN (R 4.4.1)
## codetools 0.2-20 2024-03-31 [3] CRAN (R 4.4.1)
## colorspace 2.1-1 2024-07-26 [2] CRAN (R 4.4.1)
## compiler 4.4.1 2024-09-25 [3] local
## cowplot 1.1.3 2024-01-22 [2] CRAN (R 4.4.1)
## crayon 1.5.3 2024-06-20 [2] CRAN (R 4.4.1)
## curl 5.2.3 2024-09-20 [2] CRAN (R 4.4.1)
## data.table 1.16.2 2024-10-10 [2] CRAN (R 4.4.1)
## datasets * 4.4.1 2024-09-25 [3] local
## DBI 1.2.3 2024-06-02 [2] CRAN (R 4.4.1)
## dbplyr * 2.5.0 2024-03-19 [2] CRAN (R 4.4.1)
## DelayedArray 0.32.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## digest 0.6.37 2024-08-19 [2] CRAN (R 4.4.1)
## dplyr * 1.1.4 2023-11-17 [2] CRAN (R 4.4.1)
## edgeR 4.4.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## evaluate 1.0.1 2024-10-10 [2] CRAN (R 4.4.1)
## ExperimentHub * 2.14.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## fansi 1.0.6 2023-12-08 [2] CRAN (R 4.4.1)
## farver 2.1.2 2024-05-13 [2] CRAN (R 4.4.1)
## fastmap 1.2.0 2024-05-15 [2] CRAN (R 4.4.1)
## filelock 1.0.3 2023-12-11 [2] CRAN (R 4.4.1)
## generics 0.1.3 2022-07-05 [2] CRAN (R 4.4.1)
## GenomeInfoDb * 1.42.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## GenomeInfoDbData 1.2.13 2024-10-01 [2] Bioconductor
## GenomicRanges * 1.58.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## ggbeeswarm 0.7.2 2023-04-29 [2] CRAN (R 4.4.1)
## ggplot2 * 3.5.1 2024-04-23 [2] CRAN (R 4.4.1)
## ggrastr 1.0.2 2023-06-01 [2] CRAN (R 4.4.1)
## glue 1.8.0 2024-09-30 [2] CRAN (R 4.4.1)
## graphics * 4.4.1 2024-09-25 [3] local
## grDevices * 4.4.1 2024-09-25 [3] local
## grid 4.4.1 2024-09-25 [3] local
## gridExtra 2.3 2017-09-09 [2] CRAN (R 4.4.1)
## gtable 0.3.6 2024-10-25 [2] CRAN (R 4.4.1)
## gtools 3.9.5 2023-11-20 [2] CRAN (R 4.4.1)
## HiCcompare 1.28.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## HiCExperiment * 1.6.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## HiContacts * 1.8.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## HiContactsData * 1.7.0 2024-10-26 [2] Bioconductor 3.20 (R 4.4.1)
## hms 1.1.3 2023-03-21 [2] CRAN (R 4.4.1)
## htmltools 0.5.8.1 2024-04-04 [2] CRAN (R 4.4.1)
## htmlwidgets 1.6.4 2023-12-06 [2] CRAN (R 4.4.1)
## httr 1.4.7 2023-08-15 [2] CRAN (R 4.4.1)
## InteractionSet * 1.34.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## IRanges * 2.40.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## jsonlite 1.8.9 2024-09-20 [2] CRAN (R 4.4.1)
## KEGGREST 1.46.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## KernSmooth 2.23-24 2024-05-17 [3] CRAN (R 4.4.1)
## knitr 1.48 2024-07-07 [2] CRAN (R 4.4.1)
## labeling 0.4.3 2023-08-29 [2] CRAN (R 4.4.1)
## lattice 0.22-6 2024-03-20 [3] CRAN (R 4.4.1)
## lifecycle 1.0.4 2023-11-07 [2] CRAN (R 4.4.1)
## limma 3.62.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## locfit 1.5-9.10 2024-06-24 [2] CRAN (R 4.4.1)
## magrittr 2.0.3 2022-03-30 [2] CRAN (R 4.4.1)
## MASS 7.3-61 2024-06-13 [3] CRAN (R 4.4.1)
## Matrix 1.7-1 2024-10-18 [3] CRAN (R 4.4.1)
## MatrixGenerics * 1.18.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## matrixStats * 1.4.1 2024-09-08 [2] CRAN (R 4.4.1)
## memoise 2.0.1 2021-11-26 [2] CRAN (R 4.4.1)
## methods * 4.4.1 2024-09-25 [3] local
## mgcv 1.9-1 2023-12-21 [3] CRAN (R 4.4.1)
## mime 0.12 2021-09-28 [2] CRAN (R 4.4.1)
## multiHiCcompare * 1.24.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## munsell 0.5.1 2024-04-01 [2] CRAN (R 4.4.1)
## nlme 3.1-166 2024-08-14 [3] CRAN (R 4.4.1)
## parallel 4.4.1 2024-09-25 [3] local
## pbapply 1.7-2 2023-06-27 [2] CRAN (R 4.4.1)
## pheatmap 1.0.12 2019-01-04 [2] CRAN (R 4.4.1)
## pillar 1.9.0 2023-03-22 [2] CRAN (R 4.4.1)
## pkgconfig 2.0.3 2019-09-22 [2] CRAN (R 4.4.1)
## png 0.1-8 2022-11-29 [2] CRAN (R 4.4.1)
## purrr * 1.0.2 2023-08-10 [2] CRAN (R 4.4.1)
## qqman 0.1.9 2023-08-23 [2] CRAN (R 4.4.1)
## R6 2.5.1 2021-08-19 [2] CRAN (R 4.4.1)
## rappdirs 0.3.3 2021-01-31 [2] CRAN (R 4.4.1)
## RColorBrewer 1.1-3 2022-04-03 [2] CRAN (R 4.4.1)
## Rcpp 1.0.13 2024-07-17 [2] CRAN (R 4.4.1)
## readr 2.1.5 2024-01-10 [2] CRAN (R 4.4.1)
## rhdf5 2.50.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## rhdf5filters 1.18.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## Rhdf5lib 1.28.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## rlang 1.1.4 2024-06-04 [2] CRAN (R 4.4.1)
## rmarkdown 2.28 2024-08-17 [2] CRAN (R 4.4.1)
## RSpectra 0.16-2 2024-07-18 [2] CRAN (R 4.4.1)
## RSQLite 2.3.7 2024-05-27 [2] CRAN (R 4.4.1)
## S4Arrays 1.6.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## S4Vectors * 0.44.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## scales 1.3.0 2023-11-28 [2] CRAN (R 4.4.1)
## sessioninfo 1.2.2 2021-12-06 [2] CRAN (R 4.4.1)
## SparseArray 1.6.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## splines 4.4.1 2024-09-25 [3] local
## statmod 1.5.0 2023-01-06 [2] CRAN (R 4.4.1)
## stats * 4.4.1 2024-09-25 [3] local
## stats4 * 4.4.1 2024-09-25 [3] local
## strawr 0.0.92 2024-07-16 [2] CRAN (R 4.4.1)
## stringi 1.8.4 2024-05-06 [2] CRAN (R 4.4.1)
## stringr 1.5.1 2023-11-14 [2] CRAN (R 4.4.1)
## SummarizedExperiment * 1.36.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## tibble 3.2.1 2023-03-20 [2] CRAN (R 4.4.1)
## tidyr 1.3.1 2024-01-24 [2] CRAN (R 4.4.1)
## tidyselect 1.2.1 2024-03-11 [2] CRAN (R 4.4.1)
## tools 4.4.1 2024-09-25 [3] local
## tzdb 0.4.0 2023-05-12 [2] CRAN (R 4.4.1)
## UCSC.utils 1.2.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## utf8 1.2.4 2023-10-22 [2] CRAN (R 4.4.1)
## utils * 4.4.1 2024-09-25 [3] local
## vctrs 0.6.5 2023-12-01 [2] CRAN (R 4.4.1)
## vipor 0.4.7 2023-12-18 [2] CRAN (R 4.4.1)
## vroom 1.6.5 2023-12-05 [2] CRAN (R 4.4.1)
## withr 3.0.2 2024-10-28 [2] CRAN (R 4.4.1)
## xfun 0.48 2024-10-03 [2] CRAN (R 4.4.1)
## XVector 0.46.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
## yaml 2.3.10 2024-07-26 [2] CRAN (R 4.4.1)
## zlibbioc 1.52.0 2024-10-29 [2] Bioconductor 3.20 (R 4.4.1)
##
## [1] /home/biocbuild/bbs-3.20-books/tmpdir/RtmphzOvfD/Rinst2020235c565b68
## [2] /home/biocbuild/bbs-3.20-bioc/R/site-library
## [3] /home/biocbuild/bbs-3.20-bioc/R/library
##
## ───────────────────────────────────────────────────────────────────────────