rnaseq-analysis

Identify differential expression in bulk RNA-seq count matrices using DESeq2.

1|1|Updated Feb 7, 2026
One-click install
npx skills add https://github.com/gexijin/vitiligo --skill rnaseq-analysis
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: rnaseq-analysis
Source: https://github.com/gexijin/vitiligo/tree/main/.claude/skills/rnaseq-analysis
Command: npx skills add https://github.com/gexijin/vitiligo --skill rnaseq-analysis

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires DESeq2, ggplot2, dplyr, tidyr, pheatmap, RColorBrewer, here, GEOquery, AnnotationDbi, org.Hs.eg.db, clusterProfiler, biomaRt, apeglm.

What problem does it solve?

Analyze bulk RNA-seq count data to perform QC, normalization, differential expression testing, and gene ID mapping using DESeq2, enabling reproducible end-to-end analysis workflows.

Core Features & Use Cases

  • QC, normalization, and DESeq2 modeling for bulk RNA-seq count matrices from various data formats (raw counts, featureCounts, RSEM, GEO series).
  • Differential expression results with shrunken LFC (via apeglm/normal method for paired designs) and annotated gene IDs for downstream interpretation.
  • Flexible analysis for unpaired vs paired designs, multi-factor models, and downstream visualizations (PCA, heatmaps, GO/KEGG enrichment) with a portable, here-based workflow.

Quick Start

Analyze the provided RNA-seq count matrix by performing QC, normalization, DESeq2 modeling, and annotating results.

Frequently Asked Questions about rnaseq-analysis

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

FAQPage Schema
How do I run differential expression analysis on RNA-seq count data using DESeq2?

Differential expression analysis on RNA-seq count data using DESeq2 involves applying the tool to a count matrix to perform QC, normalization, and modeling, yielding results with shrunken LFC and annotated gene IDs for downstream interpretation.

Can I analyze paired design bulk RNA-seq data with batch effects?

Yes, paired design bulk RNA-seq data with batch effects can be analyzed using DESeq2 multi-factor models. The workflow supports both unpaired and paired designs, outputting DE results with shrunken LFC via the apeglm or normal method.

What input formats are supported for RNA-seq differential expression?

RNA-seq differential expression supports count matrices from raw counts, featureCounts output, RSEM counts, or GEO-derived series like GSE####, requiring proper data formatting and the here package for portable paths.

How do I perform GO and KEGG enrichment after RNA-seq differential expression?

GO and KEGG enrichment after RNA-seq differential expression is performed using clusterProfiler on annotated gene IDs. The workflow outputs DE results with shrunken LFC, annotated genes, PCA plots, and heatmaps for downstream interpretation.

What is the best way to map gene IDs for RNA-seq DESeq2 results?

Mapping gene IDs for RNA-seq DESeq2 results is handled via AnnotationDbi, org.Hs.eg.db, and biomaRt. The workflow annotates differential expression results, providing annotated gene IDs required for downstream visualizations and enrichment analysis.

Does DESeq2 work with GEO series for bulk RNA-seq normalization?

Yes, DESeq2 works with GEO series for bulk RNA-seq normalization by processing GEO-derived series like GSE####. The workflow applies QC, normalization, and DESeq2 modeling to generate DE results with shrunken LFC and plots.