seurat-analysis

Analyze single-cell and spatial transcriptomics data with Seurat v5 on FGCZ infrastructure.

Updated May 28, 2026
One-click install
npx skills add https://github.com/cpanse/skills --skill seurat-analysis-cpanse
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: seurat-analysis
Source: https://github.com/cpanse/skills/tree/main/single-cell-spatial-general/skills/seurat-analysis
Command: npx skills add https://github.com/cpanse/skills --skill seurat-analysis-cpanse

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill includes scripts (resource) and references (resource) and assets (resource) components.

What problem does it solve? Single-cell RNA-seq analysis involves many error-prone steps—QC, normalization, clustering, integration, annotation, and reporting—and FGCZ analysts must additionally follow strict institutional conventions for plots, storage paths, and R Markdown reports. This Skill encodes those workflows and standards so analyses are reproducible and consistent. ## Core Features & Use Cases - Standard scRNA-seq pipeline: Load data from h5ad, 10x CellRanger, or qs2 objects, then run QC, normalization, PCA, UMAP, clustering, and differential expression with an automated R script. - Multi-sample integration: Integrate datasets with Harmony (recommended), CCA, RPCA, FastMNN, or scVI, including parameter tuning and troubleshooting guidance. - Advanced analyses: Compositional analysis with sccomp, pseudo-bulk PCA, MSE distance analysis, Augur cell-type prioritization, and CytoTRACE 2 potency scoring. - FGCZ-standard reporting: Generate R Markdown HTML reports with the ezRun template, polychrome cluster colors, coord_fixed UMAPs, interactive DEG tables, and 300 DPI figures. - Use Case: A bioinformatician receives CellRanger outputs for project p12345, runs the QC pipeline script, integrates samples with Harmony, annotates cell types, and delivers an FGCZ-branded HTML report saved to gStore via g-req. ## Quick Start Ask the agent to load your h5ad or CellRanger data and run the Seurat QC and clustering pipeline following FGCZ conventions, then generate an R Markdown report.

Frequently Asked Questions about seurat-analysis

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I run a standard Seurat v5 single-cell analysis pipeline?

Load data from h5ad, 10x CellRanger, or qs2 files, then run QC filtering, normalization, variable feature selection, PCA, UMAP, and graph-based clustering. The included seurat_qc_pipeline.R script automates these steps via Rscript on the FGCZ cluster.

Harmony vs CCA vs RPCA for Seurat integration—which should I use?

Harmony is recommended for most cases: it is fast, memory-efficient, works directly on PCA embeddings, and handles large datasets over 100k cells. Use CCA or RPCA for strong, complex batch effects or cross-species integration requiring anchor-based correction.

How do I load h5ad files into Seurat?

Use the anndataR package: read_h5ad("/path/to/file.h5ad", as = "Seurat") converts the h5ad file directly into a Seurat object. Alternatively, load 10x outputs with Read10X or saved objects efficiently with qs2's qs_read.

Why is my Seurat v5 analysis failing after integration?

Seurat v5 splits data into layers per sample, and some operations require unified layers. Run JoinLayers() on the RNA assay after integration before differential expression or other operations that need a single combined matrix.

Can I copy analysis results directly to /srv/gstore/projects/?

No, gStore is read-only via direct filesystem access. You must use g-req commands such as g-req copynow to transfer files, specifying only the parent directory to avoid creating nested directories.

How do I create interactive DEG tables in R Markdown reports?

Use ezInteractiveTableRmd() from the ezRun package instead of knitr::kable(). It produces sortable, searchable tables with pagination and CSV/Excel download buttons, which is the mandatory FGCZ standard for DEG and marker tables.