Welcome to online tools of bulk or scRNA-seq data.

Please wait for the page to load completely, it may take a few seconds to initialize ShinyApp.

Loading...

2.Please Select a Dataset of Interest Based on the Information Mentioned Above:


3.Initial Parameters for Bulk/Single-cell RNA-seq Data:

Identification of Differentially Expressed Genes


Powered by Limma



Plot

Loading...

Plot Download Options

Download Plot

Weighted Gene Co-expression Network Analysis



Filter Cells within Thresholds


Filter Cells

Seurat allows you to easily explore QC metrics and filter cells based on any user-defined criteria. You can visualize gene and molecule counts, plot their relationship, and exclude cells with a clear outlier number of genes detected as potential multiplets. This is not a guaranteed method to exclude cell doublets, see tutorial for more information. As an example in the tutorial, you can filter cells based on the percentage of mitochondrial genes present.

Scatter Plot

Scatter Plot is typically used to visualize feature-feature relationships. See tutorial for examples/details

Violin Plot

Note that low.thresholds and high.thresholds are used to define a 'gate'

Select thresholds to filter cells

Loading...
01,00002004006008001,000
01,00002004006008001,000
10.200.40.8
10.200.40.8

Number Of Genes & Cells


Data Normalization & Detection of Variable Genes


SCTransform function

Use this function as an alternative to the NormalizeData, FindVariableFeatures, ScaleData workflow. Results are saved in a new assay called SCT with counts being (corrected) counts, data being log1p(counts), scale.data being pearson residuals; sctransform::vst intermediate results are saved in misc slot of new assay.

We can regress out cell-cell variation in gene expression driven by batch (if applicable), cell alignment rate (as provided by Drop-seq tools for Drop-seq data), the number of detected molecules, and mitochondrial gene expression. Refer to tutorial to see an example of regressing on the number of detected molecules per cell as well as the percentage mitochondrial gene content for post-mitotic blood cells.

Variables to regress out in a second non-regularized linear regression. For example, percent.mito.


Data Normalization

After removing unwanted cells from the dataset, the next step is to normalize the data.

By default, we employ a global-scaling normalization method “LogNormalize” that normalizes the gene expression measurements for each cell by the total expression, multiplies this by a scale factor (10,000 by default), and log-transforms the result.

LogNormalize

Detection of variable genes across the single cells

Seurat calculates highly variable genes and focuses on these for downstream analysis. FindVariableGenes calculates the average expression and dispersion for each gene, places these genes into bins, and then calculates a z-score for dispersion within each bin. This helps control for the relationship between variability and average expression.

We suggest that users set these parameters to mark visual outliers on the dispersion plot, but the exact parameter settings may vary based on the data type, heterogeneity in the sample, and normalization strategy.

vst

Scaling the data and removing unwanted sources of variation

Your single cell dataset likely contains ‘uninteresting’ sources of variation. This could include not only technical noise, but batch effects, or even biological sources of variation (cell cycle stage). As suggested in Buettner et al, NBT, 2015, regressing these signals out of the analysis can improve downstream dimensionality reduction and clustering. To mitigate the effect of these signals, Seurat constructs linear models to predict gene expression based on user-defined variables. The scaled z-scored residuals of these models are stored in the scale.data slot, and are used for dimensionality reduction and clustering.

We can regress out cell-cell variation in gene expression driven by batch (if applicable), cell alignment rate (as provided by Drop-seq tools for Drop-seq data), the number of detected molecules, and mitochondrial gene expression. Refer to tutorial to see an example of regressing on the number of detected molecules per cell as well as the percentage mitochondrial gene content for post-mitotic blood cells.


Dispersion Plot

The average expression and dispersion for each gene

Loading...

Principal Component Analysis Reduction

Compute PCA

Next we perform PCA on the scaled data. By default, the variable features are used as input, but can be defined using features.

We have typically found that running dimensionality reduction on highly variable genes can improve performance.

However, with UMI data - particularly after regressing out technical variables, we often see that PCA returns similar (albeit slower) results when run on much larger subsets of genes, including the whole transcriptome.


Compute ICA (Optional)


PCA Print Output

visualizing cells and features that define the PCA

Seurat provides several useful ways of visualizing both cells and features that define the PCA, including VizDimLoadings, DimPlot, and DimHeatmap

VizDimLoadings allows to visualize top genes associated with reduction components.

DimPlot allows to visualize the cells in the 2 principal components

DimHeatmap allows for easy exploration of the primary sources of heterogeneity in a dataset, and can be useful when trying to decide which PCs to include for further downstream analyses. Both cells and features are ordered according to their PCA scores

Loading...
Loading...
Loading...

ICA Print Output

visualizing cells and features that define the PCA

Seurat provides a useful way of visualizing both cells and features that define the ICA

DimPlotIca allows to visualize the cells in the 2 principal components

Loading...


Determining Significant PCs


To overcome the extensive technical noise in any single gene for scRNA-seq data, Seurat clusters cells based on their PCA scores, with each PC essentially representing a ‘metagene’ that combines information across a correlated gene set. Determining how many PCs to include downstream is therefore an important step.

PC selection – identifying the true dimensionality of a dataset – is an important step for Seurat, but can be challenging/uncertain for the user. We therefore suggest these three approaches to consider.

The first is more supervised, exploring PCs to determine relevant sources of heterogeneity, and could be used in conjunction with GSEA for example.

The second implements a statistical test based on a random null model, but is time-consuming for large datasets, and may not return a clear PC cutoff.

The third is a heuristic that is commonly used, and can be calculated instantly.


1) PC Elbow Plot (quick)

A more ad hoc method for determining which PCs to use is to look at a plot of the standard deviations of the principle components and draw your cutoff where there is a clear elbow in the graph. This can be done with PCElbowPlot.

Loading...

Plot Download Options

Download Plot

2) JackStraw (optional)

In Macosko et al, we implemented a resampling test inspired by the jackStraw procedure. We randomly permute a subset of the data (1% by default) and rerun PCA, constructing a ‘null distribution’ of gene scores, and repeat this procedure. We identify ‘significant’ PCs as those who have a strong enrichment of low p-value genes.

The JackStrawPlot function provides a visualization tool for comparing the distribution of p-values for each PC with a uniform distribution (dashed line). ‘Significant’ PCs will show a strong enrichment of genes with low p-values (solid curve above the dashed line).

NOTE: this process will take a long time for big datasets. It can take ~10 minutes.


PCs to use:


Loading...

Plot Download Options

Download Plot

Cell Clustering

Seurat now includes an graph-based clustering approach. Importantly, the distance metric which drives the clustering analysis (based on previously identified PCs) remains the same. However, our approach to partioning the cellular distance matrix into clusters has dramatically improved. Our approach was heavily inspired by recent manuscripts which applied graph-based clustering approaches to scRNA-seq data [SNN-Cliq, Xu and Su, Bioinformatics, 2015] and CyTOF data [PhenoGraph, Levine et al., Cell, 2015].

Briefly, these methods embed cells in a graph structure - for example a K-nearest neighbor (KNN) graph, with edges drawn between cells with similar gene expression patterns, and then attempt to partition this graph into highly interconnected ‘quasi-cliques’ or ‘communities’. As in PhenoGraph, we first construct a KNN graph based on the euclidean distance in PCA space, and refine the edge weights between any two cells based on the shared overlap in their local neighborhoods (Jaccard distance). To cluster the cells, we apply modularity optimization techniques [SLM, Blondel et al., Journal of Statistical Mechanics], to iteratively group cells together, with the goal of optimizing the standard modularity function.


FindClusters paramters

The FindClusters function implements the procedure, and contains a resolution parameter that sets the ‘granularity’ of the downstream clustering, with increased values leading to a greater number of clusters. We find that setting this parameter between 0.6-1.2 typically returns good results for single cell datasets of around 3K cells. Optimal resolution often increases for larger datasets. The clusters are saved in the object@ident slot.

pca
1
2
3
4
5


Cluster tree

At this step, a cluster tree is performed on the scaled data. Only the selected dimension are used as input.

TSNE/UMAP Reduction


Cell Type Annotation


Next Step or Download





Cell Clustering Visualization


Plot of cell clustering

Plot of gene(set) expression

Plot of select cells

Table of select cells

Download


You can save the object at this point so that it can easily be loaded back in R for further analysis & exploration without having to rerun the computationally intensive steps performed above, or easily shared with collaborators.

It is also recommended that you keep it as a reference.


Generate and Download the R script to reproduce these steps in R/RStudio

Please note that you need to edit the data file(s)/directory path in the script before you run it in R/RStudio

Identification of Cluster Markers


Seurat can help you find markers that define clusters via differential expression. By default, it identifes positive and negative markers of a single cluster (specified in ident.1), compared to all other cells. FindAllMarkers automates this process for all clusters, but you can also test groups of clusters vs. each other, or against all cells.

The min.pct argument requires a gene to be detected at a minimum percentage in either of the two groups of cells, and the thresh.test argument requires a gene to be differentially expressed (on average) by some amount between the two groups. You can set both of these to 0, but with a dramatic increase in time - since this will test a large number of genes that are unlikely to be highly discriminatory. As another option to speed up these computations, max.cells.per.ident can be set. This will downsample each identity class to have no more cells than whatever this is set to. While there is generally going to be a loss in power, the speed increases can be significiant and the most highly differentially expressed genes will likely still rise to the top.

Find All Markers:

wilcox

Loading...

Visualizing Marker Expression:

Note: make sure to find cluster markers first, as the list of genes you can plot will from the output of that step

Feature Plot:

UMAP

Plot Download Options

Download Plot
Loading...

Gene Enrichment Analysis


Select a gene enrichment result


Cell clustering visualization


About Gene Expression Nebulas Tools