pydeseq2

Analyze bulk RNA-seq count matrices with DESeq2 workflows and FDR control.

Updated Jul 1, 2026
One-click install
npx skills add https://github.com/jasrajtulsi/GRAD-SCOPE --skill pydeseq2-jasrajtulsi
Or copy as Structured Prompt for Agent▼
Please help me install this Agent Skill.
Skill: pydeseq2
Source: https://github.com/jasrajtulsi/GRAD-SCOPE/tree/main/.claude/skills/pydeseq2
Command: npx skills add https://github.com/jasrajtulsi/GRAD-SCOPE --skill pydeseq2-jasrajtulsi

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires pydeseq2, pandas, numpy, matplotlib, and includes scripts (resource) and references (resource) components.

What problem does it solve?

This Skill automates bulk RNA-seq differential expression analysis so you can move from raw count tables to statistically tested gene-level results without assembling the workflow by hand.

Core Features & Use Cases

  • End-to-End DESeq2 Workflow: Fit size factors, dispersions, log fold changes, and Wald tests with Benjamini-Hochberg FDR correction.
  • Flexible Experimental Designs: Handle single-factor comparisons, batch-adjusted models, continuous covariates, and interaction terms using formulaic design strings.
  • Result Refinement and Reporting: Apply optional LFC shrinkage, export CSV and H5AD outputs, and generate volcano and MA plots for interpretation.
  • Use Cases: Compare treated versus control samples, control for batch effects, and prioritize significant genes for downstream biology or validation.

Quick Start

Use the pydeseq2 skill to analyze the attached count matrix with the provided metadata and a treated-versus-control contrast, then return significant genes and saved result files.

Frequently Asked Questions about pydeseq2

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

FAQPage Schema
How do I identify differentially expressed genes from bulk RNA-seq count data?▼

To identify differentially expressed genes from bulk RNA-seq count data, you need an integer samples × genes count matrix and sample metadata. The analysis fits size factors, dispersions, and applies Wald tests with Benjamini-Hochberg FDR correction to quantify condition effects.

Can I control for batch effects in RNA-seq differential expression analysis?▼

Yes, you can control for batch effects in RNA-seq differential expression analysis by using flexible formulaic design strings. This allows you to handle batch-adjusted models alongside single-factor comparisons, continuous covariates, and interaction terms across samples.

How do I perform LFC shrinkage and generate volcano plots for RNA-seq results?▼

You perform LFC shrinkage and generate volcano plots for RNA-seq results during the result refinement stage. After fitting log fold changes and running Wald tests, apply optional LFC shrinkage to prioritize significant genes and export visual MA and volcano plots for interpretation.

Does pydeseq2 support multi-factor experimental designs with interaction terms?▼

Yes, pydeseq2 supports multi-factor experimental designs with interaction terms. It handles treatment-versus-control comparisons, continuous covariates, and batch correction using explicit formulaic design strings and contrasts to quantify specific condition effects.

What is the best way to export statistically tested gene-level results from RNA-seq workflows?▼

The best way to export statistically tested gene-level results from RNA-seq workflows is to output CSV and H5AD files after running the analysis. This refines the data by applying Benjamini-Hochberg FDR control and optional LFC shrinkage for downstream biology validation.