1 What is Spatial Registration?

Spatial Registration is an analysis that compares the gene expression of groups in a query RNA-seq data set (typically spatially resolved RNA-seq or single cell RNA-seq) to groups in a reference spatially resolved RNA-seq data set (such annotated anatomical features).

For spatial data, this can be helpful to compare manual annotations, or annotating clusters. For scRNA-seq data it can check if a cell type might be more concentrated in one area or anatomical feature of the tissue.

The spatial annotation process correlates the \(t\)-statistics from the gene enrichment analysis between spatial features from the reference data set, with the \(t\)-statistics from the gene enrichment of features in the query data set. Pairs with high positive correlation show where similar patterns of gene expression are occurring and what anatomical feature the new spatial feature or cell population may map to.

1.1 Overview of the Spatial Registration method

  1. Perform gene set enrichment analysis between spatial features (ex. anatomical features, histological layers) on reference spatial data set. Or access existing statistics.

  2. Perform gene set enrichment analysis between features (ex. new annotations, data-driven clusters) on new query data set.

  3. Correlate the \(t\)-statistics between the reference and query features.

  4. Annotate new spatial features with the most strongly associated reference feature.

  5. Plot correlation heat map to observe patterns between the two data sets.

Spatial Registration Overview

2 How to run Spatial Registration with spatialLIBD tools

2.1 Introduction.

In this example we will utilize the human DLPFC 10x Genomics Visium dataset from Maynard, Collado-Torres et al. (Maynard, Collado-Torres, Weber et al., 2021) as the reference. This data contains manually annotated features: the six cortical layers + white matter present in the DLPFC. We will use the pre-calculated enrichment statistics for the layers, which are available from spatialLIBD.

Dotplot of sample from refernce DLPFC data, colored by annotated layers

The query dataset will be the DLPFC single nucleus RNA-seq (snRNA-seq) data from (Tran, Maynard, Spangler, Huuki, Montgomery, Sadashivaiah, Tippani, Barry, Hancock, Hicks, Kleinman, Hyde, Collado-Torres, Jaffe, and Martinowich, 2021).

We will compare the gene expression in the cell type populations of the query dataset to the annotated layers in the reference.

2.2 Important Notes

2.2.1 Required knowledge

It may be helpful to review Introduction to spatialLIBD vignette available through GitHub or Bioconductor for more information about this data set and R package.

2.2.2 Citing spatialLIBD

We hope that spatialLIBD will be useful for your research. Please use the following information to cite the package and the overall approach. Thank you!

## Citation info
citation("spatialLIBD")
#> To cite package 'spatialLIBD' in publications use:
#> 
#>   Pardo B, Spangler A, Weber LM, Hicks SC, Jaffe AE, Martinowich K, Maynard KR, Collado-Torres L (2022).
#>   "spatialLIBD: an R/Bioconductor package to visualize spatially-resolved transcriptomics data." _BMC
#>   Genomics_. doi:10.1186/s12864-022-08601-w <https://doi.org/10.1186/s12864-022-08601-w>,
#>   <https://doi.org/10.1186/s12864-022-08601-w>.
#> 
#>   Maynard KR, Collado-Torres L, Weber LM, Uytingco C, Barry BK, Williams SR, II JLC, Tran MN, Besich Z,
#>   Tippani M, Chew J, Yin Y, Kleinman JE, Hyde TM, Rao N, Hicks SC, Martinowich K, Jaffe AE (2021).
#>   "Transcriptome-scale spatial gene expression in the human dorsolateral prefrontal cortex." _Nature
#>   Neuroscience_. doi:10.1038/s41593-020-00787-0 <https://doi.org/10.1038/s41593-020-00787-0>,
#>   <https://www.nature.com/articles/s41593-020-00787-0>.
#> 
#>   Huuki-Myers LA, Spangler A, Eagles NJ, Montgomergy KD, Kwon SH, Guo B, Grant-Peters M, Divecha HR,
#>   Tippani M, Sriworarat C, Nguyen AB, Ravichandran P, Tran MN, Seyedian A, Consortium P, Hyde TM, Kleinman
#>   JE, Battle A, Page SC, Ryten M, Hicks SC, Martinowich K, Collado-Torres L, Maynard KR (2023). "Integrated
#>   single cell and unsupervised spatial transcriptomic analysis defines molecular anatomy of the human
#>   dorsolateral prefrontal cortex." _bioRxiv_. doi:10.1101/2023.02.15.528722
#>   <https://doi.org/10.1101/2023.02.15.528722>,
#>   <https://www.biorxiv.org/content/10.1101/2023.02.15.528722v1>.
#> 
#>   Kwon SH, Parthiban S, Tippani M, Divecha HR, Eagles NJ, Lobana JS, Williams SR, Mark M, Bharadwaj RA,
#>   Kleinman JE, Hyde TM, Page SC, Hicks SC, Martinowich K, Maynard KR, Collado-Torres L (2023). "Influence
#>   of Alzheimer’s disease related neuropathology on local microenvironment gene expression in the human
#>   inferior temporal cortex." _bioRxiv_. doi:10.1101/2023.04.20.537710
#>   <https://doi.org/10.1101/2023.04.20.537710>,
#>   <https://www.biorxiv.org/content/10.1101/2023.04.20.537710v1>.
#> 
#> To see these entries in BibTeX format, use 'print(<citation>, bibtex=TRUE)', 'toBibtex(.)', or set
#> 'options(citation.bibtex.max=999)'.

2.3 Setup

2.3.1 Install spatialLIBD

if (!requireNamespace("BiocManager", quietly = TRUE)) {
    install.packages("BiocManager")
}

BiocManager::install("spatialLIBD")

## Check that you have a valid Bioconductor installation
BiocManager::valid()

2.3.2 Load required packages

library("spatialLIBD")
library("SingleCellExperiment")

2.4 Download Data

2.4.1 Spatial Reference

The reference data is easily accessed through spatialLIBD. The modeling results for the annotated layers is already calculated and can be accessed with the fetch_data() function.

This data contains the results form three models (anova, enrichment, and pairwise), we will use the enrichment results for spatial registration. The tables contain the \(t\)-statistics, p-values, and gene ensembl ID and symbol.

## get reference layer enrichment statistics
layer_modeling_results <- fetch_data(type = "modeling_results")
#> 2023-11-30 10:57:28.223936 loading file /home/biocbuild/.cache/R/BiocFileCache/ad44f1d390e83_Human_DLPFC_Visium_modeling_results.Rdata%3Fdl%3D1

layer_modeling_results$enrichment[1:5, 1:5]
#>    t_stat_WM t_stat_Layer1 t_stat_Layer2 t_stat_Layer3 t_stat_Layer4
#> 1 -0.6344143    -1.0321320     0.1781501   -0.72835965    1.56703859
#> 2 -2.4758891     1.2232062    -0.8733745    1.93793650    1.33150141
#> 3 -3.0079360    -0.8564572     2.1335852    0.48741121    0.35212807
#> 4 -1.2916584    -0.9494234    -0.9485440    0.56378302   -0.11206713
#> 5  2.3175897     0.6156900     0.1127478   -0.09907566   -0.03376771

2.4.2 Query Data: snRNA-seq

For the query data set, we will use the public single nucleus RNA-seq (snRNA-seq) data from (Tran, Maynard, Spangler et al., 2021) can be accessed on github.

This data is also from postmortem human brain DLPFC, and contains gene expression data for 11k nuclei and 19 cell types.

We will use BiocFileCache() to cache this data. It is stored as a SingleCellExperiment object named sce.dlpfc.tran, and takes 1.01 GB of RAM memory to load.

# Download and save a local cache of the data available at:
# https://github.com/LieberInstitute/10xPilot_snRNAseq-human#processed-data
bfc <- BiocFileCache::BiocFileCache()
url <- paste0(
    "https://libd-snrnaseq-pilot.s3.us-east-2.amazonaws.com/",
    "SCE_DLPFC-n3_tran-etal.rda"
)
local_data <- BiocFileCache::bfcrpath(url, x = bfc)

load(local_data, verbose = TRUE)
#> Loading objects:
#>   sce.dlpfc.tran

DLPFC tissue consists of many cell types, some are quite rare and will not have enough data to complete the analysis

table(sce.dlpfc.tran$cellType)
#> 
#>      Astro    Excit_A    Excit_B    Excit_C    Excit_D    Excit_E    Excit_F    Inhib_A    Inhib_B    Inhib_C 
#>        782        529        773        524        132        187        243        333        454        365 
#>    Inhib_D    Inhib_E    Inhib_F Macrophage      Micro      Mural      Oligo        OPC      Tcell 
#>        413          7          8         10        388         18       5455        572          9

The data will be pseudo-bulked over donor x cellType, it is recommended to drop groups with < 10 nuclei (this is done automatically in the pseudobulk step).

table(sce.dlpfc.tran$donor, sce.dlpfc.tran$cellType)
#>         
#>          Astro Excit_A Excit_B Excit_C Excit_D Excit_E Excit_F Inhib_A Inhib_B Inhib_C Inhib_D Inhib_E Inhib_F
#>   donor1   371     111      75      44      22      77     102      39      98      47     119       2       0
#>   donor2   137     120     154     155      27      25      36      89     106      56      78       2       1
#>   donor6   274     298     544     325      83      85     105     205     250     262     216       3       7
#>         
#>          Macrophage Micro Mural Oligo  OPC Tcell
#>   donor1          1   152     3  2754  196     2
#>   donor2          3    92     2   517   91     2
#>   donor6          6   144    13  2184  285     5

2.5 Get Enrichment statistics for snRNA-seq data

spatialLIBD contains many functions to compute modeling_results for the query sc/snRNA-seq or spatial data.

The process includes the following steps

  1. registration_pseudobulk(): Pseudo-bulks data, filter low expressed genes, and normalize counts
  2. registration_mod(): Defines the statistical model that will be used for computing the block correlation
  3. registration_block_cor() : Computes the block correlation using the sample ID as the blocking factor, used as correlation in eBayes call
  4. registration_stats_enrichment() : Computes the gene enrichment \(t\)-statistics (one group vs. All other groups)

The function registration_wrapper() makes life easy by wrapping these functions together in to one step!

## Perform the spatial registration
sce_modeling_results <- registration_wrapper(
    sce = sce.dlpfc.tran,
    var_registration = "cellType",
    var_sample_id = "donor",
    gene_ensembl = "gene_id",
    gene_name = "gene_name"
)
#> 2023-11-30 10:57:37.306299 make pseudobulk object
#> 2023-11-30 10:57:40.867953 dropping 13 pseudo-bulked samples that are below 'min_ncells'.
#> 2023-11-30 10:57:40.932733 drop lowly expressed genes
#> 2023-11-30 10:57:41.105781 normalize expression
#> 2023-11-30 10:57:43.39106 create model matrix
#> 2023-11-30 10:57:43.427891 run duplicateCorrelation()
#> 2023-11-30 10:58:06.293711 The estimated correlation is: 0.138734774807097
#> 2023-11-30 10:58:06.298316 computing enrichment statistics
#> 2023-11-30 10:58:08.408998 extract and reformat enrichment results
#> 2023-11-30 10:58:08.614462 running the baseline pairwise model
#> 2023-11-30 10:58:08.869322 computing pairwise statistics
#> 2023-11-30 10:58:12.720189 computing F-statistics

2.6 Extract Enrichment t-statistics

## extract t-statics and rename
registration_t_stats <- sce_modeling_results$enrichment[, grep("^t_stat", colnames(sce_modeling_results$enrichment))]
colnames(registration_t_stats) <- gsub("^t_stat_", "", colnames(registration_t_stats))

## cell types x gene
dim(registration_t_stats)
#> [1] 18620    15

## check out table
registration_t_stats[1:5, 1:5]
#>                       Astro    Excit_A    Excit_B    Excit_C    Excit_D
#> ENSG00000238009 -0.71009456  0.7957792  0.0497619  0.6825793  0.5526941
#> ENSG00000237491 -4.24672326  1.7724150  1.6819367  0.9047336  2.8222782
#> ENSG00000225880  0.06152726  0.6941825  0.9819037 -0.1958094 -0.7766439
#> ENSG00000223764  7.69037575 -0.4106009 -0.4107015 -0.4106542 -0.4107933
#> ENSG00000187634 10.14422194 -0.4721603  0.2733466 -0.1397438 -1.0111055

2.7 Correlate statsics with Layer Reference

cor_layer <- layer_stat_cor(
    stats = registration_t_stats,
    modeling_results = layer_modeling_results,
    model_type = "enrichment",
    top_n = 100
)

cor_layer
#>                 WM       Layer6      Layer5       Layer4      Layer3       Layer2      Layer1
#> Oligo    0.7536847 -0.038947167 -0.22313462 -0.216143460 -0.39418956 -0.336272668 -0.04559186
#> Astro    0.2902852 -0.215372148 -0.32009814 -0.320189919 -0.24920587 -0.128017818  0.66950829
#> OPC      0.3309766 -0.076892980 -0.19246279 -0.254222683 -0.21417502 -0.073209686  0.22845154
#> Micro    0.2548264 -0.066391136 -0.14907964 -0.132236668 -0.18353524 -0.118307525  0.19395081
#> Mural    0.1652788 -0.046487411 -0.12107471 -0.190269592 -0.11924637 -0.076057486  0.24642004
#> Excit_B -0.3718199 -0.119246604 -0.26805619 -0.034391012  0.57980667  0.675251157 -0.12563319
#> Excit_C -0.5232250 -0.187359062  0.01974092  0.381935813  0.69540412  0.326841133 -0.26729177
#> Excit_A -0.3902148  0.128043074  0.63735797  0.405168995  0.06011232 -0.171933501 -0.38823893
#> Excit_E -0.3651414  0.489930498  0.13044242 -0.001907714  0.15141352  0.194625039 -0.33496664
#> Excit_F -0.2108888  0.443748447  0.35812706 -0.064476608 -0.09199288 -0.012550826 -0.30917103
#> Inhib_A -0.2497876 -0.035402374  0.08242644  0.105770587  0.12358824  0.133982729  0.01021479
#> Inhib_C -0.2100641 -0.080878595  0.05282252  0.052311336  0.09462094  0.160416496  0.06801709
#> Excit_D -0.3363905 -0.004936028  0.29430601  0.402875525  0.19153576 -0.014930006 -0.27248937
#> Inhib_B -0.1574240 -0.016415116  0.15172576  0.190096545  0.05317902 -0.004791031 -0.09307309
#> Inhib_D -0.2224211 -0.062788566  0.20260766  0.314434866  0.08963108 -0.067981519 -0.09229913

3 Explore Results

Now we can use these correlation values to learn about the cell types.

3.1 Create Heatmap of Correlations

We can see from this heatmap what layers the different cell types are associated with.

  • Oligo with WM

  • Astro with Layer 1

  • Excitatory neurons to different layers of the cortex

  • Weak associate with Inhibitory Neurons

layer_stat_cor_plot(cor_layer, max = max(cor_layer))

3.2 Annotate Cell Types by Top Correlation

We can use annotate_registered_clusters to create annotation labels for the cell types based on the correlation values.

anno <- annotate_registered_clusters(
    cor_stats_layer = cor_layer,
    confidence_threshold = 0.25,
    cutoff_merge_ratio = 0.25
)

anno
#>    cluster layer_confidence layer_label
#> 1    Oligo             good          WM
#> 2    Astro             good          L1
#> 3      OPC             good          WM
#> 4    Micro             good          WM
#> 5    Mural             poor         L1*
#> 6  Excit_B             good        L2/3
#> 7  Excit_C             good          L3
#> 8  Excit_A             good          L5
#> 9  Excit_E             good          L6
#> 10 Excit_F             good        L6/5
#> 11 Inhib_A             poor     L2/3/4*
#> 12 Inhib_C             poor         L2*
#> 13 Excit_D             good          L4
#> 14 Inhib_B             poor         L4*
#> 15 Inhib_D             good          L4

4 Reproducibility

The spatialLIBD package (Pardo, Spangler, Weber et al., 2022) was made possible thanks to:

This package was developed using biocthis.

Code for creating the vignette

## Create the vignette
library("rmarkdown")
system.time(render("guide_to_spatial_registration.Rmd", "BiocStyle::html_document"))

## Extract the R code
library("knitr")
knit("guide_to_spatial_registration.Rmd", tangle = TRUE)

Date the vignette was generated.

#> [1] "2023-11-30 10:58:13 EST"

Wallclock time spent generating the vignette.

#> Time difference of 47.071 secs

R session information.

#> ─ Session info ───────────────────────────────────────────────────────────────────────────────────────────────────────
#>  setting  value
#>  version  R version 4.3.2 Patched (2023-11-13 r85521)
#>  os       Ubuntu 22.04.3 LTS
#>  system   x86_64, linux-gnu
#>  ui       X11
#>  language (EN)
#>  collate  C
#>  ctype    en_US.UTF-8
#>  tz       America/New_York
#>  date     2023-11-30
#>  pandoc   2.7.3 @ /usr/bin/ (via rmarkdown)
#> 
#> ─ Packages ───────────────────────────────────────────────────────────────────────────────────────────────────────────
#>  package                * version   date (UTC) lib source
#>  abind                    1.4-5     2016-07-21 [2] CRAN (R 4.3.2)
#>  AnnotationDbi            1.64.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  AnnotationHub            3.10.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  attempt                  0.3.1     2020-05-03 [2] CRAN (R 4.3.2)
#>  backports                1.4.1     2021-12-13 [2] CRAN (R 4.3.2)
#>  beachmat                 2.18.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  beeswarm                 0.4.0     2021-06-01 [2] CRAN (R 4.3.2)
#>  benchmarkme              1.0.8     2022-06-12 [2] CRAN (R 4.3.2)
#>  benchmarkmeData          1.0.4     2020-04-23 [2] CRAN (R 4.3.2)
#>  bibtex                   0.5.1     2023-01-26 [2] CRAN (R 4.3.2)
#>  Biobase                * 2.62.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocFileCache          * 2.10.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocGenerics           * 0.48.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocIO                   1.12.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocManager              1.30.22   2023-08-08 [2] CRAN (R 4.3.2)
#>  BiocNeighbors            1.20.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocParallel             1.36.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocSingular             1.18.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocStyle              * 2.30.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  BiocVersion              3.18.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  Biostrings               2.70.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  bit                      4.0.5     2022-11-15 [2] CRAN (R 4.3.2)
#>  bit64                    4.0.5     2020-08-30 [2] CRAN (R 4.3.2)
#>  bitops                   1.0-7     2021-04-24 [2] CRAN (R 4.3.2)
#>  blob                     1.2.4     2023-03-17 [2] CRAN (R 4.3.2)
#>  bookdown                 0.36      2023-10-16 [2] CRAN (R 4.3.2)
#>  bslib                    0.6.1     2023-11-28 [2] CRAN (R 4.3.2)
#>  cachem                   1.0.8     2023-05-01 [2] CRAN (R 4.3.2)
#>  cli                      3.6.1     2023-03-23 [2] CRAN (R 4.3.2)
#>  codetools                0.2-19    2023-02-01 [3] CRAN (R 4.3.2)
#>  colorspace               2.1-0     2023-01-23 [2] CRAN (R 4.3.2)
#>  config                   0.3.2     2023-08-30 [2] CRAN (R 4.3.2)
#>  cowplot                  1.1.1     2020-12-30 [2] CRAN (R 4.3.2)
#>  crayon                   1.5.2     2022-09-29 [2] CRAN (R 4.3.2)
#>  curl                     5.1.0     2023-10-02 [2] CRAN (R 4.3.2)
#>  data.table               1.14.8    2023-02-17 [2] CRAN (R 4.3.2)
#>  DBI                      1.1.3     2022-06-18 [2] CRAN (R 4.3.2)
#>  dbplyr                 * 2.4.0     2023-10-26 [2] CRAN (R 4.3.2)
#>  DelayedArray             0.28.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  DelayedMatrixStats       1.24.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  digest                   0.6.33    2023-07-07 [2] CRAN (R 4.3.2)
#>  doParallel               1.0.17    2022-02-07 [2] CRAN (R 4.3.2)
#>  dotCall64                1.1-1     2023-11-28 [2] CRAN (R 4.3.2)
#>  dplyr                    1.1.4     2023-11-17 [2] CRAN (R 4.3.2)
#>  dqrng                    0.3.1     2023-08-30 [2] CRAN (R 4.3.2)
#>  DropletUtils             1.22.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  DT                       0.30      2023-10-05 [2] CRAN (R 4.3.2)
#>  edgeR                    4.0.2     2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  ellipsis                 0.3.2     2021-04-29 [2] CRAN (R 4.3.2)
#>  evaluate                 0.23      2023-11-01 [2] CRAN (R 4.3.2)
#>  ExperimentHub            2.10.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  fansi                    1.0.5     2023-10-08 [2] CRAN (R 4.3.2)
#>  farver                   2.1.1     2022-07-06 [2] CRAN (R 4.3.2)
#>  fastmap                  1.1.1     2023-02-24 [2] CRAN (R 4.3.2)
#>  fields                   15.2      2023-08-17 [2] CRAN (R 4.3.2)
#>  filelock                 1.0.2     2018-10-05 [2] CRAN (R 4.3.2)
#>  foreach                  1.5.2     2022-02-02 [2] CRAN (R 4.3.2)
#>  generics                 0.1.3     2022-07-05 [2] CRAN (R 4.3.2)
#>  GenomeInfoDb           * 1.38.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  GenomeInfoDbData         1.2.11    2023-11-14 [2] Bioconductor
#>  GenomicAlignments        1.38.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  GenomicRanges          * 1.54.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  ggbeeswarm               0.7.2     2023-04-29 [2] CRAN (R 4.3.2)
#>  ggplot2                  3.4.4     2023-10-12 [2] CRAN (R 4.3.2)
#>  ggrepel                  0.9.4     2023-10-13 [2] CRAN (R 4.3.2)
#>  glue                     1.6.2     2022-02-24 [2] CRAN (R 4.3.2)
#>  golem                    0.4.1     2023-06-05 [2] CRAN (R 4.3.2)
#>  gridExtra                2.3       2017-09-09 [2] CRAN (R 4.3.2)
#>  gtable                   0.3.4     2023-08-21 [2] CRAN (R 4.3.2)
#>  HDF5Array                1.30.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  highr                    0.10      2022-12-22 [2] CRAN (R 4.3.2)
#>  htmltools                0.5.7     2023-11-03 [2] CRAN (R 4.3.2)
#>  htmlwidgets              1.6.3     2023-11-22 [2] CRAN (R 4.3.2)
#>  httpuv                   1.6.12    2023-10-23 [2] CRAN (R 4.3.2)
#>  httr                     1.4.7     2023-08-15 [2] CRAN (R 4.3.2)
#>  interactiveDisplayBase   1.40.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  IRanges                * 2.36.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  irlba                    2.3.5.1   2022-10-03 [2] CRAN (R 4.3.2)
#>  iterators                1.0.14    2022-02-05 [2] CRAN (R 4.3.2)
#>  jquerylib                0.1.4     2021-04-26 [2] CRAN (R 4.3.2)
#>  jsonlite                 1.8.7     2023-06-29 [2] CRAN (R 4.3.2)
#>  KEGGREST                 1.42.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  knitr                    1.45      2023-10-30 [2] CRAN (R 4.3.2)
#>  labeling                 0.4.3     2023-08-29 [2] CRAN (R 4.3.2)
#>  later                    1.3.1     2023-05-02 [2] CRAN (R 4.3.2)
#>  lattice                  0.22-5    2023-10-24 [3] CRAN (R 4.3.2)
#>  lazyeval                 0.2.2     2019-03-15 [2] CRAN (R 4.3.2)
#>  lifecycle                1.0.4     2023-11-07 [2] CRAN (R 4.3.2)
#>  limma                    3.58.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  lobstr                 * 1.1.2     2022-06-22 [2] CRAN (R 4.3.2)
#>  locfit                   1.5-9.8   2023-06-11 [2] CRAN (R 4.3.2)
#>  lubridate                1.9.3     2023-09-27 [2] CRAN (R 4.3.2)
#>  magick                   2.8.1     2023-10-22 [2] CRAN (R 4.3.2)
#>  magrittr                 2.0.3     2022-03-30 [2] CRAN (R 4.3.2)
#>  maps                     3.4.1.1   2023-11-03 [2] CRAN (R 4.3.2)
#>  Matrix                   1.6-3     2023-11-14 [3] CRAN (R 4.3.2)
#>  MatrixGenerics         * 1.14.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  matrixStats            * 1.1.0     2023-11-07 [2] CRAN (R 4.3.2)
#>  memoise                  2.0.1     2021-11-26 [2] CRAN (R 4.3.2)
#>  mime                     0.12      2021-09-28 [2] CRAN (R 4.3.2)
#>  munsell                  0.5.0     2018-06-12 [2] CRAN (R 4.3.2)
#>  paletteer                1.5.0     2022-10-19 [2] CRAN (R 4.3.2)
#>  pillar                   1.9.0     2023-03-22 [2] CRAN (R 4.3.2)
#>  pkgconfig                2.0.3     2019-09-22 [2] CRAN (R 4.3.2)
#>  plotly                   4.10.3    2023-10-21 [2] CRAN (R 4.3.2)
#>  plyr                     1.8.9     2023-10-02 [2] CRAN (R 4.3.2)
#>  png                      0.1-8     2022-11-29 [2] CRAN (R 4.3.2)
#>  prettyunits              1.2.0     2023-09-24 [2] CRAN (R 4.3.2)
#>  promises                 1.2.1     2023-08-10 [2] CRAN (R 4.3.2)
#>  purrr                    1.0.2     2023-08-10 [2] CRAN (R 4.3.2)
#>  R.methodsS3              1.8.2     2022-06-13 [2] CRAN (R 4.3.2)
#>  R.oo                     1.25.0    2022-06-12 [2] CRAN (R 4.3.2)
#>  R.utils                  2.12.3    2023-11-18 [2] CRAN (R 4.3.2)
#>  R6                       2.5.1     2021-08-19 [2] CRAN (R 4.3.2)
#>  rappdirs                 0.3.3     2021-01-31 [2] CRAN (R 4.3.2)
#>  RColorBrewer             1.1-3     2022-04-03 [2] CRAN (R 4.3.2)
#>  Rcpp                     1.0.11    2023-07-06 [2] CRAN (R 4.3.2)
#>  RCurl                    1.98-1.13 2023-11-02 [2] CRAN (R 4.3.2)
#>  RefManageR             * 1.4.0     2022-09-30 [2] CRAN (R 4.3.2)
#>  rematch2                 2.1.2     2020-05-01 [2] CRAN (R 4.3.2)
#>  restfulr                 0.0.15    2022-06-16 [2] CRAN (R 4.3.2)
#>  rhdf5                    2.46.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  rhdf5filters             1.14.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  Rhdf5lib                 1.24.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  rjson                    0.2.21    2022-01-09 [2] CRAN (R 4.3.2)
#>  rlang                    1.1.2     2023-11-04 [2] CRAN (R 4.3.2)
#>  rmarkdown                2.25      2023-09-18 [2] CRAN (R 4.3.2)
#>  Rsamtools                2.18.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  RSQLite                  2.3.3     2023-11-04 [2] CRAN (R 4.3.2)
#>  rsvd                     1.0.5     2021-04-16 [2] CRAN (R 4.3.2)
#>  rtracklayer            * 1.62.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  S4Arrays                 1.2.0     2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  S4Vectors              * 0.40.2    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  sass                     0.4.7     2023-07-15 [2] CRAN (R 4.3.2)
#>  ScaledMatrix             1.10.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  scales                   1.3.0     2023-11-28 [2] CRAN (R 4.3.2)
#>  scater                   1.30.1    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  scuttle                  1.12.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  sessioninfo            * 1.2.2     2021-12-06 [2] CRAN (R 4.3.2)
#>  shiny                    1.8.0     2023-11-17 [2] CRAN (R 4.3.2)
#>  shinyWidgets             0.8.0     2023-08-30 [2] CRAN (R 4.3.2)
#>  SingleCellExperiment   * 1.24.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  spam                     2.10-0    2023-10-23 [2] CRAN (R 4.3.2)
#>  SparseArray              1.2.2     2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  sparseMatrixStats        1.14.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  SpatialExperiment      * 1.12.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  spatialLIBD            * 1.14.1    2023-11-30 [1] Bioconductor 3.18 (R 4.3.2)
#>  statmod                  1.5.0     2023-01-06 [2] CRAN (R 4.3.2)
#>  stringi                  1.8.2     2023-11-23 [2] CRAN (R 4.3.2)
#>  stringr                  1.5.1     2023-11-14 [2] CRAN (R 4.3.2)
#>  SummarizedExperiment   * 1.32.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  tibble                   3.2.1     2023-03-20 [2] CRAN (R 4.3.2)
#>  tidyr                    1.3.0     2023-01-24 [2] CRAN (R 4.3.2)
#>  tidyselect               1.2.0     2022-10-10 [2] CRAN (R 4.3.2)
#>  timechange               0.2.0     2023-01-11 [2] CRAN (R 4.3.2)
#>  utf8                     1.2.4     2023-10-22 [2] CRAN (R 4.3.2)
#>  vctrs                    0.6.4     2023-10-12 [2] CRAN (R 4.3.2)
#>  vipor                    0.4.5     2017-03-22 [2] CRAN (R 4.3.2)
#>  viridis                  0.6.4     2023-07-22 [2] CRAN (R 4.3.2)
#>  viridisLite              0.4.2     2023-05-02 [2] CRAN (R 4.3.2)
#>  withr                    2.5.2     2023-10-30 [2] CRAN (R 4.3.2)
#>  xfun                     0.41      2023-11-01 [2] CRAN (R 4.3.2)
#>  XML                      3.99-0.16 2023-11-29 [2] CRAN (R 4.3.2)
#>  xml2                     1.3.5     2023-07-06 [2] CRAN (R 4.3.2)
#>  xtable                   1.8-4     2019-04-21 [2] CRAN (R 4.3.2)
#>  XVector                  0.42.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#>  yaml                     2.3.7     2023-01-23 [2] CRAN (R 4.3.2)
#>  zlibbioc                 1.48.0    2023-11-29 [2] Bioconductor 3.18 (R 4.3.2)
#> 
#>  [1] /tmp/RtmpsA2RYx/Rinst76f1b48d43843
#>  [2] /home/biocbuild/bbs-3.18-bioc/R/site-library
#>  [3] /home/biocbuild/bbs-3.18-bioc/R/library
#> 
#> ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────

5 Bibliography

This vignette was generated using BiocStyle (Oleś, 2023) with knitr (Xie, 2023) and rmarkdown (Allaire, Xie, Dervieux et al., 2023) running behind the scenes.

Citations made with RefManageR (McLean, 2017).

[1] J. Allaire, Y. Xie, C. Dervieux, et al. rmarkdown: Dynamic Documents for R. R package version 2.25. 2023. URL: https://github.com/rstudio/rmarkdown.

[2] K. R. Maynard, L. Collado-Torres, L. M. Weber, et al. “Transcriptome-scale spatial gene expression in the human dorsolateral prefrontal cortex”. In: Nature Neuroscience (2021). DOI: 10.1038/s41593-020-00787-0. URL: https://www.nature.com/articles/s41593-020-00787-0.

[3] M. W. McLean. “RefManageR: Import and Manage BibTeX and BibLaTeX References in R”. In: The Journal of Open Source Software (2017). DOI: 10.21105/joss.00338.

[4] A. Oleś. BiocStyle: Standard styles for vignettes and other Bioconductor documents. R package version 2.30.0. 2023. URL: https://github.com/Bioconductor/BiocStyle.

[5] B. Pardo, A. Spangler, L. M. Weber, et al. “spatialLIBD: an R/Bioconductor package to visualize spatially-resolved transcriptomics data”. In: BMC Genomics (2022). DOI: 10.1186/s12864-022-08601-w. URL: https://doi.org/10.1186/s12864-022-08601-w.

[6] R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna, Austria, 2023. URL: https://www.R-project.org/.

[7] M. N. Tran, K. R. Maynard, A. Spangler, et al. “Single-nucleus transcriptome analysis reveals cell-type-specific molecular signatures across reward circuitry in the human brain”. In: Neuron (2021). DOI: 10.1016/j.neuron.2021.09.001.

[8] H. Wickham. “testthat: Get Started with Testing”. In: The R Journal 3 (2011), pp. 5–10. URL: https://journal.r-project.org/archive/2011-1/RJournal_2011-1_Wickham.pdf.

[9] H. Wickham, W. Chang, R. Flight, et al. sessioninfo: R Session Information. R package version 1.2.2. 2021. URL: https://CRAN.R-project.org/package=sessioninfo.

[10] Y. Xie. knitr: A General-Purpose Package for Dynamic Report Generation in R. R package version 1.45. 2023. URL: https://yihui.org/knitr/.