Abstract
Understanding gene regulatory networks (GRNs) is essential for deciphering biological processes and disease mechanisms. Single-cell multiome technologies now enable joint profiling of chromatin accessibility and gene expression, offering an powerful means to infer cell type-specific GRNs. However, existing methods analyze each cell type independently or aggregate data into pseudo-bulk profiles, limiting their ability to resolve rare populations and capture cellular heterogeneity. We introduce BayesCNet, a Bayesian hierarchical model that jointly infers enhancer-gene linkages across all cell types while leveraging their hierarchical relationships for information sharing. Through extensive simulations, BayesCNet consistently outperforms state-of-the-art methods, with the largest improvements in rare cell types. When applied to real datasets, BayesCNet identifies enhancer-gene linkages with higher accuracy validated by promoter-capture Hi-C data, and reconstructs cell type-specific GRNs that highlight key regulators, demonstrating its power to resolve gene regulatory programs across diverse cell types.