Human ventricular fibrillation during global ischemia and reperfusion: paradoxical changes in activation rate and wavefront complexity

Sodium current preserves electrical excitability in the heart of hibernating ground squirrel (Citellus undulatus)

Hibernating mammals are capable of maintaining normal cardiac function at low temperatures. Excitability of cardiac myocytes crucially depends on the fast sodium current (I_Na), which is decreased in hypothermia due to both depolarization of resting membrane potential and direct negative effect of low temperature. Therefore, I_Na in hibernating mammals should have specific features allowing to maintain excitability of myocardium at low temperatures. The current-voltage dependence of I_Na, its steady-state inactivation and activation and recovery from inactivation were studied in winter hibernating (WH) and summer active (SA) ground squirrels and in rats using whole-cell patch clamp at 10 °C and 20 °C. I_Na peak amplitude and the parameters of steady-state activation and inactivation curves did not differ between SA and WH ground squirrels at both temperatures. However, at both temperatures strong positive shift of activation and inactivation curves by 5-12 mV was observed in both WH and SA ground squirrels if compared to rats. This peculiarity of cardiac I_Na in ground squirrels helps to maintain excitability in conditions of depolarized resting membrane potential. The time course of I_Na recovery from inactivation at 10 °C was faster in WH than in SA ground squirrels, which could ensure normal activation of myocardium during hibernation.

Topological charge-density method of identifying phase singularities in cardiac fibrillation

Spiral waves represent the key motifs of typical self-sustained dynamical patterns in excitable systems such as cardiac tissue. The motion of phase singularities (PSs) that lies at the center of spiral waves captures many qualitative and, in some cases, quantitative features of their complex dynamics. Recent clinical studies suggested that ablating the tissue at PS locations may cure atrial fibrillation. Here, we propose a different method to determine the location of PSs. Starting from the definition of the topological charge of spiral waves, we obtain the expression of the topological charge density in a discrete case. With this expression, we can calculate the topological charge at each grid in the space directly, so as to accurately identify the position of PSs.

M/M/Infinity Birth-Death Processes - A Quantitative Representational Framework to Summarize and Explain Phase Singularity and Wavelet Dynamics in Atrial Fibrillation

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.

Insights Into the Spatiotemporal Patterns of Complexity of Ventricular Fibrillation by Multilead Analysis of Body Surface Potential Maps

BACKGROUND

Ventricular fibrillation (VF) is the main cause of sudden cardiac death, but its mechanisms are still unclear. We propose a noninvasive approach to describe the progression of VF complexity from body surface potential maps (BSPMs).

METHODS

We mapped 252 VF episodes (16 ± 10 s) with a 252-electrode vest in 110 patients (89 male, 47 ± 18 years): 50 terminated spontaneously, otherwise by electrical cardioversion (DCC). Changes in complexity were assessed between the onset ("VF start") and the end ("VF end") of VF by the nondipolar component index (N D I B S P M ), measuring the fraction of energy nonpreserved by an equivalent 3D dipole from BSPMs. Higher NDI reflected lower VF organization. We also examined other standard body surface markers of VF dynamics, including fibrillatory wave amplitude (A BSPM ), surface cycle length (BsCL BSPM ) and Shannon entropy (S h E n B S P M ). Differences between patients with and without structural heart diseases (SHD, 32 vs. NSHD, 78) were also tested at those stages. Electrocardiographic features were validated with simultaneous endocardium cycle length (CL) in a subset of 30 patients.

RESULTS

All BSPM markers measure an increase in electrical complexity during VF (p < 0.0001), and more significantly in NSHD patients. Complexity is significantly higher at the end of sustained VF episodes requiring DCC. Intraepisode intracardiac CL shortening (VF start 197 ± 24 vs. VF end 169 ± 20 ms; p < 0.0001) correlates with an increase in NDI, and decline in surface CL, f-wave amplitude, and entropy (p < 0.0001). In SHD patients VF is initially more complex than in NSHD patients (N D I B S P M , p = 0.0007; S h E n B S P M , p < 0.0001), with moderately slower (BsCL BSPM , p = 0.06), low-amplitude f-waves (A BSPM , p < 0.0001). In this population, lower NDI (p = 0.004) and slower surface CL (p = 0.008) at early stage of VF predict self-termination. In the NSHD group, a more abrupt increase in VF complexity is quantified by all BSPM parameters during sustained VF (p < 0.0001), whereas arrhythmia evolution is stable during self-terminating episodes, hinting at additional mechanisms driving VF dynamics.

CONCLUSION

Multilead BSPM analysis underlines distinct degrees of VF complexity based on substrate characteristics.

Unraveling the Underlying Arrhythmia Mechanism in Persistent Atrial Fibrillation: Results From the STARLIGHT Study

BACKGROUND

The mechanisms that initiate and sustain persistent atrial fibrillation are not well characterized. Ablation results remain significantly worse than in paroxysmal atrial fibrillation in which the mechanism is better understood and subsequent targeted therapy has been developed. The aim of this study was to characterize and quantify patterns of activation during atrial fibrillation using contact mapping.

METHODS

Patients with persistent atrial fibrillation (n=14; mean age, 61±8 years; ejection fraction, 59±10%) underwent simultaneous biatrial contact mapping with 64 electrode catheters. The atrial electrograms were transformed into phase, and subsequent spatiotemporal mapping was performed to identify phase singularities (PSs).

RESULTS

PSs were located in both atria, but we observed more PSs in the left atrium compared with the right atrium (779±302, 552±235; P=0.015). Although some PSs of duration sufficient to complete >1 rotation were detected, the maximum PS duration was only 1150 ms, and the vast majority (97%) of PSs persisted for too short a period to complete a full rotation. Although in selected patients there was evidence of PS local clustering, overall, PSs were distributed globally throughout both chambers with no clear anatomic predisposition. In a subset of patients (n=7), analysis was repeated using an alternative established atrial PS mapping technique, which confirmed our initial findings.

CONCLUSIONS

No sustained rotors or localized drivers were detected, and instead, the mechanism of arrhythmia maintenance was consistent with the multiple wavelet hypothesis, with passive activation of short-lived rotational activity.

CLINICAL TRIAL REGISTRATION

URL: https://www.clinicaltrials.gov. Unique identifier: NCT01765075.

Dynamical anchoring of distant arrhythmia sources by fibrotic regions via restructuring of the activation pattern

Rotors are functional reentry sources identified in clinically relevant cardiac arrhythmias, such as ventricular and atrial fibrillation. Ablation targeting rotor sites has resulted in arrhythmia termination. Recent clinical, experimental and modelling studies demonstrate that rotors are often anchored around fibrotic scars or regions with increased fibrosis. However, the mechanisms leading to abundance of rotors at these locations are not clear. The current study explores the hypothesis whether fibrotic scars just serve as anchoring sites for the rotors or whether there are other active processes which drive the rotors to these fibrotic regions. Rotors were induced at different distances from fibrotic scars of various sizes and degree of fibrosis. Simulations were performed in a 2D model of human ventricular tissue and in a patient-specific model of the left ventricle of a patient with remote myocardial infarction. In both the 2D and the patient-specific model we found that without fibrotic scars, the rotors were stable at the site of their initiation. However, in the presence of a scar, rotors were eventually dynamically anchored from large distances by the fibrotic scar via a process of dynamical reorganization of the excitation pattern. This process coalesces with a change from polymorphic to monomorphic ventricular tachycardia.

Rotors are waves of cardiac excitation like a tornado causing cardiac arrhythmia. Recent research shows that they are found in ventricular and atrial fibrillation. Burning (via ablation) the site of a rotor can result in the termination of the arrhythmia. Recent studies showed that rotors are often anchored to regions surrounding scar tissue, where part of the tissue still survived called fibrotic tissue. However, it is unclear why these rotors anchor to these locations. Therefore, in this work, we investigated why rotors are so abundant in fibrotic tissue with the help of computer simulations. We performed simulations in a 2D model of human ventricular tissue and in a patient-specific model of a patient with an infarction. We found that even when rotors are initially at large distances from the fibrotic region, they are attracted by this region, to finally end up at the fibrotic tissue. We called this process dynamical anchoring and explained how the process works.

Analytical approaches for myocardial fibrillation signals

Atrial and ventricular fibrillation are complex arrhythmias, and their underlying mechanisms remain widely debated and incompletely understood. This is partly because the electrical signals recorded during myocardial fibrillation are themselves complex and difficult to interpret with simple analytical tools. There are currently a number of analytical approaches to handle fibrillation data. Some of these techniques focus on mapping putative drivers of myocardial fibrillation, such as dominant frequency, organizational index, Shannon entropy and phase mapping. Other techniques focus on mapping the underlying myocardial substrate sustaining fibrillation, such as voltage mapping and complex fractionated electrogram mapping. In this review, we discuss these techniques, their application and their limitations, with reference to our experimental and clinical data. We also describe novel tools including a new algorithm to map microreentrant circuits sustaining fibrillation.

Spatio-temporal Organization During Ventricular Fibrillation in the Human Heart

In this paper, we present a novel approach to quantify the spatio-temporal organization of electrical activation during human ventricular fibrillation (VF). We propose three different methods based on correlation analysis, graph theoretical measures and hierarchical clustering. Using the proposed approach, we quantified the level of spatio-temporal organization during three episodes of VF in ten patients, recorded using multi-electrode epicardial recordings with 30 s coronary perfusion, 150 s global myocardial ischaemia and 30 s reflow. Our findings show a steady decline in spatio-temporal organization from the onset of VF with coronary perfusion. We observed transient increases in spatio-temporal organization during global myocardial ischaemia. However, the decline in spatio-temporal organization continued during reflow. Our results were consistent across all patients, and were consistent with the numbers of phase singularities. Our findings show that the complex spatio-temporal patterns can be studied using complex network analysis.

Spatiotemporal Progression of Early Human Ventricular Fibrillation

Objectives

The objective of this study was to evaluate the spatio-temporal organization and progression of human ventricular fibrillation (VF) in the left (LV) and right (RV) ventricles.

Background

Studies suggest that localized sources contribute to VF maintenance, but the evolution of VF episodes has not been quantified.

Methods

Synchrony between electrograms recorded from 25 patients with induced VF is computed and used to define the Asynchronous Index (ASI), indicating regions which are out-of-step with surrounding tissue. Computer simulations show that ASI can identify the location of VF-maintaining sources, where larger values of ASI_max correlate with more stable sources.

Results

Automated synchrony analysis shows elevated values of ASI in a majority of self-terminating episodes (LV: 8/9, RV: 7/8) and sustained episodes (LV: 11/11, RV: 12/12). The locations of ASI_max in sustained episodes co-localize with rotor cores when rotational activity is simultaneously present in phase maps (LV: 8/8, RV: 5/7, p<.05). The distribution of ASI_max differentiates self-terminating from sustained episodes (mean ASI_max = 0.60±0.14 and 0.70±0.16, respectively; p=0.01). Across sustained episodes the LV exhibits an increase in ASI_max with time.

Conclusions

Quantitative analysis identifies localized asynchronous regions that correlate with sources in VF, with sustained episodes evolving to exhibit more stable activation in the LV. This successive increase in stability indicates a stabilizing agent may be responsible for perpetuating fibrillation in a "migrate-and-capture" mechanism in the LV.

Autophagy: an adaptive physiological countermeasure to cellular senescence and ischaemia/reperfusion-associated cardiac arrhythmias

Oxidative stress placed on tissues that involved in pathogenesis of a disease activates compensatory metabolic changes, such as DNA damage repair that in turn causes intracellular accumulation of detritus and ‘proteotoxic stress’, leading to emergence of ‘senescent’ cellular phenotypes, which express high levels of inflammatory mediators, resulting in degradation of tissue function. Proteotoxic stress resulting from hyperactive inflammation following reperfusion of ischaemic tissue causes accumulation of proteinaceous debris in cells of the heart in ways that cause potentially fatal arrhythmias, in particular ventricular fibrillation (VF). An adaptive response to VF is occurrence of autophagy, an intracellular bulk degradation of damaged macromolecules and organelles that may restore cellular and tissue homoeostasis, improving chances for recovery. Nevertheless, depending on the type and intensity of stressors and inflammatory responses, autophagy may become pathological, resulting in excessive cell death. The present review examines the multilayered defences that cells have evolved to reduce proteotoxic stress by degradation of potentially toxic material beginning with endoplasmic reticulum‐associated degradation, and the unfolded protein response, which are mechanisms for removal from the endoplasmic reticulum of misfolded proteins, and then progressing through the stages of autophagy, including descriptions of autophagosomes and related vesicular structures which process material for degradation and autophagy‐associated proteins including Beclin‐1 and regulatory complexes. The physiological roles of each mode of proteotoxic defence will be examined along with consideration of how emerging understanding of autophagy, along with a newly discovered regulatory cell type called telocytes, may be used to augment existing strategies for the prevention and management of cardiovascular disease.

Human-based approaches to pharmacology and cardiology: an interdisciplinary and intersectorial workshop

Both biomedical research and clinical practice rely on complex datasets for the physiological and genetic characterization of human hearts in health and disease. Given the complexity and variety of approaches and recordings, there is now growing recognition of the need to embed computational methods in cardiovascular medicine and science for analysis, integration and prediction. This paper describes a Workshop on Computational Cardiovascular Science that created an international, interdisciplinary and inter-sectorial forum to define the next steps for a human-based approach to disease supported by computational methodologies. The main ideas highlighted were (i) a shift towards human-based methodologies, spurred by advances in new in silico, in vivo, in vitro, and ex vivo techniques and the increasing acknowledgement of the limitations of animal models. (ii) Computational approaches complement, expand, bridge, and integrate in vitro, in vivo, and ex vivo experimental and clinical data and methods, and as such they are an integral part of human-based methodologies in pharmacology and medicine. (iii) The effective implementation of multi- and interdisciplinary approaches, teams, and training combining and integrating computational methods with experimental and clinical approaches across academia, industry, and healthcare settings is a priority. (iv) The human-based cross-disciplinary approach requires experts in specific methodologies and domains, who also have the capacity to communicate and collaborate across disciplines and cross-sector environments. (v) This new translational domain for human-based cardiology and pharmacology requires new partnerships supported financially and institutionally across sectors. Institutional, organizational, and social barriers must be identified, understood and overcome in each specific setting.

Ventricular fibrillation: triggers, mechanisms and therapies

Ventricular fibrillation (VF) is a common, life-threatening arrhythmia responsible for significant morbidity and mortality. Due to challenges in safely mapping VF, a comprehensive understanding of its mechanisms remains elusive. Recent findings have provided new insights into mechanisms that sustain early VF. Notably, the central role of electrical rotors and catheter-based ablation of VF rotor substrate have been recently reported. In this article, we will review data regarding four stages of VF: initiation, transition, maintenance and evolution. We will discuss the particular mechanisms for each stage and therapies targeting these mechanisms. We also examine inherited arrhythmia syndromes, including the mechanisms and therapies specific to each. We hope that the overview of VF outlined in this work will assist other investigators in designing future therapies to interrupt this life-threatening arrhythmia.

Acute Stress Decreases but Chronic Stress Increases Myocardial Sensitivity to Ischemic Injury in Rodents

Cardiovascular disease (CVD) is the largest cause of mortality worldwide, and stress is a significant contributor to the development of CVD. The relationship between acute and chronic stress and CVD is well evidenced. Acute stress can lead to arrhythmias and ischemic injury. However, recent evidence in rodent models suggests that acute stress can decrease sensitivity to myocardial ischemia-reperfusion injury (IRI). Conversely, chronic stress is arrhythmogenic and increases sensitivity to myocardial IRI. Few studies have examined the impact of validated animal models of stress-related psychological disorders on the ischemic heart. This review examines the work that has been completed using rat models to study the effects of stress on myocardial sensitivity to ischemic injury. Utilization of animal models of stress-related psychological disorders is critical in the prevention and treatment of cardiovascular disorders in patients experiencing stress-related psychiatric conditions.

Models of ventricular arrhythmia mechanisms

The mechanisms that initiate and sustain ventricular arrhythmias in the human heart are clinically important, but hard to study experimentally. In this study, a monodomain model of electrical activation was used to examine how dynamics of electrophysiology at the cell scale influence the surface activation patterns of VF at the tissue scale. Cellular electrophysiology was described with two variants of a phenomenological model of the human ventricular epicardial action potential. The tissue geometry was an 8.0 × 8.0 × 1.2 cm 3D tissue slab with axially symmetric anisotropy. In both cases an initial re-entrant wave fragmented into multiple wavelets of activation. The model variant with steep action potential duration restitution produced much more complex activation, with a greater average number of filaments (13.79) than the variant with less steep restitution (3.08). More complex activation was associated with proportionally fewer transmural filaments, and so the average number of epicardial wavefronts and phase singularities per filament was lower. The average number of epicardial phase singularities and wavefronts for the model variant with less steep restitution were consistent with experimental observations in the human heart. This study shows that small changes in cell scale dynamics can have a large influence on the complexity of re-entrant activation in simulated 3D tissue, as well as on the features observed on the epicardial surface.

Ventricular Tachycardia and Early Fibrillation in Patients With Brugada Syndrome and Ischemic Cardiomyopathy Show Predictable Frequency-Phase Properties on the Precordial ECG Consistent With the Respective Arrhythmogenic Substrate

Supplemental Digital Content is available in the text.

Ventricular fibrillation (VF) has been proposed to be maintained by localized high-frequency sources. We tested whether spectral-phase analysis of the precordial ECG enabled identification of periodic activation patterns generated by such sources.

Analysis of cardiac fibrillation using phase mapping

The sequence of myocardial electrical activation during fibrillation is complex and changes with each cycle. Phase analysis represents the electrical activation-recovery process as an angle. Lines of equal phase converge at a phase singularity at the center of rotation of a reentrant wave, and the identification of reentry and tracking of reentrant wavefronts can be automated. We examine the basic ideas behind phase analysis. With the exciting prospect of using phase analysis of atrial electrograms to guide ablation in the human heart, we highlight several recent developments in preprocessing electrograms so that phase can be estimated reliably.

Effect of global cardiac ischemia on human ventricular fibrillation: insights from a multi-scale mechanistic model of the human heart

Acute regional ischemia in the heart can lead to cardiac arrhythmias such as ventricular fibrillation (VF), which in turn compromise cardiac output and result in secondary global cardiac ischemia. The secondary ischemia may influence the underlying arrhythmia mechanism. A recent clinical study documents the effect of global cardiac ischaemia on the mechanisms of VF. During 150 seconds of global ischemia the dominant frequency of activation decreased, while after reperfusion it increased rapidly. At the same time the complexity of epicardial excitation, measured as the number of epicardical phase singularity points, remained approximately constant during ischemia. Here we perform numerical studies based on these clinical data and propose explanations for the observed dynamics of the period and complexity of activation patterns. In particular, we study the effects on ischemia in pseudo-1D and 2D cardiac tissue models as well as in an anatomically accurate model of human heart ventricles. We demonstrate that the fall of dominant frequency in VF during secondary ischemia can be explained by an increase in extracellular potassium, while the increase during reperfusion is consistent with washout of potassium and continued activation of the ATP-dependent potassium channels. We also suggest that memory effects are responsible for the observed complexity dynamics. In addition, we present unpublished clinical results of individual patient recordings and propose a way of estimating extracellular potassium and activation of ATP-dependent potassium channels from these measurements.

Rotor stability separates sustained ventricular fibrillation from self-terminating episodes in humans

OBJECTIVES

This study mapped human ventricular fibrillation (VF) to define mechanistic differences between episodes requiring defibrillation versus those that spontaneously terminate.

BACKGROUND

VF is a leading cause of mortality; yet, episodes may also self-terminate. We hypothesized that the initial maintenance of human VF is dependent upon the formation and stability of VF rotors.

METHODS

We enrolled 26 consecutive patients (age 64 ± 10 years, n = 13 with left ventricular dysfunction) during ablation procedures for ventricular arrhythmias, using 64-electrode basket catheters in both ventricles to map VF prior to prompt defibrillation per the institutional review board-approved protocol. A total of 52 inductions were attempted, and 36 VF episodes were observed. Phase analysis was applied to identify biventricular rotors in the first 10 s or until VF terminated, whichever came first (11.4 ± 2.9 s to defibrillator charging).

RESULTS

Rotors were present in 16 of 19 patients with VF and in all patients with sustained VF. Sustained, but not self-limiting VF, was characterized by greater rotor stability: 1) rotors were present in 68 ± 17% of cycles in sustained VF versus 11 ± 18% of cycles in self-limiting VF (p < 0.001); and 2) maximum continuous rotations were greater in sustained (17 ± 11, range 7 to 48) versus self-limiting VF (1.1 ± 1.4, range 0 to 4, p < 0.001). Additionally, biventricular rotor locations in sustained VF were conserved across multiple inductions (7 of 7 patients, p = 0.025).

CONCLUSIONS

In patients with and without structural heart disease, the formation of stable rotors identifies individuals whose VF requires defibrillation from those in whom VF spontaneously self-terminates. Future work should define the mechanisms that stabilize rotors and evaluate whether rotor modulation may reduce subsequent VF risk.

Evaluation of the complexity of myocardial activation during ventricular fibrillation. An experimental study

INTRODUCTION AND OBJECTIVES

An experimental model is used to analyze the characteristics of ventricular fibrillation in situations of variable complexity, establishing relationships among the data produced by different methods for analyzing the arrhythmia.

METHODS

In 27 isolated rabbit heart preparations studied under the action of drugs (propranolol and KB-R7943) or physical procedures (stretching) that produce different degrees of change in the complexity of myocardial activation during ventricular fibrillation, use was made of spectral, morphological, and mapping techniques to process the recordings obtained with epicardial multielectrodes.

RESULTS

The complexity of ventricular fibrillation assessed by mapping techniques was related to the dominant frequency, normalized spectral energy, signal regularity index, and their corresponding coefficients of variation, as well as the area of the regions of interest identified on the basis of these parameters. In the multivariate analysis, we used as independent variables the area of the regions of interest related to the spectral energy and the coefficient of variation of the energy (complexity index=-0.005×area of the spectral energy regions -2.234×coefficient of variation of the energy+1.578; P=.0001; r=0.68).

CONCLUSIONS

The spectral and morphological indicators and, independently, those derived from the analysis of normalized energy regions of interest provide a reliable approach to the evaluation of the complexity of ventricular fibrillation as an alternative to complex mapping techniques.

Experiment-model interaction for analysis of epicardial activation during human ventricular fibrillation with global myocardial ischaemia

We describe a combined experiment-modelling framework to investigate the effects of ischaemia on the organisation of ventricular fibrillation in the human heart. In a series of experimental studies epicardial activity was recorded from 10 patients undergoing routine cardiac surgery. Ventricular fibrillation was induced by burst pacing, and recording continued during 2.5 min of global cardiac ischaemia followed by 30 s of coronary reflow. Modelling used a 2D description of human ventricular tissue. Global cardiac ischaemia was simulated by (i) decreased intracellular ATP concentration and subsequent activation of an ATP sensitive K⁺ current, (ii) elevated extracellular K⁺ concentration, and (iii) acidosis resulting in reduced magnitude of the L-type Ca²⁺ current I(Ca,L). Simulated ischaemia acted to shorten action potential duration, reduce conduction velocity, increase effective refractory period, and flatten restitution. In the model, these effects resulted in slower re-entrant activity that was qualitatively consistent with our observations in the human heart. However, the flattening of restitution also resulted in the collapse of many re-entrant waves to several stable re-entrant waves, which was different to the overall trend we observed in the experimental data. These findings highlight a potential role for other factors, such as structural or functional heterogeneity in sustaining wavebreak during human ventricular fibrillation with global myocardial ischaemia.