Islet1 is a direct transcriptional target of the homeodomain transcription factor Shox2 and rescues the Shox2-mediated bradycardia

The cardiac conduction system: History, development, and disease

The heart is the first organ to form during embryonic development, establishing the circulatory infrastructure necessary to sustain life and enable downstream organogenesis. Critical to the heart's function is its ability to initiate and propagate electrical impulses that allow for the coordinated contraction and relaxation of its chambers, and thus, the movement of blood and nutrients. Several specialized structures within the heart, collectively known as the cardiac conduction system (CCS), are responsible for this phenomenon. In this review, we discuss the discovery and scientific history of the mammalian cardiac conduction system as well as the key genes and transcription factors implicated in the formation of its major structures. We also describe known human diseases related to CCS development and explore existing challenges in the clinical context.

Cadherin-5 facilitated the differentiation of human induced pluripotent stem cells into sinoatrial node-like pacemaker cells by regulating β-catenin

Our study was conducted to investigate whether cadherin-5 (CDH5), a vascular endothelial cell adhesion glycoprotein, could facilitate the differentiation of human induced pluripotent stem cells (hiPSCs) into sinoatrial node-like pacemaker cells (SANLPCs), following previous findings of silk-fibroin hydrogel-induced direct conversion of quiescent cardiomyocytes into pacemaker cells in rats through the activation of CDH5. In this study, the differentiating hiPSCs were treated with CDH5 (40 ng/mL) between Day 5 and 7 during cardiomyocytes differentiation. The findings in the present study demonstrated that CDH5 stimulated the expression of pacemaker-specific markers while suppressing markers associated with working cardiomyocytes, resulting in an increased proportion of SANLPCs among hiPSCs-derived cardiomyocytes (hiPSC-CMs) population. Moreover, CDH5 induced typical electrophysiological characteristics resembling cardiac pacemaker cells in hiPSC-CMs. Further mechanistic investigations revealed that the enriched differentiation of hiPSCs into SANLPCs induced by CDH5 was partially reversed by iCRT14, an inhibitor of β-catenin. Therefore, based on the aforementioned findings, it could be inferred that the regulation of β-catenin by CDH5 played a crucial role in promoting the enriched differentiation of hiPSCs into SANLPCs, which presents a novel avenue for the construction of biological pacemakers in forthcoming research.

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.

Lack of the transcription factor Nfix causes tachycardia in mice sinus node and rats neonatal cardiomyocytes

AIMS

Nfix is a transcription factor belonging to the Nuclear Factor I (NFI) family comprising four members (Nfia, b, c, x). Nfix plays important roles in the development and function of several organs. In muscle development, Nfix controls the switch from embryonic to fetal myogenesis by promoting fast twitching fibres. In the adult muscle, following injury, lack of Nfix impairs regeneration, inducing higher content of slow-twitching fibres. Nfix is expressed also in the heart, but its function has been never investigated before. We studied Nfix role in this organ.

METHODS

Using Nfix-null and wild type (WT) mice we analyzed: (1) the expression pattern of Nfix during development by qPCR and (2) the functional alterations caused by its absence, by in vivo telemetry and in vitro patch clamp analysis.

RESULTS AND CONCLUSIONS

Nfix expression start in the heart from E12.5. Adult hearts of Nfix-null mice show a hearts morphology and sarcomeric proteins expression similar to WT. However, Nfix-null animals show tachycardia that derives form an intrinsic higher beating rate of the sinus node (SAN). Molecular and functional analysis revealed that sinoatrial cells of Nfix-null mice express a significantly larger L-type calcium current (Cacna1d + Cacna1c). Interestingly, downregulation of Nfix by sh-RNA in primary cultures of neonatal rat ventricular cardiomyocytes induced a similar increase in their spontaneous beating rate and in ICaL current. In conclusion, our data provide the first demonstration of a role of Nfix that, increasing the L-type calcium current, modulates heart rate.

Nr2f1a maintains atrial nkx2.5 expression to repress pacemaker identity within venous atrial cardiomyocytes of zebrafish

Maintenance of cardiomyocyte identity is vital for normal heart development and function. However, our understanding of cardiomyocyte plasticity remains incomplete. Here, we show that sustained expression of the zebrafish transcription factor Nr2f1a prevents the progressive acquisition of ventricular cardiomyocyte (VC) and pacemaker cardiomyocyte (PC) identities within distinct regions of the atrium. Transcriptomic analysis of flow-sorted atrial cardiomyocytes (ACs) from nr2f1a mutant zebrafish embryos showed increased VC marker gene expression and altered expression of core PC regulatory genes, including decreased expression of nkx2.5, a critical repressor of PC differentiation. At the arterial (outflow) pole of the atrium in nr2f1a mutants, cardiomyocytes resolve to VC identity within the expanded atrioventricular canal. However, at the venous (inflow) pole of the atrium, there is a progressive wave of AC transdifferentiation into PCs across the atrium toward the arterial pole. Restoring Nkx2.5 is sufficient to repress PC marker identity in nr2f1a mutant atria and analysis of chromatin accessibility identified an Nr2f1a-dependent nkx2.5 enhancer expressed in the atrial myocardium directly adjacent to PCs. CRISPR/Cas9-mediated deletion of the putative nkx2.5 enhancer leads to a loss of Nkx2.5-expressing ACs and expansion of a PC reporter, supporting that Nr2f1a limits PC differentiation within venous ACs via maintaining nkx2.5 expression. The Nr2f-dependent maintenance of AC identity within discrete atrial compartments may provide insights into the molecular etiology of concurrent structural congenital heart defects and associated arrhythmias.

Novel Candidate Regulators and Developmental Trajectory of Pituitary Thyrotropes

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

Comparative Study of Transcriptome in the Hearts Isolated from Mice, Rats, and Humans

The heart is a significant organ in mammalian life, and the heartbeat mechanism has been an essential focus of science. However, few studies have focused on species differences. Accordingly, challenges remain in studying genes that have universal functions across species and genes that determine species differences. Here, we analyzed transcriptome data in mouse, rat, and human atria, ventricles, and sinoatrial nodes (SA) obtained from different platforms and compared them by calculating specificity measure (SPM) values in consideration of species differences. Among the three heart regions, the species differences in SA were the greatest, and we searched for genes that determined the essential characteristics of SA, which was SHOX2 in our criteria. The SPM value of SHOX2 was prominently high across species. Similarly, by calculating SPM values, we identified 3 atrial-specific, 11 ventricular-specific, and 17 SA-specific markers. Ontology analysis identified 70 cardiac region- and species-specific ontologies. These results suggest that reanalyzing existing data by calculating SPM values may identify novel tissue-specific genes and species-dependent gene expression. This study identified the importance of SHOX2 as an SA-specific transcription factor, a novel cardiac regional marker, and species-dependent ontologies.

Preferential Expression of Ca2+-Stimulable Adenylyl Cyclase III in the Supraventricular Area, Including Arrhythmogenic Pulmonary Vein of the Rat Heart

Ectopic excitability in pulmonary veins (PVs) is the major cause of atrial fibrillation. We previously reported that the inositol trisphosphate receptor in rat PV cardiomyocytes cooperates with the Na+-Ca2+ exchanger to provoke ectopic automaticity in response to norepinephrine. Here, we focused on adenylyl cyclase (AC) as another effector of norepinephrine stimulation. RT-PCR, immunohistochemistry, and Western blotting revealed that the abundant expression of Ca2+-stimulable AC3 was restricted to the supraventricular area, including the PVs. All the other AC isotypes hardly displayed any region-specific expressions. Immunostaining of isolated cardiomyocytes showed an enriched expression of AC3 along the t-tubules in PV myocytes. The cAMP-dependent response of L-type Ca2+ currents in the PV and LA cells is strengthened by the 0.1 mM intracellular Ca2+ condition, unlike in the ventricular cells. The norepinephrine-induced automaticity of PV cardiomyocytes was reversibly suppressed by 100 µM SQ22536, an adenine-like AC inhibitor. These findings suggest that the specific expression of AC3 along t-tubules may contribute to arrhythmogenic automaticity in rat PV cardiomyocytes.

SARS-CoV-2 Infection Induces Ferroptosis of Sinoatrial Node Pacemaker Cells

BACKGROUND

Increasing evidence suggests that cardiac arrhythmias are frequent clinical features of coronavirus disease 2019 (COVID-19). Sinus node damage may lead to bradycardia. However, it is challenging to explore human sinoatrial node (SAN) pathophysiology due to difficulty in isolating and culturing human SAN cells. Embryonic stem cells (ESCs) can be a source to derive human SAN-like pacemaker cells for disease modeling.

METHODS

We used both a hamster model and human ESC (hESC)-derived SAN-like pacemaker cells to explore the impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection on the pacemaker cells of the heart. In the hamster model, quantitative real-time polymerase chain reaction and immunostaining were used to detect viral RNA and protein, respectively. We then created a dual knock-in SHOX2:GFP;MYH6:mCherry hESC reporter line to establish a highly efficient strategy to derive functional human SAN-like pacemaker cells, which was further characterized by single-cell RNA sequencing. Following exposure to SARS-CoV-2, quantitative real-time polymerase chain reaction, immunostaining, and RNA sequencing were used to confirm infection and determine the host response of hESC-SAN-like pacemaker cells. Finally, a high content chemical screen was performed to identify drugs that can inhibit SARS-CoV-2 infection, and block SARS-CoV-2-induced ferroptosis.

RESULTS

Viral RNA and spike protein were detected in SAN cells in the hearts of infected hamsters. We established an efficient strategy to derive from hESCs functional human SAN-like pacemaker cells, which express pacemaker markers and display SAN-like action potentials. Furthermore, SARS-CoV-2 infection causes dysfunction of human SAN-like pacemaker cells and induces ferroptosis. Two drug candidates, deferoxamine and imatinib, were identified from the high content screen, able to block SARS-CoV-2 infection and infection-associated ferroptosis.

CONCLUSIONS

Using a hamster model, we showed that primary pacemaker cells in the heart can be infected by SARS-CoV-2. Infection of hESC-derived functional SAN-like pacemaker cells demonstrates ferroptosis as a potential mechanism for causing cardiac arrhythmias in patients with COVID-19. Finally, we identified candidate drugs that can protect the SAN cells from SARS-CoV-2 infection.

Retinoic acid signaling in heart development: Application in the differentiation of cardiovascular lineages from human pluripotent stem cells

Retinoic acid (RA) signaling plays an important role during heart development in establishing anteroposterior polarity, formation of inflow and outflow tract progenitors, and growth of the ventricular compact wall. RA is also utilized as a key ingredient in protocols designed for generating cardiac cell types from pluripotent stem cells (PSCs). This review discusses the role of RA in cardiogenesis, currently available protocols that employ RA for differentiation of various cardiovascular lineages, and plausible transcriptional mechanisms underlying this fate specification. These insights will inform further development of desired cardiac cell types from human PSCs and their application in preclinical and clinical research.

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

RESEARCH PURPOSE

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

PRINCIPAL RESULTS

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

MAJOR CONCLUSIONS

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

Identification and Tissue-Specific Characterization of Novel SHOX-Regulated Genes in Zebrafish Highlights SOX Family Members Among Other Genes

SHOX deficiency causes a spectrum of clinical phenotypes related to skeletal dysplasia and short stature, including Léri-Weill dyschondrosteosis, Langer mesomelic dysplasia, Turner syndrome, and idiopathic short stature. SHOX controls chondrocyte proliferation and differentiation, bone maturation, and cellular growth arrest and apoptosis via transcriptional regulation of its direct target genes NPPB, FGFR3, and CTGF. However, our understanding of SHOX-related pathways is still incomplete. To elucidate the underlying molecular mechanisms and to better understand the broad phenotypic spectrum of SHOX deficiency, we aimed to identify novel SHOX targets. We analyzed differentially expressed genes in SHOX-overexpressing human fibroblasts (NHDF), and confirmed the known SHOX target genes NPPB and FGFR among the most strongly regulated genes, together with 143 novel candidates. Altogether, 23 genes were selected for further validation, first by whole-body characterization in developing shox-deficient zebrafish embryos, followed by tissue-specific expression analysis in three shox-expressing zebrafish tissues: head (including brain, pharyngeal arches, eye, and olfactory epithelium), heart, and pectoral fins. Most genes were physiologically relevant in the pectoral fins, while only few genes were also significantly regulated in head and heart tissue. Interestingly, multiple sox family members (sox5, sox6, sox8, and sox18) were significantly dysregulated in shox-deficient pectoral fins together with other genes (nppa, nppc, cdkn1a, cdkn1ca, cyp26b1, and cy26c1), highlighting an important role for these genes in shox-related growth disorders. Network-based analysis integrating data from the Ingenuity pathways revealed that most of these genes act in a common network. Our results provide novel insights into the genetic pathways and molecular events leading to the clinical manifestation of SHOX deficiency.

Genetic Complexity of Sinoatrial Node Dysfunction

The pacemaker cells of the cardiac sinoatrial node (SAN) are essential for normal cardiac automaticity. Dysfunction in cardiac pacemaking results in human sinoatrial node dysfunction (SND). SND more generally occurs in the elderly population and is associated with impaired pacemaker function causing abnormal heart rhythm. Individuals with SND have a variety of symptoms including sinus bradycardia, sinus arrest, SAN block, bradycardia/tachycardia syndrome, and syncope. Importantly, individuals with SND report chronotropic incompetence in response to stress and/or exercise. SND may be genetic or secondary to systemic or cardiovascular conditions. Current management of patients with SND is limited to the relief of arrhythmia symptoms and pacemaker implantation if indicated. Lack of effective therapeutic measures that target the underlying causes of SND renders management of these patients challenging due to its progressive nature and has highlighted a critical need to improve our understanding of its underlying mechanistic basis of SND. This review focuses on current information on the genetics underlying SND, followed by future implications of this knowledge in the management of individuals with SND.

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.

Atrial and Sinoatrial Node Development in the Zebrafish Heart

Proper development and function of the vertebrate heart is vital for embryonic and postnatal life. Many congenital heart defects in humans are associated with disruption of genes that direct the formation or maintenance of atrial and pacemaker cardiomyocytes at the venous pole of the heart. Zebrafish are an outstanding model for studying vertebrate cardiogenesis, due to the conservation of molecular mechanisms underlying early heart development, external development, and ease of genetic manipulation. Here, we discuss early developmental mechanisms that instruct appropriate formation of the venous pole in zebrafish embryos. We primarily focus on signals that determine atrial chamber size and the specialized pacemaker cells of the sinoatrial node through directing proper specification and differentiation, as well as contemporary insights into the plasticity and maintenance of cardiomyocyte identity in embryonic zebrafish hearts. Finally, we integrate how these insights into zebrafish cardiogenesis can serve as models for human atrial defects and arrhythmias.

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.

ISL1 loss-of-function variation causes familial atrial fibrillation

Atrial fibrillation (AF) represents the most frequent form of sustained cardiac rhythm disturbance, affecting approximately 1% of the general population worldwide, and confers a substantially enhanced risk of cerebral stroke, heart failure, and death. Increasing epidemiological studies have clearly demonstrated a strong genetic basis for AF, and variants in a wide range of genes, including those coding for ion channels, gap junction channels, cardiac structural proteins and transcription factors, have been identified to underlie AF. Nevertheless, the genetic pathogenesis of AF is complex and still far from completely understood. Here, whole-exome sequencing and bioinformatics analyses of a three-generation family with AF were performed, and after filtering variants by multiple metrics, we identified a heterozygous variant in the ISL1 gene (encoding a transcription factor critical for embryonic cardiogenesis and postnatal cardiac remodeling), NM_002202.2: c.481G > T; p.(Glu161*), which was validated by Sanger sequencing and segregated with autosome-dominant AF in the family with complete penetrance. The nonsense variant was absent from 284 unrelated healthy individuals used as controls. Functional assays with a dual-luciferase reporter assay system revealed that the truncating ISL1 protein lost transcriptional activation on the verified target genes MEF2C and NKX2-5. Additionally, the variant nullified the synergistic transactivation between ISL1 and TBX5 as well as GATA4, two other transcription factors that have been implicated in AF. The findings suggest ISL1 as a novel gene contributing to AF, which adds new insight to the genetic mechanisms underpinning AF, implying potential implications for genetic testing and risk stratification of the AF family members.

Reptiles as a Model System to Study Heart Development

A chambered heart is common to all vertebrates, but reptiles show unparalleled variation in ventricular septation, ranging from almost absent in tuataras to full in crocodilians. Because mammals and birds evolved independently from reptile lineages, studies on reptile development may yield insight into the evolution and development of the full ventricular septum. Compared with reptiles, mammals and birds have evolved several other adaptations, including compact chamber walls and a specialized conduction system. These adaptations appear to have evolved from precursor structures that can be studied in present-day reptiles. The increase in the number of studies on reptile heart development has been greatly facilitated by sequencing of several genomes and the availability of good staging systems. Here, we place reptiles in their phylogenetic context with a focus on features that are primitive when compared with the homologous features of mammals. Further, an outline of major developmental events is given, and variation between reptile species is discussed.

Muscarinic receptors promote pacemaker fate at the expense of secondary conduction system tissue in zebrafish

Deterioration or inborn malformations of the cardiac conduction system (CCS) interfere with proper impulse propagation in the heart and may lead to sudden cardiac death or heart failure. Patients afflicted with arrhythmia depend on antiarrhythmic medication or invasive therapy, such as pacemaker implantation. An ideal way to treat these patients would be CCS tissue restoration. This, however, requires precise knowledge regarding the molecular mechanisms underlying CCS development. Here, we aimed to identify regulators of CCS development. We performed a compound screen in zebrafish embryos and identified tolterodine, a muscarinic receptor antagonist, as a modifier of CCS development. Tolterodine provoked a lower heart rate, pericardiac edema, and arrhythmia. Blockade of muscarinic M3, but not M2, receptors induced transcriptional changes leading to amplification of sinoatrial cells and loss of atrioventricular identity. Transcriptome data from an engineered human heart muscle model provided additional evidence for the contribution of muscarinic M3 receptors during cardiac progenitor specification and differentiation. Taken together, we found that muscarinic M3 receptors control the CCS already before the heart becomes innervated. Our data indicate that muscarinic receptors maintain a delicate balance between the developing sinoatrial node and the atrioventricular canal, which is probably required to prevent the development of arrhythmia.

Transcription Factor prrx1 Promotes Brown Adipose-Derived Stem Cells Differentiation to Sinus Node-Like Cells

This study investigated whether overexpression of paired-related homeobox 1 (prrx1) can successfully induce differentiation of brown adipose-derived stem cells (BADSCs) into sinus node-like cells. The experiments were performed in two groups: adenovirus-green fluorescent protein (Ad-GFP) group and Ad-prrx1 group. After 5-7 days of adenoviral transfection, the expression levels of sinus node cell-associated pacing protein (hyperpolarization-activated cyclic nucleotide-gated potassium channel 4 [HCN4]) and ion channel (calcium channel, voltage-dependent, T type, alpha 1G subunit [Cacna1g]), as well as transcription factors (T-box 18 [TBX18], insulin gene enhancer binding protein 1 [ISL-1], paired-like homeodomain transcription factor 2 [pitx2], short stature homeobox 2 [shox2]), were detected by western blot and reverse transcription-quantitative polymerase chain reaction. Immunofluorescence assay was carried out to detect whether prrx1 was coexpressed with HCN4, TBX18, and ISL-1. Finally, whole-cell patch-clamp technique was used to record pacing current hyperpolarization-activated inward current (I_f). The isolated cells were CD90⁺, CD29⁺, and CD45^-, indicating that pure BADSCs were successfully isolated. After 5-7 days of Ad transfection into cells, the mRNA levels and protein levels of pacing-related factors (TBX18, ISL-1, HCN4, shox2, and Cacna1g) in Ad-prrx1 group were significantly higher than those in Ad-GFP group. However, the expression level of pitx2 was decreased. Immunofluorescence analysis showed that prrx1 was coexpressed with TBX18, ISL-1, and HCN4 in the Ad-prrx1 group, which did not appear in the Ad-GFP group. Whole-cell patch clamps were able to record the I_f current in the experimental group rather than in the Ad-GFP group. Overexpression of prrx1 can successfully induce sinus node-like cells.

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.

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.

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.

Insulin gene enhancer binding protein 1 induces adipose tissue‑derived stem cells to differentiate into pacemaker‑like cells

Hybrid approaches combining gene- and cell-based therapies to make biological pacemakers are a promising therapeutic avenue for bradyarrhythmia. The present study aimed to direct adipose tissue-derived stem cells (ADSCs) to differentiate specifically into cardiac pacemaker cells by overexpressing a single transcription factor, insulin gene enhancer binding protein 1 (ISL-1). In the present study, the ADSCs were transfected with ISL‑1 or mCherry fluorescent protein lentiviral vectors and co-cultured with neonatal rat ventricular cardiomyocytes (NRVMs) in vitro for 5-7 days. The feasibility of regulating the differentiation of ADSCs into pacemaker-like cells by overexpressing ISL-1 was evaluated by observation of cell morphology and beating rate, reverse transcription-quantitative polymerase chain reaction analysis, western blotting, immunofluorescence and analysis of electrophysiological activity. In conclusion, these data indicated that the overexpression of ISL-1 in ADSCs may enhance the pacemaker phenotype and automaticity in vitro, features which were significantly increased following co‑culture induction.

No effect of thymosin beta-4 on the expression of the transcription factor Islet-1 in the adult murine heart

The transcription factor Islet-1 marks a progenitor cell population of the second heart field during cardiogenesis. In the adult heart Islet-1 expression is limited to the sinoatrial node, the ventricular outflow tract, and parasympathetic ganglia. The regenerative effect in the injured mouse ventricle of thymosin beta-4 (TB4), a 43-aminoacid peptide, was associated with increased Islet-1 immunostaining, suggesting the induction of an Islet-1-positive progenitor state by TB4. Here we aimed to reassess this effect in a genetic model. Mice from the reporter mouse line Isl1-nLacZ were primed with TB4 and subsequently underwent myocardial infarction. Islet-1 expression was assessed 2, 7, and 14 days after infarction. We detected only a single Islet-1+ cell in 8 TB4 treated and infarcted hearts which located outside of the sinoatrial node, the outflow tract or cardiac ganglia (in ~2500 sections). Two cells were identified in 5 control infarcted hearts. TB4 did not induce LacZ positivity in ventricular explants cultures of Isl1-nLacZ mice nor did it affect the density of LacZ+ cells in explant cultures of nLacZ+ regions of the heart. In summary, we found no evidence that TB4 reactivates Islet-1 expression in adult mouse ventricle.

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.

Mutation of the Na+/K+-ATPase Atp1a1a.1 causes QT interval prolongation and bradycardia in zebrafish

The genetic underpinnings that orchestrate the vertebrate heart rate are not fully understood yet, but of high clinical importance, since diseases of cardiac impulse formation and propagation are common and severe human arrhythmias. To identify novel regulators of the vertebrate heart rate, we deciphered the pathogenesis of the bradycardia in the homozygous zebrafish mutant hiphop (hip) and identified a missense-mutation (N851K) in Na⁺/K⁺-ATPase α1-subunit (atp1a1a.1). N851K affects zebrafish Na⁺/K⁺-ATPase ion transport capacity, as revealed by in vitro pump current measurements. Inhibition of the Na⁺/K⁺-ATPase in vivo indicates that hip rather acts as a hypomorph than being a null allele. Consequently, reduced Na⁺/K⁺-ATPase function leads to prolonged QT interval and refractoriness in the hip mutant heart, as shown by electrocardiogram and in vivo electrical stimulation experiments. We here demonstrate for the first time that Na⁺/K⁺-ATPase plays an essential role in heart rate regulation by prolonging myocardial repolarization.

ZebraPace: An Open-Source Method for Cardiac-Rhythm Estimation in Untethered Zebrafish Larvae

For the assessment of cardiac function, heartbeat represents one key parameter. Current methods of heartbeat measurements in the zebrafish larvae usually require larval immobilization, fluorescent transgenic strains and a confocal microscope, costly commercial software for analysis, or strong programming skills if the software is open-source. Here, we present a simple yet powerful method of heartbeat analysis using untethered, unlabeled zebrafish larva using ImageJ (open-source software), which does not require programming skills. We named it as ZebraPace for Zebrafish Precise Algorithm for Cardiac-rhythm Estimation. ZebraPace works directly with AVI videos and requires no image processing steps. ZebraPace uses pixel intensity change in a grayscale video to count the number of beats. We have validated the ZebraPace method by pharmacological alterations of the heartbeat in zebrafish larvae of 48 and 72 hpf stages. We have also determined beat-to-beat interval, which relates to rhythmicity of heartbeat. The results obtained by using ZebraPace corroborates well with the heartbeat values previously reported for similarly aged larvae as determined by using specialized software. We believe that the ZebraPace method is simple, cost-effective, and easy to grasp as it involves fewer steps. It not only reduces the manual workload but also eliminates sample preparation time and researcher subjectivity.

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.

Nkx genes establish second heart field cardiomyocyte progenitors at the arterial pole and pattern the venous pole through Isl1 repression

NKX2-5 is the most commonly mutated gene associated with human congenital heart defects (CHDs), with a predilection for cardiac pole abnormalities. This homeodomain transcription factor is a central regulator of cardiac development and is expressed in both the first and second heart fields (FHF and SHF). We have previously revealed essential functions of nkx2.5 and nkx2.7, two Nkx2-5 homologs expressed in zebrafish cardiomyocytes, in maintaining ventricular identity. However, the differential roles of these genes in the specific subpopulations of the anterior (aSHF) and posterior (pSHF) SHFs have yet to be fully defined. Here, we show that Nkx genes regulate aSHF and pSHF progenitors through independent mechanisms. We demonstrate that Nkx genes restrict proliferation of aSHF progenitors in the outflow tract, delimit the number of pSHF progenitors at the venous pole and pattern the sinoatrial node acting through Isl1 repression. Moreover, optical mapping highlights the requirement for Nkx gene dose in establishing electrophysiological chamber identity and in integrating the physiological connectivity of FHF and SHF cardiomyocytes. Ultimately, our results may shed light on the discrete errors responsible for NKX2-5-dependent human CHDs of the cardiac outflow and inflow tracts.

Gata6 restricts Isl1 to the posterior of nascent hindlimb buds through Isl1 cis-regulatory modules

Isl1 is required for two processes during hindlimb development: initiation of the processes directing hindlimb development in the lateral plate mesoderm and configuring posterior hindlimb field in the nascent hindlimb buds. During these processes, Isl1 expression is restricted to the posterior mesenchyme of hindlimb buds. How this dynamic change in Isl1 expression is regulated remains unknown. We found that two evolutionarily conserved sequences, located 3' to the Isl1 gene, regulate LacZ transgene expression in the hindlimb-forming region in mouse embryos. Both sequences contain GATA binding motifs, and expression pattern analysis identified that Gata6 is expressed in the flank and the anterior portion of nascent hindlimb buds. Recent studies have shown that conditional inactivation of Gata6 in mice causes hindlimb-specific pre-axial polydactyly, indicating a role of Gata6 in anterior-posterior patterning of hindlimbs. We studied whether Gata6 restricts Isl1 in the nascent hindlimb bud through the cis-regulatory modules. In vitro experiments demonstrate that GATA6 binds to the conserved GATA motifs in the cis-regulatory modules. GATA6 repressed expression of a luciferase reporter that contains the cis-regulatory modules by synergizing with Zfpm2. Analyses of Gata6 mutant embryos showed that ISL1 levels are higher in the anterior of nascent hindlimb buds than in wild type. Moreover, we detected a greater number of Isl1-transcribing cells in the anterior of nascent hindlimb buds in Gata6 mutants. Our results support a model in which Gata6 contributes to repression of Isl1 expression in the anterior of nascent hindlimb buds.

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.

Morpho-functional characterization of the systemic venous pole of the reptile heart

Mammals evolved from reptile-like ancestors, and while the mammalian heart is driven by a distinct sinus node, a sinus node is not apparent in reptiles. We characterized the myocardial systemic venous pole, the sinus venosus, in reptiles to identify the dominant pacemaker and to assess whether the sinus venosus remodels and adopts an atrium-like phenotype as observed in mammals. Anolis lizards had an extensive sinus venosus of myocardium expressing Tbx18. A small sub-population of cells encircling the sinuatrial junction expressed Isl1, Bmp2, Tbx3, and Hcn4, homologues of genes marking the mammalian sinus node. Electrical mapping showed that hearts of Anolis lizards and Python snakes were driven from the sinuatrial junction. The electrical impulse was delayed between the sinus venosus and the right atrium, allowing the sinus venosus to contract and aid right atrial filling. In proximity of the systemic veins, the Anolis sinus venosus expressed markers of the atrial phenotype Nkx2-5 and Gja5. In conclusion, the reptile heart is driven by a pacemaker region with an expression signature similar to that of the immature sinus node of mammals. Unlike mammals, reptiles maintain a sinuatrial delay of the impulse, allowing the partly atrialized sinus venosus to function as a chamber.

Intrinsic regulation of sinoatrial node function and the zebrafish as a model of stretch effects on pacemaking

Excitation of the heart occurs in a specialised region known as the sinoatrial node (SAN). Tight regulation of SAN function is essential for the maintenance of normal heart rhythm and the response to (patho-)physiological changes. The SAN is regulated by extrinsic (central nervous system) and intrinsic (neurons, peptides, mechanics) factors. The positive chronotropic response to stretch in particular is essential for beat-by-beat adaptation to changes in hemodynamic load. Yet, the mechanism of this stretch response is unknown, due in part to the lack of an appropriate experimental model for targeted investigations. We have been investigating the zebrafish as a model for the study of intrinsic regulation of SAN function. In this paper, we first briefly review current knowledge of the principal components of extrinsic and intrinsic SAN regulation, derived primarily from experiments in mammals, followed by a description of the zebrafish as a novel experimental model for studies of intrinsic SAN regulation. This mini-review is followed by an original investigation of the response of the zebrafish isolated SAN to controlled stretch. Stretch causes an immediate and continuous increase in beating rate in the zebrafish isolated SAN. This increase reaches a maximum part way through a period of sustained stretch, with the total change dependent on the magnitude and direction of stretch. This is comparable to what occurs in isolated SAN from most mammals (including human), suggesting that the zebrafish is a novel experimental model for the study of mechanisms involved in the intrinsic regulation of SAN function by mechanical effects.

Development and Function of the Cardiac Conduction System in Health and Disease

The generation and propagation of the cardiac impulse is the central function of the cardiac conduction system (CCS). Impulse initiation occurs in nodal tissues that have high levels of automaticity, but slow conduction properties. Rapid impulse propagation is a feature of the ventricular conduction system, which is essential for synchronized contraction of the ventricular chambers. When functioning properly, the CCS produces ~2.4 billion heartbeats during a human lifetime and orchestrates the flow of cardiac impulses, designed to maximize cardiac output. Abnormal impulse initiation or propagation can result in brady- and tachy-arrhythmias, producing an array of symptoms, including syncope, heart failure or sudden cardiac death. Underlying the functional diversity of the CCS are gene regulatory networks that direct cell fate towards a nodal or a fast conduction gene program. In this review, we will discuss our current understanding of the transcriptional networks that dictate the components of the CCS, the growth factor-dependent signaling pathways that orchestrate some of these transcriptional hierarchies and the effect of aberrant transcription factor expression on mammalian conduction disease.

On the Evolution of the Cardiac Pacemaker

The rhythmic contraction of the heart is initiated and controlled by an intrinsic pacemaker system. Cardiac contractions commence at very early embryonic stages and coordination remains crucial for survival. The underlying molecular mechanisms of pacemaker cell development and function are still not fully understood. Heart form and function show high evolutionary conservation. Even in simple contractile cardiac tubes in primitive invertebrates, cardiac function is controlled by intrinsic, autonomous pacemaker cells. Understanding the evolutionary origin and development of cardiac pacemaker cells will help us outline the important pathways and factors involved. Key patterning factors, such as the homeodomain transcription factors Nkx2.5 and Shox2, and the LIM-homeodomain transcription factor Islet-1, components of the T-box (Tbx), and bone morphogenic protein (Bmp) families are well conserved. Here we compare the dominant pacemaking systems in various organisms with respect to the underlying molecular regulation. Comparative analysis of the pathways involved in patterning the pacemaker domain in an evolutionary context might help us outline a common fundamental pacemaker cell gene programme. Special focus is given to pacemaker development in zebrafish, an extensively used model for vertebrate development. Finally, we conclude with a summary of highly conserved key factors in pacemaker cell development and function.

Comparative expression analysis of Shox2-deficient embryonic stem cell-derived sinoatrial node-like cells

The homeodomain transcription factor Shox2 controls the development and function of the native cardiac pacemaker, the sinoatrial node (SAN). Moreover, SHOX2 mutations have been associated with cardiac arrhythmias in humans. For detailed examination of Shox2-dependent developmental mechanisms in SAN cells, we established a murine embryonic stem cell (ESC)-based model using Shox2 as a molecular tool. Shox2^+/+ and Shox2^-/- ESC clones were isolated and differentiated according to five different protocols in order to evaluate the most efficient enrichment of SAN-like cells. Expression analysis of cell subtype-specific marker genes revealed most efficient enrichment after CD166-based cell sorting. Comparative cardiac expression profiles of Shox2^+/+ and Shox2^-/- ESCs were examined by nCounter technology. Among other genes, we identified Nppb as a novel putative Shox2 target during differentiation in ESCs. Differential expression of Nppb could be confirmed in heart tissue of Shox2^-/- embryos. Taken together, we established an ESC-based cardiac differentiation model and successfully purified Shox2^+/+ and Shox2^-/- SAN-like cells. This now provides an excellent basis for the investigation of molecular mechanisms under physiological and pathophysiological conditions for evaluating novel therapeutic approaches.

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.

Development of the cardiac conduction system in zebrafish

The cardiac conduction system (CCS) propagates and coordinates the electrical excitation that originates from the pacemaker cells, throughout the heart, resulting in rhythmic heartbeat. Its defects result in life-threatening arrhythmias and sudden cardiac death. Understanding of the factors involved in the formation and function of the CCS remains incomplete. By transposon assisted transgenesis, we have developed enhancer trap (ET) lines of zebrafish that express fluorescent protein in the pacemaker cells at the sino-atrial node (SAN) and the atrio-ventricular region (AVR), termed CCS transgenics. This expression pattern begins at the stage when the heart undergoes looping morphogenesis at 36 h post fertilization (hpf) and is maintained into adulthood. Using the CCS transgenics, we investigated the effects of perturbation of cardiac function, as simulated by either the absence of endothelium or hemodynamic stimulation, on the cardiac conduction cells, which resulted in abnormal compaction of the SAN. To uncover the identity of the gene represented by the EGFP expression in the CCS transgenics, we mapped the transposon integration sites on the zebrafish genome to positions in close proximity to the gene encoding fibroblast growth homologous factor 2a (fhf2a). Fhf2a is represented by three transcripts, one of which is expressed in the developing heart. These transgenics are useful tools for studies of development of the CCS and cardiac disease.

The formation and function of the cardiac conduction system

The cardiac conduction system (CCS) consists of distinctive components that initiate and conduct the electrical impulse required for the coordinated contraction of the cardiac chambers. CCS development involves complex regulatory networks that act in stage-, tissue- and dose-dependent manners, and recent findings indicate that the activity of these networks is sensitive to common genetic variants associated with cardiac arrhythmias. Here, we review how these findings have provided novel insights into the regulatory mechanisms and transcriptional networks underlying CCS formation and function.

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.

Direct nkx2-5 transcriptional repression of isl1 controls cardiomyocyte subtype identity

During cardiogenesis, most myocytes arise from cardiac progenitors expressing the transcription factors Isl1 and Nkx2-5. Here, we show that a direct repression of Isl1 by Nkx2-5 is necessary for proper development of the ventricular myocardial lineage. Overexpression of Nkx2-5 in mouse embryonic stem cells (ESCs) delayed specification of cardiac progenitors and inhibited expression of Isl1 and its downstream targets in Isl1(+) precursors. Embryos deficient for Nkx2-5 in the Isl1(+) lineage failed to downregulate Isl1 protein in cardiomyocytes of the heart tube. We demonstrated that Nkx2-5 directly binds to an Isl1 enhancer and represses Isl1 transcriptional activity. Furthermore, we showed that overexpression of Isl1 does not prevent cardiac differentiation of ESCs and in Xenopus laevis embryos. Instead, it leads to enhanced specification of cardiac progenitors, earlier cardiac differentiation, and increased cardiomyocyte number. Functional and molecular characterization of Isl1-overexpressing cardiomyocytes revealed higher beating frequencies in both ESC-derived contracting areas and Xenopus Isl1-gain-of-function hearts, which associated with upregulation of nodal-specific genes and downregulation of transcripts of working myocardium. Immunocytochemistry of cardiomyocyte lineage-specific markers demonstrated a reduction of ventricular cells and an increase of cells expressing the pacemaker channel Hcn4. Finally, optical action potential imaging of single cardiomyocytes combined with pharmacological approaches proved that Isl1 overexpression in ESCs resulted in normally electrophysiologically functional cells, highly enriched in the nodal subtype at the expense of the ventricular lineage. Our findings provide an Isl1/Nkx2-5-mediated mechanism that coordinately regulates the specification of cardiac progenitors toward the different myocardial lineages and ensures proper acquisition of myocyte subtype identity.

Functional Diversification of Motor Neuron-specific Isl1 Enhancers during Evolution

Functional diversification of motor neurons has occurred in order to selectively control the movements of different body parts including head, trunk and limbs. Here we report that transcription of Isl1, a major gene necessary for motor neuron identity, is controlled by two enhancers, CREST1 (E1) and CREST2 (E2) that allow selective gene expression of Isl1 in motor neurons. Introduction of GFP reporters into the chick neural tube revealed that E1 is active in hindbrain motor neurons and spinal cord motor neurons, whereas E2 is active in the lateral motor column (LMC) of the spinal cord, which controls the limb muscles. Genome-wide ChIP-Seq analysis combined with reporter assays showed that Phox2 and the Isl1-Lhx3 complex bind to E1 and drive hindbrain and spinal cord-specific expression of Isl1, respectively. Interestingly, Lhx3 alone was sufficient to activate E1, and this may contribute to the initiation of Isl1 expression when progenitors have just developed into motor neurons. E2 was induced by onecut 1 (OC-1) factor that permits Isl1 expression in LMCm neurons. Interestingly, the core region of E1 has been conserved in evolution, even in the lamprey, a jawless vertebrate with primitive motor neurons. All E1 sequences from lamprey to mouse responded equally well to Phox2a and the Isl1-Lhx3 complex. Conversely, E2, the enhancer for limb-innervating motor neurons, was only found in tetrapod animals. This suggests that evolutionarily-conserved enhancers permit the diversification of motor neurons.

During evolution, motor neurons became specialized to control movements of different body parts including head, trunk and limbs. Here we report that two enhancers of Isl1, E1 and E2, are active together with transcription factors in motor neurons. Surprisingly, E1 and its response to transcription factors has been conserved in evolution from the lamprey to man, whereas E2 is only found in animals with limbs. Our study provides an evolutionary example of how functional diversification of motor neurons is achieved by a dynamic interplay between enhancers and transcription factors.

Transcription factor ISL1 is essential for pacemaker development and function

The sinoatrial node (SAN) maintains a rhythmic heartbeat; therefore, a better understanding of factors that drive SAN development and function is crucial to generation of potential therapies, such as biological pacemakers, for sinus arrhythmias. Here, we determined that the LIM homeodomain transcription factor ISL1 plays a key role in survival, proliferation, and function of pacemaker cells throughout development. Analysis of several Isl1 mutant mouse lines, including animals harboring an SAN-specific Isl1 deletion, revealed that ISL1 within SAN is a requirement for early embryonic viability. RNA-sequencing (RNA-seq) analyses of FACS-purified cells from ISL1-deficient SANs revealed that a number of genes critical for SAN function, including those encoding transcription factors and ion channels, were downstream of ISL1. Chromatin immunoprecipitation assays performed with anti-ISL1 antibodies and chromatin extracts from FACS-purified SAN cells demonstrated that ISL1 directly binds genomic regions within several genes required for normal pacemaker function, including subunits of the L-type calcium channel, Ank2, and Tbx3. Other genes implicated in abnormal heart rhythm in humans were also direct ISL1 targets. Together, our results demonstrate that ISL1 regulates approximately one-third of SAN-specific genes, indicate that a combination of ISL1 and other SAN transcription factors could be utilized to generate pacemaker cells, and suggest ISL1 mutations may underlie sick sinus syndrome.

Shox2 is required for the proper development of the facial motor nucleus and the establishment of the facial nerves

Background

Axons from the visceral motor neurons (vMNs) project from nuclei in the hindbrain to innervate autonomic ganglia and branchial arch-derived muscles. Although much is known about the events that govern specification of somatic motor neurons, the genetic pathways responsible for the development of vMNs are less well characterized. We know that vMNs, like all motor neurons, depend on sonic hedgehog signaling for their generation. Similarly, the paired-like homeobox 2b (Phox2b) gene, which is expressed in both proliferating progenitors and post-mitotic motor neurons, is essential for the development of vMNs. Given that our previous study identified a novel role for the short stature homeobox 2 (Shox2) gene in the hindbrain, and since SHOX2 has been shown to regulate transcription of islet 1 (Isl1), an important regulator of vMN development, we sought to determine whether Shox2 is required for the proper development of the facial motor nucleus.

Results

Using a Nestin-Cre driver, we show that elimination of Shox2 throughout the brain results in elevated cell death in the facial motor nucleus at embryonic day 12.5 (E12.5) and E14.5, which correlates with impaired axonal projection properties of vMNs. We also observed changes in the spatial expression of the vMN cell fate factors Isl1 and Phox2b, and concomitant defects in Shh and Ptch1 expression in Shox2 mutants. Furthermore, we demonstrate that elimination of Shox2 results in the loss of dorsomedial and ventromedial subnuclei by postnatal day 0 (P0), which may explain the changes in physical activity and impaired feeding/nursing behavior in Shox2 mutants.

Conclusions

Combined, our data show that Shox2 is required for development of the facial motor nucleus and its associated facial (VII) nerves, and serves as a new molecular tool to probe the genetic programs of this understudied hindbrain region.

Electronic supplementary material

The online version of this article (doi:10.1186/s12868-015-0176-0) contains supplementary material, which is available to authorized users.

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

In humans, atrial fibrillation is often triggered by ectopic pacemaking activity in the myocardium sleeves of the pulmonary vein (PV) and systemic venous return. The genetic programs that abnormally reinforce pacemaker properties at these sites and how this relates to normal sinoatrial node (SAN) development remain uncharacterized. It was noted previously that Nkx2-5, which is expressed in the PV myocardium and reinforces a chamber-like myocardial identity in the PV, is lacking in the SAN. Here we present evidence that in mice Shox2 antagonizes the transcriptional output of Nkx2-5 in the PV myocardium and in a functional Nkx2-5(+) domain within the SAN to determine cell fate. Shox2 deletion in the Nkx2-5(+) domain of the SAN caused sick sinus syndrome, associated with the loss of the pacemaker program. Explanted Shox2(+) cells from the embryonic PV myocardium exhibited pacemaker characteristics including node-like electrophysiological properties and the capability to pace surrounding Shox2(-) cells. Shox2 deletion led to Hcn4 ablation in the developing PV myocardium. Nkx2-5 hypomorphism rescued the requirement for Shox2 for the expression of genes essential for SAN development in Shox2 mutants. Similarly, the pacemaker-like phenotype induced in the PV myocardium in Nkx2-5 hypomorphs reverted back to a working myocardial phenotype when Shox2 was simultaneously deleted. A similar mechanism is also adopted in differentiated embryoid bodies. We found that Shox2 interacts with Nkx2-5 directly, and discovered a substantial genome-wide co-occupancy of Shox2, Nkx2-5 and Tbx5, further supporting a pivotal role for Shox2 in the core myogenic program orchestrating venous pole and pacemaker development.