Skip to main content
SimilarExport ZIP
100%

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

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

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.

AI-generated description

Similar vizzes

Loading thumbnail…

Use Seurat to find specific markers of each cell type

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.

FFloWuenne
89% 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
68% 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
68% match
Loading thumbnail…

Make a nice ClusterTree from the initial ClusterTree plot from Seurat

This example shows a cluster tree visualization derived from a Seurat single-cell analysis, displaying the hierarchical relationships between cell clusters as a phylogenetic-style tree. The plot renders a ladderized dendrogram where nodes are labeled with grey circles, making the branching structure easier to interpret than the default Seurat output. The code uses Seurat’s `BuildClusterTree` function to compute the tree, then plots it with base R’s `plot` and `nodelabels` functions. The visualization applies `ladderize` to reorder branches for clarity and sets custom graphical parameters like `edge.width` and `label.offset` to improve readability. The data source is the Seurat object’s cluster tree, saved as a PNG file.

FFloWuenne
66% match
Loading thumbnail…

Plot mean expression of a gene vs other parameters and highlight specific genes using ggrepel

This example uses the scater package to calculate quality metrics for single-cell RNA-seq data, then plots mean gene expression against percent dropout and percent total expression. The ggplot2 code highlights a user-defined set of genes of interest in red, while ggrepel’s `geom_text_repel` labels those genes to avoid overlapping points. Two separate scatterplots are generated from the same `gene_info` data frame, each saved as an SVG, showing how specific genes relate to global expression distributions. The analysis relies on the `newSCESet` and `calculateQCMetrics` functions from scater to derive the plotted metrics.

FFloWuenne
65% 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
64% 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
64% 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
64% match