Contents

1 Introduction

Here, we explain the way to interpret of HTML report generated by cellCellReport. If cellCellDecomp is properly finished, we can perform cellCellReport function to output the HTML report. The results can be confirmed by typing example(cellCellReport). The report will be generated in the temporary directory (it costs 5 to 10 minutes). The output directory contains some files and directories as follows.

Here, look at the index.html.

Figure1 : HTML report of cellCellReport

Figure1 : HTML report of cellCellReport

2 Interpretation of “1. About scTensor Algorithm”

In the HTML report, the 1st item describes the overview of scTensor and other CCI-related packages.

Figure2: 1. About scTensor Algorithm

Figure2: 1. About scTensor Algorithm

3 Interpretation of “2. Global statistics and plots”

The 2nd item describes all the R objects saved in reanalysis.RData, which contains the result of scTensor. This file is saved in the output directory (out.dir) specified in cellCellReport, and the user also can re-analyze the result of scTensor.

Figure3: 2. Global statistics and plots

Figure3: 2. Global statistics and plots

Using plotly package, cellCellReport generates some interactive plots. For example, in item 2.1, the number of cells in each cell type can be confirmed when the cursor moved on the box.

Figure4: 2.1 Number of cells in each celltype

Figure4: 2.1 Number of cells in each celltype

In item 2.2, the number of expressed genes in each cell type (Non-zero genes) can be confirmed when the cursor moved on the box. Figure5: 2.2 Number of expressed genes in each cell type (Non-zero genes)

In item 2.3, the two-dimensional plot user specified can be confirmed. Figure6: 2.3 Two dimensional plot of all cells

In item 2.4, the distribution of core tensor values and the value of each (Ligand-Cell-type, Receptor-Cell-type, LR-pair) pattern can be confirmed.

The red bars mean that these values are selected by the threshold (thr parameters) in cellCellReport.

Note that the thr can be specified from 0 to 100, the large thr value will generate too many HTML files (cf. 8. (Ligand-Cell, Receptor-Cell, LR-pair) Patterns) and takes a long time.

Figure7 : 2.4 Distribution of core tensor values

Figure7 : 2.4 Distribution of core tensor values

The 3-order CCI-tensor consisting of Cell_L \(\times\) Cell_R \(\times\) LR-pair (LR) are decomposed by nnTensor, in which the tensor is iteratively matricised to mode-1 (Ligand-Cell direction) and mode-2 (Receptor-Cell direction).

In each direction, NMF is performed and the strength of each directional pattern is summarized in the bar plots.

For example, in item 2.5, the distribution of mode-1 matricised tensor can be confirmed.

Figure8: 2.5 Distribution of mode-1 matricised tensor (Ligand-Cell Direction) (1/2)

Figure8: 2.5 Distribution of mode-1 matricised tensor (Ligand-Cell Direction) (1/2)

Figure9: 2.5 Distribution of mode-1 matricised tensor (Ligand-Cell Direction) (2/2)

Figure9: 2.5 Distribution of mode-1 matricised tensor (Ligand-Cell Direction) (2/2)

Likewise, in item 2.6, the distribution of mode-2 matricised tensor can be confirmed,

Figure10 : 2.6 Distribution of mode-2 matricised tensor (Receptor-Cell Direction) (1/2)

Figure10 : 2.6 Distribution of mode-2 matricised tensor (Receptor-Cell Direction) (1/2)

Figure11 : 2.6 Distribution of mode-2 matricised tensor (Receptor-Cell Direction) (2/2)

Figure11 : 2.6 Distribution of mode-2 matricised tensor (Receptor-Cell Direction) (2/2)

4 Interpretation of “3. Ligand-Cell Patterns”

In the 3rd item, using the heatmap of plotly, the user can interactively confirm the detail of Ligand-Cell Patterns extracted by nnTensor.

Figure14 : 3. Ligand-Cell Patterns

Figure14 : 3. Ligand-Cell Patterns

5 Interpretation of “4. Receptor-Cell Patterns”

Likewise, in the 4th item, the user can interactively confirm the detail of Receptor-Cell Patterns.

Figure15 : 4. Receptor-Cell Patterns

Figure15 : 4. Receptor-Cell Patterns

6 Interpretation of “5. CCI-wise Hypergraph”

In the 6th item describes, the strength between Ligand-Cell Patterns and Receptor-Cell Patterns (CCI-strength), by the summation of the core tensor with the mode-3 direction, a matrix consisting of the number of Ligand-Cell Patterns \(\times\) the number of Receptor-Cell Patterns.

Figure18 : 6. CCI-wise Hypergraph (1/2) Figure19 : 6. CCI-wise Hypergraph (2/2)

7 Interpretation of “6. Gene-wise Hypergraph”

In the 7th item, the relationship between LR-pairs, which coexpressed in any LR-pair pattern at least one time. Ligand genes are described as red nodes, receptor genes are described as blue nodes, and corresponding LR-pair patterns are described as the color of edges. Using visNetwork package, these interactions can be interactively visualized.

Figure20 : 7. Gene-wise Hypergraph

Figure20 : 7. Gene-wise Hypergraph

Under the gene-wise hypergraph, four hyperlinks are embedded.

In the 1st link, the details of the gene-wise hypergraph can be confirmed as a corresponding table in a ligand gene-centric manner. This page can work as a reverse lookup search by “Ctrl + F”; by typing the gene name of ligand that the user is interested in, the partner receptors, which are coexpressed in some LR-pair patterns, also can be found. Figure21: Details of Ligand Gene-centric Overview (selected)

In the 2nd link, the user can find all the partner receptors, even if the partner receptors are not coexpressed in any LR-pair pattern, and if they are not included in the data matrix. Figure22: Details of Ligand Gene-centric Overview (all)

Likewise, the receptor gene-centric reverse search page is embedded in the 3rd link, Figure23: Details of Receptor Gene-centric Overview (selected)

and, in the 4th link, all the partner ligand genes are included. Figure24: Details of Receptor Gene-centric Overview (all) (1/2) Figure25 : Details of Receptor Gene-centric Overview (all) (2/2)

8 Interpretation of “7. (Ligand-Cell, Receptor-Cell, LR-pair) Patterns”

In the 8th item, the details of (Ligand-Cell, Receptor-Cell, LR-pair) Patterns are ordered by the size of the core tensor, and the link of each pattern is embedded.

(Note that the number of links is dependent on the thr parameter of cellCellReport.) Figure26: 8. (Ligand-Cell, Receptor-Cell, LR-pair) Patterns

For example, the 1st link describes the details of (3,2,) Pattern, which means the relationship of 1st pattern of Ligand-Cell patterns, 1st pattern of Receptor-Cell patterns, and 5th pattern of LR-pair patterns. Figure27 : Details of (3,2,) Pattern (1/3)

In this pattern, only one LR-pair is coexpressed (INSL3 and GNG11). The hyperlinks to many databases and PubMed are also embedded. The degree of the size of the LR-pair in the LR-pair pattern is quantified as P-value and Q-value. Figure28 : Details of (3,2,) Pattern (2/3)

Under the LR-pair list, the results of many enrichment analysis are also embedded such as Gene Ontology (BP/MF/CC), Reactome, MeSH…etc. Figure29 : Details of (3,2,) Pattern (3/3)

User can confirm the detail of the result of scTensor, and perform the biological interpretation.

Session information

## R version 4.1.0 (2021-05-18)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.2 LTS
## 
## Matrix products: default
## BLAS:   /home/biocbuild/bbs-3.13-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.13-bioc/R/lib/libRlapack.so
## 
## locale:
##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=en_GB              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       
## 
## attached base packages:
## [1] parallel  stats4    stats     graphics  grDevices utils     datasets 
## [8] methods   base     
## 
## other attached packages:
##  [1] AnnotationHub_3.0.0                    
##  [2] BiocFileCache_2.0.0                    
##  [3] dbplyr_2.1.1                           
##  [4] scTGIF_1.6.0                           
##  [5] Homo.sapiens_1.3.1                     
##  [6] TxDb.Hsapiens.UCSC.hg19.knownGene_3.2.2
##  [7] org.Hs.eg.db_3.13.0                    
##  [8] GO.db_3.13.0                           
##  [9] OrganismDbi_1.34.0                     
## [10] GenomicFeatures_1.44.0                 
## [11] AnnotationDbi_1.54.0                   
## [12] LRBase.Mmu.eg.db_2.0.0                 
## [13] SingleCellExperiment_1.14.0            
## [14] SummarizedExperiment_1.22.0            
## [15] Biobase_2.52.0                         
## [16] GenomicRanges_1.44.0                   
## [17] GenomeInfoDb_1.28.0                    
## [18] IRanges_2.26.0                         
## [19] S4Vectors_0.30.0                       
## [20] BiocGenerics_0.38.0                    
## [21] MatrixGenerics_1.4.0                   
## [22] matrixStats_0.58.0                     
## [23] scTensor_2.2.0                         
## [24] RSQLite_2.2.7                          
## [25] LRBase.Hsa.eg.db_2.0.0                 
## [26] LRBaseDbi_2.2.0                        
## [27] BiocStyle_2.20.0                       
## 
## loaded via a namespace (and not attached):
##   [1] Hmisc_4.5-0                   ica_1.0-2                    
##   [3] Rsamtools_2.8.0               foreach_1.5.1                
##   [5] lmtest_0.9-38                 crayon_1.4.1                 
##   [7] spatstat.core_2.1-2           MASS_7.3-54                  
##   [9] nlme_3.1-152                  backports_1.2.1              
##  [11] GOSemSim_2.18.0               MeSHDbi_1.28.0               
##  [13] rlang_0.4.11                  XVector_0.32.0               
##  [15] ROCR_1.0-11                   irlba_2.3.3                  
##  [17] nnTensor_1.0.7                filelock_1.0.2               
##  [19] GOstats_2.58.0                BiocParallel_1.26.0          
##  [21] rjson_0.2.20                  tagcloud_0.6                 
##  [23] bit64_4.0.5                   glue_1.4.2                   
##  [25] sctransform_0.3.2             spatstat.sparse_2.0-0        
##  [27] dotCall64_1.0-1               spatstat.geom_2.1-0          
##  [29] tcltk_4.1.0                   DOSE_3.18.0                  
##  [31] tidyselect_1.1.1              SeuratObject_4.0.1           
##  [33] fitdistrplus_1.1-3            XML_3.99-0.6                 
##  [35] tidyr_1.1.3                   zoo_1.8-9                    
##  [37] GenomicAlignments_1.28.0      xtable_1.8-4                 
##  [39] magrittr_2.0.1                evaluate_0.14                
##  [41] ggplot2_3.3.3                 zlibbioc_1.38.0              
##  [43] rstudioapi_0.13               miniUI_0.1.1.1               
##  [45] bslib_0.2.5.1                 rpart_4.1-15                 
##  [47] fastmatch_1.1-0               ensembldb_2.16.0             
##  [49] treeio_1.16.0                 maps_3.3.0                   
##  [51] fields_12.3                   shiny_1.6.0                  
##  [53] xfun_0.23                     cluster_2.1.2                
##  [55] tidygraph_1.2.0               TSP_1.1-10                   
##  [57] KEGGREST_1.32.0               tibble_3.1.2                 
##  [59] interactiveDisplayBase_1.30.0 ggrepel_0.9.1                
##  [61] biovizBase_1.40.0             ape_5.5                      
##  [63] listenv_0.8.0                 dendextend_1.15.1            
##  [65] Biostrings_2.60.0             png_0.1-7                    
##  [67] future_1.21.0                 withr_2.4.2                  
##  [69] bitops_1.0-7                  ggforce_0.3.3                
##  [71] RBGL_1.68.0                   plyr_1.8.6                   
##  [73] GSEABase_1.54.0               AnnotationFilter_1.16.0      
##  [75] pillar_1.6.1                  cachem_1.0.5                 
##  [77] graphite_1.38.0               vctrs_0.3.8                  
##  [79] ellipsis_0.3.2                generics_0.1.0               
##  [81] plot3D_1.3                    MeSH.Aca.eg.db_1.15.0        
##  [83] outliers_0.14                 tools_4.1.0                  
##  [85] foreign_0.8-81                entropy_1.3.0                
##  [87] munsell_0.5.0                 tweenr_1.0.2                 
##  [89] fgsea_1.18.0                  DelayedArray_0.18.0          
##  [91] fastmap_1.1.0                 compiler_4.1.0               
##  [93] abind_1.4-5                   httpuv_1.6.1                 
##  [95] rtracklayer_1.52.0            Gviz_1.36.0                  
##  [97] plotly_4.9.3                  GenomeInfoDbData_1.2.6       
##  [99] gridExtra_2.3                 lattice_0.20-44              
## [101] deldir_0.2-10                 visNetwork_2.0.9             
## [103] AnnotationForge_1.34.0        utf8_1.2.1                   
## [105] later_1.2.0                   dplyr_1.0.6                  
## [107] jsonlite_1.7.2                concaveman_1.1.0             
## [109] scales_1.1.1                  graph_1.70.0                 
## [111] tidytree_0.3.3                pbapply_1.4-3                
## [113] genefilter_1.74.0             lazyeval_0.2.2               
## [115] promises_1.2.0.1              MeSH.db_1.15.0               
## [117] latticeExtra_0.6-29           goftest_1.2-2                
## [119] spatstat.utils_2.1-0          reticulate_1.20              
## [121] checkmate_2.0.0               rmarkdown_2.8                
## [123] cowplot_1.1.1                 schex_1.6.0                  
## [125] MeSH.Syn.eg.db_1.15.0         webshot_0.5.2                
## [127] Rtsne_0.15                    dichromat_2.0-0              
## [129] BSgenome_1.60.0               uwot_0.1.10                  
## [131] igraph_1.2.6                  survival_3.2-11              
## [133] yaml_2.2.1                    plotrix_3.8-1                
## [135] htmltools_0.5.1.1             memoise_2.0.0                
## [137] VariantAnnotation_1.38.0      rTensor_1.4.8                
## [139] BiocIO_1.2.0                  Seurat_4.0.1                 
## [141] seriation_1.2-9               graphlayouts_0.7.1           
## [143] viridisLite_0.4.0             digest_0.6.27                
## [145] assertthat_0.2.1              ReactomePA_1.36.0            
## [147] mime_0.10                     rappdirs_0.3.3               
## [149] registry_0.5-1                spam_2.6-0                   
## [151] future.apply_1.7.0            misc3d_0.9-0                 
## [153] data.table_1.14.0             blob_1.2.1                   
## [155] cummeRbund_2.34.0             splines_4.1.0                
## [157] Formula_1.2-4                 ProtGenerics_1.24.0          
## [159] RCurl_1.98-1.3                hms_1.1.0                    
## [161] colorspace_2.0-1              base64enc_0.1-3              
## [163] BiocManager_1.30.15           aplot_0.0.6                  
## [165] nnet_7.3-16                   sass_0.4.0                   
## [167] Rcpp_1.0.6                    bookdown_0.22                
## [169] RANN_2.6.1                    MeSH.PCR.db_1.15.0           
## [171] enrichplot_1.12.0             fansi_0.4.2                  
## [173] parallelly_1.25.0             R6_2.5.0                     
## [175] grid_4.1.0                    ggridges_0.5.3               
## [177] lifecycle_1.0.0               curl_4.3.1                   
## [179] MeSH.Bsu.168.eg.db_1.15.0     leiden_0.3.7                 
## [181] MeSH.AOR.db_1.15.0            meshr_1.28.0                 
## [183] jquerylib_0.1.4               DO.db_2.9                    
## [185] Matrix_1.3-3                  qvalue_2.24.0                
## [187] RcppAnnoy_0.0.18              RColorBrewer_1.1-2           
## [189] iterators_1.0.13              stringr_1.4.0                
## [191] htmlwidgets_1.5.3             polyclip_1.10-0              
## [193] biomaRt_2.48.0                purrr_0.3.4                  
## [195] shadowtext_0.0.8              reactome.db_1.76.0           
## [197] mgcv_1.8-35                   globals_0.14.0               
## [199] htmlTable_2.2.1               patchwork_1.1.1              
## [201] codetools_0.2-18              prettyunits_1.1.1            
## [203] gtable_0.3.0                  DBI_1.1.1                    
## [205] highr_0.9                     tensor_1.5                   
## [207] httr_1.4.2                    KernSmooth_2.23-20           
## [209] stringi_1.6.2                 progress_1.2.2               
## [211] msigdbr_7.4.1                 reshape2_1.4.4               
## [213] farver_2.1.0                  heatmaply_1.2.1              
## [215] annotate_1.70.0               viridis_0.6.1                
## [217] hexbin_1.28.2                 fdrtool_1.2.16               
## [219] Rgraphviz_2.36.0              magick_2.7.2                 
## [221] ggtree_3.0.0                  rvcheck_0.1.8                
## [223] restfulr_0.0.13               Category_2.58.0              
## [225] scattermore_0.7               BiocVersion_3.13.1           
## [227] bit_4.0.4                     spatstat.data_2.1-0          
## [229] scatterpie_0.1.6              jpeg_0.1-8.1                 
## [231] ggraph_2.0.5                  babelgene_21.4               
## [233] pkgconfig_2.0.3               MeSH.Hsa.eg.db_1.15.0        
## [235] knitr_1.33