Search

We describe a novel approach for simultaneously determining regional differences in action potential (AP) morphology and tissue electrophysiological properties in isolated atria. The epicardial surface of rat atrial preparations was placed in contact with a multi-electrode array (9 x 10 silver chloride electrodes, 0.1 mm diameter and 0.1 mm pitch). A glass microelectrode (100 MH) was simultaneously inserted into the endocardial surface to record intracellular AP from either of 2 regions (A, B) during pacing from 2 opposite corners of the tissue. AP duration at 80% of repolarisation and its restitution curve was significantly different only in region A (p < 0.01) when AP was initiated at different stimulation sites. Alternans in AP duration and AP amplitude, and in conduction velocity were observed during 2 separate arrhythmic episodes. This approach of combining microelectrode array and intracellular membrane potential recording may provide new insights into arrhythmogenic mechanisms in animal models of cardiovascular disease. Key words: high-density mapping, atrial action potential membrane recordings, electrophysiology, arrhythmogenic mechanisms, atrial fibrillation. Nous decrivons une nouvelle approche visant a determiner les differences regionales entre les morphologies des potentiels d'action (PA) et les proprietes electrophysiologiques du tissu dans un modele d'oreillettes isolees. Nous avons mis la surface epicardique de preparations d'oreillettes de rat en contact avec un reseau d'electrodes multiples (9 rangees de 10 electrodes en chlorure d'argent de 0,1 mm de diametre a une profondeur de 0,1 mm). Simultanement, nous avons insere une electrode de verre (100 MH) a l'interieur de la surface endocardique pour enregistrer des PA intracellulaires dans l'une ou l'autre de 2 regions (A et B) pendant une stimulation electrique a partir de 2 coins opposes de l'echantillon de tissu. A 80 % de la repolarisation, nous notions que la duree du PA et sa courbe de restitution etaient differentes dans la region A (p < 0,01) uniquement lorsque le PA etait declenche a partir de sites de stimulation differents. Nous avons observe une alternance de la duree du PA et de son amplitude, ainsi que de la vitesse de conduction au cours de 2 periodes d'arythmie distinctes. Cette approche qui consiste a associer un reseau de microelectrodes a l'enregistrement de potentiels membranaires intracellulaires pourrait permettre d'obtenir de nouvelles donnees sur les mecanismes proarythmiques dans des modeles animaux de maladie cardiovasculaire. [Traduit par la Redaction] Mots-cles : cartographie a haute densite, enregistrements de potentiels d'action auriculaires, electrophysiologie, mecanismes proarythmiques, fibrillation auriculaire., Introduction To understand the genesis and maintenance of atrial arrhythmias, detailed evaluation of the structural and electrophysiological characteristics of atrial tissue is necessary. Electrophysiological techniques generally available are electro-anatomical mapping, [...]

Aims Cardiac arrhythmia originating from the papillary muscle (PM) can trigger ventricular fibrillation (VF) and cause sudden cardiac death even in the absence of structural heart disease. Most premature ventricular contractions, however, are benign and hitherto difficult to distinguish from a potentially fatal arrhythmia. Altered repolarization characteristics are associated with electrical instability, but electrophysiological changes which precede degeneration into VF are still not fully understood. Methods and results Ventricular arrhythmia (VA) was induced by aconitine injection into PMs of healthy sheep. To investigate mechanisms of degeneration of stable VA into VF in structurally healthy hearts, endocardial high-density and epicardial mapping was performed during sinus rhythm (SR) and VA. The electrical restitution curve, modelling the relation of diastolic interval and activation recovery interval (a surrogate parameter for action potential duration), is steeper in VA than in non-arrhythmia (ventricular pacing and SR). Steeper restitution curves reflect electrical instability and propensity to degenerate into VF. Importantly, we find the parameter repolarization time in relation to cycle length (RT/CL) to differentiate self-limiting from degenerating arrhythmia with high specificity and sensitivity. Conclusion RT/CL may serve as a simple index to aid differentiation between self-limiting and electrically instable arrhythmia with the propensity to degenerate to VF. RT/CL is independent of cycle length and could easily be measured to identify electrical instability in patients. [ABSTRACT FROM AUTHOR]

Background: Charge density mapping of atrial fibrillation (AF) reveals dynamic localized rotational activation (LRA), irregular activation (LIA) and focal firing (FF). Their spatial stability, conduction characteristics and the optimal duration of mapping required to reveal these phenomena and has not been explored. Methods: Bi‐atrial mapping of AF propagation was undertaken using AcQMap (Acutus Medical) and variability of activation patterns quantified up to a duration of 30 s. The frequency of each pattern was quantified at each unique point of the chamber over two separate 30‐s recordings before ablation and R2 calculated to quantify spatial stability. Regions with the highest frequency were identified at increasing time durations and compared to the result over 30 s using Cohen's kappa. Properties of regions with the most stable patterns were assessed during sinus rhythm and extrastimulus pacing. Results: In 21 patients, 62 paired LA and RA maps were obtained. LIA was highly spatially stable with R2 between maps of 0.83 (0.71–0.88) compared to 0.39 (0.24–0.57), and 0.64 (0.54–0.73) for LRA and FF, respectively. LIA was most temporally stable with a kappa of >0.8 reached by 12 s. LRA showed greatest variability with kappa >0.8 only after 22 s. Regions of LIA were of normal voltage amplitude (1.09 mv) but showed increased conduction heterogeneity during extrastimulus pacing (p =.0480). Conclusion: Irregular activation patterns characterized by changing wavefront direction are temporally and spatially stable in contrast with LRA that is transient with least spatial stability. Focal activation appears of intermediate stability. Regions of LIA show increased heterogeneity following extrastimulus pacing and may represent fixed anatomical substrate. [ABSTRACT FROM AUTHOR]

Background: Atrial fibrillation (AF) is characterized by the repetitive regeneration of unstable rotational events, the pivot of which are known as phase singularities (PSs). The spatial concentration and distribution of PSs have not been systematically investigated using quantitative statistical approaches. Objectives: We utilized a geospatial statistical approach to determine the presence of local spatial concentration and global clustering of PSs in biatrial human AF recordings. Methods: 64‐electrode conventional basket (~5 min, n = 18 patients, persistent AF) recordings were studied. Phase maps were produced using a Hilbert‐transform based approach. PSs were characterized spatially using the following approaches: (i) local "hotspots" of high phase singularity (PS) concentration using Getis‐Ord Gi* (Z ≥ 1.96, P ≤.05) and (ii) global spatial clustering using Moran's I (inverse distance matrix). Results: Episodes of AF were analyzed from basket catheter recordings (H: 41 epochs, 120 000 s, n = 18 patients). The Getis‐Ord Gi* statistic showed local PS hotspots in 12/41 basket recordings. As a metric of spatial clustering, Moran's I showed an overall mean of 0.033 (95% CI: 0.0003‐0.065), consistent with the notion of complete spatial randomness. Conclusion: Using a systematic, quantitative geospatial statistical approach, evidence for the existence of spatial concentrations ("hotspots") of PSs were detectable in human AF, along with evidence of spatial clustering. Geospatial statistical approaches offer a new approach to map and ablate PS clusters using substrate‐based approaches. [ABSTRACT FROM AUTHOR]

Background: Myocardial slow conduction is a cornerstone of ventricular tachycardia (VT). Prolonged electrogram (EGM) duration is a useful surrogate parameter and manual annotation of EGM characteristics are widely used during catheter-based ablation of the arrhythmogenic substrate. However, this remains time-consuming and prone to inter-operator variability. We aimed to develop an algorithm for 3-D visualization of EGM duration relative to the 17-segment American Heart Association model. Methods: To calculate and visualize EGM duration, in sinus rhythm acquired high-density maps of patients with ischemic cardiomyopathy undergoing substrate-based VT ablation using a 64-mini polar basket-catheter with low noise of 0.01 mV were analyzed. Using a custom developed algorithm based on standard deviation and threshold, the relationship between EGM duration, endocardial voltage and ablation areas was studied by creating 17-segment 3-D models and 2-D polar plots. Results: 140,508 EGMs from 272 segments (n = 16 patients, 94% male, age: 66±2.4, ejection fraction: 31±2%) were studied and 3-D visualization of EGM duration was performed. Analysis of signal processing parameters revealed that a 40 ms sliding SD-window, 15% SD-threshold and >70 ms EGM duration cutoff was chosen based on diagnostic odds ratio of 12.77 to visualize rapidly prolonged EGM durations. EGMs > 70 ms matched to 99% of areas within dense scar (<0.2 mV), in 95% of zones within scar border zone (0.2–1.0 mV) and detected ablated areas having resulted in non-inducibility at the end of the procedure. Ablation targets were identified with a sensitivity of 65.6% and a specificity of 94.6% avoiding false positive labeling of prolonged EGMs in segments with healthy myocardium. Conclusion: The novel algorithm allows rapid visualization of prolonged EGM durations. This may facilitate more objective characterization of arrhythmogenic substrate in patients with ischemic cardiomyopathy. [ABSTRACT FROM AUTHOR]

Aims: Complex propagation patterns are observed in patients and models with stable atrial fibrillation (AF). The degree of this complexity is associated with AF stability. Experimental work suggests reduced wavefront turning as an important mechanism for widening of the excitable gap. The aim of this study was to investigate how sodium channel inhibition by vernakalant affects turning behaviour and propagation patterns during AF.Methods and Results: Two groups of 8 goats were instrumented with electrodes on the left atrium, and AF was maintained by burst pacing for 3 or 22 weeks. Measurements were performed at baseline and two dosages of vernakalant. Unipolar electrograms were mapped (249 electrodes/array) on the left and right atrium in an open-chest experiment. Local activation times and conduction vectors, flow lines, the number of fibrillation waves, and local re-entries were determined. At baseline, fibrillation patterns contained numerous individual fibrillation waves conducting in random directions. Vernakalant induced conduction slowing and cycle length prolongation and terminated AF in 13/15 goats. Local re-entries were strongly reduced. Local conduction vectors showed increased preferential directions and less beat-to-beat variability. Breakthroughs and waves were significantly reduced in number. Flow line curvature reduced and waves conducted more homogenously in one direction. Overall, complex propagation patterns were strongly reduced. No substantial differences in drug effects between right and left atria or between goats with different AF durations were observed.Conclusions: Destabilization of AF by vernakalant is associated with a lowering of fibrillation frequency and inhibition of complex propagation patterns, wave turning, local re-entries, and breakthroughs. [ABSTRACT FROM AUTHOR]

Rationale: A quantitative framework to summarize and explain the quasi-stationary population dynamics of unstable phase singularities (PS) and wavelets in human atrial fibrillation (AF) is at present lacking. Building on recent evidence showing that the formation and destruction of PS and wavelets in AF can be represented as renewal processes, we sought to establish such a quantitative framework, which could also potentially provide insight into the mechanisms of spontaneous AF termination. Objectives: Here, we hypothesized that the observed number of PS or wavelets in AF could be governed by a common set of renewal rate constants λ f (for PS or wavelet formation) and λ d (PS or wavelet destruction), with steady-state population dynamics modeled as an M/M/∞ birth–death process. We further hypothesized that changes to the M/M/∞ birth–death matrix would explain spontaneous AF termination. Methods and Results: AF was studied in in a multimodality, multispecies study in humans, animal experimental models (rats and sheep) and Ramirez-Nattel-Courtemanche model computer simulations. We demonstrated: (i) that λ f and λ d can be combined in a Markov M/M/∞ process to accurately model the observed average number and population distribution of PS and wavelets in all systems at different scales of mapping; and (ii) that slowing of the rate constants λ f and λ d is associated with slower mixing rates of the M/M/∞ birth–death matrix, providing an explanation for spontaneous AF termination. Conclusion: M/M/∞ birth–death processes provide an accurate quantitative representational architecture to characterize PS and wavelet population dynamics in AF, by providing governing equations to understand the regeneration of PS and wavelets during sustained AF, as well as providing insight into the mechanism of spontaneous AF termination. [ABSTRACT FROM AUTHOR]

Introduction: Respiratory sinus arrhythmia (RSA) describes heart rate (HR) changes in synchrony with respiration. It is relevant for exercise capacity and mechanistically linked with the cardiac autonomic nervous system. After pulmonary vein isolation (PVI), the current therapy of choice for patients with paroxysmal atrial fibrillation (AF), the cardiac vagal tone is often diminished. We hypothesized that RSA is modulated by PVI in patients with paroxysmal AF.Material and Methods: Respiratory sinus arrhythmia, measured by using a deep breathing test and heart rate variability parameters, was studied in 10 patients (64 ±3 years) with paroxysmal AF presenting in stable sinus rhythm for their first catheter-based PVI. Additionally, heart rate dynamics before and after PVI were studied during sympathetic/parasympathetic coactivation by using a cold-face test. All tests were performed within 24 h before and 48 h after PVI.Results: After PVI RSA (E/I difference: 7.9 ±1.0 vs. 3.5 ±0.6 bpm, p = 0.006; E/I ratio: 1.14 ±0.02 vs. 1.05 ±0.01, p = 0.003), heart rate variability (SDNN: 31 ±3 vs. 14 ±3 ms, p = 0.006; RMSSD: 17 ±2 vs. 8 ±2 ms, p = 0.002) and the HR response to sympathetic/parasympathetic coactivation (10.2 ±0.7% vs. 5.7 ±1.1%, p = 0.014) were diminished. The PVI-related changes in RSA correlated with the heart rate change during sympathetic/parasympathetic coactivation before vs. after PVI (E/I difference: r = 0.849, p = 0.002; E/I ratio: r = 0.786, p = 0.007). One patient with vagal driven arrhythmia experienced AF recurrence during follow-up (mean: 6.5 ±0.6 months).Conclusions: Respiratory sinus arrhythmia is reduced after PVI in patients with paroxysmal AF. Our findings suggest that this is related to a decrease in cardiac vagal tone. Whether and how this affects the clinical outcome including exercise capacity need to be determined. [ABSTRACT FROM AUTHOR]

Background: The potential utility of entropy (En) for atrial fibrillation (AF) mapping has been demonstrated in previous studies by multiple groups, where an association between high bipolar electrogram (EGM) entropy and the pivot of rotors has been shown. Though En is potentially attractive new approach to ablation, no studies have examined its temporal stability and specificity, which are critical to the application of entropy to clinical ablation. In the current study, we sought to objectively measure the temporal stability and specificity of bipolar EGM entropy in medium to long term recordings using three studies: i) a human basket catheter AF study, ii) a tachypaced sheep AF study and iii) a computer simulation study.Objective: To characterize the temporal dynamics and specificity of Approximate, Sample and Shannon entropy (ApEn/SampEn/ShEn) in human (H), sheep (S), and computer simulated AF.Methods: 64-electrode basket bi-atria sustained AF recordings (H:15 min; S:40 min) were separated into 5 s segments. ShEn/ApEn/SampEn were computed, and co-registered with NavX 3D maps. Temporal stability was determined in terms of: (i) global pattern stability of En and (ii) the relative stability the top 10% of En regions. To provide mechanistic insights into underlying mechanisms, stability characteristics were compared to models depicting various propagation patterns. To verify these results, cross-validation was performed across multiple En algorithms, across species, and compared with dominant frequency (DF) temporal characteristics. The specificity of En was also determined by looking at the association of En to rotors and areas of wave cross propagation.Results: Episodes of AF were analysed (H:26 epochs, 6040 s; S:15 epochs, 14,160 s). The global pattern of En was temporally unstable (CV- H:13.42% ± 4.58%; S:14.13% ± 8.13%; Friedman- H: p > 0.001; S: p > 0.001). However, within this dynamic flux, the top 10% of ApEn/SampEn/ShEn regions were relatively temporally stable (Kappa >0.6) whilst the top 10% of DF regions were unstable (Kappa <0.06). In simulated AF scenarios, the experimental data were optimally reproduced in the context of an AF pattern with stable rotating waves surrounded by wavelet breakup (Kappa: 0.610; p < 0.0001).Conclusion: En shows global temporal instability, however within this dynamic flux, the top 10% regions exhibited relative temporal stability. This suggests that high En regions may be an appealing ablation target. Despite this, high En was associated with not just the pivot of rotors but also with areas of cross propagation, which suggests the need for future work before clinical application is possible. [ABSTRACT FROM AUTHOR]

Background: QRS abnormalities may not be apparent in sinus rhythm in electrically stable cardiomyopathy patients who can have quiescent but highly arrhythmogenic substrate. Here, we test the hypothesis that differential changes in QRS construction during right‐ventricular apex pacing (RVP) as opposed to atrial pacing (AP) will identify latent substrate for ventricular arrhythmias (VA) and death. Methods: Forty patients with cardiomyopathy free of VA underwent baseline 114‐electrode body‐surface electrocardiogram during AP (100 beats per minute [bpm]) and RVP (100 and 120 bpm). The filtered‐averaged QRS at each electrode was deconstructed into individual intra‐QRS and post‐QRS ventricular myopotentials (VMP). The primary outcome was VA or death. Prognostic accuracy of VMP was validated using V1 to V6 leads in another prospective cohort of 44‐cardiomyopathy patients. Results: Twenty‐six patients were eligible for initial analysis. After 5 ± 2 years of follow‐up, eight (31%) patients had VA (VAPos) while rest were uneventful (VANeg). During AP100, VAPos patients expressed more VMP than VANeg patients (16 ± 1 vs 12 ± 1, P = 0.02). RVP100 and RVP120 in VAPos patients introduced an additional 5.5 ± 0.5 and 6.0 ± 0.5 VMP (P < 0.0001 vs AP100). The relative change with RVP120 versus AP100 in VANeg patients exceeded VAPos patients by 1.2 ± 0.5 VMP (P = 0.03). Increment in VMP count of <8 in lead‐V5 with RVP120 compared to AP100 best predicted VA (area under curve 0.81, P = 0.01). In the validation cohort, primary outcome occurred in 13 (33%) patients. Native QRS features and AP100 alone failed to predict primary outcome. Patients with increment in VMP count of <8 in lead‐V5 with RVP120 versus AP100 had 7.9‐fold increased risk of primary outcome (95% confidence interval 1.01, 61.61; P = 0.049). Conclusion: Cardiomyopathy patients at risk of VA or death perturb the QRS less than low‐risk patients with differential pacing. This functional response may be useful to identify arrhythmogenic substrate. [ABSTRACT FROM AUTHOR]

Atrial Fibrillation (AF) is the most common cardiac rhythm disorder seen in hospitals and in general practice, accounting for up to a third of arrhythmia related hospitalizations. Unfortunately, AF treatment is in practice complicated by the lack of understanding of the fundamental mechanisms underlying the arrhythmia, which makes detection of effective ablation targets particularly difficult. Various approaches to AF mapping have been explored in the hopes of better pinpointing these effective targets, such as Dominant Frequency (DF) analysis, complex fractionated electrograms (CFAE) and unipolar reconstruction (FIRM), but many of these methods have produced conflicting results or require further investigation. Exploration of AF using information theoretic-based approaches may have the potential to provide new insights into the complex system dynamics of AF, whilst also providing the benefit of being less reliant on empirically derived definitions in comparison to alternate mapping approaches. This work provides an overview of information theory and reviews its applications in AF analysis, with particular focus on AF mapping. The works discussed in this review demonstrate how understanding AF from a signal property perspective can provide new insights into the arrhythmic phenomena, which may have valuable clinical implications for AF mapping and ablation in the future. [ABSTRACT FROM AUTHOR]

Background: Current phase mapping systems for atrial fibrillation create 2-dimensional (2D) maps. This process may affect the accurate detection of rotors. We developed a 3-dimensional (3D) phase mapping technique that uses the 3D locations of basket electrodes to project phase onto patient-specific left atrial 3D surface anatomy.Objective: We sought to determine whether rotors detected in 2D phase maps were present at the corresponding time segments and anatomical locations in 3D phase maps.Methods: One-minute left atrial atrial fibrillation recordings were obtained in 14 patients using the basket catheter and analyzed off-line. Using the same phase values, 2D and 3D phase maps were created. Analysis involved determining the dominant propagation patterns in 2D phase maps and evaluating the presence of rotors detected in 2D phase maps in the corresponding 3D phase maps.Results: Using 2D phase mapping, the dominant propagation pattern was single wavefront (36.6%) followed by focal activation (34.0%), disorganized activity (23.7%), rotors (3.3%), and multiple wavefronts (2.4%). Ten transient rotors were observed in 9 of 14 patients (64%). The mean rotor duration was 1.1 ± 0.7 seconds. None of the 10 rotors observed in 2D phase maps were seen at the corresponding time segments and anatomical locations in 3D phase maps; 4 of 10 corresponded with single wavefronts in 3D phase maps, 2 of 10 with 2 simultaneous wavefronts, 1 of 10 with disorganized activity, and in 3 of 10 there was no coverage by the basket catheter at the corresponding 3D anatomical location.Conclusion: Rotors detected in 2D phase maps were not observed in the corresponding 3D phase maps. These findings may have implications for current systems that use 2D phase mapping. [ABSTRACT FROM AUTHOR]

Purpose: Protected channels of surviving myocytes in late postinfarction ventricular scar predispose to ventricular tachycardia (VT). However, only a few patients develop VT spontaneously. We studied differences in electric remodeling and protected channels in late postinfarction patients with and without spontaneous VT.Methods: Patients with ischemic cardiomyopathy (ICM) with recurrent sustained monomorphic VT (n = 22) were compared with stable ICM patients without spontaneous VT (control group; n = 5). Left ventricular mapping was performed with a 20-pole catheter. Detailed pace mapping was used to identify channels of protected conduction, and confirmed, when feasible, by entrainment. Anatomical and electrophysiological properties of VT channels and non-VT channels in VT patients and channels in controls were evaluated.Results: Seventy-three (median 3) VTs were inducible in VT patients compared to two (median 0) in controls. The VT channels in VT patients (n = 57, 3 ± 1 per patient) were lengthier (mean ± SEM 53 ± 5 vs. 33 ± 4 vs. 24 ± 8 mm), had longer S-QRS (73 ± 4 vs. 63 ± 3 vs. 44 ± 8 ms), longer conduction time (103 ± 13 vs. 33 ± 4 vs. 24 ± 8 ms), and slower conduction velocity (CV) (0.85 ± 0.21 vs. 1.39 ± 0.20 vs. 1.31 ± 0.41 m/s) than non-VT channels in VT patients (n = 183, 8 ± 6 per patient) (p ≤ 0.01) and channels in controls (n = 46, 9 ± 8 per patient) (p ≤ 0.01). Additionally, non-VT channels in VT patients had longer S-QRS (p = 0.02); however, they were similar in length, conduction time, and CV compared to channels in controls.Conclusions: Channels supporting VT are lengthier, with longer conduction times and slower CV compared to channels in patients without spontaneous VT. These observations may explain why some ICM patients have spontaneous VT and others do not. [ABSTRACT FROM AUTHOR]

Objective: To explore technical challenges of phase singularity (PS) mapping during atrial fibrillation (AF) using direct contact electrograms. Methods: AF mapping was performed in high-density epicardial recordings of human paroxysmal (PAF) or persistent (PersAF) (N = 20 pts) AF with an array of 16 × 16 electrodes placed on atrial epicardium. PS points were detected using subsets of electrodes forming rings of varying sizes. Results: PS detection using a 2 × 2 electrode ring identified 0.88 ± 1.00 PS/s in PAF group and 3.91 ± 2.51 per s in PersAF group (p < 0.001) in 2.4 × 2.4 cm mapping area. All detected PS had a short lifespan with the longest being 1100 ms (6.8 rotations). Exploration of the PS detection in a numerical model demonstrated that at least eight electrodes are required to avoid frequent false positive PS detection due to chance. Application of a detection grid consisting a double ring of electrodes (2 × 2 and 4 × 4 rings) decreased the number of false positive detections. The double ring was more resilient to electrode swapping (with just three instances of false positives versus 4380 false positives using 2 × 2 ring). Conclusions: The number of detected rotors critically depends upon the parameters of the detection algorithm, especially the number of electrodes used to detect PS. Based on our results, we recommend double ring comprised of 2 × 2 and 4 × 4 grid of electrodes for robust rotor detection. Significance: Great methodological care has to be taken before equating detected PS with rotating waves and using PS detection algorithms to guide catheter ablation of AF. [ABSTRACT FROM AUTHOR]

Aims: Loss of side-to-side electrical connections between atrial muscle bundles is thought to underlie conduction disturbances predisposing to atrial fibrillation (AF). Putatively, disruption of electrical connections occurs not only within the epicardial layer but also between the epicardial layer and the endocardial bundle network, thus impeding transmural conductions (‘breakthroughs’). However, both clinical and experimental studies have shown an enhancement of breakthroughs during later stages of AF. We tested the hypothesis that endo-epicardial uncoupling enhances endo-epicardial electrical dyssynchrony, breakthrough rate (BTR), and AF stability.Methods and Results: In a novel dual-layer computer model of the human atria, 100% connectivity between the two layers served as healthy control. Atrial structural remodelling was simulated by reducing the number of connections between the layers from 96 to 6 randomly chosen locations. With progressive elimination of connections, AF stability increased. Reduction in the number of connections from 96 to 24 resulted in an increase in endo-epicardial dyssynchrony from 6.6 ± 1.9 to 24.6 ± 1.3%, with a concomitant increase in BTR. A further reduction to 12 and 6 resulted in more pronounced endo-epicardial dyssynchrony of 34.4 ± 1.15 and 40.2 ± 0.52% but with BTR reduction. This biphasic relationship between endo-epicardial coupling and BTR was found independently from whether AF was maintained by re-entry or by ectopic focal discharges.Conclusion: Loss of endo-epicardial coupling increases AF stability. There is a biphasic relation between endo-epicardial coupling and BTR. While at high degrees of endo-epicardial connectivity, the BTR is limited by the endo-epicardial synchronicity, at low degrees of connectivity, it is limited by the number of endo-epicardial connections. [ABSTRACT FROM AUTHOR]