Pilot study to evaluate left-to-right ventricular offset in biventricular pacing-comparison of electrocardiographic imaging and ECG.

INTRODUCTION: Biventricular pacing (BiVp) improves outcomes in systolic heart failure patients with electrical dyssynchrony. BiVp is delivered from epicardial left ventricular (LV) and endocardial right ventricular (RV) electrodes. Acute electrical activation changes with different LV-RV stimulation offsets can help guide individually optimized BiVp programming. We sought to study the BiVp ventricular activation with different LV-RV offsets and compare with 12-lead ECG. METHODS: In five patients with BiVp (63 ± 17-year-old, 80% male, LV ejection fraction 27 ± 6%), we evaluated acute ventricular epicardial activation, varying LV-RV offsets in 20 ms increments from -40 to 80 ms, using electrocardiographic imaging (ECGI) to obtain absolute ventricular electrical uncoupling (VEUabs, absolute difference in average LV and average RV activation time) and total activation time (TAT). For each patient, we calculated the correlation between ECGI and corresponding ECG (3D-QRS-area and QRS duration) with different LV-RV offsets. RESULTS: The LV-RV offset to attain minimum VEUabs in individual patients ranged 20-60 ms. In all patients, a larger LV-RV offset was required to achieve minimum VEUabs (36 ± 17 ms) or 3D-QRS-area (40 ± 14 ms) than that for minimum TAT (-4 ± 9 ms) or QRS duration (-8 ± 11 ms). In individual patients, 3D-QRS-area correlated with VEUabs (r 0.65 ± 0.24) and QRS duration correlated with TAT (r 0.95 ± 0.02). Minimum VEUabs and minimum 3D-QRS-area were obtained by LV-RV offset within 20 ms of each other in all five patients. CONCLUSIONS: LV-RV electrical uncoupling, as assessed by ECGI, can be minimized by optimizing LV-RV stimulation offset. 3D-QRS-area is a surrogate to identify LV-RV offset that minimizes LV-RV uncoupling.

Technical development and feasibility of a reusable vest to integrate cardiovascular magnetic resonance with electrocardiographic imaging.

BACKGROUND: Electrocardiographic imaging (ECGI) generates electrophysiological (EP) biomarkers while cardiovascular magnetic resonance (CMR) imaging provides data about myocardial structure, function and tissue substrate. Combining this information in one examination is desirable but requires an affordable, reusable, and high-throughput solution. We therefore developed the CMR-ECGI vest and carried out this technical development study to assess its feasibility and repeatability in vivo. METHODS: CMR was prospectively performed at 3T on participants after collecting surface potentials using the locally designed and fabricated 256-lead ECGI vest. Epicardial maps were reconstructed to generate local EP parameters such as activation time (AT), repolarization time (RT) and activation recovery intervals (ARI). 20 intra- and inter-observer and 8 scan re-scan repeatability tests. RESULTS: 77 participants were recruited: 27 young healthy volunteers (HV, 38.9 ± 8.5 years, 35% male) and 50 older persons (77.0 ± 0.1 years, 52% male). CMR-ECGI was achieved in all participants using the same reusable, washable vest without complications. Intra- and inter-observer variability was low (correlation coefficients [rs] across unipolar electrograms = 0.99 and 0.98 respectively) and scan re-scan repeatability was high (rs between 0.81 and 0.93). Compared to young HV, older persons had significantly longer RT (296.8 vs 289.3 ms, p = 0.002), ARI (249.8 vs 235.1 ms, p = 0.002) and local gradients of AT, RT and ARI (0.40 vs 0.34 ms/mm, p = 0,01; 0.92 vs 0.77 ms/mm, p = 0.03; and 1.12 vs 0.92 ms/mm, p = 0.01 respectively). CONCLUSION: Our high-throughput CMR-ECGI solution is feasible and shows good reproducibility in younger and older participants. This new technology is now scalable for high throughput research to provide novel insights into arrhythmogenesis and potentially pave the way for more personalised risk stratification. CLINICAL TRIAL REGISTRATION: Title: Multimorbidity Life-Course Approach to Myocardial Health-A Cardiac Sub-Study of the MRC National Survey of Health and Development (NSHD) (MyoFit46). National Clinical Trials (NCT) number: NCT05455125. URL: https://clinicaltrials.gov/ct2/show/NCT05455125?term=MyoFit&draw=2&rank=1.

Study protocol: MyoFit46-the cardiac sub-study of the MRC National Survey of Health and Development.

BACKGROUND: The life course accumulation of overt and subclinical myocardial dysfunction contributes to older age mortality, frailty, disability and loss of independence. The Medical Research Council National Survey of Health and Development (NSHD) is the world's longest running continued surveillance birth cohort providing a unique opportunity to understand life course determinants of myocardial dysfunction as part of MyoFit46-the cardiac sub-study of the NSHD. METHODS: We aim to recruit 550 NSHD participants of approximately 75 years+ to undertake high-density surface electrocardiographic imaging (ECGI) and stress perfusion cardiovascular magnetic resonance (CMR). Through comprehensive myocardial tissue characterization and 4-dimensional flow we hope to better understand the burden of clinical and subclinical cardiovascular disease. Supercomputers will be used to combine the multi-scale ECGI and CMR datasets per participant. Rarely available, prospectively collected whole-of-life data on exposures, traditional risk factors and multimorbidity will be studied to identify risk trajectories, critical change periods, mediators and cumulative impacts on the myocardium. DISCUSSION: By combining well curated, prospectively acquired longitudinal data of the NSHD with novel CMR-ECGI data and sharing these results and associated pipelines with the CMR community, MyoFit46 seeks to transform our understanding of how early, mid and later-life risk factor trajectories interact to determine the state of cardiovascular health in older age. TRIAL REGISTRATION: Prospectively registered on ClinicalTrials.gov with trial ID: 19/LO/1774 Multimorbidity Life-Course Approach to Myocardial Health- A Cardiac Sub-Study of the MCRC National Survey of Health and Development (NSHD).

A completely automated pipeline for 3D reconstruction of human heart from 2D cine magnetic resonance slices.

Cardiac magnetic resonance (CMR) imaging is a valuable modality in the diagnosis and characterization of cardiovascular diseases, since it can identify abnormalities in structure and function of the myocardium non-invasively and without the need for ionizing radiation. However, in clinical practice, it is commonly acquired as a collection of separated and independent 2D image planes, which limits its accuracy in 3D analysis. This paper presents a completely automated pipeline for generating patient-specific 3D biventricular heart models from cine magnetic resonance (MR) slices. Our pipeline automatically selects the relevant cine MR images, segments them using a deep learning-based method to extract the heart contours, and aligns the contours in 3D space correcting possible misalignments due to breathing or subject motion first using the intensity and contours information from the cine data and next with the help of a statistical shape model. Finally, the sparse 3D representation of the contours is used to generate a smooth 3D biventricular mesh. The computational pipeline is applied and evaluated in a CMR dataset of 20 healthy subjects. Our results show an average reduction of misalignment artefacts from 1.82 ± 1.60 mm to 0.72 ± 0.73 mm over 20 subjects, in terms of distance from the final reconstructed mesh. The high-resolution 3D biventricular meshes obtained with our computational pipeline are used for simulations of electrical activation patterns, showing agreement with non-invasive electrocardiographic imaging. The automatic methodologies presented here for patient-specific MR imaging-based 3D biventricular representations contribute to the efficient realization of precision medicine, enabling the enhanced interpretability of clinical data, the digital twin vision through patient-specific image-based modelling and simulation, and augmented reality applications. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Phase I/II Trial of Electrophysiology-Guided Noninvasive Cardiac Radioablation for Ventricular Tachycardia.

BACKGROUND: Case studies have suggested the efficacy of catheter-free, electrophysiology-guided noninvasive cardiac radioablation for ventricular tachycardia (VT) using stereotactic body radiation therapy, although prospective data are lacking. METHODS: We conducted a prospective phase I/II trial of noninvasive cardiac radioablation in adults with treatment-refractory episodes of VT or cardiomyopathy related to premature ventricular contractions (PVCs). Arrhythmogenic scar regions were targeted by combining noninvasive anatomic and electric cardiac imaging with a standard stereotactic body radiation therapy workflow followed by delivery of a single fraction of 25 Gy to the target. The primary safety end point was treatment-related serious adverse events in the first 90 days. The primary efficacy end point was any reduction in VT episodes (tracked by indwelling implantable cardioverter defibrillators) or any reduction in PVC burden (as measured by a 24-hour Holter monitor) comparing the 6 months before and after treatment (with a 6-week blanking window after treatment). Health-related quality of life was assessed using the Short Form-36 questionnaire. RESULTS: Nineteen patients were enrolled (17 for VT, 2 for PVC cardiomyopathy). Median noninvasive ablation time was 15.3 minutes (range, 5.4-32.3). In the first 90 days, 2/19 patients (10.5%) developed a treatment-related serious adverse event. The median number of VT episodes was reduced from 119 (range, 4-292) to 3 (range, 0-31; P<0.001). Reduction was observed for both implantable cardioverter defibrillator shocks and antitachycardia pacing. VT episodes or PVC burden were reduced in 17/18 evaluable patients (94%). The frequency of VT episodes or PVC burden was reduced by 75% in 89% of patients. Overall survival was 89% at 6 months and 72% at 12 months. Use of dual antiarrhythmic medications decreased from 59% to 12% ( P=0.008). Quality of life improved in 5 of 9 Short Form-36 domains at 6 months. CONCLUSIONS: Noninvasive electrophysiology-guided cardiac radioablation is associated with markedly reduced ventricular arrhythmia burden with modest short-term risks, reduction in antiarrhythmic drug use, and improvement in quality of life. CLINICAL TRIAL REGISTRATION: URL: https://www.clinicaltrials.gov/ . Unique identifier: NCT02919618.

Noninvasive Cardiac Radiation for Ablation of Ventricular Tachycardia.

BACKGROUND: Recent advances have enabled noninvasive mapping of cardiac arrhythmias with electrocardiographic imaging and noninvasive delivery of precise ablative radiation with stereotactic body radiation therapy (SBRT). We combined these techniques to perform catheter-free, electrophysiology-guided, noninvasive cardiac radioablation for ventricular tachycardia. METHODS: We targeted arrhythmogenic scar regions by combining anatomical imaging with noninvasive electrocardiographic imaging during ventricular tachycardia that was induced by means of an implantable cardioverter-defibrillator (ICD). SBRT simulation, planning, and treatments were performed with the use of standard techniques. Patients were treated with a single fraction of 25 Gy while awake. Efficacy was assessed by counting episodes of ventricular tachycardia, as recorded by ICDs. Safety was assessed by means of serial cardiac and thoracic imaging. RESULTS: From April through November 2015, five patients with high-risk, refractory ventricular tachycardia underwent treatment. The mean noninvasive ablation time was 14 minutes (range, 11 to 18). During the 3 months before treatment, the patients had a combined history of 6577 episodes of ventricular tachycardia. During a 6-week postablation "blanking period" (when arrhythmias may occur owing to postablation inflammation), there were 680 episodes of ventricular tachycardia. After the 6-week blanking period, there were 4 episodes of ventricular tachycardia over the next 46 patient-months, for a reduction from baseline of 99.9%. A reduction in episodes of ventricular tachycardia occurred in all five patients. The mean left ventricular ejection fraction did not decrease with treatment. At 3 months, adjacent lung showed opacities consistent with mild inflammatory changes, which had resolved by 1 year. CONCLUSIONS: In five patients with refractory ventricular tachycardia, noninvasive treatment with electrophysiology-guided cardiac radioablation markedly reduced the burden of ventricular tachycardia. (Funded by Barnes-Jewish Hospital Foundation and others.).

The Structural Basis of IKs Ion-Channel Activation: Mechanistic Insights from Molecular Simulations.

Relating ion channel (iCh) structural dynamics to physiological function remains a challenge. Current experimental and computational techniques have limited ability to explore this relationship in atomistic detail over physiological timescales. A framework associating iCh structure to function is necessary for elucidating normal and disease mechanisms. We formulated a modeling schema that overcomes the limitations of current methods through applications of artificial intelligence machine learning. Using this approach, we studied molecular processes that underlie human IKs voltage-mediated gating. IKs malfunction underlies many debilitating and life-threatening diseases. Molecular components of IKs that underlie its electrophysiological function include KCNQ1 (a pore-forming tetramer) and KCNE1 (an auxiliary subunit). Simulations, using the IKs structure-function model, reproduced experimentally recorded saturation of gating-charge displacement at positive membrane voltages, two-step voltage sensor (VS) movement shown by fluorescence, iCh gating statistics, and current-voltage relationship. Mechanistic insights include the following: 1) pore energy profile determines iCh subconductance; 2) the entire protein structure, not limited to the pore, contributes to pore energy and channel subconductance; 3) interactions with KCNE1 result in two distinct VS movements, causing gating-charge saturation at positive membrane voltages and current activation delay; and 4) flexible coupling between VS and pore permits pore opening at lower VS positions, resulting in sequential gating. The new modeling approach is applicable to atomistic scale studies of other proteins on timescales of physiological function.

Effects of ß-subunit on gating of a potassium ion channel: Molecular simulations of cardiac IKs activation.

Dynamic conformational changes of ion channel proteins during activation gating determine their function as carriers of current. The relationship between these molecular movements and channel function over the physiological timescale of the action potential (AP) has not been fully established due to limitations of existing techniques. We constructed a library of possible cardiac IKs protein conformations and applied a combination of protein segmentation and energy linearization to study this relationship computationally. Simulations reproduced the effects of the beta-subunit (KCNE1) on the alpha-subunit (KCNQ1) dynamics and function, observed in experiments. Mechanistically, KCNE1 increased the probability of "visiting" conducting pore conformations on activation trajectories, thereby increasing IKs current. KCNE1 slowed IKs activation by impeding the voltage sensor (VS) movement and reducing its coupling to pore opening. Conformational changes along activation trajectories determined that the S4-S5 linker (S4S5L) plays an important role in these modulatory effects by KCNE1. Integration of these molecular structure-based IKs dynamics into a model of human cardiac ventricular myocyte, revealed that KCNQ1-KCNE1 interaction is essential for normal AP repolarization.

Pro-arrhythmogenic Effects of the V141M KCNQ1 Mutation in Short QT Syndrome and Its Potential Therapeutic Targets: Insights from Modeling.

Gain-of-function mutations in the pore-forming subunit of IKs channels, KCNQ1, lead to short QT syndrome (SQTS) and lethal arrhythmias. However, how mutant IKs channels cause SQTS and the possibility of IKs-specific pharmacological treatment remain unclear. V141M KCNQ1 is a SQTS associated mutation. We studied its effect on IKs gating properties and changes in the action potentials (AP) of human ventricular myocytes. Xenopus oocytes were used to study the gating mechanisms of expressed V141M KCNQ1/KCNE1 channels. Computational models were used to simulate human APs in endocardial, mid-myocardial, and epicardial ventricular myocytes with and without ß-adrenergic stimulation. V141M KCNQ1 caused a gain-of-function in IKs characterized by increased current density, faster activation, and slower deactivation leading to IKs accumulation. V141M KCNQ1 also caused a leftward shift of the conductance-voltage curve compared to wild type (WT) IKs (V1/2 = 33.6 ± 4.0 mV for WT, and 24.0 ± 1.3 mV for heterozygous V141M). A Markov model of heterozygous V141M mutant IKs was developed and incorporated into the O'Hara-Rudy model. Compared to the WT, AP simulations demonstrated marked rate-dependent shortening of AP duration (APD) for V141M, predicting a SQTS phenotype. Transmural electrical heterogeneity was enhanced in heterozygous V141M AP simulations, especially under ß-adrenergic stimulation. Computational simulations identified specific IK1 blockade as a beneficial pharmacologic target for reducing the transmural APD heterogeneity associated with V141M KCNQ1 mutation. V141M KCNQ1 mutation shortens ventricular APs and enhances transmural APD heterogeneity under ß-adrenergic stimulation. Computational simulations identified IK1 blockers as a potential antiarrhythmic drug of choice for SQTS.

Cardiac electrophysiological substrate underlying the ECG phenotype and electrogram abnormalities in Brugada syndrome patients.

BACKGROUND: Brugada syndrome (BrS) is a highly arrhythmogenic cardiac disorder, associated with an increased incidence of sudden death. Its arrhythmogenic substrate in the intact human heart remains ill-defined. METHODS AND RESULTS: Using noninvasive ECG imaging, we studied 25 BrS patients to characterize the electrophysiological substrate and 6 patients with right bundle-branch block for comparison. Seven healthy subjects provided control data. Abnormal substrate was observed exclusively in the right ventricular outflow tract with the following properties (in comparison with healthy controls; P<0.005): (1) ST-segment elevation and inverted T wave of unipolar electrograms (2.21±0.67 versus 0 mV); (2) delayed right ventricular outflow tract activation (82±18 versus 37±11 ms); (3) low-amplitude (0.47±0.16 versus 3.74±1.60 mV) and fractionated electrograms, suggesting slow discontinuous conduction; (4) prolonged recovery time (381±30 versus 311±34 ms) and activation-recovery intervals (318±32 versus 241±27 ms), indicating delayed repolarization; (5) steep repolarization gradients (Δrecovery time/Δx=96±28 versus 7±6 ms/cm, Δactivation-recovery interval/Δx=105±24 versus 7±5 ms/cm) at right ventricular outflow tract borders. With increased heart rate in 6 BrS patients, reduced ST-segment elevation and increased fractionation were observed. Unlike BrS, right bundle-branch block had delayed activation in the entire right ventricle, without ST-segment elevation, fractionation, or repolarization abnormalities on electrograms. CONCLUSIONS: The results indicate that both slow discontinuous conduction and steep dispersion of repolarization are present in the right ventricular outflow tract of BrS patients. ECG imaging could differentiate between BrS and right bundle-branch block.

Electrophysiologic Scar Substrate in Relation to VT: Noninvasive High-Resolution Mapping and Risk Assessment with ECGI.

BACKGROUND: Ischemic cardiomyopathy (ICM) can provide the substrate for ventricular tachycardia (VT). OBJECTIVE: To map noninvasively with high resolution the electrophysiologic (EP) scar substrate, identify its relationship to reentry circuits during VT, and stratify VT risk in ICM patients. METHODS: Noninvasive high-resolution epicardial mapping with electrocardiographic imaging (ECGI) was performed in 32 ICM patients (17 with clinical VT, 15 without VT). Abnormal scar EP substrate was determined based on electrogram (EGM) amplitude (as percentage of maximal peak-to-peak voltage over the entire ventricular epicardium; total scar [TS] < 30%; dense scar [DS] < 15%), fractionation, and presence of late potentials (LPs). Scar burden was defined as the ratio of the scar size to the total epicardial surface area. The VT activation pattern was mapped and correlated with the EP substrate to identify components of the reentry circuit. RESULTS: Patients with VT had higher scar burden (TS: 51.0 ± 9.3% vs 36.5 ± 5.4%, P < 0.05; DS: 29.5 ± 7.3% vs 16.8 ± 6.8%, P < 0.05) with lower normalized unipolar EGM voltage (TS: 0.107 ± 0.027 vs 0.153 ± 0.031, P < 0.05; DS: 0.073 ± 0.023 vs 0.098 ± 0.026, P < 0.05), greater prevalence of fractionated EGMs (TS: 44.1 ± 10.6% vs 26.8 ± 6.3%, P < 0.05; DS: 50.8 ± 10.8% vs 30.9 ± 7.0%, P < 0.05), and LPs (TS: 26.8 ± 10.7% vs 15.8 ± 5.3, P < 0.05). VTs were mapped in eight patients; the reentry circuits were closely related to the EP substrate. CONCLUSIONS: ECGI noninvasively identified scar EP substrate that underlies abnormal conduction in ICM patients. It identified regions within the scar that aligned with critical elements of the reentry circuit during VT. ECGI can potentially be used for VT risk stratification in ICM patients.

Electrophysiologic substrate in congenital Long QT syndrome: noninvasive mapping with electrocardiographic imaging (ECGI).

BACKGROUND: Congenital Long QT syndrome (LQTS) is an arrhythmogenic disorder that causes syncope and sudden death. Although its genetic basis has become well-understood, the mechanisms whereby mutations translate to arrhythmia susceptibility in the in situ human heart have not been fully defined. We used noninvasive ECG imaging to map the cardiac electrophysiological substrate and examine whether LQTS patients display regional heterogeneities in repolarization, a substrate that promotes arrhythmogenesis. METHODS AND RESULTS: Twenty-five subjects (9 LQT1, 9 LQT2, 5 LQT3, and 2 LQT5) with genotype and phenotype positive LQTS underwent ECG imaging. Seven normal subjects provided control. Epicardial maps of activation, recovery times, activation-recovery intervals, and repolarization dispersion were constructed. Activation was normal in all patients. However, recovery times and activation-recovery intervals were prolonged relative to control, indicating delayed repolarization and abnormally long action potential duration (312±30 ms versus 235±21 ms in control). Activation-recovery interval prolongation was spatially heterogeneous, with repolarization gradients much steeper than control (119±19 ms/cm versus 2.0±2.0 ms/cm). There was variability in steepness and distribution of repolarization gradients between and within LQTS types. Repolarization gradients were steeper in symptomatic patients (130±27 ms/cm in 12 symptomatic patients versus 98±19 ms/cm in 13 asymptomatic patients; P<0.05). CONCLUSIONS: LQTS patients display regions with steep repolarization dispersion caused by localized action potential duration prolongation. This defines a substrate for reentrant arrhythmias, not detectable by surface ECG. Steeper dispersion in symptomatic patients suggests a possible role for ECG imaging in risk stratification.

Noninvasive electrocardiographic imaging of arrhythmogenic substrates in humans.

Cardiac excitation is determined by interactions between the source of electric activation (membrane depolarization) and the load that cardiac tissue presents. This relationship is altered in pathology by remodeling processes that often create a substrate favoring the development of cardiac arrhythmias. Most studies of arrhythmia mechanisms and arrhythmogenic substrates have been conducted in animal models, which may differ in important ways from the human pathologies they are designed to represent. Electrocardiographic imaging is a noninvasive method for mapping the electric activity of the heart in humans in real-world conditions. This review summarizes results from electrocardiographic imaging studies of arrhythmogenic substrates associated with human clinical arrhythmias. Examples include heart failure, myocardial infarction scar, atrial fibrillation, and abnormal ventricular repolarization.

Determinants of beat-to-beat variability of repolarization duration in the canine ventricular myocyte: a computational analysis.

Beat-to-beat variability of repolarization duration (BVR) is an intrinsic characteristic of cardiac function and a better marker of proarrhythmia than repolarization prolongation alone. The ionic mechanisms underlying baseline BVR in physiological conditions, its rate dependence, and the factors contributing to increased BVR in pathologies remain incompletely understood. Here, we employed computer modeling to provide novel insights into the subcellular mechanisms of BVR under physiological conditions and during simulated drug-induced repolarization prolongation, mimicking long-QT syndromes type 1, 2, and 3. We developed stochastic implementations of 13 major ionic currents and fluxes in a model of canine ventricular-myocyte electrophysiology. Combined stochastic gating of these components resulted in short- and long-term variability, consistent with experimental data from isolated canine ventricular myocytes. The model indicated that the magnitude of stochastic fluctuations is rate dependent due to the rate dependence of action-potential (AP) duration (APD). This process (the "active" component) and the intrinsic nonlinear relationship between membrane current and APD ("intrinsic component") contribute to the rate dependence of BVR. We identified a major role in physiological BVR for stochastic gating of the persistent Na(+) current (INa) and rapidly activating delayed-rectifier K(+) current (IKr). Inhibition of IKr or augmentation of INa significantly increased BVR, whereas subsequent ß-adrenergic receptor stimulation reduced it, similar to experimental findings in isolated myocytes. In contrast, ß-adrenergic stimulation increased BVR in simulated long-QT syndrome type 1. In addition to stochastic channel gating, AP morphology, APD, and beat-to-beat variations in Ca(2+) were found to modulate single-cell BVR. Cell-to-cell coupling decreased BVR and this was more pronounced when a model cell with increased BVR was coupled to a model cell with normal BVR. In conclusion, our results provide new insights into the ionic mechanisms underlying BVR and suggest that BVR reflects multiple potentially proarrhythmic parameters, including increased ion-channel stochasticity, prolonged APD, and abnormal Ca(2+) handling.

Electrophysiological Characterization of Subclinical and Overt Hypertrophic Cardiomyopathy by Magnetic Resonance Imaging-Guided Electrocardiography.

BACKGROUND: Ventricular arrhythmia in hypertrophic cardiomyopathy (HCM) relates to adverse structural change and genetic status. Cardiovascular magnetic resonance (CMR)-guided electrocardiographic imaging (ECGI) noninvasively maps cardiac structural and electrophysiological (EP) properties. OBJECTIVES: The purpose of this study was to establish whether in subclinical HCM (genotype [G]+ left ventricular hypertrophy [LVH]-), ECGI detects early EP abnormality, and in overt HCM, whether the EP substrate relates to genetic status (G+/G-LVH+) and structural phenotype. METHODS: This was a prospective 211-participant CMR-ECGI multicenter study of 70 G+LVH-, 104 LVH+ (51 G+/53 G-), and 37 healthy volunteers (HVs). Local activation time (AT), corrected repolarization time, corrected activation-recovery interval, spatial gradients (GAT/GRTc), and signal fractionation were derived from 1,000 epicardial sites per participant. Maximal wall thickness and scar burden were derived from CMR. A support vector machine was built to discriminate G+LVH- from HV and low-risk HCM from those with intermediate/high-risk score or nonsustained ventricular tachycardia. RESULTS: Compared with HV, subclinical HCM showed mean AT prolongation (P = 0.008) even with normal 12-lead electrocardiograms (ECGs) (P = 0.009), and repolarization was more spatially heterogenous (GRTc: P = 0.005) (23% had normal ECGs). Corrected activation-recovery interval was prolonged in overt vs subclinical HCM (P < 0.001). Mean AT was associated with maximal wall thickness; spatial conduction heterogeneity (GAT) and fractionation were associated with scar (all P < 0.05), and G+LVH+ had more fractionation than G-LVH+ (P = 0.002). The support vector machine discriminated subclinical HCM from HV (10-fold cross-validation accuracy 80% [95% CI: 73%-85%]) and identified patients at higher risk of sudden cardiac death (accuracy 82% [95% CI: 78%-86%]). CONCLUSIONS: In the absence of LVH or 12-lead ECG abnormalities, HCM sarcomere gene mutation carriers express an aberrant EP phenotype detected by ECGI. In overt HCM, abnormalities occur more severely with adverse structural change and positive genetic status.

Ionic bases for electrical remodeling of the canine cardiac ventricle.

Emerging evidence suggests that ventricular electrical remodeling (VER) is triggered by regional myocardial strain via mechanoelectrical feedback mechanisms; however, the ionic mechanisms underlying strain-induced VER are poorly understood. To determine its ionic basis, VER induced by altered electrical activation in dogs undergoing left ventricular pacing (n = 6) were compared with unpaced controls (n = 4). Action potential (AP) durations (APDs), ionic currents, and Ca(2+) transients were measured from canine epicardial myocytes isolated from early-activated (low strain) and late-activated (high strain) left ventricular regions. VER in the early-activated region was characterized by minimal APD prolongation, but marked attenuation of the AP phase 1 notch attributed to reduced transient outward K(+) current. In contrast, VER in the late-activated region was characterized by significant APD prolongation. Despite marked APD prolongation, there was surprisingly minimal change in ion channel densities but a twofold increase in diastolic Ca(2+). Computer simulations demonstrated that changes in sarcolemmal ion channel density could only account for attenuation of the AP notch observed in the early-activated region but failed to account for APD remodeling in the late-activated region. Furthermore, these simulations identified that cytosolic Ca(2+) accounted for APD prolongation in the late-activated region by enhancing forward-mode Na(+)/Ca(2+) exchanger activity, corroborated by increased Na(+)/Ca(2+) exchanger protein expression. Finally, assessment of skinned fibers after VER identified altered myofilament Ca(2+) sensitivity in late-activated regions to be associated with increased diastolic levels of Ca(2+). In conclusion, we identified two distinct ionic mechanisms that underlie VER: 1) strain-independent changes in early-activated regions due to remodeling of sarcolemmal ion channels with no changes in Ca(2+) handling and 2) a novel and unexpected mechanism for strain-induced VER in late-activated regions in the canine arising from remodeling of sarcomeric Ca(2+) handling rather than sarcolemmal ion channels.

A model of canine purkinje cell electrophysiology and Ca(2+) cycling: rate dependence, triggered activity, and comparison to ventricular myocytes.

Purkinje cells (Pcell) are characterized by different electrophysiological properties and Ca(2+) cycling processes than ventricular myocytes (Vcell) and are frequently involved in ventricular arrhythmias. Yet, the mechanistic basis for their arrhythmic vulnerability is not completely understood. The objectives were to: (1) characterize Pcell electrophysiology, Ca(2+) cycling, and their rate dependence; (2) investigate mechanisms underlying Pcell arrhythmogenicity; and compare Pcell and Vcell electrophysiology, Ca(2+) cycling, and arrhythmic properties. We developed a new mathematical model of Pcell. The Ca(2+) subsystem includes spatial organization and receptors distribution unique to Pcell. Results were: (1) in Pcell and Vcell, Na(+) accumulation via its augmentation of repolarizing I(NaK) dominates action potential duration adaptation and, in Pcell, I(NaL) contributes additional action potential duration shortening at short cycle length; (2) steep Pcell restitution is attributable to slow recovery of I(NaL); (3) biphasic Ca(2+) transients of Pcell reflect the delay between Ca(2+) release from junctional sarcoplasmic reticulum and corbular sarcoplasmic reticulum; (4) Pcell Ca(2+) alternans, unlike Vcell, can develop without inducing action potential alternans; (5) Pcell action potential alternans develops at a shorter cycle length than Vcell, with increased subcellular heterogeneity of Ca(2+) cycling attributable to refractoriness of Ca(2+) release from corbular sarcoplasmic reticulum and junctional sarcoplasmic reticulum; (6) greater Pcell vulnerability to delayed afterdepolarizations is attributable to higher sarcoplasmic reticulum Ca(2+) content and ionic currents that reduce excitation threshold and promote triggered activity; and (7) early after depolarizations generation in Pcell is mostly attributable to reactivation of I(NaL2), whereas I(CaL) plays this role in Vcell. Steeper rate dependence of action potential and Ca(2+) transients, central peripheral heterogeneity of Ca(2+) cycling, and distinct ion channel profile underlie greater arrhythmic vulnerability of Pcell compared to Vcell.

Noninvasive mapping of ventricular activation in patients with transplanted hearts.

This is the first reported study of ventricular activation patterns after cardiac transplantation, using electrocardiographic imaging (ECGI), a noninvasive method for electrophysiologic mapping. This study of ten patients reveals that transplanted hearts have unique ventricular activation patterns in sinus rhythm, activating early in the epicardial aspect of the anterior or inferior septum, with intact right and left bundle branch conduction. They have late activation with slowing of conduction near the right ventricular (RV) basal free wall, causing a mild QRS prolongation and an rSr' pattern in lead V1 of the ECG. PVCs arise from both endocardial and epicardial locations in both ventricles.

Data-based theoretical identification of subcellular calcium compartments and estimation of calcium dynamics in cardiac myocytes.

In cardiac cells, Ca(2+) release flux (J(rel)) via ryanodine receptors (RyRs) from the sarcoplasmic reticulum (SR) has a complex effect on the action potential (AP). Coupling between J(rel) and the AP occurs via L-type Ca(2+) channels (I(Ca)) and the Na(+)/Ca(2+) exchanger (I(NCX)). We used a combined experimental and modelling approach to study interactions between J(rel), I(Ca) and I(NCX) in porcine ventricular myocytes.We tested the hypothesis that during normal uniform J(rel), the interaction between these fluxes can be represented as occurring in two myoplasmic subcompartments for Ca(2+) distribution, one (T-space) associated with RyR and enclosed by the junctional portion of the SR membrane and corresponding T-tubular portion of the sarcolemma, the other (M-space) encompassing the rest of the myoplasm. I(Ca) and I(NCX) were partitioned into subpopulations in the T-space and M-space sarcolemma. We denoted free Ca(2+) concentrations in T-space and M-space Ca(t) and Ca(m), respectively. Experiments were designed to allow separate measurements of I(Ca) and I(NCX) as a function of J(rel). Inclusion of T-space in themodel allowed us to reproduce in silico the following important experimental results: (1) hysteresis of I(NCX) dependence on Ca(m); (2) delay between peak I(NCX) and peak Ca(m) during caffeine application protocol; (3) delay between I(NCX) and Ca(m) during Ca(2+)-induced-Ca(2+)-release; (4) rapid I(Ca) inactivation (within 2 ms) due to J(rel), with magnitude graded as a function of the SR Ca(2+) content; (5) time delay between I(Ca) inactivation due to J(rel) and Ca(m). Partition of 25% NCX in T-space and 75% in M-space provided the best fit to the experimental data. Measured Ca(m) and I(Ca) or I(NCX) were used as input to the model for estimating Ca(t). The actual model-computed Ca(t), obtained by simulating specific experimental protocols, was used as a gold standard for comparison. The model predicted peak Ca(t) in the range of 6­25 µM, with time to equilibrium of Ca(t) with Ca(m) of ~350 ms. These Ca(t) values are in the range of LCC and RyR sensitivity to Ca(2+). An increase of the SR Ca(2+) load increased the time to equilibrium. The I(Ca)-based estimation method was most accurate during the ascending phase of Ca(t). The I(NCX)-based method provided a good estimate for the descending phase of Ca(t). Thus, application of both methods in combination provides the best estimate of the entire Ca(t) time course.