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

SCN10A-short gene therapy to restore conduction and protect against malignant cardiac arrhythmias

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.

A Rare Noncoding Enhancer Variant in SCN5A Contributes to the High Prevalence of Brugada Syndrome in Thailand

BACKGROUND

Brugada syndrome (BrS) is a cardiac arrhythmia disorder that causes sudden death in young adults. Rare genetic variants in the SCN5A gene encoding the Nav1.5 sodium channel and common noncoding variants at this locus are robustly associated with the condition. BrS is particularly prevalent in Southeast Asia but the underlying ancestry-specific factors remain largely unknown.

METHODS

Genome sequencing of BrS probands and population-matched controls from Thailand was performed to identify rare noncoding variants at the SCN5A-SCN10A locus that were enriched in patients with BrS. A likely causal variant was prioritized by computational methods and introduced into human induced pluripotent stem cell (hiPSC) lines using CRISPR-Cas9. The effect of the variant on SCN5A expression and Nav1.5 sodium channel current was then assessed in hiPSC-derived cardiomyocytes (hiPSC-CMs).

RESULTS

A rare noncoding variant in an SCN5A intronic enhancer region was highly enriched in patients with BrS (detected in 3.9% of cases with a case-control odds ratio of 45.2). The variant affects a nucleotide conserved across all mammalian species and predicted to disrupt a Mef2 transcription factor binding site. Heterozygous introduction of the enhancer variant in hiPSC-CMs caused significantly reduced SCN5A expression from the variant-containing allele and a 30% reduction in Nav1.5-mediated sodium current density compared with isogenic controls, confirming its pathogenicity. Patients with the variant had severe phenotypes, with 89% experiencing cardiac arrest.

CONCLUSIONS

This is the first example of a functionally validated rare noncoding variant at the SCN5A locus and highlights how genome sequencing in understudied populations can identify novel disease mechanisms. The variant partly explains the increased prevalence of BrS in this region and enables the identification of at-risk variant carriers to reduce the burden of sudden cardiac death in Thailand.

Genetic and Epigenetic Mechanisms Regulating Blood Pressure and Kidney Dysfunction

The pioneering work of Dr Lewis K. Dahl established a relationship between kidney, salt, and high blood pressure (BP), which led to the major genetic-based experimental model of hypertension. BP, a heritable quantitative trait affected by numerous biological and environmental stimuli, is a major cause of morbidity and mortality worldwide and is considered to be a primary modifiable factor in renal, cardiovascular, and cerebrovascular diseases. Genome-wide association studies have identified monogenic and polygenic variants affecting BP in humans. Single nucleotide polymorphisms identified in genome-wide association studies have quantified the heritability of BP and the effect of genetics on hypertensive phenotype. Changes in the transcriptional program of genes may represent consequential determinants of BP, so understanding the mechanisms of the disease process has become a priority in the field. At the molecular level, the onset of hypertension is associated with reprogramming of gene expression influenced by epigenomics. This review highlights the specific genetic variants, mutations, and epigenetic factors associated with high BP and how these mechanisms affect the regulation of hypertension and kidney dysfunction.

A mammalian tripartite enhancer cluster controls hypothalamic Pomc expression, food intake, and body weight

Food intake and energy balance are tightly regulated by a group of hypothalamic arcuate neurons expressing the proopiomelanocortin (POMC) gene. In mammals, arcuate-specific POMC expression is driven by two cis-acting transcriptional enhancers known as nPE1 and nPE2. Because mutant mice lacking these two enhancers still showed hypothalamic Pomc mRNA, we searched for additional elements contributing to arcuate Pomc expression. By combining molecular evolution with reporter gene expression in transgenic zebrafish and mice, here, we identified a mammalian arcuate-specific Pomc enhancer that we named nPE3, carrying several binding sites also present in nPE1 and nPE2 for transcription factors known to activate neuronal Pomc expression, such as ISL1, NKX2.1, and ERα. We found that nPE3 originated in the lineage leading to placental mammals and remained under purifying selection in all mammalian orders, although it was lost in Simiiformes (monkeys, apes, and humans) following a unique segmental deletion event. Interestingly, ablation of nPE3 from the mouse genome led to a drastic reduction (>70%) in hypothalamic Pomc mRNA during development and only moderate (<33%) in adult mice. Comparison between double (nPE1 and nPE2) and triple (nPE1, nPE2, and nPE3) enhancer mutants revealed the relative contribution of nPE3 to hypothalamic Pomc expression and its importance in the control of food intake and adiposity in male and female mice. Altogether, these results demonstrate that nPE3 integrates a tripartite cluster of partially redundant enhancers that originated upon a triple convergent evolutionary process in mammals and that is critical for hypothalamic Pomc expression and body weight homeostasis.

Role of Genetic Variation in Transcriptional Regulatory Elements in Heart Rhythm

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.

Dynamic changes in P300 enhancers and enhancer-promoter contacts control mouse cardiomyocyte maturation

Cardiomyocyte differentiation continues throughout murine gestation and into the postnatal period, driven by temporally regulated expression changes in the transcriptome. The mechanisms that regulate these developmental changes remain incompletely defined. Here, we used cardiomyocyte-specific ChIP-seq of the activate enhancer marker P300 to identify 54,920 cardiomyocyte enhancers at seven stages of murine heart development. These data were matched to cardiomyocyte gene expression profiles at the same stages and to Hi-C and H3K27ac HiChIP chromatin conformation data at fetal, neonatal, and adult stages. Regions with dynamic P300 occupancy exhibited developmentally regulated enhancer activity, as measured by massively parallel reporter assays in cardiomyocytes in vivo, and identified key transcription factor-binding motifs. These dynamic enhancers interacted with temporal changes of the 3D genome architecture to specify developmentally regulated cardiomyocyte gene expressions. Our work provides a 3D genome-mediated enhancer activity landscape of murine cardiomyocyte development.

Differential Sodium Current Remodelling Identifies Distinct Cellular Proarrhythmic Mechanisms in Paroxysmal vs Persistent Atrial Fibrillation

BACKGROUND

The cellular mechanisms underlying progression from paroxysmal to persistent atrial fibrillation (AF) are not fully understood, but alterations in (late) sodium current (I_Na) have been proposed. Human studies investigating electrophysiological changes at the paroxysmal stage of AF are sparse, with the majority employing right atrial appendage cardiomyocytes (CMs). We here investigated action potential (AP) characteristics and (late) I_Na remodelling in left atrial appendage CMs (LAA-CMs) from patients with paroxysmal and persistent AF and patients in sinus rhythm (SR), as well as the potential contribution of the "neuronal" sodium channel SCN10A/Na_V1.8.

METHODS

Peak I_Na, late I_Na and AP properties were investigated through patch-clamp analysis on single LAA-CMs, whereas quantitative polymerase chain reaction was used to assess SCN5A/SCN10A expression levels in LAA tissue.

RESULTS

In paroxysmal and persistent AF LAA-CMs, AP duration was shorter than in SR LAA-CMs. Compared with SR, peak I_Na and SCN5A expression were significantly decreased in paroxysmal AF, whereas they were restored to SR levels in persistent AF. Conversely, although late I_Na was unchanged in paroxysmal AF compared with SR, it was significantly increased in persistent AF. Peak or late Na_v1.8-based I_Na was not detected in persistent AF LAA-CMs. Similarly, expression of SCN10A was not observed in LAAs at any stage.

CONCLUSIONS

Our findings demonstrate differences in (late) I_Na remodeling in LAA-CMs from patients with paroxysmal vs persistent AF, indicating distinct cellular proarrhythmic mechanisms in different AF forms. These observations are of particular relevance when considering potential pharmacologic approaches targeting (late) I_Na in AF.

Brugada Syndrome: From Molecular Mechanisms and Genetics to Risk Stratification

Brugada syndrome (BrS) is a rare hereditary arrhythmia disorder, with a distinctive ECG pattern, correlated with an increased risk of ventricular arrhythmias and sudden cardiac death (SCD) in young adults. BrS is a complex entity in terms of mechanisms, genetics, diagnosis, arrhythmia risk stratification, and management. The main electrophysiological mechanism of BrS requires further research, with prevailing theories centered on aberrant repolarization, depolarization, and current-load match. Computational modelling, pre-clinical, and clinical research show that BrS molecular anomalies result in excitation wavelength (k) modifications, which eventually increase the risk of arrhythmia. Although a mutation in the SCN5A (Sodium Voltage-Gated Channel Alpha Subunit 5) gene was first reported almost two decades ago, BrS is still currently regarded as a Mendelian condition inherited in an autosomal dominant manner with incomplete penetrance, despite the recent developments in the field of genetics and the latest hypothesis of additional inheritance pathways proposing a more complex mode of inheritance. In spite of the extensive use of the next-generation sequencing (NGS) technique with high coverage, genetics remains unexplained in a number of clinically confirmed cases. Except for the SCN5A which encodes the cardiac sodium channel NaV1.5, susceptibility genes remain mostly unidentified. The predominance of cardiac transcription factor loci suggests that transcriptional regulation is essential to the Brugada syndrome's pathogenesis. It appears that BrS is a multifactorial disease, which is influenced by several loci, each of which is affected by the environment. The primary challenge in individuals with a BrS type 1 ECG is to identify those who are at risk for sudden death, researchers propose the use of a multiparametric clinical and instrumental strategy for risk stratification. The aim of this review is to summarize the latest findings addressing the genetic architecture of BrS and to provide novel perspectives into its molecular underpinnings and novel models of risk stratification.

Modeling incomplete penetrance in long QT syndrome type 3 through ion channel heterogeneity: an in silico population study

Many cardiac diseases are characterized by an increased late sodium current, including heart failure, hypertrophic cardiomyopathy, and inherited long QT syndrome type 3 (LQT3). The late sodium current in LQT3 is caused by a gain-of-function mutation in the voltage-gated sodium channel Nav1.5. Despite a well-defined genetic cause of LQT3, treatment remains inconsistent because of incomplete penetrance of the mutation and variability of antiarrhythmic efficacy. Here, we investigate the relationship between LQT3-associated mutation incomplete penetrance and variability in ion channel expression, simulating a population of 1,000 individuals with the O'Hara-Rudy model of the human ventricular myocyte. We first simulate healthy electrical activity (i.e., in the absence of a mutation) and then incorporate heterozygous expression for three LQT3-associated mutations (Y1795C, I1768V, and ΔKPQ), to directly compare the effects of each mutation on individuals across a diverse population. For all mutations, we find that susceptibility, defined by either the presence of an early afterdepolarization (EAD) or prolonged action potential duration (APD), primarily depends on the balance between the conductance of IKr and INa, for which individuals with a higher IKr-to-INa ratio are less susceptible. Furthermore, we find distinct differences across the population, observing individuals susceptible to zero, one, two, or all three mutations. Individuals tend to be less susceptible with an appropriate balance of repolarizing currents, typically via increased IKs or IK1. Interestingly, the more critical repolarizing current is mutation specific. We conclude that the balance between key currents plays a significant role in mutant-specific presentation of the disease phenotype in LQT3.NEW & NOTEWORTHY An in silico population approach investigates the relationship between variability in ion channel expression and gain-of-function mutations in the voltage-gated sodium channel associated with the congenital disorder long QT syndrome type 3 (LQT3). We find that ion channel variability can contribute to incomplete penetrance of the mutation, with mutant-specific differences in ion channel conductances leading to susceptibility to proarrhythmic action potential duration prolongation or early afterdepolarizations.

DeepLUCIA: predicting tissue-specific chromatin loops using Deep Learning-based Universal Chromatin Interaction Annotator

MOTIVATION

The importance of chromatin loops in gene regulation is broadly accepted. There are mainly two approaches to predict chromatin loops: transcription factor (TF) binding-dependent approach and genomic variation-based approach. However, neither of these approaches provides an adequate understanding of gene regulation in human tissues. To address this issue, we developed a deep learning-based chromatin loop prediction model called DeepLUCIA (Deep Learning-based Universal Chromatin Interaction Annotator).

RESULTS

Although DeepLUCIA does not use TF binding profile data which previous TF binding-dependent methods critically rely on, its prediction accuracies are comparable to those of the previous TF binding-dependent methods. More importantly, DeepLUCIA enables the tissue-specific chromatin loop predictions from tissue-specific epigenomes that cannot be handled by genomic variation-based approach. We demonstrated the utility of the DeepLUCIA by predicting several novel target genes of SNPs identified in genome-wide association studies targeting Brugada syndrome, COVID-19 severity, and age-related macular degeneration.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Variant Intronic Enhancer Controls SCN10A-short Expression and Heart Conduction

BACKGROUND

Genetic variants in SCN10A, encoding the neuronal voltage-gated sodium channel Na_V1.8, are strongly associated with atrial fibrillation, Brugada syndrome, cardiac conduction velocities, and heart rate. The cardiac function of SCN10A has not been resolved, however, and diverging mechanisms have been proposed. Here, we investigated the cardiac expression of SCN10A and the function of a variant-sensitive intronic enhancer previously linked to the regulation of SCN5A, encoding the major essential cardiac sodium channel Na_V1.5.

METHODS

The expression of SCN10A was investigated in mouse and human hearts. With the use of CRISPR/Cas9 genome editing, the mouse intronic enhancer was disrupted, and mutant mice were characterized by transcriptomic and electrophysiological analyses. The association of genetic variants at SCN5A-SCN10A enhancer regions and gene expression were evaluated by genome-wide association studies single-nucleotide polymorphism mapping and expression quantitative trait loci analysis.

RESULTS

We found that cardiomyocytes of the atria, sinoatrial node, and ventricular conduction system express a short transcript comprising the last 7 exons of the gene (Scn10a-short). Transcription occurs from an intronic enhancer-promoter complex, whereas full-length Scn10a transcript was undetectable in the human and mouse heart. Expression quantitative trait loci analysis revealed that the genetic variants in linkage disequilibrium with genetic variant rs6801957 in the intronic enhancer associate with SCN10A transcript levels in the heart. Genetic modification of the enhancer in the mouse genome led to reduced cardiac Scn10a-short expression in atria and ventricles, reduced cardiac sodium current in atrial cardiomyocytes, atrial conduction slowing and arrhythmia, whereas the expression of Scn5a, the presumed enhancer target gene, remained unaffected. In patch-clamp transfection experiments, expression of Scn10a-short-encoded Na_V1.8-short increased Na_V1.5-mediated sodium current. We propose that noncoding genetic variation modulates transcriptional regulation of Scn10a-short in cardiomyocytes that impacts Na_V1.5-mediated sodium current and heart rhythm.

CONCLUSIONS

Genetic variants in and around SCN10A modulate enhancer function and expression of a cardiac-specific SCN10A-short transcript. We propose that noncoding genetic variation modulates transcriptional regulation of a functional C-terminal portion of Na_V1.8 in cardiomyocytes that impacts on Na_V1.5 function, cardiac conduction velocities, and arrhythmia susceptibility.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

Genetic Variation in Enhancers Modifies Cardiomyopathy Gene Expression and Progression

BACKGROUND

Inherited cardiomyopathy associates with a range of phenotypes, mediated by genetic and nongenetic factors. Noninherited cardiomyopathy also displays varying progression and outcomes. Expression of cardiomyopathy genes is under the regulatory control of promoters and enhancers, and human genetic variation in promoters and enhancers may contribute to this variability.

METHODS

We superimposed epigenomic profiling from hearts and cardiomyocytes, including promoter-capture chromatin conformation information, to identify enhancers for 2 cardiomyopathy genes, MYH7 and LMNA. Enhancer function was validated in human cardiomyocytes derived from induced pluripotent stem cells. We also conducted a genome-wide search to ascertain genomic variation in enhancers positioned to alter cardiac expression and correlated one of these variants to cardiomyopathy progression using biobank data.

RESULTS

Multiple enhancers were identified and validated for LMNA and MYH7, including a key enhancer that regulates the switch from MYH6 expression to MYH7 expression. Deletion of this enhancer resulted in a dose-dependent increase in MYH6 and faster contractile rate in engineered heart tissues. We searched for genomic variation in enhancer sequences across the genome, with a focus on nucleotide changes that create or interrupt transcription factor binding sites. The sequence variant, rs875908, disrupts a T-Box Transcription Factor 5 binding motif and maps to an enhancer region 2 kilobases from the transcriptional start site of MYH7. Gene editing to remove the enhancer that harbors this variant markedly reduced MYH7 expression in human cardiomyocytes. Using biobank-derived data, rs875908 associated with longitudinal echocardiographic features of cardiomyopathy.

CONCLUSIONS

Enhancers regulate cardiomyopathy gene expression, and genomic variation within these enhancer regions associates with cardiomyopathic progression over time. This integrated approach identified noncoding modifiers of cardiomyopathy and is applicable to other cardiac genes.

Epigenetic Changes Governing Scn5a Expression in Denervated Skeletal Muscle

The SCN5A gene encodes the α-subunit of the voltage-gated cardiac sodium channel (Na_V1.5), a key player in cardiac action potential depolarization. Genetic variants in protein-coding regions of the human SCN5A have been largely associated with inherited cardiac arrhythmias. Increasing evidence also suggests that aberrant expression of the SCN5A gene could increase susceptibility to arrhythmogenic diseases, but the mechanisms governing SCN5A expression are not yet well understood. To gain insights into the molecular basis of SCN5A gene regulation, we used rat gastrocnemius muscle four days following denervation, a process well known to stimulate Scn5a expression. Our results show that denervation of rat skeletal muscle induces the expression of the adult cardiac Scn5a isoform. RNA-seq experiments reveal that denervation leads to significant changes in the transcriptome, with Scn5a amongst the fifty top upregulated genes. Consistent with this increase in expression, ChIP-qPCR assays show enrichment of H3K27ac and H3K4me3 and binding of the transcription factor Gata4 near the Scn5a promoter region. Also, Gata4 mRNA levels are significantly induced upon denervation. Genome-wide analysis of H3K27ac by ChIP-seq suggest that a super enhancer recently described to regulate Scn5a in cardiac tissue is activated in response to denervation. Altogether, our experiments reveal that similar mechanisms regulate the expression of Scn5a in denervated muscle and cardiac tissue, suggesting a conserved pathway for SCN5A expression among striated muscles.

Three-dimensional chromatin organization in cardiac development and disease

Recent technological advancements in the field of chromatin biology have rewritten the textbook on nuclear organization. We now appreciate that the folding of chromatin in the three-dimensional space (i.e. its 3D "architecture") is non-random, hierarchical, and highly complex. While 3D chromatin structure is partially encoded in the primary sequence and thereby broadly conserved across cell types and states, a substantial portion of the genome seems to be dynamic during development or in disease. Moreover, there is growing evidence that at least some of the 3D structure of chromatin is functionally linked to gene regulation, both being modulated by and impacting on multiple nuclear processes (including DNA replication, transcription, and RNA splicing). In recent years, these new concepts have nourished several investigations about the functional role of 3D chromatin topology dynamics in the heart during development and disease. This review aims to provide a comprehensive overview of our current understanding in this field, and to discuss how this knowledge can inform further research as well as clinical practice.

Genetic Dissection of a Super Enhancer Controlling the Nppa-Nppb Cluster in the Heart

RATIONALE

ANP (atrial natriuretic peptide) and BNP (B-type natriuretic peptide), encoded by the clustered genes Nppa and Nppb, are important prognostic, diagnostic, and therapeutic proteins in cardiac disease. The spatiotemporal expression pattern and stress-induction of the Nppa and Nppb are tightly regulated, possibly involving their coregulation by an evolutionary conserved enhancer cluster.

OBJECTIVE

To explore the physiological functions of the enhancer cluster and elucidate the genomic mechanism underlying Nppa-Nppb coregulation in vivo.

METHODS AND RESULTS

By analyzing epigenetic data we uncovered an enhancer cluster with super enhancer characteristics upstream of Nppb. Using CRISPR/Cas9 genome editing, the enhancer cluster or parts thereof, Nppb and flanking regions or the entire genomic block spanning Nppa-Nppb, respectively, were deleted from the mouse genome. The impact on gene regulation and phenotype of the respective mouse lines was investigated by transcriptomic, epigenomic, and phenotypic analyses. The enhancer cluster was essential for prenatal and postnatal ventricular expression of Nppa and Nppb but not of any other gene. Enhancer cluster-deficient mice showed enlarged hearts before and after birth, similar to Nppa-Nppb compound knockout mice we generated. Analysis of the other deletion alleles indicated the enhancer cluster engages the promoters of Nppa and Nppb in a competitive rather than a cooperative mode, resulting in increased Nppa expression when Nppb and flanking sequences were deleted. The enhancer cluster maintained its active epigenetic state and selectivity when its target genes are absent. In enhancer cluster-deficient animals, Nppa was induced but remained low in the postmyocardial infarction border zone and in the hypertrophic ventricle, involving regulatory sequences proximal to Nppa.

CONCLUSIONS

Coordinated ventricular expression of Nppa and Nppb is controlled in a competitive manner by a shared super enhancer, which is also required to augment stress-induced expression and to prevent premature hypertrophy.

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.

Molecular Determinants of Cardiac Arrhythmias

Cardiac arrhythmias are defined as electrical disorders of the pumping heart, including therein a wide range of physiopathological entities [...]

Epigenomic and Transcriptomic Dynamics During Human Heart Organogenesis

RATIONALE

There is growing evidence that common variants and rare sequence alterations in regulatory sequences can result in birth defects or predisposition to disease. Congenital heart defects are the most common birth defect and have a clear genetic component, yet only a third of cases can be attributed to structural variation in the genome or a mutation in a gene. The remaining unknown cases could be caused by alterations in regulatory sequences.

OBJECTIVE

Identify regulatory sequences and gene expression networks that are active during organogenesis of the human heart. Determine whether these sites and networks are enriched for disease-relevant genes and associated genetic variation.

METHODS AND RESULTS

We characterized ChromHMM (chromatin state) and gene expression dynamics during human heart organogenesis. We profiled 7 histone modifications in embryonic hearts from each of 9 distinct Carnegie stages (13-14, 16-21, and 23), annotated chromatin states, and compared these maps to over 100 human tissues and cell types. We also generated RNA-sequencing data, performed differential expression, and constructed weighted gene coexpression networks. We identified 177 412 heart enhancers; 12 395 had not been previously annotated as strong enhancers. We identified 92% of all functionally validated heart-positive enhancers (n=281; 7.5× enrichment; P<2.2×10-16). Integration of these data demonstrated novel heart enhancers are enriched near genes expressed more strongly in cardiac tissue and are enriched for variants associated with ECG measures and atrial fibrillation. Our gene expression network analysis identified gene modules strongly enriched for heart-related functions, regulatory control by heart-specific enhancers, and putative disease genes.

CONCLUSIONS

Well-connected hub genes with heart-specific expression targeted by embryonic heart-specific enhancers are likely disease candidates. Our functional annotations will allow for better interpretation of whole genome sequencing data in the large number of patients affected by congenital heart defects.

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

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

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation

Genome-wide association studies have uncovered over a 100 genetic loci associated with atrial fibrillation (AF), the most common arrhythmia. Many of the top AF-associated loci harbor key cardiac transcription factors, including PITX2, TBX5, PRRX1, and ZFHX3. Moreover, the vast majority of the AF-associated variants lie within noncoding regions of the genome where causal variants affect gene expression by altering the activity of transcription factors and the epigenetic state of chromatin. In this review, we discuss a transcriptional regulatory network model for AF defined by effector genes in Genome-wide association studies loci. We describe the current state of the field regarding the identification and function of AF-relevant gene regulatory networks, including variant regulatory elements, dose-sensitive transcription factor functionality, target genes, and epigenetic states. We illustrate how altered transcriptional networks may impact cardiomyocyte function and ionic currents that impact AF risk. Last, we identify the need for improved tools to identify and functionally test transcriptional components to define the links between genetic variation, epigenetic gene regulation, and atrial function.