Search

We previously reported the feasibility and efficacy of a simulation-guided clinical catheter ablation of atrial fibrillation (AF) in an in-silico AF model. We developed a highly efficient realistic AF model reflecting the patient endocardial voltage and local conduction and tested its clinical feasibility. We acquired > 500 endocardial bipolar electrograms during right atrial pacing at the beginning of the AF ablation procedures. Based on the clinical bipolar electrograms, we generated simulated voltage maps by applying fibrosis and local activation maps adjusted for the fiber orientation. The software's accuracy (CUVIA2.5) was retrospectively tested in 17 patients and feasibility prospectively in 10 during clinical AF ablation. Results: We found excellent correlations between the clinical and simulated voltage maps (R = 0.933, p < 0.001) and clinical and virtual local conduction (R = 0.958, p < 0.001). The proportion of virtual local fibrosis was 15.4, 22.2, and 36.9% in the paroxysmal AF, persistent AF, and post-pulmonary vein isolation (PVI) states, respectively. The reconstructed virtual bipolar electrogram exhibited a relatively good similarities of morphology to the local clinical bipolar electrogram (R = 0.60 ± 0.08, p < 0.001). Feasibility testing revealed an in situ procedural computing time from the clinical data acquisition to wave-dynamics analyses of 48.2 ± 4.9 min. All virtual analyses were successfully achieved during clinical PVI procedures. We developed a highly efficient, realistic, in situ procedural simulation model reflective of individual anatomy, fiber orientation, fibrosis, and electrophysiology that can be applied during AF ablation.

Although bipolar electrograms (Bi-egms) are commonly used for catheter mapping and ablation of cardiac arrhythmias, the accuracy and reproducibility of Bi-egms have not been evaluated. We aimed to clarify the influence of the catheter orientation (CO), catheter contact angle (CA), local conduction velocity (CV), scar size, and catheter type on the Bi-egm morphology using an in silico 3-dimensional realistic model of atrial fibrillation. We constructed a 3-dimensional, realistic, in silico left atrial model with activation wave propagation including bipolar catheter models. Bi-egms were obtained by computing the extracellular potentials from the distal and proximal electrodes. The amplitude and width were measured on virtual Bi-egms obtained under different conditions created by changing the CO according to the wave direction, catheter-atrial wall CA, local CV, size of the non-conductive area, and catheter type. Bipolar voltages were also compared between virtual and clinically acquired Bi-egms. Bi-egm amplitudes were lower for a perpendicular than parallel CO relative to the wave direction (p<0.001), lower for a 90° than 0° CA (p<0.001), and lower for a CV of 0.13m/s than 0.48m/s (p<0.001). Larger sized non-conductive areas were associated with a decreased bipolar amplitude (p<0.001) and increased bipolar width (p<0.001). Among three commercially available catheters (Orion, Pentaray, and Thermocool), those with more narrowly spaced and smaller electrodes produced higher voltages on the virtual Bi-egms (p<0.001). Multiple factors including the CO, CA, CV, and catheter design significantly influence the Bi-egm morphology. Universal voltage cut-off values may not be appropriate for bipolar voltage-guided substrate mapping., Competing Interests: The authors have declared that no competing interests exist.

The purpose of this study was to investigate whether the location of the hinge affects the incidence of hinge fracture during medial closing-wedge distal femoral varus osteotomy (DFVO). Twenty knees from 10 fresh-frozen human cadavers (mean age, 75 ± 17 years) were used to perform uniplanar medial closing-wedge DFVO with a 7-mm wedge. Each specimen was randomly assigned to either group A (supracondylar hinge) or group B (lateral condylar hinge). The incidence of hinge fracture and stability was compared between both groups after uniplanar medial closing-wedge DFVO. In group A, 8 of 10 knees had a lateral cortex fracture during closure of the osteotomy gap, and all fractured knees were unstable. Two knees with an intact lateral cortical hinge showed stability under manual valgus and varus forces. After intentional breakage of the lateral cortical hinge, both knees were found to be unstable under the same force. In group B, 2 of 10 knees had a lateral cortex fracture, and 8 knees had no fractures. All specimens were found to be stable under manual valgus and varus forces. After intentional breakage of the lateral cortical hinge in group B, 2 knees were unstable, while 8 knees remained stable. This study showed a significantly higher incidence of lateral cortical hinge fracture and instability in group A than in group B during closure of the osteotomy gap., Competing Interests: None., (Thieme Medical Publishers 333 Seventh Avenue, New York, NY 10001, USA.)

Background: The arrhythmogenic role of complex atrial morphology has not yet been clearly elucidated. We hypothesized that bumpy tissue geometry can induce action potential duration (APD) dispersion and wavebreak in atrial fibrillation (AF). Methods and Results: We simulated a 2D-bumpy atrial model by varying the degree of bumpiness, and 3D-left atrial (LA) models integrated by LA computed tomographic (CT) images taken from 14 patients with persistent AF. We also analyzed wave-dynamic parameters with bipolar electrograms during AF and compared them with LA-CT geometry in 30 patients with persistent AF. In the 2D-bumpy model, APD dispersion increased (P<0.001) and wavebreak occurred spontaneously when the surface bumpiness was greater, showing phase transition-like behavior (P<0.001). The bumpiness gradient 2D-model showed that spiral wave drifted in the direction of higher bumpiness, and phase singularity (PS) points were mostly located in areas with higher bumpiness. In the 3D-LA model, PS density was higher in the LA appendage (LAA) compared with other parts of the LA (P<0.05). In 30 persistent-AF patients, the surface bumpiness of LAA was 5.8-fold that of other LA parts (P<0.001), and exceeded critical bumpiness to induce wavebreak. Wave dynamics complexity parameters were consistently dominant in the LAA (P<0.001)., Conclusions: Bumpy tissue geometry promoted APD dispersion, wavebreak, and spiral wave drift in in-silico human atrial tissue, and corresponded to clinical electroanatomical maps.

Background: We previously reported that stable rotors are observed in in-silico human atrial fibrillation (AF) models, and are well represented by a dominant frequency (DF). In the current study, we hypothesized that the outcome of DF ablation is affected by conduction velocity (CV) conditions and examined this hypothesis using in-silico 3D-AF modeling., Methods: We integrated 3D CT images of left atrium obtained from 10 patients with persistent AF (80% male, 61.8±13.5 years old) into in-silico AF model. We compared AF maintenance durations (max 300s), spatiotemporal stabilities of DF, phase singularity (PS) number, life-span of PS, and AF termination or defragmentation rates after virtual DF ablation with 5 different CV conditions (0.2, 0.3, 0.4, 0.5, and 0.6m/s)., Results: 1. AF maintenance duration (p<0.001), spatiotemporal mean variance of DF (p<0.001), and the number of PS (p = 0.023) showed CV dependent bimodal patterns (highest at CV0.4m/s and lowest at CV0.6m/s) consistently. 2. After 10% highest DF ablation, AF defragmentation rates were the lowest at CV0.4m/s (37.8%), but highest at CV0.5 and 0.6m/s (all 100%, p<0.001). 3. In the episodes with AF termination or defragmentation followed by 10% highest DF ablation, baseline AF maintenance duration was shorter (p<0.001), spatiotemporal mean variance of DF was lower (p = 0.014), and the number of PS was lower (p = 0.004) than those with failed AF defragmentation after DF ablation., Conclusion: Virtual ablation of DF, which may indicate AF driver, was more likely to terminate or defragment AF with spatiotemporally stable DF, but not likely to do so in long-lasting and sustained AF conditions, depending on CV.

Objective: Radiofrequency catheter ablation for persistent atrial fibrillation (PeAF) still has a substantial recurrence rate. This study aims to investigate whether an AF ablation lesion set chosen using in-silico ablation (V-ABL) is clinically feasible and more effective than an empirically chosen ablation lesion set (Em-ABL) in patients with PeAF. Methods: We prospectively included 108 patients with antiarrhythmic drug-resistant PeAF (77.8% men, age 60.8 ± 9.9 years), and randomly assigned them to the V-ABL ( n = 53) and Em-ABL ( n = 55) groups. Five different in-silico ablation lesion sets [1 pulmonary vein isolation (PVI), 3 linear ablations, and 1 electrogram-guided ablation] were compared using heart-CT integrated AF modeling. We evaluated the feasibility, safety, and efficacy of V-ABL compared with that of Em-ABL. Results: The pre-procedural computing time for five different ablation strategies was 166 ± 11 min. In the Em-ABL group, the earliest terminating blinded in-silico lesion set matched with the Em-ABL lesion set in 21.8%. V-ABL was not inferior to Em-ABL in terms of procedure time ( p = 0.403), ablation time ( p = 0.510), and major complication rate ( p = 0.900). During 12.6 ± 3.8 months of follow-up, the clinical recurrence rate was 14.0% in the V-ABL group and 18.9% in the Em-ABL group ( p = 0.538). In Em-ABL group, clinical recurrence rate was significantly lower after PVI+posterior box+anterior linear ablation, which showed the most frequent termination during in-silico ablation (log-rank p = 0.027). Conclusions: V-ABL was feasible in clinical practice, not inferior to Em-ABL, and predicts the most effective ablation lesion set in patients who underwent PeAF ablation.

Background/aims: Atrial anatomy and thickness may affect the electrical wave-dynamics of atrial fibrillation (AF). We explored the relationship between left atrial (LA) wall thickness (LAWT) or LA geometry and AF wave-dynamics., Methods: We included 15 patients with persistent AF (age, 62.3 ± 11.9 years) who underwent AF catheter ablation. We measured the LAWT, LA endocardial curvature, and SD-curvature (surface bumpiness) from preprocedural computed tomography images. We compared those anatomical characteristics with electrophysiologic parameters such as dominant frequency (DF), Shannon entropy (ShEn), or complex fractionated atrial electrogram (CFAE)-cycle length (CL), calculated from intracardiac bipolar electrograms (300-500 points, 5 s), acquired during ablation procedures., Results: 1. LAWT (excluding fat) varied widely among patients, locations, and types of AF. LAWT was inversely correlated with LA volume (r = -0.565, p = 0.028) and positively correlated with SD-curvature (r = 0.272, p < 0.001). 2. LAWT was positively correlated with ShEn (r = 0.233, p < 0.001) and negatively correlated with CFAE-CL (r = -0.107, p = 0.038). 3. In the multivariate linear regression analyses for AF wave-dynamics parameters, DF (β = -0.29 [95% CI -0.44--0.14], p < 0.001), ShEn (β = 0.19 [95% CI 0.12-0.25], p < 0.001), and CFAE-CL (β = 7.49 [95% CI 0.65-14.34], p = 0.032) were independently associated with LAWT., Conclusion: Regional LAWT is associated with LA structural features, and has significant correlations with the wave-dynamics parameters associated with electrical wavebreaks or rotors in patients with persistent AF.

Background: The role of the autonomic nervous system (ANS) on atrial fibrillation (AF) is difficult to demonstrate in the intact human left atrium (LA) due to technical limitations of the current electrophysiological mapping technique. We examined the effects of the ANS on the initiation and maintenance of AF by employing a realistic in silico human left atrium (LA) model integrated with a model of ganglionated plexi (GPs)., Methods: We incorporated the morphology of the GP and parasympathetic nerves in a three-dimensional (3D) realistic LA model. For the model of ionic currents, we used a human atrial model. GPs were stimulated by increasing the IK[ACh], and sympathetic nerve stimulation was conducted through a homogeneous increase in the ICa-L. ANS-induced wave-dynamics changes were evaluated in a model that integrated a patient's LA geometry, and we repeated simulation studies using LA geometries from 10 different patients., Results: The two-dimensional model of pulmonary vein (PV) cells exhibited late phase 3 early afterdepolarization-like activity under 0.05μM acetylcholine (ACh) stimulation. In the 3D simulation model, PV tachycardia was induced, which degenerated to AF via GP (0.05μM ACh) and sympathetic (7.0×ICa-L) stimulations. Under sustained AF, local reentries were observed at the LA-PV junction. We also observed that GP stimulation reduced the complex fractionated atrial electrogram (CFAE)-cycle length (CL, p<0.01) and the life span of phase singularities (p<0.01). GP stimulation also increased the overlap area of the GP and CFAE areas (CFAE-CL≤120ms, p<0.01). When 3 patterns of virtual ablations were applied to the 3D AF models, circumferential PV isolation including the GP was the most effective in terminating AF., Conclusion: Cardiac ANS stimulations demonstrated triggered activity, automaticity, and local reentries at the LA-PV junction, as well as co-localized GP and CFAE areas in the 3D in silico GP model of the LA.

Background: The point of phase singularity (PS) is considered to represent a spiral wave core or a rotor in cardiac fibrillation. Computational efficiency is important for detection of PS in clinical electrophysiology. We developed a novel algorithm for highly efficient and robust detection of PS., Methods: In contrast to the conventional method, which calculates PS based on the line integral of the phase around a PS point equal to ±2π (the Iyer-Gray method), the proposed algorithm (the location-centric method) looks for the phase discontinuity point at which PS actually occurs. We tested the efficiency and robustness of these two methods in a two-dimensional mathematical model of atrial fibrillation (AF), with and without remodeling of ionic currents., Results: 1. There was a significant association, in terms of the Hausdorff distance (3.30 ± 0.0 mm), between the PS points measured using the Iyer-Gray and location-centric methods, with almost identical PS trajectories generated by the two methods. 2. For the condition of electrical remodeling of AF (0.3 × ICaL), the PS points calculated by the two methods were satisfactorily co-localized (with the Hausdorff distance of 1.64 ± 0.09 mm). 3. The proposed location-centric method was substantially more efficient than the Iyer-Gray method, with a 28.6-fold and 28.2-fold shorter run times for the control and remodeling scenarios, respectively., Conclusion: We propose a new location-centric method for calculating PS, which is robust and more efficient compared with the conventionally used method., Competing Interests: The authors have declared that no competing interests exist.

Although 3D-complex fractionated atrial electrogram (CFAE) mapping is useful in radiofrequency catheter ablation for persistent atrial fibrillation (AF), the directions and configuration of the bipolar electrodes may affect the electrogram. This study aimed to compare the spatial reproducibility of CFAE by changing the catheter orientations and electrode distance in an in-silico left atrium (LA). We conducted this study by importing the heart CT image of a patient with AF into a 3D-homogeneous human LA model. Electrogram morphology, CFAE-cycle lengths (CLs) were compared for 16 different orientations of a virtual bipolar conventional catheter (conv-cath: size 3.5 mm, inter-electrode distance 4.75 mm). Additionally, the spatial correlations of CFAE-CLs and the percentage of consistent sites with CFAE-CL<120 ms were analyzed. The results from the conv-cath were compared with that obtained using a mini catheter (mini-cath: size 1 mm, inter-electrode distance 2.5 mm). Depending on the catheter orientation, the electrogram morphology and CFAE-CLs varied (conv-cath: 11.5±0.7% variation, mini-cath: 7.1±1.2% variation), however the mini-cath produced less variation of CFAE-CL than conv-cath (p<0.001). There were moderate spatial correlations among CFAE-CL measured at 16 orientations (conv-cath: r=0.3055±0.2194 vs. mini-cath: 0.6074±0.0733, p<0.001). Additionally, the ratio of consistent CFAE sites was higher for mini catheter than conventional one (38.3±4.6% vs. 22.3±1.4%, p<0.05). Electrograms and CFAE distribution are affected by catheter orientation and electrode configuration in the in-silico LA model. However, there was moderate spatial consistency of CFAE areas, and narrowly spaced bipolar catheters were less influenced by catheter direction than conventional catheters., Competing Interests: The authors declare no conflicts of interest.

Background: We previously reported that stable rotors were observed in in-silico human atrial fibrillation (AF) models, and were well represented by dominant frequency (DF). We explored the spatiotemporal stability of DF sites in 3D-AF models imported from patient CT images of the left atrium (LA)., Methods: We integrated 3-D CT images of the LA obtained from ten patients with persistent AF (male 80%, 61.8 ± 13.5 years old) into an in-silico AF model. After induction, we obtained 6 seconds of AF simulation data for DF analyses in 30 second intervals (T1-T9). The LA was divided into ten sections. Spatiotemporal changes and variations in the temporal consistency of DF were evaluated at each section of the LA. The high DF area was defined as the area with the highest 10% DF., Results: 1. There was no spatial consistency in the high DF distribution at each LA section during T1-T9 except in one patient (p = 0.027). 2. Coefficients of variation for the high DF area were highly different among the ten LA sections (p < 0.001), and they were significantly higher in the four pulmonary vein (PV) areas, the LA appendage, and the peri-mitral area than in the other LA sections (p < 0.001). 3. When we conducted virtual ablation of 10%, 15%, and 20% of the highest DF areas (n = 270 cases), AF was changed to atrial tachycardia (AT) or terminated at a rate of 40%, 57%, and 76%, respectively., Conclusions: Spatiotemporal consistency of the DF area was observed in 10% of AF patients, and high DF areas were temporally variable. Virtual ablation of DF is moderately effective in AF termination and AF changing into AT.

Persistent atrial fibrillation (PeAF) in humans is characterized by shortening of action potential duration (APD) and attenuation of APD rate-adaptation. However, the quantitative influences of particular ionic current alterations on rate-dependent APD changes, and effects on patterns of reentry in atrial tissue, have not been systematically investigated. Using mathematical models of human atrial cells and tissue and performing parameter sensitivity analysis, we evaluated the quantitative contributions to action potential (AP) shortening and APD rate-adaptation of ionic current remodeling seen with PeAF. Ionic remodeling in PeAF was simulated by reducing L-type Ca2+ channel current (ICaL), increasing inward rectifier K+ current (IK1) and modulating five other ionic currents. Parameter sensitivity analysis, which quantified how each ionic current influenced APD in control and PeAF conditions, identified interesting results, including a negative effect of Na+/Ca2+ exchange on APD only in the PeAF condition. At high pacing rate (2 Hz), electrical remodeling in IK1 alone accounts for the APD reduction of PeAF, but at slow pacing rate (0.5 Hz) both electrical remodeling in ICaL alone (-70%) and IK1 alone (+100%) contribute equally to the APD reduction. Furthermore, AP rate-adaptation was affected by IKur in control and by INaCa in the PeAF condition. In a 2D tissue model, a large reduction (-70%) of ICaL becomes a dominant factor leading to a stable spiral wave in PeAF. Our study provides a quantitative and unifying understanding of the roles of ionic current remodeling in determining rate-dependent APD changes at the cellular level and spatial reentry patterns in tissue.

Background: Although rotors have been considered among the drivers of atrial fibrillation (AF), the rotor definition is inconsistent. We evaluated the nature of rotors in 2D and 3D in- silico models of persistent AF (PeAF) by analyzing phase singularity (PS), dominant frequency (DF), Shannon entropy (ShEn), and complex fractionated atrial electrogram cycle length (CFAE-CL) and their ablation., Methods: Mother rotor was spatiotemporally defined as stationary reentries with a meandering tip remaining within half the wavelength and lasting longer than 5 s. We generated 2D- and 3D-maps of the PS, DF, ShEn, and CFAE-CL during AF. The spatial correlations and ablation outcomes targeting each parameter were analyzed., Results: 1. In the 2D PeAF model, we observed a mother rotor that matched relatively well with DF (>9 Hz, 71.0%, p<0.001), ShEn (upper 2.5%, 33.2%, p<0.001), and CFAE-CL (lower 2.5%, 23.7%, p<0.001). 2. The 3D-PeAF model also showed mother rotors that had spatial correlations with DF (>5.5 Hz, 39.7%, p<0.001), ShEn (upper 8.5%, 15.1%, p <0.001), and CFAE (lower 8.5%, 8.0%, p = 0.002). 3. In both the 2D and 3D models, virtual ablation targeting the upper 5% of the DF terminated AF within 20 s, but not the ablations based on long-lasting PS, high ShEn area, or lower CFAE-CL area., Conclusion: Mother rotors were observed in both 2D and 3D human AF models. Rotor locations were well represented by DF, and their virtual ablation altered wave dynamics and terminated AF.