Genetic Regulation of Sinoatrial Node Development and Pacemaker Program in the Venous Pole

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.

Segregation of morphogenetic regulatory function of Shox2 from its cell fate guardian role in sinoatrial node development

Shox2 plays a vital role in the morphogenesis and physiological function of the sinoatrial node (SAN), the primary cardiac pacemaker, manifested by the formation of a hypoplastic SAN and failed differentiation of pacemaker cells in Shox2 mutants. Shox2 and Nkx2-5 are co-expressed in the developing SAN and regulate the fate of the pacemaker cells through a Shox2-Nkx2-5 antagonistic mechanism. Here we show that simultaneous inactivation of Nkx2-5 in the SAN of Shox2 mutants (dKO) rescued the pacemaking cell fate but not the hypoplastic defects, indicating uncoupling of SAN cell fate determination and morphogenesis. Single-cell RNA-seq revealed that the presumptive SAN cells of Shox2-/- mutants failed to activate pacemaking program but remained in a progenitor state preceding working myocardium, while both wildtype and dKO SAN cells displayed normal pacemaking cell fate with similar cellular state. Shox2 thus acts as a safeguard but not a determinant to ensure the pacemaking cell fate through the Shox2-Nkx2-5 antagonistic mechanism, which is segregated from its morphogenetic regulatory function in SAN development.

Local tissue mechanics control cardiac pacemaker cell embryonic patterning

Cardiac pacemaker cells (CPCs) initiate the electric impulses that drive the rhythmic beating of the heart. CPCs reside in a heterogeneous, ECM-rich microenvironment termed the sinoatrial node (SAN). Surprisingly, little is known regarding the biochemical composition or mechanical properties of the SAN, and how the unique structural characteristics present in this region of the heart influence CPC function remains poorly understood. Here, we have identified that SAN development involves the construction of a "soft" macromolecular ECM that specifically encapsulates CPCs. In addition, we demonstrate that subjecting embryonic CPCs to substrate stiffnesses higher than those measured in vivo results in loss of coherent electrical oscillation and dysregulation of the HCN4 and NCX1 ion channels required for CPC automaticity. Collectively, these data indicate that local mechanics play a critical role in maintaining the embryonic CPC function while also quantitatively defining the range of material properties that are optimal for embryonic CPC maturation.

Decoding the PITX2-controlled genetic network in atrial fibrillation

Atrial fibrillation (AF), the most common sustained cardiac arrhythmia and a major risk factor for stroke, often arises through ectopic electrical impulses derived from the pulmonary veins (PVs). Sequence variants in enhancers controlling expression of the transcription factor PITX2, which is expressed in the cardiomyocytes (CMs) of the PV and left atrium (LA), have been implicated in AF predisposition. Single nuclei multiomic profiling of RNA and analysis of chromatin accessibility combined with spectral clustering uncovered distinct PV- and LA-enriched CM cell states. Pitx2-mutant PV and LA CMs exhibited gene expression changes consistent with cardiac dysfunction through cell type-distinct, PITX2-directed, cis-regulatory grammars controlling target gene expression. The perturbed network targets in each CM were enriched in distinct human AF predisposition genes, suggesting combinatorial risk for AF genesis. Our data further reveal that PV and LA Pitx2-mutant CMs signal to endothelial and endocardial cells through BMP10 signaling with pathogenic potential. This work provides a multiomic framework for interrogating the basis of AF predisposition in the PVs of humans.

Altered microRNA and mRNA profiles during heart failure in the human sinoatrial node

Heart failure (HF) is frequently accompanied with the sinoatrial node (SAN) dysfunction, which causes tachy-brady arrhythmias and increased mortality. MicroRNA (miR) alterations are associated with HF progression. However, the transcriptome of HF human SAN, and its role in HF-associated remodeling of ion channels, transporters, and receptors responsible for SAN automaticity and conduction impairments is unknown. We conducted comprehensive high-throughput transcriptomic analysis of pure human SAN primary pacemaker tissue and neighboring right atrial tissue from human transplanted HF hearts (n = 10) and non-failing (nHF) donor hearts (n = 9), using next-generation sequencing. Overall, 47 miRs and 832 mRNAs related to multiple signaling pathways, including cardiac diseases, tachy-brady arrhythmias and fibrosis, were significantly altered in HF SAN. Of the altered miRs, 27 are predicted to regulate mRNAs of major ion channels and neurotransmitter receptors which are involved in SAN automaticity (e.g. HCN1, HCN4, SLC8A1) and intranodal conduction (e.g. SCN5A, SCN8A) or both (e.g. KCNJ3, KCNJ5). Luciferase reporter assays were used to validate interactions of miRs with predicted mRNA targets. In conclusion, our study provides a profile of altered miRs in HF human SAN, and a novel transcriptome blueprint to identify molecular targets for SAN dysfunction and arrhythmia treatments in HF.

PDGFRα-Signaling Is Dispensable for the Development of the Sinoatrial Node After Its Fate Commitment

Palate-derived growth factor receptor α (Pdgfrα) signaling has been reported to play important roles in the cardiac development. A previous study utilizing Pdgfrα conventional knockout mice reported hypoplasia of the sinus venous myocardium including the sinoatrial node (SAN) accompanied by increased expression of Nkx2.5. This mouse line embryos die by E11.5 due to embryonic lethality, rendering them difficult to investigate the details. To elucidate the underlying mechanism, in this study, we revisited this observation by generation of specific ablation of Pdgfrα in the SAN by Shox2-Cre at E9.5, using a Shox2-Cre;Pdgfrα ^flox/flox conditional mouse line. Surprisingly, we found that resultant homozygous mutant mice did not exhibit any malformation in SAN morphology as compared to their wild-type littermates. Further analysis revealed the normal cardiac function in adult mutant mice assessed by the record of heart rate and electrocardiogram and unaltered expression of Nkx2.5 in the E13.5 SAN of Pdgfrα conditional knockout mice. Our results unambiguously demonstrate that Pdgfrα is dispensable for SAN development after its fate commitment in mice.

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.

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.

Comparative study of hyperpolarization-activated currents in pulmonary vein cardiomyocytes isolated from rat, guinea pig, and rabbit

Pulmonary vein (PV) cardiomyocytes have the potential to generate spontaneous activity, in contrast to working myocytes of atria. Different electrophysiological properties underlie the potential automaticity of PV cardiomyocytes, one being the hyperpolarization-activated inward current (I_h), which facilitates the slow diastolic depolarization. In the present study, we examined pharmacological characteristics of the I_h of PV cardiomyocytes in rat, guinea pig and rabbit. The results showed that guinea pig and rat PV cardiomyocytes possessed sizeable amplitudes of the I_h, and the I_h of guinea pig was suppressed by Cs⁺, a blocker of the hyperpolarization-activated cation current. However, the I_h of rat was not suppressed by Cs⁺, but by Cd²⁺, a blocker of the Cl⁻ current. The current density of the I_h of rabbit PV cardiomyocytes was significantly smaller than those of other species. This suggests that the ion channels that carry the I_h of PV cardiomyocytes differ among the animal species.

Association of Dysplastic Coronoid Process with Long-Face Morphology

Vertical malocclusion is a developmental condition, resulting from complex interactions among multiple etiological factors during the growth period. As a tricky dentofacial deformity clinically, long-face (LF) morphology is characterized by excessive vertical facial growth with severe disarrangement of jaws and teeth. Since the improvement of LF patients on facial profile and occlusion is often difficult and lacks long-term stability, it becomes important to unravel the etiology of LF pattern formation for early prevention and treatment. In the current studies, we identified a transgenic mouse model that exhibited a dysplastic coronoid process and LF morphology. Although the mutant mice exhibited jaw structures and occlusion comparable to controls at birth, they all acquired typical LF morphology with anterior open bite during postnatal growth, resembling clinical features of the selected skeletal class III patients. Since the coronoid process provides an insertion site for the temporalis attachment, we examined the initial development and differentiation of the temporalis and found identical results in both control and mutant mice before E17.5 when the temporal muscle makes attachment to the coronoid process. However, thereafter, we observed altered orientation and reduced size of the cross-sectional area of the temporalis in mutant mice, which persisted to the weaning stage. Biomechanical analysis and simulation modeling further support the idea that altered morphology of the coronoid process may impair the efficiency of the vertical temporalis contraction and appears to correlate with LF formation. Consistently, we present evidence that a dysplastic mandibular coronoid process was also seen in some human patients with skeletal III LF morphology. Taken together, the results presented in this study establish an association of the craniofacial bony structures with vertical patterning, which will have implications in earlier prediction for clinical precaution and intervention.

Differences in Expression of Genes Involved in Bone Development and Morphogenesis in the Walls of Internal Thoracic Artery and Saphenous Vein Conduits May Provide Markers Useful for Evaluation Graft Patency

Coronary artery bypass grafting (CABG) is one of the most efficient procedures for patients with advanced coronary artery disease. From all the blood vessels with the potential to be used in this procedure, the internal thoracic artery (ITA) and the saphenous vein (SV) are the most commonly applied as aortocoronary conduits. Nevertheless, in order to evaluate the graft patency and efficiency effectively, basic knowledge should be constantly expanding at the molecular level as well, as the understanding of predictive factors is still limited. In this study, we have employed the expressive microarray approach, validated with Real-Time Quantitative Polymerase Chain Reaction (RT-qPCR), to analyze the transcriptome of both venous and arterial grafts. Searching for potential molecular factors, we analyzed differentially expressed gene ontologies involved in bone development and morphogenesis, for the possibility of discovery of new markers for the evaluation of ITA and SV segment quality. Among three ontological groups of interest—“endochondral bone morphogenesis”, “ossification”, and “skeletal system development”—we found six genes common to all of them. BMP6, SHOX2, COL13A1, CSGALNACT1, RUNX2, and STC1 showed differential expression patterns in both analyzed vessels. STC1 and COL13A1 were upregulated in ITA samples, whereas others were upregulated in SV. With regard to the Runx2 protein function in osteogenic phenotype regulation, the RUNX2 gene seems to be of paramount importance in assessing the potential of ITA, SV, and other vessels used in the CABG procedure. Overall, the presented study provided valuable insight into the molecular background of conduit characterization, and thus indicated genes that may be the target of subsequent studies, also at the protein level. Moreover, it has been suggested that RUNX2 may be recognized as a molecular marker of osteogenic changes in human blood vessels.

Nkx2-5 defines a subpopulation of pacemaker cells and is essential for the physiological function of the sinoatrial node in mice

The sinoatrial node (SAN), the primary cardiac pacemaker, consists of a head domain and a junction/tail domain that exhibit different functional properties. However, the underlying molecular mechanism defining these two pacemaker domains remains elusive. Nkx2-5 is a key transcription factor essential for the formation of the working myocardium, but it was generally thought to be detrimental to SAN development. However, Nkx2-5 is expressed in the developing SAN junction, suggesting a role for Nkx2-5 in SAN junction development and function. In this study, we present unambiguous evidence that SAN junction cells exhibit unique action potential configurations intermediate to those manifested by the SAN head and the surrounding atrial cells, suggesting a specific role for the junction cells in impulse generation and in SAN-atrial exit conduction. Single-cell RNA-seq analyses support this concept. Although Nkx2-5 inactivation in the SAN junction did not cause a malformed SAN at birth, the mutant mice manifested sinus node dysfunction. Thus, Nkx2-5 defines a population of pacemaker cells in the transitional zone. Despite Nkx2-5 being dispensable for SAN morphogenesis during embryogenesis, its deletion hampers atrial activation by the pacemaker.

Functional Characterization of Rare Variants in the SHOX2 Gene Identified in Sinus Node Dysfunction and Atrial Fibrillation

Sinus node dysfunction (SND) and atrial fibrillation (AF) often coexist; however, the molecular mechanisms linking both conditions remain elusive. Mutations in the homeobox-containing SHOX2 gene have been recently associated with early-onset and familial AF. Shox2 is a key regulator of sinus node development, and its deficiency leads to bradycardia, as demonstrated in animal models. To provide an extended SHOX2 gene analysis in patients with distinct arrhythmias, we investigated SHOX2 as a susceptibility gene for SND and AF by screening 98 SND patients and 450 individuals with AF. The functional relevance of the novel mutations was investigated in vivo and in vitro, together with the previously reported p.H283Q variant. A heterozygous missense mutation (p.P33R) was identified in the SND cohort and four heterozygous variants (p.G77D, p.L129=, p.L130F, p.A293=) in the AF cohort. Overexpression of the pathogenic predicted mutations in zebrafish revealed pericardial edema for p.G77D and the positive control p.H283Q, whereas the p.P33R and p.A293= variants showed no effect. In addition, a dominant-negative effect with reduced heart rates was detected for p.G77D and p.H283Q. In vitro reporter assays demonstrated for both missense variants p.P33R and p.G77D significantly impaired transactivation activity, similar to the described p.H283Q variant. Also, a reduced Bmp4 target gene expression was revealed in zebrafish hearts upon overexpression of the p.P33R mutant. This study associates additional rare variants in the SHOX2 gene implicated in the susceptibility to distinct arrhythmias and allows frequency estimations in the AF cohort (3/990). We also demonstrate for the first time a genetic link between SND and AF involving SHOX2. Moreover, our data highlight the importance of functional investigations of rare variants.

Shox2: The Role in Differentiation and Development of Cardiac Conduction System

The formation and conduction of electrocardiosignals and the synchronous contraction of atria and ventricles with rhythmicity are both triggered and regulated by the cardiac conduction system (CCS). Defect of this system will lead to various types of cardiac arrhythmias. In recent years, the research progress of molecular genetics and developmental biology revealed a clearer understanding of differentiation and development of the cardiac conduction system and their regulatory mechanisms. Short stature homeobox 2 (Shox2) transcription factor, encoded by Shox2 gene in the mouse, is crucial in the formation and differentiation of the sinoatrial node (SAN). Shox2 drives embryonic development processes and is widely expressed in the appendicular skeleton, palate, temporomandibular joints, and heart. Mutations of Shox2 can lead to dysembryoplasia and abnormal phenotypes, including bradycardiac arrhythmia. In this review, we provide a summary of the latest research progress on the regulatory effects of the Shox2 gene in differentiation and development processes of the cardiac conduction system, hoping to deepen the knowledge and understanding of this systematic process based on the cardiac conduction system. Overall, the Shox2 gene is intimately involved in the differentiation and development of cardiac conduction system, especially sinoatrial node. We also summarize the current information about human SHOX2. This review article provides a new direction in biological pacemaker therapies.

A SHOX2 loss-of-function mutation underlying familial atrial fibrillation

Atrial fibrillation (AF), as the most common sustained cardiac arrhythmia, is associated with substantially increased morbidity and mortality. Aggregating evidence demonstrates that genetic defects play a crucial role in the pathogenesis of AF, especially in familial AF. Nevertheless, AF is of pronounced genetic heterogeneity, and in an overwhelming majority of cases the genetic determinants underlying AF remain elusive. In the current study, 162 unrelated patients with familial AF and 238 unrelated healthy individuals served as controls were recruited. The coding exons and splicing junction sites of the SHOX2 gene, which encodes a homeobox-containing transcription factor essential for proper development and function of the cardiac conduction system, were sequenced in all study participants. The functional effect of the mutant SHOX2 protein was characterized with a dual-luciferase reporter assay system. As a result, a novel heterozygous SHOX2 mutation, c.580C>T or p.R194X, was identified in an index patient, which was absent from the 476 control chromosomes. Genetic analysis of the proband's pedigree revealed that the nonsense mutation co-segregated with AF in the family with complete penetrance. Functional assays demonstrated that the mutant SHOX2 protein had no transcriptional activity compared with its wild-type counterpart. In conclusion, this is the first report on the association of SHOX2 loss-of-function mutation with enhanced susceptibility to familial AF, which provides novel insight into the molecular mechanism underpinning AF, suggesting potential implications for genetic counseling and individualized management of AF patients.

Development of the cardiac pacemaker

The sinoatrial node (SAN) is the dominant pacemaker of the heart. Abnormalities in SAN formation and function can cause sinus arrhythmia, including sick sinus syndrome and sudden death. A better understanding of genes and signaling pathways that regulate SAN development and function is essential to develop more effective treatment to sinus arrhythmia, including biological pacemakers. In this review, we briefly summarize the key processes of SAN morphogenesis during development, and focus on the transcriptional network that drives SAN development.

A unique stylopod patterning mechanism by Shox2-controlled osteogenesis

Vertebrate appendage patterning is programmed by Hox-TALE factor-bound regulatory elements. However, it remains unclear which cell lineages are commissioned by Hox-TALE factors to generate regional specific patterns and whether other Hox-TALE co-factors exist. In this study, we investigated the transcriptional mechanisms controlled by the Shox2 transcriptional regulator in limb patterning. Harnessing an osteogenic lineage-specific Shox2 inactivation approach we show that despite widespread Shox2 expression in multiple cell lineages, lack of the stylopod observed upon Shox2 deficiency is a specific result of Shox2 loss of function in the osteogenic lineage. ChIP-Seq revealed robust interaction of Shox2 with cis-regulatory enhancers clustering around skeletogenic genes that are also bound by Hox-TALE factors, supporting a lineage autonomous function of Shox2 in osteogenic lineage fate determination and skeleton patterning. Pbx ChIP-Seq further allowed the genome-wide identification of cis-regulatory modules exhibiting co-occupancy of Pbx, Meis and Shox2 transcriptional regulators. Integrative analysis of ChIP-Seq and RNA-Seq data and transgenic enhancer assays indicate that Shox2 patterns the stylopod as a repressor via interaction with enhancers active in the proximal limb mesenchyme and antagonizes the repressive function of TALE factors in osteogenesis.

Coding and non-coding variants in the SHOX2 gene in patients with early-onset atrial fibrillation

Atrial fibrillation (AF) is the most prevalent cardiac arrhythmia with a strong genetic component. Molecular pathways involving the homeodomain transcription factor Shox2 control the development and function of the cardiac conduction system in mouse and zebrafish. Here we report the analysis of human SHOX2 as a potential susceptibility gene for early-onset AF. To identify causal variants and define the underlying mechanisms, results from 378 patients with early-onset AF before the age of 60 years were analyzed and compared to 1870 controls or reference datasets. We identified two missense mutations (p.G81E, p.H283Q), that were predicted as damaging. Transactivation studies using SHOX2 targets and phenotypic rescue experiments in zebrafish demonstrated that the p.H283Q mutation severely affects SHOX2 pacemaker function. We also demonstrate an association between a 3'UTR variant c.*28T>C of SHOX2 and AF (p = 0.00515). Patients carrying this variant present significantly longer PR intervals. Mechanistically, this variant creates a functional binding site for hsa-miR-92b-5p. Circulating hsa-miR-92b-5p plasma levels were significantly altered in AF patients carrying the 3'UTR variant (p = 0.0095). Finally, we demonstrate significantly reduced SHOX2 expression levels in right atrial appendages of AF patients compared to patients with sinus rhythm. Together, these results suggest a genetic contribution of SHOX2 in early-onset AF.