Transcriptional regulation of the cardiac conduction system

Molecular Regulation of Cardiac Conduction System Development

PURPOSE OF REVIEW

The cardiac conduction system, composed of pacemaker cells and conducting cardiomyocytes, orchestrates the propagation of electrical activity to synchronize heartbeats. The conduction system plays a crucial role in the development of cardiac arrhythmias. In the embryo, the cells of the conduction system derive from the same cardiac progenitors as the contractile cardiomyocytes and and the key question is how this choice is made during development.

RECENT FINDINGS

This review focuses on recent advances in developmental biology using the mouse as animal model to better understand the cellular origin and molecular regulations that control morphogenesis of the cardiac conduction system, including the latest findings in single-cell transcriptomics. The conducting cell fate is acquired during development starting with pacemaking activity and last with the formation of a complex fast-conducting network. Cardiac conduction system morphogenesis is controlled by complex transcriptional and gene regulatory networks that differ in the components of the cardiac conduction system.

Revisiting Cardiac Biology in the Era of Single Cell and Spatial Omics

Throughout our lifetime, each beat of the heart requires the coordinated action of multiple cardiac cell types. Understanding cardiac cell biology, its intricate microenvironments, and the mechanisms that govern their function in health and disease are crucial to designing novel therapeutical and behavioral interventions. Recent advances in single-cell and spatial omics technologies have significantly propelled this understanding, offering novel insights into the cellular diversity and function and the complex interactions of cardiac tissue. This review provides a comprehensive overview of the cellular landscape of the heart, bridging the gap between suspension-based and emerging in situ approaches, focusing on the experimental and computational challenges, comparative analyses of mouse and human cardiac systems, and the rising contextualization of cardiac cells within their niches. As we explore the heart at this unprecedented resolution, integrating insights from both mouse and human studies will pave the way for novel diagnostic tools and therapeutic interventions, ultimately improving outcomes for patients with cardiovascular diseases.

A GABAergic system in atrioventricular node pacemaker cells controls electrical conduction between the atria and ventricles

Physiologically, the atria contract first, followed by the ventricles, which is the prerequisite for normal blood circulation. The above phenomenon of atrioventricular sequential contraction results from the characteristically slow conduction of electrical excitation of the atrioventricular node (AVN) between the atria and the ventricles. However, it is not clear what controls the conduction of electrical excitation within AVNs. Here, we find that AVN pacemaker cells (AVNPCs) possess an intact intrinsic GABAergic system, which plays a key role in electrical conduction from the atria to the ventricles. First, along with the discovery of abundant GABA-containing vesicles under the surface membranes of AVNPCs, key elements of the GABAergic system, including GABA metabolic enzymes, GABA receptors, and GABA transporters, were identified in AVNPCs. Second, GABA synchronously elicited GABA-gated currents in AVNPCs, which significantly weakened the excitability of AVNPCs. Third, the key molecular elements of the GABAergic system markedly modulated the conductivity of electrical excitation in the AVN. Fourth, GABAA receptor deficiency in AVNPCs accelerated atrioventricular conduction, which impaired the AVN's protective potential against rapid ventricular frequency responses, increased susceptibility to lethal ventricular arrhythmias, and decreased the cardiac contractile function. Finally, interventions targeting the GABAergic system effectively prevented the occurrence and development of atrioventricular block. In summary, the endogenous GABAergic system in AVNPCs determines the slow conduction of electrical excitation within AVNs, thereby ensuring sequential atrioventricular contraction. The endogenous GABAergic system shows promise as a novel intervention target for cardiac arrhythmias.

TBX3 transfection and nodal signal pathway inhibition promote differentiation of adipose mesenchymal stem cell to cardiac pacemaker-like cells

BACKGROUND

Mesenchymal stem cells (MSCs) are known as one of the best candidate cells to produce cardiac pacemaker-like cells (CPLCs). Upregulation of TBX3 transcription factor and inhibition of the nodal signal pathway have a significant role in the formation of cardiac pacemaker cells such as sinoatrial and atrioventricular nodes, which initiate the heartbeat and control the rhythm of heart contractions. This study aimed to confirm the effects of transfection of TBX3 transcription factor and inhibition of the nodal signal pathway on differentiating adipose-derived MSCs (AD-MSCs) to CPLCs. AD-MSCs were characterized using flow cytometry and three-lineage differentiation staining.

METHODS

The transfection of TBX3 plasmid was carried out using lipofectamine, and inhibition of the nodal signal pathway was done using the small-molecule SB431542. The morphology of the cells was observed using a light microscope. Pacemaker-specific markers, including TBX3, Cx30, HCN4, HCN1, HCN3, and KCNN4, were evaluated using the qRT-PCR method. For protein level, TBX3 and Cx30 were evaluated using ELISA and immunofluorescence staining. The electrophysiology of cells was evaluated using a patch clamp.

RESULTS

The TBX3 expression in the TBX3, SM, and TBX + SM groups significantly higher (p < 0.05) compared to the control group and cardiomyocytes. The expression of Cx40 and Cx43 genes were lower in TBX3, SM, TBX + SM groups. In contrast, Cx30 gene showed higher expression in TBX3 group. The expression HCN1, HCN3, and HCN4 genes are higher in TBX3 group.

CONCLUSION

The transfection of TBX3 and inhibition of the nodal signal pathway by small-molecule SB431542 enhanced differentiation of AD-MSCs to CPLCs.

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.

Canonical Wnt signaling directs the generation of functional human PSC-derived atrioventricular canal cardiomyocytes in bioprinted cardiac tissues

The creation of a functional 3D bioprinted human heart remains challenging, largely due to the lack of some crucial cardiac cell types, including the atrioventricular canal (AVC) cardiomyocytes, which are essential to slow down the electrical impulse between the atrium and ventricle. By utilizing single-cell RNA sequencing analysis and a 3D bioprinting technology, we discover that stage-specific activation of canonical Wnt signaling creates functional AVC cardiomyocytes derived from human pluripotent stem cells. These cardiomyocytes display morphological characteristics and express molecular markers of AVC cardiomyocytes, including transcription factors TBX2 and MSX2. When bioprinted in prefabricated cardiac tissues, these cardiomyocytes successfully delay the electrical impulse, demonstrating their capability of functioning as the AVC cardiomyocytes in vitro. Thus, these findings not only identify canonical Wnt signaling as a key regulator of the AVC cardiomyocyte differentiation in vitro, but, more importantly, provide a critical cellular source for the biofabrication of a functional human heart.

Cardiac Development at a Single-Cell Resolution

Mammalian cardiac development is a complex, multistage process. Though traditional lineage tracing studies have characterized the broad trajectories of cardiac progenitors, the advent and rapid optimization of single-cell RNA sequencing methods have yielded an ever-expanding toolkit for characterizing heterogeneous cell populations in the developing heart. Importantly, they have allowed for a robust profiling of the spatiotemporal transcriptomic landscape of the human and mouse heart, revealing the diversity of cardiac cells-myocyte and non-myocyte-over the course of development. These studies have yielded insights into novel cardiac progenitor populations, chamber-specific developmental signatures, the gene regulatory networks governing cardiac development, and, thus, the etiologies of congenital heart diseases. Furthermore, single-cell RNA sequencing has allowed for the exquisite characterization of distinct cardiac populations such as the hard-to-capture cardiac conduction system and the intracardiac immune population. Therefore, single-cell profiling has also resulted in new insights into the regulation of cardiac regeneration and injury repair. Single-cell multiomics approaches combining transcriptomics, genomics, and epigenomics may uncover an even more comprehensive atlas of human cardiac biology. Single-cell analyses of the developing and adult mammalian heart offer an unprecedented look into the fundamental mechanisms of cardiac development and the complex diseases that may arise from it.

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.

Dbh+ catecholaminergic cardiomyocytes contribute to the structure and function of the cardiac conduction system in murine heart

The heterogeneity of functional cardiomyocytes arises during heart development, which is essential to the complex and highly coordinated cardiac physiological function. Yet the biological and physiological identities and the origin of the specialized cardiomyocyte populations have not been fully comprehended. Here we report a previously unrecognised population of cardiomyocytes expressing Dbhgene encoding dopamine beta-hydroxylase in murine heart. We determined how these myocytes are distributed across the heart by utilising advanced single-cell and spatial transcriptomic analyses, genetic fate mapping and molecular imaging with computational reconstruction. We demonstrated that they form the key functional components of the cardiac conduction system by using optogenetic electrophysiology and conditional cardiomyocyte Dbh gene deletion models. We revealed their close relationship with sympathetic innervation during cardiac conduction system formation. Our study thus provides new insights into the development and heterogeneity of the mammalian cardiac conduction system by revealing a new cardiomyocyte population with potential catecholaminergic endocrine function.

Single-cell RNA sequencing of murine hearts for studying the development of the cardiac conduction system

The development of the cardiac conduction system (CCS) is essential for correct heart function. However, critical details on the cell types populating the CCS in the mammalian heart during the development remain to be resolved. Using single-cell RNA sequencing, we generated a large dataset of transcriptomes of ~0.5 million individual cells isolated from murine hearts at six successive developmental corresponding to the early, middle and late stages of heart development. The dataset provides a powerful library for studying the development of the heart's CCS and other cardiac components. Our initial analysis identified distinct cell types between 20 to 26 cell types across different stages, of which ten are involved in forming the CCS. Our dataset allows researchers to reuse the datasets for data mining and a wide range of analyses. Collectively, our data add valuable transcriptomic resources for further study of cardiac development, such as gene expression, transcriptional regulation and functional gene activity in developing hearts, particularly the CCS.

Efficient generation of TBX3+ atrioventricular conduction-like cardiomyocytes from human pluripotent stem cells

Atrioventricular conduction cardiomyocytes (AVCCs) regulate the rate and rhythm of heart contractions. Dysfunction due to aging or disease can cause atrioventricular (AV) block, interrupting electrical impulses from the atria to the ventricles. Generation of functional atrioventricular conduction like cardiomyocytes (AVCLCs) from human pluripotent stem cells (hPSCs) provides a promising approach to repair damaged atrioventricular conduction tissue by cell transplantation. In this study, we put forward the generation of AVCLCs from hPSCs by stage-specific manipulation of the retinoic acid (RA), WNT, and bone morphogenetic protein (BMP) signaling pathways. These cells express AVCC-specific markers, including the transcription factors TBX3, MSX2 and NKX2.5, display functional electrophysiological characteristics and present low conduction velocity (0.07 ± 0.02 m/s). Our findings provide new insights into the understanding of the development of the atrioventricular conduction system and propose a strategy for the treatment of severe atrioventricular conduction block by cell transplantation in future.

Spatially resolved multiomics of human cardiac niches

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

An evolutionarily conserved pacemaker role for HCN ion channels in smooth muscle

Although hyperpolarization-activated cation (HCN) ion channels are well established to underlie cardiac pacemaker activity, their role in smooth muscle organs remains controversial. HCN-expressing cells are localized to renal pelvic smooth muscle (RPSM) pacemaker tissues of the murine upper urinary tract and HCN channel conductance is required for peristalsis. To date, however, the Ih pacemaker current conducted by HCN channels has never been detected in these cells, raising questions on the identity of RPSM pacemakers. Indeed, the RPSM pacemaker mechanisms of the unique multicalyceal upper urinary tract exhibited by humans remains unknown. Here, we developed immunopanning purification protocols and demonstrate that 96% of isolated HCN+ cells exhibit Ih . Single-molecule STORM to whole-tissue imaging showed HCN+ cells express single HCN channels on their plasma membrane and integrate into the muscular syncytium. By contrast, PDGFR-α+ cells exhibiting the morphology of ICC gut pacemakers were shown to be vascular mural cells. Translational studies in the homologous human and porcine multicalyceal upper urinary tracts showed that contractions and pacemaker depolarizations originate in proximal calyceal RPSM. Critically, HCN+ cells were shown to integrate into calyceal RPSM pacemaker tissues, and HCN channel block abolished electrical pacemaker activity and peristalsis of the multicalyceal upper urinary tract. Cumulatively, these studies demonstrate that HCN ion channels play a broad, evolutionarily conserved pacemaker role in both cardiac and smooth muscle organs and have implications for channelopathies as putative aetiologies of smooth muscle disorders. KEY POINTS: Pacemakers trigger contractions of involuntary muscles. Hyperpolarization-activated cation (HCN) ion channels underpin cardiac pacemaker activity, but their role in smooth muscle organs remains controversial. Renal pelvic smooth muscle (RPSM) pacemakers trigger contractions that propel waste away from the kidney. HCN+ cells localize to murine RPSM pacemaker tissue and HCN channel conductance is required for peristalsis. The HCN (Ih ) current has never been detected in RPSM cells, raising doubt whether HCN+ cells are bona fide pacemakers. Moreover, the pacemaker mechanisms of the unique multicalyceal RPSM of higher order mammals remains unknown. In total, 97% of purified HCN+ RPSM cells exhibit Ih . HCN+ cells integrate into the RPSM musculature, and pacemaker tissue peristalsis is dependent on HCN channels. Translational studies in human and swine demonstrate HCN channels are conserved in the multicalyceal RPSM and that HCN channels underlie pacemaker activity that drives peristalsis. These studies provide insight into putative channelopathies that can underlie smooth muscle dysfunction.

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.

An atrial fibrillation-associated regulatory region modulates cardiac Tbx5 levels and arrhythmia susceptibility

Heart development and rhythm control are highly Tbx5 dosage-sensitive. TBX5 haploinsufficiency causes congenital conduction disorders, whereas increased expression levels of TBX5 in human heart samples has been associated with atrial fibrillation (AF). We deleted the conserved mouse orthologues of two independent AF-associated genomic regions in the Tbx5 locus, one intronic (RE(int)) and one downstream (RE(down)) of Tbx5. In both lines, we observed a modest (30%) increase of Tbx5 in the postnatal atria. To gain insight into the effects of slight dosage increase in vivo, we investigated the atrial transcriptional, epigenetic and electrophysiological properties of both lines. Increased atrial Tbx5 expression was associated with induction of genes involved in development, ion transport and conduction, with increased susceptibility to atrial arrhythmias, and increased action potential duration of atrial cardiomyocytes. We identified an AF-associated variant in the human RE(int) that increases its transcriptional activity. Expression of the AF-associated transcription factor Prrx1 was induced in Tbx5^RE(int)KO cardiomyocytes. We found that some of the transcriptional and functional changes in the atria caused by increased Tbx5 expression were normalized when reducing cardiac Prrx1 expression in Tbx5^RE(int)KO mice, indicating an interaction between these two AF genes. We conclude that modest increases in expression of dose-dependent transcription factors, caused by common regulatory variants, significantly impact on the cardiac gene regulatory network and disease susceptibility.

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.

Direct reprogramming of cardiomyocytes into cardiac Purkinje-like cells

Currently, there are no treatments that ameliorate cardiac cell death, the underlying basis of cardiovascular disease. An unexplored cell type in cardiac regeneration is cardiac Purkinje cells; specialized cells from the cardiac conduction system (CCS) responsible for propagating electrical signals. Purkinje cells have tremendous potential as a regenerative treatment because they may intrinsically integrate with the CCS of a recipient myocardium, resulting in more efficient electrical conduction in diseased hearts. This study is the first to demonstrate an effective protocol for the direct reprogramming of human cardiomyocytes into cardiac Purkinje-like cells using small molecules. The cells generated were genetically and functionally similar to native cardiac Purkinje cells, where expression of key cardiac Purkinje genes such as CNTN2, ETV1, PCP4, IRX3, SCN5a, HCN2 and the conduction of electrical signals with increased velocity was observed. This study may help to advance the quest to finding an optimized cell therapy for heart regeneration.

Molecular and Functional Remodeling of Superior and Inferior SAN in a Rat Model of HCM

BACKGROUND

Recently, our laboratory presented functional and molecular evidence for the presence of 2 competing sinoatrial node (SAN) pacemakers in healthy human and rat hearts. Anatomically localized near the superior vena cava and inferior vena cava, the superior and inferior SANs (sSAN and iSAN, respectively) preferentially control fast and slow normal heart rates. However, only 1 dominant pacemaker, primarily the sSAN, was functional in the failing rat heart with hypertrophic cardiomyopathy.

OBJECTIVES

This study aimed to determine the transcriptional basis of functional silencing of 1 of 2 dominant pacemakers in failing rat hearts.

METHODS

Ascending aortic constriction was performed on 1-week-old male Sprague-Dawley rat pups to induce left ventricular hypertrophy and heart failure. The dominant pacemaker was anatomically mapped in adult (10-12 weeks old) healthy and failing rat hearts using optical mapping in isolated right atrial tissue preparations. RNA sequencing was used to identify regional sSAN/iSAN gene expression differences between healthy and failing rat hearts.

RESULTS

In all failing rat hearts optically mapped in this study (n = 4), only the sSAN pacemaker was functional, while the iSAN was silent. Compared to healthy rat hearts, a total of 3,640 genes were downregulated, and 4,518 genes were upregulated in failing rat hearts. The functional quiescence of the iSAN in these failing rat hearts may be explained by their downregulation of sodium, potassium, and calcium ion channels as well as their downregulation of specific structural genes, including ankyrin, titin, and myosin heavy chain. Moreover, the iSAN showed predominant downregulation of several key transcription factors such as Tbx5, Tbx3, Shox2, and Smad9.

CONCLUSIONS

Pressure-overload-induced heart failure resulted in significant downregulation of critical transcription factors, ion channels, and structural transcripts of the iSAN, which could explain the functional silencing of the iSAN in failing rat hearts.

In vivo visualization and molecular targeting of the cardiac conduction system

Accidental injury to the cardiac conduction system (CCS), a network of specialized cells embedded within the heart and indistinguishable from the surrounding heart muscle tissue, is a major complication in cardiac surgeries. Here, we addressed this unmet need by engineering targeted antibody-dye conjugates directed against the CCS, allowing for the visualization of the CCS in vivo following a single intravenous injection in mice. These optical imaging tools showed high sensitivity, specificity, and resolution, with no adverse effects on CCS function. Further, with the goal of creating a viable prototype for human use, we generated a fully human monoclonal Fab that similarly targets the CCS with high specificity. We demonstrate that, when conjugated to an alternative cargo, this Fab can also be used to modulate CCS biology in vivo, providing a proof of principle for targeted cardiac therapeutics. Finally, in performing differential gene expression analyses of the entire murine CCS at single-cell resolution, we uncovered and validated a suite of additional cell surface markers that can be used to molecularly target the distinct subcomponents of the CCS, each prone to distinct life-threatening arrhythmias. These findings lay the foundation for translational approaches targeting the CCS for visualization and therapy in cardiothoracic surgery, cardiac imaging, and arrhythmia management.

A single cell transcriptional roadmap of human pacemaker cell differentiation

Each heartbeat is triggered by the sinoatrial node (SAN), the primary pacemaker of the heart. Studies in animal models have revealed that pacemaker cells share a common progenitor with the (pro)epicardium, and that the pacemaker cardiomyocytes further diversify into ‘transitional’, ‘tail’, and ‘head’ subtypes. However, the underlying molecular mechanisms, especially of human pacemaker cell development, are poorly understood. Here, we performed single cell RNA sequencing (scRNA-seq) and trajectory inference on human induced pluripotent stem cells (hiPSCs) differentiating to SAN-like cardiomyocytes (SANCMs) to construct a roadmap of transcriptional changes and lineage decisions. In differentiated SANCM, we identified distinct clusters that closely resemble different subpopulations of the in vivo SAN. Moreover, the presence of a side population of proepicardial cells suggested their shared ontogeny with SANCM, as also reported in vivo. Our results demonstrate that the divergence of SANCM and proepicardial lineages is determined by WNT signaling. Furthermore, we uncovered roles for TGFβ and WNT signaling in the branching of transitional and head SANCM subtypes, respectively. These findings provide new insights into the molecular processes involved in human pacemaker cell differentiation, opening new avenues for complex disease modeling in vitro and inform approaches for cell therapy-based regeneration of the SAN.

Pregestational diabetes alters cardiac structure and function of neonatal rats through developmental plasticity

Pregestational diabetes (PGDM) leads to developmental impairment, especially cardiac dysfunction, in their offspring. The hyperglycemic microenvironment inside the uterus alters the cardiac plasticity characterized by electrical and structural remodeling of the heart. The altered expression of several transcription factors due to hyperglycemia during fetal development might be responsible for molecular defects and phenotypic changes in the heart. The molecular mechanism of the developmental defects in the heart due to PGDM remains unclear. To understand the molecular defects in the 2-days old neonatal rats, streptozotocin-induced diabetic female rats were bred with healthy male rats. We collected 2-day-old hearts from the neonates and identified the molecular basis for phenotypic changes. Neonates from diabetic mothers showed altered electrocardiography and echocardiography parameters. Transcriptomic profiling of the RNA-seq data revealed that several altered genes were associated with heart development, myocardial fibrosis, cardiac conduction, and cell proliferation. Histopathology data showed the presence of focal cardiac fibrosis and increased cell proliferation in neonates from diabetic mothers. Thus, our results provide a comprehensive map of the cellular events and molecular pathways perturbed in the neonatal heart during PGDM. All of the molecular and structural changes lead to developmental plasticity in neonatal rat hearts and develop cardiac anomalies in their early life.

Application of Three-Dimensional Culture Method in the Cardiac Conduction System Research

Congenital heart defects (CHD) are the most common type of birth defects. Several human case studies and genetically altered animal models have identified abnormalities in the development of ventricular conduction system (VCS) in the heart. While cell-based therapies hold promise for treating CHDs, translational efforts are limited by the lack of suitable in vitro models for feasibility and safety studies. A better understanding of cell differentiation pathways can lead to development of cell-based therapies for individuals living with CHD/VCS disorders. Here, we describe a new and reproducible 3-D cell culture method for studying cardiac cell lineage differentiation in vitro. We used primary ventricular cells isolated from embryonic day 11.5 (E11.5) mouse embryos, which can differentiate into multiple cardiac cell types including VCS cells. We compared 3-D cultures with three types of basement membrane extracts (BME) for their abilities to support E11.5 ventricular cell differentiation. In addition, the effects of atrial natriuretic peptide (ANP) and an inhibitor for its high affinity receptor were tested on cell differentiation in 3-D cultures. Following the cell culture, protocols for immunofluorescence imaging, cell extraction and protein isolation from the 3-D culture matrix and in-cell western methods are described. Further, these approaches can be used to study the effects of various ligands and genetic interventions on VCS cell development. We propose that these methodologies may also be extended for differentiation studies using other sources of stem cells such as induced pluripotent stem cells.

Genetic analysis of right heart structure and function in 40,000 people

Congenital heart diseases often involve maldevelopment of the evolutionarily recent right heart chamber. To gain insight into right heart structure and function, we fine-tuned deep learning models to recognize the right atrium, right ventricle and pulmonary artery, measuring right heart structures in 40,000 individuals from the UK Biobank with magnetic resonance imaging. Genome-wide association studies identified 130 distinct loci associated with at least one right heart measurement, of which 72 were not associated with left heart structures. Loci were found near genes previously linked with congenital heart disease, including NKX2-5, TBX5/TBX3, WNT9B and GATA4. A genome-wide polygenic predictor of right ventricular ejection fraction was associated with incident dilated cardiomyopathy (hazard ratio, 1.33 per standard deviation; P = 7.1 × 10-13) and remained significant after accounting for a left ventricular polygenic score. Harnessing deep learning to perform large-scale cardiac phenotyping, our results yield insights into the genetic determinants of right heart structure and function.

Rare and Common Variants Uncover the Role of the Atria in Coarctation of the Aorta

Coarctation of the aorta (CoA) and bicuspid aortic valve (BAV) often cooccur and are genetically linked congenital heart defects (CHD). While CoA is thought to have a hemodynamic origin from ventricular dysfunction, we provide evidence pointing to atrial hemodynamics based on investigating the genetic etiology of CoA. Previous studies have shown a rare MYH6 variant in an Icelandic cohort, and two common deletions in the protocadherin α cluster (PCDHA delCNVs) are significantly associated with CoA and BAV. Here, analysis of a non-Icelandic white CHD cohort (n = 166) recovered rare MYH6 variants in 10.9% of CoA and 32.7% of BAV/CoA patients, yielding odds ratios of 18.6 (p = 2.5 × 10⁻⁷) and 20.5 (p = 7.4 × 10⁻⁵) for the respective association of MYH6 variants with CoA and BAV/CoA. In combination with the PCHDA delCNVs, they accounted for a third of CoA cases. Gene expression datasets for the human and mouse embryonic heart showed that both genes are predominantly expressed in the atria, not the ventricle. Moreover, cis-eQTLs analysis showed the PCHDA delCNV is associated with reduced atrial expression of PCHDA10, a gene in the delCNV interval. Together, these findings showed that PCDHA/MYH6 variants account for a substantial fraction of CoA cases. An atrial rather than ventricular hemodynamic model for CoA is indicated, consistent with the known early atrial functional dominance of the human embryonic heart.

Getting to the heart of rhythm: A century of progress

The human heart beats over eighty thousand times a day, and the average person's heart may have beaten up to 3 billion times by the age of 80. During the early stages of pregnancy, the heart beat provides the first visual and auditory sign of life of the foetus. Conversely, the first audible sound that the foetus is likely to hear is the heart beat of the mother. How fitting then, that at the "birth" Physiological Reviews the very first article published in 1921 written by Eyster and Meek addressed "The origin and conduction of the heart beat".¹ In their insightful review, the authors discussed the landmark discoveries made from the mid-19^th century on the electrical function of the heart. Now, a hundred years later, at the start of the next century of Physiological Reviews, an update on the huge progress made in the "exciting" field of cardiac electrophysiology is warranted. Guided by a number of excellent reviews published in Physiological Reviews since 1921 as well as a large body of literature, an overview of the important advancements made on the topic is provided here.

GOing Forward With the Cardiac Conduction System Using Gene Ontology

The cardiac conduction system (CCS) comprises critical components responsible for the initiation, propagation, and coordination of the action potential. Aberrant CCS development can cause conduction abnormalities, including sick sinus syndrome, accessory pathways, and atrioventricular and bundle branch blocks. Gene Ontology (GO; http://geneontology.org/) is an invaluable global bioinformatics resource which provides structured, computable knowledge describing the functions of gene products. Many gene products are known to be involved in CCS development; however, this information is not comprehensively captured by GO. To address the needs of the heart development research community, this study aimed to describe the specific roles of proteins reported in the literature to be involved with CCS development and/or function. 14 proteins were prioritized for GO annotation which led to the curation of 15 peer-reviewed primary experimental articles using carefully selected GO terms. 152 descriptive GO annotations, including those describing sinoatrial node and atrioventricular node development were created and submitted to the GO Consortium database. A functional enrichment analysis of 35 key CCS development proteins confirmed that this work has improved the in-silico interpretation of this CCS dataset. This work may improve future investigations of the CCS with application of high-throughput methods such as genome-wide association studies analysis, proteomics, and transcriptomics.

Genome-wide association analyses identify new Brugada syndrome risk loci and highlight a new mechanism of sodium channel regulation in disease susceptibility

Brugada syndrome (BrS) is a cardiac arrhythmia disorder associated with sudden death in young adults. With the exception of SCN5A, encoding the cardiac sodium channel Na_V1.5, susceptibility genes remain largely unknown. Here we performed a genome-wide association meta-analysis comprising 2,820 unrelated cases with BrS and 10,001 controls, and identified 21 association signals at 12 loci (10 new). Single nucleotide polymorphism (SNP)-heritability estimates indicate a strong polygenic influence. Polygenic risk score analyses based on the 21 susceptibility variants demonstrate varying cumulative contribution of common risk alleles among different patient subgroups, as well as genetic associations with cardiac electrical traits and disorders in the general population. The predominance of cardiac transcription factor loci indicates that transcriptional regulation is a key feature of BrS pathogenesis. Furthermore, functional studies conducted on MAPRE2, encoding the microtubule plus-end binding protein EB2, point to microtubule-related trafficking effects on Na_V1.5 expression as a new underlying molecular mechanism. Taken together, these findings broaden our understanding of the genetic architecture of BrS and provide new insights into its molecular underpinnings.

Patient-specific TBX5-G125R Variant Induces Profound Transcriptional Deregulation and Atrial Dysfunction

Background: The pathogenic missense variant p.G125R in TBX5 causes Holt-Oram syndrome (HOS; hand-heart syndrome) and early onset of atrial fibrillation. Revealing how an altered key developmental transcription factor modulates cardiac physiology in vivo will provide unique insights into the mechanisms underlying atrial fibrillation in these patients. Methods: We analyzed electrocardiograms (ECGs) of an extended family pedigree of HOS patients. Next, we introduced the TBX5-p.G125R variant in the mouse genome (Tbx5^G125R) and performed electrophysiological analyses (ECG, optical mapping, patch clamp, intracellular calcium measurements), transcriptomics (single nuclei and tissue RNA sequencing) and epigenetic profiling (ATAC-sequencing, H3K27ac CUT&RUN-sequencing). Results: We discovered high incidence of atrial extra systoles and atrioventricular conduction disturbances in HOS patients. Tbx5^G125R/+ mice were morphologically unaffected and displayed variable RR intervals, atrial extra systoles and susceptibility to atrial fibrillation, reminiscent of TBX5-p.G125R patients. Atrial conduction velocity was not affected but systolic and diastolic intracellular calcium concentrations were decreased and action potentials prolonged in isolated cardiomyocytes of Tbx5^G125R/+ mice compared to controls. Transcriptional profiling of atria revealed most profound transcriptional changes in cardiomyocytes versus other cell types, and identified over a thousand coding and non-coding transcripts that were differentially expressed. Epigenetic profiling uncovered thousands of TBX5-p.G125R sensitive putative regulatory elements (including enhancers) that gained accessibility in atrial cardiomyocytes. The majority of sites with increased accessibility were occupied by Tbx5. The small group of sites with reduced accessibility was enriched for DNA binding motifs of members of the SP- and KLF families of transcription factors. These data show that Tbx5-p.G125R induces changes in regulatory element activity, altered transcriptional regulation and changed cardiomyocyte behavior, possibly caused by altered DNA binding and cooperativity properties. Conclusions: Our data reveal how a disease-causing missense variant in TBX5 induces profound changes in the atrial transcriptional regulatory network and epigenetic state in vivo, leading to arrhythmia reminiscent of those seen in human TBX5-p.G125R variant carriers.

Modeling Human Cardiac Arrhythmias: Insights from Zebrafish

Cardiac arrhythmia, or irregular heart rhythm, is associated with morbidity and mortality and is described as one of the most important future public health challenges. Therefore, developing new models of cardiac arrhythmia is critical for understanding disease mechanisms, determining genetic underpinnings, and developing new therapeutic strategies. In the last few decades, the zebrafish has emerged as an attractive model to reproduce in vivo human cardiac pathologies, including arrhythmias. Here, we highlight the contribution of zebrafish to the field and discuss the available cardiac arrhythmia models. Further, we outline techniques to assess potential heart rhythm defects in larval and adult zebrafish. As genetic tools in zebrafish continue to bloom, this model will be crucial for functional genomics studies and to develop personalized anti-arrhythmic therapies.

Tbx18 promoted the conversion of human-induced pluripotent stem cell-derived cardiomyocytes into sinoatrial node-like pacemaker cells

Sinoatrial node (SAN) pacemaker cells originate from T-box transcription factor 18 (Tbx18)-expressing progenitor cells. The present study aimed to investigate whether overexpression of human transcription factor Tbx18 could reprogram human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) into SAN-like pacemaker cells (SANLPCs) in vitro. In the study, hiPSCs were first differentiated into hiPSC-CMs through regulating the Wnt/β-catenin pathway, then purified hiPSC-CMs were transfected by Tbx18 adenovirus (Tbx18-CMs group) or green fluorescent protein (GFP) adenovirus (GFP-CMs group). The beating frequency of the Tbx18-CMs group was significantly higher than that of the hiPSC-CMs group and GFP-CMs group. Compared with the other two groups, the expression levels of hyperpolarization-activated cyclic nucleotide-gated potassium channel isoform 4, connexin-45 in the Tbx18-CMs group were markedly upregulated, while the expressions of transcription factor NKX2.5, CX43 were significantly downregulated. Whole-cell patch-clamp results illustrated that action potential and "funny" current (I_f ) similar to SAN pacemaker cells could be recorded in the Tbx18-CMs group. In conclusion, this present study demonstrated that overexpression of Tbx18 promoted the conversion of hiPSC-CMs into SANLPCs.

Further insights into the molecular complexity of the human sinus node - The role of 'novel' transcription factors and microRNAs

RESEARCH PURPOSE

The sinus node (SN) is the heart's primary pacemaker. Key ion channels (mainly the funny channel, HCN4) and Ca2+-handling proteins in the SN are responsible for its function. Transcription factors (TFs) regulate gene expression through inhibition or activation and microRNAs (miRs) do this through inhibition. There is high expression of macrophages and mast cells within the SN connective tissue. 'Novel'/unexplored TFs and miRs in the regulation of ion channels and immune cells in the SN are not well understood. Using RNAseq and bioinformatics, the expression profile and predicted interaction of key TFs and cell markers with key miRs in the adult human SN vs. right atrial tissue (RA) were determined.

PRINCIPAL RESULTS

68 and 60 TFs significantly more or less expressed in the SN vs. RA respectively. Among those more expressed were ISL1 and TBX3 (involved in embryonic development of the SN) and 'novel' RUNX1-2, CEBPA, GLI1-2 and SOX2. These TFs were predicted to regulate HCN4 expression in the SN. Markers for different cells: fibroblasts (COL1A1), fat (FABP4), macrophages (CSF1R and CD209), natural killer (GZMA) and mast (TPSAB1) were significantly more expressed in the SN vs. RA. Interestingly, RUNX1-3, CEBPA and GLI1 also regulate expression of these cells. MiR-486-3p inhibits HCN4 and markers involved in immune response.

MAJOR CONCLUSIONS

In conclusion, RUNX1-2, CSF1R, TPSAB1, COL1A1 and HCN4 are highly expressed in the SN but not miR-486-3p. Their complex interactions can be used to treat SN dysfunction such as bradycardia. Interestingly, another research group recently reported miR-486-3p is upregulated in blood samples from severe COVID-19 patients who suffer from bradycardia.

New Insights into the Development and Morphogenesis of the Cardiac Purkinje Fiber Network: Linking Architecture and Function

The rapid propagation of electrical activity through the ventricular conduction system (VCS) controls spatiotemporal contraction of the ventricles. Cardiac conduction defects or arrhythmias in humans are often associated with mutations in key cardiac transcription factors that have been shown to play important roles in VCS morphogenesis in mice. Understanding of the mechanisms of VCS development is thus crucial to decipher the etiology of conduction disturbances in adults. During embryogenesis, the VCS, consisting of the His bundle, bundle branches, and the distal Purkinje network, originates from two independent progenitor populations in the primary ring and the ventricular trabeculae. Differentiation into fast-conducting cardiomyocytes occurs progressively as ventricles develop to form a unique electrical pathway at late fetal stages. The objectives of this review are to highlight the structure–function relationship between VCS morphogenesis and conduction defects and to discuss recent data on the origin and development of the VCS with a focus on the distal Purkinje fiber network.

Human model of IRX5 mutations reveals key role for this transcription factor in ventricular conduction

AIMS

Several inherited arrhythmic diseases have been linked to single gene mutations in cardiac ion channels and interacting proteins. However, the mechanisms underlying most arrhythmias, are thought to involve altered regulation of the expression of multiple effectors. In this study, we aimed to examine the role of a transcription factor (TF) belonging to the Iroquois homeobox family, IRX5, in cardiac electrical function.

METHODS AND RESULTS

Using human cardiac tissues, transcriptomic correlative analyses between IRX5 and genes involved in cardiac electrical activity showed that in human ventricular compartment, IRX5 expression strongly correlated to the expression of major actors of cardiac conduction, including the sodium channel, Nav1.5, and Connexin 40 (Cx40). We then generated human-induced pluripotent stem cells (hiPSCs) derived from two Hamamy syndrome-affected patients carrying distinct homozygous loss-of-function mutations in IRX5 gene. Cardiomyocytes derived from these hiPSCs showed impaired cardiac gene expression programme, including misregulation in the control of Nav1.5 and Cx40 expression. In accordance with the prolonged QRS interval observed in Hamamy syndrome patients, a slower ventricular action potential depolarization due to sodium current reduction was observed on electrophysiological analyses performed on patient-derived cardiomyocytes, confirming the functional role of IRX5 in electrical conduction. Finally, a cardiac TF complex was newly identified, composed by IRX5 and GATA4, in which IRX5 potentiated GATA4-induction of SCN5A expression.

CONCLUSION

Altogether, this work unveils a key role for IRX5 in the regulation of human ventricular depolarization and cardiac electrical conduction, providing therefore new insights into our understanding of cardiac diseases.

Sex differences in the origin of Purkinje ectopy-initiated idiopathic ventricular fibrillation

BACKGROUND

Purkinje ectopics (PurkEs) are major triggers of idiopathic ventricular fibrillation (VF). Identifying clinical factors associated with specific PurkE characteristics could yield insights into the mechanisms of Purkinje-mediated arrhythmogenicity.

OBJECTIVE

The purpose of this study was to examine the associations of clinical, environmental, and genetic factors with PurkE origin in patients with PurkE-initiated idiopathic VF.

METHODS

Consecutive patients with PurkE-initiated idiopathic VF from 4 arrhythmia referral centers were included. We evaluated demographic characteristics, medical history, clinical circumstances associated with index VF events, and electrophysiological characteristics of PurkEs. An electrophysiology study was performed in most patients to confirm the Purkinje origin.

RESULTS

Eighty-three patients were included (mean age 38 ± 14 years; 44 [53%] women), of whom 32 had a history of syncope. Forty-four patients had VF at rest. PurkEs originated from the right ventricle (RV) in 41 patients (49%), from the left ventricle (LV) in 36 (44%), and from both ventricles in 6 (7%). Seasonal and circadian distributions of VF episodes were similar according to PurkE origin. The clinical characteristics of patients with RV vs LV PurkE origins were similar, except for sex. RV PurkEs were more frequent in men than in women (76% vs 24%), whereas LV and biventricular PurkEs were more frequent in women (81% vs 19% and 83% vs 17%, respectively) (P < .0001).

CONCLUSION

PurkEs triggering idiopathic VF originate dominantly from the RV in men and from the LV or both ventricles in women, adding to other sex-related arrhythmias such as Brugada syndrome or long QT syndrome. Sex-based factors influencing Purkinje arrhythmogenicity warrant 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.

Hedgehog signalling controls sinoatrial node development and atrioventricular cushion formation

Smoothened is a key receptor of the hedgehog pathway, but the roles of Smoothened in cardiac development remain incompletely understood. In this study, we found that the conditional knockout of Smoothened from the mesoderm impaired the development of the venous pole of the heart and resulted in hypoplasia of the atrium/inflow tract (IFT) and a low heart rate. The blockage of Smoothened led to reduced expression of genes critical for sinoatrial node (SAN) development in the IFT. In a cardiac cell culture model, we identified a Gli2–Tbx5–Hcn4 pathway that controls SAN development. In the mutant embryos, the endocardial-to-mesenchymal transition (EndMT) in the atrioventricular cushion failed, and Bmp signalling was downregulated. The addition of Bmp2 rescued the EndMT in mutant explant cultures. Furthermore, we analysed Gli2⁺ scRNAseq and Tbx5^−/− RNAseq data and explored the potential genes downstream of hedgehog signalling in posterior second heart field derivatives. In conclusion, our study reveals that Smoothened-mediated hedgehog signalling controls posterior cardiac progenitor commitment, which suggests that the mutation of Smoothened might be involved in the aetiology of congenital heart diseases related to the cardiac conduction system and heart valves.

Understanding the Adult Mammalian Heart at Single-Cell RNA-Seq Resolution

During the last decade, extensive efforts have been made to comprehend cardiac cell genetic and functional diversity. Such knowledge allows for the definition of the cardiac cellular interactome as a reasonable strategy to increase our understanding of the normal and pathologic heart. Previous experimental approaches including cell lineage tracing, flow cytometry, and bulk RNA-Seq have often tackled the analysis of cardiac cell diversity as based on the assumption that cell types can be identified by the expression of a single gene. More recently, however, the emergence of single-cell RNA-Seq technology has led us to explore the diversity of individual cells, enabling the cardiovascular research community to redefine cardiac cell subpopulations and identify relevant ones, and even novel cell types, through their cell-specific transcriptomic signatures in an unbiased manner. These findings are changing our understanding of cell composition and in consequence the identification of potential therapeutic targets for different cardiac diseases. In this review, we provide an overview of the continuously changing cardiac cellular landscape, traveling from the pre-single-cell RNA-Seq times to the single cell-RNA-Seq revolution, and discuss the utilities and limitations of this technology.

Tissue clearing and imaging methods for cardiovascular development

Tissue imaging in 3D using visible light is limited and various clearing techniques were developed to increase imaging depth, but none provides universal solution for all tissues at all developmental stages. In this review, we focus on different tissue clearing methods for 3D imaging of heart and vasculature, based on chemical composition (solvent-based, simple immersion, hyperhydration, and hydrogel embedding techniques). We discuss in detail compatibility of various tissue clearing techniques with visualization methods: fluorescence preservation, immunohistochemistry, nuclear staining, and fluorescent dyes vascular perfusion. We also discuss myocardium visualization using autofluorescence, tissue shrinking, and expansion. Then we overview imaging methods used to study cardiovascular system and live imaging. We discuss heart and vessels segmentation methods and image analysis. The review covers the whole process of cardiovascular system 3D imaging, starting from tissue clearing and its compatibility with various visualization methods to the types of imaging methods and resulting image analysis.

Assembly of the Cardiac Pacemaking Complex: Electrogenic Principles of Sinoatrial Node Morphogenesis

Cardiac pacemaker cells located in the sinoatrial node initiate the electrical impulses that drive rhythmic contraction of the heart. The sinoatrial node accounts for only a small proportion of the total mass of the heart yet must produce a stimulus of sufficient strength to stimulate the entire volume of downstream cardiac tissue. This requires balancing a delicate set of electrical interactions both within the sinoatrial node and with the downstream working myocardium. Understanding the fundamental features of these interactions is critical for defining vulnerabilities that arise in human arrhythmic disease and may provide insight towards the design and implementation of the next generation of potential cellular-based cardiac therapeutics. Here, we discuss physiological conditions that influence electrical impulse generation and propagation in the sinoatrial node and describe developmental events that construct the tissue-level architecture that appears necessary for sinoatrial node function.

Bioengineering Clinically Relevant Cardiomyocytes and Cardiac Tissues from Pluripotent Stem Cells

The regenerative capacity of cardiomyocytes is insufficient to functionally recover damaged tissue, and as such, ischaemic heart disease forms the largest proportion of cardiovascular associated deaths. Human-induced pluripotent stem cells (hiPSCs) have enormous potential for developing patient specific cardiomyocytes for modelling heart disease, patient-based cardiac toxicity testing and potentially replacement therapy. However, traditional protocols for hiPSC-derived cardiomyocytes yield mixed populations of atrial, ventricular and nodal-like cells with immature cardiac properties. New insights gleaned from embryonic heart development have progressed the precise production of subtype-specific hiPSC-derived cardiomyocytes; however, their physiological immaturity severely limits their utility as model systems and their use for drug screening and cell therapy. The long-entrenched challenges in this field are being addressed by innovative bioengingeering technologies that incorporate biophysical, biochemical and more recently biomimetic electrical cues, with the latter having the potential to be used to both direct hiPSC differentiation and augment maturation and the function of derived cardiomyocytes and cardiac tissues by mimicking endogenous electric fields.

Network-driven discovery yields new insight into Shox2-dependent cardiac rhythm control

The homeodomain transcription factor SHOX2 is involved in the development and function of the heart's primary pacemaker, the sinoatrial node (SAN), and has been associated with cardiac conduction-related diseases such as atrial fibrillation and sinus node dysfunction. To shed light on Shox2-dependent genetic processes involved in these diseases, we established a murine embryonic stem cell (ESC) cardiac differentiation model to investigate Shox2 pathways in SAN-like cardiomyocytes. Differential RNA-seq-based expression profiling of Shox2^+/+ and Shox2^-/- ESCs revealed 94 dysregulated transcripts in Shox2^-/- ESC-derived SAN-like cells. Of these, 15 putative Shox2 target genes were selected for further validation based on comparative expression analysis with SAN- and right atria-enriched genes. Network-based analyses, integrating data from the Mouse Organogenesis Cell Atlas and the Ingenuity pathways, as well as validation in mouse and zebrafish models confirmed a regulatory role for the novel identified Shox2 target genes including Cav1, Fkbp10, Igfbp5, Mcf2l and Nr2f2. Our results indicate that genetic networks involving SHOX2 may contribute to conduction traits through the regulation of these genes.

Sinoatrial node pacemaker cells share dominant biological properties with glutamatergic neurons

Activation of the heart normally begins in the sinoatrial node (SAN). Electrical impulses spontaneously released by SAN pacemaker cells (SANPCs) trigger the contraction of the heart. However, the cellular nature of SANPCs remains controversial. Here, we report that SANPCs exhibit glutamatergic neuron-like properties. By comparing the single-cell transcriptome of SANPCs with that of cells from primary visual cortex in mouse, we found that SANPCs co-clustered with cortical neurons. Tissue and cellular imaging confirmed that SANPCs contained key elements of glutamatergic neurotransmitter system, expressing genes encoding glutamate synthesis pathway (Gls), ionotropic and metabotropic glutamate receptors (Grina, Gria3, Grm1 and Grm5), and glutamate transporters (Slc17a7). SANPCs highly expressed cell markers of glutamatergic neurons (Snap25 and Slc17a7), whereas Gad1, a marker of GABAergic neurons, was negative. Functional studies revealed that inhibition of glutamate receptors or transporters reduced spontaneous pacing frequency of isolated SAN tissues and spontaneous Ca²⁺ transients frequency in single SANPC. Collectively, our work suggests that SANPCs share dominant biological properties with glutamatergic neurons, and the glutamatergic neurotransmitter system may act as an intrinsic regulation module of heart rhythm, which provides a potential intervention target for pacemaker cell-associated arrhythmias.

Cardiac Development: A Glimpse on Its Translational Contributions

Cardiac development is a complex developmental process that is initiated soon after gastrulation, as two sets of precardiac mesodermal precursors are symmetrically located and subsequently fused at the embryonic midline forming the cardiac straight tube. Thereafter, the cardiac straight tube invariably bends to the right, configuring the first sign of morphological left–right asymmetry and soon thereafter the atrial and ventricular chambers are formed, expanded and progressively septated. As a consequence of all these morphogenetic processes, the fetal heart acquired a four-chambered structure having distinct inlet and outlet connections and a specialized conduction system capable of directing the electrical impulse within the fully formed heart. Over the last decades, our understanding of the morphogenetic, cellular, and molecular pathways involved in cardiac development has exponentially grown. Multiples aspects of the initial discoveries during heart formation has served as guiding tools to understand the etiology of cardiac congenital anomalies and adult cardiac pathology, as well as to enlighten novels approaches to heal the damaged heart. In this review we provide an overview of the complex cellular and molecular pathways driving heart morphogenesis and how those discoveries have provided new roads into the genetic, clinical and therapeutic management of the diseased hearts.

Cellular and molecular landscape of mammalian sinoatrial node revealed by single-cell RNA sequencing

Bioelectrical impulses intrinsically generated within the sinoatrial node (SAN) trigger the contraction of the heart in mammals. Though discovered over a century ago, the molecular and cellular features of the SAN that underpin its critical function in the heart are uncharted territory. Here, we identify four distinct transcriptional clusters by single-cell RNA sequencing in the mouse SAN. Functional analysis of differentially expressed genes identifies a core cell cluster enriched in the electrogenic genes. The similar cellular features are also observed in the SAN from both rabbit and cynomolgus monkey. Notably, Vsnl1, a core cell cluster marker in mouse, is abundantly expressed in SAN, but is barely detectable in atrium or ventricle, suggesting that Vsnl1 is a potential SAN marker. Importantly, deficiency of Vsnl1 not only reduces the beating rate of human induced pluripotent stem cell - derived cardiomyocytes (hiPSC-CMs) but also the heart rate of mice. Furthermore, weighted gene co-expression network analysis (WGCNA) unveiled the core gene regulation network governing the function of the SAN in mice. Overall, these findings reveal the whole transcriptome profiling of the SAN at single-cell resolution, representing an advance toward understanding of both the biology and the pathology of SAN.

The spontaneous bioelectrical activity of pacemaker cells in sinoatrial node (SAN) triggers the heartbeats. Here, the authors perform single-cell RNA sequencing in the mouse SAN and identify molecular and cellular features of the SAN conserved in rabbit and cynomolgus monkey, identifying a new potential SAN marker.

GATA6 is a regulator of sinus node development and heart rhythm

The sinus node (SAN) is the primary pacemaker of the human heart, and abnormalities in its structure or function cause sick sinus syndrome, the most common reason for electronic pacemaker implantation. Here we report that transcription factor GATA6, whose mutations in humans are linked to arrhythmia, is highly expressed in the SAN and its haploinsufficiency in mice results in hypoplastic SANs and rhythm abnormalities. Cell-specific deletion reveals a requirement for GATA6 in various SAN lineages. Mechanistically, GATA6 directly activates key regulators of the SAN genetic program in conduction and nonconduction cells, such as TBX3 and EDN1, respectively. The data identify GATA6 as an important regulator of the SAN and provide a molecular basis for understanding the conduction abnormalities associated with GATA6 mutations in humans. They also suggest that GATA6 may be a potential modifier of the cardiac pacemaker.

Genome-Wide Analysis Identifies an Essential Human TBX3 Pacemaker Enhancer

RATIONALE

The development and function of the pacemaker cardiomyocytes of the sinoatrial node (SAN), the leading pacemaker of the heart, are tightly controlled by a conserved network of transcription factors, including TBX3 (T-box transcription factor 3), ISL1 (ISL LIM homeobox 1), and SHOX2 (short stature homeobox 2). Yet, the regulatory DNA elements (REs) controlling target gene expression in the SAN pacemaker cells have remained undefined.

OBJECTIVE

Identification of the regulatory landscape of human SAN-like pacemaker cells and functional assessment of SAN-specific REs potentially involved in pacemaker cell gene regulation.

METHODS AND RESULTS

We performed Assay for Transposase-Accessible Chromatin using sequencing on human pluripotent stem cell-derived SAN-like pacemaker cells and ventricle-like cells and identified thousands of putative REs specific for either human cell type. We validated pacemaker cell-specific elements in the SHOX2 and TBX3 loci. CRISPR-mediated homozygous deletion of the mouse ortholog of a noncoding region with candidate pacemaker-specific REs in the SHOX2 locus resulted in selective loss of Shox2 expression from the developing SAN and embryonic lethality. Putative pacemaker-specific REs were identified up to 1 Mbp upstream of TBX3 in a region close to MED13L harboring variants associated with heart rate recovery after exercise. The orthologous region was deleted in mice, which resulted in selective loss of expression of Tbx3 from the SAN and (cardiac) ganglia and in neonatal lethality. Expression of Tbx3 was maintained in other tissues including the atrioventricular conduction system, lungs, and liver. Heterozygous adult mice showed increased SAN recovery times after pacing. The human REs harboring the associated variants robustly drove expression in the SAN of transgenic mouse embryos.

CONCLUSIONS

We provided a genome-wide collection of candidate human pacemaker-specific REs, including the loci of SHOX2, TBX3, and ISL1, and identified a link between human genetic variants influencing heart rate recovery after exercise and a variant RE with highly conserved function, driving SAN expression of TBX3.

Development of the Cardiac Conduction System

The cardiac conduction system initiates and propagates each heartbeat. Specialized conducting cells are a well-conserved phenomenon across vertebrate evolution, although mammalian and avian species harbor specific components unique to organisms with four-chamber hearts. Early histological studies in mammals provided evidence for a dominant pacemaker within the right atrium and clarified the existence of the specialized muscular axis responsible for atrioventricular conduction. Building on these seminal observations, contemporary genetic techniques in a multitude of model organisms has characterized the developmental ontogeny, gene regulatory networks, and functional importance of individual anatomical compartments within the cardiac conduction system. This review describes in detail the transcriptional and regulatory networks that act during cardiac conduction system development and homeostasis with a particular emphasis on networks implicated in human electrical variation by large genome-wide association studies. We conclude with a discussion of the clinical implications of these studies and describe some future directions.

Evidence of Superior and Inferior Sinoatrial Nodes in the Mammalian Heart

OBJECTIVES

This study sought to investigate the shift of leading pacemaker locations in healthy and failing mammalian hearts over the entire range of physiological heart rates (HRs), and to molecularly characterize spatial regions of spontaneous activity.

BACKGROUND

A normal heartbeat originates as an action potential in a group of pacemaker cells known as the sinoatrial node (SAN), located near the superior vena cava. HRs and the anatomical site of origin of pacemaker activity in the adult heart are known to dynamically change in response to various physiological inputs, yet the mechanism of this pacemaker shift is not well understood.

METHODS

Optical mapping was applied to ex vivo rat and human isolated right atrial tissues, and HRs were modulated with acetylcholine and isoproterenol. RNA sequencing was performed on tissue areas that elicited spontaneous activity, and comparisons were made to neighboring myocardial tissues.

RESULTS

Functional and molecular evidence identified and confirmed the presence of 2 competing right atrial pacemakers localized near the superior vena cava and the inferior vena cava—the superior SAN (sSAN) and inferior SAN (iSAN), respectively—which preferentially control the fast and slow HRs. Both of these regions were evident in non-failing rat and human hearts and maintained spontaneous activity in the rat heart when physically separated from one another. Molecular analysis of these 2 pacemaker regions revealed unique but similar transcriptional profiles, suggesting iSAN dominance when the sSAN is silent.

CONCLUSIONS

The presence of 2 spatially distinct dominant pacemakers, sSAN and iSAN, in the mammalian heart clarifies previous identification of migrating pacemakers and corresponding changes in P-wave morphology in mammalian species.

Nkx2-5 defines distinct scaffold and recruitment phases during formation of the murine cardiac Purkinje fiber network

The ventricular conduction system coordinates heartbeats by rapid propagation of electrical activity through the Purkinje fiber (PF) network. PFs share common progenitors with contractile cardiomyocytes, yet the mechanisms of segregation and network morphogenesis are poorly understood. Here, we apply genetic fate mapping and temporal clonal analysis to identify murine cardiomyocytes committed to the PF lineage as early as E7.5. We find that a polyclonal PF network emerges by progressive recruitment of conductive precursors to this scaffold from a pool of bipotent progenitors. At late fetal stages, the segregation of conductive cells increases during a phase of rapid recruitment to build the definitive PF network through a non-cell autonomous mechanism. We also show that PF differentiation is impaired in Nkx2-5 haploinsufficient embryos leading to failure to extend the scaffold. In particular, late fetal recruitment fails, resulting in PF hypoplasia and persistence of bipotent progenitors. Our results identify how transcription factor dosage regulates cell fate divergence during distinct phases of PF network morphogenesis.

Here, the authors apply genetic fate mapping and temporal clonal analysis to study progenitor recruitment and network morphogenesis of murine cardiac Purkinje fibers. Additionally, they characterize how transcription factor dosage regulates cell fate divergence during distinct phases of this process.