Evaluation and optimization of novel extraction algorithms for the automatic detection of atrial activations recorded within the pulmonary veins during atrial fibrillation.

BACKGROUND AND OBJECTIVE: The automated detection of atrial activations (AAs) recorded from intracardiac electrograms (IEGMs) during atrial fibrillation (AF) is challenging considering their various amplitudes, morphologies and cycle length. Activation time estimation is further complicated by the constant changes in the IEGM active zones in complex and/or fractionated signals. We propose a new method which provides reliable automatic extraction of intracardiac AAs recorded within the pulmonary veins during AF and an accurate estimation of their local activation times. METHODS: First, two recently developed algorithms were evaluated and optimized on 118 recordings of pulmonary vein IEGM taken from 35 patients undergoing ablation of persistent AF. The adaptive mathematical morphology algorithm (AMM) uses an adaptive structuring element to extract AAs based on their morphological features. The relative-energy algorithm (Rel-En) uses short- and long-term energies to enhance and detect the AAs in the IEGM signals. Second, following the AA extraction, the signal amplitude was weighted using statistics of the AA sequences in order to reduce over- and undersensing of the algorithms. The detection capacity of our algorithms was compared with manually annotated activations and with two previously developed algorithms based on the Teager-Kaiser energy operator and the AF cycle length iteration, respectively. Finally, a method based on the barycenter was developed to reduce artificial variations in the activation annotations of complex IEGM signals. RESULTS: The best detection was achieved using Rel-En, yielding a false negative rate of 0.76% and a false positive rate of only 0.12% (total error rate 0.88%) against expert annotation. The post-processing further reduced the total error rate of the Rel-En algorithm by 70% (yielding to a final total error rate of 0.28%). CONCLUSION: The proposed method shows reliable detection and robust temporal annotation of AAs recorded within pulmonary veins in AF. The method has low computational cost and high robustness for automatic detection of AAs, which makes it a suitable approach for online use in a procedural context.

Stochastic pacing reveals the propensity to cardiac action potential alternans and uncovers its underlying dynamics.

KEY POINTS: Beat-to-beat alternation (alternans) of the cardiac action potential duration is known to precipitate life-threatening arrhythmias and can be driven by the kinetics of voltage-gated membrane currents or by instabilities in intracellular calcium fluxes. To prevent alternans and associated arrhythmias, suitable markers must be developed to quantify the susceptibility to alternans; previous theoretical studies showed that the eigenvalue of the alternating eigenmode represents an ideal marker of alternans. Using rabbit ventricular myocytes, we show that this eigenvalue can be estimated in practice by pacing these cells at intervals varying stochastically. We also show that stochastic pacing permits the estimation of further markers distinguishing between voltage-driven and calcium-driven alternans. Our study opens the perspective to use stochastic pacing during clinical investigations and in patients with implanted pacing devices to determine the susceptibility to, and the type of alternans, which are both important to guide preventive or therapeutic measures. ABSTRACT: Alternans of the cardiac action potential (AP) duration (APD) is a well-known arrhythmogenic mechanism. APD depends on several preceding diastolic intervals (DIs) and APDs, which complicates the prediction of alternans. Previous theoretical studies pinpointed a marker called λalt that directly quantifies how an alternating perturbation persists over successive APs. When the propensity to alternans increases, λalt decreases from 0 to -1. Our aim was to quantify λalt experimentally using stochastic pacing and to examine whether stochastic pacing allows discriminating between voltage-driven and Ca(2+) -driven alternans. APs were recorded in rabbit ventricular myocytes paced at cycle lengths (CLs) decreasing progressively and incorporating stochastic variations. Fitting APD with a function of two previous APDs and CLs permitted us to estimate λalt along with additional markers characterizing whether the dependence of APD on previous DIs or CLs is strong (typical for voltage-driven alternans) or weak (Ca(2+) -driven alternans). During the recordings, λalt gradually decreased from around 0 towards -1. Intermittent alternans appeared when λalt reached -0.8 and was followed by sustained alternans. The additional markers detected that alternans was Ca(2+) driven in control experiments and voltage driven in the presence of ryanodine. This distinction could be made even before alternans was manifest (specificity/sensitivity >80% for -0.4 > λalt  > -0.5). These observations were confirmed in a mathematical model of a rabbit ventricular myocyte. In conclusion, stochastic pacing allows the practical estimation of λalt to reveal the onset of alternans and distinguishes between voltage-driven and Ca(2+) -driven mechanisms, which is important since these two mechanisms may precipitate arrhythmias in different manners.

Slow conduction in mixed cultured strands of primary ventricular cells and stem cell-derived cardiomyocytes.

Modern concepts for the treatment of myocardial diseases focus on novel cell therapeutic strategies involving stem cell-derived cardiomyocytes (SCMs). However, functional integration of SCMs requires similar electrophysiological properties as primary cardiomyocytes (PCMs) and the ability to establish intercellular connections with host myocytes in order to contribute to the electrical and mechanical activity of the heart. The aim of this project was to investigate the properties of cardiac conduction in a co-culture approach using SCMs and PCMs in cultured cell strands. Murine embryonic SCMs were pooled with fetal ventricular cells and seeded in predefined proportions on microelectrode arrays to form patterned strands of mixed cells. Conduction velocity (CV) was measured during steady state pacing. SCM excitability was estimated from action potentials measured in single cells using the patch clamp technique. Experiments were complemented with computer simulations of conduction using a detailed model of cellular architecture in mixed cell strands. CV was significantly lower in strands composed purely of SCMs (5.5 ± 1.5 cm/s, n = 11) as compared to PCMs (34.9 ± 2.9 cm/s, n = 21) at similar refractoriness (100% SCMs: 122 ± 25 ms, n = 9; 100% PCMs: 139 ± 67 ms, n = 14). In mixed strands combining both cell types, CV was higher than in pure SCMs strands, but always lower than in 100% PCM strands. Computer simulations demonstrated that both intercellular coupling and electrical excitability limit CV. These data provide evidence that in cultures of murine ventricular cardiomyocytes, SCMs cannot restore CV to control levels resulting in slow conduction, which may lead to reentry circuits and arrhythmias.

Nonlinear behaviour of conduction and block in cardiac tissue with heterogeneous expression of connexin 43.

Altered gap junctional coupling potentiates slow conduction and arrhythmias. To better understand how heterogeneous connexin expression affects conduction at the cellular scale, we investigated conduction in tissue consisting of two cardiomyocyte populations expressing different connexin levels. Conduction was mapped using microelectrode arrays in cultured strands of foetal murine ventricular myocytes with predefined contents of connexin 43 knockout (Cx43KO) cells. Corresponding computer simulations were run in randomly generated two-dimensional tissues mimicking the cellular architecture of the strands. In the cultures, the relationship between conduction velocity (CV) and Cx43KO cell content was nonlinear. CV first decreased significantly when Cx43KO content was increased from 0 to 50%. When the Cx43KO content was ≥60%, CV became comparable to that in 100% Cx43KO strands. Co-culturing Cx43KO and wild-type cells also resulted in significantly more heterogeneous conduction patterns and in frequent conduction blocks. The simulations replicated this behaviour of conduction. For Cx43KO contents of 10-50%, conduction was slowed due to wavefront meandering between Cx43KO cells. For Cx43KO contents ≥60%, clusters of remaining wild-type cells acted as electrical loads that impaired conduction. For Cx43KO contents of 40-60%, conduction exhibited fractal characteristics, was prone to block, and was more sensitive to changes in ion currents compared to homogeneous tissue. In conclusion, conduction velocity and stability behave in a nonlinear manner when cardiomyocytes expressing different connexin amounts are combined. This behaviour results from heterogeneous current-to-load relationships at the cellular level. Such behaviour is likely to be arrhythmogenic in various clinical contexts in which gap junctional coupling is heterogeneous.

Measures of spatiotemporal organization differentiate persistent from long-standing atrial fibrillation.

AIMS: This study presents an automatic diagnostic method for the discrimination between persistent and long-standing atrial fibrillation (AF) based on the surface electrocardiogram (ECG). METHODS AND RESULTS: Standard 12-lead ECG recordings were acquired in 53 patients with either persistent (N= 20) or long-standing AF (N= 33), the latter including both long-standing persistent and permanent AF. A combined frequency analysis of multiple ECG leads followed by the computation of standard complexity measures provided a method for the quantification of spatiotemporal AF organization. All possible pairs of precordial ECG leads were analysed by this method and resulting organization measures were used for automatic classification of persistent and long-standing AF signals. Correct classification rates of 84.9% were obtained, with a predictive value for long-standing AF of 93.1%. Spatiotemporal organization as measured in lateral precordial leads V5 and V6 was shown to be significantly lower during long-standing AF than persistent AF, suggesting that time-related alterations in left atrial electrical activity can be detected in the ECG. CONCLUSION: Accurate discrimination between persistent and long-standing AF based on standard surface recordings was demonstrated. This information could contribute to optimize the management of sustained AF, permitting appropriate therapeutic decisions and thereby providing substantial clinical cost savings.

Adaptive tracking of EEG oscillations.

Neuronal oscillations are an important aspect of EEG recordings. These oscillations are supposed to be involved in several cognitive mechanisms. For instance, oscillatory activity is considered a key component for the top-down control of perception. However, measuring this activity and its influence requires precise extraction of frequency components. This processing is not straightforward. Particularly, difficulties with extracting oscillations arise due to their time-varying characteristics. Moreover, when phase information is needed, it is of the utmost importance to extract narrow-band signals. This paper presents a novel method using adaptive filters for tracking and extracting these time-varying oscillations. This scheme is designed to maximize the oscillatory behavior at the output of the adaptive filter. It is then capable of tracking an oscillation and describing its temporal evolution even during low amplitude time segments. Moreover, this method can be extended in order to track several oscillations simultaneously and to use multiple signals. These two extensions are particularly relevant in the framework of EEG data processing, where oscillations are active at the same time in different frequency bands and signals are recorded with multiple sensors. The presented tracking scheme is first tested with synthetic signals in order to highlight its capabilities. Then it is applied to data recorded during a visual shape discrimination experiment for assessing its usefulness during EEG processing and in detecting functionally relevant changes. This method is an interesting additional processing step for providing alternative information compared to classical time-frequency analyses and for improving the detection and analysis of cross-frequency couplings.

Phase-rectified signal averaging used to estimate the dominant frequencies in ECG signals during atrial fibrillation.

Phase-rectified signal averaging (PRSA) is a technique recently introduced to enhance quasi-periodic signal components. An important parameter that can be extracted from surface ECG is the dominant frequency (DF) of atrial fibrillation (AF). AF signal components are always highly contaminated by the ventricular complexes, and the cancellation of these components is never perfect. The remaining artifacts tend to induce erroneous DF estimates. In this paper, we report on the use of PRSA in the context of noninvasive AF classification procedures for improving DF estimation. The potential of PRSA is demonstrated by experiments both on synthetic and clinical ECG signals.

Propagation velocity kinetics and repolarization alternans in a free-behaving sheep model of pacing-induced atrial fibrillation.

AIMS: Experimental models have reported conflicting results regarding the role of dispersion of repolarization in promoting atrial fibrillation (AF). Repolarization alternans, a beat-to-beat alternation in action potential duration, enhances dispersion of repolarization when propagation velocity is involved. METHODS AND RESULTS: In this work, original electrophysiological parameters were analysed to study AF susceptibility in a chronic sheep model of pacing-induced AF. Two pacemakers were implanted, each with a single right atrial lead. Right atrial depolarization and repolarization waves were documented at 2-week intervals. A significant and gradual decrease in the propagation velocity at all pacing rates and in the right atrial effective refractory period (ERP) was observed during the weeks of burst pacing before sustained AF developed when compared with baseline conditions. Right atrial repolarization alternans was observed, but because of the development of 2/1 atrioventricular block with far-field ventricular interference, its threshold could not be precisely measured. Non-sustained AF was not observed at baseline, but appeared during the electrical remodelling in association with a decrease in both ERP and propagation velocity. CONCLUSION: We report here on the feasibility of measuring ERP, atrial repolarization alternans, and propagation velocity kinetics and their potential in predicting susceptibility to AF in a free-behaving model of pacing-induced AF using the standard pacemaker technology.

Phase-rectified signal averaging: a useful tool for the estimation of the dominant frequency in ECG signals during atrial fibrillation.

Atrial fibrillation (AF) is the most common type of human cardiac arrhythmia. An important parameter that can be extracted from surface electrocardiogram (ECG) during AF is the dominant frequency (DF) of AF. Unfortunately, AF signal components are always highly contaminated by the ventricular QRST complexes, and the cancellation of these components is never perfect. The remaining artifacts tend to induce DF overestimates. In this paper we report on the use of phase-rectified signal analysis, a technique introduced recently to enhance quasi-periodic signal components, for improving DF estimation. The potential of phase-rectified analysis is demonstrated through experiments both on synthetic and clinical ECG signals.