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

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.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.

Potential causal association of a prolonged PR interval and clinical recurrence of atrial fibrillation after catheter ablation: a Mendelian randomization analysis

A prolonged PR interval predicts atrial fibrillation (AF) recurrence after catheter ablation. We investigated the causal association between the PR interval and AF clinical recurrence by a Mendelian randomization. We prospectively included 1722 patients with AF (73.2% male, 58.6 ± 10.8 years old, 71.3% paroxysmal AF) who underwent catheter ablation into a genome-wide association study (GWAS). We searched for the genetic associations between the PR interval and AF recurrence by analyzing 44 single nucleotide polymorphisms (SNPs) already known to be associated with the PR interval, and investigated the Mendelian randomization. Based on the quartile analysis, the highest quartile of the PR interval was associated with an increased risk of AF recurrence compared with the lowest quartile (Hazard ratio (HR) = 1.91, 95% CI = 1.51-2.42, P = 8.41 × 10^-8) during 35.7 ± 28.5 months of follow-up. Among 44 SNPs known to be associated with the PR interval, two SNPs had significant associations with the PR interval (P < 0.001 for each SNP). CAV1 (HR = 1.15, 95% CI = 1.02-1.31, P = 0.024) was associated with clinical recurrence of AF. A Mendelian randomization analysis demonstrated a significant association with CAV1 (HR = 1.04, 95% CI = 1.01-1.07, P = 0.006). A prolonged PR interval was a risk factor for an AF recurrence, and the PR interval had a potentially causal association with an AF clinical recurrence after catheter ablation at the genetic level.

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

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

Allele-specific NKX2-5 binding underlies multiple genetic associations with human electrocardiographic traits

The cardiac transcription factor (TF) gene NKX2-5 has been associated with electrocardiographic (EKG) traits through genome-wide association studies (GWASs), but the extent to which differential binding of NKX2-5 at common regulatory variants contributes to these traits has not yet been studied. We analyzed transcriptomic and epigenomic data from induced pluripotent stem cell-derived cardiomyocytes from seven related individuals, and identified ~2,000 single-nucleotide variants associated with allele-specific effects (ASE-SNVs) on NKX2-5 binding. NKX2-5 ASE-SNVs were enriched for altered TF motifs, for heart-specific expression quantitative trait loci and for EKG GWAS signals. Using fine-mapping combined with epigenomic data from induced pluripotent stem cell-derived cardiomyocytes, we prioritized candidate causal variants for EKG traits, many of which were NKX2-5 ASE-SNVs. Experimentally characterizing two NKX2-5 ASE-SNVs (rs3807989 and rs590041) showed that they modulate the expression of target genes via differential protein binding in cardiac cells, indicating that they are functional variants underlying EKG GWAS signals. Our results show that differential NKX2-5 binding at numerous regulatory variants across the genome contributes to EKG phenotypes.

Variation in a Left Ventricle-Specific Hand1 Enhancer Impairs GATA Transcription Factor Binding and Disrupts Conduction System Development and Function

RATIONALE

The ventricular conduction system (VCS) rapidly propagates electrical impulses through the working myocardium of the ventricles to coordinate chamber contraction. GWAS (Genome-wide association studies) have associated nucleotide polymorphisms, most are located within regulatory intergenic or intronic sequences, with variation in VCS function. Two highly correlated polymorphisms (r²>0.99) associated with VCS functional variation (rs13165478 and rs13185595) occur 5' to the gene encoding the basic helix-loop-helix transcription factor HAND1 (heart- and neural crest derivatives-expressed protein 1).

OBJECTIVE

Here, we test the hypothesis that these polymorphisms influence HAND1 transcription thereby influencing VCS development and function.

METHODS AND RESULTS

We employed transgenic mouse models to identify an enhancer that is sufficient for left ventricle (LV) cis-regulatory activity. Two evolutionarily conserved GATA transcription factor cis-binding elements within this enhancer are bound by GATA4 and are necessary for cis-regulatory activity, as shown by in vitro DNA binding assays. CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9-mediated deletion of this enhancer dramatically reduces Hand1 expression solely within the LV but does not phenocopy previously published mouse models of cardiac Hand1 loss-of-function. Electrophysiological and morphological analyses reveals that mice homozygous for this deleted enhancer display a morphologically abnormal VCS and a conduction system phenotype consistent with right bundle branch block. Using 1000 Genomes Project data, we identify 3 additional single nucleotide polymorphisms (SNPs), located within the Hand1 LV enhancer, that compose a haplotype with rs13165478 and rs13185595. One of these SNPs, rs10054375, overlaps with a critical GATA cis-regulatory element within the Hand1 LV enhancer. This SNP, when tested in electrophoretic mobility shift assays, disrupts GATA4 DNA-binding. Modeling 2 of these SNPs in mice causes diminished Hand1 expression and mice present with abnormal VCS function.

CONCLUSIONS

Together, these findings reveal that SNP rs10054375, which is located within a necessary and sufficient LV-specific Hand1 enhancer, exhibits reduces GATA DNA-binding in electrophoretic mobility shift assay, and this enhancer in total, is required for VCS development and function in mice and perhaps humans.

Novel missense variant in TTN cosegregating with familial atrioventricular block

BACKGROUND

Cardiovascular diseases are the most common cause of death globally. In which atrioventricular block (AVB) is a common disorder with genetic causes, but the responsible genes have not been fully identified yet. To determine the underlying causative genes involved in cardiac AVB, here we report a three-generation Chinese family with severe autosomal dominant cardiac AVB that has been ruled out as being caused by known genes mutations.

METHODS

Whole-exome sequencing was performed in five affected family members across three generations, and co-segregation analysis was validated on other members of this family.

RESULTS

Whole-exome sequencing and subsequent co-segregation validation identified a novel germline heterozygous point missense mutation, c.49287C > A (p.N16429K), in the titin (TTN, NM_001267550.2) gene in all 5 affected family members but not in the unaffected family members, neither in the large population according to the Genome Aggregation Database (https://gnomad.broadinstitute.org/). The point mutation is predicted to be functionally deleterious by in-silico software tools. Our finding was further supported by the conservative analysis across species.

CONCLUSION

Based on this study, TTN was identified as a potential novel candidate gene for autosomal dominant AVB; this study expands the mutational spectrum of TTN gene and is the first to implicate TTN mutations as AVB disease causing in a Chinese pedigree.

A newly identified missense mutation in CLCA2 is associated with autosomal dominant cardiac conduction block

BACKGROUND

Progressive cardiac conduction defect (PCCD), also known as Lenegre-Lev disease, is one of the most common heart conduction abnormalities. Previous studies have screened for known mutation sites that cause heart block in a 68-person family with a history of PCCD, revealed no mutations.

OBJECTIVE

To screen pathogenic genes of the PCCD family and to study the function of the gene mutations related to heart block diseases.

METHODS

Whole exome sequencing (WES) was performed on two PCCD patients and one non-PCCD family member to find the related pathogenic gene. After family co-segregation and preliminary functional analysis, we identified the mutant gene CLCA2. To study the function of this gene, we constructed mutant-gene mice using CRISPR-Cas9 technology, and electrocardiogram monitoring was performed after genotype verification.

RESULTS

The CLCA2 c.G1725T mutation was identified and co-segregated with the phenotype. The analysis showed that the CLCA2 c.G1725T mutation is harmful and mainly affects protein glycosylation. Immunofluorescence staining revealed that CLCA2 was highly expressed in the sinoatrial node (SAN) tissues. Electrocardiogram monitoring of the mice revealed that CLCA2 point mutations induced mild conduction block and ectopic pacemakers.

CONCLUSION

Our findings indicate that a novel heterozygous missense mutation c.G1725T of the CLCA2 gene may be associated with heart block disease and the mutation in this gene may lead to sinus node lesions and conduction blocking.

GWAS of QRS duration identifies new loci specific to Hispanic/Latino populations

BACKGROUND

The electrocardiographically quantified QRS duration measures ventricular depolarization and conduction. QRS prolongation has been associated with poor heart failure prognosis and cardiovascular mortality, including sudden death. While previous genome-wide association studies (GWAS) have identified 32 QRS SNPs across 26 loci among European, African, and Asian-descent populations, the genetics of QRS among Hispanics/Latinos has not been previously explored.

METHODS

We performed a GWAS of QRS duration among Hispanic/Latino ancestry populations (n = 15,124) from four studies using 1000 Genomes imputed genotype data (adjusted for age, sex, global ancestry, clinical and study-specific covariates). Study-specific results were combined using fixed-effects, inverse variance-weighted meta-analysis.

RESULTS

We identified six loci associated with QRS (P<5x10-8), including two novel loci: MYOCD, a nuclear protein expressed in the heart, and SYT1, an integral membrane protein. The top SNP in the MYOCD locus, intronic SNP rs16946539, was found in Hispanics/Latinos with a minor allele frequency (MAF) of 0.04, but is monomorphic in European and African descent populations. The most significant QRS duration association was with intronic SNP rs3922344 (P = 1.19x10-24) in SCN5A/SCN10A. Three other previously identified loci, CDKN1A, VTI1A, and HAND1, also exceeded the GWAS significance threshold among Hispanics/Latinos. A total of 27 of 32 previously identified QRS duration SNPs were shown to generalize in Hispanics/Latinos.

CONCLUSIONS

Our QRS duration GWAS, the first in Hispanic/Latino populations, identified two new loci, underscoring the utility of extending large scale genomic studies to currently under-examined populations.

Genome-wide association meta-analysis of 30,000 samples identifies seven novel loci for quantitative ECG traits

Genome-wide association studies (GWAS) of quantitative electrocardiographic (ECG) traits in large consortia have identified more than 130 loci associated with QT interval, QRS duration, PR interval, and heart rate (RR interval). In the current study, we meta-analyzed genome-wide association results from 30,000 mostly Dutch samples on four ECG traits: PR interval, QRS duration, QT interval, and RR interval. SNP genotype data was imputed using the Genome of the Netherlands reference panel encompassing 19 million SNPs, including millions of rare SNPs (minor allele frequency < 5%). In addition to many known loci, we identified seven novel locus-trait associations: KCND3, NR3C1, and PLN for PR interval, KCNE1, SGIP1, and NFKB1 for QT interval, and ATP2A2 for QRS duration, of which six were successfully replicated. At these seven loci, we performed conditional analyses and annotated significant SNPs (in exons and regulatory regions), demonstrating involvement of cardiac-related pathways and regulation of nearby genes.

An expanded proteome of cardiac t-tubules

BACKGROUND

Transverse tubules (t-tubules) are important structural elements, derived from sarcolemma, found on all striated myocytes. These specialized organelles create a scaffold for many proteins crucial to the effective propagation of signal in cardiac excitation-contraction coupling. The full protein composition of this region is unknown.

METHODS

We characterized the t-tubule subproteome using 52,033 immunohistochemical images covering 13,203 proteins from the Human Protein Atlas (HPA) cardiac tissue microarrays. We used HPASubC, a suite of Python tools, to rapidly review and classify each image for a specific t-tubule staining pattern. The tools Gene Cards, String 11, and Gene Ontology Consortium as well as literature searches were used to understand pathways and relationships between the proteins.

RESULTS

There were 96 likely t-tubule proteins identified by HPASubC. Of these, 12 were matrisome proteins and 3 were mitochondrial proteins. A separate literature search identified 50 known t-tubule proteins. A comparison of the 2 lists revealed only 17 proteins in common, including 8 of the matrisome proteins. String11 revealed that 94 of 127 combined t-tubule proteins generated a single interconnected network.

CONCLUSION

Using HPASubC and the HPA, we identified 78 novel, putative t-tubule proteins and validated 17 within the literature. This expands and improves our knowledge of this important subcellular structure of the cardiac myocyte. This information can be used to identify new structural targets involved in excitation-contraction coupling that may be altered in disease.