TBX5 drives Scn5a expression to regulate cardiac conduction system function

14-3-3ε/YWHAE regulates the transcriptional expression of cardiac sodium channel NaV1.5

BACKGROUND

Cardiac voltage-gated sodium channel alpha subunit 5 (NaV1.5) encoded by SCN5A is associated with arrhythmia disorders. However, the molecular mechanism underlying NaV1.5 expression remains to be fully elucidated. Previous studies have reported that the 14-3-3 family acts as an adaptor involved in regulating kinetic characteristics of cardiac ion channels.

OBJECTIVE

The purpose of this study was to establish 14-3-3ε/YWHAE, a member of the 14-3-3 family, as a crucial regulator of NaV1.5 and to explore the potential role of 14-3-3ε in the heart.

METHODS

Western blotting, patch clamping, real-time reverse transcription-polymerase chain reaction, RNA immunoprecipitation, electrocardiogram recording, echocardiography, and histologic analysis were performed.

RESULTS

YWHAE overexpression significantly reduced the expression level of SCN5A mRNA and sodium current density, whereas YWHAE knockdown significantly increased SCN5A mRNA expression and sodium current density in HEK293/NaV1.5 and H9c2 cells. Similar results were observed in mice injected with adeno-associated virus serotype 9-mediated YWHAE knockdown. The effect of 14-3-3ε on NaV1.5 expression was abrogated by knockdown of TBX5, a transcription factor of NaV1.5. An interaction between 14-3-3ε protein and TBX5 mRNA was identified, and YWHAE overexpression significantly decreased TBX5 mRNA stability without affecting SCN5A mRNA stability. In addition, mice subjected to adeno-associated virus serotype 9-mediated YWHAE knockdown exhibited shorter R-R intervals and higher prevalence of premature ventricular contractions.

CONCLUSION

Our data unveil a novel regulatory mechanism of NaV1.5 by 14-3-3ε and highlight the significance of 14-3-3ε in transcriptional regulation of NaV1.5 expression and cardiac arrhythmias.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network

BACKGROUND

The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice.

METHODS

FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription.

RESULTS

FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency.

CONCLUSIONS

Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

TBX5 genetic variants and SCD-CAD susceptibility: insights from Chinese Han cohorts

Background

The prevention and prediction of sudden cardiac death (SCD) present persistent challenges, prompting exploration into common genetic variations for potential insights. T-box 5 (TBX5), a critical cardiac transcription factor, plays a pivotal role in cardiovascular development and function. This study systematically examined variants within the 500-bp region downstream of the TBX5 gene, focusing on their potential impact on susceptibility to SCD associated with coronary artery disease (SCD-CAD) in four different Chinese Han populations.

Methods

In a comprehensive case-control analysis, we explored the association between rs11278315 and SCD-CAD susceptibility using a cohort of 553 controls and 201 SCD-CAD cases. Dual luciferase reporter assays and genotype-phenotype correlation studies using human cardiac tissue samples as well as integrated in silicon analysis were applied to explore the underlining mechanism.

Result

Binary logistic regression results underscored a significantly reduced risk of SCD-CAD in individuals harboring the deletion allele (odds ratio = 0.70, 95% CI [0.55-0.88], p = 0.0019). Consistent with the lower transcriptional activity of the deletion allele observed in dual luciferase reporter assays, genotype-phenotype correlation studies on human cardiac tissue samples affirmed lower expression levels associated with the deletion allele at both mRNA and protein levels. Furthermore, our investigation revealed intriguing insights into the role of rs11278315 in TBX5 alternative splicing, which may contribute to alterations in its ultimate functional effects, as suggested by sQTL analysis. Gene ontology analysis and functional annotation further underscored the potential involvement of TBX5 in alternative splicing and cardiac-related transcriptional regulation.

Conclusions

In summary, our current dataset points to a plausible correlation between rs11278315 and susceptibility to SCD-CAD, emphasizing the potential of rs11278315 as a genetic risk marker for aiding in molecular diagnosis and risk stratification of SCD-CAD.

Hypoxia Attenuates Pressure Overload-Induced Heart Failure

BACKGROUND

Alveolar hypoxia is protective in the context of cardiovascular and ischemic heart disease; however, the underlying mechanisms are incompletely understood. The present study sought to test the hypothesis that hypoxia is cardioprotective in left ventricular pressure overload (LVPO)-induced heart failure. We furthermore aimed to test that overlapping mechanisms promote cardiac recovery in heart failure patients following left ventricular assist device-mediated mechanical unloading and circulatory support.

METHODS AND RESULTS

We established a novel murine model of combined chronic alveolar hypoxia and LVPO following transverse aortic constriction (HxTAC). The HxTAC model is resistant to cardiac hypertrophy and the development of heart failure. The cardioprotective mechanisms identified in our HxTAC model include increased activation of HIF (hypoxia-inducible factor)-1α-mediated angiogenesis, attenuated induction of genes associated with pathological remodeling, and preserved metabolic gene expression as identified by RNA sequencing. Furthermore, LVPO decreased Tbx5 and increased Hsd11b1 mRNA expression under normoxic conditions, which was attenuated under hypoxic conditions and may induce additional hypoxia-mediated cardioprotective effects. Analysis of samples from patients with advanced heart failure that demonstrated left ventricular assist device-mediated myocardial recovery revealed a similar expression pattern for TBX5 and HSD11B1 as observed in HxTAC hearts.

CONCLUSIONS

Hypoxia attenuates LVPO-induced heart failure. Cardioprotective pathways identified in the HxTAC model might also contribute to cardiac recovery following left ventricular assist device support. These data highlight the potential of our novel HxTAC model to identify hypoxia-mediated cardioprotective mechanisms and therapeutic targets that attenuate LVPO-induced heart failure and mediate cardiac recovery following mechanical circulatory support.

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.

Tbx5 maintains atrial identity in post-natal cardiomyocytes by regulating an atrial-specific enhancer network

Understanding how the atrial and ventricular heart chambers maintain distinct identities is a prerequisite for treating chamber-specific diseases. Here, we selectively knocked out (KO) the transcription factor Tbx5 in the atrial working myocardium to evaluate its requirement for atrial identity. Atrial Tbx5 inactivation downregulated atrial cardiomyocyte (aCM) selective gene expression. Using concurrent single nucleus transcriptome and open chromatin profiling, genomic accessibility differences were identified between control and Tbx5 KO aCMs, revealing that 69% of the control-enriched ATAC regions were bound by TBX5. Genes associated with these regions were downregulated in KO aCMs, suggesting they function as TBX5-dependent enhancers. Comparing enhancer chromatin looping using H3K27ac HiChIP identified 510 chromatin loops sensitive to TBX5 dosage, and 74.8% of control-enriched loops contained anchors in control-enriched ATAC regions. Together, these data demonstrate TBX5 maintains the atrial gene expression program by binding to and preserving the tissue-specific chromatin architecture of atrial enhancers.

Tbx5 overexpression in embryoid bodies increases TAK1 expression but does not enhance the differentiation of sinoatrial node cardiomyocytes

Genetic studies place Tbx5 at the apex of the sinoatrial node (SAN) transcriptional program. To understand its role in SAN differentiation, clonal embryonic stem (ES) cell lines were made that conditionally overexpress Tbx5, Tbx3, Tbx18, Shox2, Islet-1, and MAP3k7/TAK1. Cardiac cells differentiated using embryoid bodies (EBs). EBs overexpressing Tbx5, Islet1, and TAK1 beat faster than cardiac cells differentiated from control ES cell lines, suggesting possible roles in SAN differentiation. Tbx5 overexpressing EBs showed increased expression of TAK1, but cardiomyocytes did not differentiate as SAN cells. EBs showed no change in the expression of the SAN transcription factors Shox2 and Islet1 and decreased expression of the SAN channel protein HCN4. EBs constitutively overexpressing TAK1 direct cardiac differentiation to the SAN fate but have reduced phosphorylation of its targets, p38 and Jnk. This opens the possibility that blocking the phosphorylation of TAK1 targets may have the same impact as forced overexpression. To test this, we treated EBs with 5z-7-Oxozeanol (OXO), an inhibitor of TAK1 phosphorylation. Like TAK1 overexpressing cardiac cells, cardiomyocytes differentiated in the presence of OXO beat faster and showed increased expression of SAN genes (Shox2, HCN4, and Islet1). This suggests that activation of the SAN transcriptional network can be accomplished by blocking the phosphorylation of TAK1.

Tbx5 maintains atrial identity by regulating an atrial enhancer network

Understanding how the atrial and ventricular chambers of the heart maintain their distinct identity is a prerequisite for treating chamber-specific diseases. Here, we selectively inactivated the transcription factor Tbx5 in the atrial working myocardium of the neonatal mouse heart to show that it is required to maintain atrial identity. Atrial Tbx5 inactivation downregulated highly chamber specific genes such as Myl7 and Nppa , and conversely, increased the expression of ventricular identity genes including Myl2 . Using combined single nucleus transcriptome and open chromatin profiling, we assessed genomic accessibility changes underlying the altered atrial identity expression program, identifying 1846 genomic loci with greater accessibility in control atrial cardiomyocytes compared to KO aCMs. 69% of the control-enriched ATAC regions were bound by TBX5, demonstrating a role for TBX5 in maintaining atrial genomic accessibility. These regions were associated with genes that had higher expression in control aCMs compared to KO aCMs, suggesting they act as TBX5-dependent enhancers. We tested this hypothesis by analyzing enhancer chromatin looping using HiChIP and found 510 chromatin loops that were sensitive to TBX5 dosage. Of the loops enriched in control aCMs, 73.7% contained anchors in control-enriched ATAC regions. Together, these data demonstrate a genomic role for TBX5 in maintaining the atrial gene expression program by binding to atrial enhancers and preserving tissue-specific chromatin architecture of atrial enhancers.

Return of the Tbx5; lineage-tracing reveals ventricular cardiomyocyte-like precursors in the injured adult mammalian heart

The single curative measure for heart failure patients is a heart transplantation, which is limited due to a shortage of donors, the need for immunosuppression and economic costs. Therefore, there is an urgent unmet need for identifying cell populations capable of cardiac regeneration that we will be able to trace and monitor. Injury to the adult mammalian cardiac muscle, often leads to a heart attack through the irreversible loss of a large number of cardiomyocytes, due to an idle regenerative capability. Recent reports in zebrafish indicate that Tbx5a is a vital transcription factor for cardiomyocyte regeneration. Preclinical data underscore the cardioprotective role of Tbx5 upon heart failure. Data from our earlier murine developmental studies have identified a prominent unipotent Tbx5-expressing embryonic cardiac precursor cell population able to form cardiomyocytes, in vivo, in vitro and ex vivo. Using a developmental approach to an adult heart injury model and by employing a lineage-tracing mouse model as well as the use of single-cell RNA-seq technology, we identify a Tbx5-expressing ventricular cardiomyocyte-like precursor population, in the injured adult mammalian heart. The transcriptional profile of that precursor cell population is closer to that of neonatal than embryonic cardiomyocyte precursors. Tbx5, a cardinal cardiac development transcription factor, lies in the center of a ventricular adult precursor cell population, which seems to be affected by neurohormonal spatiotemporal cues. The identification of a Tbx5-specific cardiomyocyte precursor-like cell population, which is capable of dedifferentiating and potentially deploying a cardiomyocyte regenerative program, provides a clear target cell population for translationally-relevant heart interventional studies.

Sinus node dysfunction: current understanding and future directions

The sinoatrial node (SAN) is the primary pacemaker of the heart. Normal SAN function is crucial in maintaining proper cardiac rhythm and contraction. Sinus node dysfunction (SND) is due to abnormalities within the SAN, which can affect the heartbeat frequency, regularity, and the propagation of electrical pulses through the cardiac conduction system. As a result, SND often increases the risk of cardiac arrhythmias. SND is most commonly seen as a disease of the elderly given the role of degenerative fibrosis as well as other age-dependent changes in its pathogenesis. Despite the prevalence of SND, current treatment is limited to pacemaker implantation, which is associated with substantial medical costs and complications. Emerging evidence has identified various genetic abnormalities that can cause SND, shedding light on the molecular underpinnings of SND. Identification of these molecular mechanisms and pathways implicated in the pathogenesis of SND is hoped to identify novel therapeutic targets for the development of more effective therapies for this disease. In this review article, we examine the anatomy of the SAN and the pathophysiology and epidemiology of SND. We then discuss in detail the most common genetic mutations correlated with SND and provide our perspectives on future research and therapeutic opportunities in this field.

Global chromatin landscapes identify candidate noncoding modifiers of cardiac rhythm

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

Deciphering Transcriptional Networks during Human Cardiac Development

Human heart development is governed by transcription factor (TF) networks controlling dynamic and temporal gene expression alterations. Therefore, to comprehensively characterize these transcriptional regulations, day-to-day transcriptomic profiles were generated throughout the directed cardiac differentiation, starting from three distinct human- induced pluripotent stem cell lines from healthy donors (32 days). We applied an expression-based correlation score to the chronological expression profiles of the TF genes, and clustered them into 12 sequential gene expression waves. We then identified a regulatory network of more than 23,000 activation and inhibition links between 216 TFs. Within this network, we observed previously unknown inferred transcriptional activations linking IRX3 and IRX5 TFs to three master cardiac TFs: GATA4, NKX2-5 and TBX5. Luciferase and co-immunoprecipitation assays demonstrated that these five TFs could (1) activate each other's expression; (2) interact physically as multiprotein complexes; and (3) together, finely regulate the expression of SCN5A, encoding the major cardiac sodium channel. Altogether, these results unveiled thousands of interactions between TFs, generating multiple robust hypotheses governing human cardiac development.

Enhancer RNA (eRNA) in Human Diseases

Enhancer RNAs (eRNAs), a class of non-coding RNAs (ncRNAs) transcribed from enhancer regions, serve as a type of critical regulatory element in gene expression. There is increasing evidence demonstrating that the aberrant expression of eRNAs can be broadly detected in various human diseases. Some studies also revealed the potential clinical utility of eRNAs in these diseases. In this review, we summarized the recent studies regarding the pathological mechanisms of eRNAs as well as their potential utility across human diseases, including cancers, neurodegenerative disorders, cardiovascular diseases and metabolic diseases. It could help us to understand how eRNAs are engaged in the processes of diseases and to obtain better insight of eRNAs in diagnosis, prognosis or therapy. The studies we reviewed here indicate the enormous therapeutic potency of eRNAs across human diseases.

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

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

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.

The role of the cardiomyocyte circadian clocks in ion channel regulation and cardiac electrophysiology

Daily variations in cardiac electrophysiology and the incidence for different types of arrhythmias reflect ≈24 h changes in the environment, behaviour and internal circadian rhythms. This article focuses on studies that use animal models to separate the impact that circadian rhythms, as well as changes in the environment and behaviour, have on 24 h rhythms in heart rate and ventricular repolarization. Circadian rhythms are initiated at the cellular level by circadian clocks, transcription-translation feedback loops that cycle with a periodicity of 24 h. Several studies now show that the circadian clock in cardiomyocytes regulates the expression of cardiac ion channels by multiple mechanisms; underlies time-of-day changes in sinoatrial node excitability/intrinsic heart rate; and limits the duration of the ventricular action potential waveform. However, the 24 h rhythms in heart rate and ventricular repolarization are primarily driven by autonomic signalling. A functional role for the cardiomyocyte circadian clock appears to buffer the heart against perturbations. For example, the cardiomyocyte circadian clock limits QT-interval prolongation (especially at slower heart rates), and it may facilitate the realignment of the 24 h rhythm in heart rate to abrupt changes in the light cycle. Additional studies show that modifying rhythmic behaviours (including feeding behaviour) can dramatically impact the 24 h rhythms in heart rate and ventricular repolarization. If these mechanisms are conserved, these studies suggest that targeting endogenous circadian mechanisms in the heart, as well as modifying the timing of certain rhythmic behaviours, could emerge as therapeutic strategies to support heart function against perturbations and regulate 24 h rhythms in cardiac electrophysiology.

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

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

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

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

Genomic and Non-Genomic Regulatory Mechanisms of the Cardiac Sodium Channel in Cardiac Arrhythmias

Nav1.5 is the predominant cardiac sodium channel subtype, encoded by the SCN5A gene, which is involved in the initiation and conduction of action potentials throughout the heart. Along its biosynthesis process, Nav1.5 undergoes strict genomic and non-genomic regulatory and quality control steps that allow only newly synthesized channels to reach their final membrane destination and carry out their electrophysiological role. These regulatory pathways are ensured by distinct interacting proteins that accompany the nascent Nav1.5 protein along with different subcellular organelles. Defects on a large number of these pathways have a tremendous impact on Nav1.5 functionality and are thus intimately linked to cardiac arrhythmias. In the present review, we provide current state-of-the-art information on the molecular events that regulate SCN5A/Nav1.5 and the cardiac channelopathies associated with defects in these pathways.

Genomic enhancers in cardiac development and disease

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

MicroRNA is a potential target for therapies to improve the physiological function of skeletal muscle after trauma

MicroRNAs can regulate the function of ion channels in many organs. Based on our previous study we propose that miR-142a-39, which is highly expressed in denervated skeletal muscle, might affect cell excitability through similar mechanisms. In this study, we overexpressed or knocked down miR-142a-3p in C2C12 cells using a lentivirus method. After 7 days of differentiation culture, whole-cell currents were recorded. The results showed that overexpression of miR-142a-3p reduced the cell membrane capacitance, increased potassium current density and decreased calcium current density. Knockdown of miR-142a-3p reduced sodium ion channel current density. The results showed that change in miR-142a-3p expression affected the ion channel currents in C2C12 cells, suggesting its possible roles in muscle cell electrophysiology. This study was approved by the Animal Ethics Committee of Peking University in July 2020 (approval No. LA2017128).

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.

Genome-wide association studies of cardiac electrical phenotypes

The genetic basis of cardiac electrical phenotypes has in the last 25 years been the subject of intense investigation. While in the first years, such efforts were dominated by the study of familial arrhythmia syndromes, in recent years, large consortia of investigators have successfully pursued genome-wide association studies (GWAS) for the identification of single-nucleotide polymorphisms that govern inter-individual variability in electrocardiographic parameters in the general population. We here provide a review of GWAS conducted on cardiac electrical phenotypes in the last 14 years and discuss the implications of these discoveries for our understanding of the genetic basis of disease susceptibility and variability in disease severity. Furthermore, we review functional follow-up studies that have been conducted on GWAS loci associated with cardiac electrical phenotypes and highlight the challenges and opportunities offered by such studies.

Preclinical evidence for the therapeutic value of TBX5 normalization in arrhythmia control

AIMS

Arrhythmias and sudden cardiac death (SCD) occur commonly in patients with heart failure. We found T-box 5 (TBX5) dysregulated in ventricular myocardium from heart failure patients and thus we hypothesized that TBX5 reduction contributes to arrhythmia development in these patients. To understand the underlying mechanisms, we aimed to reveal the ventricular TBX5-dependent transcriptional network and further test the therapeutic potential of TBX5 level normalization in mice with documented arrhythmias.

METHODS AND RESULTS

We used a mouse model of TBX5 conditional deletion in ventricular cardiomyocytes. Ventricular (v) TBX5 loss in mice resulted in mild cardiac dysfunction and arrhythmias and was associated with a high mortality rate (60%) due to SCD. Upon angiotensin stimulation, vTbx5KO mice showed exacerbated cardiac remodelling and dysfunction suggesting a cardioprotective role of TBX5. RNA-sequencing of a ventricular-specific TBX5KO mouse and TBX5 chromatin immunoprecipitation was used to dissect TBX5 transcriptional network in cardiac ventricular tissue. Overall, we identified 47 transcripts expressed under the control of TBX5, which may have contributed to the fatal arrhythmias in vTbx5KO mice. These included transcripts encoding for proteins implicated in cardiac conduction and contraction (Gja1, Kcnj5, Kcng2, Cacna1g, Chrm2), in cytoskeleton organization (Fstl4, Pdlim4, Emilin2, Cmya5), and cardiac protection upon stress (Fhl2, Gpr22, Fgf16). Interestingly, after TBX5 loss and arrhythmia development in vTbx5KO mice, TBX5 protein-level normalization by systemic adeno-associated-virus (AAV) 9 application, re-established TBX5-dependent transcriptome. Consequently, cardiac dysfunction was ameliorated and the propensity of arrhythmia occurrence was reduced.

CONCLUSIONS

This study uncovers a novel cardioprotective role of TBX5 in the adult heart and provides preclinical evidence for the therapeutic value of TBX5 protein normalization in the control of arrhythmia.

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

PURPOSE OF REVIEW

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

RECENT FINDINGS

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

Cardiomyocyte Deletion of Bmal1 Exacerbates QT- and RR-Interval Prolongation in Scn5a +/ΔKPQ Mice

Circadian rhythms are generated by cell autonomous circadian clocks that perform a ubiquitous cellular time-keeping function and cell type-specific functions important for normal physiology. Studies show inducing the deletion of the core circadian clock transcription factor Bmal1 in adult mouse cardiomyocytes disrupts cardiac circadian clock function, cardiac ion channel expression, slows heart rate, and prolongs the QT-interval at slow heart rates. This study determined how inducing the deletion of Bmal1 in adult cardiomyocytes impacted the in vivo electrophysiological phenotype of a knock-in mouse model for the arrhythmogenic long QT syndrome (Scn5a^+/ΔKPQ). Electrocardiographic telemetry showed inducing the deletion of Bmal1 in the cardiomyocytes of mice with or without the ΔKPQ-Scn5a mutation increased the QT-interval at RR-intervals that were ≥130 ms. Inducing the deletion of Bmal1 in the cardiomyocytes of mice with or without the ΔKPQ-Scn5a mutation also increased the day/night rhythm-adjusted mean in the RR-interval, but it did not change the period, phase or amplitude. Compared to mice without the ΔKPQ-Scn5a mutation, mice with the ΔKPQ-Scn5a mutation had reduced heart rate variability (HRV) during the peak of the day/night rhythm in the RR-interval. Inducing the deletion of Bmal1 in cardiomyocytes did not affect HRV in mice without the ΔKPQ-Scn5a mutation, but it did increase HRV in mice with the ΔKPQ-Scn5a mutation. The data demonstrate that deleting Bmal1 in cardiomyocytes exacerbates QT- and RR-interval prolongation in mice with the ΔKPQ-Scn5a mutation.

Tbx5 variants disrupt Nav1.5 function differently in patients diagnosed with Brugada or Long QT Syndrome

AIMS

The transcription factor Tbx5 controls cardiogenesis and drives Scn5a expression in mice. We have identified two variants in TBX5 encoding p.D111Y and p.F206L Tbx5, respectively, in two unrelated patients with structurally normal hearts diagnosed with Long QT (LQTS) and Brugada (BrS) Syndrome. Here we characterized the consequences of each variant to unravel the underlying disease mechanisms.

METHODS AND RESULTS

We combined clinical analysis with in vivo and in vitro electrophysiological and molecular techniques in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), HL-1 cells, and cardiomyocytes from mice trans-expressing human wildtype (WT) or mutant proteins. Tbx5 increased transcription of SCN5A encoding cardiac Nav1.5 channels, while repressing CAMK2D and SPTBN4 genes encoding Ca-calmodulin kinase IIδ (CaMKIIδ) and βIV-spectrin, respectively. These effects significantly increased Na current (INa) in hiPSC-CMs and in cardiomyocytes from mice trans-expressing Tbx5. Consequently, action potential (AP) amplitudes increased and QRS interval narrowed in the mouse electrocardiogram. p.F206L Tbx5 bound to the SCN5A promoter failed to transactivate it, thus precluding the pro-transcriptional effect of WT Tbx5. Therefore, p.F206L markedly decreased INa in hiPSC-CM, HL-1 cells, and mouse cardiomyocytes. The INa decrease in p.F206L trans-expressing mice translated into QRS widening and increased flecainide sensitivity. p.D111Y Tbx5 increased SCN5A expression but failed to repress CAMK2D and SPTBN4. The increased CaMKIIδ and βIV-spectrin significantly augmented the late component of INa (INaL) which, in turn, significantly prolonged AP duration in both hiPSC-CMs and mouse cardiomyocytes. Ranolazine, a selective INaL inhibitor, eliminated the QT and QTc intervals prolongation seen in p.D111Y trans-expressing mice.

CONCLUSIONS

In addition to peak INa, Tbx5 critically regulates INaL and the duration of repolarization in human cardiomyocytes. Our original results suggest that TBX5 variants associate with and modulate the intensity of the electrical phenotype in LQTS and BrS patients.

Cholecystokinin peptide signaling is regulated by a TBX5-MEF2 axis in the heart

The procholecystokinin (proCCK) gene encodes a secreted peptide known to regulate the digestive, endocrine, and nervous systems. Though recently proposed as a biomarker for heart dysfunction, its physiological role in both the embryonic and adult heart is poorly understood, and there are no reports of tissue-specific regulators of cholecystokinin signaling in the heart or other tissues. In the present study, mRNA of proCCK was observed in cardiac tissues during mouse embryonic development, establishing proCCK as an early marker of differentiated cardiomyocytes which is later restricted to anatomical subdomains of the neonatal heart. Three-dimensional analysis of the expression of proCCK and CCKAR/CCKBR receptors was performed using in situ hybridization and optical projection tomography, illustrating chamber-specific expression patterns in the postnatal heart. Transcription factor motif analyses indicated developmental cardiac transcription factors TBX5 and MEF2C as upstream regulators of proCCK, and this regulatory activity was confirmed in reporter gene assays. proCCK mRNA levels were also measured in the infarcted heart and in response to cyclic mechanical stretch and endothelin-1, indicating dynamic transcriptional regulation which might be leveraged for improved biomarker development. Functional analyses of exogenous cholecystokinin octapeptide (CCK-8) administration were performed in differentiating mouse embryonic stem cells (mESCs), and the results suggest that CCK-8 does not act as a differentiation modulator of cardiomyocyte subtypes. Collectively, these findings indicate that proCCK is regulated at the transcriptional level by TBX5-MEF2 and neurohormonal signaling, informing use of proCCK as a biomarker and future strategies for upstream manipulation of cholecystokinin signaling in the heart and other tissues.

Elevated EZH2 in ischemic heart disease epigenetically mediates suppression of NaV1.5 expression

Suppression of the cardiac sodium channel Na_V1.5 leads to fatal arrhythmias in ischemic heart disease (IHD). However, the transcriptional regulation of Na_V1.5 in cardiac ischemia is still unclear. Our studies are aimed to investigate the expression of enhancer of zeste homolog 2 (EZH2) in IHD and regulation of cardiac Na_V1.5 expression by EZH2. Human heart tissue was obtained from IHD and non-failing heart (NFH) patients; mouse heart tissue was obtained from the peri-infarct zone of hearts with myocardial infarction (MI) and hearts with a sham procedure. Protein and mRNA expression were measured by immunoblotting, immunostaining, and qRT-PCR. Protein-DNA binding and promoter activity were analyzed by ChIP-qPCR and luciferase assays, respectively. Na⁺ channel activity was assessed by whole-cell patch clamp recordings. EZH2 and H3K27me3 were increased while Na_V1.5 expression was reduced in IHD hearts and in mouse MI hearts compared to the controls. Reduced Na_V1.5 and increased EZH2 mRNA levels were observed in mouse MI hearts. A selective EZH2 inhibitor, GSK126 decreased H3K27me3 and elevated Na_V1.5 in HL-1 cells. Silencing of EZH2 expression decreased H3K27me3 and increased Na_V1.5 in these cells. EZH2 and H3K27me3 were enriched in the promoter regions of Scn5a and were decreased by treatment with EZH2 siRNA. GSK126 inhibited the enrichment of H3K27me3 in the Scn5a promoter and enhanced Scn5a transcriptional activity. GSK126 significantly increased Na⁺ channel activity. Taken together, EZH2 is increased in ischemic hearts and epigenetically suppresses Scn5a transcription by H3K27me3, leading to decreased Na_V1.5 expression and Na⁺ channel activity underlying the pathogenesis of arrhythmias.

TBX5-encoded T-box transcription factor 5 variant T223M is associated with long QT syndrome and pediatric sudden cardiac death

Long QT syndrome (LQTS) is a genetic disease resulting in a prolonged QT interval on a resting electrocardiogram, predisposing affected individuals to polymorphic ventricular tachycardia and sudden death. Although a number of genes have been implicated in this disease, nearly one in four individuals exhibiting the LQTS phenotype are genotype-negative. Whole-exome sequencing identified a missense T223M variant in TBX5 that cosegregates with prolonged QT interval in a family with otherwise genotype-negative LQTS and sudden death. The TBX5-T223M variant was absent among large ostensibly healthy populations (gnomAD) and predicted to be pathogenic by in silico modeling based on Panther, PolyPhen-2, Provean, SIFT, SNAP2, and PredictSNP prediction tools. The variant was located in a highly conserved region of TBX5 predicted to be part of the DNA-binding interface. A luciferase assay identified a 57.5% reduction in the ability of TBX5-T223M to drive expression at the atrial natriuretic factor promotor compared to wildtype TBX5 in vitro. We conclude that the variant is pathogenic in this family, and we put TBX5 forward as a disease susceptibility allele for genotype-negative LQTS. The identification of this familial variant may serve as a basis for the identification of previously unknown mechanisms of LQTS with broader implications for cardiac electrophysiology.

Life Cycle of the Cardiac Voltage-Gated Sodium Channel NaV1.5

Cardiac voltage-gated sodium channel Na_V1.5, encoded by SCN5A, is crucial for the upstroke of action potential and excitation of cardiomyocytes. Na_V1.5 undergoes complex processes before it reaches the target membrane microdomains and performs normal functions. A variety of protein partners are needed to achieve the balance between SCN5A transcription and mRNA decay, endoplasmic reticulum retention and export, Golgi apparatus retention and export, selective anchoring and degradation, activation, and inactivation of sodium currents. Subtle alterations can impair Na_V1.5 in terms of expression or function, eventually leading to Na_V1.5-associated diseases such as lethal arrhythmias and cardiomyopathy.

Electrical stimulation applied during differentiation drives the hiPSC-CMs towards a mature cardiac conduction-like cells

Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) resemble fetal cardiomyocytes and electrical stimulation (ES) has been explored to mature the differentiated cells. Here, we hypothesize that ES applied at the beginning of the differentiation process, triggers both differentiation of the hiPSC-CMs into a specialized conduction system (CS) phenotype and cell maturation. We applied ES for 15 days starting on day 0 of the differentiation process and found an increased expression of transcription factors and proteins associated with the development and function of CS including Irx3, Nkx2.5 and contactin 2, Hcn4 and Scn5a, respectively. We also found activation of intercalated disc proteins (Nrap and β-catenin). We detected ES-induced CM maturation as indicated by increased Tnni1 and Tnni3 expression. Confocal micrographs showed a shift towards expression of the gap junction protein connexin 40 in ES hiPSC-CM compared to the more dominant expression of connexin 43 in controls. Finally, analysis of functional parameters revealed that ES hiPSC-CMs exhibited faster action potential (AP) depolarization, longer intracellular Ca2+ transients, and slower AP duration at 90% of repolarization, resembling fast conducting fibers. Altogether, we provided evidence that ES during the differentiation of hiPSC to cardiomyocytes lead to development of cardiac conduction-like cells with more mature cytoarchitecture. Thus, hiPSC-CMs exposed to ES during differentiation can be instrumental to develop CS cells for cardiac disease modelling, screening individual drugs on a precison medicine type platform and support the development of novel therapeutics for arrhythmias.

Epigenetic Analyses of Human Left Atrial Tissue Identifies Gene Networks Underlying Atrial Fibrillation

BACKGROUND

Atrial fibrillation (AF) often arises from structural abnormalities in the left atria (LA). Annotation of the noncoding genome in human LA is limited, as are effects on gene expression and chromatin architecture. Many AF-associated genetic variants reside in noncoding regions; this knowledge gap impairs efforts to understand the molecular mechanisms of AF and cardiac conduction phenotypes.

METHODS

We generated a model of the LA noncoding genome by profiling 7 histone post-translational modifications (active: H3K4me3, H3K4me2, H3K4me1, H3K27ac, H3K36me3; repressive: H3K27me3, H3K9me3), CTCF binding, and gene expression in samples from 5 individuals without structural heart disease or AF. We used MACS2 to identify peak regions (P<0.01), applied a Markov model to classify regulatory elements, and annotated this model with matched gene expression data. We intersected chromatin states with expression quantitative trait locus, DNA methylation, and HiC chromatin interaction data from LA and left ventricle. Finally, we integrated genome-wide association data for AF and electrocardiographic traits to link disease-related variants to genes.

RESULTS

Our model identified 21 epigenetic states, encompassing regulatory motifs, such as promoters, enhancers, and repressed regions. Genes were regulated by proximal chromatin states; repressive states were associated with a significant reduction in gene expression (P<2×10-16). Chromatin states were differentially methylated, promoters were less methylated than repressed regions (P<2×10-16). We identified over 15 000 LA-specific enhancers, defined by homeobox family motifs, and annotated several cardiovascular disease susceptibility loci. Intersecting AF and PR genome-wide association studies loci with long-range chromatin conformation data identified a gene interaction network dominated by NKX2-5, TBX3, ZFHX3, and SYNPO2L.

CONCLUSIONS

Profiling the noncoding genome provides new insights into the gene expression and chromatin regulation in human LA tissue. These findings enabled identification of a gene network underlying AF; our experimental and analytic approach can be extended to identify molecular mechanisms for other cardiac diseases and traits.

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.

Role of Non-Coding Variants in Brugada Syndrome

Brugada syndrome (BrS) is an inherited electrical heart disease associated with a high risk of sudden cardiac death (SCD). The genetic characterization of BrS has always been challenging. Although several cardiac ion channel genes have been associated with BrS, SCN5A is the only gene that presents definitive evidence for causality to be used for clinical diagnosis of BrS. However, more than 65% of diagnosed cases cannot be explained by variants in SCN5A or other genes. Therefore, in an important number of BrS cases, the underlying mechanisms are still elusive. Common variants, mostly located in non-coding regions, have emerged as potential modulators of the disease by affecting different regulatory mechanisms, including transcription factors (TFs), three-dimensional organization of the genome, or non-coding RNAs (ncRNAs). These common variants have been hypothesized to modulate the interindividual susceptibility of the disease, which could explain incomplete penetrance of BrS observed within families. Altogether, the study of both common and rare variants in parallel is becoming increasingly important to better understand the genetic basis underlying BrS. In this review, we aim to describe the challenges of studying non-coding variants associated with disease, re-examine the studies that have linked non-coding variants with BrS, and provide further evidence for the relevance of regulatory elements in understanding this cardiac disorder.

Biphasic Electrical Pulse by a Micropillar Electrode Array Enhances Maturation and Drug Response of Reprogrammed Cardiac Spheroids

Direct reprogramming is an efficient strategy to produce cardiac lineage cells necessary for cardiac tissue engineering and drug testing for cardiac toxicity. However, functional maturation of reprogrammed cardiomyocytes, which is of great importance for their regenerative potential and drug response, still remains challenging. In this study, we propose a novel electrode platform to promote direct cardiac reprogramming and improve the functionality of reprogrammed cardiac cells. Nonviral cardiac reprogramming was improved via a three-dimensional spheroid culture of chemically induced cardiomyocytes exposed to a small-molecule cocktail. A micropillar electrode array providing biphasic electrical pulses mimicking the heartbeat further enhanced maturation and electrophysiological properties of reprogrammed cardiac spheroids, leading to proper responses and increased sensitivity to drugs. On the basis of our results, we conclude that our device may have a wider application in the generation of functional cardiac cells for regenerative medicine and screening of novel drugs.

T-box transcription factor 3 governs a transcriptional program for the function of the mouse atrioventricular conduction system

Genome-wide association studies have identified noncoding variants near TBX3 that are associated with PR interval and QRS duration, suggesting that subtle changes in TBX3 expression affect atrioventricular conduction system function. To explore whether and to what extent the atrioventricular conduction system is affected by Tbx3 dose reduction, we first characterized electrophysiological properties and morphology of heterozygous Tbx3 mutant (Tbx3+/-) mouse hearts. We found PR interval shortening and prolonged QRS duration, as well as atrioventricular bundle hypoplasia after birth in heterozygous mice. The atrioventricular node size was unaffected. Transcriptomic analysis of atrioventricular nodes isolated by laser capture microdissection revealed hundreds of deregulated genes in Tbx3+/- mutants. Notably, Tbx3+/- atrioventricular nodes showed increased expression of working myocardial gene programs (mitochondrial and metabolic processes, muscle contractility) and reduced expression of pacemaker gene programs (neuronal, Wnt signaling, calcium/ion channel activity). By integrating chromatin accessibility profiles (ATAC sequencing) of atrioventricular tissue and other epigenetic data, we identified Tbx3-dependent atrioventricular regulatory DNA elements (REs) on a genome-wide scale. We used transgenic reporter assays to determine the functionality of candidate REs near Ryr2, an up-regulated chamber-enriched gene, and in Cacna1g, a down-regulated conduction system-specific gene. Using genome editing to delete candidate REs, we showed that a strong intronic bipartite RE selectively governs Cacna1g expression in the conduction system in vivo. Our data provide insights into the multifactorial Tbx3-dependent transcriptional network that regulates the structure and function of the cardiac conduction system, which may underlie the differences in PR duration and QRS interval between individuals carrying variants in the TBX3 locus.

miR-182-5p is an evolutionarily conserved Tbx5 effector that impacts cardiac development and electrical activity in zebrafish

To dissect the TBX5 regulatory circuit, we focused on microRNAs (miRNAs) that collectively contribute to make TBX5 a pivotal cardiac regulator. We profiled miRNAs in hearts isolated from wild-type, CRE, Tbx5lox/+and Tbx5del/+ mice using a Next Generation Sequencing (NGS) approach. TBX5 deficiency in cardiomyocytes increased the expression of the miR-183 cluster family that is controlled by Kruppel-like factor 4, a transcription factor repressed by TBX5. MiR-182-5p, the most highly expressed miRNA of this family, was functionally analyzed in zebrafish. Transient overexpression of miR-182-5p affected heart morphology, calcium handling and the onset of arrhythmias as detected by ECG tracings. Accordingly, several calcium channel proteins identified as putative miR-182-5p targets were downregulated in miR-182-5p overexpressing hearts. In stable zebrafish transgenic lines, we demonstrated that selective miRNA-182-5p upregulation contributes to arrhythmias. Moreover, cardiac-specific down-regulation of miR-182-5p rescued cardiac defects in a zebrafish model of Holt-Oram syndrome. In conclusion, miR-182-5p exerts an evolutionarily conserved role as a TBX5 effector in the onset of cardiac propensity for arrhythmia, and constitutes a relevant target for mediating the relationship between TBX5, arrhythmia and heart development.

Identifying the key regulators that promote cell-cycle activity in the hearts of early neonatal pigs after myocardial injury

Mammalian cardiomyocytes exit the cell cycle shortly after birth. As a result, an occurrence of coronary occlusion-induced myocardial infarction often results in heart failure, postinfarction LV dilatation, or death, and represents one of the most significant public health morbidities worldwide. Interestingly however, the hearts of neonatal pigs have been shown to regenerate following an acute myocardial infarction (MI) occuring on postnatal day 1 (P1); a recovery period which is accompanied by an increased expression of markers for cell-cycle activity, and suggests that early postnatal myocardial regeneration may be driven in part by the MI-induced proliferation of pre-existing cardiomyocytes. In this study, we identified signaling pathways known to regulate the cell cycle, and determined of these, the pathways persistently upregulated in response to MI injury. We identified five pathways (mitogen associated protein kinase [MAPK], Hippo, cyclic [cAMP], Janus kinase/signal transducers and activators of transcription [JAK-STAT], and Ras) which were comprehensively upregulated in cardiac tissues collected on day 7 (P7) and/or P28 of the P1 injury hearts. Several of the initiating master regulators (e.g., CSF1/CSF1R, TGFB, and NPPA) and terminal effector molecules (e.g., ATF4, FOS, RELA/B, ITGB2, CCND1/2/3, PIM1, RAF1, MTOR, NKF1B) in these pathways were persistently upregulated at day 7 through day 28, suggesting there exists at least some degree of regenerative activity up to 4 weeks following MI at P1. Our observations provide a list of key regulators to be examined in future studies targeting cell-cycle activity as an avenue for myocardial regeneration.

Transcriptional Patterning of the Ventricular Cardiac Conduction System

RATIONALE

The heartbeat is organized by the cardiac conduction system (CCS), a specialized network of cardiomyocytes. Patterning of the CCS into atrial node versus ventricular conduction system (VCS) components with distinct physiology is essential for the normal heartbeat. Distinct node versus VCS physiology has been recognized for more than a century, but the molecular basis of this regional patterning is not well understood.

OBJECTIVE

To study the genetic and genomic mechanisms underlying node versus VCS distinction and investigate rhythm consequences of failed VCS patterning.

METHODS AND RESULTS

Using mouse genetics, we found that the balance between T-box transcriptional activator, Tbx5, and T-box transcriptional repressor, Tbx3, determined the molecular and functional output of VCS myocytes. Adult VCS-specific removal of Tbx5 or overexpression of Tbx3 re-patterned the fast VCS into slow, nodal-like cells based on molecular and functional criteria. In these cases, gene expression profiling showed diminished expression of genes required for VCS-specific fast conduction but maintenance of expression of genes required for nodal slow conduction physiology. Action potentials of Tbx5-deficient VCS myocytes adopted nodal-specific characteristics, including increased action potential duration and cellular automaticity. Removal of Tbx5 in vivo precipitated inappropriate depolarizations in the atrioventricular (His)-bundle associated with lethal ventricular arrhythmias. TBX5 bound and directly activated cis-regulatory elements at fast conduction channel genes required for fast physiological characteristics of the VCS action potential, defining the identity of the adult VCS.

CONCLUSIONS

The CCS is patterned entirely as a slow, nodal ground state, with a T-box dependent, physiologically dominant, fast conduction network driven specifically in the VCS. Disruption of the fast VCS gene regulatory network allowed nodal physiology to emerge, providing a plausible molecular mechanism for some lethal ventricular arrhythmias.

Trait-associated noncoding variant regions affect TBX3 regulation and cardiac conduction

Genome-wide association studies have implicated common genomic variants in the gene desert upstream of TBX3 in cardiac conduction velocity. Whether these noncoding variants affect expression of TBX3 or neighboring genes and how they affect cardiac conduction is not understood. Here, we use high-throughput STARR-seq to test the entire 1.3 Mb human and mouse TBX3 locus, including two cardiac conduction-associated variant regions, for regulatory function. We identified multiple accessible and functional regulatory DNA elements that harbor variants affecting their activity. Both variant regions drove gene expression in the cardiac conduction tissue in transgenic reporter mice. Genomic deletion from the mouse genome of one of the regions caused increased cardiac expression of only Tbx3, PR interval shortening and increased QRS duration. Combined, our findings address the mechanistic link between trait-associated variants in the gene desert, TBX3 regulation and cardiac conduction.

ZO-1 Regulates Intercalated Disc Composition and Atrioventricular Node Conduction

RATIONALE

ZO-1 (Zona occludens 1), encoded by the tight junction protein 1 (TJP1) gene, is a regulator of paracellular permeability in epithelia and endothelia. ZO-1 interacts with the actin cytoskeleton, gap, and adherens junction proteins and localizes to intercalated discs in cardiomyocytes. However, the contribution of ZO-1 to cardiac physiology remains poorly defined.

OBJECTIVE

We aim to determine the role of ZO-1 in cardiac function.

METHODS AND RESULTS

Inducible cardiomyocyte-specific Tjp1 deletion mice (Tjp1fl/fl; Myh6Cre/Esr1*) were generated by crossing the Tjp1 floxed mice and Myh6Cre/Esr1* transgenic mice. Tamoxifen-induced loss of ZO-1 led to atrioventricular (AV) block without changes in heart rate, as measured by ECG and ex vivo optical mapping. Mice with tamoxifen-induced conduction system-specific deletion of Tjp1 (Tjp1fl/fl; Hcn4CreERt2) developed AV block while tamoxifen-induced conduction system deletion of Tjp1 distal to the AV node (Tjp1fl/fl; Kcne1CreERt2) did not demonstrate conduction defects. Western blot and immunostaining analyses of AV nodes showed that ZO-1 loss decreased Cx (connexin) 40 expression and intercalated disc localization. Consistent with the mouse model study, immunohistochemical staining showed that ZO-1 is abundantly expressed in the human AV node and colocalizes with Cx40. Ventricular conduction was not altered despite decreased localization of ZO-1 and Cx43 at the ventricular intercalated disc and modestly decreased left ventricular ejection fraction, suggesting ZO-1 is differentially required for AV node and ventricular conduction.

CONCLUSIONS

ZO-1 is a key protein responsible for maintaining appropriate AV node conduction through maintaining gap junction protein localization.

Joint reconstruction of cis-regulatory interaction networks across multiple tissues using single-cell chromatin accessibility data

The rapid accumulation of single-cell chromatin accessibility data offers a unique opportunity to investigate common and specific regulatory mechanisms across different cell types. However, existing methods for cis-regulatory network reconstruction using single-cell chromatin accessibility data were only designed for cells belonging to one cell type, and resulting networks may be incomparable directly due to diverse cell numbers of different cell types. Here, we adopt a computational method to jointly reconstruct cis-regulatory interaction maps (JRIM) of multiple cell populations based on patterns of co-accessibility in single-cell data. We applied JRIM to explore common and specific regulatory interactions across multiple tissues from single-cell ATAC-seq dataset containing ~80 000 cells across 13 mouse tissues. Reconstructed common interactions among 13 tissues indeed relate to basic biological functions, and individual cis-regulatory networks show strong tissue specificity and functional relevance. More importantly, tissue-specific regulatory interactions are mediated by coordination of histone modifications and tissue-related TFs, and many of them may reveal novel regulatory mechanisms.

Atrial fibrillation risk loci interact to modulate Ca2+-dependent atrial rhythm homeostasis

Atrial fibrillation (AF), defined by disorganized atrial cardiac rhythm, is the most prevalent cardiac arrhythmia worldwide. Recent genetic studies have highlighted a major heritable component and identified numerous loci associated with AF risk, including the cardiogenic transcription factor genes TBX5, GATA4, and NKX2-5. We report that Tbx5 and Gata4 interact with opposite signs for atrial rhythm controls compared with cardiac development. Using mouse genetics, we found that AF pathophysiology caused by Tbx5 haploinsufficiency, including atrial arrhythmia susceptibility, prolonged action potential duration, and ectopic cardiomyocyte depolarizations, were all rescued by Gata4 haploinsufficiency. In contrast, Nkx2-5 haploinsufficiency showed no combinatorial effect. The molecular basis of the TBX5/GATA4 interaction included normalization of intra-cardiomyocyte calcium flux and expression of calcium channel genes Atp2a2 and Ryr2. Furthermore, GATA4 and TBX5 showed antagonistic interactions on an Ryr2 enhancer. Atrial rhythm instability caused by Tbx5 haploinsufficiency was rescued by a decreased dose of phospholamban, a sarco/endoplasmic reticulum Ca2+-ATPase inhibitor, consistent with a role for decreased sarcoplasmic reticulum calcium flux in Tbx5-dependent AF susceptibility. This work defines a link between Tbx5 dose, sarcoplasmic reticulum calcium flux, and AF propensity. The unexpected interactions between Tbx5 and Gata4 in atrial rhythm control suggest that evaluating specific interactions between genetic risk loci will be necessary for ascertaining personalized risk from genetic association data.

Direct and Indirect Suppression of Scn5a Gene Expression Mediates Cardiac Na+ Channel Inhibition by Wnt Signalling

BACKGROUND

Myocardial infarction and heart failure are associated with reduced voltage-gated Na⁺ current (I_Na) that promotes arrhythmias and sudden deaths. We have previously shown that the Wnt/β-catenin signalling (Wnt signalling), which is active in heart disease, reduces cardiac I_Na, suggesting that Wnt signalling may be a potential therapeutic target. However, because Wnt signalling is required for the homeostasis of many noncardiac tissues, administration of Wnt inhibitors to heart patients would cause significant side effects. The present study aims to elucidate the molecular mechanisms of cardiac I_Na inhibition by Wnt, which would identify cardiac-specific therapeutic targets.

METHODS

Wnt signalling was activated in neonatal rat ventricular myocytes by Wnt3a protein. Adenovirus expressing Wnt3a was injected into the adult rat ventricle. CRISPR/Cas9 and chromatin immunoprecipitation were used for mechanistic studies.

RESULTS

Wnt signalling activation in neonatal rat ventricular myocytes reduced Na_v1.5 protein and Scn5a mRNA, but increased Tbx3, a known suppressor of Scn5a. Chromatin immunoprecipitation showed that Wnt signalling inhibits Scn5a expression through downstream mediator (TCF4) binding to both Tbx3 and Scn5a promoters. Overexpression or knockdown of Tbx3 directly modified Na_v1.5 and I_Na, whereas CRISPR/Cas9-induced mutations at TCF4 binding sites within the Scn5a promoter attenuated Wnt inhibition of Scn5a and Na_v1.5. In adult rat hearts, adenovirus expressing Wnt3a reduced Na_v1.5, increased QRS duration in electrocardiogram, and increased the susceptibility to ventricular tachycardia.

CONCLUSIONS

Wnt signalling inhibits the Na⁺ channel by direct and indirect (via Tbx3) suppression of Scn5a transcription. Strategies to block TCF4 binding to the Tbx3 and Scn5a promoters would represent novel strategies for cardiac-specific inhibition of the Wnt pathway to rescue I_Na and prevent sudden cardiac deaths.

Genetic risk and atrial fibrillation in patients with heart failure

AIMS

To study the association between an atrial fibrillation (AF) genetic risk score with prevalent AF and all-cause mortality in patients with heart failure.

METHODS AND RESULTS

An AF genetic risk score was calculated in 3759 European ancestry individuals (1783 with sinus rhythm, 1976 with AF) from the BIOlogy Study to TAilored Treatment in Chronic Heart Failure (BIOSTAT-CHF) by summing 97 single nucleotide polymorphism (SNP) alleles (ranging from 0-2) weighted by the natural logarithm of the relative SNP risk from the latest AF genome-wide association study. Further, we assessed AF risk variance explained by additive SNP variation, and performance of clinical or genetic risk factors, and the combination in classifying AF prevalence. AF was classified as AF or atrial flutter (AFL) at baseline electrocardiogram and/or a history of AF or AFL. The genetic risk score was associated with AF after multivariable adjustment. Odds ratio for AF prevalence per 1-unit increase genetic risk score was 2.12 (95% confidence interval 1.84-2.45, P = 2.15 × 10^-24 ) in the total cohort, 2.08 (1.72-2.50, P = 1.30 × 10^-14 ) in heart failure with reduced ejection fraction (HFrEF) and 2.02 (1.37-2.99, P = 4.37 × 10^-4 ) in heart failure with preserved ejection fraction (HFpEF). AF-associated loci explained 22.9% of overall AF SNP heritability. Addition of the genetic risk score to clinical risk factors increased the C-index by 2.2% to 0.721.

CONCLUSIONS

The AF genetic risk score was associated with increased AF prevalence in HFrEF and HFpEF. Genetic variation accounted for 22.9% of overall AF SNP heritability. Addition of genetic risk to clinical risk improved model performance in classifying AF prevalence.

β-Adrenergic stimuli and rotating suspension culture enhance conversion of human adipogenic mesenchymal stem cells into highly conductive cardiac progenitors

Clinical trials using human adipogenic mesenchymal stem cells (hAdMSCs) for the treatment of cardiac diseases have shown improvement in cardiac function and were proven safe. However, hAdMSCs do not convert efficiently into cardiomyocytes (CMs) or vasculature. Thus, reprogramming hAdMSCs into myocyte progenitors may fare better in future investigations. To reprogramme hAdMSCs into electrically conductive cardiac progenitor cells, we pioneered a three-step reprogramming strategy that uses proven MESP1/ETS2 transcription factors, β-adrenergic and hypoxic signalling induced in three-dimensional (3D) cardiospheres. In Stage 1, ETS2 and MESP1 activated NNKX2.5, TBX5, MEF2C, dHAND, and GATA4 during the conversion of hAdMSCs into cardiac progenitor cells. Next, in Stage 2, β2AR activation repositioned cardiac progenitors into de novo immature conductive cardiac cells, along with the appearance of RYR2, CAV2.1, CAV3.1, NAV1.5, SERCA2, and CX45 gene transcripts and displayed action potentials. In Stage 3, electrical conduction that was fostered by 3D cardiospheres formed in a Synthecon®, Inc. rotating bioreactor induced the appearance of hypoxic genes: HIF-1α/β, PCG 1α/β, and NOS2, which coincided with the robust activation of adult contractile genes including MLC2v, TNNT2, and TNNI3, ion channel genes, and the appearance of hyperpolarization-activated and cyclic nucleotide-gated channels (HCN1-4). Conduction velocities doubled to ~200 mm/s after hypoxia and doubled yet again after dissociation of the 3D cell clusters to ~400 mm/s. By comparison, normal conduction velocities within working ventricular myocytes in the whole heart range from 0.5 to 1 m/s. Epinephrine stimulation of stage 3 cardiac cells in patches resulted in an increase in amplitude of the electrical wave, indicative of conductive cardiac cells. Our efficient protocol that converted hAdMSCs into highly conductive cardiac progenitors demonstrated the potential utilization of stage 3 cells for tissue engineering applications for cardiac repair.

An enhancer cluster controls gene activity and topology of the SCN5A-SCN10A locus in vivo

Mutations and variations in and around SCN5A, encoding the major cardiac sodium channel, influence impulse conduction and are associated with a broad spectrum of arrhythmia disorders. Here, we identify an evolutionary conserved regulatory cluster with super enhancer characteristics downstream of SCN5A, which drives localized cardiac expression and contains conduction velocity-associated variants. We use genome editing to create a series of deletions in the mouse genome and show that the enhancer cluster controls the conformation of a >0.5 Mb genomic region harboring multiple interacting gene promoters and enhancers. We find that this cluster and its individual components are selectively required for cardiac Scn5a expression, normal cardiac conduction and normal embryonic development. Our studies reveal physiological roles of an enhancer cluster in the SCN5A-SCN10A locus, show that it controls the chromatin architecture of the locus and Scn5a expression, and suggest that genetic variants affecting its activity may influence cardiac function.

A small-molecule LF3 abrogates β-catenin/TCF4-mediated suppression of NaV1.5 expression in HL-1 cardiomyocytes

Increased nuclear β-catenin interacting with T-cell factor 4 (TCF4) affects the expression of target genes including SCN5A in ischemic heart disease, which is characterized by frequent ventricular tachycardia/fibrillation. A complex of β-catenin and TCF4 inhibits cardiac Na+ channel activity by reducing NaV1.5 expression through suppressing SCN5A promoter activity in HL-1 cardiomyocytes. LF3, a 4-thioureido-benzenesulfonamide derivative and an inhibitor of β-catenin/TCF4 interaction, has been shown to block the self-renewal capacity of cancer stem cells. We performed studies to determine if LF3 can reverse suppressive effects of β-catenin/TCF4 signaling on the expression of NaV1.5 in HL-1 cardiomyocytes. Western blotting and real-time qRT-PCR analyses showed that 10 μM LF3 significantly increased the expression of NaV1.5 but it did not alter β-catenin and TCF4 expression. Subcellular fractionation analysis demonstrated that LF3 significantly increased the levels of NaV1.5 in both membrane and cytoplasm. Whole-cell patch-clamp recordings revealed that Na+ currents were significantly increased with no changes in the steady-state parameters, activation and inactivation time constants and recovery from inactivation of Na+ channel in HL-1 cells treated with LF3. Immunoprecipitation exhibited that LF3 blocked the interaction of β-catenin and TCF4. Luciferase reporter assays performed in HEK 293 cells and HL-1 revealed that LF3 increased the SCN5A promoter activity in HL-1 cells and prevented β-catenin suppressive effect on SCN5A promoter activity in HEK 293 cells. Taken together, we conclude that LF3, an inhibitor of β-catenin/TCF4 interaction, elevates NaV1.5 expression, leading to increase Na+ channel activity in HL-1 cardiomyocytes.