Building rustscenic, a Rust + PyO3 rewrite of SCENIC+, with the Kuan-lin Huang Lab at Icahn Mount Sinai.
rustscenic (v0.4.0 release): full SCENIC+ pipeline for single-cell regulatory network analysis. GRN, AUCell, topics, cisTarget, enhancer links, eRegulons.
Benchmarks vs the reference Python stack on identical input (v0.3.x measurements; v0.4.x cross-dataset sweep underway):
- AUCell 88× faster than pyscenic at per-cell Pearson 0.99 (0.21 s vs 18.6 s on 10x Multiome, 10k cells × 1,457 regulons)
- GRN 1.78× faster than pinned arboreto 0.1.6 on PBMC (per-TF Spearman 0.63)
- ~6.3× less memory at 100k cells × 20k genes × 4 stages (6.3 GB vs >40 GB reported for scenicplus)
- 5 runtime dependencies (numpy, pandas, pyarrow, scipy, anndata) vs 40+ for the reference stack
- 9/9 cortex TFs recovered on mouse brain E18 multiome (biological validation)
- 7 merged scverse PRs across scanpy and PyDESeq2; open AnnData concat API PR
- 1,773 DE genes (Zenodo DOI): DESeq2 on SARS-CoV-2 nasopharyngeal RNA-seq, n = 60 primary, 99.8% concordant with n = 484 sensitivity (repo)
- r = 0.954 on 5-fold CV: PyTorch deconvolution of 484 bulk RNA-seq samples into 14 airway cell types, upper bound with no batch effects (repo)
- Live deploy: SafetyNett, AI safety-netting prototype for NHS GP workflows, OpenClaw Clinical Hackathon (repo)



