PR interval genome-wide association meta-analysis identifies 50 loci associated with atrial and atrioventricular electrical activity

SCN10A gene polymorphism is associated with pain sensitivity and postoperative analgesic effects in patients undergoing gynecological laparoscopy

BACKGROUND

Postoperative pain intensity is influenced by various factors, including genetic variations. The SCN10A gene encodes the Nav1.8 sodium channel protein, which is crucial for pain signal transmission in peripheral sensory neurons.

OBJECTIVES

This study aims to investigate the relationship between genetic mutations in the SCN10A gene (rs6795970) and postoperative analgesic effects following gynecological laparoscopic surgery.

METHODS

Two hundred female patients undergoing gynecological laparoscopic surgery under general anesthesia were included. pain sensitivity was evaluated using the catastrophizing scale and pain sensitivity questionnaire (PSQ). Patients received patient-controlled intravenous analgesia with sufentanil and dexmedetomidine for 48 h post-surgery. Postoperative pain indicators, such as visual analog scale (VAS) scores, Ramsay scores, and side effects were recorded. SCN10A rs6795970 mutations were identified using MassARRAY SNP typing technology, and patients were categoried into homozygous mutant (AA), wild type (GG), and heterozygous mutation (GA) groups for analysis.

RESULTS

Patients in the AA group had higher scores on the pain Catastrophizing Scale, PSQ-total, PSQ-minor, and PSQ-moderate compared to GA and GG groups (P < 0.05). VAS scores at 4, 6, and 12 h post-operation were higher in the AA group than the GG group (P < 0.05). Ramsay scores were lower in AA patients at 2 and 4 h post-operation compared to GA and GG groups (P < 0.05). The AA group exhibited more self-control analgesic pump compressions within the first 24 h post-surgery, quicker depletion of analgesics in the pump, and lower patient satisfaction with pain relief compared to GA and GG groups (P < 0.05).

CONCLUSIONS

Female patients with homozygous SCN10A mutations may experience higher preoperative pain scores and increased sensitivity to postoperative pain following gynecological laparoscopic surgery with intravenous patient-controlled analgesia.

TRIAL REGISTRATION

www.chictr.org.cn , registration number: ChiCTR2200062425.

Atrial cardiomyopathy revisited-evolution of a concept: a clinical consensus statement of the European Heart Rhythm Association (EHRA) of the ESC, the Heart Rhythm Society (HRS), the Asian Pacific Heart Rhythm Society (APHRS), and the Latin American Heart Rhythm Society (LAHRS)

AIMS

The concept of "atrial cardiomyopathy" (AtCM) had been percolating through the literature since its first mention in 1972. Since then, publications using the term were sporadic until the decision was made to convene an expert working group with representation from four multinational arrhythmia organizations to prepare a consensus document on atrial cardiomyopathy in 2016 (EHRA/HRS/APHRS/SOLAECE expert consensus on atrial cardiomyopathies: definition, characterization, and clinical implication). Subsequently, publications on AtCM have increased progressively.

METHODS AND RESULTS

The present consensus document elaborates the 2016 AtCM document further to implement a simple AtCM staging system (AtCM stages 1-3) by integrating biomarkers, atrial geometry, and electrophysiological changes. However, the proposed AtCM staging needs clinical validation. Importantly, it is clearly stated that the presence of AtCM might serve as a substrate for the development of atrial fibrillation (AF) and AF may accelerates AtCM substantially, but AtCM per se needs to be viewed as a separate entity.

CONCLUSION

Thus, the present document serves as a clinical consensus statement of the European Heart Rhythm Association (EHRA) of the ESC, the Heart Rhythm Society (HRS), the Asian Pacific Heart Rhythm Society (APHRS), and the Latin American Heart Rhythm Society (LAHRS) to contribute to the evolution of the AtCM concept.

Advancing drug development for atrial fibrillation by prioritising findings from human genetic association studies

BACKGROUND

Drug development for atrial fibrillation (AF) has failed to yield new approved compounds. We sought to identify and prioritise potential druggable targets with support from human genetics, by integrating the available evidence with bioinformatics sources relevant for AF drug development.

METHODS

Genetic hits for AF and related traits were identified through structured search of MEDLINE. Genes derived from each paper were cross-referenced with the OpenTargets platform for drug interactions. Confirmation/validation was demonstrated through structured searches and review of evidence on MEDLINE and ClinialTrials.gov for each drug and its association with AF.

FINDINGS

613 unique drugs were identified, with 21 already included in AF Guidelines. Cardiovascular drugs from classes not currently used for AF (e.g. ranolazine and carperitide) and anti-inflammatory drugs (e.g. dexamethasone and mehylprednisolone) had evidence of potential benefit. Further targets were considered druggable but remain open for drug development.

INTERPRETATION

Our systematic approach, combining evidence from different bioinformatics platforms, identified drug repurposing opportunities and druggable targets for AF.

FUNDING

KK is supported by Barts Charity grant G-002089 and is mentored on the AFGen 2023-24 Fellowship funded by the AFGen NIH/NHLBI grant R01HL092577. RP is supported by the UCL BHF Research Accelerator AA/18/6/34223 and NIHR grant NIHR129463. AFS is supported by the BHF grants PG/18/5033837, PG/22/10989 and UCL BHF Accelerator AA/18/6/34223 as well as the UK Research and Innovation (UKRI) under the UK government's Horizon Europe funding guarantee EP/Z000211/1 and by the UKRI-NIHR grant MR/V033867/1 for the Multimorbidity Mechanism and Therapeutics Research Collaboration. AF is supported by UCL BHF Accelerator AA/18/6/34223. CF is supported by UCL BHF Accelerator AA/18/6/34223.

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.

PR prolongation as a predictor of atrial fibrillation onset: A state-of-the-art review

First-degree atrioventricular block (1-AVB), characterized by a PR interval exceeding 200 milliseconds, has traditionally been perceived as a benign cardiac condition. Recently, this perception has been challenged by investigations that indicate a potential association between PR prolongation and an elevated risk of atrial fibrillation (AF). To consolidate these findings, we performed a comprehensive review to assess the available evidence indicating a relationship between these two conditions. We searched MEDLINE and EMBASE databases as well as manually searched references of retrieved articles. We selected 18 cohort studies/meta-analyses involving general and special populations. Consistent findings across expansive cohort studies reveal that incremental increases in the PR interval may serve as an independent risk factor for AF. However, our analyses underscore the need for further research into the association between 1-AVB, defined by a specified PR interval cutoff, and the risk of AF.

Improvement in Atrioventricular Conduction Using Cardioneuroablation Performed Immediately after Pulmonary Vein Isolation

In patients with atrial fibrillation (AF) recurrences after pulmonary vein isolation (PVI), concomitant treatment using anti arrhythmic drugs (AADs) can lead to clinical success. However, patients with atrioventricular (AV) block may not be good candidates for concomitant AAD therapy due to the risk of further worsening of conduction abnormalities. Cardioneuroablation (CNA), as an adjunct to PVI, may offer a solution to this problem. We present a case of a 74-year-old male with paroxysmal AF and first degree AV block in whom CNA following PVI led to PR normalization. The presented case describes an example of CNA utilization in patients with AF undergoing PVI who have concomitant problems with AV conduction and shows that CNA can be sometimes useful in older patients with functional AV block.

Development of the Cardiac Conduction System

The electrical impulses that coordinate the sequential, rhythmic contractions of the atria and ventricles are initiated and tightly regulated by the specialized tissues of the cardiac conduction system. In the mature heart, these impulses are generated by the pacemaker cardiomyocytes of the sinoatrial node, propagated through the atria to the atrioventricular node where they are delayed and then rapidly propagated to the atrioventricular bundle, right and left bundle branches, and finally, the peripheral ventricular conduction system. Each of these specialized components arise by complex patterning events during embryonic development. This chapter addresses the origins and transcriptional networks and signaling pathways that drive the development and maintain the function of the cardiac conduction system.

Pericardial Delivery of SDF-1α Puerarin Hydrogel Promotes Heart Repair and Electrical Coupling

The stromal-derived factor 1α/chemokine receptor 4 (SDF-1α/CXCR4) axis contributes to myocardial protection after myocardial infarction (MI) by recruiting endogenous stem cells into the ischemic tissue. However, excessive inflammatory macrophages are also recruited simultaneously, aggravating myocardial damage. More seriously, the increased inflammation contributes to abnormal cardiomyocyte electrical coupling, leading to inhomogeneities in ventricular conduction and retarded conduction velocity. It is highly desirable to selectively recruit the stem cells but block the inflammation. In this work, SDF-1α-encapsulated Puerarin (PUE) hydrogel (SDF-1α@PUE) is capable of enhancing endogenous stem cell homing and simultaneously polarizing the recruited monocyte/macrophages into a repairing phenotype. Flow cytometry analysis of the treated heart tissue shows that endogenous bone marrow mesenchymal stem cells, hemopoietic stem cells, and immune cells are recruited while SDF-1α@PUE efficiently polarizes the recruited monocytes/macrophages into the M2 type. These macrophages influence the preservation of connexin 43 (Cx43) expression which modulates intercellular coupling and improves electrical conduction. Furthermore, by taking advantage of the improved "soil", the recruited stem cells mediate an improved cardiac function by preventing deterioration, promoting neovascular architecture, and reducing infarct size. These findings demonstrate a promising therapeutic platform for MI that not only facilitates heart regeneration but also reduces the risk of cardiac arrhythmias.

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.

Phenotypic and Genetic Factors Associated with Absence of Cardiomyopathy Symptoms in PLN:c.40_42delAGA Carriers

The c.40_42delAGA variant in the phospholamban gene (PLN) has been associated with dilated and arrhythmogenic cardiomyopathy, with up to 70% of carriers experiencing a major cardiac event by age 70. However, there are carriers who remain asymptomatic at older ages. To understand the mechanisms behind this incomplete penetrance, we evaluated potential phenotypic and genetic modifiers in 74 PLN:c.40_42delAGA carriers identified in 36,339 participants of the Lifelines population cohort. Asymptomatic carriers (N = 48) showed shorter QRS duration (- 5.73 ms, q value = 0.001) compared to asymptomatic non-carriers, an effect we could replicate in two different independent cohorts. Furthermore, symptomatic carriers showed a higher correlation (rPearson = 0.17) between polygenic predisposition to higher QRS (PGSQRS) and QRS (p value = 1.98 × 10-8), suggesting that the effect of the genetic variation on cardiac rhythm might be increased in symptomatic carriers. Our results allow for improved clinical interpretation for asymptomatic carriers, while our approach could guide future studies on genetic diseases with incomplete penetrance.

Genetic insights into resting heart rate and its role in cardiovascular disease

Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.

Genome-wide association studies: utility and limitations for research in physiology

Physiological systems are subject to interindividual variation encoded by genetics. Genome-wide association studies (GWAS) operate by surveying thousands of genetic variants from a substantial number of individuals and assessing their association to a trait of interest, be it a physiological variable, a molecular phenotype (e.g. gene expression), or even a disease or condition. Through a myriad of methods, GWAS downstream analyses then explore the functional consequences of each variant and attempt to ascertain a causal relationship to the phenotype of interest, as well as to delve into its links to other traits. This type of investigation allows mechanistic insights into physiological functions, pathological disturbances and shared biological processes between traits (i.e. pleiotropy). An exciting example is the discovery of a new thyroid hormone transporter (SLC17A4) and hormone metabolising enzyme (AADAT) from a GWAS on free thyroxine levels. Therefore, GWAS have substantially contributed with insights into physiology and have been shown to be useful in unveiling the genetic control underlying complex traits and pathological conditions; they will continue to do so with global collaborations and advances in genotyping technology. Finally, the increasing number of trans-ancestry GWAS and initiatives to include ancestry diversity in genomics will boost the power for discoveries, making them also applicable to non-European populations.

Multi-omic molecular profiling and network biology for precision anaesthesiology: a narrative review

Technological advancement, data democratisation, and decreasing costs have led to a revolution in molecular biology in which the entire set of DNA, RNA, proteins, and various other molecules - the 'multi-omic' profile - can be measured in humans. Sequencing 1 million bases of human DNA now costs US$0.01, and emerging technologies soon promise to reduce the cost of sequencing the whole genome to US$100. These trends have made it feasible to sample the multi-omic profile of millions of people, much of which is publicly available for medical research. Can anaesthesiologists use these data to improve patient care? This narrative review brings together a rapidly growing literature in multi-omic profiling across numerous fields that points to the future of precision anaesthesiology. Here, we discuss how DNA, RNA, proteins, and other molecules interact in molecular networks that can be used for preoperative risk stratification, intraoperative optimisation, and postoperative monitoring. This literature provides evidence for four fundamental insights: (1) Clinically similar patients have different molecular profiles and, as a consequence, different outcomes. (2) Vast, publicly available, and rapidly growing molecular datasets have been generated in chronic disease patients and can be repurposed to estimate perioperative risk. (3) Multi-omic networks are altered in the perioperative period and influence postoperative outcomes. (4) Multi-omic networks can serve as empirical, molecular measurements of a successful postoperative course. With this burgeoning universe of molecular data, the anaesthesiologist-of-the-future will tailor their clinical management to an individual's multi-omic profile to optimise postoperative outcomes and long-term health.

Atrioventricular node dysfunction in pressure overload-induced heart failure—Involvement of the immune system and transcriptomic remodelling

Heart failure is associated with atrioventricular (AV) node dysfunction, and AV node dysfunction in the setting of heart failure is associated with an increased risk of mortality and heart failure hospitalisation. This study aims to understand the causes of AV node dysfunction in heart failure by studying changes in the whole nodal transcriptome. The mouse transverse aortic constriction model of pressure overload-induced heart failure was studied; functional changes were assessed using electrocardiography and echocardiography and the transcriptome of the AV node was quantified using RNAseq. Heart failure was associated with a significant increase in the PR interval, indicating a slowing of AV node conduction and AV node dysfunction, and significant changes in 3,077 transcripts (5.6% of the transcriptome). Many systems were affected: transcripts supporting AV node conduction were downregulated and there were changes in transcripts identified by GWAS as determinants of the PR interval. In addition, there was evidence of remodelling of the sarcomere, a shift from fatty acid to glucose metabolism, remodelling of the extracellular matrix, and remodelling of the transcription and translation machinery. There was evidence of the causes of this widespread remodelling of the AV node: evidence of dysregulation of multiple intracellular signalling pathways, dysregulation of 109 protein kinases and 148 transcription factors, and an immune response with a proliferation of neutrophils, monocytes, macrophages and B lymphocytes and a dysregulation of 40 cytokines. In conclusion, inflammation and a widespread transcriptional remodelling of the AV node underlies AV node dysfunction in heart failure.

Fine mapping spatiotemporal mechanisms of genetic variants underlying cardiac traits and disease

The causal variants and genes underlying thousands of cardiac GWAS signals have yet to be identified. Here, we leverage spatiotemporal information on 966 RNA-seq cardiac samples and perform an expression quantitative trait locus (eQTL) analysis detecting eQTLs considering both eGenes and eIsoforms. We identify 2,578 eQTLs associated with a specific developmental stage-, tissue- and/or cell type. Colocalization between eQTL and GWAS signals of five cardiac traits identified variants with high posterior probabilities for being causal in 210 GWAS loci. Pulse pressure GWAS loci are enriched for colocalization with fetal- and smooth muscle- eQTLs; pulse rate with adult- and cardiac muscle- eQTLs; and atrial fibrillation with cardiac muscle- eQTLs. Fine mapping identifies 79 credible sets with five or fewer SNPs, of which 15 were associated with spatiotemporal eQTLs. Our study shows that many cardiac GWAS variants impact traits and disease in a developmental stage-, tissue- and/or cell type-specific fashion.

Clinical characteristics and electrophysiologic properties of SCN5A variants in fever-induced Brugada syndrome

BACKGROUND

Brugada syndrome (BrS) is a severe inherited arrhythmia syndrome that can be unmasked by fever.

METHODS

A multicentre clinical analysis was performed in 261 patients diagnosed with fever-induced BrS, including 198 (75.9%) and 27 (10.3%) patients who received next-generation genetic sequencing and epicardial arrhythmogenic substrate (AS) mapping, respectively.

FINDINGS

In fever-induced BrS patients, pathogenic or likely pathogenic (P/LP) SCN5A variant carriers developed fever-induced BrS at a younger age, and more often in females and those of Caucasian descent. They exhibited significant electrophysical abnormalities, including a larger epicardial AS area, and more prolonged abnormal epicardial electrograms. During a median follow-up of 50.5 months (quartiles 32.5-81.5 months) after the diagnosis, major cardiac events (MCE) occurred in 27 (14.4%) patients. Patients with P/LP SCN5A variants had a higher ratio of MCE compared with the rest. Additionally, history of syncope, QRS duration, and Tpe interval could also predict an increased risk for future MCE according to univariate analysis. Multivariate analysis indicated that only P/LP SCN5A variants were independent significant predictors of MCE. Computational structural modelling showed that most variants are destabilizing, suggesting that Nav1.5 structure destabilization caused by SCN5A missense variants may contribute to fever-induced BrS.

INTERPRETATION

In our cohort, P/LP SCN5A variant carriers with fever-induced BrS are more prevalent among patients of Caucasian descent, females, and younger patients. These patients exhibit aggressive electrophysiological abnormalities and worse outcome, which warrants closer monitoring and more urgent management of fever.

FUNDING

The current work was supported by the National Natural Science Foundation Project of China (Nos. 82270332 & 81670304), The Fundamental Research Funds for the Central Universities of China - Independent Research Project of Wuhan University (No. 2042022kf1217) from China; the National Institutes of Health of USA [NIH R56 (HL47678), NIH R01 (HL138103), and NIH R01 (HL152201)], the W. W. Smith Charitable Trust and the Wistar and Martha Morris Fund, Sharpe-Strumia Research Foundation, the American Heart Association Postdoctoral Fellowship (20POST35220002) from United States; the Netherlands CardioVascular Research Initiative: the Dutch Heart Foundation, Dutch Federation of University Medical Centers, the Netherlands Organization for Health Research and Development, and the Royal Netherlands Academy of Sciences (PREDICT2) from the Netherlands.

Genetic architecture of natural variation of cardiac performance from flies to humans

Deciphering the genetic architecture of human cardiac disorders is of fundamental importance but their underlying complexity is a major hurdle. We investigated the natural variation of cardiac performance in the sequenced inbred lines of the Drosophila Genetic Reference Panel (DGRP). Genome-wide associations studies (GWAS) identified genetic networks associated with natural variation of cardiac traits which were used to gain insights as to the molecular and cellular processes affected. Non-coding variants that we identified were used to map potential regulatory non-coding regions, which in turn were employed to predict transcription factors (TFs) binding sites. Cognate TFs, many of which themselves bear polymorphisms associated with variations of cardiac performance, were also validated by heart-specific knockdown. Additionally, we showed that the natural variations associated with variability in cardiac performance affect a set of genes overlapping those associated with average traits but through different variants in the same genes. Furthermore, we showed that phenotypic variability was also associated with natural variation of gene regulatory networks. More importantly, we documented correlations between genes associated with cardiac phenotypes in both flies and humans, which supports a conserved genetic architecture regulating adult cardiac function from arthropods to mammals. Specifically, roles for PAX9 and EGR2 in the regulation of the cardiac rhythm were established in both models, illustrating that the characteristics of natural variations in cardiac function identified in Drosophila can accelerate discovery in humans.

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

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

Partial and complete loss of myosin binding protein H-like cause cardiac conduction defects

A premature truncation of MYBPHL in humans and a loss of Mybphl in mice is associated with dilated cardiomyopathy, atrial and ventricular arrhythmias, and atrial enlargement. MYBPHL encodes myosin binding protein H-like (MyBP-HL). Prior work in mice indirectly identified Mybphl expression in the atria and in small puncta throughout the ventricle. Because of its genetic association with human and mouse cardiac conduction system disease, we evaluated the anatomical localization of MyBP-HL and the consequences of loss of MyBP-HL on conduction system function. Immunofluorescence microscopy of normal adult mouse ventricles identified MyBP-HL-positive ventricular cardiomyocytes that co-localized with the ventricular conduction system marker contactin-2 near the atrioventricular node and in a subset of Purkinje fibers. Mybphl heterozygous ventricles had a marked reduction of MyBP-HL-positive cells compared to controls. Lightsheet microscopy of normal perinatal day 5 mouse hearts showed enrichment of MyBP-HL-positive cells within and immediately adjacent to the contactin-2-positive ventricular conduction system, but this association was not apparent in Mybphl heterozygous hearts. Surface telemetry of Mybphl-null mice revealed atrioventricular block and atrial bigeminy, while intracardiac pacing revealed a shorter atrial relative refractory period and atrial tachycardia. Calcium transient analysis of isolated Mybphl-null atrial cardiomyocytes demonstrated an increased heterogeneity of calcium release and faster rates of calcium release compared to wild type controls. Super-resolution microscopy of Mybphl heterozygous and homozygous null atrial cardiomyocytes showed ryanodine receptor disorganization compared to wild type controls. Abnormal calcium release, shorter atrial refractory period, and atrial dilation seen in Mybphl null, but not wild type control hearts, agree with the observed atrial arrhythmias, bigeminy, and atrial tachycardia, whereas the proximity of MyBP-HL-positive cells with the ventricular conduction system provides insight into how a predominantly atrial expressed gene contributes to ventricular arrhythmias and ventricular dysfunction.

Knockout of Sorbin And SH3 Domain Containing 2 (Sorbs2) in Cardiomyocytes Leads to Dilated Cardiomyopathy in Mice

Background

Sorbin and SH3 domain containing 2 (Sorbs2) protein is a cytoskeletal adaptor with an emerging role in cardiac biology and disease; yet, its potential relevance to adult‐onset cardiomyopathies remains underexplored. Sorbs2 global knockout mice display lethal arrhythmogenic cardiomyopathy; however, the causative mechanisms remain unclear. Herein, we examine Sorbs2 dysregulation in heart failure, characterize novel Sorbs2 cardiomyocyte‐specific knockout mice (Sorbs2‐cKO), and explore associations between Sorbs2 genetic variations and human cardiovascular disease.

Methods and Results

Bioinformatic analyses show myocardial Sorbs2 mRNA is consistently upregulated in humans with adult‐onset cardiomyopathies and in heart failure models. We generated Sorbs2‐cKO mice and report that they develop progressive systolic dysfunction and enlarged cardiac chambers, and they die with congestive heart failure at about 1 year old. After 3 months, Sorbs2‐cKO mice begin to show atrial enlargement and P‐wave anomalies, without dysregulation of action potential–associated ion channel and gap junction protein expressions. After 6 months, Sorbs2‐cKO mice exhibit impaired contractility in dobutamine‐treated hearts and skinned myofibers, without dysregulation of contractile protein expressions. From our comprehensive survey of potential mechanisms, we found that within 4 months, Sorbs2‐cKO hearts have defective microtubule polymerization and compensatory upregulation of structural cytoskeletal and adapter proteins, suggesting that this early intracellular structural remodeling is responsible for contractile dysfunction. Finally, we identified genetic variants that associate with decreased Sorbs2 expression and human cardiac phenotypes, including conduction abnormalities, atrial enlargement, and dilated cardiomyopathy, consistent with Sorbs2‐cKO mice phenotypes.

Conclusions

Our studies show that Sorbs2 is essential for maintaining structural integrity in cardiomyocytes, likely through strengthening the interactions between microtubules and other cytoskeletal proteins at cross‐link sites.

Familial risk of atrioventricular block in first-degree relatives

Objective

Rare cases of genetically inherited atrioventricular block (AVB) have been reported; however, the heredity of AVB remains unknown. We aimed to assess the heredity of AVB.

Design, setting and participants

Using data from the Danish Civil Registration Registry, we established a nationwide cohort of individuals with parental links. Data were merged with information from the Danish Pacemaker and Implantable Cardioverter Defibrillator Registry, containing information on all pacemaker implantations performed in Denmark during the study period, to identify patients who received a first-time pacemaker because of AVB.

Results

A total of 4 648 204 individuals had parental links and a total of 26 880 consecutive patients received a first-time pacemaker due to AVB. Overall, the adjusted rate ratio (RR) of pacemaker implantation due to AVB was 2.1 (95% CI 1.8 to 2.5) if a father, mother or sibling had AVB compared with the risk in the general population. The adjusted RR was 2.2 (1.7–2.9) for offspring of mothers with AVB, 1.9 (1.5–2.4) for offspring of fathers with AVB and 3.5 (2.3–5.4) for siblings to a patient with AVB. The risk increased inversely proportionally with the age of the index case at the time of pacemaker implantation. The corresponding adjusted RRs were 15.8 (4.8–52.3) and 10.0 (3.3–30.4) if a mother or father, respectively, had a pacemaker implantation before 50 years.

Conclusion and relevance

First-degree relatives to a patient with AVB carry an increased risk of AVB with the risk being strongly inversely associated with the age of the index case at pacemaker implantation. These findings indicate a genetic component in the development of AVB in families with an early-onset disease.

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 NaV1.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 NaV1.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.

Common variants in SCN10A gene associated with Brugada syndrome

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. 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. In conclusion, 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.

Inherited and Acquired Rhythm Disturbances in Sick Sinus Syndrome, Brugada Syndrome, and Atrial Fibrillation: Lessons from Preclinical Modeling

Rhythm disturbances are life-threatening cardiovascular diseases, accounting for many deaths annually worldwide. Abnormal electrical activity might arise in a structurally normal heart in response to specific triggers or as a consequence of cardiac tissue alterations, in both cases with catastrophic consequences on heart global functioning. Preclinical modeling by recapitulating human pathophysiology of rhythm disturbances is fundamental to increase the comprehension of these diseases and propose effective strategies for their prevention, diagnosis, and clinical management. In silico, in vivo, and in vitro models found variable application to dissect many congenital and acquired rhythm disturbances. In the copious list of rhythm disturbances, diseases of the conduction system, as sick sinus syndrome, Brugada syndrome, and atrial fibrillation, have found extensive preclinical modeling. In addition, the electrical remodeling as a result of other cardiovascular diseases has also been investigated in models of hypertrophic cardiomyopathy, cardiac fibrosis, as well as arrhythmias induced by other non-cardiac pathologies, stress, and drug cardiotoxicity. This review aims to offer a critical overview on the effective ability of in silico bioinformatic tools, in vivo animal studies, in vitro models to provide insights on human heart rhythm pathophysiology in case of sick sinus syndrome, Brugada syndrome, and atrial fibrillation and advance their safe and successful translation into the cardiology arena.

Mapping gene and gene pathways associated with coronary artery disease: a CARDIoGRAM exome and multi-ancestry UK biobank analysis

Coronary artery disease (CAD) genome-wide association studies typically focus on single nucleotide variants (SNVs), and many potentially associated SNVs fail to reach the GWAS significance threshold. We performed gene and pathway-based association (GBA) tests on publicly available Coronary ARtery DIsease Genome wide Replication and Meta-analysis consortium Exome (n = 120,575) and multi ancestry pan UK Biobank study (n = 442,574) summary data using versatile gene-based association study (VEGAS2) and Multi-marker analysis of genomic annotation (MAGMA) to identify novel genes and pathways associated with CAD. We included only exonic SNVs and excluded regulatory regions. VEGAS2 and MAGMA ranked genes and pathways based on aggregated SNV test statistics. We used Bonferroni corrected gene and pathway significance threshold at 3.0 × 10^-6 and 1.0 × 10^-5, respectively. We also report the top one percent of ranked genes and pathways. We identified 17 top enriched genes with four genes (PCSK9, FAM177, LPL, ARGEF26), reaching statistical significance (p ≤ 3.0 × 10^-6) using both GBA tests in two GWAS studies. In addition, our analyses identified ten genes (DUSP13, KCNJ11, CD300LF/RAB37, SLCO1B1, LRRFIP1, QSER1, UBR2, MOB3C, MST1R, and ABCC8) with previously unreported associations with CAD, although none of the single SNV associations within the genes were genome-wide significant. Among the top 1% non-lipid pathways, we detected pathways regulating coagulation, inflammation, neuronal aging, and wound healing.

Rare Coding Variants Associated With Electrocardiographic Intervals Identify Monogenic Arrhythmia Susceptibility Genes: A Multi-Ancestry Analysis

BACKGROUND

Alterations in electrocardiographic (ECG) intervals are well-known markers for arrhythmia and sudden cardiac death (SCD) risk. While the genetics of arrhythmia syndromes have been studied, relations between electrocardiographic intervals and rare genetic variation at a population level are poorly understood.

METHODS

Using a discovery sample of 29 000 individuals with whole-genome sequencing from Trans-Omics in Precision Medicine and replication in nearly 100 000 with whole-exome sequencing from the UK Biobank and MyCode, we examined associations between low-frequency and rare coding variants with 5 routinely measured electrocardiographic traits (RR, P-wave, PR, and QRS intervals and corrected QT interval).

RESULTS

We found that rare variants associated with population-based electrocardiographic intervals identify established monogenic SCD genes (KCNQ1, KCNH2, and SCN5A), a controversial monogenic SCD gene (KCNE1), and novel genes (PAM and MFGE8) involved in cardiac conduction. Loss-of-function and pathogenic SCN5A variants, carried by 0.1% of individuals, were associated with a nearly 6-fold increased odds of the first-degree atrioventricular block (P=8.4×10^-5). Similar variants in KCNQ1 and KCNH2 (0.2% of individuals) were associated with a 23-fold increased odds of marked corrected QT interval prolongation (P=4×10^-25), a marker of SCD risk. Incomplete penetrance of such deleterious variation was common as over 70% of carriers had normal electrocardiographic intervals.

CONCLUSIONS

Our findings indicate that large-scale high-depth sequence data and electrocardiographic analysis identifies monogenic arrhythmia susceptibility genes and rare variants with large effects. Known pathogenic variation in conventional arrhythmia and SCD genes exhibited incomplete penetrance and accounted for only a small fraction of marked electrocardiographic interval prolongation.

A Mendelian randomization analysis of the relationship between cardioembolic risk factors and ischemic stroke

Observational studies have shown that several risk factors are associated with cardioembolic stroke. However, whether such associations reflect causality remains unknown. We aimed to determine whether established and provisional cardioembolic risk factors are causally associated with cardioembolic stroke. Genetic instruments for atrial fibrillation (AF), myocardial infarction (MI), electrocardiogram (ECG) indices and N-terminal pro-brain natriuretic peptide (NT-pro BNP) were obtained from large genetic consortiums. Summarized data of ischemic stroke and its subtypes were extracted from the MEGASTROKE consortium. Causal estimates were calculated by applying inverse-variance weighted analysis, weighted median analysis, simple median analysis and Mendelian randomization (MR)-Egger regression. Genetically predicted AF was significantly associated with higher odds of ischemic stroke (odds ratio (OR): 1.20, 95% confidence intervals (CI): 1.16-1.24, P = 6.53 × 10^-30) and cardioembolic stroke (OR: 1.95, 95% CI: 1.85-2.06, P = 8.81 × 10^-125). Suggestive associations were found between genetically determined resting heart rate and higher odds of ischemic stroke (OR: 1.01, 95% CI: 1.00-1.02, P = 0.005), large-artery atherosclerotic stroke (OR: 1.02, 95% CI: 1.00-1.04, P = 0.026) and cardioembolic stroke (OR: 1.02, 95% CI: 1.00-1.04, P = 0.028). There was no causal association of P-wave terminal force in the precordial lead V1 (PTFVI), P-wave duration (PWD), NT-pro BNP or PR interval with ischemic stroke or any subtype.

Genome-wide association study of neck circumference identifies sex-specific loci independent of generalized adiposity

BACKGROUND/OBJECTIVES

Neck circumference, an index of upper airway fat, has been suggested to be an important measure of body-fat distribution with unique associations with health outcomes such as obstructive sleep apnea and metabolic disease. This study aims to study the genetic bases of neck circumference.

METHODS

We conducted a multi-ethnic genome-wide association study of neck circumference, adjusted and unadjusted for BMI, in up to 15,090 European Ancestry (EA) and African American (AA) individuals. Because sexually dimorphic associations have been observed for anthropometric traits, we conducted both sex-combined and sex-specific analysis.

RESULTS

We identified rs227724 near the Noggin (NOG) gene as a possible quantitative locus for neck circumference in men (N = 8831, P = 1.74 × 10-9) but not in women (P = 0.08). The association was replicated in men (N = 1554, P = 0.045) in an independent dataset. This locus was previously reported to be associated with human height and with self-reported snoring. We also identified rs13087058 on chromosome 3 as a suggestive locus in sex-combined analysis (N = 15090, P = 2.94 × 10-7; replication P =0.049). This locus was also associated with electrocardiogram-assessed PR interval and is a cis-expression quantitative locus for the PDZ Domain-containing ring finger 2 (PDZRN3) gene. Both NOG and PDZRN3 interact with members of transforming growth factor-beta superfamily signaling proteins.

CONCLUSIONS

Our study suggests that neck circumference may have unique genetic basis independent of BMI.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

Genome-first approach to rare EYA4 variants and cardio-auditory phenotypes in adults

While newborns and children with hearing loss are routinely offered genetic testing, adults are rarely clinically tested for a genetic etiology. One clinically actionable result from genetic testing in children is the discovery of variants in syndromic hearing loss genes. EYA4 is a known hearing loss gene which is also involved in important pathways in cardiac tissue. The pleiotropic effects of rare EYA4 variants are poorly understood and their prevalence in a large cohort has not been previously reported. We investigated cardio-auditory phenotypes in 11,451 individuals in a large biobank using a rare variant, genome-first approach to EYA4. We filtered 256 EYA4 variants carried by 6737 participants to 26 rare and predicted deleterious variants carried by 42 heterozygotes. We aggregated predicted deleterious EYA4 gene variants into a combined variable (i.e. "gene burden") and performed association studies across phenotypes compared to wildtype controls. We validated findings with replication in three independent cohorts and human tissue expression data. EYA4 gene burden was significantly associated with audiometric-proven HL (p = [Formula: see text], Mobitz Type II AV block (p = [Formula: see text]) and the syndromic presentation of HL and primary cardiomyopathy (p = 0.0194). Analyses on audiogram, echocardiogram, and electrocardiogram data validated these associations. Prior reports have focused on identifying variants in families with severe or syndromic phenotypes. In contrast, we found, using a genotype-first approach, that gene burden in EYA4 is associated with more subtle cardio-auditory phenotypes in an adult medical biobank population, including cardiac conduction disorders which have not been previously reported. We show the value of using a focused approach to uncover human disease related to pleiotropic gene variants and suggest a role for genetic testing in adults presenting with hearing loss.

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.

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.

Network-driven discovery yields new insight into Shox2-dependent cardiac rhythm control

The homeodomain transcription factor SHOX2 is involved in the development and function of the heart's primary pacemaker, the sinoatrial node (SAN), and has been associated with cardiac conduction-related diseases such as atrial fibrillation and sinus node dysfunction. To shed light on Shox2-dependent genetic processes involved in these diseases, we established a murine embryonic stem cell (ESC) cardiac differentiation model to investigate Shox2 pathways in SAN-like cardiomyocytes. Differential RNA-seq-based expression profiling of Shox2^+/+ and Shox2^-/- ESCs revealed 94 dysregulated transcripts in Shox2^-/- ESC-derived SAN-like cells. Of these, 15 putative Shox2 target genes were selected for further validation based on comparative expression analysis with SAN- and right atria-enriched genes. Network-based analyses, integrating data from the Mouse Organogenesis Cell Atlas and the Ingenuity pathways, as well as validation in mouse and zebrafish models confirmed a regulatory role for the novel identified Shox2 target genes including Cav1, Fkbp10, Igfbp5, Mcf2l and Nr2f2. Our results indicate that genetic networks involving SHOX2 may contribute to conduction traits through the regulation of these genes.

A genome-wide association and polygenic risk score study on abnormal electrocardiogram in a Chinese population

Electrocardiography is a common and widely-performed medical examination based on the measurement and evaluation of electrocardiogram (ECG) to assess the up-to-date cardiac rhythms and thus suggest the health conditions of cardiovascular system and on a larger level the individual’s wellness. Abnormal ECG assessment from the detection of abnormal heart rhythms may have clinical implications including blood clots in formation, ongoing heart attack, coronary artery blockage, etc. Past genetic-phenotypic research focused primarily on the physical parameters of ECG but not the medical evaluation. To unbiasedly uncover the underlying links of genetic variants with normal vs. abnormal ECG assessment, a genome-wide association study (GWAS) is carried out in a 1006-participant cohort of Chinese population effectively genotyped for 243487 single nucleotide polymorphisms (SNPs). Both age and sex are influential factors, and six novel SNPs are identified for potential association with abnormal ECG. With the selected SNPs, a polygenic risk score (PRS) differentiates the case–control subgroups, and correlates well with increased risk of abnormal ECG. The findings are reproduced in an independent validation cohort. The derived PRS may function as a potential biomarker for prospectively screening the high-risk subgroup of heart issues in the Chinese population.

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

AIMS

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

METHODS AND RESULTS

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

CONCLUSIONS

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

GATA6 is a regulator of sinus node development and heart rhythm

The sinus node (SAN) is the primary pacemaker of the human heart, and abnormalities in its structure or function cause sick sinus syndrome, the most common reason for electronic pacemaker implantation. Here we report that transcription factor GATA6, whose mutations in humans are linked to arrhythmia, is highly expressed in the SAN and its haploinsufficiency in mice results in hypoplastic SANs and rhythm abnormalities. Cell-specific deletion reveals a requirement for GATA6 in various SAN lineages. Mechanistically, GATA6 directly activates key regulators of the SAN genetic program in conduction and nonconduction cells, such as TBX3 and EDN1, respectively. The data identify GATA6 as an important regulator of the SAN and provide a molecular basis for understanding the conduction abnormalities associated with GATA6 mutations in humans. They also suggest that GATA6 may be a potential modifier of the cardiac pacemaker.

Genome-Wide Analysis Identifies an Essential Human TBX3 Pacemaker Enhancer

RATIONALE

The development and function of the pacemaker cardiomyocytes of the sinoatrial node (SAN), the leading pacemaker of the heart, are tightly controlled by a conserved network of transcription factors, including TBX3 (T-box transcription factor 3), ISL1 (ISL LIM homeobox 1), and SHOX2 (short stature homeobox 2). Yet, the regulatory DNA elements (REs) controlling target gene expression in the SAN pacemaker cells have remained undefined.

OBJECTIVE

Identification of the regulatory landscape of human SAN-like pacemaker cells and functional assessment of SAN-specific REs potentially involved in pacemaker cell gene regulation.

METHODS AND RESULTS

We performed Assay for Transposase-Accessible Chromatin using sequencing on human pluripotent stem cell-derived SAN-like pacemaker cells and ventricle-like cells and identified thousands of putative REs specific for either human cell type. We validated pacemaker cell-specific elements in the SHOX2 and TBX3 loci. CRISPR-mediated homozygous deletion of the mouse ortholog of a noncoding region with candidate pacemaker-specific REs in the SHOX2 locus resulted in selective loss of Shox2 expression from the developing SAN and embryonic lethality. Putative pacemaker-specific REs were identified up to 1 Mbp upstream of TBX3 in a region close to MED13L harboring variants associated with heart rate recovery after exercise. The orthologous region was deleted in mice, which resulted in selective loss of expression of Tbx3 from the SAN and (cardiac) ganglia and in neonatal lethality. Expression of Tbx3 was maintained in other tissues including the atrioventricular conduction system, lungs, and liver. Heterozygous adult mice showed increased SAN recovery times after pacing. The human REs harboring the associated variants robustly drove expression in the SAN of transgenic mouse embryos.

CONCLUSIONS

We provided a genome-wide collection of candidate human pacemaker-specific REs, including the loci of SHOX2, TBX3, and ISL1, and identified a link between human genetic variants influencing heart rate recovery after exercise and a variant RE with highly conserved function, driving SAN expression of TBX3.

Transcriptomic Heterogeneity of Alzheimer's Disease Associated with Lipid Genetic Risk

Alzheimer's disease (AD) is a multifactorial disease that affects more than 5 million Americans. Multiple pathways might be involved in the AD pathogenesis. The implication of lipid genetic susceptibility on brain gene expression is yet to be investigated. The current study included 192 brain samples from AD patients who were enrolled in the ROSMAP study. The samples were genotyped and imputed to the HRC Reference Panel. Lipid polygenetic risk score was constructed from the weighted sum of genetic variants associated with low-density lipoprotein cholesterol (LDL-C). The gene expression was profiled by RNA sequencing, and the  association of gene expression with lipid polygenetic risk scores was tested by linear regression models adjusted for age, sex and APOE e4 alleles. Three genes were found to associate with lipid polygenetic risk scores, including HMCN2 (P = 3.6 × 10-7), PDLIM5 (P = 1.2 × 10-6), and FHL5 (P = 2.0 × 10-6). Network analysis revealed multiple related pathways, including dopaminergic synapse (P = 4.5 × 10-5), circadian entrainment (P = 1.1 × 10-4), and cholinergic synapse (P = 2.3 × 10-4). Our study underscores the importance of lipid regulation and metabolism to AD heterogeneity.

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.

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.

Genetic Determinants of Electrocardiographic P-Wave Duration and Relation to Atrial Fibrillation

BACKGROUND

The P-wave duration (PWD) is an electrocardiographic measurement that represents cardiac conduction in the atria. Shortened or prolonged PWD is associated with atrial fibrillation (AF). We used exome-chip data to examine the associations between common and rare variants with PWD.

METHODS

Fifteen studies comprising 64 440 individuals (56 943 European, 5681 African, 1186 Hispanic, 630 Asian) and ≈230 000 variants were used to examine associations with maximum PWD across the 12-lead ECG. Meta-analyses summarized association results for common variants; gene-based burden and sequence kernel association tests examined low-frequency variant-PWD associations. Additionally, we examined the associations between PWD loci and AF using previous AF genome-wide association studies.

RESULTS

We identified 21 common and low-frequency genetic loci (14 novel) associated with maximum PWD, including several AF loci (TTN, CAND2, SCN10A, PITX2, CAV1, SYNPO2L, SOX5, TBX5, MYH6, RPL3L). The top variants at known sarcomere genes (TTN, MYH6) were associated with longer PWD and increased AF risk. However, top variants at other loci (eg, PITX2 and SCN10A) were associated with longer PWD but lower AF risk.

CONCLUSIONS

Our results highlight multiple novel genetic loci associated with PWD, and underscore the shared mechanisms of atrial conduction and AF. Prolonged PWD may be an endophenotype for several different genetic mechanisms of AF.

The Genetic Makeup of the Electrocardiogram

The electrocardiogram (ECG) is one of the most useful non-invasive diagnostic tests for a wide array of cardiac disorders. Traditional approaches to analyzing ECGs focus on individual segments. Here, we performed comprehensive deep phenotyping of 77,190 ECGs in the UK Biobank across the complete cycle of cardiac conduction, resulting in 500 spatial-temporal datapoints, across 10 million genetic variants. In addition to characterizing polygenic risk scores for the traditional ECG segments, we identified over 300 genetic loci that are statistically associated with the high-dimensional representation of the ECG. We established the genetic ECG signature for dilated cardiomyopathy, associated the BAG3, HSPB7/CLCNKA, PRKCA, TMEM43, and OBSCN loci with disease risk and confirmed this association in an independent cohort. In total, our work demonstrates that a high-dimensional analysis of the entire ECG provides unique opportunities for studying cardiac biology and disease and furthering drug development. A record of this paper's transparent peer review process is included in the Supplemental Information.

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

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

Transcriptional and Cellular Diversity of the Human Heart

BACKGROUND

The human heart requires a complex ensemble of specialized cell types to perform its essential function. A greater knowledge of the intricate cellular milieu of the heart is critical to increase our understanding of cardiac homeostasis and pathology. As recent advances in low-input RNA sequencing have allowed definitions of cellular transcriptomes at single-cell resolution at scale, we have applied these approaches to assess the cellular and transcriptional diversity of the nonfailing human heart.

METHODS

Microfluidic encapsulation and barcoding was used to perform single nuclear RNA sequencing with samples from 7 human donors, selected for their absence of overt cardiac disease. Individual nuclear transcriptomes were then clustered based on transcriptional profiles of highly variable genes. These clusters were used as the basis for between-chamber and between-sex differential gene expression analyses and intersection with genetic and pharmacologic data.

RESULTS

We sequenced the transcriptomes of 287 269 single cardiac nuclei, revealing 9 major cell types and 20 subclusters of cell types within the human heart. Cellular subclasses include 2 distinct groups of resident macrophages, 4 endothelial subtypes, and 2 fibroblast subsets. Comparisons of cellular transcriptomes by cardiac chamber or sex reveal diversity not only in cardiomyocyte transcriptional programs but also in subtypes involved in extracellular matrix remodeling and vascularization. Using genetic association data, we identified strong enrichment for the role of cell subtypes in cardiac traits and diseases. Intersection of our data set with genes on cardiac clinical testing panels and the druggable genome reveals striking patterns of cellular specificity.

CONCLUSIONS

Using large-scale single nuclei RNA sequencing, we defined the transcriptional and cellular diversity in the normal human heart. Our identification of discrete cell subtypes and differentially expressed genes within the heart will ultimately facilitate the development of new therapeutics for cardiovascular diseases.

Multi-Ethnic Genome-Wide Association Study of Decomposed Cardioelectric Phenotypes Illustrates Strategies to Identify and Characterize Evidence of Shared Genetic Effects for Complex Traits

BACKGROUND

We examined how expanding electrocardiographic trait genome-wide association studies to include ancestrally diverse populations, prioritize more precise phenotypic measures, and evaluate evidence for shared genetic effects enabled the detection and characterization of loci.

METHODS

We decomposed 10 seconds, 12-lead electrocardiograms from 34 668 multi-ethnic participants (15% Black; 30% Hispanic/Latino) into 6 contiguous, physiologically distinct (P wave, PR segment, QRS interval, ST segment, T wave, and TP segment) and 2 composite, conventional (PR interval and QT interval) interval scale traits and conducted multivariable-adjusted, trait-specific univariate genome-wide association studies using 1000-G imputed single-nucleotide polymorphisms. Evidence of shared genetic effects was evaluated by aggregating meta-analyzed univariate results across the 6 continuous electrocardiographic traits using the combined phenotype adaptive sum of powered scores test.

RESULTS

We identified 6 novels (CD36, PITX2, EMB, ZNF592, YPEL2, and BC043580) and 87 known loci (adaptive sum of powered score test P<5×10-9). Lead single-nucleotide polymorphism rs3211938 at CD36 was common in Blacks (minor allele frequency=10%), near monomorphic in European Americans, and had effects on the QT interval and TP segment that ranked among the largest reported to date for common variants. The other 5 novel loci were observed when evaluating the contiguous but not the composite electrocardiographic traits. Combined phenotype testing did not identify novel electrocardiographic loci unapparent using traditional univariate approaches, although this approach did assist with the characterization of known loci.

CONCLUSIONS

Despite including one-third as many participants as published electrocardiographic trait genome-wide association studies, our study identified 6 novel loci, emphasizing the importance of ancestral diversity and phenotype resolution in this era of ever-growing genome-wide association studies.

Transethnic Genome-Wide Association Study Provides Insights in the Genetic Architecture and Heritability of Long QT Syndrome

BACKGROUND

Long QT syndrome (LQTS) is a rare genetic disorder and a major preventable cause of sudden cardiac death in the young. A causal rare genetic variant with large effect size is identified in up to 80% of probands (genotype positive) and cascade family screening shows incomplete penetrance of genetic variants. Furthermore, a proportion of cases meeting diagnostic criteria for LQTS remain genetically elusive despite genetic testing of established genes (genotype negative). These observations raise the possibility that common genetic variants with small effect size contribute to the clinical picture of LQTS. This study aimed to characterize and quantify the contribution of common genetic variation to LQTS disease susceptibility.

METHODS

We conducted genome-wide association studies followed by transethnic meta-analysis in 1656 unrelated patients with LQTS of European or Japanese ancestry and 9890 controls to identify susceptibility single nucleotide polymorphisms. We estimated the common variant heritability of LQTS and tested the genetic correlation between LQTS susceptibility and other cardiac traits. Furthermore, we tested the aggregate effect of the 68 single nucleotide polymorphisms previously associated with the QT-interval in the general population using a polygenic risk score.

RESULTS

Genome-wide association analysis identified 3 loci associated with LQTS at genome-wide statistical significance (P<5×10-8) near NOS1AP, KCNQ1, and KLF12, and 1 missense variant in KCNE1(p.Asp85Asn) at the suggestive threshold (P<10-6). Heritability analyses showed that ≈15% of variance in overall LQTS susceptibility was attributable to common genetic variation (h2SNP 0.148; standard error 0.019). LQTS susceptibility showed a strong genome-wide genetic correlation with the QT-interval in the general population (rg=0.40; P=3.2×10-3). The polygenic risk score comprising common variants previously associated with the QT-interval in the general population was greater in LQTS cases compared with controls (P<10-13), and it is notable that, among patients with LQTS, this polygenic risk score was greater in patients who were genotype negative compared with those who were genotype positive (P<0.005).

CONCLUSIONS

This work establishes an important role for common genetic variation in susceptibility to LQTS. We demonstrate overlap between genetic control of the QT-interval in the general population and genetic factors contributing to LQTS susceptibility. Using polygenic risk score analyses aggregating common genetic variants that modulate the QT-interval in the general population, we provide evidence for a polygenic architecture in genotype negative LQTS.

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

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

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.