Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways.

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries.

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

Rare variants in KDR, encoding VEGF Receptor 2, are associated with tetralogy of Fallot.

PURPOSE: Rare genetic variants in KDR, encoding the vascular endothelial growth factor receptor 2 (VEGFR2), have been reported in patients with tetralogy of Fallot (TOF). However, their role in disease causality and pathogenesis remains unclear. METHODS: We conducted exome sequencing in a familial case of TOF and large-scale genetic studies, including burden testing, in >1,500 patients with TOF. We studied gene-targeted mice and conducted cell-based assays to explore the role of KDR genetic variation in the etiology of TOF. RESULTS: Exome sequencing in a family with two siblings affected by TOF revealed biallelic missense variants in KDR. Studies in knock-in mice and in HEK 293T cells identified embryonic lethality for one variant when occurring in the homozygous state, and a significantly reduced VEGFR2 phosphorylation for both variants. Rare variant burden analysis conducted in a set of 1,569 patients of European descent with TOF identified a 46-fold enrichment of protein-truncating variants (PTVs) in TOF cases compared to controls (P = 7 × 10-11). CONCLUSION: Rare KDR variants, in particular PTVs, strongly associate with TOF, likely in the setting of different inheritance patterns. Supported by genetic and in vivo and in vitro functional analysis, we propose loss-of-function of VEGFR2 as one of the mechanisms involved in the pathogenesis of TOF.

Evidence-Based Assessment of Genes in Dilated Cardiomyopathy.

BACKGROUND: Each of the cardiomyopathies, classically categorized as hypertrophic cardiomyopathy, dilated cardiomyopathy (DCM), and arrhythmogenic right ventricular cardiomyopathy, has a signature genetic theme. Hypertrophic cardiomyopathy and arrhythmogenic right ventricular cardiomyopathy are largely understood as genetic diseases of sarcomere or desmosome proteins, respectively. In contrast, >250 genes spanning >10 gene ontologies have been implicated in DCM, representing a complex and diverse genetic architecture. To clarify this, a systematic curation of evidence to establish the relationship of genes with DCM was conducted. METHODS: An international panel with clinical and scientific expertise in DCM genetics evaluated evidence supporting monogenic relationships of genes with idiopathic DCM. The panel used the Clinical Genome Resource semiquantitative gene-disease clinical validity classification framework with modifications for DCM genetics to classify genes into categories on the basis of the strength of currently available evidence. Representation of DCM genes on clinically available genetic testing panels was evaluated. RESULTS: Fifty-one genes with human genetic evidence were curated. Twelve genes (23%) from 8 gene ontologies were classified as having definitive (BAG3, DES, FLNC, LMNA, MYH7, PLN, RBM20, SCN5A, TNNC1, TNNT2, TTN) or strong (DSP) evidence. Seven genes (14%; ACTC1, ACTN2, JPH2, NEXN, TNNI3, TPM1, VCL) including 2 additional ontologies were classified as moderate evidence; these genes are likely to emerge as strong or definitive with additional evidence. Of these 19 genes, 6 were similarly classified for hypertrophic cardiomyopathy and 3 for arrhythmogenic right ventricular cardiomyopathy. Of the remaining 32 genes (63%), 25 (49%) had limited evidence, 4 (8%) were disputed, 2 (4%) had no disease relationship, and 1 (2%) was supported by animal model data only. Of the 16 evaluated clinical genetic testing panels, most definitive genes were included, but panels also included numerous genes with minimal human evidence. CONCLUSIONS: In the curation of 51 genes, 19 had high evidence (12 definitive/strong, 7 moderate). It is notable that these 19 genes explain only a minority of cases, leaving the remainder of DCM genetic architecture incompletely addressed. Clinical genetic testing panels include most high-evidence genes; however, genes lacking robust evidence are also commonly included. We recommend that high-evidence DCM genes be used for clinical practice and that caution be exercised in the interpretation of variants in variable-evidence DCM genes.

Biallelic loss-of-function variants in PLD1 cause congenital right-sided cardiac valve defects and neonatal cardiomyopathy.

Congenital heart disease is the most common type of birth defect, accounting for one-third of all congenital anomalies. Using whole-exome sequencing of 2718 patients with congenital heart disease and a search in GeneMatcher, we identified 30 patients from 21 unrelated families of different ancestries with biallelic phospholipase D1 (PLD1) variants who presented predominantly with congenital cardiac valve defects. We also associated recessive PLD1 variants with isolated neonatal cardiomyopathy. Furthermore, we established that p.I668F is a founder variant among Ashkenazi Jews (allele frequency of ~2%) and describe the phenotypic spectrum of PLD1-associated congenital heart defects. PLD1 missense variants were overrepresented in regions of the protein critical for catalytic activity, and, correspondingly, we observed a strong reduction in enzymatic activity for most of the mutant proteins in an enzymatic assay. Finally, we demonstrate that PLD1 inhibition decreased endothelial-mesenchymal transition, an established pivotal early step in valvulogenesis. In conclusion, our study provides a more detailed understanding of disease mechanisms and phenotypic expression associated with PLD1 loss of function.

Shared genetic pathways contribute to risk of hypertrophic and dilated cardiomyopathies with opposite directions of effect.

The heart muscle diseases hypertrophic (HCM) and dilated (DCM) cardiomyopathies are leading causes of sudden death and heart failure in young, otherwise healthy, individuals. We conducted genome-wide association studies and multi-trait analyses in HCM (1,733 cases), DCM (5,521 cases) and nine left ventricular (LV) traits (19,260 UK Biobank participants with structurally normal hearts). We identified 16 loci associated with HCM, 13 with DCM and 23 with LV traits. We show strong genetic correlations between LV traits and cardiomyopathies, with opposing effects in HCM and DCM. Two-sample Mendelian randomization supports a causal association linking increased LV contractility with HCM risk. A polygenic risk score explains a significant portion of phenotypic variability in carriers of HCM-causing rare variants. Our findings thus provide evidence that polygenic risk score may account for variability in Mendelian diseases. More broadly, we provide insights into how genetic pathways may lead to distinct disorders through opposing genetic effects.

Enhancing rare variant interpretation in inherited arrhythmias through quantitative analysis of consortium disease cohorts and population controls.

PURPOSE: Stringent variant interpretation guidelines can lead to high rates of variants of uncertain significance (VUS) for genetically heterogeneous disease like long QT syndrome (LQTS) and Brugada syndrome (BrS). Quantitative and disease-specific customization of American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) guidelines can address this false negative rate. METHODS: We compared rare variant frequencies from 1847 LQTS (KCNQ1/KCNH2/SCN5A) and 3335 BrS (SCN5A) cases from the International LQTS/BrS Genetics Consortia to population-specific gnomAD data and developed disease-specific criteria for ACMG/AMP evidence classes-rarity (PM2/BS1 rules) and case enrichment of individual (PS4) and domain-specific (PM1) variants. RESULTS: Rare SCN5A variant prevalence differed between European (20.8%) and Japanese (8.9%) BrS patients (p = 5.7 × 10-18) and diagnosis with spontaneous (28.7%) versus induced (15.8%) Brugada type 1 electrocardiogram (ECG) (p = 1.3 × 10-13). Ion channel transmembrane regions and specific N-terminus (KCNH2) and C-terminus (KCNQ1/KCNH2) domains were characterized by high enrichment of case variants and >95% probability of pathogenicity. Applying the customized rules, 17.4% of European BrS and 74.8% of European LQTS cases had (likely) pathogenic variants, compared with estimated diagnostic yields (case excess over gnomAD) of 19.2%/82.1%, reducing VUS prevalence to close to background rare variant frequency. CONCLUSION: Large case-control data sets enable quantitative implementation of ACMG/AMP guidelines and increased sensitivity for inherited arrhythmia genetic testing.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals.

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

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.

Genome-wide association studies of cardiac electrical phenotypes.

The genetic basis of cardiac electrical phenotypes has in the last 25 years been the subject of intense investigation. While in the first years, such efforts were dominated by the study of familial arrhythmia syndromes, in recent years, large consortia of investigators have successfully pursued genome-wide association studies (GWAS) for the identification of single-nucleotide polymorphisms that govern inter-individual variability in electrocardiographic parameters in the general population. We here provide a review of GWAS conducted on cardiac electrical phenotypes in the last 14 years and discuss the implications of these discoveries for our understanding of the genetic basis of disease susceptibility and variability in disease severity. Furthermore, we review functional follow-up studies that have been conducted on GWAS loci associated with cardiac electrical phenotypes and highlight the challenges and opportunities offered by such studies.

The yield of postmortem genetic testing in sudden death cases with structural findings at autopsy.

Sudden cardiac death (SCD) is often associated with structural abnormalities of the heart during autopsy. This study sought to compare the diagnostic yield of postmortem genetic testing in (1) cases with structural findings of uncertain significance at autopsy to (2) cases with autopsy findings diagnostic of cardiomyopathy. We evaluated 57 SCD cases with structural findings at cardiac autopsy. Next-generation sequencing using a panel of 77 primary electrical disorder and cardiomyopathy genes was performed. Pathogenic and likely pathogenic variants were classified using American College of Medical Genetics (ACMG) consensus guidelines. In 29 cases (51%) autopsy findings of uncertain significance were identified whereas in 28 cases (49%) a diagnosis of cardiomyopathy was established. We identified a pathogenic or likely pathogenic variant in 10 cases (18%); in 1 (3%) case with non-specific autopsy findings compared with 9 (32%) cases with autopsy findings diagnostic of cardiomyopathy (p = 0.0054). The yield of genetic testing in SCD cases with autopsy findings consistent with cardiomyopathy is comparable with the yield in cardiomyopathy patients that are alive. Genetic testing in cases with findings of uncertain significance offers lower clinical utility than in cardiomyopathy, with lower yields than detected previously. This highlights the need for stringent evaluation of variant pathogenicity.

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 (PRSPR), QRS duration (PRSQRS), and Brugada syndrome (PRSBrS) 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 PRSPR, baseline PR, and female sex are associated with more pronounced PR slope, while PRSQRS and age are positively associated with QRS slope (P < 0.01 for all). PRSBrS, 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. PRSBrS, 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.

Blood Pressure-Associated Genetic Variants in the Natriuretic Peptide Receptor 1 Gene Modulate Guanylate Cyclase Activity.

BACKGROUND: Human genetic variation in the NPR1 (natriuretic peptide receptor 1 gene, encoding NPR-A, atrial natriuretic peptide receptor 1) was recently shown to affect blood pressure (BP). NPR-A catalyzes the intracellular conversion of guanosine triphosphate to cGMP (cyclic 3',5'-guanosine monophosphate) on binding of ANP, BNP (atrial or brain natriuretic peptide). Increased levels of cGMP decrease BP by inducing natriuresis, diuresis, and vasodilation. METHODS: We performed a meta-analysis of low-frequency and rare NPR1 variants for BP association in up to 491 584 unrelated individuals. To examine whether the identified BP-associated variants affect NPR-A function, the cGMP response to ANP and BNP was measured in cells expressing wild-type NPR1 and cells expressing the NPR1 variants. RESULTS: In this study, we identified BP associations of 3 amino acid altering variants of NPR1. The minor alleles of rs35479618 (p.E967K, gnomAD non-Finnish European allele frequency 0.017) and rs116245325 (p.L1034F, allele frequency 0.0007) were associated with higher BP (P=4.0×10-25 and P=9.9×10-8, respectively), while the minor allele of rs61757359 (p.G541S, allele frequency 0.003) was associated with lower BP (P=1.8×10-9). Cells transiently expressing 967K or 1034F NPR-A displayed decreased cGMP production in response to ANP and BNP (all P<10-6), while cells expressing 541S NPR-A produced more cGMP compared with cells expressing wild-type NPR-A (P≤4.13×10-5 for ANP and P≤4.24×10-3 for BNP). CONCLUSIONS: In summary, the loss or gain of guanylate cyclase activity for these NPR1 allelic variants could explain the higher or lower BP observed for carriers in large population-based studies.

GATA6 mutations: Characterization of two novel patients and a comprehensive overview of the GATA6 genotypic and phenotypic spectrum.

The first human mutations in GATA6 were described in a cohort of patients with persistent truncus arteriosus, and the phenotypic spectrum has expanded since then. This study underscores the broad phenotypic spectrum by presenting two patients with de novo GATA6 mutations, both exhibiting complex cardiac defects, pancreatic, and other abnormalities. Furthermore, we provided a detailed overview of all published human genetic variation in/near GATA6 published to date and the associated phenotypes (n = 78). We conclude that the most common phenotypes associated with a mutation in GATA6 were structural cardiac and pancreatic abnormalities, with a penetrance of 87 and 60%, respectively. Other common malformations were gallbladder agenesis, congenital diaphragmatic hernia, and neurocognitive abnormalities, mostly developmental delay. Fifty-eight percent of the mutations were de novo, and these patients more often had an anomaly of intracardiac connections, an anomaly of the great arteries, and hypothyroidism, compared with those with inherited mutations. Functional studies mostly support loss-of-function as the pathophysiological mechanism. In conclusion, GATA6 mutations give a wide range of phenotypic defects, most frequently malformations of the heart and pancreas. This highlights the importance of detailed clinical evaluation of identified carriers to evaluate their full phenotypic spectrum.

Calmodulin mutations and life-threatening cardiac arrhythmias: insights from the International Calmodulinopathy Registry.

AIMS: Calmodulinopathies are rare life-threatening arrhythmia syndromes which affect mostly young individuals and are, caused by mutations in any of the three genes (CALM 1-3) that encode identical calmodulin proteins. We established the International Calmodulinopathy Registry (ICalmR) to understand the natural history, clinical features, and response to therapy of patients with a CALM-mediated arrhythmia syndrome. METHODS AND RESULTS: A dedicated Case Report File was created to collect demographic, clinical, and genetic information. ICalmR has enrolled 74 subjects, with a variant in the CALM1 (n = 36), CALM2 (n = 23), or CALM3 (n = 15) genes. Sixty-four (86.5%) were symptomatic and the 10-year cumulative mortality was 27%. The two prevalent phenotypes are long QT syndrome (LQTS; CALM-LQTS, n = 36, 49%) and catecholaminergic polymorphic ventricular tachycardia (CPVT; CALM-CPVT, n = 21, 28%). CALM-LQTS patients have extremely prolonged QTc intervals (594 ± 73 ms), high prevalence (78%) of life-threatening arrhythmias with median age at onset of 1.5 years [interquartile range (IQR) 0.1-5.5 years] and poor response to therapies. Most electrocardiograms (ECGs) show late onset peaked T waves. All CALM-CPVT patients were symptomatic with median age of onset of 6.0 years (IQR 3.0-8.5 years). Basal ECG frequently shows prominent U waves. Other CALM-related phenotypes are idiopathic ventricular fibrillation (IVF, n = 7), sudden unexplained death (SUD, n = 4), overlapping features of CPVT/LQTS (n = 3), and predominant neurological phenotype (n = 1). Cardiac structural abnormalities and neurological features were present in 18 and 13 patients, respectively. CONCLUSION: Calmodulinopathies are largely characterized by adrenergically-induced life-threatening arrhythmias. Available therapies are disquietingly insufficient, especially in CALM-LQTS. Combination therapy with drugs, sympathectomy, and devices should be considered.

Homozygous frameshift mutations in FAT1 cause a syndrome characterized by colobomatous-microphthalmia, ptosis, nephropathy and syndactyly.

A failure in optic fissure fusion during development can lead to blinding malformations of the eye. Here, we report a syndrome characterized by facial dysmorphism, colobomatous microphthalmia, ptosis and syndactyly with or without nephropathy, associated with homozygous frameshift mutations in FAT1. We show that Fat1 knockout mice and zebrafish embryos homozygous for truncating fat1a mutations exhibit completely penetrant coloboma, recapitulating the most consistent developmental defect observed in affected individuals. In human retinal pigment epithelium (RPE) cells, the primary site for the fusion of optic fissure margins, FAT1 is localized at earliest cell-cell junctions, consistent with a role in facilitating optic fissure fusion during vertebrate eye development. Our findings establish FAT1 as a gene with pleiotropic effects in human, in that frameshift mutations cause a severe multi-system disorder whereas recessive missense mutations had been previously associated with isolated glomerulotubular nephropathy.

Trans-ethnic association study of blood pressure determinants in over 750,000 individuals.

In this trans-ethnic multi-omic study, we reinterpret the genetic architecture of blood pressure to identify genes, tissues, phenomes and medication contexts of blood pressure homeostasis. We discovered 208 novel common blood pressure SNPs and 53 rare variants in genome-wide association studies of systolic, diastolic and pulse pressure in up to 776,078 participants from the Million Veteran Program (MVP) and collaborating studies, with analysis of the blood pressure clinical phenome in MVP. Our transcriptome-wide association study detected 4,043 blood pressure associations with genetically predicted gene expression of 840 genes in 45 tissues, and mouse renal single-cell RNA sequencing identified upregulated blood pressure genes in kidney tubule cells.