rosetta-bioc 0.2.1

Pythonic API for R/Bioconductor statistical methods — calls validated R code, returns pandas DataFrames.

🪨 rosetta

Python interface to R/Bioconductor — pandas in, pandas out, .report() when you're done.

pip install rosetta-bioc

30-second demo

import rosetta as rb

# DESeq2 differential expression — one call, pandas out
results = rb.deseq2(counts_df, metadata_df, design="~ condition")
results.report()

DESeq2 Results Summary
──────────────────────────────
Total genes tested:      12,000
Significant (padj<0.05): 843 (7.0%)
  ↑ Upregulated:         428
  ↓ Downregulated:       415
LFC range:               [-4.71, 3.50]

That's it. No R code. No rpy2 boilerplate. No type conversion. Just results.

What it wraps

R Package	Python	What it does
DESeq2	rb.deseq2()	Differential expression (negative binomial)
edgeR	rb.edger()	Quasi-likelihood differential expression
limma	rb.limma_voom()	Linear models + TREAT significance
clusterProfiler	rb.enrich_go()	GO/KEGG/Reactome pathway enrichment
phyloseq	rb.phyloseq()	Microbiome diversity analysis
Seurat	rb.seurat()	Single-cell RNA-seq

All functions return a RosettaDataFrame (pandas DataFrame subclass) with a .report() method.

Not a toy — full design support

• Multi-factor designs: design="~ batch + condition", interaction terms, blocking factors
• LFC thresholds: proper hypothesis testing via lfcThreshold (not post-hoc filtering)
• Shrinkage: apeglm, ashr, normal — via lfc_shrink()
• Contrasts: contrast=["genotype", "mutant", "wildtype"]
• QC/normalization/outliers: DESeq2's size factors, Cook's distance, independent filtering all run normally — Rosetta doesn't hide the fitted object
• Weights, correlations: limma-voom with duplicateCorrelation, sample weights — everything the R function accepts, Rosetta passes through

Show me the R code

Don't trust a black box? Turn on codegen to see exactly what's running:

import rosetta as rb
rb.codegen.enable()

dds = rb.wrappers.deseq2.run_deseq2(counts, meta, design="~ batch + condition")
res = rb.wrappers.deseq2.get_results(dds, lfc_threshold=1.0)

  R> library(DESeq2)
  R> dds <- DESeqDataSetFromMatrix(countData=counts, colData=metadata, design=~ batch + condition)
  R> dds <- DESeq(dds)
  R> res <- results(dds, alpha=0.1, lfcThreshold=1.0)

rb.codegen.last() returns it as a string — paste into R to reproduce independently.

Modular DESeq2 API

For more control, use the step-by-step interface:

from rosetta.wrappers.deseq2 import run_deseq2, get_results, lfc_shrink

dds = run_deseq2(counts_df, metadata_df, design="~ condition")
res = get_results(dds, contrast=["condition", "treated", "control"], alpha=0.05)
shrunk = lfc_shrink(dds, coef="condition_treated_vs_control", type="apeglm")

res.report()
shrunk.report()

Enrichment analysis

import rosetta as rb

# Over-representation analysis
go_results = rb.enrich_go(gene_list, org_db="org.Hs.eg.db", ont="BP")
go_results.report()

# KEGG pathways
kegg = rb.enrich_kegg(gene_list, organism="hsa")
kegg.report()

Setup

Python side:

pip install rosetta-bioc

R side (one-time):

Rscript install.R

Or manually:

BiocManager::install(c("DESeq2", "edgeR", "limma", "clusterProfiler"))

Posit Cloud: See docs/posit-cloud.md for zero-config setup.

Requirements

• Python 3.9+
• R 4.0+ with Bioconductor
• rpy2 ≥ 3.5

Philosophy

1. Rosetta calls R — it doesn't reimplement it. All statistics run in the original, validated R packages.
2. Pandas in, pandas out. No R objects leak into your Python workflow.
3. Fail early, fail clearly. Input validation happens in Python before crossing the R boundary.
4. .report() everything. Results should be immediately interpretable without manual inspection.
5. Show your work. codegen prints the equivalent R code so you can verify, reproduce, or learn.

Contributing

See CONTRIBUTING.md. Good first issues are labeled — start with Issue #1: report() enhancements.

Acknowledgments

Built on rpy2 and the extraordinary R/Bioconductor ecosystem. All credit for the statistical methods goes to the original R package authors.

Supported by:

• Google Summer of Code 2026 — funding Catherine's development work
• JPMorgan Chase — startup banking and advisory through their Innovation Economy program
• AWS — quantum computing infrastructure via Amazon Braket
• Nodes Bio, Inc. — project lead, CI/hosting, and engineering

GSoC 2026 · MIT License · Nodes Bio