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.

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.

Systematic large-scale assessment of the genetic architecture of left ventricular noncompaction reveals diverse etiologies.

PURPOSE: To characterize the genetic architecture of left ventricular noncompaction (LVNC) and investigate the extent to which it may represent a distinct pathology or a secondary phenotype associated with other cardiac diseases. METHODS: We performed rare variant association analysis with 840 LVNC cases and 125,748 gnomAD population controls, and compared results to similar analyses on dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). RESULTS: We observed substantial genetic overlap indicating that LVNC often represents a phenotypic variation of DCM or HCM. In contrast, truncating variants in MYH7, ACTN2, and PRDM16 were uniquely associated with LVNC and may reflect a distinct LVNC etiology. In particular, MYH7 truncating variants (MYH7tv), generally considered nonpathogenic for cardiomyopathies, were 20-fold enriched in LVNC cases over controls. MYH7tv heterozygotes identified in the UK Biobank and healthy volunteer cohorts also displayed significantly greater noncompaction compared with matched controls. RYR2 exon deletions and HCN4 transmembrane variants were also enriched in LVNC, supporting prior reports of association with arrhythmogenic LVNC phenotypes. CONCLUSION: LVNC is characterized by substantial genetic overlap with DCM/HCM but is also associated with distinct noncompaction and arrhythmia etiologies. These results will enable enhanced application of LVNC genetic testing and help to distinguish pathological from physiological noncompaction.

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.

Functional modulation of atrio-ventricular conduction by enhanced late sodium current and calcium-dependent mechanisms in Scn5a1798insD/+ mice.

AIMS: SCN5A mutations are associated with arrhythmia syndromes, including Brugada syndrome, long QT syndrome type 3 (LQT3), and cardiac conduction disease. Long QT syndrome type 3 patients display atrio-ventricular (AV) conduction slowing which may contribute to arrhythmogenesis. We here investigated the as yet unknown underlying mechanisms. METHODS AND RESULTS: We assessed electrophysiological and molecular alterations underlying AV-conduction abnormalities in mice carrying the Scn5a1798insD/+ mutation. Langendorff-perfused Scn5a1798insD/+ hearts showed prolonged AV-conduction compared to wild type (WT) without changes in atrial and His-ventricular (HV) conduction. The late sodium current (INa,L) inhibitor ranolazine (RAN) normalized AV-conduction in Scn5a1798insD/+ mice, likely by preventing the mutation-induced increase in intracellular sodium ([Na+]i) and calcium ([Ca2+]i) concentrations. Indeed, further enhancement of [Na+]i and [Ca2+]i by the Na+/K+-ATPase inhibitor ouabain caused excessive increase in AV-conduction time in Scn5a1798insD/+ hearts. Scn5a1798insD/+ mice from the 129P2 strain displayed more severe AV-conduction abnormalities than FVB/N-Scn5a1798insD/+ mice, in line with their larger mutation-induced INa,L. Transverse aortic constriction (TAC) caused excessive prolongation of AV-conduction in FVB/N-Scn5a1798insD/+ mice (while HV-intervals remained unchanged), which was prevented by chronic RAN treatment. Scn5a1798insD/+-TAC hearts showed decreased mRNA levels of conduction genes in the AV-nodal region, but no structural changes in the AV-node or His bundle. In Scn5a1798insD/+-TAC mice deficient for the transcription factor Nfatc2 (effector of the calcium-calcineurin pathway), AV-conduction and conduction gene expression were restored to WT levels. CONCLUSIONS: Our findings indicate a detrimental role for enhanced INa,L and consequent calcium dysregulation on AV-conduction in Scn5a1798insD/+ mice, providing evidence for a functional mechanism underlying AV-conduction disturbances secondary to gain-of-function SCN5A mutations.

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.

GNB5 Mutations Cause an Autosomal-Recessive Multisystem Syndrome with Sinus Bradycardia and Cognitive Disability.

GNB5 encodes the G protein ß subunit 5 and is involved in inhibitory G protein signaling. Here, we report mutations in GNB5 that are associated with heart-rate disturbance, eye disease, intellectual disability, gastric problems, hypotonia, and seizures in nine individuals from six families. We observed an association between the nature of the variants and clinical severity; individuals with loss-of-function alleles had more severe symptoms, including substantial developmental delay, speech defects, severe hypotonia, pathological gastro-esophageal reflux, retinal disease, and sinus-node dysfunction, whereas related heterozygotes harboring missense variants presented with a clinically milder phenotype. Zebrafish gnb5 knockouts recapitulated the phenotypic spectrum of affected individuals, including cardiac, neurological, and ophthalmological abnormalities, supporting a direct role of GNB5 in the control of heart rate, hypotonia, and vision.

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.

The Brugada Syndrome Susceptibility Gene HEY2 Modulates Cardiac Transmural Ion Channel Patterning and Electrical Heterogeneity.

RATIONALE: Genome-wide association studies previously identified an association of rs9388451 at chromosome 6q22.3 (near HEY2) with Brugada syndrome. The causal gene and underlying mechanism remain unresolved. OBJECTIVE: We used an integrative approach entailing transcriptomic studies in human hearts and electrophysiological studies in Hey2+/- (Hey2 heterozygous knockout) mice to dissect the underpinnings of the 6q22.31 association with Brugada syndrome. METHODS AND RESULTS: We queried expression quantitative trait locus data acquired in 190 human left ventricular samples from the genotype-tissue expression consortium for cis-expression quantitative trait locus effects of rs9388451, which revealed an association between Brugada syndrome risk allele dosage and HEY2 expression (ß=+0.159; P=0.0036). In the same transcriptomic data, we conducted genome-wide coexpression analysis for HEY2, which uncovered KCNIP2, encoding the ß-subunit of the channel underlying the transient outward current (Ito), as the transcript most robustly correlating with HEY2 expression (ß=+1.47; P=2×10-34). Transcript abundance of Hey2 and the Ito subunits Kcnip2 and Kcnd2, assessed by quantitative reverse transcription-polymerase chain reaction, was higher in subepicardium versus subendocardium in both left and right ventricles, with lower levels in Hey2+/- mice compared with wild type. Surface ECG measurements showed less prominent J waves in Hey2+/- mice compared with wild-type. In wild-type mice, patch-clamp electrophysiological studies on cardiomyocytes from right ventricle demonstrated a shorter action potential duration and a lower Vmax in subepicardium compared with subendocardium cardiomyocytes, which was paralleled by a higher Ito and a lower sodium current (INa) density in subepicardium versus subendocardium. These transmural differences were diminished in Hey2+/- mice because of changes in subepicardial cardiomyocytes. CONCLUSIONS: This study uncovers a role of HEY2 in the normal transmural electrophysiological gradient in the ventricle and provides compelling evidence that genetic variation at 6q22.31 (rs9388451) is associated with Brugada syndrome through a HEY2-dependent alteration of ion channel expression across the cardiac ventricular wall.

Sudden cardiac arrest associated with use of a non-cardiac drug that reduces cardiac excitability: evidence from bench, bedside, and community.

AIMS: Non-cardiac drugs that impair cardiac repolarization (electrocardiographic QT prolongation) are associated with an increased sudden cardiac arrest (SCA) risk. Emerging evidence suggests that non-cardiac drugs that impair cardiac depolarization and excitability (electrocardiographic QRS prolongation) also increase the risk for SCA. Nortriptyline, which blocks the SCN5A-encoded cardiac sodium channel, may exemplify such drugs. We aimed to study whether nortriptyline increases the risk for SCA, and to establish the underlying mechanisms. METHODS AND RESULTS: We studied QRS durations during rest/exercise in an index patient who experienced ventricular tachycardia during exercise while using nortriptyline, and compared them with those of 55 controls with/without nortriptyline and 24 controls with Brugada syndrome (BrS) without nortriptyline, who carried an SCN5A mutation. We performed molecular-genetic (exon-trapping) and functional (patch-clamp) experiments to unravel the mechanisms of QRS prolongation by nortriptyline and the SCN5A mutation found in the index patient. We conducted a prospective community-based study among 944 victims of ECG-documented SCA and 4354-matched controls to determine the risk for SCA associated with nortriptyline use. Multiple mechanisms may act in concert to increase the risk for SCA during nortriptyline use. Pharmacological (nortriptyline), genetic (loss-of-function SCN5A mutation), and/or functional (sodium channel inactivation at fast heart rates) factors conspire to reduce the cardiac sodium current and increase the risk for SCA. Nortriptyline use in the community was associated with a 4.5-fold increase in the risk for SCA [adjusted OR: 4.5 (95% CI: 1.1-19.5)], particularly when other sodium channel-blocking factors were present. CONCLUSIONS: Nortriptyline increases the risk for SCA in the general population, particularly in the presence of genetic and/or non-genetic factors that decrease cardiac excitability by blocking the cardiac sodium channel.

Biallelic variants in POPDC2 cause a novel autosomal recessive syndrome presenting with cardiac conduction defects and variable hypertrophic cardiomyopathy.

POPDC2 encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss of Popdc2 in mice results in sinus pauses and bradycardia and morpholino knockdown of popdc2 in zebrafish results in atrioventricular (AV) block. We identified bi-allelic variants in POPDC2 in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identified POPDC2 variants are predicted to diminish the ability of POPDC2 to bind cAMP. In in vitro electrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1-POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants in POPDC1 and POPDC2 and presence of sinus node disease in POPDC2, but not in POPDC1 related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence for POPDC2 as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functional POPDC2 display sinus and AV node dysfunction.

Genetic Burden of TNNI3K in Diagnostic Testing of Patients With Dilated Cardiomyopathy and Supraventricular Arrhythmias.

BACKGROUND: Genetic variants in TNNI3K (troponin-I interacting kinase) have previously been associated with dilated cardiomyopathy (DCM), cardiac conduction disease, and supraventricular tachycardias. However, the link between TNNI3K variants and these cardiac phenotypes shows a lack of consensus concerning phenotype and protein function. METHODS: We describe a systematic retrospective study of a cohort of patients undergoing genetic testing for cardiac arrhythmias and cardiomyopathy including TNNI3K. We further performed burden testing of TNNI3K in the UK Biobank. For 2 novel TNNI3K variants, we tested cosegregation. TNNI3K kinase function was estimated by TNNI3K autophosphorylation assays. RESULTS: We demonstrate enrichment of rare coding TNNI3K variants in DCM patients in the Amsterdam cohort. In the UK Biobank, we observed an association between TNNI3K missense (but not loss-of-function) variants and DCM and atrial fibrillation. Furthermore, we demonstrate genetic segregation for 2 rare variants, TNNI3K-p.Ile512Thr and TNNI3K-p.His592Tyr, with phenotypes consisting of DCM, cardiac conduction disease, and supraventricular tachycardia, together with increased autophosphorylation. In contrast, TNNI3K-p.Arg556_Asn590del, a likely benign variant, demonstrated depleted autophosphorylation. CONCLUSIONS: Our findings demonstrate an increased burden of rare coding TNNI3K variants in cardiac patients with DCM. Furthermore, we present 2 novel likely pathogenic TNNI3K variants with increased autophosphorylation, suggesting that enhanced autophosphorylation is likely to drive pathogenicity.

Haplotype-sharing analysis implicates chromosome 7q36 harboring DPP6 in familial idiopathic ventricular fibrillation.

Idiopathic Ventricular Fibrillation (IVF) is defined as spontaneous VF without any known structural or electrical heart disease. A family history is present in up to 20% of probands with the disorder, suggesting that at least a subset of IVF is hereditary. A genome-wide haplotype-sharing analysis was performed for identification of the responsible gene in three distantly related families in which multiple individuals died suddenly or were successfully resuscitated at young age. We identified a haplotype, on chromosome 7q36, that was conserved in these three families and was also shared by 7 of 42 independent IVF patients. The shared chromosomal segment harbors part of the DPP6 gene, which encodes a putative component of the transient outward current in the heart. We demonstrated a 20-fold increase in DPP6 mRNA levels in the myocardium of carriers as compared to controls. Clinical evaluation of 84 risk-haplotype carriers and 71 noncarriers revealed no ECG or structural parameters indicative of cardiac disease. Penetrance of IVF was high; 50% of risk-haplotype carriers experienced (aborted) sudden cardiac death before the age of 58 years. We propose DPP6 as a gene for IVF and increased DPP6 expression as the likely pathogenetic mechanism.

Chronically elevated branched chain amino acid levels are pro-arrhythmic.

AIMS: Cardiac arrhythmias comprise a major health and economic burden and are associated with significant morbidity and mortality, including cardiac failure, stroke, and sudden cardiac death (SCD). Development of efficient preventive and therapeutic strategies is hampered by incomplete knowledge of disease mechanisms and pathways. Our aim is to identify novel mechanisms underlying cardiac arrhythmia and SCD using an unbiased approach. METHODS AND RESULTS: We employed a phenotype-driven N-ethyl-N-nitrosourea mutagenesis screen and identified a mouse line with a high incidence of sudden death at young age (6-9 weeks) in the absence of prior symptoms. Affected mice were found to be homozygous for the nonsense mutation Bcat2p.Q300*/p.Q300* in the Bcat2 gene encoding branched chain amino acid transaminase 2. At the age of 4-5 weeks, Bcat2p.Q300*/p.Q300* mice displayed drastic increase of plasma levels of branch chain amino acids (BCAAs-leucine, isoleucine, valine) due to the incomplete catabolism of BCAAs, in addition to inducible arrhythmias ex vivo as well as cardiac conduction and repolarization disturbances. In line with these findings, plasma BCAA levels were positively correlated to electrocardiogram indices of conduction and repolarization in the German community-based KORA F4 Study. Isolated cardiomyocytes from Bcat2p.Q300*/p.Q300* mice revealed action potential (AP) prolongation, pro-arrhythmic events (early and late afterdepolarizations, triggered APs), and dysregulated calcium homeostasis. Incubation of human pluripotent stem cell-derived cardiomyocytes with elevated concentration of BCAAs induced similar calcium dysregulation and pro-arrhythmic events which were prevented by rapamycin, demonstrating the crucial involvement of mTOR pathway activation. CONCLUSIONS: Our findings identify for the first time a causative link between elevated BCAAs and arrhythmia, which has implications for arrhythmogenesis in conditions associated with BCAA metabolism dysregulation such as diabetes, metabolic syndrome, and heart failure.

Quantitative trait loci for electrocardiographic parameters and arrhythmia in the mouse.

Cardiac arrhythmias associated with sudden death are influenced by multiple biological pathways and are modulated by numerous genetic and environmental factors. Elevated heart rate and prolonged ECG indices of conduction and repolarization have been associated with risk of sudden death. Insight into the genetic underpinnings of these parameters thus provides an important means to the dissection of the genetic components modulating risk of sudden cardiac death. In this study we mapped quantitative trait loci (QTL) modulating heart rate, ECG indices of conduction and repolarization, and susceptibility to arrhythmia, in a conduction disease-sensitized F(2) mouse population. Heart rate, P-duration, PR-, QRS- and QT-interval were measured at baseline (n=502) and after flecainide administration (n=370) in mutant F(2) progeny (F(2)-MUT) resulting from the FVB/NJ-Scn5a1798(insD/+) X 129P2-Scn5a1798(insD/+) mouse cross. Episodes of sinus arrhythmia and ventricular tachyarrhythmia occurring post-flecainide were treated as binary traits. F(2)-MUT mice were genotyped using a genome-wide 768 single nucleotide polymorphism (SNP) panel. Interval mapping uncovered multiple QTL for ECG parameters and arrhythmia. A sex-interacting scan identified QTL displaying sex-dependency, and a two-dimensional QTL scan unmasked locus-locus (epistasis) interactions influencing ECG traits. A number of QTL coincided at specific chromosomal locations, suggesting pleiotropic effects at these loci. Through transcript profiling in myocardium from the parental mouse strains we identified genes co-localizing at the identified QTL that constitute highly relevant candidates for the observed effects. The detection of QTL influencing ECG indices and arrhythmia is an essential step towards identifying genetic networks for sudden, arrhythmic, cardiac death.

Genetically determined differences in sodium current characteristics modulate conduction disease severity in mice with cardiac sodium channelopathy.

Conduction slowing of the electric impulse that drives the heartbeat may evoke lethal cardiac arrhythmias. Mutations in SCN5A, which encodes the pore-forming cardiac sodium channel alpha subunit, are associated with familial arrhythmia syndromes based on conduction slowing. However, disease severity among mutation carriers is highly variable. We hypothesized that genetic modifiers underlie the variability in conduction slowing and disease severity. With the aim of identifying such modifiers, we studied the Scn5a(1798insD/+) mutation in 2 distinct mouse strains, FVB/N and 129P2. In 129P2 mice, the mutation resulted in more severe conduction slowing particularly in the right ventricle (RV) compared to FVB/N. Pan-genomic mRNA expression profiling in the 2 mouse strains uncovered a drastic reduction in mRNA encoding the sodium channel auxiliary subunit beta4 (Scn4b) in 129P2 mice compared to FVB/N. This corresponded to low to undetectable beta4 protein levels in 129P2 ventricular tissue, whereas abundant beta4 protein was detected in FVB/N. Sodium current measurements in isolated myocytes from the 2 mouse strains indicated that sodium channel activation in myocytes from 129P2 mice occurred at more positive potentials compared to FVB/N. Using computer simulations, this difference in activation kinetics was predicted to explain the observed differences in conduction disease severity between the 2 strains. In conclusion, genetically determined differences in sodium current characteristics on the myocyte level modulate disease severity in cardiac sodium channelopathies. In particular, the sodium channel subunit beta4 (SCN4B) may constitute a potential genetic modifier of conduction and cardiac sodium channel disease.

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.

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.