Multi-scale gene regulatory networks for resolving transcription factor and autoimmune variant impacts on human B cell fate dynamics

Bifurcation of activated human B cells into plasmablast (PB) and germinal center precursor (preGC) fates underlies protective and autoimmune antibody responses, yet gene regulatory networks (GRNs) governing the alternative trajectories remain poorly defined. Using temporal single-cell multiomics, we assembled state-specific human B cell GRNs spanning four scales: transcription factor (TF)-to-fate, TF-to-gene, cis-regulatory element (CRE)-to-gene and nucleotide-to-CRE. Applying the framework to in vitro differentiated and tonsil B cells revealed concordant regulatory architectures and impacts of in silico TF perturbations. CRISPR perturbations validated many TF-to-fate and TF-to-gene predictions and revealed a mutually repressive BATF-IRF4/PRDM1 network module. The GRNs were used to predict and interpret effects of autoimmune disease variants, uncovering partitioning of disease risk at distinct TF motifs and B cell states. Predicted functional variants were independently supported by chromatin and expression QTLs and reporter assays. A web application enables exploration of B cell regulatory determinants of vaccine responses and autoimmune diseases.