A New Efficient Method for Detecting Phase Singularity in Cardiac Fibrillation.
Ontology highlight
ABSTRACT: 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.
SUBMITTER: Lee YS
PROVIDER: S-EPMC5131933 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The mechanisms underlying atrial fibrillation (AF), the most common sustained cardiac rhythm disturbance, remain elusive. Atrial fibrosis plays an important role in the development of AF and rotor dynamics. Both electrical wavelength (WL) and the degree of atrial fibrosis change as AF progresses. However, their combined effect on rotor core location remains unknown. The aim of this study was to analyze the effects of WL change on rotor core location in both fibrotic and non-fibrotic atria. Three patient specific fibrosis distributions (total fibrosis content: 16.6, 22.8, and 19.2%) obtained from clinical imaging data of persistent AF patients were incorporated in a bilayer atrial computational model. Fibrotic effects were modeled as myocyte-fibroblast coupling + conductivity remodeling; structural remodeling; ionic current changes + conductivity remodeling; and combinations of these methods. To change WL, action potential duration (APD) was varied from 120 to 240ms, representing the range of clinically observed AF cycle length, by modifying the inward rectifier potassium current (IK1) conductance between 80 and 140% of the original value. Phase singularities (PSs) were computed to identify rotor core locations. Our results show that IK1 conductance variation resulted in a decrease of APD and WL across the atria. For large WL in the absence of fibrosis, PSs anchored to regions with high APD gradient at the center of the left atrium (LA) anterior wall and near the junctions of the inferior pulmonary veins (PVs) with the LA. Decreasing the WL induced more PSs, whose distribution became less clustered. With fibrosis, PS locations depended on the fibrosis distribution and the fibrosis implementation method. The proportion of PSs in fibrotic areas and along the borders varied with both WL and fibrosis modeling method: for patient one, this was 4.2-14.9% as IK1 varied for the structural remodeling representation, but 12.3-88.4% using the combination of structural remodeling with myocyte-fibroblast coupling. The degree and distribution of fibrosis and the choice of implementation technique had a larger effect on PS locations than the WL variation. Thus, distinguishing the fibrotic mechanisms present in a patient is important for interpreting clinical fibrosis maps to create personalized models.
| S-EPMC6139329 | biostudies-literature
Project description:PurposeRecent investigations failed to reproduce the positive rotor-guided ablation outcomes shown by initial studies for treating persistent atrial fibrillation (persAF). Phase singularity (PS) is an important feature for AF driver detection, but algorithms for automated PS identification differ. We aim to investigate the performance of four different techniques for automated PS detection.Methods2048-channel virtual electrogram (VEGM) and electrocardiogram signals were collected for 30 s from 10 patients undergoing persAF ablation. QRST-subtraction was performed and VEGMs were processed using sinusoidal wavelet reconstruction. The phase was obtained using Hilbert transform. PSs were detected using four algorithms: (1) 2D image processing based and neighbor-indexing algorithm; (2) 3D neighbor-indexing algorithm; (3) 2D kernel convolutional algorithm estimating topological charge; (4) topological charge estimation on 3D mesh. PS annotations were compared using the structural similarity index (SSIM) and Pearson's correlation coefficient (CORR). Optimized parameters to improve detection accuracy were found for all four algorithms using F β score and 10-fold cross-validation compared with manual annotation. Local clustering with density-based spatial clustering of applications with noise (DBSCAN) was proposed to improve algorithms 3 and 4.ResultsThe PS density maps created by each algorithm with default parameters were poorly correlated. Phase gradient threshold and search radius (or kernels) were shown to affect PS detections. The processing times for the algorithms were significantly different (p < 0.0001). The F β scores for algorithms 1, 2, 3, 3 + DBSCAN, 4 and 4 + DBSCAN were 0.547, 0.645, 0.742, 0.828, 0.656, and 0.831. Algorithm 4 + DBSCAN achieved the best classification performance with acceptable processing time (2.0 ± 0.3 s).ConclusionAF driver identification is dependent on the PS detection algorithms and their parameters, which could explain some of the inconsistencies in rotor-guided ablation outcomes in different studies. For 3D triangulated meshes, algorithm 4 + DBSCAN with optimal parameters was the best solution for real-time, automated PS detection due to accuracy and speed. Similarly, algorithm 3 + DBSCAN with optimal parameters is preferred for uniform 2D meshes. Such algorithms - and parameters - should be preferred in future clinical studies for identifying AF drivers and minimizing methodological heterogeneities. This would facilitate comparisons in rotor-guided ablation outcomes in future works.
| S-EPMC7386053 | biostudies-literature
Project description:Optical phase singularities are zeros of a scalar light field. The most systematically studied class of singular fields is vortices: beams with helical wavefronts and a linear (1D) singularity along the optical axis. Beyond these common and stable 1D topologies, we show that a broader family of zero-dimensional (point) and two-dimensional (sheet) singularities can be engineered. We realize sheet singularities by maximizing the field phase gradient at the desired positions. These sheets, owning to their precise alignment requirements, would otherwise only be observed in rare scenarios with high symmetry. Furthermore, by applying an analogous procedure to the full vectorial electric field, we can engineer paraxial transverse polarization singularity sheets. As validation, we experimentally realize phase and polarization singularity sheets with heart-shaped cross-sections using metasurfaces. Singularity engineering of the dark enables new degrees of freedom for light-matter interaction and can inspire similar field topologies beyond optics, from electron beams to acoustics.
| S-EPMC8263812 | biostudies-literature
Project description:R-wave singularity (RWS) measures the intermittence or discontinuousness of R waves. It has been broadly used in QRS (QRS complex of electrocardiogram) detection, electrocardiogram (ECG) beats classification, etc. In this article, we novelly developed RWS to the analysis of QRS morphology as the measurement of ventricular dyssynchrony and tested the hypothesis that RWS could enhance the discrimination between control and acute myocardial infarction (AMI) patients. Holter ECG recordings were obtained from the Telemetric and Holter ECG Warehouse database, among which database Normal was extracted as normal controls (n = 202) and database AMI (n = 93) as typical subjects of autonomic nervous system dysfunction and cardiac electrical dyssynchrony with high risk for cardiac arrhythmias and sudden cardiac death. Experimental results demonstrate that RWS measured by Lipschitz exponent calculated from 5-min Holter recordings was significantly less negative in early AMI and late AMI than that in Normal subjects for overall, elderly, and elderly male groups, which suggested the heterogeneous depolarization of the ventricular myocardium during AMI. Receiver operating characteristic curve analyses show that combined with heart rate variability parameters, Lipschitz exponent provides higher accuracy in distinguishing between the patients with AMI and healthy control subjects for overall, elderly, elderly male, and elderly female groups. In summary, our study demonstrates the significance of using RWS to probe the cardiac electrical dyssynchrony for AMI. Lipschitz exponent may be valuable and complementary for existing cardiac resynchronization therapy and autonomic nervous system assessment.
| S-EPMC7783454 | biostudies-literature
Project description:BackgroundAlthough 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.MethodsMother 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.Results1. 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.ConclusionMother 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.
| S-EPMC4766081 | biostudies-literature
Project description:Multiple genetic factors and their interactive effects are speculated to contribute to complex diseases. Detecting such genetic interactive effects, i.e., epistatic interactions, however, remains a significant challenge in large-scale association studies.We have developed a new method, named SNPInterForest, for identifying epistatic interactions by extending an ensemble learning technique called random forest. Random forest is a predictive method that has been proposed for use in discovering single-nucleotide polymorphisms (SNPs), which are most predictive of the disease status in association studies. However, it is less sensitive to SNPs with little marginal effect. Furthermore, it does not natively exhibit information on interaction patterns of susceptibility SNPs. We extended the random forest framework to overcome the above limitations by means of (i) modifying the construction of the random forest and (ii) implementing a procedure for extracting interaction patterns from the constructed random forest. The performance of the proposed method was evaluated by simulated data under a wide spectrum of disease models. SNPInterForest performed very well in successfully identifying pure epistatic interactions with high precision and was still more than capable of concurrently identifying multiple interactions under the existence of genetic heterogeneity. It was also performed on real GWAS data of rheumatoid arthritis from the Wellcome Trust Case Control Consortium (WTCCC), and novel potential interactions were reported.SNPInterForest, offering an efficient means to detect epistatic interactions without statistical analyses, is promising for practical use as a way to reveal the epistatic interactions involved in common complex diseases.
| S-EPMC3260223 | biostudies-literature
Project description:RationaleA 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.ObjectivesHere, 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 resultsAF 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.ConclusionM/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.
| S-EPMC7841497 | biostudies-literature
Project description:Copy number variants (CNVs) play an important role in the etiology of many diseases such as cancers and psychiatric disorders. Due to a modest marginal effect size or the rarity of the CNVs, collapsing rare CNVs together and collectively evaluating their effect serves as a key approach to evaluating the collective effect of rare CNVs on disease risk. While a plethora of powerful collapsing methods are available for sequence variants (e.g., SNPs) in association analysis, these methods cannot be directly applied to rare CNVs due to the CNV-specific challenges, i.e., the multi-faceted nature of CNV polymorphisms (e.g., CNVs vary in size, type, dosage, and details of gene disruption), and etiological heterogeneity (e.g., heterogeneous effects of duplications and deletions that occur within a locus or in different loci). Existing CNV collapsing analysis methods (a.k.a. the burden test) tend to have suboptimal performance due to the fact that these methods often ignore heterogeneity and evaluate only the marginal effects of a CNV feature. We introduce CCRET, a random effects test for collapsing rare CNVs when searching for disease associations. CCRET is applicable to variants measured on a multi-categorical scale, collectively modeling the effects of multiple CNV features, and is robust to etiological heterogeneity. Multiple confounders can be simultaneously corrected. To evaluate the performance of CCRET, we conducted extensive simulations and analyzed large-scale schizophrenia datasets. We show that CCRET has powerful and robust performance under multiple types of etiological heterogeneity, and has performance comparable to or better than existing methods when there is no heterogeneity.
| S-EPMC4592002 | biostudies-literature
Project description:Social learning has been documented in a wide diversity of animals. In free-living animals, however, it has been difficult to discern whether animals learn socially by observing other group members or asocially by acquiring a new behaviour independently. We addressed this challenge by developing network-based diffusion analysis (NBDA), which analyses the spread of traits through animal groups and takes into account that social network structure directs social learning opportunities. NBDA fits agent-based models of social and asocial learning to the observed data using maximum-likelihood estimation. The underlying learning mechanism can then be identified using model selection based on the Akaike information criterion. We tested our method with artificially created learning data that are based on a real-world co-feeding network of macaques. NBDA is better able to discriminate between social and asocial learning in comparison with diffusion curve analysis, the main method that was previously applied in this context. NBDA thus offers a new, more reliable statistical test of learning mechanisms. In addition, it can be used to address a wide range of questions related to social learning, such as identifying behavioural strategies used by animals when deciding whom to copy.
| S-EPMC2674490 | biostudies-literature
Project description:Assessing the compatibility between gene family phylogenies is a crucial and often computationally demanding step in many phylogenomic analyses. Here, we describe the Evolutionary Similarity Index (IES), a means to assess shared evolution between gene families using a weighted orthogonal distance regression model applied to sequence distances. The utilization of pairwise distance matrices circumvents comparisons between gene tree topologies, which are inherently uncertain and sensitive to evolutionary model choice, phylogenetic reconstruction artifacts, and other sources of error. Furthermore, IES enables the many-to-many pairing of multiple copies between similarly evolving gene families. This is done by selecting non-overlapping pairs of copies, one from each assessed family, and yielding the least sum of squared residuals. Analyses of simulated gene family data sets show that IES's accuracy is on par with popular tree-based methods while also less susceptible to noise introduced by sequence alignment and evolutionary model fitting. Applying IES to an empirical data set of 1,322 genes from 42 archaeal genomes identified eight major clusters of gene families with compatible evolutionary trends. The most cohesive cluster consisted of 62 genes with compatible evolutionary signal, which occur as both single-copy and multiple homologs per genome; phylogenetic analysis of concatenated alignments from this cluster produced a tree closely matching previously published species trees for Archaea. Four other clusters are mainly composed of accessory genes with limited distribution among Archaea and enriched toward specific metabolic functions. Pairwise evolutionary distances obtained from these accessory gene clusters suggest patterns of interphyla horizontal gene transfer. An IES implementation is available at https://github.com/lthiberiol/evolSimIndex.
| S-EPMC8483891 | biostudies-literature