Population-scale skeletal muscle single-nucleus multi-omic profiling reveals extensive context specific genetic regulation
Arushi Varshney,Nandini Manickam,Peter Orchard,Adelaide Tovar,Christa Ventresca,Zhenhao Zhang,Fan Feng,Joseph Mears,Michael R Erdos,Narisu Narisu,Kirsten Nishino,Vivek Rai,Heather Stringham,Anne U Jackson,Tricia Tamsen,Chao Gao,Mao Yang,Olivia I Koues,Joshua D Welch,Charles F Burant,L. Keoki Williams,Chris Jenkinson,Ralph A DeFronzo,Luke Norton,Jouko Saramies,Timo A Lakka,Markku Laakso,Jaakko Tuomilehto,Karen L Mohlke,Jacob O Kitzman,Heikki A Koistinen,Jie Liu,Michael Boehnke,Francis S Collins,Laura J Scott,Stephen CJ Parker
DOI: https://doi.org/10.1101/2023.12.15.571696
2024-12-17
Abstract:Skeletal muscle, the largest human organ by weight, is relevant in several polygenic metabolic traits and diseases including type 2 diabetes (T2D). Identifying genetic mechanisms underlying these traits requires pinpointing cell types, regulatory elements, target genes, and causal variants. Here, we use genetic multiplexing to generate population-scale single nucleus (sn) chromatin accessibility (snATAC-seq) and transcriptome (snRNA-seq) maps across 287 frozen human skeletal muscle biopsies representing nearly half a million nuclei. We identify 13 cell types and integrate genetic variation to discover >7,000 expression quantitative trait loci (eQTL) and >100,000 chromatin accessibility QTLs (caQTL) across cell types. Learning patterns of e/caQTL sharing across cell types increased precision of effect estimates. We identify high-resolution cell-states and context-specific e/caQTL with significant genotype by context interaction. We identify nearly 2,000 eGenes colocalized with caQTL and construct causal directional maps for chromatin accessibility and gene expression. Almost 3,500 genome-wide association study (GWAS) signals across 38 relevant traits colocalize with sn-e/caQTL, most in a cell-specific manner. These signals typically colocalize with caQTL and not eQTL, highlighting the importance of population-scale chromatin profiling for GWAS functional studies. Finally, our GWAS-caQTL colocalization data reveal distinct cell-specific regulatory paradigms. Our results illuminate the genetic regulatory architecture of human skeletal muscle at high resolution epigenomic, transcriptomic, and cell-state scales and serve as a template for population-scale multi-omic mapping in complex tissues and traits.
Biology