Phil Barnett

Background and Aims: Life-threatening arrhythmias are a well-established consequence of reduced cardiac sodium current (INa). Gene therapy approaches to increase INa have demonstrated potential benefits to prevent arrhythmias. However, the development of such therapies is hampered by the large size of sodium channels. In this study, SCN10A-short (S10s), a short transcript encoding the carboxy-terminal domain of the human neuronal sodium channel, was evaluated as a gene therapy target to increase INa and prevent arrhythmias. Methods: Adeno-associated viral vector overexpressing S10s was injected into wild type and Scn5a-haploinsufficient mice on which patch-clamp studies, optical mapping, electrocardiogram analyses, and ischaemia reperfusion were performed. In vitro and in silico studies were conducted to further explore the effect of S10s gene therapy in the context of human hearts. Results: Cardiac S10s overexpression increased cellular INa, maximal action potential upstroke velocity, and action potential amplitude in Scn5a-haploinsufficient cardiomyocytes. S10s gene therapy rescues conduction slowing in Scn5a-haploinsufficient mice and prevented ventricular tachycardia induced by ischaemia-reperfusion in wild type mice. S10s overexpression increased maximal action potential upstroke velocity in human inducible pluripotent stem cell-derived cardiomyocytes and prevented inducible arrhythmias in simulated human heart models. Conclusions: S10s gene therapy may be effective to treat cardiac conduction abnormalities and associated arrhythmias.

Genetic predisposition to cardiac arrhythmias has been a field of intense investigation. Research initially focused on rare hereditary arrhythmias, but over the last two decades, the role of genetic variation (single nucleotide polymorphisms) in heart rate, rhythm, and arrhythmias has been taken into consideration as well. In particular, genome-wide association studies have identified hundreds of genomic loci associated with quantitative electrocardiographic traits, atrial fibrillation, and less common arrhythmias such as Brugada syndrome. A significant number of associated variants have been found to systematically localize in non-coding regulatory elements that control the tissue-specific and temporal transcription of genes encoding transcription factors, ion channels, and other proteins. However, the identification of causal variants and the mechanism underlying their impact on phenotype has proven difficult due to the complex tissue-specific, time-resolved, condition-dependent, and combinatorial function of regulatory elements, as well as their modest conservation across different model species. In this review, we discuss research efforts aimed at identifying and characterizing-trait-associated variant regulatory elements and the molecular mechanisms underlying their impact on heart rate or rhythm.

Background: Congenital heart diseases (CHD) are the most common congenital malformations in newborns and remain the leading cause of mortality among infants under one year old. Molecular diagnosis is crucial to evaluate the recurrence risk and to address future prenatal diagnosis. Here, we describe two families with various forms of inherited non-syndromic CHD and the genetic work-up and resultant findings. Methods: Next-generation sequencing (NGS) was employed in both families to uncover the genetic cause. In addition, we performed functional analysis to investigate the consequences of the identified variants in vitro. Results: NGS identified possible causative variants in both families in the protein kinase domain of the TGFBR1 gene. These variants occurred on the same amino acid, but resulted in differently substituted amino acids (p.R398C/p.R398H). Both variants co-segregate with the disease, are extremely rare or unique, and occur in an evolutionary highly conserved domain of the protein. Furthermore, both variants demonstrated a significantly altered TGFBR1-smad signaling activity. Clinical investigation revealed that none of the carriers had (signs of) aortopathy. Conclusion: In conclusion, we describe two families, with various forms of inherited non-syndromic CHD without aortopathies, associated with unique/rare variants in TGFBR1 that display altered TGF-beta signaling. These findings highlight involvement of TGFBR1 in CHD, and warrant consideration of potential causative TGFBR1 variants also in CHD patients without aortopathies.

Rationale: Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored. Objective: We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA. Methods and Results: We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10 ^-10 , OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10 ^-5 ). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A , which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart. Conclusions: This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A . Genomic and functional data support a causal role of WNT5A at the locus.

Myocardial infarction causes ventricular muscle loss and formation of scar tissue. The surviving myocardium in the border zone, located adjacent to the infarct, undergoes profound changes in function, structure and composition. How and to what extent these changes of border zone cardiomyocytes are regulated epigenetically is not fully understood. Here, we obtained transcriptomes of PCM-1-sorted mouse cardiomyocyte nuclei of healthy left ventricle and 7 days post myocardial infarction border zone tissue. We validated previously observed downregulation of genes involved in fatty acid metabolism, oxidative phosphorylation and mitochondrial function in border zone-derived cardiomyocytes, and observed a modest induction of genes involved in glycolysis, including Slc2a1 (Glut1) and Pfkp. To gain insight into the underlying epigenetic regulatory mechanisms, we performed H3K27ac profiling of healthy and border zone cardiomyocyte nuclei. We confirmed the switch from Mef2- to AP-1 chromatin association in border zone cardiomyocytes, and observed, in addition, an enrichment of PPAR/RXR binding motifs in the sites with reduced H3K27ac signal. We detected downregulation and accompanying epigenetic state changes at several key PPAR target genes including Ppargc1a (PGC-1α), Cpt2, Ech1, Fabpc3 and Vldrl in border zone cardiomyocytes. These data indicate that changes in epigenetic state and gene regulation underlie the maintained metabolic switch in border zone cardiomyocytes.

Atrial natriuretic factor and brain natriuretic peptide are two important biomarkers in clinical cardiology. These two natriuretic peptide hormones are encoded by the paralogous genes Nppa and Nppb, which are evolutionary conserved. Both genes are predominantly expressed by the heart muscle during the embryonic and fetal stages, and in particular Nppa expression is strongly reduced in the ventricles after birth. Upon cardiac stress, Nppa and Nppb are strongly upregulated in the ventricular myocardium. Much is known about the molecular and physiological ques inducing Nppa and Nppb expression; however, the transcriptional regulatory mechanisms of the Nppa–Nppb cluster in vivo has proven to be quite complex and is not well understood. In this review, we will provide recent insights into the dynamic and complex regulation of Nppa and Nppb during heart development and hypertrophy, and the association of this gene cluster with the cardiomyocyte-intrinsic program of heart regeneration.

ChIP-seq has become a major tool for the genome-wide identification of transcription factor binding or histone modification sites. Most peak-calling algorithms require input control datasets to model the occurrence of background reads to account for local sequencing and GC bias. However, the GC-content of reads in Input-seq datasets deviates significantly from that in ChIP-seq datasets. Moreover, we observed that a commonly used peak calling program performed equally well when the use of a simulated uniform background set was compared to an Input-seq dataset. This contradicts the assumption that input control datasets are necessary to fatefully reflect the background read distribution. Because the GC-content of the abundant single reads in ChIP-seq datasets is similar to those of randomly sampled regions we designed a peak-calling algorithm with a background model based on overlapping single reads. The application, OccuPeak, uses the abundant low frequency tags present in each ChIP-seq dataset to model the background, thereby avoiding the need for additional datasets. Analysis of the performance of OccuPeak showed robust model parameters. Its measure of peak significance, the excess ratio, is only dependent on the tag density of a peak and the global noise levels. Compared to the commonly used peak-calling applications MACS and CisGenome, OccuPeak had the highest sensitivity in an enhancer identification benchmark test, and performed similar in an overlap tests of transcription factor occupation with DNase I hypersensitive sites and H3K27ac sites. Moreover, peaks called by OccuPeak were significantly enriched with cardiac disease-associated SNPs. OccuPeak runs as a standalone application and does not require extensive tweaking of parameters, making its use straightforward and user friendly. Availability: http://occupeak.hfrc.nl

Variants in SCN10A, which encodes a voltage-gated sodium channel, are associated with alterations of cardiac conduction parameters and the cardiac rhythm disorder Brugada syndrome; however, it is unclear how SCN10A variants promote dysfunctional cardiac conduction. Here we showed by high-resolution 4C-seq analysis of the Scn10a-Scn5a locus in murine heart tissue that a cardiac enhancer located in Scn10a, encompassing SCN10A functional variant rs6801957, interacts with the promoter of Scn5a, a sodium channel–encoding gene that is critical for cardiac conduction. We observed that SCN5A transcript levels were several orders of magnitude higher than SCN10A transcript levels in both adult human and mouse heart tissue. Analysis of BAC transgenic mouse strains harboring an engineered deletion of the enhancer within Scn10a revealed that the enhancer was essential for Scn5a expression in cardiac tissue. Furthermore, the common SCN10A variant rs6801957 modulated Scn5a expression in the heart. In humans, the SCN10A variant rs6801957, which correlated with slowed conduction, was associated with reduced SCN5A expression. These observations establish a genomic mechanism for how a common genetic variation at SCN10A influences cardiac physiology and predisposes to arrhythmia.