A common Shox2-Nkx2-5 antagonistic mechanism primes the pacemaker cell fate in the pulmonary vein myocardium and sinoatrial node

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.

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.

Development of ventricular trabeculae affects electrical conduction in the early endothermic heart

BACKGROUND

The ventricular trabeculae play a role, among others, in the impulse spreading in ectothermic hearts. Despite the morphological similarity with the early developing hearts of endotherms, this trabecular function in mammalian and avian embryos was poorly addressed.

RESULTS

We simulated impulse propagation inside the looping ventricle and revealed delayed apical activation in the heart with inhibited trabecular growth. This finding was corroborated by direct imaging of the endocardial surface showing early activation within the trabeculae implying preferential spreading of depolarization along with them. Targeting two crucial pathways of trabecular formation (Neuregulin/ErbB and Nkx2.5), we showed that trabecular development is also essential for proper conduction patterning. Persistence of the slow isotropic conduction likely contributed to the pumping failure in the trabeculae-deficient hearts.

CONCLUSIONS

Our results showed the essential role of trabeculae in intraventricular impulse spreading and conduction patterning in the early endothermic heart. Lack of trabeculae leads to the failure of conduction parameters differentiation resulting in primitive ventricular activation with consequent impact on the cardiac pumping function.

Development of the Cardiac Conduction System

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

The method of sinus node-like pacemaker cells from human induced pluripotent stem cells by BMP and Wnt signaling

The embryonic development of sinus nodes (SAN) is co-regulated by multiple signaling pathways. Among these, the bone morphogenetic protein (BMP) and Wnt signaling pathways are involved in the development of SAN. In this study, the effects of BMP and Wnt signaling on the differentiation of SAN-like pacemaker cells (SANLPCs) were investigated. Human induced pluripotent stem cells (hiPSCs) were divided into four groups: control, BMP4, CHIR-3, and BMP4 + CHIR (CHIR: a Wnt signaling activator). The samples were tested at day (D) 15 of differentiation. The final protocol for the activation of BMP signaling at D0-D3 and reactivation of Wnt signaling at D5-D7 in the differentiation of hiPSCs were determined. The results showed that the mRNA levels of pacemaker markers (TBX18, SHOX2, TBX3, HCN4, and HCN1) were higher in the BMP4 + CHIR group than in the control group, and working myocardial genes were downregulated. The immunofluorescence assay revealed that the expression of SHOX2 and HCN4 increased in the BMP4 + CHIR group compared to that in the other groups. In addition, the results of patch clamps revealed that a funny current of higher density and typical SAN action potentials were recorded, except in the control group, in which the L-type calcium current was higher in the BMP4 + CHIR group than in the other groups. Finally, the proportion of SANLPCs (cTnT+ NKX2.5-) was further enhanced by the combination of BMP4 and CHIR treatment. In summary, the combination of BMP and Wnt signaling promotes the differentiation of SANLPCs from hiPSCs.

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.

Trends in research on sick sinus syndrome: A bibliometric analysis from 2000 to 2022

Sick sinus syndrome (SSS) is a refractory arrhythmia disease caused by the pathological changes of sinoatrial node and its adjacent tissues. 2,251 publications related to SSS were retrieved from Web of Science database from 2000 to 2022 and analyzed by using VOS viewer and CiteSpace software. The results showed the United States dominated the field, followed by Japan, Germany, and China. SSS was closely related to risk factors such as atrial fibrillation and aging. Sick sinus syndrome, atrial fibrillation and sinus node dysfunction were the top three keywords that had the strongest correlation with the study. Pacemaker implantation, differentiation and mutation are research hotspots currently. Clinical studies on SSS found that sick sinus syndrome, atrial fibrillation, and pacemakers were the top three keywords that had the largest nodes and the highest frequency. In the field of basic applied research and basic research, atrial fibrillation and pacemaker cells were the focus of research. In conclusion, bibliometric analysis provided valuable information for the prevention, treatment and future research trends of SSS.

α1-adrenergic receptors accompanied by GATA4 expression are related to proarrhythmic conduction and automaticity in rat interatrial septum

The development of interatrial septum (IAS) is a complicated process, which continues during postnatal life. The hypertrophic signals in developing heart are mediated among others by α-adrenergic pathways. These facts suggest the presence of specific electrophysiological features in developing IAS. This study was aimed to investigate the electrical activity in the tissue preparations of IAS from rat heart in normal conditions and under stimulation of adrenoreceptors. Intracellular recording of electrical activity revealed less negative level of resting membrane potential in IAS if compared to myocardium of left atrium. In normal conditions, non-paced IAS preparations were quiescent, but noradrenaline (10^-5 M) and phenylephrine (10^-5 M) induced spontaneous action potentials, which could be abolished by α1-blocker prazosin (10^-5 M), but not β1-blocker atenolol (10^-5 M). Optical mapping showed drastic phenylephrine-induced slowing of conduction in adult rat IAS. The α1-dependent ectopic automaticity of IAS myocardium might be explained by immunohistochemical data indicating the presence of transcription factor GATA4 and abundant α1A-adrenoreceptors in myocytes from adult rat IAS. An elevated sensitivity to adrenergic stimulation due to involvement of α1-adrenergic pathways may underlie increased proarrhythmic potential of adult IAS at least in rats.

Chemically defined and small molecules-based generation of sinoatrial node-like cells

Background

Existing methods for in vitro differentiation of human pluripotent stem cells (hPSCs) into sinoatrial node-like cells (SANLCs) require complex and undefined medium constituents. This might hinder the elucidation of the molecular mechanisms involved in cardiac subtype specification and prevent translational application. In our study, we aimed to establish a chemically defined differentiation methods to generate SANLCs effectively and stably.

Methods

We induced human embryonic stem cells (hESCs)/induced PSCs (hiPSCs) to pan-cardiomyocytes by temporal modulation of the WNT/β-catenin (WNT) signaling pathway with GSK3 inhibitor and WNT inhibitor. During cardiac mesoderm stage of the differentiation process, signaling of WNT, retinoid acid (RA), and fibroblast growth factor (FGF) was manipulated by three specific molecules. Moreover, metabolic selection was designed to improve the enrichment of SANLCs. Finally, RT-PCR, immunofluorescence, flow cytometry, and whole cell patch clamp were used to identify the SANLCs.

Results

WNT, RA, and FGF signaling promote the differentiation of hPSCs into SANLCs in a concentration- and time window-sensitive manner, respectively. Synergetic modulation of WNT, FGF, and RA signaling pathways enhance the pacemaker phenotype and improve the differentiation efficiency of SANLCs (up to 45%). Moreover, the purification based on lactate metabolism and glucose starvation further reached approximately 50% of SANLCs. Finally, the electrophysiological data demonstrate that cells differentiated with the proposed protocol produce a considerable number of SANLCs that display typical electrophysiological characteristics of pacemaker cells in vitro.

Conclusion

We provide an optimized and chemically defined protocol to generate SANLCs by combined modulation of WNT, RA, and FGF signaling pathways and metabolic selection by lactate enrichment and glucose starvation. This chemically defined method for generating SANLCs might provide a platform for disease modeling, drug discovery, predictive toxicology, and biological pacemaker construction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-022-02834-y.

A dual SHOX2:GFP; MYH6:mCherry knockin hESC reporter line for derivation of human SAN-like cells

The sinoatrial node (SAN) is the primary pacemaker of the heart. The human SAN is poorly understood due to limited primary tissue access and limitations in robust in vitro derivation methods. We developed a dual SHOX2:GFP; MYH6:mCherry knockin human embryonic stem cell (hESC) reporter line, which allows the identification and purification of SAN-like cells. Using this line, we performed several rounds of chemical screens and developed an efficient strategy to generate and purify hESC-derived SAN-like cells (hESC-SAN). The derived hESC-SAN cells display molecular and electrophysiological characteristics of bona fide nodal cells, which allowed exploration of their transcriptional profile at single-cell level. In sum, our dual reporter system facilitated an effective strategy for deriving human SAN-like cells, which can potentially be used for future disease modeling and drug discovery.

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

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

Implementing Biological Pacemakers: Design Criteria for Successful

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

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

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

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.

Adherens junction engagement regulates functional patterning of the cardiac pacemaker cell lineage

Cardiac pacemaker cells (CPCs) rhythmically initiate the electrical impulses that drive heart contraction. CPCs display the highest rate of spontaneous depolarization in the heart despite being subjected to inhibitory electrochemical conditions that should theoretically suppress their activity. While several models have been proposed to explain this apparent paradox, the actual molecular mechanisms that allow CPCs to overcome electrogenic barriers to their function remain poorly understood. Here, we have traced CPC development at single-cell resolution and uncovered a series of cytoarchitectural patterning events that are critical for proper pacemaking. Specifically, our data reveal that CPCs dynamically modulate adherens junction (AJ) engagement to control characteristics including surface area, volume, and gap junctional coupling. This allows CPCs to adopt a structural configuration that supports their overall excitability. Thus, our data have identified a direct role for local cellular mechanics in patterning critical morphological features that are necessary for CPC electrical activity.

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.

The susceptibility of the rat pulmonary and caval vein myocardium to the catecholamine-induced ectopy changes oppositely in postnatal development

KEY POINTS

The developmental changes of the caval (SVC) and pulmonary vein (PV) myocardium electrophysiology are traced throughout postnatal ontogenesis. The myocardium in SVC as well as in PV demonstrate age-dependent differences in the ability to maintain resting membrane potential, to manifest automaticity in a form of ectopic action potentials in basal condition and in responses to the adrenergic stimulation. Electrophysiological characteristics of two distinct types of thoracic vein myocardium change in an opposite manner during early postnatal ontogenesis with increased proarrhythmicity of pulmonary and decreased automaticity in caval veins. Predisposition of PV cardiac tissue to proarrhythmycity develops during ontogenesis in time correlation with the establishment of sympathetic innervation of the tissue. The electrophysiological properties of caval vein cardiac tissue shift from a pacemaker-like phenotype to atrial phenotype in accompaniment with sympathetic nerve growth and adrenergic receptor expression changes.

ABSTRACT

The thoracic vein myocardium is considered as a main source for atrial fibrillation initiation due to its high susceptibility to ectopic activity. The mechanism by which and when pulmonary (PV) and superior vena cava (SVC) became proarrhythmic during postnatal ontogenesis is still unknown. In this study, we traced postnatal changes of electrophysiology in a correlation with the sympathetic innervation and adrenergic receptor distribution to reveal developmental differences in proarrhythmicity occurrence in PV and SVC myocardium. A standard microelectrode technique was used to assess the changes in ability to maintain resting membrane potential (RMP), generate spontaneous action potentials (SAP) and adrenergically induced ectopy in multicellular SVC and PV preparations of rats of different postnatal ages. Immunofluorescence imaging was used to trace postnatal changes in sympathetic innervation, β1- and α1A-adrenergic receptor (AR) distribution. We revealed that the ability to generate SAP and susceptibility to adrenergic stimulation changes during postnatal ontogenesis in an opposite manner in PV and SVC myocardium. While SAP occurrence decreases with age in SVC myocardium, it significantly increases in PV cardiac tissue. PV myocardium starts to demonstrate RMP instability and proarrhythmic activity from the 14th day of postnatal life which correlates with the appearance of the sympathetic innervation of the thoracic veins. In addition, postnatal attenuation of SVC myocardium automaticity occurs concomitantly with sympathetic innervation establishment and increase in β1-ARs, but not α1A-AR levels. Our results support the contention that SVC and PV myocardium electrophysiology change during postnatal development, resulting in higher PV proarrhythmicity in adults.

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.

The Transcription Factor Shox2 Shapes Neuron Firing Properties and Suppresses Seizures by Regulation of Key Ion Channels in Thalamocortical Neurons

Thalamocortical neurons (TCNs) play a critical role in the maintenance of thalamocortical oscillations, dysregulation of which can result in certain types of seizures. Precise control over firing rates of TCNs is foundational to these oscillations, yet the transcriptional mechanisms that constrain these firing rates remain elusive. We hypothesized that Shox2 is a transcriptional regulator of ion channels important for TCN function and that loss of Shox2 alters firing frequency and activity, ultimately perturbing thalamocortical oscillations into an epilepsy-prone state. In this study, we used RNA sequencing and quantitative PCR of control and Shox2 knockout mice to determine Shox2-affected genes and revealed a network of ion channel genes important for neuronal firing properties. Protein regulation was confirmed by Western blotting, and electrophysiological recordings showed that Shox2 KO impacted the firing properties of a subpopulation of TCNs. Computational modeling showed that disruption of these conductances in a manner similar to Shox2's effects modulated frequency of oscillations and could convert sleep spindles to near spike and wave activity, which are a hallmark for absence epilepsy. Finally, Shox2 KO mice were more susceptible to pilocarpine-induced seizures. Overall, these results reveal Shox2 as a transcription factor important for TCN function in adult mouse thalamus.

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.

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.

Precise Correction of Heterozygous SHOX2 Mutations in hiPSCs Derived from Patients with Atrial Fibrillation via Genome Editing and Sib Selection

Patient-specific human induced pluripotent stem cells (hiPSCs) offer unprecedented opportunities for the investigation of multigenic disease, personalized medicine, and stem cell therapy. For heterogeneous diseases such as atrial fibrillation (AF), however, precise correction of the associated mutation is crucial. Here, we generated and corrected hiPSC lines from two AF patients carrying different heterozygous SHOX2 mutations. We developed a strategy for the scarless correction of heterozygous mutations, based on stochastic enrichment by sib selection, followed by allele quantification via digital PCR and next-generation sequencing to detect isogenic subpopulations. This allowed enriching edited cells 8- to 20-fold. The method does not require antibiotic selection or cell sorting and can be easily combined with base-and-prime editing approaches. Our strategy helps to overcome low efficiencies of homology-dependent repair in hiPSCs and facilitates the generation of isogenic control lines that represent the gold standard for modeling complex diseases in vitro.

NODAL inhibition promotes differentiation of pacemaker-like cardiomyocytes from human induced pluripotent stem cells

Directed cardiomyogenesis from human induced pluripotent stem cells (hiPSCs) has been greatly improved in the last decade but directed differentiation to pacemaking cardiomyocytes (CMs) remains incompletely understood. In this study, we demonstrated that inhibition of NODAL signaling by a specific NODAL inhibitor (SB431542) in the cardiac mesoderm differentiation stage downregulated PITX2c, a transcription factor that is known to inhibit the formation of the sinoatrial node in the left atrium during cardiac development. The resulting hiPSC-CMs were smaller in cell size, expressed higher pro-pacemaking transcription factors, TBX3 and TBX18, and exhibited pacemaking-like electrophysiological characteristics compared to control hiPSC-CMs differentiated from established Wnt-based protocol. The pacemaker-like subtype increased up to 2.4-fold in hiPSC-CMs differentiated with the addition of SB431542 relative to the control. Hence, Nodal inhibition in the cardiac mesoderm stage promoted pacemaker-like CM differentiation from hiPSCs. Improving the yield of human pacemaker-like CMs is a critical first step in the development of functional human cell-based biopacemakers.

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.

Chamber-specific transcriptional responses in atrial fibrillation

Atrial fibrillation (AF) is the most common cardiac arrhythmia, yet the molecular signature of the vulnerable atrial substrate is not well understood. Here, we delineated a distinct transcriptional signature in right versus left atrial cardiomyocytes (CMs) at baseline and identified chamber-specific gene expression changes in patients with a history of AF in the setting of end-stage heart failure (AF+HF) that are not present in heart failure alone (HF). We observed that human left atrial (LA) CMs exhibited Notch pathway activation and increased ploidy in AF+HF but not in HF alone. Transient activation of Notch signaling within adult CMs in a murine genetic model is sufficient to increase ploidy in both atrial chambers. Notch activation within LA CMs generated a transcriptomic fingerprint resembling AF, with dysregulation of transcription factor and ion channel genes, including Pitx2, Tbx5, Kcnh2, Kcnq1, and Kcnip2. Notch activation also produced distinct cellular electrophysiologic responses in LA versus right atrial CMs, prolonging the action potential duration (APD) without altering the upstroke velocity in the left atrium and reducing the maximal upstroke velocity without altering the APD in the right atrium. Our results support a shared human/murine model of increased Notch pathway activity predisposing to AF.

Distinct transcriptional changes occur in human left versus right atrial cardiomyocytes in atrial fibrillation, including Notch pathway activation, which alters electric properties and ploidy in murine models.

Conjugated activation of myocardial-specific transcription of Gja5 by a pair of Nkx2-5-Shox2 co-responsive elements

The sinoatrial node (SAN) is the primary pacemaker in the heart. During cardiogenesis, Shox2 and Nkx2-5 are co-expressed in the junction domain of the SAN and regulate pacemaker cell fate through a Shox2-Nkx2-5 antagonism. Cx40 is a marker of working myocardium and an Nkx2-5 transcriptional output antagonized by Shox2, but the underlying regulatory mechanisms remain elusive. Here we characterized a bona fide myocardial-specific Gja5 (coding gene of Cx40) distal enhancer consisting of a pair of Nkx2-5 and Shox2 co-bound elements in the regulatory region of Gja5. Transgenic reporter assays revealed that neither element alone, but the conjugation of both elements together, drives myocardial-specific transcription. Genetic analyses confirmed that the activation of this enhancer depends on Nkx2-5 but is inhibited by Shox2 in vivo, and its presence is essential for Gja5 expression in the myocardium but not the endothelial cells of the heart. Furthermore, chromatin conformation analysis showed an Nkx2-5-dependent loop formation between these two elements and the Gja5 promoter in vivo, indicating that Nkx2-5 bridges the conjugated activation of this enhancer by pairing the two elements to the Gja5 promoter.

Relations between plasma microRNAs, echocardiographic markers of atrial remodeling, and atrial fibrillation: Data from the Framingham Offspring study

BACKGROUND

Circulating microRNAs may reflect or influence pathological cardiac remodeling and contribute to atrial fibrillation (AF).

OBJECTIVE

The purpose of this study was to identify candidate plasma microRNAs that are associated with echocardiographic phenotypes of atrial remodeling, and incident and prevalent AF in a community-based cohort.

METHODS

We analyzed left atrial function index (LAFI) of 1788 Framingham Offspring 8 participants. We quantified expression of 339 plasma microRNAs. We examined associations between microRNA levels with LAFI and prevalent and incident AF. We constructed pathway analysis of microRNAs' predicted gene targets to identify molecular processes involved in adverse atrial remodeling in AF.

RESULTS

The mean age of the participants was 66 ± 9 years, and 54% were women. Five percent of participants had prevalent AF at the initial examination and 9% (n = 157) developed AF over a median 8.6 years of follow-up (IQR 8.1-9.2 years). Plasma microRNAs were associated with LAFI (N = 73, p<0.0001). Six of these plasma microRNAs were significantly associated with incident AF, including 4 also associated with prevalent AF (microRNAs 106b, 26a-5p, 484, 20a-5p). These microRNAs are predicted to regulate genes involved in cardiac hypertrophy, inflammation, and myocardial fibrosis.

CONCLUSIONS

Circulating microRNAs 106b, 26a-5p, 484, 20a-5p are associated with atrial remodeling and AF.

The Joint Analysis of Multi-Omics Data Revealed the Methylation-Expression Regulations in Atrial Fibrillation

Atrial fibrillation (AF) is one of the most prevalent heart rhythm disorder. The causes of AF include age, male sex, diabetes, hypertension, valve disease, and systolic/diastolic dysfunction. But on molecular level, its mechanisms are largely unknown. In this study, we collected 10 patients with persistent atrial fibrillation, 10 patients with paroxymal atrial fibrillation and 10 healthy individuals and did Methylation EPICBead Chip and RNA sequencing. By analyzing the methylation and gene expression data using machine learning based feature selection method Boruta, we identified the key genes that were strongly associated with AF and found their interconnections. The results suggested that the methylation of KIF15 may regulate the expression of PSMC3, TINAG, and NUDT6. The identified AF associated methylation-expression regulations may help understand the molecular mechanisms of AF from a multi-omics perspective.

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.

The local repolarization heterogeneity in the murine pulmonary veins myocardium contributes to the spatial distribution of the adrenergically induced ectopic foci

An atrial tachyarrhythmias is predominantly triggered by a proarrhythmic activity originate from the pulmonary veins (PV) myocardial sleeves; sympathetic or adrenergic stimulation facilitates PV proarrhythmia. In the present study the electrophysiological inhomogeneity, spatiotemporal characteristics of the adrenergically induced ectopic firing and sympathetic nerves distribution have been investigated in a murine PV myocardium to clarify mechanisms of adrenergic PV ectopy. Electrically paced murine PV demonstrate atrial-like pattern of conduction and atrial-like action potentials (AP) with longest duration in the mouth of PV. The application of norepinephrine (NE), agonists of α- and β-adrenergic receptors (ARs) or intracardiac nerves stimulation induced spontaneous AP in a form of periodical bursts or continuous firing. NE- or ARs agonists-induced SAP originated from unifocal ectopic foci with predominant localization in the region surrounding PV mouth, but not in the distal portions of a murine PV myocardium. A higher level of catecholamine content and catecholamine fiber network density was revealed in the PV myocardial sleeves relative to LA appendage. However, no significant local variation of catecholamine content and fiber density was observed in the murine PV. In conclusion, PV mouth region appear to be a most susceptible to adrenergic proarrhythmia in mice. Intrinsic spatial heterogeneity of AP duration can be considered as a factor influencing localization of the ectopic foci in PV.

Nkx2.5 insufficiency leads to atrial electrical remodeling through Wnt signaling in HL-1 cells

Homeobox protein Nxk-2.5 (Nkx2.5) is a homeobox transcription factor that promotes chamber-like myocardial gene expression. Data from a previous genome-wide association study suggested that Nkx2.5 may be associated with the genetic variation that underlies atrial fibrillation (AF). Nkx2.5 loss of function has been demonstrated to be associated with an increasing susceptibility of familial AF. Therefore, the aim of the present study was to investigate the effect of Nkx2.5 loss of function on electrophysiological substrates in HL-1 cells. To the best of our knowledge, the results demonstrated for the first time that Nkx2.5 expression was significantly decreased in a rat model exhibiting AF. The effect of silencing Nkx2.5 was assessed following transfection with adenoviral vectors with specific NKX2.5-shRNA. The effect of Nkx2.5 silencing on potassium/sodium hyperpolarization-activated cyclic nucleotide-gated channel 4 (HCN4), gap junction alpha-5 protein (Cx40), calcium handling proteins and protein Wnt-11 (Wnt11) expression levels was also assessed in HL-1 cells. The results revealed that Nkx2.5 silencing increased HCN4 expression, decreased Cx40 expression and disrupted the expression of calcium handling proteins. Additionally, Wnt11 signal protein expression was decreased following Nkx2.5 silencing. The results of the present study demonstrated that Nkx2.5 served as a transcriptional regulator of the electrophysiological substrates associated with AF.

Shox2 regulates osteogenic differentiation and pattern formation during hard palate development in mice

During mammalian palatogenesis, cranial neural crest-derived mesenchymal cells undergo osteogenic differentiation and form the hard palate, which is divided into palatine process of the maxilla and the palatine. However, it remains unknown whether these bony structures originate from the same cell lineage and how the hard palate is patterned at the molecular level. Using mice, here we report that deficiency in Shox2 (short stature homeobox 2), a transcriptional regulator whose expression is restricted to the anterior palatal mesenchyme, leads to a defective palatine process of the maxilla but does not affect the palatine. Shox2 overexpression in palatal mesenchyme resulted in a hyperplastic palatine process of the maxilla and a hypoplastic palatine. RNA sequencing and assay for transposase-accessible chromatin-sequencing analyses revealed that Shox2 controls the expression of pattern specification and skeletogenic genes associated with accessible chromatin in the anterior palate. This highlighted a lineage-autonomous function of Shox2 in patterning and osteogenesis of the hard palate. H3K27ac ChIP-Seq and transient transgenic enhancer assays revealed that Shox2 binds distal-acting cis-regulatory elements in an anterior palate-specific manner. Our results suggest that the palatine process of the maxilla and palatine arise from different cell lineages and differ in ossification mechanisms. Shox2 evidently controls osteogenesis of a cell lineage and contributes to the palatine process of the maxilla by interacting with distal cis-regulatory elements to regulate skeletogenic gene expression and to pattern the hard palate. Genome-wide Shox2 occupancy in the developing palate may provide a marker for identifying active anterior palate-specific gene enhancers.

Transcription Factor TBX18 Reprograms Vascular Smooth Muscle Cells of Ascending Aorta to Pacemaker-Like Cells

Vascular smooth muscle cells (VSMCs) of ascending aorta and TBX18⁺ sinus node both originated from the second heart field. The study explored whether ascending aortic smooth muscle cells in vitro could be reprogrammed into pacemaker-like cells with human TBX18. In the study, VSMCs were infected with TBX18, and then cocultured with neonatal rat ventricular cardiomyocytes (NRVMs) in vitro. By overexpressing TBX18, the transfected VSMCs expressed high levels of hyperpolarization-activated cyclic nucleotide-gated channel 4 (HCN4), insulin gene enhancer binding protein 1, and human dwarf homeobox gene SHOX2, cardiac troponin I, and low level of connexin 43. In addition, funny current (I_f) was recorded by patch clamp appeared the time and voltage dependence in TBX18 group, which the amplitude of I_f density was from -5.164 ± 0.662 pA/pF to -0.765 ± 0.358 pA/pF (n = 14). Furthermore, TBX18-transfected VSMCs coupled with NRVMs showed typical action potential of pacemaker-like cells and the beating rate was faster (178.00 ± 7.55 bpm, p < 0.05) compared with other groups. In conclusion, our study indicated that transcription factor TBX18 could reprogram VSMCs into pacemaker-like cells in vitro.

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.

Learn from Your Elders: Developmental Biology Lessons to Guide Maturation of Stem Cell-Derived Cardiomyocytes

Human pluripotent stem cells (hPSCs) offer a multifaceted platform to study cardiac developmental biology, understand disease mechanisms, and develop novel therapies. Remarkable progress over the last two decades has led to methods to obtain highly pure hPSC-derived cardiomyocytes (hPSC-CMs) with reasonable ease and scalability. Nevertheless, a major bottleneck for the translational application of hPSC-CMs is their immature phenotype, resembling that of early fetal cardiomyocytes. Overall, bona fide maturation of hPSC-CMs represents one of the most significant goals facing the field today. Developmental biology studies have been pivotal in understanding the mechanisms to differentiate hPSC-CMs. Similarly, evaluation of developmental cues such as electrical and mechanical activities or neurohormonal and metabolic stimulations revealed the importance of these pathways in cardiomyocyte physiological maturation. Those signals cooperate and dictate the size and the performance of the developing heart. Likewise, this orchestra of stimuli is important in promoting hPSC-CM maturation, as demonstrated by current in vitro maturation approaches. Different shades of adult-like phenotype are achieved by prolonging the time in culture, electromechanical stimulation, patterned substrates, microRNA manipulation, neurohormonal or metabolic stimulation, and generation of human-engineered heart tissue (hEHT). However, mirroring this extremely dynamic environment is challenging, and reproducibility and scalability of these approaches represent the major obstacles for an efficient production of mature hPSC-CMs. For this reason, understanding the pattern behind the mechanisms elicited during the late gestational and early postnatal stages not only will provide new insights into postnatal development but also potentially offer new scalable and efficient approaches to mature hPSC-CMs.

Smad7 Regulates Dental Epithelial Proliferation during Tooth Development

Tooth morphogenesis involves dynamic changes in shape and size as it proceeds through the bud, cap, and bell stages. This process requires exact regulation of cell proliferation and differentiation. Smad7, a general antagonist against transforming growth factor-β (TGF-β) signaling, is necessary for maintaining homeostasis and proper functionality in many organs. While TGF-β signaling is widely involved in tooth morphogenesis, the precise role of Smad7 in tooth development remains unknown. In this study, we showed that Smad7 is expressed in the developing mouse molars with a high level in the dental epithelium but a moderate to weak level in the dental mesenchyme. Smad7 deficiency led to a profound decrease in tooth size primarily due to a severely compromised cell proliferation capability in the dental epithelium. Consistent with the tooth shrinkage phenotype, RNA sequencing (RNA-seq) analysis revealed that Smad7 ablation downregulated genes referred to epithelial cell proliferation and cell cycle G1/S phase transition, whereas the upregulated genes were involved in responding to TGF-β signaling and cell cycle arrest. Among these genes, the expression of Cdkn1a (encoding p21), a negative cell proliferation regulator, was remarkably elevated in parallel with the diminution of Ccnd1 encoding the crucial cell cycle regulator cyclin D1 in the dental epithelium. Meanwhile, the expression level of p-Smad2/3 was ectopically elevated in the developing tooth germ of Smad7 null mice, indicating the hyperactivation of the canonical TGF-β signaling. These effects were reversed by addition of TGF-β signaling inhibitor in cell cultures of Smad7^-/- molar tooth germs, with rescued expression of cyclin D1 and cell proliferation rate. In sum, our studies demonstrate that Smad7 functions primarily as a positive regulator of cell proliferation via inhibition of the canonical TGF-β signaling during dental epithelium development and highlight a crucial role for Smad7 in regulating tooth size.

Transcriptional regulation of the cardiac conduction system

The rate and rhythm of heart muscle contractions are coordinated by the cardiac conduction system (CCS), a generic term for a collection of different specialized muscular tissues within the heart. The CCS components initiate the electrical impulse at the sinoatrial node, propagate it from atria to ventricles via the atrioventricular node and bundle branches, and distribute it to the ventricular muscle mass via the Purkinje fibre network. The CCS thereby controls the rate and rhythm of alternating contractions of the atria and ventricles. CCS function is well conserved across vertebrates from fish to mammals, although particular specialized aspects of CCS function are found only in endotherms (mammals and birds). The development and homeostasis of the CCS involves transcriptional and regulatory networks that act in an embryonic-stage-dependent, tissue-dependent, and dose-dependent manner. This Review describes emerging data from animal studies, stem cell models, and genome-wide association studies that have provided novel insights into the transcriptional networks underlying CCS formation and function. How these insights can be applied to develop disease models and therapies is also discussed.

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.

The lncRNA Hand2os1/Uph locus orchestrates heart development through regulation of precise expression of Hand2

Exploration and dissection of potential actions and effects of long noncoding RNA (lncRNA) in animals remain challenging. Here, using multiple knockout mouse models and single cell RNA sequencing, we demonstrate that the divergent lncRNA Hand2os1/Uph has a key complex modulatory effect on the expression of its neighboring gene HAND2 and subsequently on heart development and function. Short deletion of the Hand2os1 promoter in mouse diminishes Hand2os1 transcription to ∼8-32%, but fails to affect HAND2 expression and yields no discernable heart phenotypes. Interestingly, full-length deletion of Hand2os1 in mouse causes moderate yet prevalent upregulation of HAND2 in hundreds of cardiac cells, leading to profound biological consequences, including dysregulated cardiac gene programs, congenital heart defects and perinatal lethality. We propose that the Hand2os1 locus dampens HAND2 expression to restrain cardiomyocyte proliferation, thereby orchestrating a balanced development of cardiac cell lineages. This study highlights the regulatory complexity of the lncRNA Hand2os1 on HAND2 expression, emphasizing the need for complementary genetic and single cell approaches to delineate the function and primary molecular effects of an lncRNA in animals.

Geometrical Patterning and Constituent Cell Heterogeneity Facilitate Electrical Conduction Disturbances in a Human Induced Pluripotent Stem Cell-Based Platform: An In vitro Disease Model of Atrial Arrhythmias

Ectopic foci from pulmonary veins (PVs) comprise the main trigger associated with the initiation of atrial fibrillation (AF). An abrupt anatomical narrow-to-wide transition, modeled as in vitro geometrical patterning with similar configuration in the present study, is located at the junction of PVs and the left atrium (LA). Complex cellular composition, i.e., constituent cell heterogeneity, is also observed in PVs and the PVs-LA junction. High frequency triggers accompanied with anatomical irregularity and constituent cell heterogeneity provoke impaired conduction, a prerequisite for AF genesis. However, few experiments investigating the effects of these factors on electrophysiological properties using human-based cardiomyocytes (CMs) with atrial properties have been reported. The aim of the current study was to estimate whether geometrical patterning and constituent cell heterogeneity under high frequency stimuli undergo conduction disturbance utilizing an in vitro two-dimensional (2D) monolayer preparation consisting of atrial-like CMs derived from human induced pluripotent stem cells (hiPSCs) and atrial fibroblasts (Fbs). We induced hiPSCs into atrial-like CMs using a directed cardiac differentiation protocol with the addition of all-trans retinoic acid (ATRA). The atrial-like hiPSC-derived CMs (hiPSC-CMs) and atrial Fbs were transferred in defined ratios (CMs/Fbs: 100%/0% or 70%/30%) on manually fabricated plates with or without geometrical patterning imitating the PVs-LA junction. High frequency field stimulation emulating repetitive ectopic foci originated in PVs were delivered, and the electrical propagation was assessed by optical mapping. We generated high purity CMs with or without the ATRA application. ATRA-treated hiPSC-CMs exhibited significantly higher atrial-specific properties by immunofluorescence staining, gene expression patterns, and optical action potential parameters than those of ATRA-untreated hiPSC-CMs. Electrical stimuli at a higher frequency preferentially induced impaired electrical conduction on atrial-like hiPSC-CMs monolayer preparations with an abrupt geometrical transition than on those with uniform geometry. Additionally, the application of human atrial Fbs to the geometrically patterned atrial-like hiPSC-CMs tended to further deteriorate the integrity of electrical conduction compared with those using the atrial-like hiPSC-CM alone preparations. Thus, geometrical narrow-to-wide patterning under high frequency stimuli preferentially jeopardized electrical conduction within in vitro atrial-like hiPSC-CM monolayers. Constituent cell heterogeneity represented by atrial Fbs also contributed to the further deterioration of conduction stability.

Transcriptome analysis of mouse and human sinoatrial node cells reveals a conserved genetic program

The rate of contraction of the heart relies on proper development and function of the sinoatrial node, which consists of a small heterogeneous cell population, including Tbx3⁺ pacemaker cells. Here, we have isolated and characterized the Tbx3⁺ cells from Tbx3 ^+/Venus knock-in mice. We studied electrophysiological parameters during development and found that Venus-labeled cells are genuine Tbx3⁺ pacemaker cells. We analyzed the transcriptomes of late fetal FACS-purified Tbx3⁺ sinoatrial nodal cells and Nppb-Katushka⁺ atrial and ventricular chamber cardiomyocytes, and identified a sinoatrial node-enriched gene program, including key nodal transcription factors, BMP signaling and Smoc2, the disruption of which in mice did not affect heart rhythm. We also obtained the transcriptomes of the sinoatrial node region, including pacemaker and other cell types, and right atrium of human fetuses, and found a gene program including TBX3, SHOX2, ISL1 and HOX family members, and BMP and NOTCH signaling components conserved between human and mouse. We conclude that a conserved gene program characterizes the sinoatrial node region and that the Tbx3 ^+/Venus allele provides a reliable tool for visualizing the sinoatrial node, and studying its development and function.

Bone morphogenetic protein 4 promotes the differentiation of Tbx18-positive epicardial progenitor cells to pacemaker-like cells

Clarifying the mechanisms via which pacemaker- like cells are generated is critical for identifying novel targets for arrhythmia-associated disorders and constructing pacemakers with the ability to adapt to physiological requirements. T-box 18 (Tbx18)⁺ epicardial progenitor cells (EPCs) have the potential to differentiate into pacemaker cells. Although bone morphogenetic protein 4 (Bmp4) is likely to contribute, its role and regulatory mechanisms in the differentiation of Tbx18⁺ EPCs into pacemaker-like cells have remained to be fully elucidated. In the present study, the association between Bmp4, GATA binding protein 4 (Gata4) and hyperpolarization- activated cyclic nucleotide gated potassium channel 4 (Hcn4) to regulate NK2 homeobox 5 (Nkx2.5), which is known to be required for the differentiation of Tbx18⁺ EPCs into pacemaker-like cells, was assessed. Tbx18⁺ EPCs were isolated from Tbx18:Cre/Rosa26R^{enhanced yellow fluorescence protein (EYFP)} murine embryos at embryonic day 11.5 and divided into the following four treatment groups: Control, Bmp4, Bmp4+LDN193189 (a Bmp inhibitor) and LDN193189. In vitro Bmp4 promoted the expression of Hcn4 in Tbx18⁺ EPCs via lineage tracing of Tbx18:Cre/Rosa26R^EYFP mice, which was likely due to upregulation of Gata4 expression. Gata4 knockdown experiments were then performed using the following five treatment groups: Control, control small interfering RNA (siRNA), Bmp4, Bmp4+siRNA targeting Gata4 (siGata4) and siGata4 group. Knockdown of Gata4 caused a downregulation of Hcn4 and an upregulation of Nkx2.5, but had no effect on Bmp4 expression. In conclusion, it was indicated that in Tbx18⁺ EPCs, the expression of Nkx2.5 was regulated by Bmp4 via Gata4. Taken together, these results provide important information on regulatory networks of pacemaker cell differentiation and may serve as a basis for further studies.

RHOA-ROCK signalling is necessary for lateralization and differentiation of the developing sinoatrial node

Aims

RHOA-ROCK signalling regulates cell migration, proliferation, differentiation, and transcription. RHOA is expressed in the developing cardiac conduction system in chicken and mice. In early development, the entire sinus venosus myocardium, including both the transient left-sided and the definitive sinoatrial node (SAN), has pacemaker potential. Later, pacemaker potential is restricted to the right-sided SAN. Disruption of RHOA expression in adult mice causes arrhythmias including bradycardia and atrial fibrillation, the mechanism of which is unknown but presumed to affect the SAN. The aim of this study is to assess the role of RHOA-ROCK signalling in SAN development in the chicken heart.

Methods and results

ROCK signalling was inhibited chemically in embryonic chicken hearts using Y-27632. This prolonged the immature state of the sinus venosus myocardium, evidenced by up-regulation of the transcription factor ISL1, wide distribution of pacemaker potential, and significantly reduced heart rate. Furthermore ROCK inhibition caused aberrant expression of typical SAN genes: ROCK1, ROCK2, SHOX2, TBX3, TBX5, ISL1, HCN4, CX40, CAV3.1, and NKX2.5 and left-right asymmetry genes: PITX2C and NODAL. Anatomical abnormalities in pulmonary vein development were also observed. Patch clamp electrophysiology confirmed the immature phenotype of the SAN cells and a residual left-sided sinus venosus myocardium pacemaker-like potential.

Conclusions

RHOA-ROCK signalling is involved in establishing the right-sided SAN as the definitive pacemaker of the heart and restricts typical pacemaker gene expression to the right side of the sinus venosus myocardium.

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.

Comparative developmental biology of the cardiac inflow tract

The vertebrate heart receives the blood through the cardiac inflow tract. This area has experienced profound changes along the evolution of vertebrates; changes that have a reflection in the cardiac ontogeny. The development of the inflow tract involves dynamic changes due to the progressive addition of tissue derived from the secondary heart field. The inflow tract is the site where oxygenated blood coming from lungs is received separately from the systemic return, where the cardiac pacemaker is established and where the proepicardium develops. Differential cell migration towards the inflow tract breaks the symmetry of the primary heart tube and determines the direction of the cardiac looping. In air-breathing vertebrates, an inflow tract reorganization is essential to keep separate blood flows from systemic and pulmonary returns. Finally, the sinus venosus endocardium has recently been recognized as playing a role in the constitution of the coronary vasculature. Due to this developmental complexity, congenital anomalies of the inflow tract can cause severe cardiac diseases. We aimed to review the recent literature on the cellular and molecular mechanisms that regulate the morphogenesis of the cardiac inflow tract, together with comparative and evolutionary details, thus providing a basis for a better understanding of these mechanisms.

Spatially resolved RNA-sequencing of the embryonic heart identifies a role for Wnt/β-catenin signaling in autonomic control of heart rate

Development of specialized cells and structures in the heart is regulated by spatially -restricted molecular pathways. Disruptions in these pathways can cause severe congenital cardiac malformations or functional defects. To better understand these pathways and how they regulate cardiac development we used tomo-seq, combining high-throughput RNA-sequencing with tissue-sectioning, to establish a genome-wide expression dataset with high spatial resolution for the developing zebrafish heart. Analysis of the dataset revealed over 1100 genes differentially expressed in sub-compartments. Pacemaker cells in the sinoatrial region induce heart contractions, but little is known about the mechanisms underlying their development. Using our transcriptome map, we identified spatially restricted Wnt/β-catenin signaling activity in pacemaker cells, which was controlled by Islet-1 activity. Moreover, Wnt/β-catenin signaling controls heart rate by regulating pacemaker cellular response to parasympathetic stimuli. Thus, this high-resolution transcriptome map incorporating all cell types in the embryonic heart can expose spatially restricted molecular pathways critical for specific cardiac functions.