Skip to main content
SimilarExport ZIP
100%

Use Seurat to find specific markers of each cell type

✓ Published0🌍 Public
FFloWuenne
Last edited Oct 30, 2017
Created on Oct 30, 2017

This example demonstrates how to identify specific marker genes for each cell type in a Seurat single-cell RNA-seq analysis. It shows a workflow that compares every cell type against all others using the `FindMarkers` function with a ROC test, then filters the results to retain only markers expressed in fewer than 50% of cells in all other populations. The code uses the `tidyverse` and `Seurat` libraries to perform pairwise comparisons, calculate the difference in expression percentages, and summarize the average AUC values. The final output is a curated list of the top 10 markers per cell type, emphasizing those with high discriminatory power.

AI-generated description

Similar vizzes

Loading thumbnail…

Use Seurat to find specific markers of each cell type (modified from Dan Skelly)

This example shows a dot plot of the top five specific marker genes for each cell type in a single-cell RNA-sequencing dataset, with dot size and color encoding expression level and percentage of cells expressing each gene. The visualization relies on Seurat’s `FindMarkers` function with a ROC test to perform all pairwise comparisons between cell types, filtering for markers expressed in fewer than 50% of cells in other clusters and retaining those with a mean AUC above 0.65. The code uses tidyverse functions like `mutate`, `filter`, `summarize`, and `group_by` to process the marker data, and Seurat’s `DotPlot` to render the final figure. Marker lists are written to TSV files, and the plot is saved as a PNG.

FFloWuenne
89% match
Loading thumbnail…

Extract shiny clustering data for full dataset

This example prepares a clustered single-cell dataset for interactive visualization by extracting tSNE coordinates, cluster identities, normalized expression, and marker genes from a Seurat object. It shows how the full dataset’s clustering structure is assembled into a custom S4 object for downstream use. The code uses R with `dplyr` and `tidyr` to filter the top 50 markers per cluster by `avg_logFC`, and merges tSNE embeddings with cell identities. It saves the processed data as an RDS file for a Shiny app, relying on Seurat’s data structures and base R matrix operations.

FFloWuenne
65% match
Loading thumbnail…

Make a table of the number of cells for each group in the seurat object

The example builds a data frame of cell counts for each cluster in a Seurat object, broken down by treatment group. It first shows raw cell numbers per group per cluster, then adds normalized proportions relative to each group’s total cell count, making group compositions comparable despite unequal sample sizes. The code uses `Seurat` functions like `SubsetData` and the `@ident` and `@meta.data` slots, alongside `dplyr`’s `select` and `nrow` to count group memberships. It iterates over clusters, appending rows via `rbind`, and finally writes the result to a tab-separated file with `write.table`. The output table directly reveals per-cluster group distributions and their normalized percentages, aiding differential abundance analysis.

FFloWuenne
63% match
Loading thumbnail…

Create a cellview Rds object from a seurat expression object (updated for Seurat version 2) [MOUSE]

This example demonstrates how to prepare a Seurat single-cell expression object for use with the CellView visualization tool, specifically for mouse data. It shows the workflow of rerunning t-SNE with three dimensions, converting gene symbols to Ensembl IDs using `bitr` from clusterProfiler, and constructing the necessary `log2cpm`, `featuredata`, and `tsne.data` matrices. The code uses Seurat’s `RunTSNE`, the `Matrix` and `dplyr` packages for data manipulation, and ultimately saves the combined object as an Rds file for downstream visualization.

FFloWuenne
61% match
Loading thumbnail…

Create shiny cluster modules for renamed datasets

This example demonstrates how to process and store single-cell RNA sequencing data across multiple developmental time points, creating Shiny-ready cluster modules for each stage. The code defines an S4 class, `clustering_module`, that bundles t-SNE embeddings, normalized expression data, metadata, and a marker list for each time point from E14.5 to P7. It loads Seurat objects, extracts the t-SNE coordinates and cell identities, filters markers by expression difference using `dplyr`’s `top_n`, and saves each time point as a compressed RDS file. The workflow uses `Seurat` for data handling and `tidyr` for reshaping, preparing interactive visualizations for the Shiny server.

FFloWuenne
61% match
Loading thumbnail…

Cell cycle assignment for single cell data using scran and cyclone with gene symbols instead of ensembl IDs

This example demonstrates how to assign cell cycle phases to single-cell RNA sequencing data by converting the pretrained marker gene pairs in the scran package from Ensembl IDs to gene symbols. It shows the resulting phase assignments being added as metadata to a Seurat object for downstream analysis. The code uses the scran cyclone function, the org.Hs.eg.db database for ID conversion via bitr, and dplyr for data manipulation. The visualization implicitly relies on the Seurat object’s structure to store the new cell cycle labels, though no plot is generated in the provided script.

FFloWuenne
61% match
Loading thumbnail…

Create an Rds file for the full heart maturation dataset

This example demonstrates how to convert a Seurat single-cell RNA sequencing object into a compact RDS file for a Shiny visualization app. It shows the creation of a custom S4 class that bundles tSNE coordinates, normalized expression data, cell metadata, and the top 50 differentially expressed markers per cluster. Using the dplyr and tidyr libraries for data manipulation, the script extracts the tSNE embeddings and cluster identities from the Seurat object, then merges them with metadata. The resulting clustering_module object is saved via saveRDS, enabling efficient downstream interactive exploration of heart maturation data.

FFloWuenne
60% match
Loading thumbnail…

Plot VlnPlot for markers with minimum axis for joined marker plot as in (Macosko et al for retina)

This example uses Seurat’s `VlnPlot` to draw violin plots for a set of candidate marker genes (GAPDH, NANOG, ACTB) across cell populations, with a fixed y-max of 6 and no axis labels or ticks to focus on expression distributions. The code loops through each gene, builds a violin plot, and then removes all axis lines and text while flipping coordinates so the plots are horizontal. It demonstrates how to customize Seurat’s plotting defaults—using `point.size.use = -1` to suppress points, `remove.legend = TRUE`, and theme elements to create a minimal, clean visualization for comparing marker expression across clusters. The data source is an updated Seurat object (`expression_seurat.updated`).

FFloWuenne
60% match