Luis E Solano, Negin Rahimzadeh, Zechuan Shi, Vivek Swarup
UNLABELLED: Single-nucleus multiome assays jointly profile gene expression and chromatin accessibility, yet their analysis typically requires bespoke chaining of modality-specific tools, creating barriers to reproducibility, scalability, and regulatory interpretation. We present FORGE (Flow Orchestrated Regulatory Genomics Engine), a configurable workflow that automates standalone snRNA-seq and snATAC-seq analysis, integrates the pair through complementary linear and nonlinear latent-variable models, and carries them through regulatory-network inference and differential testing. We evaluated FORGE on four human and mouse datasets spanning blood, brain, and kidney and two multiome chemistries, including a twelve-sample CRND8 Alzheimer's disease cohort. We report cross-modal agreement alongside missing-modality reconstruction and an accounting of computational cost. In the Alzheimer's cohort, FORGE nominated a glial Mef2c-associated program defensible across expression, co-accessibility, footprinting, and eRegulon evidence. In human PBMC, FORGE's layered evidence models also provide nuanced interpretations that largely corroborate previously published regulatory links while also proposing an additional CD83 myeloid module.
MOTIVATION: Single-nucleus resolution multiome (snMultiome) assays concurrently profile gene expression and chromatin accessibility in the same nucleus. Yet regulatory inference from such analyses are difficult to scale, audit, and reproduce; moreover, as a field, snMultiomics and its' toolset remains far from standardized. To address these challenges we developed FORGE, a configureable Nextflow workflow that carries paired data from raw counts and fragments through regulatory network inference with a comprehensive differential testing suite. Execution is containerized, tracks provenance, robust to interruption, optimized for cluster-based compute environments, and allows for nuanced customization of specific processes.