Cellular heterogeneity within the sinoatrial node (SAN) is functionally important but has been difficult to model in vitro , presenting a major obstacle to studies of heart rate regulation and arrhythmias. Here we describe a scalable method to derive sinoatrial node pacemaker cardiomyocytes (PCs) from human induced pluripotent stem cells that recapitulates differentiation into distinct PC subtypes, including SAN Head, SAN Tail, transitional zone cells, and sinus venosus myocardium. Single cell (sc) RNA-sequencing, sc-ATAC-sequencing, and trajectory analyses were used to define epigenetic and transcriptomic signatures of each cell type, and to identify novel transcriptional pathways important for PC subtype differentiation. Integration of our multi-omics datasets with genome wide association studies uncovered cell type-specific regulatory elements that associated with heart rate regulation and susceptibility to atrial fibrillation. Taken together, these datasets validate a novel, robust, and realistic in vitro platform that will enable deeper mechanistic exploration of human cardiac automaticity and arrhythmia.

Single Cell Multi-Omics of an iPSC Model of Human Sinoatrial Node Development Reveals Genetic Determinants of Heart Rate and Arrhythmia Susceptibility