Genome-Wide Analysis Identifies an Essential Human TBX3 Pacemaker Enhancer

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Transcriptome and open chromatin analysis reveals the process of myocardial cell development and key pathogenic target proteins in Long QT syndrome type 7

OBJECTIVE

Long QT syndrome type 7 (Andersen-Tawil syndrome, ATS), which is caused by KCNJ2 gene mutation, often leads to ventricular arrhythmia, periodic paralysis and skeletal malformations. The development, differentiation and electrophysiological maturation of cardiomyocytes (CMs) changes promote the pathophysiology of Long QT syndrome type 7(LQT7). We aimed to specifically reproduce the ATS disease phenotype and study the pathogenic mechanism.

METHODS AND RESULTS

We established a cardiac cell model derived from human induced pluripotent stem cells (hiPSCs) to the phenotypes and electrophysiological function, and the establishment of a human myocardial cell model that specifically reproduces the symptoms of ATS provides a reliable platform for exploring the mechanism of this disease or potential drugs. The spontaneous pulsation rate of myocardial cells in the mutation group was significantly lower than that in the repair CRISPR group, the action potential duration was prolonged, and the Kir2.1 current of the inward rectifier potassium ion channel was decreased, which is consistent with the clinical symptoms of ATS patients. Only ZNF528, a chromatin-accessible TF related to pathogenicity, was continuously regulated beginning from the cardiac mesodermal precursor cell stage (day 4), and continued to be expressed at low levels, which was identified by WGCNA method and verified with ATAC-seq data in the mutation group. Subsequently, it indicated that seven pathways were downregulated (all p < 0.05) by used single sample Gene Set Enrichment Analysis to evaluate the overall regulation of potassium-related pathways enriched in the transcriptome and proteome of late mature CMs. Among them, the three pathways (GO: 0008076, GO: 1990573 and GO: 0030007) containing the mutated gene KCNJ2 is involved that are related to the whole process by which a potassium ion enters the cell via the inward rectifier potassium channel to exert its effect were inhibited. The other four pathways are related to regulation of the potassium transmembrane pathway and sodium:potassium exchange ATPase (p < 0.05). ZNF528 small interfering (si)-RNA was applied to hiPSC-derived cardiomyocytes for CRISPR group to explore changes in potassium ion currents and growth and development related target protein levels that affect disease phenotype. Three consistently downregulated proteins (KCNJ2, CTTN and ATP1B1) associated with pathogenicity were verificated through correlation and intersection analysis.

CONCLUSION

This study uncovers TFs and target proteins related to electrophysiology and developmental pathogenicity in ATS myocardial cells, obtaining novel targets for potential therapeutic candidate development that does not rely on gene editing.

Development of the Cardiac Conduction System

The electrical impulses that coordinate the sequential, rhythmic contractions of the atria and ventricles are initiated and tightly regulated by the specialized tissues of the cardiac conduction system. In the mature heart, these impulses are generated by the pacemaker cardiomyocytes of the sinoatrial node, propagated through the atria to the atrioventricular node where they are delayed and then rapidly propagated to the atrioventricular bundle, right and left bundle branches, and finally, the peripheral ventricular conduction system. Each of these specialized components arise by complex patterning events during embryonic development. This chapter addresses the origins and transcriptional networks and signaling pathways that drive the development and maintain the function of the cardiac conduction system.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Novel Candidate Regulators and Developmental Trajectory of Pituitary Thyrotropes

The pituitary gland regulates growth, metabolism, reproduction, the stress response, uterine contractions, lactation, and water retention. It secretes hormones in response to hypothalamic input, end organ feedback, and diurnal cues. The mechanisms by which pituitary stem cells are recruited to proliferate, maintain quiescence, or differentiate into specific cell types, especially thyrotropes, are not well understood. We used single-cell RNA sequencing in juvenile P7 mouse pituitary cells to identify novel factors in pituitary cell populations, with a focus on thyrotropes and rare subtypes. We first observed cells coexpressing markers of both thyrotropes and gonadotropes, such as Pou1f1 and Nr5a1. This was validated in vivo by both immunohistochemistry and lineage tracing of thyrotropes derived from Nr5a1-Cre; mTmG mice and demonstrates that Nr5a1-progenitors give rise to a proportion of thyrotropes during development. Our data set also identifies novel factors expressed in pars distalis and pars tuberalis thyrotropes, including the Shox2b isoform in all thyrotropes and Sox14 specifically in Pou1f1-negative pars tuberalis thyrotropes. We have therefore used single-cell transcriptomics to determine a novel developmental trajectory for thyrotropes and potential novel regulators of thyrotrope populations.

Coronary Artery Disease Risk Gene PRDM16 is Preferentially Expressed in Vascular Smooth Muscle Cells and a Potential Novel Regulator of Smooth Muscle Homeostasis

Objective

Vascular smooth muscle cells (VSMCs) are the primary contractile component of blood vessels and can undergo phenotypic switching from a contractile to a synthetic phenotype in vascular diseases such as coronary artery disease (CAD). This process leads to decreased expression of SMC lineage genes and increased proliferative, migratory and secretory abilities that drive disease progression. Super-enhancers (SE) and occupied transcription factors are believed to drive expression of genes that maintain cell identify and homeostasis. The goal of this study is to identify novel regulator of VSMC homeostasis by screening for SE-regulated transcription factors in arterial tissues.

Approach and Results

We characterized human artery SEs by analyzing the enhancer histone mark H3K27ac ChIP-seq data of multiple arterial tissues. We unexpectedly discovered the transcription factor PRDM16, a GWAS identified CAD risk gene with previously well-documented roles in brown adipocytes but with an unknown function in vascular disease progression, is enriched with artery-specific SEs. Further analysis of public bulk RNA-seq and scRNA-seq datasets, as well as qRT-PCR and Western blotting analysis, demonstrated that PRDM16 is preferentially expressed in arterial tissues and in contractile VSMCs but not in visceral SMCs, and down-regulated in phenotypically modulated VSMCs. To explore the function of Prdm16 in vivo, we generated Prdm16 SMC-specific knockout mice and performed histological and bulk RNA-Seq analysis of aortic tissues. SMC-deficiency of Prdm16 does not affect the aortic morphology but significantly alters expression of many CAD risk genes and genes involved in VSMC phenotypic modulation. Specifically, Prdm16 negatively regulates the expression of Tgfb2 that encodes for an upstream ligand of TGF-β signaling pathway, potentially through binding to the promoter region of Tgfb2 . These transcriptomic changes likely disrupt VSMC homeostasis and predispose VSMCs to a disease state.

Conclusions

Our results suggest that the CAD risk gene PRDM16 is preferentially expressed in VSMCs and is a novel regulator of VSMC homeostasis. Future studies are warranted to investigate its role in VSMCs under pathological conditions such as atherosclerosis.

Global chromatin landscapes identify candidate noncoding modifiers of cardiac rhythm

Comprehensive cis-regulatory landscapes are essential for accurate enhancer prediction and disease variant mapping. Although cis-regulatory element (CRE) resources exist for most tissues and organs, many rare - yet functionally important - cell types remain overlooked. Despite representing only a small fraction of the heart's cellular biomass, the cardiac conduction system (CCS) unfailingly coordinates every life-sustaining heartbeat. To globally profile the mouse CCS cis-regulatory landscape, we genetically tagged CCS component-specific nuclei for comprehensive assay for transposase-accessible chromatin-sequencing (ATAC-Seq) analysis. Thus, we established a global CCS-enriched CRE database, referred to as CCS-ATAC, as a key resource for studying CCS-wide and component-specific regulatory functions. Using transcription factor (TF) motifs to construct CCS component-specific gene regulatory networks (GRNs), we identified and independently confirmed several specific TF sub-networks. Highlighting the functional importance of CCS-ATAC, we also validated numerous CCS-enriched enhancer elements and suggested gene targets based on CCS single-cell RNA-Seq data. Furthermore, we leveraged CCS-ATAC to improve annotation of existing human variants related to cardiac rhythm and nominated a potential enhancer-target pair that was dysregulated by a specific SNP. Collectively, our results established a CCS-regulatory compendium, identified novel CCS enhancer elements, and illuminated potential functional associations between human genomic variants and CCS component-specific CREs.

Hippo-Yap Signaling Maintains Sinoatrial Node Homeostasis

BACKGROUND

The sinoatrial node (SAN) functions as the pacemaker of the heart, initiating rhythmic heartbeats. Despite its importance, the SAN is one of the most poorly understood cardiac entities because of its small size and complex composition and function. The Hippo signaling pathway is a molecular signaling pathway fundamental to heart development and regeneration. Although abnormalities of the Hippo pathway are associated with cardiac arrhythmias in human patients, the role of this pathway in the SAN is unknown.

METHODS

We investigated key regulators of the Hippo pathway in SAN pacemaker cells by conditionally inactivating the Hippo signaling kinases Lats1 and Lats2 using the tamoxifen-inducible, cardiac conduction system-specific Cre driver Hcn4CreERT2 with Lats1 and Lats2 conditional knockout alleles. In addition, the Hippo-signaling effectors Yap and Taz were conditionally inactivated in the SAN. To determine the function of Hippo signaling in the SAN and other cardiac conduction system components, we conducted a series of physiological and molecular experiments, including telemetry ECG recording, echocardiography, Masson Trichrome staining, calcium imaging, immunostaining, RNAscope, cleavage under targets and tagmentation sequencing using antibodies against Yap1 or H3K4me3, quantitative real-time polymerase chain reaction, and Western blotting. We also performed comprehensive bioinformatics analyses of various datasets.

RESULTS

We found that Lats1/2 inactivation caused severe sinus node dysfunction. Compared with the controls, Lats1/2 conditional knockout mutants exhibited dysregulated calcium handling and increased fibrosis in the SAN, indicating that Lats1/2 function through both cell-autonomous and non-cell-autonomous mechanisms. It is notable that the Lats1/2 conditional knockout phenotype was rescued by genetic deletion of Yap and Taz in the cardiac conduction system. These rescued mice had normal sinus rhythm and reduced fibrosis of the SAN, indicating that Lats1/2 function through Yap and Taz. Cleavage Under Targets and Tagmentation sequencing data showed that Yap potentially regulates genes critical for calcium homeostasis such as Ryr2 and genes encoding paracrine factors important in intercellular communication and fibrosis induction such as Tgfb1 and Tgfb3. Consistent with this, Lats1/2 conditional knockout mutants had decreased Ryr2 expression and increased Tgfb1 and Tgfb3 expression compared with control mice.

CONCLUSIONS

We reveal, for the first time to our knowledge, that the canonical Hippo-Yap pathway plays a pivotal role in maintaining SAN homeostasis.

Association of T-box gene polymorphisms with the risk of Wolff-Parkinson-White syndrome in a Han Chinese population

Abnormal development of the atrioventricular ring can lead to the formation of a bypass pathway and the occurrence of Wolff-Parkinson-White (WPW) syndrome. The genetic mechanism underlying the sporadic form of WPW syndrome remains unclear. Existing evidence suggests that both T-box transcription factor 3 (TBX3) and T-box transcription factor 2 (TBX2) genes participate in regulating annulus fibrosus formation and atrioventricular canal development. Thus, we aimed to examine whether single-nucleotide polymorphisms (SNPs) in the TBX3 and TBX2 genes confer susceptibility to WPW syndrome in a Han Chinese Population. We applied a SNaPshot SNP assay to analyze 5 selected tagSNPs of TBX3 and TBX2 in 230 patients with sporadic WPW syndrome and 231 sex- and age-matched controls. Haplotype analysis was performed using Haploview software. Allele C of TBX3 rs1061657 was associated with a higher risk of WPW syndrome (odds ratio [OR] = 1.41, 95% confidence interval [CI]: 1.08-1.83, P = .011) and left-sided accessory pathways (OR = 1.40, 95% CI: 1.07-1.84, P = .016). However, allele C of TBX3 rs8853 was likely to reduce these risks (OR = 0.71, 95% CI: 0.54-0.92, P = .011; OR = 0.70, 95% CI: 0.53-0.92, P = .011, respectively). The data revealed no association between TBX3 rs77412687, TBX3 rs2242442, or TBX2 rs75743672 and WPW syndrome. TBX3 rs1061657 and rs8853 are significantly associated with sporadic WPW syndrome among a Han Chinese population. To verify our results, larger sample sizes are required in future studies.

Regulation of cardiac ion channels by transcription factors: Looking for new opportunities of druggable targets for the treatment of arrhythmias

Cardiac electrical activity is governed by different ion channels that generate action potentials. Acquired or inherited abnormalities in the expression and/or function of ion channels usually result in electrophysiological changes that can cause cardiac arrhythmias. Transcription factors (TFs) control gene transcription by binding to specific DNA sequences adjacent to target genes. Linkage analysis, candidate-gene screening within families, and genome-wide association studies have linked rare and common genetic variants in the genes encoding TFs with genetically-determined cardiac arrhythmias. Besides its critical role in cardiac development, recent data demonstrated that they control cardiac electrical activity through the direct regulation of the expression and function of cardiac ion channels in adult hearts. This narrative review summarizes some studies showing functional data on regulation of the main human atrial and ventricular Na⁺, Ca²⁺, and K⁺ channels by cardiac TFs such as Pitx2c, Tbx20, Tbx5, Zfhx3, among others. The results have improved our understanding of the mechanisms regulating cardiac electrical activity and may open new avenues for therapeutic interventions in cardiac acquired or inherited arrhythmias through the identification of TFs as potential drug targets. Even though TFs have for a long time been considered as 'undruggable' targets, advances in structural biology have led to the identification of unique pockets in TFs amenable to be targeted with small-molecule drugs or peptides that are emerging as novel therapeutic drugs.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries

RATIONALE

Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored.

OBJECTIVE

We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA.

METHODS AND RESULTS

We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart.

CONCLUSIONS

This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

Genomic enhancers in cardiac development and disease

The Human Genome Project marked a major milestone in the scientific community as it unravelled the ~3 billion bases that are central to crucial aspects of human life. Despite this achievement, it only scratched the surface of understanding how each nucleotide matters, both individually and as part of a larger unit. Beyond the coding genome, which comprises only ~2% of the whole genome, scientists have realized that large portions of the genome, not known to code for any protein, were crucial for regulating the coding genes. These large portions of the genome comprise the 'non-coding genome'. The history of gene regulation mediated by proteins that bind to the regulatory non-coding genome dates back many decades to the 1960s. However, the original definition of 'enhancers' was first used in the early 1980s. In this Review, we summarize benchmark studies that have mapped the role of cardiac enhancers in disease and development. We highlight instances in which enhancer-localized genetic variants explain the missing link to cardiac pathogenesis. Finally, we inspire readers to consider the next phase of exploring enhancer-based gene therapy for cardiovascular disease.

Inherited and Acquired Rhythm Disturbances in Sick Sinus Syndrome, Brugada Syndrome, and Atrial Fibrillation: Lessons from Preclinical Modeling

Rhythm disturbances are life-threatening cardiovascular diseases, accounting for many deaths annually worldwide. Abnormal electrical activity might arise in a structurally normal heart in response to specific triggers or as a consequence of cardiac tissue alterations, in both cases with catastrophic consequences on heart global functioning. Preclinical modeling by recapitulating human pathophysiology of rhythm disturbances is fundamental to increase the comprehension of these diseases and propose effective strategies for their prevention, diagnosis, and clinical management. In silico, in vivo, and in vitro models found variable application to dissect many congenital and acquired rhythm disturbances. In the copious list of rhythm disturbances, diseases of the conduction system, as sick sinus syndrome, Brugada syndrome, and atrial fibrillation, have found extensive preclinical modeling. In addition, the electrical remodeling as a result of other cardiovascular diseases has also been investigated in models of hypertrophic cardiomyopathy, cardiac fibrosis, as well as arrhythmias induced by other non-cardiac pathologies, stress, and drug cardiotoxicity. This review aims to offer a critical overview on the effective ability of in silico bioinformatic tools, in vivo animal studies, in vitro models to provide insights on human heart rhythm pathophysiology in case of sick sinus syndrome, Brugada syndrome, and atrial fibrillation and advance their safe and successful translation into the cardiology arena.

Implementing Biological Pacemakers: Design Criteria for Successful

Each heartbeat that pumps blood throughout the body is initiated by an electrical impulse generated in the sinoatrial node (SAN). However, a number of disease conditions can hamper the ability of the SAN's pacemaker cells to generate consistent action potentials and maintain an orderly conduction path, leading to arrhythmias. For symptomatic patients, current treatments rely on implantation of an electronic pacing device. However, complications inherent to the indwelling hardware give pause to categorical use of device therapy for a subset of populations, including pediatric patients or those with temporary pacing needs. Cellular-based biological pacemakers, derived in vitro or in situ, could function as a therapeutic alternative to current electronic pacemakers. Understanding how biological pacemakers measure up to the SAN would facilitate defining and demonstrating its advantages over current treatments. In this review, we discuss recent approaches to creating biological pacemakers and delineate design criteria to guide future progress based on insights from basic biology of the SAN. We emphasize the need for long-term efficacy in vivo via maintenance of relevant proteins, source-sink balance, a niche reflective of the native SAN microenvironment, and chronotropic competence. With a focus on such criteria, combined with delivery methods tailored for disease indications, clinical implementation will be attainable.

Bioengineering the Cardiac Conduction System: Advances in Cellular, Gene, and Tissue Engineering for Heart Rhythm Regeneration

Heart rhythm disturbances caused by different etiologies may affect pediatric and adult patients with life-threatening consequences. When pharmacological therapy is ineffective in treating the disturbances, the implantation of electronic devices to control and/or restore normal heart pacing is a unique clinical management option. Although these artificial devices are life-saving, they display many limitations; not least, they do not have any capability to adapt to somatic growth or respond to neuroautonomic physiological changes. A biological pacemaker could offer a new clinical solution for restoring heart rhythms in the conditions of disorder in the cardiac conduction system. Several experimental approaches, such as cell-based, gene-based approaches, and the combination of both, for the generation of biological pacemakers are currently established and widely studied. Pacemaker bioengineering is also emerging as a technology to regenerate nodal tissues. This review analyzes and summarizes the strategies applied so far for the development of biological pacemakers, and discusses current translational challenges toward the first-in-human clinical application.

Transcriptional and Epigenetic Landscape of Cardiac Pacemaker Cells: Insights Into Cellular Specialization in the Sinoatrial Node

Cardiac pacemaker cells differentiate and functionally specialize early in embryonic development through activation of critical gene regulatory networks. In general, cellular specification and differentiation require that combinations of cell type-specific transcriptional regulators activate expression of key effector genes by binding to DNA regulatory elements including enhancers and promoters. However, because genomic DNA is tightly packaged by histones that must be covalently modified in order to render DNA regulatory elements and promoters accessible for transcription, the process of development and differentiation is intimately connected to the epigenetic regulation of chromatin accessibility. Although the difficulty of obtaining sufficient quantities of pure populations of pacemaker cells has limited progress in this field, the advent of low-input genomic technologies has the potential to catalyze a rapid growth of knowledge in this important area. The goal of this review is to outline the key transcriptional networks that control pacemaker cell development, with particular attention to our emerging understanding of how chromatin accessibility is modified and regulated during pacemaker cell differentiation. In addition, we will discuss the relevance of these findings to adult sinus node function, sinus node diseases, and origins of genetic variation in heart rhythm. Lastly, we will outline the current challenges facing this field and promising directions for future investigation.

Molecular Profiling of the Cardiac Conduction System: the Dawn of a New Era

PURPOSE OF REVIEW

Recent technological advances have led to an increased ability to define the gene expression profile of the cardiac conduction system (CCS). Here, we review the most salient studies to emerge in recent years and discuss existing gaps in our knowledge as well as future areas of investigation.

RECENT FINDINGS

Molecular profiling of the CCS spans several decades. However, the advent of high-throughput sequencing strategies has allowed for the discovery of unique transcriptional programs of the many diverse CCS cell types. The CCS, a diverse structure with significant inter- and intra-component cellular heterogeneity, is essential to the normal function of the heart. Progress in transcriptomic profiling has improved the resolution and depth of characterization of these unique and clinically relevant CCS cell types. Future studies leveraging this big data will play a crucial role in improving our understanding of CCS development and function as well as translating these findings into tangible translational tools for the improved detection, prevention, and treatment of cardiac arrhythmias.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

ATAC-Seq Reveals an Isl1 Enhancer That Regulates Sinoatrial Node Development and Function

RATIONALE

Cardiac pacemaker cells (PCs) in the sinoatrial node (SAN) have a distinct gene expression program that allows them to fire automatically and initiate the heartbeat. Although critical SAN transcription factors, including Isl1 (Islet-1), Tbx3 (T-box transcription factor 3), and Shox2 (short-stature homeobox protein 2), have been identified, the cis-regulatory architecture that governs PC-specific gene expression is not understood, and discrete enhancers required for gene regulation in the SAN have not been identified.

OBJECTIVE

To define the epigenetic profile of PCs using comparative ATAC-seq (assay for transposase-accessible chromatin with sequencing) and to identify novel enhancers involved in SAN gene regulation, development, and function.

METHODS AND RESULTS

We used ATAC-seq on sorted neonatal mouse SAN to compare regions of accessible chromatin in PCs and right atrial cardiomyocytes. PC-enriched assay for transposase-accessible chromatin peaks, representing candidate SAN regulatory elements, were located near established SAN genes and were enriched for distinct sets of TF (transcription factor) binding sites. Among several novel SAN enhancers that were experimentally validated using transgenic mice, we identified a 2.9-kb regulatory element at the Isl1 locus that was active specifically in the cardiac inflow at embryonic day 8.5 and throughout later SAN development and maturation. Deletion of this enhancer from the genome of mice resulted in SAN hypoplasia and sinus arrhythmias. The mouse SAN enhancer also directed reporter activity to the inflow tract in developing zebrafish hearts, demonstrating deep conservation of its upstream regulatory network. Finally, single nucleotide polymorphisms in the human genome that occur near the region syntenic to the mouse enhancer exhibit significant associations with resting heart rate in human populations.

CONCLUSIONS

(1) PCs have distinct regions of accessible chromatin that correlate with their gene expression profile and contain novel SAN enhancers, (2) cis-regulation of Isl1 specifically in the SAN depends upon a conserved SAN enhancer that regulates PC development and SAN function, and (3) a corresponding human ISL1 enhancer may regulate human SAN function.