scellrun 0.6.0

Opinionated, report-first CLI for single-cell and multi-omics analysis. Sane defaults baked in, defensible deliverables out.

scellrun

scellrun lowers the bar to running a defensible single-cell analysis. One command, a publication-quality report. You don't need to learn scanpy first.

# 1. one-time: a clean Python environment so scellrun's deps don't collide with anything else on your machine
conda create -n scellrun python=3.11 -y
conda activate scellrun

# 2. install
pip install scellrun

# 3. run the full pipeline in one shot — qc → integrate → markers → annotate → report
scellrun analyze data.h5ad --tissue "OA cartilage"
# → opens with a single index.html link you can drop into a browser

Don't have a .h5ad? Cellranger output works directly:

scellrun scrna convert path/to/cellranger_outs -o data.h5ad
scellrun analyze data.h5ad --tissue "OA cartilage"

Want a Chinese report? Add --lang zh.

Status: v0.1.0, early alpha. APIs will change without warning. Don't pin against a tag yet.

Who this is for

• A clinician or rotating student looking at scRNA-seq data for the first time
• A postdoc on a deadline who doesn't want to write QC boilerplate
• A bioinformatics core that needs every project to look the same in a report
• An LLM coding agent that should reach for a real tool instead of re-deriving thresholds

If you already have your own pipeline, you don't need scellrun. If you don't, or you're tired of re-litigating the same decisions every project, this is for you.

Why this exists

Tools like scanpy, pyDESeq2, decoupler, SCENIC and nf-core are excellent and intentionally unopinionated. Every new project re-litigates the same decisions:

• Which mt% threshold? Doublet filter before or after batch correction?
• The composite score correlates with disease — Spearman ρ on each of 18 clinical features, or one well-chosen dichotomy?
• What does the QC report actually need to look like for a reviewer to stop asking?

scellrun answers these once, in code, by encoding the working practice of a clinician + bioinformatics team. Defaults aren't "neutral" — they're a real position, made for real reasons, with the rationale rendered into every report.

You can override anything. But if you don't override, you get a defensible analysis on day one.

Where this fits

Tool	What it is
Anthropic scientific skills	Prompt scaffolds that teach an LLM how to call scanpy etc.
Bioconda recipes	Packaging
nf-core pipelines	General-purpose, infrastructure-heavy, vendor-neutral workflows
scellrun	An opinionated, report-first CLI optimized for low barrier to entry

v0.1 scope

• scellrun scrna qc <h5ad> — single-cell QC (mt%, ribo%, hb%, n_genes, doublets) with defensible thresholds. Outputs an HTML report with rationale text, a per-cell CSV, and an annotated qc.h5ad ready for the v0.2 integration step. Cells are flagged (obs.scellrun_qc_pass), never silently dropped.

That's it. See ROADMAP.md for v0.2-v0.5 (integrate, markers, annotate, multi-stage report) and the planned scellrun run --steps ... pipeline mode.

Install

For users — clean conda env (recommended)

conda create -n scellrun python=3.11 -y
conda activate scellrun
pip install scellrun

If you don't have conda, miniconda takes about 2 minutes to install. Or use uv for a faster venv:

uv venv .scellrun-env && source .scellrun-env/bin/activate
uv pip install scellrun

Either way, scellrun ends up in its own environment — your other Python projects (an old scanpy, Seurat-via-rpy2, etc.) won't be touched.

For contributors — editable from source

git clone https://github.com/drstrangerujn/scellrun.git
cd scellrun
conda create -n scellrun-dev python=3.11 -y && conda activate scellrun-dev
pip install -e ".[dev]"
pytest -q          # 25 tests should pass

Profiles

Different tissue domains have different defaults. v0.1 ships two profiles:

• default — fresh-tissue 10x v3 chemistry, joint-tissue-aware mt% ceiling
• joint-disease — tighter hb% for avascular cartilage; ships the Fan 2024 chondrocyte 11-subtype panel + 15-group broad celltype panel for v0.4 annotation

scellrun profiles list
scellrun profiles show joint-disease   # prints thresholds + panels

Adding a profile = one Python file under src/scellrun/profiles/. If your community has working practice, contribute a profile.

For LLM agents

skills/scellrun/SKILL.md is a portable instruction document that teaches Claude Code, Hermes, Codex, or any markdown-aware agent harness when and how to invoke scellrun. Install it by symlinking into your agent's skills directory — see skills/README.md for one-liner install commands per agent.

License

MIT — see LICENSE.

Acknowledgements

Defaults reflect the working practice of a clinician + bioinformatics team that has shipped these analyses for the OARSI / MSK community. The R "AIO" pipeline that prefigured scellrun's design is documented in ROADMAP.md.

Built with assistance from Claude (Anthropic).