bio-workflows-crispr-screen-pipeline

Analyzes pooled and single-cell CRISPR screens from FASTQ to consensus hit genes.

Updated Sep 15, 2026
One-click install
npx skills add https://github.com/mrsonord2240/openscience-specialists --skill bio-workflows-crispr-screen-pipeline-mrsonord2240
Or copy as Structured Prompt for Agent▼
Please help me install this Agent Skill.
Skill: bio-workflows-crispr-screen-pipeline
Source: https://github.com/mrsonord2240/openscience-specialists/tree/main/specialists/crispr-screen-analyst/versions/1.0.0/package/skills/bio-workflows-crispr-screen-pipeline
Command: npx skills add https://github.com/mrsonord2240/openscience-specialists --skill bio-workflows-crispr-screen-pipeline-mrsonord2240

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires pandas, numpy, matplotlib, pertpy, scanpy, anndata, scipy, seaborn, scikit-learn, biopython, statsmodels.

What problem does it solve? Pooled CRISPR screen analysis involves many fragile decisions—baseline choice, control classes, copy-number artifacts, and method selection—where a wrong-but-silent commitment invalidates every downstream hit with no error thrown. This Skill orchestrates the entire pipeline so each step uses the correct reference, QC gate, and hit-calling method for the experimental design. ## Core Features & Use Cases - End-to-end orchestration: Runs guide counting with MAGeCK, six-stage QC (plasmid Gini, replicate Pearson, CEGv2 PR-AUC), copy-number correction with CRISPRcleanR or Chronos, and batch covariates via MAGeCK MLE. - Design-matched hit calling: Selects among MAGeCK RRA/MLE, BAGEL2, drugZ, JACKS, and Chronos based on screen type, then builds tier-based consensus (Tier 1 = 3-method agreement) across methods. - Specialized branches: Hands off to dedicated workflows for Perturb-seq, combinatorial paralog screens, base-editor and prime-editor variant screens, and in vivo bottleneck-aware screens. - Use Case: Given FASTQ files from a Brunello dropout screen with vehicle and drug arms, the Skill counts guides, verifies QC gates, runs drugZ and MAGeCK in parallel, and outputs tier-2 consensus hits at FDR < 0.05. ## Quick Start Ask the agent to run the end-to-end CRISPR screen pipeline on your pooled screen FASTQ files, specifying the baseline, library control classes, and screen type so it can count guides, run QC, and call consensus hits.

Frequently Asked Questions about bio-workflows-crispr-screen-pipeline

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

FAQPage Schema
How do I analyze a pooled CRISPR screen from FASTQ to hit genes?▼

Run mageck count with your library CSV and sample labels to build a count matrix, pass six QC gates (plasmid Gini <0.1, replicate Pearson >=0.8, CEGv2 PR-AUC >0.7), then call hits with the design-matched method such as MAGeCK RRA, BAGEL2, or drugZ, and consolidate with tier-based consensus.

MAGeCK RRA vs BAGEL2 vs drugZ: which hit-calling method should I use?▼

Match the method to the design: MAGeCK RRA for two-condition essentiality, MAGeCK MLE for time course or multi-condition, drugZ for drug-versus-vehicle screens, JACKS for multi-screen joint analysis, and Chronos for cancer cell-line panels. Running two methods enables tier-based consensus.

Why do amplified regions like ERBB2 and MYC appear as false essential genes?▼

Multiple simultaneous Cas9 cuts at amplified loci trigger gene-independent DNA-damage arrest, so amplified regions look essential regardless of gene function. Apply CRISPRcleanR or Chronos copy-number correction before hit calling and verify abs(rho(LFC,CN)) < 0.1 afterward.

Should I use vehicle or Day 0 as the control for a drug screen?▼

Use vehicle as the control for drug screens, not Day 0. Comparing drug-treated samples to Day 0 conflates the drug effect with normal proliferation, rescaling every downstream hit.

Why are my BAGEL2 results different between identical reruns?▼

BAGEL2's bf step uses a clock-based default seed, so unseeded reruns produce different Bayes factors and can flip hit calls. Always pass a fixed -s seed to make results reproducible.

What does it mean when CEGv2 PR-AUC falls below 0.7?▼

A CEGv2 PR-AUC below 0.7 means known essential genes did not deplete, indicating the screen failed selection. No novel hit from such a screen is trustworthy regardless of its p-value; diagnose Cas9 activity, timepoint, or library issues before proceeding.