arboreto-grn-inference

npx skills add BioTender-max/awesome-bio-agent-skills --skill arboreto-grn-inference --agent claude-code

Same command for any agent — swap --agent for codex, cursor, copilot.

# Arboreto GRN Inference ## Overview Arboreto infers gene regulatory networks (GRNs) from gene expression data using parallelized tree-based regression. For each target gene, it trains a regression model with all other genes (or a specified TF list) as features and emits TF-target-importance triplets. It provides two interchangeable algorithms -- GRNBoost2 (gradient boosting, fast) and GENIE3 (Random Forest, classic) -- sharing identical input/output formats. Computation is Dask-parallelized, scaling from laptop cores to HPC clusters. ## When to Use - Inferring transcription factor-to-target gene regulatory relationships from bulk RNA-seq expression data - Building gene regulatory networks from single-cell RNA-seq count matrices (cells as rows, genes as columns) - Generating the adjacency matrix (Step 1) of the pySCENIC regulatory analysis pipeline - Comparing regulatory network structure across experimental conditions (e.g., control vs treatment) - Producing consensus regulatory networks by running inference across multiple random seeds - Validating GRN results by comparing GRNBoost2 and GENIE3 outputs on the same dataset - For downstream regulon identification and activity scorin

pip install arboreto distributed networkx matplotlib
On head node: start scheduler
dask-scheduler
On each compute node: start workers
dask-worker tcp://10.118.224.134:8786 --nprocs 4 --nthreads 1 --memory-limit 16GB