Bulk RNA-seq differential expression with omicverse

Perform bulk RNA-seq differential expression analysis with omicverse.

1.2k|145|Updated Mar 22, 2021
One-click install
npx skills add https://github.com/Starlitnightly/omicverse --skill bulk-rna-seq-differential-expression-with-omicverse
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: Bulk RNA-seq differential expression with omicverse
Source: https://github.com/Starlitnightly/omicverse/tree/main/.claude/skills/bulk-deg-analysis
Command: npx skills add https://github.com/Starlitnightly/omicverse --skill bulk-rna-seq-differential-expression-with-omicverse

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires omicverse, scanpy, matplotlib, pandas, gseapy, and includes references (resource) components.

What problem does it solves? Analyzing bulk RNA-seq data for differentially expressed genes (DEGs) involves multiple steps from raw counts to statistical testing and pathway enrichment, which can be complex and time-consuming. This Skill automates the entire pipeline, providing a streamlined workflow for robust DEG analysis.

Core Features & Use Cases

  • Gene ID Mapping: Convert gene identifiers (e.g., Ensembl IDs) to gene symbols for consistent analysis.
  • Differential Expression Testing: Perform DESeq2-style normalization and statistical testing (t-test, edgeR, limma).
  • Visualization & Filtering: Generate volcano plots and boxplots, with dynamic fold-change and significance thresholding.
  • Pathway Enrichment: Conduct comprehensive pathway enrichment analysis on significant DEGs to uncover biological insights.
  • Use Case: Analyze a bulk RNA-seq dataset comparing treated vs. control samples, identify significant DEGs, visualize them with volcano plots, and then perform GO pathway enrichment to understand affected biological processes.

Quick Start

Run differential expression analysis on my bulk RNA-seq counts using a t-test, then generate a volcano plot and perform WikiPathways enrichment.

Frequently Asked Questions about Bulk RNA-seq differential expression with omicverse

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

FAQPage Schema
How do I perform differential expression analysis on bulk RNA-seq data?

Differential expression analysis identifies genes with significant expression changes between sample groups. This Skill automates the complete workflow: normalize raw counts, apply statistical testing (t-test, edgeR, limma), filter results by fold-change and p-value thresholds, and generate volcano plots—all within omicverse's integrated pipeline across multiple reference genomes.

Can I map gene identifiers and perform DESeq2-style normalization in a single workflow?

Yes. This Skill converts gene identifiers (e.g., Ensembl IDs to gene symbols) and performs DESeq2-style normalization and statistical testing in one streamlined process. You provide raw counts, sample labels, and treatment/control groupings; the workflow handles mapping, normalization, and hypothesis testing automatically.

What visualization and filtering options are available for DEG results?

The Skill generates volcano plots and boxplots with dynamic fold-change and significance thresholds, allowing you to filter and visualize differentially expressed genes interactively. Visualization helps identify and communicate significant DEGs alongside statistical results.

Does this support pathway enrichment analysis on DEGs?

Yes. After identifying significant DEGs, the Skill performs comprehensive pathway enrichment analysis using gseapy, supporting WikiPathways and other pathway libraries. This reveals affected biological processes and functional categories beyond individual gene results.

What input data and dependencies do I need to get started?

You need a Python environment with omicverse, scanpy, matplotlib, pandas, and gseapy installed; a raw gene-count matrix; gene ID annotation pairs in genesets/; and sample metadata with treatment versus control labels. Optional: pathway libraries for enrichment analysis.

Can I use this Skill for comparing treated versus control samples across multiple reference genomes?

Yes. The Skill is applicable to users with raw gene-count matrices and supports gene ID mapping across multiple reference genomes. Provide your counts, sample labels, and treatment/control groupings; the workflow produces DEG results tables, plots, and enrichment summaries for each genome context.