Transcriptional response of the heart to vagus nerve stimulation

Heart failure is a major clinical problem, with treatments involving medication, devices, and emerging neuromodulation therapies such as vagus nerve stimulation (VNS). Considering the ongoing interest in using VNS to treat cardiovascular disease, it is important to understand the genetic and molecular changes developing in the heart in response to this form of autonomic neuromodulation. This experimental animal (rat) study investigated the immediate transcriptional response of the ventricular myocardium to selective stimulation of vagal efferent activity using an optogenetic approach. Vagal preganglionic neurons in the dorsal motor nucleus of the vagus nerve were genetically targeted to express light-sensitive chimeric channelrhodopsin variant ChIEF and stimulated using light. RNA sequencing of the left ventricular myocardium identified 294 differentially expressed genes (false discovery rate < 0.05). Qiagen Ingenuity Pathway Analysis (IPA) highlighted 118 canonical pathways that were significantly modulated by vagal activity, of which 14 had a z score of ≥2/≤-2, including EIF-2, IL-2, integrin, and NFAT-regulated cardiac hypertrophy. IPA revealed the effect of efferent vagus stimulation on protein synthesis, autophagy, fibrosis, autonomic signaling, inflammation, and hypertrophy. IPA further predicted that the identified differentially expressed genes were the targets of 50 upstream regulators, including transcription factors (e.g., MYC and NRF1) and microRNAs (e.g., miR-335-3p and miR-338-3p). These data demonstrate that the vagus nerve has a major impact on the myocardial expression of genes involved in the regulation of key biological pathways. The transcriptional response of the ventricular myocardium induced by stimulation of vagal efferents is consistent with the beneficial effect of maintained/increased vagal activity on the heart.NEW & NOTEWORTHY This experimental animal study investigated the immediate transcriptional response of the ventricular myocardium to selective stimulation of vagal efferent activity. Vagal stimulation induced significant transcriptional changes in the heart involving the pathways controlling autonomic signaling, inflammation, fibrosis, and hypertrophy. This study provides the first direct evidence that myocardial gene expression is modulated by the activity of the autonomic nervous system.

Prognostic Significance of Different Ventricular Ectopic Burdens During Submaximal Exercise in Asymptomatic UK Biobank Subjects

BACKGROUND

The consequences of exercise-induced premature ventricular contractions (PVCs) in asymptomatic individuals remain unclear. This study aimed to assess the association between PVC burdens during submaximal exercise and major adverse cardiovascular events (MI/HF/LTVA: myocardial infarction [MI], heart failure [HF], and life-threatening ventricular arrhythmia [LTVA]), and all-cause mortality. Additional end points were MI, LTVA, HF, and cardiovascular mortality.

METHODS

A neural network was developed to count PVCs from ECGs recorded during exercise (6 minutes) and recovery (1 minute) in 48 315 asymptomatic participants from UK Biobank. Associations were estimated using multivariable Cox proportional hazard models. Explorative studies were conducted in subgroups with cardiovascular magnetic resonance imaging data (n=6290) and NT-proBNP (N-terminal Pro-B-type natriuretic peptide) levels (n=4607) to examine whether PVC burden was associated with subclinical cardiomyopathy.

RESULTS

Mean age was 56.8±8.2 years; 51.1% of the participants were female; and median follow-up was 12.6 years. Low PVC counts during exercise and recovery were both associated with MI/HF/LTVA risk, independently of clinical factors: adjusted hazard ratio (HR), 1.2 (1-5 exercise PVCs, P<0.001) and HR, 1.3 (1-5 recovery PVCs, P<0.001). Risks were higher with increasing PVC count: HR, 1.8 (>20 exercise PVCs, P<0.001) and HR, 1.6 (>5 recovery PVCs, P<0.001). A similar trend was observed for all-cause mortality, although associations were only significant for high PVC burdens: HRs, 1.6 (>20 exercise PVCs, P<0.001) and 1.5 (>5 recovery PVCs, P<0.001). Complex PVC rhythms were associated with higher risk compared with PVC count alone. PVCs were also associated with incident HF, LTVA, and cardiovascular mortality, but not MI. In the explorative studies, high PVC burden was associated with larger left ventricular volumes, lower ejection fraction, and higher levels of NT-proBNP compared with participants without PVCs.

CONCLUSIONS

In this cohort of middle-aged and older adults, PVC count during submaximal exercise and recovery were both associated with MI/HF/LTVA, all-cause mortality, HF, LTVAs, and cardiovascular mortality, independent of clinical and exercise test factors, indicating an incremental increase in risk as PVC count rises. Complex PVC rhythms were associated with higher risk compared with PVC count alone. Underlying mechanisms may include the presence of subclinical cardiomyopathy.

Long-term association of ultra-short heart rate variability with cardiovascular events

Heart rate variability (HRV) is a cardiac autonomic marker with predictive value in cardiac patients. Ultra-short HRV (usHRV) can be measured at scale using standard and wearable ECGs, but its association with cardiovascular events in the general population is undetermined. We aimed to validate usHRV measured using ≤ 15-s ECGs (using RMSSD, SDSD and PHF indices) and investigate its association with atrial fibrillation, major adverse cardiac events, stroke and mortality in individuals without cardiovascular disease. In the National Survey for Health and Development (n = 1337 participants), agreement between 15-s and 6-min HRV, assessed with correlation analysis and Bland-Altman plots, was very good for RMSSD and SDSD and good for PHF. In the UK Biobank (n = 51,628 participants, 64% male, median age 58), after a median follow-up of 11.5 (11.4-11.7) years, incidence of outcomes ranged between 1.7% and 4.3%. Non-linear Cox regression analysis showed that reduced usHRV from 15-, 10- and 5-s ECGs was associated with all outcomes. Individuals with low usHRV (< 20th percentile) had hazard ratios for outcomes between 1.16 and 1.29, p < 0.05, with respect to the reference group. In conclusion, usHRV from ≤ 15-s ECGs correlates with standard short-term HRV and predicts increased risk of cardiovascular events in a large population-representative cohort.

Immediate and sustained increases in the activity of vagal preganglionic neurons during exercise and after exercise training

AIMS

The brain controls the heart by dynamic recruitment and withdrawal of cardiac parasympathetic (vagal) and sympathetic activity. Autonomic control is essential for the development of cardiovascular responses during exercise, however, the patterns of changes in the activity of the two autonomic limbs, and their functional interactions in orchestrating physiological responses during exercise, are not fully understood. The aim of this study was to characterize changes in vagal parasympathetic drive in response to exercise and exercise training by directly recording the electrical activity of vagal preganglionic neurons in experimental animals (rats).

METHODS AND RESULTS

Single unit recordings were made using carbon-fibre microelectrodes from the populations of vagal preganglionic neurons of the nucleus ambiguus (NA) and the dorsal vagal motor nucleus of the brainstem. It was found that (i) vagal preganglionic neurons of the NA and the dorsal vagal motor nucleus are strongly activated during bouts of acute exercise, and (ii) exercise training markedly increases the resting activity of both populations of vagal preganglionic neurons and augments the excitatory responses of NA neurons during exercise.

CONCLUSIONS

These data show that central vagal drive increases during exercise and provide the first direct neurophysiological evidence that exercise training increases vagal tone. The data argue against the notion of exercise-induced central vagal withdrawal during exercise. We propose that robust increases in the activity of vagal preganglionic neurons during bouts of exercise underlie activity-dependent plasticity, leading to higher resting vagal tone that confers multiple health benefits associated with regular exercise.

Hypoxia Promotes Atrial Tachyarrhythmias via Opening of ATP-Sensitive Potassium Channels

BACKGROUND

Hypoxia-ischemia predisposes to atrial arrhythmia. Atrial ATP-sensitive potassium channel (KATP) modulation during hypoxia has not been explored. We investigated the effects of hypoxia on atrial electrophysiology in mice with global deletion of KATP pore-forming subunits.

METHODS

Whole heart KATP RNA expression was probed. Whole-cell KATP current and action potentials were recorded in isolated wild-type (WT), Kir6.1 global knockout (6.1-gKO), and Kir6.2 global knockout (6.2-gKO) murine atrial myocytes. Langendorff-perfused hearts were assessed for atrial effective refractory period (ERP), conduction velocity, wavefront path length (WFPL), and arrhymogenicity under normoxia/hypoxia using a microelectrode array and programmed electrical stimulation. Heart histology was assessed.

RESULTS

Expression patterns were essentially identical for all KATP subunit RNA across human heart, whereas in mouse, Kir6.1 and SUR2 (sulphonylurea receptor subunit) were higher in ventricle than atrium, and Kir6.2 and SUR1 were higher in atrium. Compared with WT, 6.2-gKO atrial myocytes had reduced tolbutamide-sensitive current and action potentials were more depolarized with slower upstroke and reduced peak amplitude. Action potential duration was prolonged in 6.1-gKO atrial myocytes, absent of changes in other ion channel gene expression or atrial myocyte hypertrophy. In Langendorff-perfused hearts, baseline atrial ERP was prolonged and conduction velocity reduced in both KATP knockout mice compared with WT, without histological fibrosis. Compared with baseline, hypoxia led to conduction velocity slowing, stable ERP, and WFPL shortening in WT and 6.1-gKO hearts, whereas WFPL was stable in 6.2-gKO hearts due to ERP prolongation with conduction velocity slowing. Tolbutamide reversed hypoxia-induced WFPL shortening in WT and 6.1-gKO hearts through ERP prolongation. Atrial tachyarrhythmias inducible with programmed electrical stimulation during hypoxia in WT and 6.1-gKO mice correlated with WFPL shortening. Spontaneous arrhythmia was not seen.

CONCLUSIONS

KATP block/absence leads to cellular and tissue level atrial electrophysiological modification. Kir6.2 global knockout prevents hypoxia-induced atrial WFPL shortening and atrial arrhythmogenicity to programmed electrical stimulation. This mechanism could be explored translationally to treat ischemically driven atrial arrhythmia.

Retinoic acid signaling modulation guides in vitro specification of human heart field-specific progenitor pools

Cardiogenesis relies on the precise spatiotemporal coordination of multiple progenitor populations. Understanding the specification and differentiation of these distinct progenitor pools during human embryonic development is crucial for advancing our knowledge of congenital cardiac malformations and designing new regenerative therapies. By combining genetic labelling, single-cell transcriptomics, and ex vivo human-mouse embryonic chimeras we uncovered that modulation of retinoic acid signaling instructs human pluripotent stem cells to form heart field-specific progenitors with distinct fate potentials. In addition to the classical first and second heart fields, we observed the appearance of juxta-cardiac field progenitors giving rise to both myocardial and epicardial cells. Applying these findings to stem-cell based disease modelling we identified specific transcriptional dysregulation in first and second heart field progenitors derived from stem cells of patients with hypoplastic left heart syndrome. This highlights the suitability of our in vitro differentiation platform for studying human cardiac development and disease.

Genetic architecture of spatial electrical biomarkers for cardiac arrhythmia and relationship with cardiovascular disease

The 3-dimensional spatial and 2-dimensional frontal QRS-T angles are measures derived from the vectorcardiogram. They are independent risk predictors for arrhythmia, but the underlying biology is unknown. Using multi-ancestry genome-wide association studies we identify 61 (58 previously unreported) loci for the spatial QRS-T angle (N = 118,780) and 11 for the frontal QRS-T angle (N = 159,715). Seven out of the 61 spatial QRS-T angle loci have not been reported for other electrocardiographic measures. Enrichments are observed in pathways related to cardiac and vascular development, muscle contraction, and hypertrophy. Pairwise genome-wide association studies with classical ECG traits identify shared genetic influences with PR interval and QRS duration. Phenome-wide scanning indicate associations with atrial fibrillation, atrioventricular block and arterial embolism and genetically determined QRS-T angle measures are associated with fascicular and bundle branch block (and also atrioventricular block for the frontal QRS-T angle). We identify potential biology involved in the QRS-T angle and their genetic relationships with cardiovascular traits and diseases, may inform future research and risk prediction.

Premature atrial and ventricular contractions detected on wearable-format electrocardiograms and prediction of cardiovascular events

Aims

Wearable devices are transforming the electrocardiogram (ECG) into a ubiquitous medical test. This study assesses the association between premature ventricular and atrial contractions (PVCs and PACs) detected on wearable-format ECGs (15 s single lead) and cardiovascular outcomes in individuals without cardiovascular disease (CVD).

Methods and results

Premature atrial contractions and PVCs were identified in 15 s single-lead ECGs from N = 54 016 UK Biobank participants (median age, interquartile range, age 58, 50-63 years, 54% female). Cox regression models adjusted for traditional risk factors were used to determine associations with atrial fibrillation (AF), heart failure (HF), myocardial infarction (MI), stroke, life-threatening ventricular arrhythmias (LTVAs), and mortality over a period of 11.5 (11.4-11.7) years. The strongest associations were found between PVCs (prevalence 2.2%) and HF (hazard ratio, HR, 95% confidence interval = 2.09, 1.58-2.78) and between PACs (prevalence 1.9%) and AF (HR = 2.52, 2.11-3.01), with shorter prematurity further increasing risk. Premature ventricular contractions and PACs were also associated with LTVA (P < 0.05). Associations with MI, stroke, and mortality were significant only in unadjusted models. In a separate UK Biobank sub-study sample [UKB-2, N = 29,324, age 64, 58-60 years, 54% female, follow-up 3.5 (2.6-4.8) years] used for independent validation, after adjusting for risk factors, PACs were associated with AF (HR = 1.80, 1.12-2.89) and PVCs with HF (HR = 2.32, 1.28-4.22).

Conclusion

In middle-aged individuals without CVD, premature contractions identified in 15 s single-lead ECGs are strongly associated with an increased risk of AF and HF. These data warrant further investigation to assess the role of wearable ECGs for early cardiovascular risk stratification.

Omega-3 polyunsaturated fatty acid-induced vasodilation in mouse aorta and mesenteric arteries is not mediated by ATP-sensitive potassium channels

There is strong evidence that the omega-3 polyunsaturated fatty acids (n-3 PUFAs) docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) have cardioprotective effects. n-3 PUFAs cause vasodilation in hypertensive patients, in part controlled by increased membrane conductance to potassium. As KATP channels play a major role in vascular tone regulation and are involved in hypertension, we aimed to verify whether n-3 PUFA-mediated vasodilation involved the opening of KATP channels. We used a murine model in which the KATP channel pore subunit, Kir6.1, is deleted in vascular smooth muscle. The vasomotor response of preconstricted arteries to physiologically relevant concentrations of DHA and EPA was measured using wire myography, using the channel blocker PNU-37883A. The effect of n-3 PUFAs on potassium currents in wild-type native smooth muscle cells was investigated using whole-cell patch clamping. DHA and EPA induced vasodilation in mouse aorta and mesenteric arteries; relaxations in the aorta were sensitive to KATP blockade with PNU-37883A. Endothelium removal didn't affect relaxation to EPA and caused a small but significant inhibition of relaxation to DHA. In the knock-out model, relaxations to DHA and EPA were unaffected by channel knockdown but were still inhibited by PNU-37883A, indicating that the action of PNU-37883A on relaxation may not reflect inhibition of KATP. In native aortic smooth muscle cells DHA failed to activate KATP currents. We conclude that DHA and EPA cause vasodilation in mouse aorta and mesenteric arteries. Relaxations in blocker-treated arteries from knock-out mice demonstrate that KATP channels are not involved in the n-3 PUFA-induced relaxation.

Role of cardiac mitofusins in cardiac conduction following simulated ischemia-reperfusion

Mitochondrial dysfunction induced by acute cardiac ischemia-reperfusion (IR), may increase susceptibility to arrhythmias by perturbing energetics, oxidative stress production and calcium homeostasis. Although changes in mitochondrial morphology are known to impact on mitochondrial function, their role in cardiac arrhythmogenesis is not known. To assess action potential duration (APD) in cardiomyocytes from the Mitofusins-1/2 (Mfn1/Mfn2)-double-knockout (Mfn-DKO) compared to wild-type (WT) mice, optical-electrophysiology was conducted. To measure conduction velocity (CV) in atrial and ventricular tissue from the Mfn-DKO and WT mice, at both baseline and following simulated acute IR, multi-electrode array (MEA) was employed. Intracellular localization of connexin-43 (Cx43) at baseline was evaluated by immunohistochemistry, while Cx-43 phosphorylation was assessed by Western-blotting. Mfn-DKO cardiomyocytes demonstrated an increased APD. At baseline, CV was significantly lower in the left ventricle of the Mfn-DKO mice. CV decreased with simulated-ischemia and returned to baseline levels during simulated-reperfusion in WT but not in atria of Mfn-DKO mice. Mfn-DKO hearts displayed increased Cx43 lateralization, although phosphorylation of Cx43 at Ser-368 did not differ. In summary, Mfn-DKO mice have increased APD and reduced CV at baseline and impaired alterations in CV following cardiac IR. These findings were associated with increased Cx43 lateralization, suggesting that the mitofusins may impact on post-MI cardiac-arrhythmogenesis.

CRISPR-mediated correction of skeletal muscle Ca2+ handling in a novel DMD patient-derived pluripotent stem cell model

Mutations in the dystrophin gene cause the most common and currently incurable Duchenne muscular dystrophy (DMD) characterized by progressive muscle wasting. Although abnormal Ca²⁺ handling is a pathological feature of DMD, mechanisms underlying defective Ca²⁺ homeostasis remain unclear. Here we generate a novel DMD patient-derived pluripotent stem cell (PSC) model of skeletal muscle with an isogenic control using clustered regularly interspaced short palindromic repeat (CRISPR)-mediated precise gene correction. Transcriptome analysis identifies dysregulated gene sets in the absence of dystrophin, including genes involved in Ca²⁺ handling, excitation-contraction coupling and muscle contraction. Specifically, analysis of intracellular Ca²⁺ transients and mathematical modeling of Ca²⁺ dynamics reveal significantly reduced cytosolic Ca²⁺ clearance rates in DMD-PSC derived myotubes. Pharmacological assays demonstrate Ca²⁺ flux in myotubes is determined by both intracellular and extracellular sources. DMD-PSC derived myotubes display significantly reduced velocity of contractility. Compared with a non-isogenic wildtype PSC line, these pathophysiological defects could be rescued by CRISPR-mediated precise gene correction. Our study provides new insights into abnormal Ca²⁺ homeostasis in DMD and suggests that Ca²⁺ signaling pathways amenable to pharmacological modulation are potential therapeutic targets. Importantly, we have established a human physiology-relevant in vitro model enabling rapid pre-clinical testing of potential therapies for DMD.

Slowing Heart Rate Protects Against Pathological Cardiac Hypertrophy

We aimed to determine the pathophysiological impact of heart rate (HR) slowing on cardiac function. We have recently developed a murine model in which it is possible to conditionally delete the stimulatory heterotrimeric G-protein (Gα_s) in the sinoatrial (SA) node after the addition of tamoxifen using cre-loxP technology. The addition of tamoxifen leads to bradycardia. We used this approach to examine the physiological and pathophysiological effects of HR slowing. We first looked at the impact on exercise performance by running the mice on a treadmill. After the addition of tamoxifen, mice with conditional deletion of Gα_s in the SA node ran a shorter distance at a slower speed. Littermate controls preserved their exercise capacity after tamoxifen. Results consistent with impaired cardiac capacity in the mutants were also obtained with a dobutamine echocardiographic stress test. We then examined if HR reduction influenced pathological cardiac hypertrophy using two models: ligation of the left anterior descending coronary artery for myocardial infarction and abdominal aortic banding for hypertensive heart disease. In littermate controls, both procedures resulted in cardiac hypertrophy. However, induction of HR reduction prior to surgical intervention significantly ameliorated the hypertrophy. In order to assess potential protein kinase pathways that may be activated in the left ventricle by relative bradycardia, we used a phospho-antibody array and this revealed selective activation of phosphoinositide-3 kinase. In conclusion, HR reduction protects against pathological cardiac hypertrophy but limits physiological exercise capacity.

Prediction of Coronary Artery Disease and Major Adverse Cardiovascular Events Using Clinical and Genetic Risk Scores for Cardiovascular Risk Factors

BACKGROUND

Coronary artery disease (CAD) and major adverse cardiovascular events (MACE) are the leading causes of death in the general population, but risk stratification remains suboptimal. CAD genetic risk scores (GRSs) predict risk independently from clinical tools, like QRISK3. We assessed the added value of GRSs for a variety of cardiovascular traits (CV GRSs) for predicting CAD and MACE and tested their early-life screening potential by comparing against the CAD GRS only.

METHODS

We used data from 379 581 participants in the UK Biobank without known cardiovascular conditions (follow-up, 11.3 years; 3.3% CAD cases and 5.2% MACE cases). In a training subset (50%) we built 3 scores: QRISK3; QRISK3 and an established CAD GRS; and QRISK3, the CAD GRS and the CV GRSs. In an independent subset (50%), we evaluated each score's performance using the concordance index, odds ratio and net reclassification index. We then repeated the analyses without considering QRISK3.

RESULTS

For CAD, the combination of QRISK3 and the CAD GRS had a better performance than QRISK3 alone (concordance index, 0.766 versus 0.753; odds ratio, 5.47 versus 4.82; net reclassification index, 7.7%). Adding the CV GRSs did not significantly improve risk stratification. When only looking at genetic information, the combination of CV GRSs and the CAD GRS had a better performance than the CAD GRS alone (concordance index, 0.637 versus 0.625; odds ratio, 2.17 versus 2.07; net reclassification index, 3.3%). Similar results were obtained for MACE.

CONCLUSIONS

In individuals without known cardiovascular disease, the inclusion of CV GRSs to a clinical tool and an established CAD GRS does not improve CAD or MACE risk stratification. However, their combination only with the CAD GRS increases prediction performance indicating potential use in early-life screening before the advanced development of conventional cardiovascular risk factors.

ECG T-Wave Morphologic Variations Predict Ventricular Arrhythmic Risk in Low- and Moderate-Risk Populations

Background Early identification of individuals at risk of sudden cardiac death (SCD) remains a major challenge. The ECG is a simple, common test, with potential for large-scale application. We developed and tested the predictive value of a novel index quantifying T-wave morphologic variations with respect to a normal reference (TMV), which only requires one beat and a single-lead ECG. Methods and Results We obtained reference T-wave morphologies from 23 962 participants in the UK Biobank study. With Cox models, we determined the association between TMV and life-threatening ventricular arrhythmia in an independent data set from UK Biobank study without a history of cardiovascular events (N=51 794; median follow-up of 122 months) and SCD in patients with coronary artery disease from ARTEMIS (N=1872; median follow-up of 60 months). In UK Biobank study, 220 (0.4%) individuals developed life-threatening ventricular arrhythmias. TMV was significantly associated with life-threatening ventricular arrhythmias (hazard ratio [HR] of 1.13 per SD increase [95% CI, 1.03-1.24]; P=0.009). In ARTEMIS, 34 (1.8%) individuals reached the primary end point. Patients with TMV ≥5 had an HR for SCD of 2.86 (95% CI, 1.40-5.84; P=0.004) with respect to those with TMV <5, independently from QRS duration, corrected QT interval, and left ventricular ejection fraction. TMV was not significantly associated with death from a cause other than SCD. Conclusions TMV identifies individuals at life-threatening ventricular arrhythmia and SCD risk using a single-beat single-lead ECG, enabling inexpensive, quick, and safe risk assessment in large populations.

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.

Abstract P2005: Loss Of Full-length Mylk3 Causes Dilated Cardiomyopathy Via A MYL2-independent Mechanism

Introduction: We recently discovered C57BL/6N mice develop dilated cardiomyopathy (DCM), while the closely related C57BL/6J do not. We identified a variant (Chr8:85365179A>T, c.-5T>A) in myosin light chain kinase 3 ( Mylk3 ) in the C57BL/6N substrain that abolishes translation of the predominant transcript (ENSMUST00000034133). MYLK3 phosphorylates many sarcomere proteins, in particular myosin light chain 2 (MYL2), vital for cardiac development and function.

Hypothesis: The MYLK3 variant is sufficient to drive the development of DCM.

Results: We generated a novel isogenic mouse line carrying the Mylk3 c.-5T>A variant on the C57BL/6J background (C57BL/6J Mylk3.6N ). We found symptoms of DCM in male mutant mice at just 4 weeks old, with increased left ventricular diameter and volume, thinning of the posterior wall and a reduced ejection fraction. These data suggest the Mylk3 variant alone is sufficient to cause DCM. Surprisingly, we found no difference in phosphorylation of myosin light chain 2 (MYL2) by Phos-Tag® blotting at any age. This was unexpected, as previous reports in Mylk3 -/- mice found complete loss of MYL2 phosphorylation with a similar profile of DCM. However, in our mouse model, loss of MYLK3 is caused by a single point mutation in exon 1 which ablates the primary isoform (ENSMUST00000034133; Q3UIZ8) but western blotting with a C-terminal antibody revealed expression of a shorter isoform that is unaffected by the variant due to a distinct N-terminus (ENSMUST00000122452; D3Z630). There was no difference between the isoforms in their ability to phosphorylate MYL2 suggesting the shorter isoform can compensate. MYLK3 may also interact with ACTN2 and phosphorylate TNNI3. We found no difference in the expression or localisation of ACTN2 in our tissues, although co-IP in cells showed the primary isoform binds more readily to ACTN2, however, there was a significant reduction in TNNI3 S23/24 phosphorylation at all ages in C57BL/6J Mylk3.6N mice. This could affect cross-bridge cycling kinetics by altering the function of TNNI3 within the tropomyosin/troponin complex.

Conclusions: Loss of full-length MYLK3 causes DCM via a MYL2-independent mechanism, possibly due to lack of phosphorylation of alternative target(s), such as TNNI3.

Whole genome and transcriptome analysis enhances precision cancer treatment options

BACKGROUND

Recent advances are enabling delivery of precision genomic medicine to cancer clinics. While the majority of approaches profile panels of selected genes or hotspot regions, comprehensive data provided by whole genome and transcriptome sequencing and analysis (WGTA) presents an opportunity to align a much larger proportion of patients to therapies.

PATIENTS AND METHODS

Samples from 570 patients with advanced or metastatic cancer of diverse types enrolled in the Personalized OncoGenomics (POG) program underwent WGTA. DNA-based data, including mutations, copy number, and mutation signatures, were combined with RNA-based data, including gene expression and fusions, to generate comprehensive WGTA profiles. A multidisciplinary molecular tumour board used WGTA profiles to identify and prioritize clinically actionable alterations and inform therapy. Patient responses to WGTA-informed therapies were collected.

RESULTS

Clinically actionable targets were identified for 83% of patients, 37% of whom received WGTA-informed treatments. RNA expression data were particularly informative, contributing to 67% of WGTA-informed treatments; 25% of treatments were informed by RNA expression alone. Of a total 248 WGTA-informed treatments, 46% resulted in clinical benefit. RNA expression data were comparable to DNA-based mutation and copy number data in aligning to clinically beneficial treatments. Genome signatures also guided therapeutics including platinum, PARP inhibitors, and immunotherapies. Patients accessed WGTA-informed treatments through clinical trials (19%), off-label use (35%), and as standard therapies (46%) including those which would not otherwise have been the next choice of therapy, demonstrating the utility of genomic information to direct use of chemotherapies as well as targeted therapies.

CONCLUSIONS

Integrating RNA expression and genome data illuminated treatment options that resulted in 46% of treated patients experiencing positive clinical benefit, supporting the use of comprehensive WGTA profiling in clinical cancer care.

CLINICAL TRIAL NUMBER

NCT02155621.

Smooth muscle ATP-sensitive potassium channels mediate migraine-relevant hypersensitivity in mouse models

BACKGROUND

Opening of KATP channels by systemic levcromakalim treatment triggers attacks in migraine patients and hypersensitivity to von Frey stimulation in a mouse model. Blocking of these channels is effective in several preclinical migraine models. It is unknown in what tissue and cell type KATP-induced migraine attacks are initiated and which KATP channel subtype is targeted.

METHODS

In mouse models, we administered levcromakalim intracerebroventricularly, intraperitoneally and intraplantarily and compared the nociceptive responses by von Frey and hotplate tests. Mice with a conditional loss-of-function mutation in the smooth muscle KATP channel subunit Kir6.1 were given levcromakalim and GTN and examined with von Frey filaments. Arteries were tested for their ability to dilate ex vivo. mRNA expression, western blotting and immunohistochemical stainings were made to identify relevant target tissue for migraine induced by KATP channel opening.

RESULTS

Systemic administration of levcromakalim induced hypersensitivity but central and local administration provided antinociception respectively no effect. The Kir6.1 smooth muscle knockout mouse was protected from both GTN and levcromakalim induced hypersensitivity, and their arteries had impaired dilatory response to the latter. mRNA and protein expression studies showed that trigeminal ganglia did not have significant KATP channel expression of any subtype, whereas brain arteries and dura mater primarily expressed the Kir6.1 + SUR2B subtype.

CONCLUSION

Hypersensitivity provoked by GTN and levcromakalim in mice is dependent on functional smooth muscle KATP channels of extracerebral origin. These results suggest a vascular contribution to hypersensitivity induced by migraine triggers.

Genomic and pleiotropic analyses of resting QT interval identifies novel loci and overlap with atrial electrical disorders

The resting QT interval, an electrocardiographic (ECG) measure of ventricular myocardial repolarization, is a heritable risk marker of cardiovascular mortality, but the mechanisms remain incompletely understood. Previously reported candidate genes have provided insights into the regulatory mechanisms of the QT interval. However, there are still important knowledge gaps. We aimed to gain new insights by (i) providing new candidate genes, (ii) identifying pleiotropic associations with other cardiovascular traits, and (iii) scanning for sexually dimorphic genetic effects. We conducted a genome-wide association analysis for resting QT interval with ~9.8 million variants in 52 107 individuals of European ancestry without known cardiovascular disease from the UK Biobank. We identified 40 loci, 13 of which were novel, including 2 potential sex-specific loci, explaining ~11% of the trait variance. Candidate genes at novel loci were involved in myocardial structure and arrhythmogenic cardiomyopathy. Investigation of pleiotropic effects of QT interval variants using phenome-wide association analyses in 302 000 unrelated individuals from the UK Biobank and pairwise genome-wide comparisons with other ECG and cardiac imaging traits revealed genetic overlap with atrial electrical pathology. These findings provide novel insights into how abnormal myocardial repolarization and increased cardiovascular mortality may be linked.

Early inflammation precedes cardiac fibrosis and heart failure in desmoglein 2 murine model of arrhythmogenic cardiomyopathy

The study of a desmoglein 2 murine model of arrhythmogenic cardiomyopathy revealed cardiac inflammation as a key early event leading to fibrosis. Arrhythmogenic cardiomyopathy (AC) is an inherited heart muscle disorder leading to ventricular arrhythmias and heart failure due to abnormalities in the cardiac desmosome. We examined how loss of desmoglein 2 (Dsg2) in the young murine heart leads to development of AC. Apoptosis was an early cellular phenotype, and RNA sequencing analysis revealed early activation of inflammatory-associated pathways in Dsg2-null (Dsg2^-/-) hearts at postnatal day 14 (2 weeks) that were absent in the fibrotic heart of adult mice (10 weeks). This included upregulation of iRhom2/ADAM17 and its associated pro-inflammatory cytokines and receptors such as TNFα, IL6R and IL-6. Furthermore, genes linked to specific macrophage populations were also upregulated. This suggests cardiomyocyte stress triggers an early immune response to clear apoptotic cells allowing tissue remodelling later on in the fibrotic heart. Our analysis at the early disease stage suggests cardiac inflammation is an important response and may be one of the mechanisms responsible for AC disease progression.

Genetic variants in TRPM7 associated with unexplained stillbirth modify ion channel function

Stillbirth is the loss of a fetus after 22 weeks of gestation, of which almost half go completely unexplained despite post-mortem. We recently sequenced 35 arrhythmia-associated genes from 70 unexplained stillbirth cases. Our hypothesis was that deleterious mutations in channelopathy genes may have a functional effect in utero that may be pro-arrhythmic in the developing fetus. We observed four heterozygous, nonsynonymous variants in transient receptor potential melastatin 7 (TRPM7), a ubiquitously expressed ion channel known to regulate cardiac development and repolarization in mice.

We used site-directed mutagenesis and single-cell patch-clamp to analyze the functional effect of the four stillbirth mutants on TRPM7 ion channel function in heterologous cells. We also used cardiomyocytes derived from human pluripotent stem cells to model the contribution of TRPM7 to action potential morphology.

Our results show that two TRPM7 variants, p.G179V and p.T860M, lead to a marked reduction in ion channel conductance. This observation was underpinned by a lack of measurable TRPM7 protein expression, which in the case of p.T860M was due to rapid proteasomal degradation. We also report that human hiPSC-derived cardiomyocytes possess measurable TRPM7 currents; however, siRNA knockdown did not directly affect action potential morphology.

TRPM7 variants found in the unexplained stillbirth population adversely affect ion channel function and this may precipitate fatal arrhythmia in utero.

Analysing electrocardiographic traits and predicting cardiac risk in UK biobank

The electrocardiogram (ECG) is a commonly used clinical tool that reflects cardiac excitability and disease. Many parameters are can be measured and with the improvement of methodology can now be quantified in an automated fashion, with accuracy and at scale. Furthermore, these measurements can be heritable and thus genome wide association studies inform the underpinning biological mechanisms. In this review we describe how we have used the resources in UK Biobank to undertake such work. In particular, we focus on a substudy uniquely describing the response to exercise performed at scale with accompanying genetic information.

Genetically Determined Serum Calcium Levels and Markers of Ventricular Repolarization: A Mendelian Randomization Study in the UK Biobank

BACKGROUND

ECG markers of ventricular depolarization and repolarization are associated with an increased risk of arrhythmia and sudden cardiac death. Our prior work indicated lower serum calcium concentrations are associated with longer QT and JT intervals in the general population. Here, we investigate whether serum calcium is a causal risk factor for changes in ECG measures using Mendelian randomization (MR).

METHODS

Independent lead variants from a newly performed genome-wide association study for serum calcium in >300 000 European-ancestry participants from UK Biobank were used as instrumental variables. Two-sample MR analyses were performed to approximate the causal effect of serum calcium on QT, JT, and QRS intervals using an inverse-weighted method in 76 226 participants not contributing to the serum calcium genome-wide association study. Sensitivity analyses including MR-Egger, weighted-median estimator, and MR pleiotropy residual sum and outlier were performed to test for the presence of horizontal pleiotropy.

RESULTS

Two hundred five independent lead calcium-associated variants were used as instrumental variables for MR. A decrease of 0.1 mmol/L serum calcium was associated with longer QT (3.01 ms [95% CI, 2.03 to 3.99]) and JT (2.89 ms [1.91 to 3.87]) intervals. A weak association was observed for QRS duration (secondary analyses only). Results were concordant in all sensitivity analyses.

CONCLUSIONS

These analyses support a causal effect of serum calcium levels on ventricular repolarization, in a middle-aged population of European-ancestry where serum calcium concentrations are likely stable and chronic. Modulation of calcium concentration may, therefore, directly influence cardiovascular disease risk.

Enhancing mutant IKS channel activity by increasing endogenous PIP2 levels and its interaction with PKA signalling pathway

Funding Acknowledgements

Type of funding sources: Public grant(s) – National budget only. Main funding source(s): British Heart Foundation

Phosphatidylinositol-4,5-biphosphate (PIP2) is implicated in the regulation and modulation of the IKS channel. The channel is formed at the plasma membrane by the co-assembly of KCNQ1 and KCNE1. Patients with Congenital Long QT 1(LQT1) syndrome are predisposed to Polymorphic VT due to mutations in KCNQ1, leading to impaired channel activity.

We initially transfect Human Embryonic Kidney (HEK) cells with a mammalian vector expressing KCNQ1 gene tagged with green fluorescent protein, along with KCNE1 to form the wild type (WT) IKS channel. The cells were also transfected with a constitutively active PI(4)P 5-kinase(PIP5K), which converts the phospholipid Phosphatidylinositol-4-phosphate to PIP2, therefore increasing endogenous levels of PIP2. To ensure the enzyme remains localised at the plasma membrane we attached it to CFP-FKBP and we co-transfected the cells with Lyn11-FRB construct that tethers to the plasma membrane. When these cells were perfused with Rapamycin it induced chemical dimerization of CF-PIP5K to lyn11. We utilised an inactive PIP5K as a control. Mutants were created with a site directed mutagenesis kit.

In the presence CF-PIP5K, whole cell voltage clamp recordings demonstrated a 2.5 fold statistically significant increase in WT channel activity (at +80mV,p &lt; 0.001), when compared to unaltered PIP2 conditions. Heterozygous Serine566phe and Phe340del mutants had statistically significant reduction in current density compared to wild type in basal conditions. When these mutants were expressed with the active CF-PIP5K, Serine566phe and Phe340 had a 2.97 and 3.30 fold increase in current density, respectively (p &lt;0.05). Homozygous Mutants D242N and T247in also showed statistically significant channel activity.

We substituted serine with alanine at site 27 and 92(S27A/S92A) to generate a mutant known to disrupt cAMP mediated upregulation, there was a statistical 3.3 fold (80 + mV) increase in current density when co-expressed with CF-PIP5K. We then substituted serine with aspartic acid (S27D/S92D) to create a Phosphomimetic mutation, this mutant reproduces the effects of sympathetic mediated augmentation of IKS channel. In the presence of enhanced PIP2 levels, the S27D/S92D failed to demonstrate a statistical increase in current, implying the channel is at its maximum activity and hence we failed to observe any further modulation.

We then proceeded to interrogate how PIP2 interacts with sympathetic signalling system. Pseudojanin(PJ) causes depletion of PIP2 hence perturbing channel activity. When PJ was expressed with KCNQ1 and KCNE1 we observed an 80% reduction in channel activity at +80mV(P &lt;0.001). When we perfused these cells with isoprenaline the channel activity was restored to normal.

Here we illustrate how increasing PIP2 levels can revive IKS channel activity in mutant genotype, therefore supporting evidence of its capabilities as a potential therapeutic tool. This modulation is independent of the PKA-cAMP system. Abstract Figure. Current Increment

The role of resistance to inhibitors of cholinesterase 8b in the control of heart rate

We have assessed the role of ric-b8 in the control of heart rate after the gene was implicated in a recent genome-wide association study of resting heart rate. We developed a novel murine model in which it was possible to conditionally delete ric-8b in the sinoatrial (SA) node after the addition of tamoxifen. Despite this, we were unable to obtain homozygotes and thus studied heterozygotes. Haploinsufficiency of ric-8b in the sinoatrial node induced by the addition of tamoxifen in adult animals leads to mice with a reduced heart rate. However, other electrocardiographic intervals (e.g., PR and QRS) were normal, and there was no apparent arrhythmia such as heart block. The positive chronotropic response to isoprenaline was abrogated, whereas the response to carbachol was unchanged. The pacemaker current I_f (funny current) has an important role in regulating heart rate, and its function is modulated by both isoprenaline and carbachol. Using a heterologous system expressing HCN4, we show that ric-8b can modulate the HCN4 current. Overexpression of ric-8b led to larger HCN4 currents, whereas silencing ric-8b led to smaller currents. Ric-8b modulates heart rate responses in vivo likely via its actions on the stimulatory G-protein.

No Clinically Relevant Effect of Heart Rate Increase and Heart Rate Recovery During Exercise on Cardiovascular Disease: A Mendelian Randomization Analysis

BACKGROUND

Reduced heart rate (HR) increase (HRI), recovery (HRR), and higher resting HR are associated with cardiovascular (CV) disease, but causal inferences have not been deduced. We investigated causal effects of HRI, HRR, and resting HR on CV risk, all-cause mortality (ACM), atrial fibrillation (AF), coronary artery disease (CAD), and ischemic stroke (IS) using Mendelian Randomization.

METHODS

11 variants for HRI, 11 for HRR, and two sets of 46 and 414 variants for resting HR were obtained from four genome-wide association studies (GWASs) on UK Biobank. We performed a lookup on GWASs for CV risk and ACM in UK Biobank (N = 375,367, 5.4% cases and N = 393,165, 4.4% cases, respectively). For CAD, AF, and IS, we used publicly available summary statistics. We used a random-effects inverse-variance weighted (IVW) method and sensitivity analyses to estimate causality.

RESULTS

IVW showed a nominally significant effect of HRI on CV events (odds ratio [OR] = 1.0012, P = 4.11 × 10-2) and on CAD and AF. Regarding HRR, IVW was not significant for any outcome. The IVW method indicated statistically significant associations of resting HR with AF (OR = 0.9825, P = 9.8 × 10-6), supported by all sensitivity analyses, and a nominally significant association with IS (OR = 0.9926, P = 9.82 × 10-3).

CONCLUSION

Our findings suggest no strong evidence of an association between HRI and HRR and any outcome and confirm prior work reporting a highly significant effect of resting HR on AF. Future research is required to explore HRI and HRR associations further using more powerful predictors, when available.

A Method to Minimise the Impact of ECG Marker Inaccuracies on the Spatial QRS-T angle: Evaluation on 1,512 Manually Annotated ECGs

•

Inaccuracies of QRS and T-wave markers significantly impact QRS-Ta estimation.

•

These errors influence the classification of clinically relevant abnormal values.

•

Our algorithm provides robust measurements in the presence of inaccurate VCG markers.

•

We present for the first time, the distribution of the QRS-Ta in a large cohort.

The spatial QRS-T angle (QRS-Ta) derived from the vectorcardiogram (VCG) is a strong risk predictor for ventricular arrhythmia and sudden cardiac death with potential use for mass screening. Accurate QRS-Ta estimation in the presence of ECG delineation errors is crucial for its deployment as a prognostic test. Our study assessed the effect of inaccurate QRS and T-wave marker placement on QRS-Ta estimation and proposes a robust method for its calculation.

Reference QRS-Ta measurements were derived from 1,512 VCGs manually annotated by three expert reviewers. We systematically changed onset and offset timings of QRS and T-wave markers to simulate inaccurate placement. The QRS-Ta was recalculated using a standard approach and our proposed algorithm, which limits the impact of VCG marker inaccuracies by defining the vector origin as an interval preceding QRS-onset and redefines the beginning and end of QRS and T-wave loops.

Using the standard approach, mean absolute errors (MAE) in peak QRS-Ta were >40% and sensitivity and precision in the detection of abnormality (>105°) were <80% and <65% respectively, when QRS-onset was delayed or QRS-offset anticipated >15 ms. Using our proposed algorithm, MAE for peak QRS-Ta were reduced to <4% and sensitivity and precision of abnormality were >94% for inaccuracies up to ±15 ms. Similar results were obtained for mean QRS-Ta.

In conclusion, inaccuracies of QRS and T-wave markers can significantly influence the QRS-Ta. Our proposed algorithm provides robust QRS-Ta measurements in the presence of inaccurate VCG annotation, enabling its use in large datasets.

Cardiac phenotype in ATP1A3-related syndromes: A multicenter cohort study

OBJECTIVE

To define the risks and consequences of cardiac abnormalities in ATP1A3-related syndromes.

METHODS

Patients meeting clinical diagnostic criteria for rapid-onset dystonia-parkinsonism (RDP), alternating hemiplegia of childhood (AHC), and cerebellar ataxia, areflexia, pes cavus, optic atrophy, and sensorineural hearing loss (CAPOS) with ATP1A3 genetic analysis and at least 1 cardiac assessment were included. We evaluated the cardiac phenotype in an Atp1a3 knock-in mouse (Mashl^+/-) to determine the sequence of events in seizure-related cardiac death.

RESULTS

Ninety-eight patients with AHC, 9 with RDP, and 3 with CAPOS (63 female, mean age 17 years) were included. Resting ECG abnormalities were found in 52 of 87 (60%) with AHC, 2 of 3 (67%) with CAPOS, and 6 of 9 (67%) with RDP. Serial ECGs showed dynamic changes in 10 of 18 patients with AHC. The first Holter ECG was abnormal in 24 of 65 (37%) cases with AHC and RDP with either repolarization or conduction abnormalities. Echocardiography was normal. Cardiac intervention was required in 3 of 98 (≈3%) patients with AHC. In the mouse model, resting ECGs showed intracardiac conduction delay; during induced seizures, heart block or complete sinus arrest led to death.

CONCLUSIONS

We found increased prevalence of ECG dynamic abnormalities in all ATP1A3-related syndromes, with a risk of life-threatening cardiac rhythm abnormalities equivalent to that in established cardiac channelopathies (≈3%). Sudden cardiac death due to conduction abnormality emerged as a seizure-related outcome in murine Atp1a3-related disease. ATP1A3-related syndromes are cardiac diseases and neurologic diseases. We provide guidance to identify patients potentially at higher risk of sudden cardiac death who may benefit from insertion of a pacemaker or implantable cardioverter-defibrillator.

Endothelial ATP-Sensitive Potassium Channel Protects Against the Development of Hypertension and Atherosclerosis

In the endothelium, ATP-sensitive potassium (KATP) channels are thought to couple cellular metabolism with membrane excitability, calcium entry, and endothelial mediator release. We hypothesized that endothelial KATP channels have a broad role protecting against high blood pressure and atherosclerosis. Endothelial-specific Kir6.1 KO mice (eKO) and eKO mice on an apolipoprotein E KO background were generated (A-eKO) to investigate the role of KATP channels in the endothelium. Basal blood pressure was not elevated in eKO mice. However, when challenged with a high-salt diet and the eNOS inhibitor L-NAME, eKO mice became more hypertensive than their littermate controls. In aorta, NO release at least partly contributes to the endothelium-dependent vasorelaxation induced by pinacidil. In A-eKO mice atherosclerotic plaque density was significantly greater than in their littermate controls when challenged with a high-fat diet, particularly in the aortic arch region. Levels of endothelial dysfunction markers were higher in eKO compared with WT mice; however, these were not significant for A-eKO mice compared with their littermate controls. Furthermore, decreased vascular reactivity was observed in the mesenteric arteries of A-eKO mice, but not in aorta when on a high-fat diet. Our data support a role for endothelial Kir6.1-containing KATP channels in the endothelial protection against environmental stressors: the maintenance of blood pressure homeostasis in response to high salt and endothelial integrity when challenged with a high-fat diet.

The contribution of non-coding regulatory elements to cardiovascular disease

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

Genetic Basis and Prognostic Value of Exercise QT Dynamics

Supplemental Digital Content is available in the text.

Background:

Abnormal QT interval responses to heart rate (QT dynamics) is an independent risk predictor for cardiovascular disease in patients, but its genetic basis and prognostic value in a population-based cohort have not been investigated.

Methods:

QT dynamics during exercise and recovery were derived in 56 643 individuals from UK Biobank without a history of cardiovascular events. Genome-wide association studies were conducted to identify genetic variants and bioinformatics analyses were performed to prioritize candidate genes. The prognostic value of QT dynamics was evaluated for cardiovascular events (death or hospitalization) and all-cause mortality.

Results:

Heritability of QT dynamics during exercise and recovery were 10.7% and 5.4%, respectively. Genome-wide association studies identified 20 loci, of which 4 loci included genes implicated in mendelian long-QT syndrome. Five loci did not overlap with previously reported resting QT interval loci; candidate genes included KCNQ4 and KIAA1755. Genetic risk scores were not associated with cardiovascular events in 357 882 unrelated individuals from UK Biobank. We also did not observe associations of QT dynamics during exercise and recovery with cardiovascular events. Increased QT dynamics during recovery was significantly associated with all-cause mortality in the univariate Cox regression analysis (hazard ratio, 1.09 [95% CI, 1.05–1.13], P=2.28×10^-5), but the association was not significant after adjusting for clinical risk factors.

Conclusions:

QT interval dynamics during exercise and recovery are heritable markers but do not carry independent prognostic information for clinical outcomes in the UK Biobank, a population-based cohort. Their prognostic importance may relate to cardiovascular disease cohorts where structural heart disease or ischemia may influence repolarization dynamics. The strong overlap between QT dynamics and resting QT interval loci suggests common biological pathways; however, nonoverlapping loci suggests alternative mechanisms may exist that underlie QT interval dynamics.

Common Genetic Variants Modulate the Electrocardiographic Tpeak-to-Tend Interval

Sudden cardiac death is responsible for half of all deaths from cardiovascular disease. The analysis of the electrophysiological substrate for arrhythmias is crucial for optimal risk stratification. A prolonged T-peak-to-Tend (Tpe) interval on the electrocardiogram is an independent predictor of increased arrhythmic risk, and Tpe changes with heart rate are even stronger predictors. However, our understanding of the electrophysiological mechanisms supporting these risk factors is limited. We conducted genome-wide association studies (GWASs) for resting Tpe and Tpe response to exercise and recovery in ∼30,000 individuals, followed by replication in independent samples (∼42,000 for resting Tpe and ∼22,000 for Tpe response to exercise and recovery), all from UK Biobank. Fifteen and one single-nucleotide variants for resting Tpe and Tpe response to exercise, respectively, were formally replicated. In a full dataset GWAS, 13 further loci for resting Tpe, 1 for Tpe response to exercise and 1 for Tpe response to exercise were genome-wide significant (p ≤ 5 × 10^-8). Sex-specific analyses indicated seven additional loci. In total, we identify 32 loci for resting Tpe, 3 for Tpe response to exercise and 3 for Tpe response to recovery modulating ventricular repolarization, as well as cardiac conduction and contraction. Our findings shed light on the genetic basis of resting Tpe and Tpe response to exercise and recovery, unveiling plausible candidate genes and biological mechanisms underlying ventricular excitability.

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.

Mylk3 null C57BL/6N mice develop cardiomyopathy, whereas Nnt null C57BL/6J mice do not

The C57BL/6J and C57BL/6N mice have well-documented phenotypic and genotypic differences, including the infamous nicotinamide nucleotide transhydrogenase (Nnt) null mutation in the C57BL/6J substrain, which has been linked to cardiovascular traits in mice and cardiomyopathy in humans. To assess whether Nnt loss alone causes a cardiovascular phenotype, we investigated the C57BL/6N, C57BL/6J mice and a C57BL/6J-BAC transgenic rescuing NNT expression, at 3, 12, and 18 mo. We identified a modest dilated cardiomyopathy in the C57BL/6N mice, absent in the two B6J substrains. Immunofluorescent staining of cardiomyocytes revealed eccentric hypertrophy in these mice, with defects in sarcomere organisation. RNAseq analysis identified differential expression of a number of cardiac remodelling genes commonly associated with cardiac disease segregating with the phenotype. Variant calling from RNAseq data identified a myosin light chain kinase 3 (Mylk3) mutation in C57BL/6N mice, which abolishes MYLK3 protein expression. These results indicate the C57BL/6J Nnt-null mice do not develop cardiomyopathy; however, we identified a null mutation in Mylk3 as a credible cause of the cardiomyopathy phenotype in the C57BL/6N.

Whole-cell and Perforated Patch-clamp Recordings from Acutely-isolated Murine Sino-atrial Node Cells

Cardiac pacemaker cells of the sino-atrial node are responsible for the initiation of the heart beat and express an array of ion channels. The patch-clamp technique is the gold standard method for investigating the function of ion channels expressed in electrically active cells. Conventional whole-cell and perforated patch-clamp techniques can be used to investigate ionic currents in the voltage-clamp mode and changes in membrane potential (e.g., action potential) in the current-clamp mode. Here, we provide details of protocols used to measure spontaneous and triggered action potentials and whole-cell funny current If (HCN4) in single cardiomyocytes isolated from the mouse sino-atrial node (SAN).

Acute Isolation of Cells from Murine Sino-atrial Node

The cardiac conduction system allows the synchronized propagation of electrical activity through heart muscle. This is initiated by the spontaneous activity of the specialized pacemaker cells of the sino-atrial node (SAN). The SAN region underlies automaticity in mammals and therefore has a crucial role in the pathogenesis of cardiac disorders such as arrhythmia. Isolation of SAN tissue and SAN cells is critical to advance our understanding of SAN structure and function in health and disease. Initially, isolation of SAN tissue and SAN cells was carried out in the rabbit owing to its larger size and similar electrical properties to human. This protocol was optimized by Mangoni and Nargeot (2001) for use in mice to take advantage of advancements in transgenic models. Here, we provide a step-by-step guide to dissecting the SAN tissue and isolating pacemaker cardiomyocytes from mouse hearts using an enzyme digestion approach.

CRISPR/Cas9 editing in human pluripotent stem cell-cardiomyocytes highlights arrhythmias, hypocontractility, and energy depletion as potential therapeutic targets for hypertrophic cardiomyopathy

Aims

Sarcomeric gene mutations frequently underlie hypertrophic cardiomyopathy (HCM), a prevalent and complex condition leading to left ventricle thickening and heart dysfunction. We evaluated isogenic genome-edited human pluripotent stem cell-cardiomyocytes (hPSC-CM) for their validity to model, and add clarity to, HCM.

Methods and results

CRISPR/Cas9 editing produced 11 variants of the HCM-causing mutation c.C9123T-MYH7 [(p.R453C-β-myosin heavy chain (MHC)] in 3 independent hPSC lines. Isogenic sets were differentiated to hPSC-CMs for high-throughput, non-subjective molecular and functional assessment using 12 approaches in 2D monolayers and/or 3D engineered heart tissues. Although immature, edited hPSC-CMs exhibited the main hallmarks of HCM (hypertrophy, multi-nucleation, hypertrophic marker expression, sarcomeric disarray). Functional evaluation supported the energy depletion model due to higher metabolic respiration activity, accompanied by abnormalities in calcium handling, arrhythmias, and contraction force. Partial phenotypic rescue was achieved with ranolazine but not omecamtiv mecarbil, while RNAseq highlighted potentially novel molecular targets.

Conclusion

Our holistic and comprehensive approach showed that energy depletion affected core cardiomyocyte functionality. The engineered R453C-βMHC-mutation triggered compensatory responses in hPSC-CMs, causing increased ATP production and αMHC to energy-efficient βMHC switching. We showed that pharmacological rescue of arrhythmias was possible, while MHY7: MYH6 and mutant: wild-type MYH7 ratios may be diagnostic, and previously undescribed lncRNAs and gene modifiers are suggestive of new mechanisms.

Cardiovascular Predictive Value and Genetic Basis of Ventricular Repolarization Dynamics

BACKGROUND

Early prediction of cardiovascular risk in the general population remains an important issue. The T-wave morphology restitution (TMR), an ECG marker quantifying ventricular repolarization dynamics, is strongly associated with cardiovascular mortality in patients with heart failure. Our aim was to evaluate the cardiovascular prognostic value of TMR in a UK middle-aged population and identify any genetic contribution.

METHODS

We analyzed ECG recordings from 55 222 individuals from a UK middle-aged population undergoing an exercise stress test in UK Biobank (UKB). TMR was used to measure ventricular repolarization dynamics, exposed in this cohort by exercise (TMR during exercise, TMR^ex) and recovery from exercise (TMR during recovery, TMR^rec). The primary end point was cardiovascular events; secondary end points were all-cause mortality, ventricular arrhythmias, and atrial fibrillation with median follow-up of 7 years. Genome-wide association studies for TMR^ex and TMR^rec were performed, and genetic risk scores were derived and tested for association in independent samples from the full UKB cohort (N=360 631).

RESULTS

A total of 1743 (3.2%) individuals in UKB who underwent the exercise stress test had a cardiovascular event, and TMR^rec was significantly associated with cardiovascular events (hazard ratio, 1.11; P=5×10^-7), independent of clinical variables and other ECG markers. TMR^rec was also associated with all-cause mortality (hazard ratio, 1.10) and ventricular arrhythmias (hazard ratio, 1.16). We identified 12 genetic loci in total for TMR^ex and TMR^rec, of which 9 are associated with another ECG marker. Individuals in the top 20% of the TMR^rec genetic risk score were significantly more likely to have a cardiovascular event in the full UKB cohort (18 997, 5.3%) than individuals in the bottom 20% (hazard ratio, 1.07; P=6×10^-3).

CONCLUSIONS

TMR and TMR genetic risk scores are significantly associated with cardiovascular risk in a UK middle-aged population, supporting the hypothesis that increased spatio-temporal heterogeneity of ventricular repolarization is a substrate for cardiovascular risk and the validity of TMR as a cardiovascular risk predictor.

A critical role for the ATP-sensitive potassium channel subunit KIR6.1 in the control of cerebral blood flow

KIR6.1 (KCNJ8) is a subunit of ATP sensitive potassium channel (KATP) that plays an important role in the control of peripheral vascular tone and is highly expressed in brain contractile cells (vascular smooth muscle cells and pericytes). This study determined the effect of global deletion of the KIR6.1 subunit on cerebral blood flow, neurovascular coupling and cerebral oxygenation in mice. In KIR6.1 deficient mice resting cerebral blood flow and brain parenchymal partial pressure of oxygen (PO2) were found to be markedly lower compared to that in their wildtype littermates. However, cortical blood oxygen level dependent responses triggered by visual stimuli were not affected in conditions of KIR6.1 deficiency. These data suggest that KATP channels containing KIR6.1 subunit are critically important for the maintenance of normal cerebral perfusion and parenchymal PO2 but play no significant role in the mechanisms underlying functional changes in brain blood flow.