How to measure propagation velocity in cardiac tissue: a simulation study

Brain-to-heart cholinergic synapse-calcium signaling mediates ischemic stroke-induced atrial fibrillation

Background: Stroke-related cardiovascular diseases have attracted considerable attention, with atrial fibrillation (AF) being among the most frequent complications. Despite increasing clinical evidence, experimental models of stroke-induced AF are still lacking, hindering mechanistic discoveries and the development of adequate therapeutics targeting this stroke-heart syndrome (SHS). This study aims to create a rat model of ischemic stroke-induced AF (ISIAF) and to explore the efficacy and mechanism of Wenxin Keli (WK), an antiarrhythmic Chinese medicine. Method: The middle cerebral artery occlusion/reperfusion model was adapted to create subacute brain ischemia in rats with normal cardiac function. Invasive electrophysiologic studies and ex vivo optical mapping were performed to evaluate the altered electrophysiological parameters and Ca2+ handling properties. RNA-seq analysis, RT-PCR, and immunohistochemistry (IHC) with immunofluorescence (IF) were employed to assess the SHS model and elucidate the mechanisms of ISIAF and the effects of WK. UPLC/Q-TOF-MS, molecular docking, and whole-cell patch recordings were used to identify the active components of WK for SHS. Results: Ischemic stroke aggravated atrial electrical instability, altered action potential duration (APD), Ca2+ transient duration (CaT), conduction heterogeneity, and spatially discordant alternans in SHS rat hearts. These abnormalities were alleviated by WK. RNA-seq analysis revealed that M3-mediated cholinergic synapse signaling and L-type calcium channel (LTCCs)-mediated Ca2+ signaling play prominent roles in ISIAF development and its reversal by WK. UPLC/Q-TOF-MS analysis identified 19 WK components as the main components in plasma after WK treatment. Molecular docking screening identified Dioscin as the major active component of WK. WK and Dioscin reduced ICa-L in a concentration-dependent manner with a half-maximal inhibitory concentration of 24.254 ± 2.051 mg/mL and 8.666 ± 0.777 µmol/L, respectively. Conclusion: This study established an experimental model of ISIAF capable of characterizing clinically relevant atrial electrophysiological changes post-cerebral ischemia. Molecular mechanistic studies revealed that the cholinergic-calcium signaling pathway is central to this brain-heart syndrome. Ischemic stroke-induced atrial fibrillation is partially reversible by the Chinese medicine Wenxin Keli, which acts via regulation of the cholinergic-calcium signaling pathway, with its active component Dioscin directly binding to IKM3 and inhibiting ICa-L.

The importance of electrode-tissue proximity in creating pulsed field ablation lesions: insights from a sub-acute preclinical model

BACKGROUND

We sought to evaluate the anatomic and functional lesion development over time at different atrial sites immediately following delivery of pulsed field ablation (PFA).

METHODS

Using a porcine model, PFA ablations were performed in the superior vena cava (SVC), right atrial lateral wall (RA), left atrial appendage (LAA), and right superior pulmonary vein (RSPV) using four different PFA profiles. Mapping was done sequentially in 5-20-min increments up to 280-min post lesion delivery for low voltage area (LVA) assessment and conduction velocity. Lesion characteristics were noted with voltage mapping immediately post ablation and at the serial time points.

RESULTS

In 9 animals, 33 sites were ablated. None of the four different profiles across all sites showed any statistical difference on acute lesion formation or persistence. Higher tissue contact was observed in the SVC and RSPV and lower tissue contact was observed in the LAA and RA locations. Higher contact areas were noted to have higher density electroanatomic low voltage area (LVA) (12/14 vs 5/18, p = 0.01) and larger lesions on gross pathology (2 /14 vs 6/16, p = 0.01) compared to lower contact areas. Lesion regression occurred in 16/33 sites. Sustained lesions were significantly more prevalent in higher versus lower contact sites (65% vs 38%, p = 0.037).

CONCLUSION

The development of significant and durable lesions for PFA in a porcine model appears to be dependent on tissue proximity and contact.

Atrial conduction velocity mapping: clinical tools, algorithms and approaches for understanding the arrhythmogenic substrate

Characterizing patient-specific atrial conduction properties is important for understanding arrhythmia drivers, for predicting potential arrhythmia pathways, and for personalising treatment approaches. One metric that characterizes the health of the myocardial substrate is atrial conduction velocity, which describes the speed and direction of propagation of the electrical wavefront through the myocardium. Atrial conduction velocity mapping algorithms are under continuous development in research laboratories and in industry. In this review article, we give a broad overview of different categories of currently published methods for calculating CV, and give insight into their different advantages and disadvantages overall. We classify techniques into local, global, and inverse methods, and discuss these techniques with respect to their faithfulness to the biophysics, incorporation of uncertainty quantification, and their ability to take account of the atrial manifold.

Circle Method for Robust Estimation of Local Conduction Velocity High-Density Maps From Optical Mapping Data: Characterization of Radiofrequency Ablation Sites

Conduction velocity (CV) slowing is associated with atrial fibrillation (AF) and reentrant ventricular tachycardia (VT). Clinical electroanatomical mapping systems used to localize AF or VT sources as ablation targets remain limited by the number of measuring electrodes and signal processing methods to generate high-density local activation time (LAT) and CV maps of heterogeneous atrial or trabeculated ventricular endocardium. The morphology and amplitude of bipolar electrograms depend on the direction of propagating electrical wavefront, making identification of low-amplitude signal sources commonly associated with fibrotic area difficulty. In comparison, unipolar electrograms are not sensitive to wavefront direction, but measurements are susceptible to distal activity. This study proposes a method for local CV calculation from optical mapping measurements, termed the circle method (CM). The local CV is obtained as a weighted sum of CV values calculated along different chords spanning a circle of predefined radius centered at a CV measurement location. As a distinct maximum in LAT differences is along the chord normal to the propagating wavefront, the method is adaptive to the propagating wavefront direction changes, suitable for electrical conductivity characterization of heterogeneous myocardium. In numerical simulations, CM was validated characterizing modeled ablated areas as zones of distinct CV slowing. Experimentally, CM was used to characterize lesions created by radiofrequency ablation (RFA) on isolated hearts of rats, guinea pig, and explanted human hearts. To infer the depth of RFA-created lesions, excitation light bands of different penetration depths were used, and a beat-to-beat CV difference analysis was performed to identify CV alternans. Despite being limited to laboratory research, studies based on CM with optical mapping may lead to new translational insights into better-guided ablation therapies.

Cardiac Conduction Velocity, Remodeling and Arrhythmogenesis

Cardiac electrophysiological disorders, in particular arrhythmias, are a key cause of morbidity and mortality throughout the world. There are two basic requirements for arrhythmogenesis: an underlying substrate and a trigger. Altered conduction velocity (CV) provides a key substrate for arrhythmogenesis, with slowed CV increasing the probability of re-entrant arrhythmias by reducing the length scale over which re-entry can occur. In this review, we examine methods to measure cardiac CV in vivo and ex vivo, discuss underlying determinants of CV, and address how pathological variations alter CV, potentially increasing arrhythmogenic risk. Finally, we will highlight future directions both for methodologies to measure CV and for possible treatments to restore normal CV.

A practical guide for investigating cardiac physiology using living myocardial slices

Ex vivo multicellular preparations are essential tools to study tissue physiology. Among them, the recent methodological and technological developments in living myocardial slices (LMS) are attracting increasing interest by the cardiac research field. Despite this, this research model remains poorly perceived and utilized by most research laboratories. Here, we provide a practical guide on how to use LMS to interrogate multiple aspects of cardiac function, structure and biochemistry. We discuss issues that should be considered to conduct successful experiments, including experimental design, sample preparation, data collection and analysis. We describe how laboratory setups can be adapted to accommodate and interrogate this multicellular research model. These adaptations can often be achieved at a reasonable cost with off-the-shelf components and operated reliably using well-established protocols and freely available software, which is essential to broaden the utilization of this method. We will also highlight how current measurements can be improved to further enhance data quality and reliability to ensure inter-laboratory reproducibility. Finally, we summarize the most promising biomedical applications and envision how living myocardial slices can lead to further breakthroughs.

Activation During Sinus Rhythm in Ventricles With Healed Infarction: Differentiation Between Arrhythmogenic and Nonarrhythmogenic Scar

BACKGROUND

In infarct-related ventricular tachycardia (VT), the circuit often corresponds to a location characterized by activation slowing during sinus rhythm (SR). However, the relationship between activation slowing during SR and vulnerability for reentry and correlation to components of the VT circuit are unknown. This study examined the relationship between activation slowing during SR and vulnerability for reentry and correlated these areas with components of the circuit.

METHODS

In a porcine model of healed infarction, the spatial distribution of endocardial activation velocity was compared between SR and VT. Isthmus sites were defined using activation and entrainment mapping as areas exhibiting diastolic activity within the circuit while bystanders were defined as areas displaying diastolic activity outside the circuit.

RESULTS

Of 15 swine, 9 had inducible VT (5.2±3.0 per animal) while in 6 swine VT could not be induced despite stimulation from 4 RV and LV sites at 2 drive trains with 6 extra-stimuli down to refractoriness. Infarcts with VT had a greater magnitude of activation slowing during SR. A minimal endocardial activation velocity cutoff ≤0.1 m/s differentiated inducible from noninducible infarctions (P=0.015). Regions of maximal endocardial slowing during SR corresponded to the VT isthmus (area under curve=0.84 95% CI, 0.78-0.90) while bystander sites exhibited near-normal activation during SR. VT circuits were complex with 41.7% exhibiting discontinuous propagation with intramural bridges of slow conduction and delayed quasi-simultaneous endocardial activation. Regions forming the VT isthmus borders had faster activation during SR while regions forming the inner isthmus were activated faster during VT.

CONCLUSIONS

Endocardial activation slowing during SR may differentiate infarctions vulnerable for VT from those less vulnerable for VT. Sites of slow activation during SR correspond to sites forming the VT isthmus but not to bystander sites.

ElectroMap: High-throughput open-source software for analysis and mapping of cardiac electrophysiology

The ability to record and analyse electrical behaviour across the heart using optical and electrode mapping has revolutionised cardiac research. However, wider uptake of these technologies is constrained by the lack of multi-functional and robustly characterised analysis and mapping software. We present ElectroMap, an adaptable, high-throughput, open-source software for processing, analysis and mapping of complex electrophysiology datasets from diverse experimental models and acquisition modalities. Key innovation is development of standalone module for quantification of conduction velocity, employing multiple methodologies, currently not widely available to researchers. ElectroMap has also been designed to support multiple methodologies for accurate calculation of activation, repolarisation, arrhythmia detection, calcium handling and beat-to-beat heterogeneity. ElectroMap implements automated signal segmentation, ensemble averaging and integrates optogenetic approaches. Here we employ ElectroMap for analysis, mapping and detection of pro-arrhythmic phenomena in silico, in cellulo, animal model and in vivo patient datasets. We anticipate that ElectroMap will accelerate innovative cardiac research and enhance the uptake, application and interpretation of mapping technologies leading to novel approaches for arrhythmia prevention.

A technique for measuring anisotropy in atrial conduction to estimate conduction velocity and atrial fibre direction

BACKGROUND

Cardiac conduction properties exhibit large variability, and affect patient-specific arrhythmia mechanisms. However, it is challenging to clinically measure conduction velocity (CV), anisotropy and fibre direction. Our aim is to develop a technique to estimate conduction anisotropy and fibre direction from clinically available electrical recordings.

METHODS

We developed and validated automated algorithms for estimating cardiac CV anisotropy, from any distribution of recording locations on the atrial surface. The first algorithm is for elliptical wavefront fitting to a single activation map (method 1), which works well close to the pacing location, but decreases in accuracy further from the pacing location (due to spatial heterogeneity in the conductivity and fibre fields). As such, we developed a second methodology for measuring local conduction anisotropy, using data from two or three activation maps (method 2: ellipse fitting to wavefront propagation velocity vectors from multiple activation maps).

RESULTS

Ellipse fitting to CV vectors from two activation maps (method 2) leads to an improved estimation of longitudinal and transverse CV compared to method 1, but fibre direction estimation is still relatively poor. Using three activation maps with method 2 provides accurate estimation, with approximately 70% of atrial fibres estimated within 20∘. We applied the technique to clinical activation maps to demonstrate the presence of heterogeneous conduction anisotropy, and then tested the effects of this conduction anisotropy on predicted arrhythmia dynamics using computational simulation.

CONCLUSIONS

We have developed novel algorithms for calculating CV and measuring the direction dependency of atrial activation to estimate atrial fibre direction, without the need for specialised pacing protocols, using clinically available electrical recordings.

Automated isochronal late activation mapping to identify deceleration zones: Rationale and methodology of a practical electroanatomic mapping approach for ventricular tachycardia ablation

Sinus rhythm surrogates for critical isthmus sites are highly desirable because the vast majority of VT is hemodynamically unstable. While many ablation strategies to decrease the arrhythmogenicity of scar have been shown to be effective, the predominant method for electroanatomic mapping relies on a voltage-based depiction of scar and abnormal electrograms. A functional prioritization of slow conduction, distinct from late activation, is feasible in clinical practice with the creation of isochronal late activation maps. Regions of slow conduction are easily visualized with isochronal displays of baseline intrinsic rhythm activation and deceleration zones, where isochrones crowd, have been observed to have strong correlation with successful ablation sites. Automated annotation of the offset of local electrograms was developed to create the propagation maps to incorporate electrogram width and completion of local activation. Simple conduction velocity estimates where three isochrones are seen within a 1 cm radium confirm that deceleration zones harbor conduction velocity of <0.6 m/s. We present a practical methodology of analyzing electroanatomic substrate in a voltage-independent manner with correlation to reentrant VT. Non-linear 3D transmyocardial conduction limits the validity of conduction velocity estimates that assume planar and tangential conduction and we show an example of a patient with 3D isthmus boundaries with an activation gap on the epicardial surface during tachycardia.

Ccoffinn: Automated Wave Tracking in Cultured Cardiac Monolayers

Cardiac arrhythmias are one of the most frequent causes of death worldwide. A popular biological model used to study arrhythmogenesis is the cultured cardiac cell monolayer, which provides a good trade-off between physiological relevance and experimental access. Excitation wave patterns are imaged using high-bandwidth detectors, producing large data sets that are typically analyzed manually. To make such analysis less time consuming and less subjective, we have designed and implemented a toolkit for segmentation and tracking of cardiac waves in optical mapping recordings. The toolkit is optimized for high-resolution detectors to accommodate the growing availability of inexpensive high-resolution detectors for life science imaging applications (e.g., scientific CMOS cameras). The software extracts key features of propagating waves, such as wavefront speed and entropy. The methods have been validated using synthetic data, and real-world examples are provided, showing a difference in conduction velocity between two different types of cardiac cell cultures.

High-Resolution Mapping of Postinfarction Reentrant Ventricular Tachycardia: Electrophysiological Characterization of the Circuit

BACKGROUND

In vivo description of ventricular tachycardia (VT) circuits is limited by insufficient spatiotemporal resolution. We used a novel high-resolution mapping technology to characterize the electrophysiological properties of the postinfarction reentrant VT circuit.

METHODS

In 15 swine, myocardial infarction was induced by left anterior descending artery balloon occlusion. Animals were studied 6 to 8 weeks after myocardial infarction. Activation mapping of VTs was performed by using the Rhythmia mapping system. Activation time was based on a combination of bipolar and unipolar electrograms. The response to overdrive pacing from different zones of the circuit was examined.

RESULTS

A total of 56 monomorphic VTs were induced (3.8±2.1 per animal). Among these, 21 (37.5%) were hemodynamically stable and allowed mapping of the circuit. Isthmuses were 16.4±7.2 mm long and 7.4±2.8 mm wide. Conduction velocities were slowest at the inward curvature into the isthmus entrance (0.28±0.2 m/s), slightly faster at the outward curvature exit (0.40±0.3 m/s) and nearly normal at the central isthmus (0.62±0.2 m/s). In 3 animals, 2 VT morphologies with opposite axes sharing the same isthmus were mapped. Conduction velocities within the shared isthmus were dependent on the activation vector, consistently slower at the proximal curvature. Overdrive pacing from isthmus sites determined by activation mapping was consistent with entrainment criteria for isthmus. However, dimensions of the isthmus defined by entrainment exceeded dimensions of the isthmus measured by activation mapping by 32±18%.

CONCLUSIONS

In postinfarction reentrant VT, conduction velocities are slowest at the proximal and distal curvatures. Entrainment mapping overestimates the true size of the isthmus. High-resolution activation mapping of VT may better guide ablation therapy.

Feasibility of a semi-automated method for cardiac conduction velocity analysis of high-resolution activation maps

Myocardial conduction velocity is important for the genesis of arrhythmias. In the normal heart, conduction is primarily dependent on fiber direction (anisotropy) and may be discontinuous at sites with tissue heterogeneities (trabeculated or fibrotic tissue). We present a semi-automated method for the accurate measurement of conduction velocity based on high-resolution activation mapping following central stimulation. The method was applied to activation maps created from myocardium from man, sheep and mouse with anisotropic and discontinuous conduction. Advantages of the presented method over existing methods are discussed.

Dynamics of propagation of premature impulses in structurally remodeled infarcted myocardium: a computational analysis

Initiation of cardiac arrhythmias typically follows one or more premature impulses either occurring spontaneously or applied externally. In this study, we characterize the dynamics of propagation of single (S2) and double premature impulses (S3), and the mechanisms of block of premature impulses at structural heterogeneities caused by remodeling of gap junctional conductance (Gj) in infarcted myocardium. Using a sub-cellular computer model of infarcted tissue, we found that |I_Na,max|, prematurity (coupling interval with the previous impulse), and conduction velocity (CV) of premature impulses change dynamically as they propagate away from the site of initiation. There are fundamental differences between the dynamics of propagation of S2 and S3 premature impulses: for S2 impulses |I_Na,max| recovers fast, prematurity decreases and CV increases as propagation proceeds; for S3 impulses low values of |I_Na,max| persist, prematurity could increase, and CV could decrease as impulses propagate away from the site of initiation. As a consequence it is more likely that S3 impulses block at sites of structural heterogeneities causing source/sink mismatch than S2 impulses block. Whether premature impulses block at Gj heterogeneities or not is also determined by the values of Gj (and the space constant λ) in the regions proximal and distal to the heterogeneity: when λ in the direction of propagation increases >40%, premature impulses could block. The maximum slope of CV restitution curves for S2 impulses is larger than for S3 impulses. In conclusion: (1) The dynamics of propagation of premature impulses make more likely that S3 impulses block at sites of structural heterogeneities than S2 impulses block; (2) Structural heterogeneities causing an increase in λ (or CV) of >40% could result in block of premature impulses; (3) A decrease in the maximum slope of CV restitution curves of propagating premature impulses is indicative of an increased potential for block at structural heterogeneities.