ABSTRACT
Atrial and ventricular fibrillation (AF/VF) are characterized by the repetitive regeneration of topological defects known as phase singularities (PSs). The effect of PS interactions has not been previously studied in human AF and VF. We hypothesized that PS population size would influence the rate of PS formation and destruction in human AF and VF, due to increased inter-defect interaction. PS population statistics were studied in computational simulations (Aliev-Panfilov), human AF and human VF. The influence of inter-PS interactions was evaluated by comparison between directly modeled discrete-time Markov chain (DTMC) transition matrices of the PS population changes, and M/M/∞ birth-death transition matrices of PS dynamics, which assumes that PS formations and destructions are effectively statistically independent events. Across all systems examined, PS population changes differed from those expected with M/M/∞. In human AF and VF, the formation rates decreased slightly with PS population when modeled with the DTMC, compared with the static formation rate expected through M/M/∞, suggesting new formations were being inhibited. In human AF and VF, the destruction rates increased with PS population for both models, with the DTMC rate increase exceeding the M/M/∞ estimates, indicating that PS were being destroyed faster as the PS population grew. In human AF and VF, the change in PS formation and destruction rates as the population increased differed between the two models. This indicates that the presence of additional PS influenced the likelihood of new PS formation and destruction, consistent with the notion of self-inhibitory inter-PS interactions.
Subject(s)
Atrial Fibrillation , Ventricular Fibrillation , Humans , Heart Atria , Markov Chains , ProbabilityABSTRACT
Heart failure (HF) is one of the leading causes of death worldwide. HF is associated with substantial microstructural remodeling, which is linked to changes in left ventricular geometry and impaired cardiac function. The role of myocardial remodeling in altering the mechanics of failing hearts remains unclear. Structurally based constitutive modeling provides an approach to improve understanding of the relationship between biomechanical function and tissue organization in cardiac muscle during HF. In this study, we used cardiac magnetic resonance imaging and extended-volume confocal microscopy to quantify the remodeling of left ventricular geometry and myocardial microstructure of healthy and spontaneously hypertensive rat hearts at the ages of 12 and 24 months. Passive cardiac mechanical function was characterized using left ventricular pressure-volume compliance measurements. We have developed a, to our knowledge, new structurally based biomechanical constitutive equation built on parameters quantified directly from collagen distributions observed in confocal images of the myocardium. Three-dimensional left ventricular finite element models were constructed from subject-specific in vivo magnetic resonance imaging data. The structurally based constitutive equation was integrated into geometrically subject-specific finite element models of the hearts and used to investigate the underlying mechanisms of ventricular dysfunction during HF. Using a single pair of material parameters for all hearts, we were able to produce compliance curves that reproduced all of the experimental compliance measurements. The value of this study is not limited to reproducing the mechanical behavior of healthy and diseased hearts, but it also provides important insights into the structure-function relationship of diseased myocardium that will help pave the way toward more effective treatments for HF.
Subject(s)
Heart Failure/pathology , Models, Cardiovascular , Animals , Disease Progression , Heart Failure/complications , Heart Failure/physiopathology , Myocardium/pathology , Pressure , Rats , Ventricular Dysfunction, Left/complicationsABSTRACT
Although magnetic resonance elastography (MRE) has been used to estimate isotropic stiffness in the heart, myocardium is known to have anisotropic properties. This study investigated the determinability of global transversely isotropic material parameters using MRE and finite-element modeling (FEM). A FEM-based material parameter identification method, using a displacement-matching objective function, was evaluated in a gel phantom and simulations of a left ventricular (LV) geometry with a histology-derived fiber field. Material parameter estimation was performed in the presence of Gaussian noise. Parameter sweeps were analyzed and characteristics of the Hessian matrix at the optimal solution were used to evaluate the determinability of each constitutive parameter. Four out of five material stiffness parameters (Young's modulii E1 and E3 , shear modulus G13 and damping coefficient s), which describe a transversely isotropic linear elastic material, were well determined from the MRE displacement field using an iterative FEM inversion method. However, the remaining parameter, Poisson's ratio, was less identifiable. In conclusion, Young's modulii, shear modulii and damping can theoretically be well determined from MRE data, but Poisson's ratio is not as well determined and could be set to a reasonable value for biological tissue (close to 0.5).
Subject(s)
Elasticity Imaging Techniques , Magnetic Resonance Imaging , Algorithms , Anisotropy , Computer Simulation , Phantoms, ImagingABSTRACT
INTRODUCTION: The influence of levator-ani muscles on second-stage labor is poorly understood. The ability of these muscles to stretch without damage may affect birth outcomes, but little is known about material properties, effects of pregnancy and/or ethnicity on levator-ani stiffness. There are strong associations between muscle damage and subsequent pelvic floor disorders. This study aimed to quantify levator-ani muscle stiffness during the third trimester of pregnancy and postpartum in European and Polynesian women. Associations between stiffness, obstetric variables, and the risk of intrapartum levator-ani injury (avulsion) were investigated. MATERIAL AND METHODS: This was a prospective observational pilot study. A total of 167 (106 European and 61 Polynesian) nulliparous women were recruited antenatally; 129 returned postnatally. Participants were assessed between 36 and 38 weeks' gestation and three to five months postpartum. Assessments included pelvic floor ultrasound, elastometry testing, and validated questionnaires on pelvic floor function. Logistic regression, Student t-, Chi-square and Mann-Whitney tests were used as appropriate. RESULTS: There are significant differences between antenatal and postnatal muscle stiffness measurements (p < 0.01). Stiffness was significantly higher in the European cohort (p = 0.03). There were more avulsion injuries in European (20%) than in Polynesian (9%) women. There were no significant differences in antenatal stiffness between women with and without avulsion, but change in stiffness (antenatal to postnatal) was significantly less in the avulsion group. There were no associations between stiffness, and other obstetric variables, epidural anesthesia seemed protective (p = 0.03). CONCLUSIONS: Quantification of levator-ani muscle stiffness is feasible. Muscle stiffness is significantly different before and after birth.
Subject(s)
Anal Canal/injuries , Muscle Contraction/physiology , Obstetric Labor Complications/physiopathology , Pelvic Floor Disorders/physiopathology , Pelvic Floor/physiopathology , Europe , Female , Humans , New Zealand , Pilot Projects , Postpartum Period , Pregnancy , Prospective StudiesABSTRACT
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.
Subject(s)
Heart/physiopathology , Models, Cardiovascular , Myocardial Ischemia/physiopathology , Myocardium/pathology , Ventricular Fibrillation/physiopathology , Computer Simulation , Humans , Hyperkalemia , Hypoxia , Imaging, Three-Dimensional , Myocardial Ischemia/pathology , Ventricular Fibrillation/pathologyABSTRACT
Vascularized biological tissue has been shown to increase in stiffness with increased perfusion pressure. The interaction between blood in the vasculature and other tissue components can be modeled with a poroelastic, biphasic approach. The ability of this model to reproduce the pressure-driven stiffening behavior exhibited by some tissues depends on the choice of the mechanical constitutive relation, defined by the Helmholtz free energy density of the skeleton. We analyzed the behavior of a number of isotropic poroelastic constitutive relations by applying a swelling pressure, followed by homogeneous uniaxial or simple-shear deformation. Our results demonstrate that a strain-stiffening constitutive relation is required for a material to show pressure-driven stiffening, and that the strain-stiffening terms must be volume-dependent.
Subject(s)
Blood Circulation , Elasticity , Models, Biological , Pressure , Porosity , Stress, MechanicalABSTRACT
BACKGROUND: Predicting breast tissue motion using biomechanical models can provide navigational guidance during breast cancer treatment procedures. These models typically do not account for changes in posture between procedures. Difference in shoulder position can alter the shape of the pectoral muscles and breast. A greater understanding of the differences in the shoulder orientation between prone and supine could improve the accuracy of breast biomechanical models. METHODS: 19 landmarks were placed on the sternum, clavicle, scapula, and humerus of the shoulder girdle in prone and supine breast MRIs (N = 10). These landmarks were used in an optimization framework to fit subject-specific skeletal models and compare joint angles of the shoulder girdle between these positions. FINDINGS: The mean Euclidean distance between joint locations from the fitted skeletal model and the manually identified joint locations was 15.7 mm ± 2.7 mm. Significant differences were observed between prone and supine. Compared to supine position, the shoulder girdle in the prone position had the lateral end of the clavicle in more anterior translation (i.e., scapula more protracted) (P < 0.05), the scapula in more protraction (P < 0.01), the scapula in more upward rotation (associated with humerus elevation) (P < 0.05); and the humerus more elevated (P < 0.05) for both the left and right sides. INTERPRETATION: Shoulder girdle orientation was found to be different between prone and supine. These differences would affect the shape of multiple pectoral muscles, which would affect breast shape and the accuracy of biomechanical models.
Subject(s)
Shoulder Joint , Shoulder , Humans , Shoulder/diagnostic imaging , Shoulder/physiology , Supine Position , Shoulder Joint/diagnostic imaging , Shoulder Joint/physiology , Range of Motion, Articular/physiology , Biomechanical Phenomena , Scapula/diagnostic imaging , Scapula/physiology , Rotation , Magnetic Resonance ImagingABSTRACT
BACKGROUND: Atrial fibrillation (AF) and ventricular fibrillation (VF) episodes exhibit varying durations, with some spontaneously ending quickly while others persist. A quantitative framework to explain episode durations remains elusive. We hypothesized that observable self-terminating AF and VF episode lengths, whereby durations are known, would conform with a power law based on the ratio of system size and correlation length ([Formula: see text]. METHODS: Using data from computer simulations (2-dimensional sheet and 3-dimensional left-atrial), human ischemic VF recordings (256-electrode sock, n=12 patients), and human AF recordings (64-electrode basket-catheter, n=9 patients; 16-electrode high definition-grid catheter, n=42 patients), conformance with a power law was assessed using the Akaike information criterion, Bayesian information criterion, coefficient of determination (R2, significance=P<0.05) and maximum likelihood estimation. We analyzed fibrillatory episode durations and [Formula: see text], computed by taking the ratio between system size ([Formula: see text], chamber/simulation size) and correlation length (xi, estimated from pairwise correlation coefficients over electrode/node distance). RESULTS: In all computer models, the relationship between episode durations and [Formula: see text] was conformant with a power law (Aliev-Panfilov R2: 0.90, P<0.001; Courtemanche R2: 0.91, P<0.001; Luo-Rudy R2: 0.61, P<0.001). Observable clinical AF/VF durations were also conformant with a power law relationship (VF R2: 0.86, P<0.001; AF basket R2: 0.91, P<0.001; AF grid R2: 0.92, P<0.001). [Formula: see text] also differentiated between self-terminating and sustained episodes of AF and VF (P<0.001; all systems), as well as paroxysmal versus persistent AF (P<0.001). In comparison, other electrogram metrics showed no statistically significant differences (dominant frequency, Shannon Entropy, mean voltage, peak-peak voltage; P>0.05). CONCLUSIONS: Observable fibrillation episode durations are conformant with a power law based on system size and correlation length.
Subject(s)
Atrial Fibrillation , Ventricular Fibrillation , Humans , Ventricular Fibrillation/physiopathology , Ventricular Fibrillation/diagnosis , Atrial Fibrillation/physiopathology , Atrial Fibrillation/diagnosis , Time Factors , Male , Female , Action Potentials , Computer Simulation , Heart Rate , Models, Cardiovascular , Middle Aged , Heart Conduction System/physiopathology , Electrophysiologic Techniques, Cardiac , Aged , Bayes TheoremABSTRACT
BACKGROUND: Echocardiography is widely used to evaluate left ventricular (LV) diastolic function in patients suspected of heart failure. For patients in sinus rhythm, a combination of several echocardiographic parameters can differentiate between normal and elevated LV filling pressure with good accuracy. However, there is no established echocardiographic approach for the evaluation of LV filling pressure in patients with atrial fibrillation. The objective of the present study was to determine if a combination of several echocardiographic and clinical parameters may be used to evaluate LV filling pressure in patients with atrial fibrillation. RESULTS: In a multicentre study of 148 atrial fibrillation patients, several echocardiographic parameters were tested against invasively measured LV filling pressure as the reference method. No single parameter had sufficiently strong association with LV filling pressure to be recommended for clinical use. Based on univariate regression analysis in the present study, and evidence from existing literature, we developed a two-step algorithm for differentiation between normal and elevated LV filling pressure, defining values ≥ 15 mmHg as elevated. The parameters in the first step included the ratio between mitral early flow velocity and septal mitral annular velocity (septal E/e'), mitral E velocity, deceleration time of E, and peak tricuspid regurgitation velocity. Patients who could not be classified in the first step were tested in a second step by applying supplementary parameters, which included left atrial reservoir strain, pulmonary venous systolic/diastolic velocity ratio, and body mass index. This two-step algorithm classified patients as having either normal or elevated LV filling pressure with 75% accuracy and with 85% feasibility. Accuracy in EF ≥ 50% and EF < 50% was similar (75% and 76%). CONCLUSIONS: In patients with atrial fibrillation, no single echocardiographic parameter was sufficiently reliable to be used clinically to identify elevated LV filling pressure. An algorithm that combined several echocardiographic parameters and body mass index, however, was able to classify patients as having normal or elevated LV filling pressure with moderate accuracy and high feasibility.
ABSTRACT
Changes in cardiovascular hemodynamics are closely related to the development of aortic regurgitation (AR), a type of valvular heart disease. Metrics derived from blood flows are used to indicate AR onset and evaluate its severity. These metrics can be non-invasively obtained using four-dimensional (4D) flow magnetic resonance imaging (MRI), where accuracy is primarily dependent on spatial resolution. However, insufficient resolution often results from limitations in 4D flow MRI and complex aortic regurgitation hemodynamics. To address this, computational fluid dynamics simulations were transformed into synthetic 4D flow MRI data and used to train a variety of neural networks. These networks generated super-resolution, full-field phase images with an upsample factor of 4. Results showed decreased velocity error, high structural similarity scores, and improved learning capabilities from previous work. Further validation was performed on two sets of in vivo 4D flow MRI data and demonstrated success in de-noising flow images. This approach presents an opportunity to comprehensively analyse AR hemodynamics in a non-invasive manner.
Subject(s)
Aortic Valve Insufficiency , Deep Learning , Humans , Aortic Valve Insufficiency/diagnostic imaging , Blood Flow Velocity/physiology , Hydrodynamics , Predictive Value of Tests , Magnetic Resonance Imaging/methods , Hemodynamics , Imaging, Three-Dimensional/methodsABSTRACT
BACKGROUND: Left ventricular (LV) global longitudinal strain (GLS) has been proposed as an early imaging biomarker of cardiac mechanical dysfunction. OBJECTIVE: To assess the impact of angiotensin-converting enzyme (ACE) inhibitor treatment of hypertensive heart disease on LV GLS and mechanical function. METHODS: The spontaneously hypertensive rat (SHR) model of hypertensive heart disease ( n â=â38) was studied. A subset of SHRs received quinapril (TSHR, n â=â16) from 3âmonths (mo). Wistar Kyoto rats (WKY, n â=â13) were used as controls. Tagged cardiac MRI was performed using a 4.7 T Varian preclinical scanner. RESULTS: The SHRs had significantly lower LV ejection fraction (EF) than the WKYs at 3 mo (53.0â±â1.7% vs. 69.6â±â2.1%, P â<â0.05), 14 mo (57.0â±â2.5% vs. 74.4â±â2.9%, P â<â0.05) and 24 mo (50.1â±â2.4% vs. 67.0â±â2.0%, P â<â0.01). At 24 mo, ACE inhibitor treatment was associated with significantly greater LV EF in TSHRs compared to untreated SHRs (64.2â±â3.4% vs. 50.1â±â2.4%, P â<â0.01). Peak GLS magnitude was significantly lower in SHRs compared with WKYs at 14âmonths (7.5%â±â0.4% vs. 9.9â±â0.8%, P â<â0.05). At 24âmonths, Peak GLS magnitude was significantly lower in SHRs compared with both WKYs (6.5â±â0.4% vs. 9.7â±â1.0%, P â<â0.01) and TSHRs (6.5â±â0.4% vs. 9.6â±â0.6%, P â<â0.05). CONCLUSIONS: ACE inhibitor treatment curtails the decline in global longitudinal strain in hypertensive rats, with the treatment group exhibiting significantly greater LV EF and GLS magnitude at 24 mo compared with untreated SHRs.
Subject(s)
Heart Diseases , Hypertension , Rats , Animals , Quinapril , Rats, Inbred WKY , Global Longitudinal Strain , Hypertension/drug therapy , Angiotensin-Converting Enzyme Inhibitors/pharmacology , Angiotensin-Converting Enzyme Inhibitors/therapeutic use , Rats, Inbred SHR , Blood PressureABSTRACT
Objective. Early diagnosis of heart problems is essential for improving patient prognosis.Approach. We created a non-contact imaging system that calculates the vessel-induced deformation of the skin to estimate the carotid artery pressure displacement waveforms. We present a clinical study of the system in patients (n= 27) with no underlying condition, aortic stenosis (AS), or mitral regurgitation (MR).Main results. Displacement waveforms were compared to aortic catheter pressures in the same patients. The morphologies of the pressure and displacement waveforms were found to be similar, and pulse wave analysis metrics, such as our modified reflection indices (RI) and waveform duration proportions, showed no significant differences. Compared with the control group, AS patients displayed a greater proportion of time to peak (p= 0.026 andp= 0.047 for catheter and displacement, respectively), whereas augmentation index (AIx)was greater for the displacement waveform only (p= 0.030). The modified RI for MR (p= 0.047 andp= 0.004 for catheter and displacement, respectively) was lower than in the controls. AS and MR were also significantly different for the proportion of time to peak (p= 0.018 for the catheter measurements), RI (p= 0.045 andp= 0.002 for the catheter and displacement, respectively), and AIx (p= 0.005 for the displacement waveform).Significance. These findings demonstrate the ability of our system to provide insights into cardiac conditions and support further development as a diagnostic/telehealth-based screening tool.
Subject(s)
Aortic Valve Stenosis , Mitral Valve Insufficiency , Humans , Mitral Valve Insufficiency/diagnostic imaging , Carotid Arteries , Aortic Valve Stenosis/diagnostic imaging , Aorta , Blood PressureABSTRACT
Clinical translation of personalised computational physiology workflows and digital twins can revolutionise healthcare by providing a better understanding of an individual's physiological processes and any changes that could lead to serious health consequences. However, the lack of common infrastructure for developing these workflows and digital twins has hampered the realisation of this vision. The Auckland Bioengineering Institute's 12 LABOURS project aims to address these challenges by developing a Digital Twin Platform to enable researchers to develop and personalise computational physiology models to an individual's health data in clinical workflows. This will allow clinical trials to be more efficiently conducted to demonstrate the efficacy of these personalised clinical workflows. We present a demonstration of the platform's capabilities using publicly available data and an existing automated computational physiology workflow developed to assist clinicians with diagnosing and treating breast cancer. We also demonstrate how the platform facilitates the discovery and exploration of data and the presentation of workflow results as part of clinical reports through a web portal. Future developments will involve integrating the platform with health systems and remote-monitoring devices such as wearables and implantables to support home-based healthcare. Integrating outputs from multiple workflows that are applied to the same individual's health data will also enable the generation of their personalised digital twin.Clinical Relevance- The proposed 12 LABOURS Digital Twin Platform will enable researchers to 1) more efficiently conduct clinical trials to assess the efficacy of their computational physiology workflows and support the clinical translation of their research; 2) reuse primary and derived data from these workflows to generate novel workflows; and 3) generate personalised digital twins by integrating the outputs of different computational physiology workflows.
Subject(s)
Computational Biology , Software , Computational Biology/methods , WorkflowABSTRACT
Cardiovascular imaging studies provide a multitude of structural and functional data to better understand disease mechanisms. While pooling data across studies enables more powerful and broader applications, performing quantitative comparisons across datasets with varying acquisition or analysis methods is problematic due to inherent measurement biases specific to each protocol. We show how dynamic time warping and partial least squares regression can be applied to effectively map between left ventricular geometries derived from different imaging modalities and analysis protocols to account for such differences. To demonstrate this method, paired real-time 3D echocardiography (3DE) and cardiac magnetic resonance (CMR) sequences from 138 subjects were used to construct a mapping function between the two modalities to correct for biases in left ventricular clinical cardiac indices, as well as regional shape. Leave-one-out cross-validation revealed a significant reduction in mean bias, narrower limits of agreement, and higher intraclass correlation coefficients for all functional indices between CMR and 3DE geometries after spatiotemporal mapping. Meanwhile, average root mean squared errors between surface coordinates of 3DE and CMR geometries across the cardiac cycle decreased from 7 ± 1 to 4 ± 1 mm for the total study population. Our generalised method for mapping between time-varying cardiac geometries obtained using different acquisition and analysis protocols enables the pooling of data between modalities and the potential for smaller studies to leverage large population databases for quantitative comparisons.
Subject(s)
Echocardiography, Three-Dimensional , Humans , Echocardiography, Three-Dimensional/methods , Magnetic Resonance Imaging , Bias , Heart Ventricles/diagnostic imaging , Reproducibility of Results , Ventricular Function, Left , Stroke VolumeABSTRACT
Our study methodology is motivated from three disparate needs: one, imaging studies have existed in silo and study organs but not across organ systems; two, there are gaps in our understanding of paediatric structure and function; three, lack of representative data in New Zealand. Our research aims to address these issues in part, through the combination of magnetic resonance imaging, advanced image processing algorithms and computational modelling. Our study demonstrated the need to take an organ-system approach and scan multiple organs on the same child. We have pilot tested an imaging protocol to be minimally disruptive to the children and demonstrated state-of-the-art image processing and personalized computational models using the imaging data. Our imaging protocol spans brain, lungs, heart, muscle, bones, abdominal and vascular systems. Our initial set of results demonstrated child-specific measurements on one dataset. This work is novel and interesting as we have run multiple computational physiology workflows to generate personalized computational models. Our proposed work is the first step towards achieving the integration of imaging and modelling improving our understanding of the human body in paediatric health and disease.
ABSTRACT
Computational physiological models continue to increase in complexity, however, the task of efficiently calibrating the model to available clinical data remains a significant challenge. One part of this challenge is associated with long calibration times, which present a barrier for the routine application of model-based prediction in clinical practice. Another aspect of this challenge is the limited available data for the unique calibration of complex models. Therefore, to calibrate a patient-specific model, it may be beneficial to verify that task-specific model predictions have acceptable uncertainty, rather than requiring all parameters to be uniquely identified. We have developed a pipeline that reduces the set of fitting parameters to make them structurally identifiable and to improve the efficiency of a subsequent Markov Chain Monte Carlo (MCMC) analysis. MCMC was used to find the optimal parameter values and to determine the confidence interval of a task-specific prediction. This approach was demonstrated on numerical experiments where a lumped parameter model of the cardiovascular system was calibrated to brachial artery cuff pressure, echocardiogram volume measurements, and synthetic cerebral blood flow data that approximates what can be obtained from 4D-flow MRI data. This pipeline provides a cerebral arterial pressure prediction that may be useful for determining the risk of hemorrhagic stroke. For a set of three patients, this pipeline successfully reduced the parameter set of a cardiovascular system model from 12 parameters to 8-10 structurally identifiable parameters. This enabled a significant ( > 4 × ) efficiency improvement in determining confidence intervals on predictions of pressure compared to performing a naive MCMC analysis with the full parameter set. This demonstrates the potential that the proposed pipeline has in helping address one of the key challenges preventing clinical application of such models. Additionally, for each patient, the MCMC approach yielded a 95% confidence interval on systolic blood pressure prediction in the middle cerebral artery smaller than ±10 mmHg (±1.3 kPa). The proposed pipeline exploits available high-performance computing parallelism to allow straightforward automation for general models and arbitrary data sets, enabling automated calibration of a parameter set that is specific to the available clinical data with minimal user interaction.
ABSTRACT
BACKGROUND: Ventricular fibrillation (VF) is characterized by multiple wavelets and rotors. No equation to predict the number of rotors and wavelets observed during fibrillation has been validated in human VF. OBJECTIVE: The purpose of this study was to test the hypothesis that a single equation derived from a Markov M/M/∞ birth-death process could predict the number of rotors and wavelets occurring in human clinical VF. METHODS: Epicardial induced VF (256-electrode) recordings obtained from patients undergoing cardiac surgery were studied (12 patients; 62 epochs). Rate constants for phase singularity (PS) (which occur at the pivot points of rotors) and wavefront (WF) formation and destruction were derived by fitting distributions to PS and WF interformation and lifetimes. These rate constants were combined in an M/M/∞ governing equation to predict the number of PS and WF in VF episodes. Observed distributions were compared to those predicted by the M/M/∞ equation. RESULTS: The M/M/∞ equation accurately predicted average PS and WF number and population distribution, demonstrated in all epochs. Self-terminating episodes of VF were distinguished from VF episodes requiring termination by a trend toward slower PS destruction, slower rates of PS formation, and a slower mixing rate of the VF process, indicated by larger values of the second largest eigenvalue modulus of the M/M/∞ birth-death matrix. The longest-lasting PS (associated with rotors) had shorter interactivation time intervals compared to shorter-lasting PS lasting <150 ms (â¼1 PS rotation in human VF). CONCLUSION: The M/M/∞ equation explains the number of wavelets and rotors observed, supporting a paradigm of VF based on statistical fibrillatory dynamics.
Subject(s)
Death, Sudden, Cardiac/etiology , Ventricular Fibrillation/physiopathology , Cardiac Surgical Procedures , Epicardial Mapping , Female , Heart Conduction System/physiopathology , Humans , Male , Markov Chains , Models, CardiovascularABSTRACT
Background and Objective: Renewal theory is a statistical approach to model the formation and destruction of phase singularities (PS), which occur at the pivots of spiral waves. A common issue arising during observation of renewal processes is an inspection paradox, due to oversampling of longer events. The objective of this study was to characterise the effect of a potential inspection paradox on the perception of PS lifetimes in cardiac fibrillation. Methods: A multisystem, multi-modality study was performed, examining computational simulations (Aliev-Panfilov (APV) model, Courtmanche-Nattel model), experimentally acquired optical mapping Atrial and Ventricular Fibrillation (AF/VF) data, and clinically acquired human AF and VF. Distributions of all PS lifetimes across full epochs of AF, VF, or computational simulations, were compared with distributions formed from lifetimes of PS existing at 10,000 simulated commencement timepoints. Results: In all systems, an inspection paradox led towards oversampling of PS with longer lifetimes. In APV computational simulations there was a mean PS lifetime shift of +84.9% (95% CI, ± 0.3%) (p < 0.001 for observed vs overall), in Courtmanche-Nattel simulations of AF +692.9% (95% CI, ±57.7%) (p < 0.001), in optically mapped rat AF +374.6% (95% CI, ± 88.5%) (p = 0.052), in human AF mapped with basket catheters +129.2% (95% CI, ±4.1%) (p < 0.05), human AF-HD grid catheters 150.8% (95% CI, ± 9.0%) (p < 0.001), in optically mapped rat VF +171.3% (95% CI, ±15.6%) (p < 0.001), in human epicardial VF 153.5% (95% CI, ±15.7%) (p < 0.001). Conclusion: Visual inspection of phase movies has the potential to systematically oversample longer lasting PS, due to an inspection paradox. An inspection paradox is minimised by consideration of the overall distribution of PS lifetimes.
ABSTRACT
Segmentation of the left ventricle (LV) in echocardiography is an important task for the quantification of volume and mass in heart disease. Continuing advances in echocardiography have extended imaging capabilities into the 3D domain, subsequently overcoming the geometric assumptions associated with conventional 2D acquisitions. Nevertheless, the analysis of 3D echocardiography (3DE) poses several challenges associated with limited spatial resolution, poor contrast-to-noise ratio, complex noise characteristics, and image anisotropy. To develop automated methods for 3DE analysis, a sufficiently large, labeled dataset is typically required. However, ground truth segmentations have historically been difficult to obtain due to the high inter-observer variability associated with manual analysis. We address this lack of expert consensus by registering labels derived from higher-resolution subject-specific cardiac magnetic resonance (CMR) images, producing 536 annotated 3DE images from 143 human subjects (10 of which were excluded). This heterogeneous population consists of healthy controls and patients with cardiac disease, across a range of demographics. To demonstrate the utility of such a dataset, a state-of-the-art, self-configuring deep learning network for semantic segmentation was employed for automated 3DE analysis. Using the proposed dataset for training, the network produced measurement biases of -9 ± 16 ml, -1 ± 10 ml, -2 ± 5 %, and 5 ± 23 g, for end-diastolic volume, end-systolic volume, ejection fraction, and mass, respectively, outperforming an expert human observer in terms of accuracy as well as scan-rescan reproducibility. As part of the Cardiac Atlas Project, we present here a large, publicly available 3DE dataset with ground truth labels that leverage the higher resolution and contrast of CMR, to provide a new benchmark for automated 3DE analysis. Such an approach not only reduces the effect of observer-specific bias present in manual 3DE annotations, but also enables the development of analysis techniques which exhibit better agreement with CMR compared to conventional methods. This represents an important step for enabling more efficient and accurate diagnostic and prognostic information to be obtained from echocardiography.