SeuratExplorer

Why build this R package

Currently, there is still no good tools for visualising the analysis results from Seurat, when the bioinformatics analyst hands over the results to the user, if the user does not have any R language foundation, it is still difficult to retrieve the results and re-analysis on their own, and this R package is designed to help such users to visualize and explore the anaysis results. The only thing to do for such users is to configure R and Rstudio on their own computers, and then install SeuratExplorer, without any other operations, an optional way is to upload the Seurat object file to a server which has been deployed with shinyserver and SeuratExplorer.

Essentially, what SeuratExplorer done is just to perform visual operations for command line tools from Seurat or other packages.

A live demo webserver

Upload an Rds or qs2 file, with file size no more than 5GB. Open, you can download a mini demo data from github here.

Installation

You can use codes bellow to install the latest version of SeuratExplorer:

# install dependency
if (!require("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("ComplexHeatmap")

# install SeuratExplorer
if(!require(devtools)){install.packages("devtools")}
install_github("fentouxungui/SeuratExplorer")

Run App

library(SeuratExplorer)
launchSeuratExplorer()

Introduction

Load data

• support Seurat object saved as .rds or .qs2 file.

• support data processed by Seurat V5 and older versions. it may takes a while to update Seurat object when loading data.

Cell Metadata

• support download cell metadata in csv format, which can be used to for further analysis.

Dimensional Reduction Plot

• support options for Dimension Reductions

• support options for Cluster Resolution

• support split plots

• support adjust the height/width ratio of the plot

• support options for showing cluster label

• support adjust label size

• support adjust point size

• support download plot in pdf format, what you see is what you get

Example plots:

Feature Plot

• support display multiple genes simultaneous, genes names are case-insensitive. Tips: paste multiple genes from excel

• support options for Dimension Reductions

• support split plots

• support change colors for the lowest expression and highest expression

• support adjust the height/width ratio of the plot

• support adjust point size

• support download plot in pdf format, what you see is what you get

Example plots:

Violin Plot

• support display multiple genes simultaneous, genes names are case-insensitive. Tips: paste multiple genes from excel

• support options for Cluster Resolution

• support split plots

• support stack and flip plot, and color mapping selection.

• support adjust point size and transparency

• support adjust font size on x and y axis

• support adjust the height/width ratio of the plot

• support download plot in pdf format, what you see is what you get

Example plots:

Dot Plot

• support display multiple genes simultaneous, genes names are case-insensitive. Tips: paste multiple genes from excel

• support options for Cluster Resolution and subset clusters

• support split plots

• support cluster clusters

• support rotate axis and flip coordinate

• support adjust point size and transparency

• support adjust font size on x and y axis

• support adjust the height/width ratio of the plot

• support download plot in pdf format, what you see is what you get

Example plots:

Heatmap for cell level expression

• support display multiple genes simultaneous, genes names are case-insensitive. Tips: paste multiple genes from excel

• support options for Cluster Resolution and reorder clusters

• support adjust font size and rotation angle of cluster label, and flip coordinate

• support adjust the height of group bar

• support adjust the gap size between groups

• support adjust the font size of gene names

• support adjust the height/width ratio of the plot

• support download plot in pdf format, what you see is what you get

Example plots:

Heatmap for group averaged expression

• support display multiple genes simultaneous, genes names are case-insensitive. Tips: paste multiple genes from excel

• support options for Cluster Resolution and reorder clusters

• support adjust font size and rotation angle of cluster label

• support adjust the font size of gene names

• support adjust the height/width ratio of the plot

• support download plot in pdf format, what you see is what you get

Example plots:

Ridge Plot

• support display multiple genes simultaneous, genes names are case-insensitive. Tips: paste multiple genes from excel

• support options for Cluster Resolution and reorder clusters

• support adjust column numbers

• support stack plot and color mapping

• support adjust font size on x and y axis

• support adjust the height/width ratio of the plot

• support download plot in pdf format, what you see is what you get

Example plots:

Plot Cell Percentage

• support facet

• support adjust the height/width ratio of the plot

• support download plot in pdf format, what you see is what you get

Example plots:

Find Cluster Markers and DEGs Analysis

This usually takes longer, please wait patiently.Please save the results before start a new analysis, the old results will be overwritten by the new results, the results can be downloaded as csv format.

Support two ways

• support find markers for all clusters

• support calculate DEGs for self-defined two groups, you can subset cells before calculate DEGs between two groups, default use all cells of two groups.

You can modify part calculation parameters before a analysis.

Screen shots:

Output description

FindMarkers(object, …)

A data.frame with a ranked list of putative markers as rows, and associated statistics as columns (p-values, ROC score, etc., depending on the test used (test.use)). The following columns are always present:

avg_logFC: log fold-chage of the average expression between the two groups. Positive values indicate that the gene is more highly expressed in the first group

pct.1: The percentage of cells where the gene is detected in the first group

pct.2: The percentage of cells where the gene is detected in the second group

p_val_adj: Adjusted p-value, based on bonferroni correction using all genes in the dataset

Top Expressed Features

Highly expressed genes can reflect the main functions of cells, there two ways to do this. the first - Find Top Genes by Cell could find gene only high express in a few cells, while the second - Find Top Genes by Accumulated UMI counts is biased to find the highly expressed genes in most cells by accumulated UMI counts.

1. Find Top Genes by Cell

How?

Step1: for each cell, find genes that has high UMI percentage, for example, if a cell has 10000 UMIs, and the UMI percentage cutoff is set to 0.01, then all genes that has more than 10000 * 0.01 = 100 UMIs is thought to be the highly expressed genes for this cell.

Step2: summary those genes for each cluster, firstly get all highly expressed genes in a cluster, some genes may has less cells, then for each gene, count cells in which this genes is highly expressed, and also calculate the mean and median UMI percentage in those highly expressed cells.

Output description

• celltype: the cluster name which is define by Choose A Cluster Resolution

• total.cells: total cell in this cluster

• Gene: this Gene is highly expressed in at least 1 cell in this cluster

• total.pos.cells: how many cells express this gene

• total.UMI.pct: (all UMIs of this gene)/(total UMIs of this cluster)

• cut.Cells: how many cells highly express this gene

• cut.pct.mean: in those highly expressed cells, the mean expression percentage

• cut.pct.median: in those highly expressed cells, the median expression percentage

2. Find Top Genes by Accumulated UMI counts

for each cluster, calculate the top n highly expressed genes by accumulated UMI counts. if a cluster has less than 3 cells, this cluster will be escaped.

Output description

• CellType: the cluster name which is define by Choose A Cluster Resolution

• total.cells: total cell in this cluster

• Gene: the top n highly expressed genes

• total.pos.cells: how many cells express this gene

• NormalizedAccumulatedUMICounts: (total accumulated UMI counts) / (total cells of this cluster)

• PCT: (total accumulated UMI counts) / (total UMIs of this cluster)

Feature Summary

Summary interested features by cluster, such as the positive cell percentage and mean/median expression level.

Output description

• celltype: the cluster name which is define by Choose A Cluster Resolution

• TotalCells: total cell in this cluster

• Gene: the input genes

• PCT: the percentage of how many cells express this gene in this cluster

• Expr.mean: the mean normalized expression in this cluster

• Expr.median: the median normalized expression in this cluster

Feature Correlation Analysis

Can calculate the correlation value of gene pairs within cells from a cluster, support pearson & spearman methods.

3 ways to do

• Find Top Correlated Gene Pairs: to find top 1000 correlated gene pairs

• Find Correlated Genes for A Gene: to find the most correlated genes for input genes

• Calculate Correlation for A Gene List: to calculate the correlation value for each pair of the input genes

Output description

• GeneA: the first gene in a Gene pair

• GeneB: the second gene in a Gene pair

• correlation: the correlation value

if nothing return, this is because the input genes has very low expression level, very low expressed genes will be removed before analysis.

Key related packages

• Hla-Lab/SeuratExplorer: An interactive R shiny application for exploring scRNAseq data processed in Seurat

• junjunlab/scRNAtoolVis: Some useful function to make your scRNA-seq plot more beautiful

• rstudio/shiny-server: Shiny Server is a server program that makes Shiny applications available over the web

Session Info

#> R version 4.4.1 (2024-06-14 ucrt)
#> Platform: x86_64-w64-mingw32/x64
#> Running under: Windows 11 x64 (build 22631)
#> 
#> Matrix products: default
#> 
#> 
#> locale:
#> [1] LC_COLLATE=Chinese (Simplified)_China.utf8 
#> [2] LC_CTYPE=Chinese (Simplified)_China.utf8   
#> [3] LC_MONETARY=Chinese (Simplified)_China.utf8
#> [4] LC_NUMERIC=C                               
#> [5] LC_TIME=Chinese (Simplified)_China.utf8    
#> 
#> time zone: Asia/Shanghai
#> tzcode source: internal
#> 
#> attached base packages:
#> [1] stats     graphics  grDevices utils     datasets  methods   base     
#> 
#> loaded via a namespace (and not attached):
#>  [1] compiler_4.4.1    fastmap_1.2.0     cli_3.6.3         tools_4.4.1      
#>  [5] htmltools_0.5.8.1 rstudioapi_0.16.0 yaml_2.3.8        rmarkdown_2.27   
#>  [9] highr_0.11        knitr_1.47        xfun_0.45         digest_0.6.36    
#> [13] rlang_1.1.4       evaluate_0.24.0