A common genetic variant within SCN10A modulates cardiac SCN5A expression

Brugada syndrome in Japan and Europe: a genome-wide association study reveals shared genetic architecture and new risk loci

BACKGROUND AND AIMS

Brugada syndrome (BrS) is an inherited arrhythmia with a higher disease prevalence and more lethal arrhythmic events in Asians than in Europeans. Genome-wide association studies (GWAS) have revealed its polygenic architecture mainly in European populations. The aim of this study was to identify novel BrS-associated loci and to compare allelic effects across ancestries.

METHODS

A GWAS was conducted in Japanese participants, involving 940 cases and 1634 controls, followed by a cross-ancestry meta-analysis of Japanese and European GWAS (total of 3760 cases and 11 635 controls). The novel loci were characterized by fine-mapping, gene expression, and splicing quantitative trait associations in the human heart.

RESULTS

The Japanese-specific GWAS identified one novel locus near ZSCAN20 (P = 1.0 × 10-8), and the cross-ancestry meta-analysis identified 17 association signals, including six novel loci. The effect directions of the 17 lead variants were consistent (94.1%; P for sign test = 2.7 × 10-4), and their allelic effects were highly correlated across ancestries (Pearson's R = .91; P = 2.9 × 10-7). The genetic risk score derived from the BrS GWAS of European ancestry was significantly associated with the risk of BrS in the Japanese population [odds ratio 2.12 (95% confidence interval 1.94-2.31); P = 1.2 × 10-61], suggesting a shared genetic architecture across ancestries. Functional characterization revealed that a lead variant in CAMK2D promotes alternative splicing, resulting in an isoform switch of calmodulin kinase II-δ, favouring a pro-inflammatory/pro-death pathway.

CONCLUSIONS

This study demonstrates novel susceptibility loci implicating potentially novel pathogenesis underlying BrS. Despite differences in clinical expressivity and epidemiology, the polygenic architecture of BrS was substantially shared across ancestries.

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.

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.

Noncoding RNAs and Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes in Cardiac Arrhythmic Brugada Syndrome

Hundreds of thousands of people die each year as a result of sudden cardiac death, and many are due to heart rhythm disorders. One of the major causes of these arrhythmic events is Brugada syndrome, a cardiac channelopathy that results in abnormal cardiac conduction, severe life-threatening arrhythmias, and, on many occasions, death. This disorder has been associated with mutations and dysfunction of about two dozen genes; however, the majority of the patients do not have a definite cause for the diagnosis of Brugada Syndrome. The protein-coding genes represent only a very small fraction of the mammalian genome, and the majority of the noncoding regions of the genome are actively transcribed. Studies have shown that most of the loci associated with electrophysiological traits are located in noncoding regulatory regions and are expected to affect gene expression dosage and cardiac ion channel function. Noncoding RNAs serve an expanding number of regulatory and other functional roles within the cells, including but not limited to transcriptional, post-transcriptional, and epigenetic regulation. The major noncoding RNAs found in Brugada Syndrome include microRNAs; however, others such as long noncoding RNAs are also identified. They contribute to pathogenesis by interacting with ion channels and/or are detectable as clinical biomarkers. Stem cells have received significant attention in the recent past, and can be differentiated into many different cell types including those in the heart. In addition to contractile and relaxational properties, BrS-relevant electrophysiological phenotypes are also demonstrated in cardiomyocytes differentiated from stem cells induced from adult human cells. In this review, we discuss the current understanding of noncoding regions of the genome and their RNA biology in Brugada Syndrome. We also delve into the role of stem cells, especially human induced pluripotent stem cell-derived cardiac differentiated cells, in the investigation of Brugada syndrome in preclinical and clinical studies.

Brugada Syndrome: More than a Monogenic Channelopathy

Brugada syndrome (BrS) is an inherited cardiac channelopathy first diagnosed in 1992 but still considered a challenging disease in terms of diagnosis, arrhythmia risk prediction, pathophysiology and management. Despite about 20% of individuals carrying pathogenic variants in the SCN5A gene, the identification of a polygenic origin for BrS and the potential role of common genetic variants provide the basis for applying polygenic risk scores for individual risk prediction. The pathophysiological mechanisms are still unclear, and the initial thinking of this syndrome as a primary electrical disease is evolving towards a partly structural disease. This review focuses on the main scientific advancements in the identification of biomarkers for diagnosis, risk stratification, pathophysiology and therapy of BrS. A comprehensive model that integrates clinical and genetic factors, comorbidities, age and gender, and perhaps environmental influences may provide the opportunity to enhance patients' quality of life and improve the therapeutic approach.

Transcriptional Dysregulation Underlies Both Monogenic Arrhythmia Syndrome and Common Modifiers of Cardiac Repolarization

BACKGROUND

Brugada syndrome (BrS) is an inherited arrhythmia syndrome caused by loss-of-function variants in the cardiac sodium channel gene SCN5A (sodium voltage-gated channel alpha subunit 5) in ≈20% of subjects. We identified a family with 4 individuals diagnosed with BrS harboring the rare G145R missense variant in the cardiac transcription factor TBX5 (T-box transcription factor 5) and no SCN5A variant.

METHODS

We generated induced pluripotent stem cells (iPSCs) from 2 members of a family carrying TBX5-G145R and diagnosed with Brugada syndrome. After differentiation to iPSC-derived cardiomyocytes (iPSC-CMs), electrophysiologic characteristics were assessed by voltage- and current-clamp experiments (n=9 to 21 cells per group) and transcriptional differences by RNA sequencing (n=3 samples per group), and compared with iPSC-CMs in which G145R was corrected by CRISPR/Cas9 approaches. The role of platelet-derived growth factor (PDGF)/phosphoinositide 3-kinase (PI3K) pathway was elucidated by small molecule perturbation. The rate-corrected QT (QTc) interval association with serum PDGF was tested in the Framingham Heart Study cohort (n=1893 individuals).

RESULTS

TBX5-G145R reduced transcriptional activity and caused multiple electrophysiologic abnormalities, including decreased peak and enhanced "late" cardiac sodium current (INa), which were entirely corrected by editing G145R to wild-type. Transcriptional profiling and functional assays in genome-unedited and -edited iPSC-CMs showed direct SCN5A down-regulation caused decreased peak INa, and that reduced PDGF receptor (PDGFRA [platelet-derived growth factor receptor α]) expression and blunted signal transduction to PI3K was implicated in enhanced late INa. Tbx5 regulation of the PDGF axis increased arrhythmia risk due to disruption of PDGF signaling and was conserved in murine model systems. PDGF receptor blockade markedly prolonged normal iPSC-CM action potentials and plasma levels of PDGF in the Framingham Heart Study were inversely correlated with the QTc interval (P<0.001).

CONCLUSIONS

These results not only establish decreased SCN5A transcription by the TBX5 variant as a cause of BrS, but also reveal a new general transcriptional mechanism of arrhythmogenesis of enhanced late sodium current caused by reduced PDGF receptor-mediated PI3K signaling.

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.

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.

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.

Modulating the activity of human nociceptors with a SCN10A promoter-specific viral vector tool

Despite the high prevalence of chronic pain as a disease in our society, there is a lack of effective treatment options for patients living with this condition. Gene therapies using recombinant AAVs are a direct method to selectively express genes of interest in target cells with the potential of, in the case of nociceptors, reducing neuronal firing in pain conditions. We designed a recombinant AAV vector expressing cargos whose expression was driven by a portion of the SCN10A (Na_V1.8) promoter, which is predominantly active in nociceptors. We validated its specificity for nociceptors in mouse and human dorsal root ganglia and showed that it can drive the expression of functional proteins. Our viral vector and promoter package drove the expression of both excitatory or inhibitory DREADDs in primary human DRG cultures and in whole cell electrophysiology experiments, increased or decreased neuronal firing, respectively. Taken together, we present a novel viral tool that drives expression of cargo specifically in human nociceptors. This will allow for future specific studies of human nociceptor properties as well as pave the way for potential future gene therapies for chronic pain.

Current challenges in understanding the role of enhancers in disease

Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

The Genetics and Epigenetics of Ventricular Arrhythmias in Patients Without Structural Heart Disease

Ventricular arrhythmia without structural heart disease is an arrhythmic disorder that occurs in structurally normal heart and no transient or reversible arrhythmia factors, such as electrolyte disorders and myocardial ischemia. Ventricular arrhythmias without structural heart disease can be induced by multiple factors, including genetics and environment, which involve different genetic and epigenetic regulation. Familial genetic analysis reveals that cardiac ion-channel disorder and dysfunctional calcium handling are two major causes of this type of heart disease. Genome-wide association studies have identified some genetic susceptibility loci associated with ventricular tachycardia and ventricular fibrillation, yet relatively few loci associated with no structural heart disease. The effects of epigenetics on the ventricular arrhythmias susceptibility genes, involving non-coding RNAs, DNA methylation and other regulatory mechanisms, are gradually being revealed. This article aims to review the knowledge of ventricular arrhythmia without structural heart disease in genetics, and summarizes the current state of epigenetic regulation.

Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity

Disease-associated single-nucleotide polymorphisms (SNPs) generally do not implicate target genes, as most disease SNPs are regulatory. Many SNP-to-gene (S2G) linking strategies have been developed to link regulatory SNPs to the genes that they regulate in cis. Here, we developed a heritability-based framework for evaluating and combining different S2G strategies to optimize their informativeness for common disease risk. Our optimal combined S2G strategy (cS2G) included seven constituent S2G strategies and achieved a precision of 0.75 and a recall of 0.33, more than doubling the recall of any individual strategy. We applied cS2G to fine-mapping results for 49 UK Biobank diseases/traits to predict 5,095 causal SNP-gene-disease triplets (with S2G-derived functional interpretation) with high confidence. We further applied cS2G to provide an empirical assessment of disease omnigenicity; we determined that the top 1% of genes explained roughly half of the SNP heritability linked to all genes and that gene-level architectures vary with variant allele frequency.

Brugada Syndrome

Brugada syndrome (BrS) is an "inherited" condition characterized by predisposition to syncope and cardiac arrest, predominantly during sleep. The prevalence is ∼1:2,000, and is more commonly diagnosed in young to middle-aged males, although patient sex does not appear to impact prognosis. Despite the perception of BrS being an inherited arrhythmia syndrome, most cases are not associated with a single causative gene variant. Electrocardiogram (ECG) findings support variable extent of depolarization and repolarization changes, with coved ST-segment elevation ≥2 mm and a negative T-wave in the right precordial leads. These ECG changes are often intermittent, and may be provoked by fever or sodium channel blocker challenge. Growing evidence from cardiac imaging, epicardial ablation, and pathology studies suggests the presence of an epicardial arrhythmic substrate within the right ventricular outflow tract. Risk stratification aims to identify those who are at increased risk of sudden cardiac death, with well-established factors being the presence of spontaneous ECG changes and a history of cardiac arrest or cardiogenic syncope. Current management involves conservative measures in asymptomatic patients, including fever management and drug avoidance. Symptomatic patients typically undergo implantable cardioverter defibrillator insertion, with quinidine and epicardial ablation used for patients with recurrent arrhythmia. This review summarizes our current understanding of BrS and provides clinicians with a practical approach to diagnosis and management.

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.

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.

Asthma-associated genetic variants induce IL33 differential expression through an enhancer-blocking regulatory region

Genome-wide association studies (GWAS) have implicated the IL33 locus in asthma, but the underlying mechanisms remain unclear. Here, we identify a 5 kb region within the GWAS-defined segment that acts as an enhancer-blocking element in vivo and in vitro. Chromatin conformation capture showed that this 5 kb region loops to the IL33 promoter, potentially regulating its expression. We show that the asthma-associated single nucleotide polymorphism (SNP) rs1888909, located within the 5 kb region, is associated with IL33 gene expression in human airway epithelial cells and IL-33 protein expression in human plasma, potentially through differential binding of OCT-1 (POU2F1) to the asthma-risk allele. Our data demonstrate that asthma-associated variants at the IL33 locus mediate allele-specific regulatory activity and IL33 expression, providing a mechanism through which a regulatory SNP contributes to genetic risk of asthma.

Susceptibility to asthma and severity of symptoms are regulated by a number of different genomic regions. Here the authors characterise a 5kb regulatory region and demonstrate genetic and topological regulation of IL33 and association with disease in different human cohorts.

Early-Onset Atrial Fibrillation and the Prevalence of Rare Variants in Cardiomyopathy and Arrhythmia Genes

Question

In patients diagnosed with atrial fibrillation before 66 years of age, what is the prevalence of disease-associated variants in susceptibility genes for inherited cardiomyopathy and arrhythmia syndromes?

Findings

In this cohort study, among 1293 participants who underwent whole genome sequencing, disease-associated rare variants in cardiomyopathy and arrhythmia genes were identified in 10.1% of participants younger than 66 years and 16.8% of those younger than 30 years. Disease-associated rare variants were more prevalent in genes associated with inherited cardiomyopathy syndromes than inherited arrhythmia syndromes.

Meaning

The results of this study suggest that genetic testing in patients with early-onset atrial fibrillation identifies pathogenic variants associated with more serious inherited cardiomyopathy and arrhythmia syndromes.

Importance

Early-onset atrial fibrillation (AF) can be the initial manifestation of a more serious underlying inherited cardiomyopathy or arrhythmia syndrome.

Objective

To examine the results of genetic testing for early-onset AF.

Design, Setting, and Participants

This prospective, observational cohort study enrolled participants from an academic medical center who had AF diagnosed before 66 years of age and underwent whole genome sequencing through the National Heart, Lung, and Blood Institute’s Trans-Omics for Precision Medicine program. Participants were enrolled from November 23, 1999, to June 2, 2015. Data analysis was performed from October 24, 2020, to March 11, 2021.

Exposures

Rare variants identified in a panel of 145 genes that are included on cardiomyopathy and arrhythmia panels used by commercial clinical genetic testing laboratories.

Main Outcomes and Measures

Sequencing data were analyzed using an automated process followed by manual review by a panel of independent, blinded reviewers. The primary outcome was classification of rare variants using American College of Medical Genetics and Genomics criteria: benign, likely benign, variant of undetermined significance, likely pathogenic, or pathogenic. Disease-associated variants were defined as pathogenic/likely pathogenic variants in genes associated with autosomal dominant or X-linked dominant disorders.

Results

Among 1293 participants (934 [72.2%] male; median [interquartile range] age at enrollment, 56 [48-61] years; median [interquartile range] age at AF diagnosis, 50 [41-56] years), genetic testing identified 131 participants (10.1%) with a disease-associated variant, 812 (62.8%) with a variant of undetermined significance, 92 (7.1%) as heterozygous carriers for an autosomal recessive disorder, and 258 (20.0%) with no suspicious variant. The likelihood of a disease-associated variant was highest in participants with AF diagnosed before the age of 30 years (20 of 119 [16.8%; 95% CI, 10.0%-23.6%]) and lowest after the age of 60 years (8 of 112 [7.1%; 95% CI, 2.4%-11.9%]). Disease-associated variants were more often associated with inherited cardiomyopathy syndromes compared with inherited arrhythmias. The most common genes were TTN (n = 38), MYH7 (n = 18), MYH6 (n = 10), LMNA (n = 9), and KCNQ1 (n = 8).

Conclusions and Relevance

In this cohort study, genetic testing identified a disease-associated variant in 10% of patients with early-onset AF (the percentage was higher if diagnosed before the age of 30 years and lower if diagnosed after the age of 60 years). Most pathogenic/likely pathogenic variants are in genes associated with cardiomyopathy. These results support the use of genetic testing in early-onset AF.

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.

Common variants in SCN10A gene associated with Brugada syndrome

BACKGROUND

Genome-wide association studies indicate that SCN10A plays an important role in cardiac electrophysiology. Common and rare SCN10A variants are suggested to contribute to Brugada Syndrome (BrS), an inherited channelopathy resulting from genetic-determined loss-of-function in cardiac sodium channel. This study sought to characterize the role of SCN10A common variants in BrS.

METHODS AND RESULTS

Clinical and genetic analyses were performed in 197 patients diagnosed with BrS. Baseline ECG parameters were evaluated in patients carrying each of four common variants associated with BrS. Cellular electrophysiological study was performed in SCN5A-SCN10A co-transfected TSA201 cells to investigate the possible electrophysiological characteristics of the allele of rs6795970, which displayed the most significant association with BrS. Four SCN10A common variants (rs7630989, rs57326399, rs6795970, rs12632942) displayed significant association with BrS susceptibility. There were no evident associations between baseline ECG parameters in BrS patients and the different genotypes of the four variants. Rs6795970 (V1073) was strongly associated with a risk for BrS, which suggests the different electrophysiological characters between these two alleles. Functional study showed a positive shift in steady-state activation (V1/2: -62.2 ± 2.6 vs. -53.5 ± 1.6 for A1073 and V1073 group, respectively; P < 0.05) and slower recovery from inactivation in mutant SCN5A-SCN10A co-transfected cells with, which contribute to the slow conduction in BrS patients with rs6795970.

CONCLUSIONS

SCN10A common variants are associated with increased susceptibility to BrS. An allele rs6795970 (V1073) increases the risk for BrS. The electrophysiological changes in a positive shift in steady-state activation and slower recovery from inactivation by SCN10A-V1073 contribute to this variant associated BrS.

Non-Coding RNAs in the Cardiac Action Potential and Their Impact on Arrhythmogenic Cardiac Diseases

Cardiac arrhythmias are prevalent among humans across all age ranges, affecting millions of people worldwide. While cardiac arrhythmias vary widely in their clinical presentation, they possess shared complex electrophysiologic properties at cellular level that have not been fully studied. Over the last decade, our current understanding of the functional roles of non-coding RNAs have progressively increased. microRNAs represent the most studied type of small ncRNAs and it has been demonstrated that miRNAs play essential roles in multiple biological contexts, including normal development and diseases. In this review, we provide a comprehensive analysis of the functional contribution of non-coding RNAs, primarily microRNAs, to the normal configuration of the cardiac action potential, as well as their association to distinct types of arrhythmogenic cardiac diseases.

Genetic predictors of sick sinus syndrome

Sick sinus syndrome (SSS) encompasses a group of conduction disorders characterized by the inability of sinoatrial node to perform its pacemaker function. Our aim was to identify genetic predictors of SSS in a prospective cohort of patients admitted to the clinic for pacemaker implantation using single-locus and multilocus approaches. We performed genotyping for polymorphic markers of CLCNKA (rs10927887), SCN10A (rs6795970), FNDC3B (rs9647379), MIR146A (rs2910164), SYT10 (rs7980799), MYH6 (rs365990), and KCNE1 (rs1805127) genes in the group of 284 patients with SSS and 243 healthy individuals. Associations between the studied loci and SSS were tested using logistic regression under recessive genetic model using sex and age as covariates. Multilocus analysis was performed using Markov chain Monte Carlo method implemented in the APSampler program. Correction for multiple testing was performed using Benjamini-Hochberg procedure. We detected an individual association between KCNE1 rs1805127*A allele and SSS in the total study group (OR 0.43, P_FDR = 0.028) and in the subgroup of patients with 2nd or 3rd degree sinoatrial block (OR 0.17, P_FDR = 0.033), and identified seven allelic patterns associated with the disease. SCN10A rs6795970*T and MIR146A rs2910164*C alleles were present in all seven combinations associated with SSS. The highest risk of SSS was conferred by the combination SCN10A rs6795970*T+FNDC3B rs9647379*C+MIR146A rs2910164*C+SYT10 rs7980799*C+KCNE1 rs1805127*G (OR 2.98, CI 1.77-5.00, P = 1.27 × 10^-5, P_FDR = 0.022). Our findings suggest that KCNE1 rs1805127 polymorphism may play a role in susceptibility to sinoatrial node dysfunction, particularly presenting as 2nd or 3rd degree sinoatrial block, and the risk-modifying effect of other studied loci is better detected using multilocus approach.

WEVar: a novel statistical learning framework for predicting noncoding regulatory variants

Understanding the functional consequence of noncoding variants is of great interest. Though genome-wide association studies or quantitative trait locus analyses have identified variants associated with traits or molecular phenotypes, most of them are located in the noncoding regions, making the identification of causal variants a particular challenge. Existing computational approaches developed for prioritizing noncoding variants produce inconsistent and even conflicting results. To address these challenges, we propose a novel statistical learning framework, which directly integrates the precomputed functional scores from representative scoring methods. It will maximize the usage of integrated methods by automatically learning the relative contribution of each method and produce an ensemble score as the final prediction. The framework consists of two modes. The first 'context-free' mode is trained using curated causal regulatory variants from a wide range of context and is applicable to predict regulatory variants of unknown and diverse context. The second 'context-dependent' mode further improves the prediction when the training and testing variants are from the same context. By evaluating the framework via both simulation and empirical studies, we demonstrate that it outperforms integrated scoring methods and the ensemble score successfully prioritizes experimentally validated regulatory variants in multiple risk loci.

Cardiac cell type-specific gene regulatory programs and disease risk association

Misregulated gene expression in human hearts can result in cardiovascular diseases that are leading causes of mortality worldwide. However, the limited information on the genomic location of candidate cis-regulatory elements (cCREs) such as enhancers and promoters in distinct cardiac cell types has restricted the understanding of these diseases. Here, we defined >287,000 cCREs in the four chambers of the human heart at single-cell resolution, which revealed cCREs and candidate transcription factors associated with cardiac cell types in a region-dependent manner and during heart failure. We further found cardiovascular disease-associated genetic variants enriched within these cCREs including 38 candidate causal atrial fibrillation variants localized to cardiomyocyte cCREs. Additional functional studies revealed that two of these variants affect a cCRE controlling KCNH2/HERG expression and action potential repolarization. Overall, this atlas of human cardiac cCREs provides the foundation for illuminating cell type-specific gene regulation in human hearts during health and disease.

Assessment of the causal relevance of ECG parameters for risk of atrial fibrillation: A mendelian randomisation study

BACKGROUND

Atrial electrical and structural remodelling in older individuals with cardiovascular risk factors has been associated with changes in surface electrocardiographic (ECG) parameters (e.g., prolongation of the PR interval) and higher risks of atrial fibrillation (AF). However, it has been difficult to establish whether altered ECG parameters are the cause or a consequence of the myocardial substrate leading to AF. This study aimed to examine the potential causal relevance of ECG parameters on risk of AF using mendelian randomisation (MR).

METHODS AND FINDINGS

Weighted genetic scores explaining lifelong differences in P-wave duration, PR interval, and QT interval were constructed, and associations between these ECG scores and risk of AF were estimated among 278,792 UK Biobank participants (mean age: 57 years at recruitment; 19,132 AF cases). The independent genetic variants contributing to each of the separate ECG scores, and their corresponding weights, were based on published genome-wide association studies. In UK Biobank, genetic scores representing a 5 ms longer P-wave duration or PR interval were significantly associated with lower risks of AF (odds ratio [OR] 0.91; 95% confidence interval [CI]: 0.87-0.96, P = 2 × 10-4 and OR 0.94; 95% CI: 0.93-0.96, P = 2 × 10-19, respectively), while longer QT interval was not significantly associated with AF. These effects were independently replicated among a further 17,931 AF cases from the AFGen Consortium. Investigation of potential mechanistic pathways showed that differences in ECG parameters associated with specific ion channel genes had effects on risk of AF consistent with the overall scores, while the overall scores were not associated with changes in left atrial size. Limitations of the study included the inherent assumptions of MR, restriction to individuals of European ancestry, and possible restriction of results to the normal ECG ranges represented in UK Biobank.

CONCLUSIONS

In UK Biobank, we observed evidence suggesting a causal relationship between lifelong differences in ECG parameters (particularly PR interval) that reflect longer atrial conduction times and a lower risk of AF. These findings, which appear to be independent of atrial size and concomitant cardiovascular comorbidity, support the relevance of varying mechanisms underpinning AF and indicate that more individualised treatment strategies warrant consideration.

The physiological function of different voltage-gated sodium channels in pain

Evidence from human genetic pain disorders shows that voltage-gated sodium channel α-subtypes Nav1.7, Nav1.8 and Nav1.9 are important in the peripheral signalling of pain. Nav1.7 is of particular interest because individuals with Nav1.7 loss-of-function mutations are congenitally insensitive to acute and chronic pain, and there is considerable hope that phenocopying these effects with a pharmacological antagonist will produce a new class of analgesic drug. However, studies in these rare individuals do not reveal how and where voltage-gated sodium channels contribute to pain signalling, which is of critical importance for drug development. More than a decade of research utilizing rodent genetic models and pharmacological tools to study voltage-gated sodium channels in pain has begun to unravel the role of different subtypes. Here, we review the contribution of individual channel subtypes in three key physiological processes necessary for transmission of sensory information to the CNS: transduction of stimuli at peripheral nerve terminals, axonal transmission of action potentials and neurotransmitter release from central terminals. These data suggest that drugs seeking to recapitulate the analgesic effects of loss of function of Nav1.7 will need to be brain-penetrant - which most of those developed to date are not.

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.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

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.

A cellular model of Brugada syndrome with SCN10A variants using human-induced pluripotent stem cell-derived cardiomyocytes

AIMS

Brugada syndrome (BrS) is associated with a pronounced risk to develop sudden cardiac death (SCD). Up to 21% of patients are related to mutations in SCN5A. Studies identified SCN10A as a contributor of BrS. However, the investigation of the human cellular phenotype of BrS in the presence of SCN10A mutations remains lacking. The objective of this study was to establish a cellular model of BrS in presence of SCN10A mutations using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs).

METHODS AND RESULTS

Dermal fibroblasts obtained from a BrS patient suffering from SCD harbouring the SCN10A double variants (c.3803G>A and c.3749G>A) and three independent healthy control subjects were reprogrammed to hiPSCs. Human-induced pluripotent stem cells were differentiated into cardiomyocytes (hiPSC-CMs).The hiPSC-CMs from the BrS patient showed a significantly reduced peak sodium channel current (INa) and a significantly reduced ATX II (sea anemone toxin, an enhancer of late INa) sensitive as well as A-887826 (a blocker of SCN10A channel) sensitive late sodium channel current (INa) when compared with the healthy control hiPSC-CMs, indicating loss-of-function of sodium channels. Consistent with reduced INa the action potential amplitude and upstroke velocity (Vmax) were significantly reduced, which may contribute to arrhythmogenesis of BrS. Moreover, Ajmaline effects on action potentials were stronger in BrS-hiPSC-CMs than in healthy control cells. This is in agreement with the higher susceptibility of patients to sodium channel blocking drugs in unmasking BrS.

CONCLUSION

Patient-specific hiPSC-CMs are able to recapitulate single-cell phenotype features of BrS with SCN10A mutations and may provide novel opportunities to further elucidate the cellular disease mechanism.

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.

Comparable clinical characteristics in Brugada syndrome patients harboring SCN5A or novel SCN10A variants

AIMS

The Brugada syndrome (BrS) is an inherited disease associated with an increased risk of sudden cardiac death. Often, the genetic cause remains undetected. Perhaps due at least in part because the NaV1.8 protein is expressed more in both the central and peripheral nervous systems than in the heart, the SCN10A gene is not included in diagnostic arrhythmia/sudden death panels in the vast majority of cardiogenetics centres.

METHODS AND RESULTS

Clinical characteristics were assessed in patients harboring either SCN5A or novel SCN10A variants. Genetic testing was performed using Next Generation Sequencing on genomic DNA. Clinical characteristics, including the arrhythmogenic substrate, in BrS patients harboring novel SCN10A variants and SCN5A variants are comparable. Clinical characteristics, including gender, age, personal history of cardiac arrest/syncope, spontaneous BrS electrocardiogram pattern, family history of sudden death, and arrhythmic substrate are not significantly different between probands harboring SCN10A or SCN5A variants.

CONCLUSION

Future studies are warranted to further characterize the role of these specific SCN10A variants.

Predicting cardiac electrical response to sodium-channel blockade and Brugada syndrome using polygenic risk scores

Aims

Sodium-channel blockers (SCBs) are associated with arrhythmia, but variability of cardiac electrical response remains unexplained. We sought to identify predictors of ajmaline-induced PR and QRS changes and Type I Brugada syndrome (BrS) electrocardiogram (ECG).

Methods and results

In 1368 patients that underwent ajmaline infusion for suspected BrS, we performed measurements of 26 721 ECGs, dose–response mixed modelling and genotyping. We calculated polygenic risk scores (PRS) for PR interval (PRS_PR), QRS duration (PRS_QRS), and Brugada syndrome (PRS_BrS) derived from published genome-wide association studies and used regression analysis to identify predictors of ajmaline dose related PR change (slope) and QRS slope. We derived and validated using bootstrapping a predictive model for ajmaline-induced Type I BrS ECG. Higher PRS_PR, baseline PR, and female sex are associated with more pronounced PR slope, while PRS_QRS and age are positively associated with QRS slope (P < 0.01 for all). PRS_BrS, baseline QRS duration, presence of Type II or III BrS ECG at baseline, and family history of BrS are independently associated with the occurrence of a Type I BrS ECG, with good predictive accuracy (optimism-corrected C-statistic 0.74).

Conclusion

We show for the first time that genetic factors underlie the variability of cardiac electrical response to SCB. PRS_BrS, family history, and a baseline ECG can predict the development of a diagnostic drug-induced Type I BrS ECG with clinically relevant accuracy. These findings could lead to the use of PRS in the diagnosis of BrS and, if confirmed in population studies, to identify patients at risk for toxicity when given SCB.

Update on Genetic Basis of Brugada Syndrome: Monogenic, Polygenic or Oligogenic?

Brugada syndrome is a rare inherited arrhythmogenic disease leading to ventricular fibrillation and high risk of sudden death. In 1998, this syndrome was linked with a genetic variant with an autosomal dominant pattern of inheritance. To date, rare variants identified in more than 40 genes have been potentially associated with this disease. Variants in regulatory regions, combinations of common variants and other genetic alterations are also proposed as potential origins of Brugada syndrome, suggesting a polygenic or oligogenic inheritance pattern. However, most of these genetic alterations remain of questionable causality; indeed, rare pathogenic variants in the SCN5A gene are the only established cause of Brugada syndrome. Comprehensive analysis of all reported genetic alterations identified the origin of disease in no more than 40% of diagnosed cases. Therefore, identifying the cause of this rare arrhythmogenic disease in the many families without a genetic diagnosis is a major current challenge in Brugada syndrome. Additional challenges are interpretation/classification of variants and translation of genetic data into clinical practice. Further studies focused on unraveling the pathophysiological mechanisms underlying the disease are needed. Here we provide an update on the genetic basis of Brugada syndrome.

Identification and Characterization of a Transcribed Distal Enhancer Involved in Cardiac Kcnh2 Regulation

The human ether-a-go-go-related gene KCNH2 encodes the voltage-gated potassium channel underlying I_Kr, a current critical for the repolarization phase of the cardiac action potential. Mutations in KCNH2 that cause a reduction of the repolarizing current can result in cardiac arrhythmias associated with long-QT syndrome. Here, we investigate the regulation of KCNH2 and identify multiple active enhancers. A transcribed enhancer ∼85 kbp downstream of Kcnh2 physically contacts the promoters of two Kcnh2 isoforms in a cardiac-specific manner in vivo. Knockdown of its ncRNA transcript results in reduced expression of Kcnh2b and two neighboring mRNAs, Nos3 and Abcb8, in vitro. Genomic deletion of the enhancer, including the ncRNA transcription start site, from the mouse genome causes a modest downregulation of both Kcnh2a and Kcnh2b in the ventricles. These findings establish that the regulation of Kcnh2a and Kcnh2b is governed by a complex regulatory landscape that involves multiple partially redundantly acting enhancers.

Susceptibility to Ventricular Arrhythmias Resulting from Mutations in FKBP1B, PXDNL, and SCN9A Evaluated in hiPSC Cardiomyocytes

BACKGROUND

We report an inherited cardiac arrhythmia syndrome consisting of Brugada and Early Repolarization Syndrome associated with variants in SCN9A, PXDNL, and FKBP1B. The proband inherited the 3 mutations and exhibited palpitations and arrhythmia-mediated syncope, whereas the parents and sister, who carried one or two of the mutations, were asymptomatic.

METHODS AND RESULTS

We assessed the functional impact of these mutations in induced pluripotent stem cell cardiomyocytes (hiPSC-CMs) derived from the proband and an unaffected family member. Current and voltage clamp recordings, as well as confocal microscopy analysis of Ca2+ transients, were evaluated in hiPSC-CMs from the proband and compared these results with hiPSC-CMs from undiseased controls. Genetic analysis using next-generation DNA sequencing revealed heterozygous mutations in SCN9A, PXDNL, and FKBP1B in the proband. The proband displayed right bundle branch block and exhibited episodes of syncope. The father carried a mutation in FKBP1B, whereas the mother and sister carried the SCN9A mutation. None of the 3 family members screened developed cardiac events. Action potential recordings from control hiPSC-CM showed spontaneous activity and a low upstroke velocity. In contrast, the hiPSC-CM from the proband showed irregular spontaneous activity. Confocal microscopy of the hiPSC-CM of the proband revealed low fluorescence intensity Ca2+ transients that were episodic in nature. Patch-clamp measurements in hiPSC-CM showed no difference in I Na but reduced I Ca in the proband compared with control. Coexpression of PXDNL-R391Q with SCN5A-WT displayed lower I Na density compared to PXDNL-WT. In addition, coexpression of PXDNL-R391Q with KCND3-WT displayed significantly higher I to density compared to PXDNL-WT.

CONCLUSION

SCN9A, PXDNL, and FKBP1B variants appeared to alter spontaneous activity in hiPSC-CM. Only the proband carrying all 3 mutations displayed the ERS/BrS phenotype, whereas one nor two mutations alone did not produce the clinical phenotype. Our results suggest a polygenic cause of the BrS/ERS arrhythmic phenotype due to mutations in these three gene variants caused a very significant loss of function of I Na and I Ca and gain of function of I to.

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.

The contribution of non-coding regulatory elements to cardiovascular disease

Cardiovascular disease collectively accounts for a quarter of deaths worldwide. Genome-wide association studies across a range of cardiovascular traits and pathologies have highlighted the prevalence of common non-coding genetic variants within candidate loci. Here, we review genetic, epigenomic and molecular approaches to investigate the contribution of non-coding regulatory elements in cardiovascular biology. We then discuss recent insights on the emerging role of non-coding variation in predisposition to cardiovascular disease, with a focus on novel mechanistic examples from functional genomics studies. Lastly, we consider the clinical significance of these findings at present, and some of the current challenges facing the field.

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.

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.

From Genome-Wide Association Studies to Cardiac Electrophysiology: Through the Maze of Biological Complexity

Genome Wide Association Studies (GWAS) have provided an enormous amount of data on genomic loci associated with cardiac electrophysiology and arrhythmias. Clinical relevance, however, remains unclear since GWAS do not provide a mechanistic explanation for this association. Determining the electrophysiological relevance of variants for arrhythmias would aid development of risk stratification models for patients with arrhythmias. In this review, we give an overview of genetic variants related to ECG intervals and arrhythmogenic pathologies and discuss how these variants may influence cardiac electrophysiology and the occurrence of arrhythmias.

Identification of Functional Variant Enhancers Associated With Atrial Fibrillation

RATIONALE

Genome-wide association studies have identified a large number of common variants (single-nucleotide polymorphisms) associated with atrial fibrillation (AF). These variants are located mainly in noncoding regions of the genome and likely include variants that modulate the function of transcriptional regulatory elements (REs) such as enhancers. However, the actual REs modulated by variants and the target genes of such REs remain to be identified. Thus, the biological mechanisms by which genetic variation promotes AF has thus far remained largely unexplored.

OBJECTIVE

To identify REs in genome-wide association study loci that are influenced by AF-associated variants.

METHODS AND RESULTS

We screened 2.45 Mbp of human genomic DNA containing 12 strongly AF-associated loci for RE activity using self-transcribing active regulatory region sequencing and a recently generated monoclonal line of conditionally immortalized rat atrial myocytes. We identified 444 potential REs, 55 of which contain AF-associated variants (P<10^-8). Subsequently, using an adaptation of the self-transcribing active regulatory region sequencing approach, we identified 24 variant REs with allele-specific regulatory activity. By mining available chromatin conformation data, the possible target genes of these REs were mapped. To define the physiological function and target genes of such REs, we deleted the orthologue of an RE containing noncoding variants in the Hcn4 (potassium/sodium hyperpolarization-activated cyclic nucleotide-gated channel 4) locus of the mouse genome. Mice heterozygous for the RE deletion showed bradycardia, sinus node dysfunction, and selective loss of Hcn4 expression.

CONCLUSIONS

We have identified REs at multiple genetic loci for AF and found that loss of an RE at the HCN4 locus results in sinus node dysfunction and reduced gene expression. Our approach can be broadly applied to facilitate the identification of human disease-relevant REs and target genes at cardiovascular genome-wide association studies loci.

The evolution of gene-guided management of inherited arrhythmia syndromes: Peering beyond monogenic paradigms towards comprehensive genomic risk scores

Inherited arrhythmia syndromes have traditionally been viewed as monogenic forms of disease whose pathophysiology is driven by a single highly penetrant rare genetic variant. Although an accurate depiction of a proportion of genetic variants, the variable penetrance frequently noted in genotype positive families and the presence of sporadic genotype negative cases have long highlighted a more nuanced truth being operative. Coupled with our more recent recognition that many rare variants implicated in inherited arrhythmia syndromes possess unexpectedly high allele frequencies within the general population, these observations have contributed to the realization that a spectrum of pathogenicity exists among clinically relevant genetic variants. Notably, variable mutation pathogenicity and corresponding variable degrees of penetrance emphasize a limitation of contemporary guidelines, which attempt to dichotomize genetic variants as pathogenic or benign. Recognition of the existence of low and intermediate penetrant variants insufficient to be causative for disease in isolation has served to emphasize the importance of additional genetic, clinical, and environmental factors in the pathogenesis of rare inherited arrhythmia syndromes. Despite being rare, it has also become increasingly evident that common genetic variants play critical roles in both heritable channelopathies and cardiomyopathies and in aggregate may even be the primary drivers in certain instances, such as genotype negative Brugada syndrome. Our growing realization that the genetic substrates of inherited arrhythmia syndromes have intricacies that extend beyond traditionally perceived monogenic paradigms has highlighted a potential value of leveraging more comprehensive genomic risk scores for predicting disease development and arrhythmic risk.