Introduction

In the single cell World, which includes flow cytometry, mass cytometry, single-cell RNA-seq (scRNA-seq), and others, there is a need to improve data visualisation and to bring analysis capabilities to researchers even from non-technical backgrounds. scDataviz (Blighe 2020) attempts to fit into this space, while also catering for advanced users. Additonally, due to the way that scDataviz is designed, which is based on SingleCellExperiment (Lun and Risso 2020), it has a ‘plug and play’ feel, and immediately lends itself as flexibile and compatibile with studies that go beyond scDataviz. Finally, the graphics in scDataviz are generated via the ggplot (Wickham 2016) engine, which means that users can ‘add on’ features to these with ease.

This package just provides some additional functions for dataviz and clustering, and provides another way of identifying cell-types in clusters. It is not strictly intended as a standalone analysis package. For a comprehensive high-dimensional cytometry workflow, it is recommended to check out the work by Nowicka et al. CyTOF workflow: differential discovery in high-throughput high-dimensional cytometry datasets. For a more comprehensive scRNA-seq workflow, please check out OSCA and Analysis of single cell RNA-seq data.

Installation

  if (!requireNamespace('BiocManager', quietly = TRUE))
    install.packages('BiocManager')

  BiocManager::install('scDataviz')

  devtools::install_github('kevinblighe/scDataviz')

  library(scDataviz)

Tutorial 1: CyTOF FCS data

Here, we will utilise some of the flow cytometry data from Deep phenotyping detects a pathological CD4+ T-cell complosome signature in systemic sclerosis.

This can normally be downloadedd via git clone from your command prompt:

  git clone https://github.com/kevinblighe/scDataviz_data/ ;

In a practical situation, we would normally read in this data from the raw FCS files and then QC filter, normalise, and transform them. This can be achieved via the processFCS function, which, by default, also removes variables based on low variance and downsamples [randomly] your data to 100000 variables. The user can change these via the downsample and downsampleVar parameters. An example (not run) is given below:

  filelist <- list.files(
    path = "scDataviz_data/FCS/",
    pattern = "*.fcs|*.FCS",
    full.names = TRUE)
  filelist

  metadata <- data.frame(
    sample = gsub('\\ [A-Za-z0-9]*\\.fcs$', '',
      gsub('scDataviz_data\\/FCS\\/\\/', '', filelist)),
    group = c(rep('Healthy', 7), rep('Disease', 11)),
    treatment = gsub('\\.fcs$', '',
      gsub('scDataviz_data\\/FCS\\/\\/[A-Z0-9]*\\ ', '', filelist)),
    row.names = filelist,
    stringsAsFactors = FALSE)
  metadata

  inclusions <- c('Yb171Di','Nd144Di','Nd145Di',
    'Er168Di','Tm169Di','Sm154Di','Yb173Di','Yb174Di',
    'Lu175Di','Nd143Di')

  markernames <- c('Foxp3','C3aR','CD4',
    'CD46','CD25','CD3','Granzyme B','CD55',
    'CD279','CD45RA')

  names(markernames) <- inclusions
  markernames

  exclusions <- c('Time','Event_length','BCKG190Di',
    'Center','Offset','Width','Residual')

  sce <- processFCS(
    files = filelist,
    metadata = metadata,
    transformation = TRUE,
    transFun = function (x) asinh(x),
    asinhFactor = 5,
    downsample = 10000,
    downsampleVar = 0.7,
    colsRetain = inclusions,
    colsDiscard = exclusions,
    newColnames = markernames)

In flow and mass cytometry, getting the correct marker names in the FCS files can be surprisingly difficult. In many cases, from experience, a facility may label the markers by their metals, such as Iridium (Ir), Ruthenium (Ru), Terbium (Tb), et cetera - this is the case for the data used in this tutorial. The true marker names may be held as pData encoded within each FCS, accessible via:

  library(flowCore)
  pData(parameters(
    read.FCS(filelist[[4]], transformation = FALSE, emptyValue = FALSE)))

Whatever the case, it is important to sort out marker naming issues prior to the experiment being conducted in order to avoid any confusion.

For this vignette, due to the fact that the raw FCS data is > 500 megabytes, we will work with a smaller pre-prepared dataset that has been downsampled to 10000 cells using the above code. This data comes included with the package.

Load the pre-prepared complosome data.

  load(system.file('extdata/', 'complosome.rdata', package = 'scDataviz'))

One can also create a new SingleCellExperiment object manually using any type of data, including any data from scRNA-seq produced elsewhere. Import functions for data deriving from other sources is covered in Tutorials 2 and 3 in this vignette. All functions in scDataviz additionally accept data-frames or matrices on their own, de-necessitating the reliance on the SingleCellExperiment class.

Perform principal component analysis (PCA)

We can use the PCAtools (Blighe and Lun 2020) package for the purpose of performing PCA.

  library(PCAtools)
  p <- pca(assay(sce, 'scaled'), metadata = metadata(sce))

  biplot(p,
    x = 'PC1', y = 'PC2',
    lab = NULL,
    xlim = c(min(p$rotated[,'PC1'])-1, max(p$rotated[,'PC1'])+1),
    ylim = c(min(p$rotated[,'PC2'])-1, max(p$rotated[,'PC2'])+1),
    pointSize = 1.0,
    colby = 'treatment',
    legendPosition = 'right',
    title = 'PCA applied to CyTOF data',
    caption = paste0('10000 cells randomly selected after ',
      'having filtered for low variance'))

Perform PCA

We can add the rotated component loadings as a new reduced dimensional component to our dataset.

  reducedDim(sce, 'PCA') <- p$rotated

For more functionality via PCAtools, check the vignette: PCAtools: everything Principal Component Analysis

  sce <- performUMAP(sce)

  config <- umap::umap.defaults
  config$min_dist <- 0.5
  performUMAP(sce, config = config)

  elbow <- findElbowPoint(p$variance)
  horn <- parallelPCA(assay(sce, 'scaled'))

  elbow

## PC3 
##   3

  horn$n

## [1] 1

  sce <- performUMAP(sce, reducedDim = 'PCA', dims = c(1:5))

Create a contour plot of the UMAP layout

This and the remaining sections in this tutorial are about producing great visualisations of the data and attempting to make sense of it, while not fully overlapping with functionalioty provided by other programs that operate in tis space.

With the contour plot, we are essentially looking at celluar density. It can provide for a beautiful viusualisation in a manuscript while also serving as a useful QC tool: if the density is ‘scrunched up’ into a single area in the plot space, then there are likely issues with your input data distribution. We want to see well-separated, high density ‘islands’, or, at least, gradual gradients that blend into one another across high density ‘peaks’.

  ggout1 <- contourPlot(sce,
    reducedDim = 'UMAP',
    bins = 150,
    subtitle = 'UMAP performed on expression values',
    legendLabSize = 18,
    axisLabSize = 22,
    titleLabSize = 22,
    subtitleLabSize = 18,
    captionLabSize = 18)

  ggout2 <- contourPlot(sce,
    reducedDim = 'UMAP_PCA',
    bins = 150,
    subtitle = 'UMAP performed on PC eigenvectors',
    legendLabSize = 18,
    axisLabSize = 22,
    titleLabSize = 22,
    subtitleLabSize = 18,
    captionLabSize = 18)

  cowplot::plot_grid(ggout1, ggout2,
    labels = c('A','B'),
    ncol = 2, align = "l", label_size = 24)

Create a contour plot of the UMAP layout

Show marker expression across the layout

Here, we randomly select some markers and then plot their expression profiles across the UMAP layouts.

  markers <- sample(rownames(sce), 6)
  markers

## [1] "CD55"   "CD25"   "Foxp3"  "C3aR"   "CD46"   "CD45RA"

  ggout1 <- markerExpression(sce,
    markers = markers,
    subtitle = 'UMAP performed on expression values',
    nrow = 1, ncol = 6,
    legendKeyHeight = 1.0,
    legendLabSize = 18,
    stripLabSize = 22,
    axisLabSize = 22,
    titleLabSize = 22,
    subtitleLabSize = 18,
    captionLabSize = 18)

  ggout2 <- markerExpression(sce,
    markers = markers,
    reducedDim = 'UMAP_PCA',
    subtitle = 'UMAP performed on PC eigenvectors',
    nrow = 1, ncol = 6,
    col = c('white', 'darkblue'),
    legendKeyHeight = 1.0,
    legendLabSize = 18,
    stripLabSize = 22,
    axisLabSize = 22,
    titleLabSize = 22,
    subtitleLabSize = 18,
    captionLabSize = 18)

  cowplot::plot_grid(ggout1, ggout2,
    labels = c('A','B'),
    nrow = 2, align = "l", label_size = 24)

Show marker expression across the layout

Shade cells by metadata

Shading cells by metadata can be useful for identifying any batch effects, but also useful for visualising, e.g., differences across treatments.

First, let’s take a look inside the metadata that we have.

  head(metadata(sce))

##       sample   group treatment
## cell1    P00 Disease    Unstim
## cell2    P00 Disease    Unstim
## cell3    P04 Disease      CD46
## cell4    P03 Disease      CD46
## cell5    P08 Disease    Unstim
## cell6    P00 Disease      CD46

  levels(metadata(sce)$group)

## [1] "Healthy" "Disease"

  levels(metadata(sce)$treatment)

## [1] "CD46"   "Unstim" "CD3"

  ggout1 <- metadataPlot(sce,
    colby = 'group',
    colkey = c(Healthy = 'royalblue', Disease = 'red2'),
    title = 'Disease status',
    subtitle = 'UMAP performed on expression values',
    legendLabSize = 16,
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 16,
    captionLabSize = 16)

  ggout2 <- metadataPlot(sce,
    reducedDim = 'UMAP_PCA',
    colby = 'group',
    colkey = c(Healthy = 'royalblue', Disease = 'red2'),
    title = 'Disease status',
    subtitle = 'UMAP performed on PC eigenvectors',
    legendLabSize = 16,
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 16,
    captionLabSize = 16)

  ggout3 <- metadataPlot(sce,
    colby = 'treatment',
    title = 'Treatment type',
    subtitle = 'UMAP performed on expression values',
    legendLabSize = 16,
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 16,
    captionLabSize = 16)

  ggout4 <- metadataPlot(sce,
    reducedDim = 'UMAP_PCA',
    colby = 'treatment',
    title = 'Treatment type',
    subtitle = 'UMAP performed on PC eigenvectors',
    legendLabSize = 16,
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 16,
    captionLabSize = 16)

  cowplot::plot_grid(ggout1, ggout3, ggout2, ggout4,
    labels = c('A','B','C','D'),
    nrow = 2, ncol = 2, align = "l", label_size = 24)

Shade cells by metadata

Find ideal clusters in the UMAP layout via k-nearest neighbours

This function utilises the k nearest neighbours (k-NN) approach from Seurat, which works quite well on flow cytometry and CyTOF UMAP layouts, from my experience.

  sce <- clusKNN(sce,
    k.param = 20,
    prune.SNN = 1/15,
    resolution = 0.01,
    algorithm = 2,
    verbose = FALSE)

  sce <- clusKNN(sce,
    reducedDim = 'UMAP_PCA',
    clusterAssignName = 'Cluster_PCA',
    k.param = 20,
    prune.SNN = 1/15,
    resolution = 0.01,
    algorithm = 2,
    verbose = FALSE)

  ggout1 <- plotClusters(sce,
    clusterColname = 'Cluster',
    labSize = 7.0,
    subtitle = 'UMAP performed on expression values',
    caption = paste0('Note: clusters / communities identified via',
      '\nLouvain algorithm with multilevel refinement'),
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 16,
    captionLabSize = 16)

  ggout2 <- plotClusters(sce,
    clusterColname = 'Cluster_PCA',
    reducedDim = 'UMAP_PCA',
    labSize = 7.0,
    subtitle = 'UMAP performed on PC eigenvectors',
    caption = paste0('Note: clusters / communities identified via',
      '\nLouvain algorithm with multilevel refinement'),
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 16,
    captionLabSize = 16)

  cowplot::plot_grid(ggout1, ggout2,
    labels = c('A','B'),
    ncol = 2, align = "l", label_size = 24)

Find ideal clusters in the UMAP layout via k-nearest neighbours

Plot marker expression per identified cluster

  markerExpressionPerCluster(sce,
    caption = 'Cluster assignments based on UMAP performed on expression values',
    stripLabSize = 22,
    axisLabSize = 22,
    titleLabSize = 22,
    subtitleLabSize = 18,
    captionLabSize = 18)

  clusters <- unique(metadata(sce)[['Cluster_PCA']])
  clusters

## [1] 3 0 1 6 4 2 5 8 7

  markers <- sample(rownames(sce), 5)
  markers

## [1] "CD25"       "Granzyme B" "CD46"       "CD45RA"     "CD279"

  markerExpressionPerCluster(sce,
    clusters = clusters,
    clusterAssign = metadata(sce)[['Cluster_PCA']],
    markers = markers,
    nrow = 2, ncol = 5,
    caption = 'Cluster assignments based on UMAP performed on PC eigenvectors',
    stripLabSize = 22,
    axisLabSize = 22,
    titleLabSize = 22,
    subtitleLabSize = 18,
    captionLabSize = 18)

Plot marker expression per identified cluster2

Try all markers across a single cluster:

  cluster <- sample(unique(metadata(sce)[['Cluster']]), 1)
  cluster

## [1] 7

  markerExpressionPerCluster(sce,
    clusters = cluster,
    markers = rownames(sce),
    stripLabSize = 20,
    axisLabSize = 20,
    titleLabSize = 20,
    subtitleLabSize = 14,
    captionLabSize = 12)

Plot marker expression per identified cluster3

Determine enriched markers in each cluster and plot the expression signature

This method also calculates metacluster abundances across a chosen phenotype. The function returns a data-frame, which can then be exported to do other analyses.

Disease vs Healthy metacluster abundances

  markerEnrichment(sce,
    method = 'quantile',
    studyvarID = 'group')

Cluster	nCells	TotalCells	PercentCells	NegMarkers	PosMarkers	PerCent_HD00	PerCent_HD01	PerCent_HD262	PerCent_P00	PerCent_P02	PerCent_P03	PerCent_P04	PerCent_P08	nCell_Healthy	nCell_Disease
0	3308	10000	33.08	NA	CD25+	0.0604595	9.9758162	26.299879	23.5489722	6.4087062	8.2527207	18.1680774	7.2853688	1202	2106
1	1927	10000	19.27	CD25-CD279-	CD3+CD45RA+	0.0000000	0.1037883	61.961598	0.6746238	0.2075765	0.0518941	0.0518941	36.9486248	1196	731
2	1513	10000	15.13	NA	CD25+	5.2875083	0.9914078	14.144085	2.7759418	19.7620621	5.7501652	46.4639788	4.8248513	309	1204
3	1303	10000	13.03	Granzyme B-CD279-	CD3+	15.5794321	7.9815810	2.455871	63.7759018	3.8372985	3.6838066	1.9186493	0.7674597	339	964
4	960	10000	9.60	NA	CD3+	16.5625000	26.0416667	7.604167	28.0208333	3.2291667	10.0000000	6.0416667	2.5000000	482	478
5	560	10000	5.60	CD25-	NA	0.3571429	15.0000000	36.785714	9.8214286	0.1785714	8.9285714	0.1785714	28.7500000	292	268
6	320	10000	3.20	CD46-CD279-	CD3+	0.0000000	39.6875000	1.875000	51.2500000	0.3125000	5.3125000	0.6250000	0.9375000	133	187
7	109	10000	1.09	CD46-	CD25+	0.9174312	0.0000000	0.000000	0.0000000	67.8899083	27.5229358	3.6697248	0.0000000	1	108

.

Treatment type metacluster abundances

  markerEnrichment(sce,
    sampleAbundances = FALSE,
    method = 'quantile',
    studyvarID = 'treatment')

Cluster	nCells	TotalCells	PercentCells	NegMarkers	PosMarkers	nCell_CD46	nCell_Unstim	nCell_CD3
0	3308	10000	33.08	NA	CD25+	3283	2	23
1	1927	10000	19.27	CD25-CD279-	CD3+CD45RA+	4	1919	4
2	1513	10000	15.13	NA	CD25+	20	409	1084
3	1303	10000	13.03	Granzyme B-CD279-	CD3+	5	1175	123
4	960	10000	9.60	NA	CD3+	4	770	186
5	560	10000	5.60	CD25-	NA	166	28	366
6	320	10000	3.20	CD46-CD279-	CD3+	308	10	2
7	109	10000	1.09	CD46-	CD25+	107	0	2

.

Expression signature

The expression signature is a quick way to visualise which markers are more or less expressed in each identified cluster of cells.

  plotSignatures(sce,
    labCex = 1.2,
    legendCex = 1.2,
    labDegree = 40)

Determine enriched markers in each cluster and plot the expression signature

Tutorial 2: Import from Seurat

Due to the fact that scDataviz is based on SingleCellExperiment, it has increased interoperability with other packages, including the popular Seurat (Stuart et al. 2018). Taking the data produced from the Seurat Tutorial on Peripheral Blood Mononuclear Cells (PBMCs), we can convert this to a SingleCellExperiment object recognisable by scDataviz via as.SingleCellExperiment().

When deriving from the Seurat route, be sure to manually assign the metadata slot, which is required for some functions. Also be sure to modify the default values for assay, reducedDim, and dimColnames, as these are assigned differently in Seurat.

  sce <- as.SingleCellExperiment(pbmc)

  metadata(sce) <- data.frame(colData(sce))

  markerExpression(sce,
    assay = 'logcounts',
    reducedDim = 'UMAP',
    dimColnames = c('UMAP_1','UMAP_2'),
    markers = c('CD79A', 'Cd79B', 'MS4A1'))

For markerEnrichment(), a typical command using an ex-Seurat object could be:

  markerEnrichment(sce,
    assay = 'logcounts',
    method = 'quantile',
    sampleAbundances = TRUE,
    sampleID = 'orig.ident',
    studyvarID = 'ident',
    clusterAssign = as.character(colData(sce)[['seurat_clusters']]))

Tutorial 3: Import any numerical data

scDataviz will work with any numerical data, too. Here, we show a quick example of how one can import a data-matrix of randomly-generated numbers that follow a negative binomial distribution, comprising 2500 cells and 20 markers:

  mat <- jitter(matrix(
    MASS::rnegbin(rexp(50000, rate=.1), theta = 4.5),
    ncol = 20))
  colnames(mat) <- paste0('CD', 1:ncol(mat))
  rownames(mat) <- paste0('cell', 1:nrow(mat))

  metadata <- data.frame(
    group = rep('A', nrow(mat)),
    row.names = rownames(mat),
    stringsAsFactors = FALSE)
  head(metadata)

##       group
## cell1     A
## cell2     A
## cell3     A
## cell4     A
## cell5     A
## cell6     A

  sce <- importData(mat,
    assayname = 'normcounts',
    metadata = metadata)
  sce

## class: SingleCellExperiment 
## dim: 20 2500 
## metadata(1): group
## assays(1): normcounts
## rownames(20): CD1 CD2 ... CD19 CD20
## rowData names(0):
## colnames(2500): cell1 cell2 ... cell2499 cell2500
## colData names(0):
## reducedDimNames(0):
## mainExpName: NULL
## altExpNames(0):

This will also work without any assigned metadata; however, having no metadata limits the functionality of the package.

  sce <- importData(mat,
    assayname = 'normcounts',
    metadata = NULL)
  sce

## class: SingleCellExperiment 
## dim: 20 2500 
## metadata(0):
## assays(1): normcounts
## rownames(20): CD1 CD2 ... CD19 CD20
## rowData names(0):
## colnames(2500): cell1 cell2 ... cell2499 cell2500
## colData names(0):
## reducedDimNames(0):
## mainExpName: NULL
## altExpNames(0):

Acknowledgments

• Jessica Timms
• James Opzoomer
• Shahram Kordasti
• Marcel Ramos (Bioconductor)
• Lori Shepherd (Bioconductor)
• Bioinformatics CRO
• Henrik Bengtsson

Session info

sessionInfo()

## 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] PCAtools_2.19.0             ggrepel_0.9.6              
##  [3] ggplot2_3.5.1               scDataviz_1.17.0           
##  [5] SingleCellExperiment_1.29.1 SummarizedExperiment_1.37.0
##  [7] Biobase_2.67.0              GenomicRanges_1.59.1       
##  [9] GenomeInfoDb_1.43.2         IRanges_2.41.1             
## [11] MatrixGenerics_1.19.0       matrixStats_1.4.1          
## [13] S4Vectors_0.45.2            BiocGenerics_0.53.3        
## [15] generics_0.1.3              kableExtra_1.4.0           
## [17] knitr_1.49                 
## 
## loaded via a namespace (and not attached):
##   [1] RcppAnnoy_0.0.22          splines_4.4.2            
##   [3] later_1.4.1               tibble_3.2.1             
##   [5] polyclip_1.10-7           fastDummies_1.7.4        
##   [7] lifecycle_1.0.4           globals_0.16.3           
##   [9] lattice_0.22-6            MASS_7.3-61              
##  [11] magrittr_2.0.3            plotly_4.10.4            
##  [13] sass_0.4.9                rmarkdown_2.29           
##  [15] jquerylib_0.1.4           yaml_2.3.10              
##  [17] httpuv_1.6.15             Seurat_5.1.0             
##  [19] sctransform_0.4.1         spam_2.11-0              
##  [21] askpass_1.2.1             sp_2.1-4                 
##  [23] spatstat.sparse_3.1-0     reticulate_1.40.0        
##  [25] cowplot_1.1.3             pbapply_1.7-2            
##  [27] buildtools_1.0.0          RColorBrewer_1.1-3       
##  [29] abind_1.4-8               zlibbioc_1.52.0          
##  [31] Rtsne_0.17                purrr_1.0.2              
##  [33] GenomeInfoDbData_1.2.13   irlba_2.3.5.1            
##  [35] listenv_0.9.1             spatstat.utils_3.1-1     
##  [37] maketools_1.3.1           umap_0.2.10.0            
##  [39] goftest_1.2-3             RSpectra_0.16-2          
##  [41] spatstat.random_3.3-2     dqrng_0.4.1              
##  [43] fitdistrplus_1.2-1        parallelly_1.39.0        
##  [45] svglite_2.1.3             DelayedMatrixStats_1.29.0
##  [47] leiden_0.4.3.1            codetools_0.2-20         
##  [49] DelayedArray_0.33.2       xml2_1.3.6               
##  [51] tidyselect_1.2.1          UCSC.utils_1.3.0         
##  [53] farver_2.1.2              ScaledMatrix_1.15.0      
##  [55] spatstat.explore_3.3-3    jsonlite_1.8.9           
##  [57] progressr_0.15.1          ggridges_0.5.6           
##  [59] survival_3.7-0            systemfonts_1.1.0        
##  [61] tools_4.4.2               ica_1.0-3                
##  [63] Rcpp_1.0.13-1             glue_1.8.0               
##  [65] gridExtra_2.3             SparseArray_1.7.2        
##  [67] xfun_0.49                 dplyr_1.1.4              
##  [69] withr_3.0.2               fastmap_1.2.0            
##  [71] fansi_1.0.6               openssl_2.2.2            
##  [73] digest_0.6.37             rsvd_1.0.5               
##  [75] R6_2.5.1                  mime_0.12                
##  [77] colorspace_2.1-1          scattermore_1.2          
##  [79] tensor_1.5                spatstat.data_3.1-4      
##  [81] utf8_1.2.4                tidyr_1.3.1              
##  [83] data.table_1.16.2         httr_1.4.7               
##  [85] htmlwidgets_1.6.4         S4Arrays_1.7.1           
##  [87] uwot_0.2.2                pkgconfig_2.0.3          
##  [89] gtable_0.3.6              RProtoBufLib_2.19.0      
##  [91] lmtest_0.9-40             XVector_0.47.0           
##  [93] sys_3.4.3                 htmltools_0.5.8.1        
##  [95] dotCall64_1.2             SeuratObject_5.0.2       
##  [97] scales_1.3.0              png_0.1-8                
##  [99] spatstat.univar_3.1-1     corrplot_0.95            
## [101] rstudioapi_0.17.1         reshape2_1.4.4           
## [103] nlme_3.1-166              cachem_1.1.0             
## [105] zoo_1.8-12                stringr_1.5.1            
## [107] KernSmooth_2.23-24        parallel_4.4.2           
## [109] miniUI_0.1.1.1            pillar_1.9.0             
## [111] grid_4.4.2                vctrs_0.6.5              
## [113] RANN_2.6.2                promises_1.3.2           
## [115] BiocSingular_1.23.0       cytolib_2.19.0           
## [117] beachmat_2.23.2           xtable_1.8-4             
## [119] cluster_2.1.6             evaluate_1.0.1           
## [121] isoband_0.2.7             cli_3.6.3                
## [123] compiler_4.4.2            rlang_1.1.4              
## [125] crayon_1.5.3              future.apply_1.11.3      
## [127] labeling_0.4.3            flowCore_2.19.0          
## [129] plyr_1.8.9                stringi_1.8.4            
## [131] viridisLite_0.4.2         deldir_2.0-4             
## [133] BiocParallel_1.41.0       munsell_0.5.1            
## [135] lazyeval_0.2.2            spatstat.geom_3.3-4      
## [137] Matrix_1.7-1              RcppHNSW_0.6.0           
## [139] patchwork_1.3.0           sparseMatrixStats_1.19.0 
## [141] future_1.34.0             shiny_1.9.1              
## [143] ROCR_1.0-11               igraph_2.1.1             
## [145] bslib_0.8.0

References

Blighe (2020)

Blighe and Lun (2020)

Lun and Risso (2020)

Stuart et al. (2018)

Wickham (2016)