# Identifying cellular neighborhood with SPIAT

Yuzhou Feng

#### 2022-12-20

#### Package

SPIAT 1.0.4

    library(SPIAT)

# 1 Cellular neighborhood

The aggregation of cells can result in ‘cellular neighbourhoods’. A neighbourhood is defined as a group of cells that cluster together. These can be homotypic, containing cells of a single class (e.g. immune cells), or heterotypic (e.g. a mixture of tumour and immune cells).

Function `identify_neighborhoods()` identifies cellular neighbourhoods. Users can select a subset of cell types of interest if desired. SPIAT includes three algorithms for the detection of neighbourhoods.

- Hierarchical Clustering algorithm: Euclidean distances between cells are calculated, and pairs of cells with a distance less than a specified radius are considered to be ‘interacting’, with the rest being ‘non-interacting’. Hierarchical clustering is then used to separate the clusters. Larger radii will result in the merging of individual clusters.
- dbscan
- phenograph

For Hierarchical Clustering algorithm and dbscan, users need to specify a radius that defines the distance for an interaction. We suggest users to test different radii and select the one that generates intuitive clusters upon visualisation. Cells not assigned to clusters are assigned as `Cluster_NA` in the output table. The argument `min_neighborhood_size` specifies the threshold of a neighborhood size to be considered as a neighborhood. Smaller neighbourhoods will be outputted, but will not be assigned a number.

Rphenograph uses the number of nearest neighbours to detect clusters. This number should be specified by `min_neighborhood_size` argument. We also encourage users to test different values.

For this part of the tutorial, we will use the image `image_no_markers` simulated with the `spaSim` package. This image contains “Tumour”, “Immune”, “Immune1” and “Immune2” cells without marker intensities.

    data("image_no_markers")

    plot_cell_categories(
      image_no_markers, c("Tumour", "Immune","Immune1","Immune2","Others"),
      c("red","blue","darkgreen", "brown","lightgray"), "Cell.Type")

Users are recommended to test out different radii and then visualise the clustering results. To aid in this process, users can use the `average_minimum_distance()` function, which calculates the average minimum distance between all cells in an image, and can be used as a starting point.

    average_minimum_distance(image_no_markers)

    ## [1] 17.01336

We then identify the cellular neighbourhoods using our hierarchical algorithm with a radius of 50, and with a minimum neighbourhood size of 100. Cells assigned to neighbourhoods smaller than 100 will be assigned to the “Cluster_NA” neighbourhood.

    clusters <- identify_neighborhoods(
      image_no_markers, method = "hierarchical", min_neighborhood_size = 100,
      cell_types_of_interest = c("Immune", "Immune1", "Immune2"), radius = 50,
      feature_colname = "Cell.Type")

This plot shows clusters of “Immune”, “Immune1” and “Immune2” cells. Each number and colour corresponds to a distinct cluster. Black cells correspond to ‘free’, un-clustered cells.

We can visualise the cell composition of neighborhoods. To do this, we can use `composition_of_neighborhoods()` to obtain the percentages of cells with a specific marker within each neighborhood and the number of cells in the neighborhood.

In this example we select cellular neighbourhoods with at least 5 cells.

    neighorhoods_vis <-
      composition_of_neighborhoods(clusters, feature_colname = "Cell.Type")
    neighorhoods_vis <-
      neighorhoods_vis[neighorhoods_vis$Total_number_of_cells >=5,]

Finally, we can use `plot_composition_heatmap()` to produce a heatmap showing the marker percentages within each cluster, which can be used to classify the derived neighbourhoods.

    plot_composition_heatmap(neighorhoods_vis, feature_colname="Cell.Type")

This plot shows that Cluster_1 and Cluster_2 contain all three types of immune cells. Cluster_3 does not have Immune1 cells. Cluster_1 and Cluster_2 are more similar to the free cells (cells not assigned to clusters) in their composition than Cluster_3.

# 2 Average Nearest Neighbour Index (ANNI)

We can test for the presence of neighbourhoods using ANNI. We can calculate the ANNI with the function `average_nearest_neighbor_index()`, which takes one cell type of interest (e.g. `Cluster_1` under `Neighborhood` column of `clusters` object) or a combinations of cell types (e.g. `Immune1` and `Immune2` cells under `Cell.Type` column of `image_no_markers` object) and outputs whether there is a clear neighbourhood (clustered) or unclear (dispersed/random), along with a P value for the estimate.

Here show the examples for both one cell type and multiple cell types.

    average_nearest_neighbor_index(clusters, reference_celltypes=c("Cluster_1"),
                                   feature_colname="Neighborhood", p_val = 0.05)

    ## $ANN_index
    ## [1] 0.4696608
    ##
    ## $pattern
    ## [1] "Clustered"
    ##
    ## $`p-value`
    ## [1] 2.623132e-68

    average_nearest_neighbor_index(
      image_no_markers, reference_celltypes=c("Immune", "Immune1" , "Immune2"),
      feature_colname="Cell.Type", p_val = 0.05)

    ## $ANN_index
    ## [1] 0.9968575
    ##
    ## $pattern
    ## [1] "Random"
    ##
    ## $`p-value`
    ## [1] 0.4000806

`p_val` is the cutoff to determine if a pattern is significant or not. If the p value of ANNI is larger than the threshold, the pattern will be “Random”. Although we give a default p value cutoff of 5e-6, we suggest the users to define their own cutoff based on the images and how they define the patterns “Clustered” and “Dispersed”.

# 3 You can access the vignettes for other modules of SPIAT here:

- Overview of SPIAT
- Data reading and formatting
- Quality control and visualisation
- Basic analysis
- Cell colocalisation
- Spatial heterogeneity
- Tissue structure

# 4 Reproducibility

    sessionInfo()

    ## R version 4.2.2 (2022-10-31 ucrt)
    ## Platform: x86_64-w64-mingw32/x64 (64-bit)
    ## Running under: Windows Server x64 (build 20348)
    ##
    ## Matrix products: default
    ##
    ## locale:
    ## [1] LC_COLLATE=C
    ## [2] LC_CTYPE=English_United States.utf8
    ## [3] LC_MONETARY=English_United States.utf8
    ## [4] LC_NUMERIC=C
    ## [5] LC_TIME=English_United States.utf8
    ##
    ## attached base packages:
    ## [1] stats4    stats     graphics  grDevices utils     datasets  methods
    ## [8] base
    ##
    ## other attached packages:
    ##  [1] SPIAT_1.0.4                 SpatialExperiment_1.8.0
    ##  [3] SingleCellExperiment_1.20.0 SummarizedExperiment_1.28.0
    ##  [5] Biobase_2.58.0              GenomicRanges_1.50.2
    ##  [7] GenomeInfoDb_1.34.4         IRanges_2.32.0
    ##  [9] S4Vectors_0.36.1            BiocGenerics_0.44.0
    ## [11] MatrixGenerics_1.10.0       matrixStats_0.63.0
    ## [13] BiocStyle_2.26.0
    ##
    ## loaded via a namespace (and not attached):
    ##   [1] colorspace_2.0-3          rjson_0.2.21
    ##   [3] deldir_1.0-6              ggridges_0.5.4
    ##   [5] circlize_0.4.15           scuttle_1.8.3
    ##   [7] XVector_0.38.0            GlobalOptions_0.1.2
    ##   [9] clue_0.3-63               spatstat.data_3.0-0
    ##  [11] farver_2.1.1              ggrepel_0.9.2
    ##  [13] bit64_4.0.5               fansi_1.0.3
    ##  [15] codetools_0.2-18          R.methodsS3_1.8.2
    ##  [17] sparseMatrixStats_1.10.0  doParallel_1.0.17
    ##  [19] cachem_1.0.6              knitr_1.41
    ##  [21] polyclip_1.10-4           mmand_1.6.2
    ##  [23] jsonlite_1.8.4            apcluster_1.4.10
    ##  [25] Cairo_1.6-0               cluster_2.1.4
    ##  [27] png_0.1-8                 R.oo_1.25.0
    ##  [29] pheatmap_1.0.12           HDF5Array_1.26.0
    ##  [31] spatstat.sparse_3.0-0     BiocManager_1.30.19
    ##  [33] compiler_4.2.2            dqrng_0.3.0
    ##  [35] assertthat_0.2.1          Matrix_1.5-3
    ##  [37] fastmap_1.1.0             limma_3.54.0
    ##  [39] cli_3.5.0                 htmltools_0.5.4
    ##  [41] tools_4.2.2               gtable_0.3.1
    ##  [43] glue_1.6.2                GenomeInfoDbData_1.2.9
    ##  [45] RANN_2.6.1                reshape2_1.4.4
    ##  [47] dplyr_1.0.10              Rcpp_1.0.9
    ##  [49] jquerylib_0.1.4           vctrs_0.5.1
    ##  [51] rhdf5filters_1.10.0       spatstat.explore_3.0-5
    ##  [53] nlme_3.1-161              iterators_1.0.14
    ##  [55] DelayedMatrixStats_1.20.0 spatstat.random_3.0-1
    ##  [57] xfun_0.35                 stringr_1.5.0
    ##  [59] beachmat_2.14.0           lifecycle_1.0.3
    ##  [61] gtools_3.9.4              goftest_1.2-3
    ##  [63] edgeR_3.40.1              zlibbioc_1.44.0
    ##  [65] scales_1.2.1              vroom_1.6.0
    ##  [67] spatstat.utils_3.0-1      parallel_4.2.2
    ##  [69] rhdf5_2.42.0              RColorBrewer_1.1-3
    ##  [71] ComplexHeatmap_2.14.0     yaml_2.3.6
    ##  [73] gridExtra_2.3             ggplot2_3.4.0
    ##  [75] sass_0.4.4                stringi_1.7.8
    ##  [77] highr_0.9                 foreach_1.5.2
    ##  [79] BiocParallel_1.32.4       shape_1.4.6
    ##  [81] rlang_1.0.6               pkgconfig_2.0.3
    ##  [83] bitops_1.0-7              pracma_2.4.2
    ##  [85] evaluate_0.19             lattice_0.20-45
    ##  [87] archive_1.1.5             tensor_1.5
    ##  [89] Rhdf5lib_1.20.0           labeling_0.4.2
    ##  [91] cowplot_1.1.1             bit_4.0.5
    ##  [93] tidyselect_1.2.0          plyr_1.8.8
    ##  [95] magrittr_2.0.3            bookdown_0.31
    ##  [97] R6_2.5.1                  magick_2.7.3
    ##  [99] generics_0.1.3            DelayedArray_0.24.0
    ## [101] DBI_1.1.3                 pillar_1.8.1
    ## [103] withr_2.5.0               abind_1.4-5
    ## [105] RCurl_1.98-1.9            tibble_3.1.8
    ## [107] crayon_1.5.2              DropletUtils_1.18.1
    ## [109] utf8_1.2.2                spatstat.geom_3.0-3
    ## [111] dittoSeq_1.10.0           tzdb_0.3.0
    ## [113] rmarkdown_2.19            GetoptLong_1.0.5
    ## [115] locfit_1.5-9.6            grid_4.2.2
    ## [117] digest_0.6.31             dbscan_1.1-11
    ## [119] R.utils_2.12.2            munsell_0.5.0
    ## [121] viridisLite_0.4.1         bslib_0.4.2

# 5 Author Contributions

AT, YF, TY, ML, JZ, VO, MD are authors of the package code. MD and YF wrote the vignette. AT, YF and TY designed the package.