In vivo diffusion MRI of the human heart using a 300 mT/m gradient system.

PURPOSE: This work reports for the first time on the implementation and application of cardiac diffusion-weighted MRI on a Connectom MR scanner with a maximum gradient strength of 300 mT/m. It evaluates the benefits of the increased gradient performance for the investigation of the myocardial microstructure. METHODS: Cardiac diffusion-weighted imaging (DWI) experiments were performed on 10 healthy volunteers using a spin-echo sequence with up to second- and third-order motion compensation ( M 2 $$ {M}_2 $$ and M 3 $$ {M}_3 $$ ) and b = 100 , 450 $$ b=100,450 $$ , and 1000 s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ (twice the b max $$ {b}_{\mathrm{max}} $$ commonly used on clinical scanners). Mean diffusivity (MD), fractional anisotropy (FA), helix angle (HA), and secondary eigenvector angle (E2A) were calculated for b = [100, 450] s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ and b = [100, 1000] s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ for both M 2 $$ {M}_2 $$ and M 3 $$ {M}_3 $$ . RESULTS: The MD values with M 3 $$ {M}_3 $$ are slightly higher than with M 2 $$ {M}_2 $$ with Δ MD = 0 . 05 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 5 ) $$ \Delta \mathrm{MD}=0.05\pm 0.05\kern0.3em \left[\times 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-5\right) $$ for b max = 450 s / mm 2 $$ {b}_{\mathrm{max}}=450\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ and Δ MD = 0 . 03 ± 0 . 03 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 4 ) $$ \Delta \mathrm{MD}=0.03\pm 0.03\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-4\right) $$ for b max = 1000 s / mm 2 $$ {b}_{\mathrm{max}}=1000\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ . A reduction in MD is observed by increasing the b max $$ {b}_{\mathrm{max}} $$ from 450 to 1000 s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ ( Δ MD = 0 . 06 ± 0 . 04 [ × 1 0 - 3 mm 2 / s ] ( p = 1 . 6 e - 9 ) $$ \Delta \mathrm{MD}=0.06\pm 0.04\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1.6e-9\right) $$ for M 2 $$ {M}_2 $$ and Δ MD = 0 . 08 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 1 e - 9 ) $$ \Delta \mathrm{MD}=0.08\pm 0.05\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1e-9\right) $$ for M 3 $$ {M}_3 $$ ). The difference between FA, E2A, and HA was not significant in different schemes ( p > 0 . 05 $$ p>0.05 $$ ). CONCLUSION: This work demonstrates cardiac DWI in vivo with higher b-value and higher order of motion compensated diffusion gradient waveforms than is commonly used. Increasing the motion compensation order from M 2 $$ {M}_2 $$ to M 3 $$ {M}_3 $$ and the maximum b-value from 450 to 1000 s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ affected the MD values but FA and the angular metrics (HA and E2A) remained unchanged. Our work paves the way for cardiac DWI on the next-generation MR scanners with high-performance gradient systems.

Measuring cardiomyocyte cellular characteristics in cardiac hypertrophy using diffusion-weighted MRI.

PURPOSE: This paper presents a hierarchical modeling approach for estimating cardiomyocyte major and minor diameters and intracellular volume fraction (ICV) using diffusion-weighted MRI (DWI) data in ex vivo mouse hearts. METHODS: DWI data were acquired on two healthy controls and two hearts 3 weeks post transverse aortic constriction (TAC) using a bespoke diffusion scheme with multiple diffusion times ( Δ $$ \Delta $$ ), q-shells and diffusion encoding directions. Firstly, a bi-exponential tensor model was fitted separately at each diffusion time to disentangle the dependence on diffusion times from diffusion weightings, that is, b-values. The slow-diffusing component was attributed to the restricted diffusion inside cardiomyocytes. ICV was then extrapolated at Δ = 0 $$ \Delta =0 $$ using linear regression. Secondly, given the secondary and the tertiary diffusion eigenvalue measurements for the slow-diffusing component obtained at different diffusion times, major and minor diameters were estimated assuming a cylinder model with an elliptical cross-section (ECS). High-resolution three-dimensional synchrotron X-ray imaging (SRI) data from the same specimen was utilized to evaluate the biophysical parameters. RESULTS: Estimated parameters using DWI data were (control 1/control 2 vs. TAC 1/TAC 2): major diameter-17.4 µ $$ \mu $$ m/18.0 µ $$ \mu $$ m versus 19.2 µ $$ \mu $$ m/19.0 µ $$ \mu $$ m; minor diameter-10.2 µ $$ \mu $$ m/9.4 µ $$ \mu $$ m versus 12.8 µ $$ \mu $$ m/13.4 µ $$ \mu $$ m; and ICV-62%/62% versus 68%/47%. These findings were consistent with SRI measurements. CONCLUSION: The proposed method allowed for accurate estimation of biophysical parameters suggesting cardiomyocyte diameters as sensitive biomarkers of hypertrophy in the heart.

Quantifying the impact of shape uncertainty on predicted arrhythmias.

BACKGROUND: Personalised computer models are increasingly used to diagnose cardiac arrhythmias and tailor treatment. Patient-specific models of the left atrium are often derived from pre-procedural imaging of anatomy and fibrosis. These images contain noise that can affect simulation predictions. There are few computationally tractable methods for propagating uncertainties from images to clinical predictions. METHOD: We describe the left atrium anatomy using our Bayesian shape model that captures anatomical uncertainty in medical images and has been validated on 63 independent clinical images. This algorithm describes the left atrium anatomy using Nmodes=15 principal components, capturing 95% of the shape variance and calculated from 70 clinical cardiac magnetic resonance (CMR) images. Latent variables encode shape uncertainty: we evaluate their posterior distribution for each new anatomy. We assume a normally distributed prior. We use the unscented transform to sample from the posterior shape distribution. For each sample, we assign the local material properties of the tissue using the projection of late gadolinium enhancement CMR (LGE-CMR) onto the anatomy to estimate local fibrosis. To test which activation patterns an atrium can sustain, we perform an arrhythmia simulation for each sample. We consider 34 possible outcomes (31 macro-re-entries, functional re-entry, atrial fibrillation, and non-sustained arrhythmia). For each sample, we determine the outcome by comparing pre- and post-ablation activation patterns following a cross-field stimulus. RESULTS: We create patient-specific atrial electrophysiology models of ten patients. We validate the mean and standard deviation maps from the unscented transform with the same statistics obtained with 12,000 Monte Carlo (ground truth) samples. We found discrepancies <3% and <2% for the mean and standard deviation for fibrosis burden and activation time, respectively. For each patient case, we then compare the predicted outcome from a model built on the clinical data (deterministic approach) with the probability distribution obtained from the simulated samples. We found that the deterministic approach did not predict the most likely outcome in 80% of the cases. Finally, we estimate the influence of each source of uncertainty independently. Fixing the anatomy to the posterior mean and maintaining uncertainty in fibrosis reduced the prediction of self-terminating arrhythmias from ≃14% to ≃7%. Keeping the fibrosis fixed to the sample mean while retaining uncertainty in shape decreased the prediction of substrate-driven arrhythmias from ≃33% to ≃18% and increased the prediction of macro-re-entries from ≃54% to ≃68%. CONCLUSIONS: We presented a novel method for propagating shape uncertainty in atrial models through to uncertainty in numerical simulations. The algorithm takes advantage of the unscented transform to compute the output distribution of the outcomes. We validated the unscented transform as a viable sampling strategy to deal with anatomy uncertainty. We then showed that the prediction computed with a deterministic model does not always coincide with the most likely outcome. Finally, we found that shape uncertainty affects the predictions of macro-re-entries, while fibrosis uncertainty affects the predictions of functional re-entries.

Calibrating cardiac electrophysiology models using latent Gaussian processes on atrial manifolds.

Models of electrical excitation and recovery in the heart have become increasingly detailed, but have yet to be used routinely in the clinical setting to guide personalized intervention in patients. One of the main challenges is calibrating models from the limited measurements that can be made in a patient during a standard clinical procedure. In this work, we propose a novel framework for the probabilistic calibration of electrophysiology parameters on the left atrium of the heart using local measurements of cardiac excitability. Parameter fields are represented as Gaussian processes on manifolds and are linked to measurements via surrogate functions that map from local parameter values to measurements. The posterior distribution of parameter fields is then obtained. We show that our method can recover parameter fields used to generate localised synthetic measurements of effective refractory period. Our methodology is applicable to other measurement types collected with clinical protocols, and more generally for calibration where model parameters vary over a manifold.

Atrial conduction velocity mapping: clinical tools, algorithms and approaches for understanding the arrhythmogenic substrate.

Characterizing patient-specific atrial conduction properties is important for understanding arrhythmia drivers, for predicting potential arrhythmia pathways, and for personalising treatment approaches. One metric that characterizes the health of the myocardial substrate is atrial conduction velocity, which describes the speed and direction of propagation of the electrical wavefront through the myocardium. Atrial conduction velocity mapping algorithms are under continuous development in research laboratories and in industry. In this review article, we give a broad overview of different categories of currently published methods for calculating CV, and give insight into their different advantages and disadvantages overall. We classify techniques into local, global, and inverse methods, and discuss these techniques with respect to their faithfulness to the biophysics, incorporation of uncertainty quantification, and their ability to take account of the atrial manifold.

Corrigendum: Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators.

[This corrects the article DOI: 10.3389/fphys.2021.693015.].

Bayesian Calibration of Electrophysiology Models Using Restitution Curve Emulators.

Calibration of cardiac electrophysiology models is a fundamental aspect of model personalization for predicting the outcomes of cardiac therapies, simulation testing of device performance for a range of phenotypes, and for fundamental research into cardiac function. Restitution curves provide information on tissue function and can be measured using clinically feasible measurement protocols. We introduce novel "restitution curve emulators" as probabilistic models for performing model exploration, sensitivity analysis, and Bayesian calibration to noisy data. These emulators are built by decomposing restitution curves using principal component analysis and modeling the resulting coordinates with respect to model parameters using Gaussian processes. Restitution curve emulators can be used to study parameter identifiability via sensitivity analysis of restitution curve components and rapid inference of the posterior distribution of model parameters given noisy measurements. Posterior uncertainty about parameters is critical for making predictions from calibrated models, since many parameter settings can be consistent with measured data and yet produce very different model behaviors under conditions not effectively probed by the measurement protocols. Restitution curve emulators are therefore promising probabilistic tools for calibrating electrophysiology models.

Gaussian process manifold interpolation for probabilistic atrial activation maps and uncertain conduction velocity.

In patients with atrial fibrillation, local activation time (LAT) maps are routinely used for characterizing patient pathophysiology. The gradient of LAT maps can be used to calculate conduction velocity (CV), which directly relates to material conductivity and may provide an important measure of atrial substrate properties. Including uncertainty in CV calculations would help with interpreting the reliability of these measurements. Here, we build upon a recent insight into reduced-rank Gaussian processes (GPs) to perform probabilistic interpolation of uncertain LAT directly on human atrial manifolds. Our Gaussian process manifold interpolation (GPMI) method accounts for the topology of the atrium, and allows for calculation of statistics for predicted CV. We demonstrate our method on two clinical cases, and perform validation against a simulated ground truth. CV uncertainty depends on data density, wave propagation direction and CV magnitude. GPMI is suitable for probabilistic interpolation of other uncertain quantities on non-Euclidean manifolds. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Sensitivity and Uncertainty Analysis of Two Human Atrial Cardiac Cell Models Using Gaussian Process Emulators.

Biophysically detailed cardiac cell models reconstruct the action potential and calcium dynamics of cardiac myocytes. They aim to capture the biophysics of current flow through ion channels, pumps, and exchangers in the cell membrane, and are highly detailed. However, the relationship between model parameters and model outputs is difficult to establish because the models are both complex and non-linear. The consequences of uncertainty and variability in model parameters are therefore difficult to determine without undertaking large numbers of model evaluations. The aim of the present study was to demonstrate how sensitivity and uncertainty analysis using Gaussian process emulators can be used for a systematic and quantitive analysis of biophysically detailed cardiac cell models. We selected the Courtemanche and Maleckar models of the human atrial action potential for analysis because these models describe a similar set of currents, with different formulations. In our approach Gaussian processes emulate the main features of the action potential and calcium transient. The emulators were trained with a set of design data comprising samples from parameter space and corresponding model outputs, initially obtained from 300 model evaluations. Variance based sensitivity indices were calculated using the emulators, and first order and total effect indices were calculated for each combination of parameter and output. The differences between the first order and total effect indices indicated that the effect of interactions between parameters was small. A second set of emulators were then trained using a new set of design data with a subset of the model parameters with a sensitivity index of more than 0.1 (10%). This second stage analysis enabled comparison of mechanisms in the two models. The second stage sensitivity indices enabled the relationship between the L-type Ca 2+ current and the action potential plateau to be quantified in each model. Our quantitative analysis predicted that changes in maximum conductance of the ultra-rapid K + channel I Kur would have opposite effects on action potential duration in the two models, and this prediction was confirmed by additional simulations. This study has demonstrated that Gaussian process emulators are an effective tool for sensitivity and uncertainty analysis of biophysically detailed cardiac cell models.

Quantifying atrial anatomy uncertainty from clinical data and its impact on electro-physiology simulation predictions.

Patient-specific computational models of structure and function are increasingly being used to diagnose disease and predict how a patient will respond to therapy. Models of anatomy are often derived after segmentation of clinical images or from mapping systems which are affected by image artefacts, resolution and contrast. Quantifying the impact of uncertain anatomy on model predictions is important, as models are increasingly used in clinical practice where decisions need to be made regardless of image quality. We use a Bayesian probabilistic approach to estimate the anatomy and to quantify the uncertainty about the shape of the left atrium derived from Cardiac Magnetic Resonance images. We show that we can quantify uncertain shape, encode uncertainty about the left atrial shape due to imaging artefacts, and quantify the effect of uncertain shape on simulations of left atrial activation times.

Probabilistic Interpolation of Uncertain Local Activation Times on Human Atrial Manifolds.

OBJECTIVE: Local activation time (LAT) mapping of the atria is important for targeted treatment of atrial arrhythmias, but current methods do not interpolate on the atrial manifold and neglect uncertainties associated with LAT observations. In this paper, we describe novel methods to, first, quantify uncertainties in LAT arising from bipolar electrogram analysis and assignment of electrode recordings to the anatomical mesh, second, interpolate uncertain LAT measurements directly on left atrial manifolds to obtain complete probabilistic activation maps, and finally, interpolate LAT jointly across both the manifold and different S1-S2 pacing protocols. METHODS: A modified center of mass approach was used to process bipolar electrograms, yielding a LAT estimate and error distribution from the electrogram morphology. An error distribution for assigning measurements to the anatomical mesh was estimated. Probabilistic LAT maps were produced by interpolating on a left atrial manifold using Gaussian Markov random fields, taking into account observation errors and characterizing LAT predictions by their mean and standard deviation. This approach was extended to interpolate across S1-S2 pacing protocols. RESULTS: We evaluated our approach using recordings from three patients undergoing atrial ablation. Cross-validation showed consistent and accurate prediction of LAT observations both at different locations on the left atrium and for different S1-S2 intervals. SIGNIFICANCE: Interpolation of scalar and vector fields across anatomical structures from point measurements is a challenging problem in biomedical engineering, compounded by uncertainties in measurements and meshes. New methods and approaches are required, and in this paper, we have demonstrated an effective method for probabilistic interpolation of uncertain LAT.

Fitting two human atrial cell models to experimental data using Bayesian history matching.

Cardiac cell models are potentially valuable tools for applications such as quantitative safety pharmacology, but have many parameters. Action potentials in real cardiac cells also vary from beat to beat, and from one cell to another. Calibrating cardiac cell models to experimental observations is difficult, because the parameter space is large and high-dimensional. In this study we have demonstrated the use of history matching to calibrate the maximum conductance of ion channels and exchangers in two detailed models of the human atrial action potential against measurements of action potential biomarkers. History matching is an approach developed in other modelling communities, based on constructing fast-running Gaussian process emulators of the model. Emulators were constructed from a small number of model runs (around 102), and then run many times (>106) at low computational cost, each time with a different set of model parameters. Emulator outputs were compared with experimental biomarkers using an implausibility measure, which took into account experimental variance as well as emulator variance. By repeating this process, the region of non-implausible parameter space was iteratively reduced. Both cardiac cell models were successfully calibrated to experimental datasets, resulting in sets of parameters that could be sampled to produce variable action potentials. However, model parameters did not occupy a small range of values. Instead, the history matching process exposed inputs that can co-vary across a wide range and still be consistent with a particular biomarker. We also found correlations between some biomarkers, indicating a need for better descriptors of action potential shape.

Quantifying the effect of uncertainty in input parameters in a simplified bidomain model of partial thickness ischaemia.

Reduced blood flow in the coronary arteries can lead to damaged heart tissue (myocardial ischaemia). Although one method for detecting myocardial ischaemia involves changes in the ST segment of the electrocardiogram, the relationship between these changes and subendocardial ischaemia is not fully understood. In this study, we modelled ST-segment epicardial potentials in a slab model of cardiac ventricular tissue, with a central ischaemic region, using the bidomain model, which considers conduction longitudinal, transverse and normal to the cardiac fibres. We systematically quantified the effect of uncertainty on the input parameters, fibre rotation angle, ischaemic depth, blood conductivity and six bidomain conductivities, on outputs that characterise the epicardial potential distribution. We found that three typical types of epicardial potential distributions (one minimum over the central ischaemic region, a tripole of minima, and two minima flanking a central maximum) could all occur for a wide range of ischaemic depths. In addition, the positions of the minima were affected by both the fibre rotation angle and the ischaemic depth, but not by changes in the conductivity values. We also showed that the magnitude of ST depression is affected only by changes in the longitudinal and normal conductivities, but not by the transverse conductivities.

Pattern formation in polymer blend thin films: surface roughening couples to phase separation.

We introduce a model for thin films of multicomponent fluids that includes lateral and vertical phase separation, preferential component attraction at both surfaces, and surface roughening. We apply our model to thin films of binary polymer blends, and use simulations of different surface-blend interaction regimes to investigate pattern formation. We demonstrate that surface roughening couples to phase separation. For films undergoing lateral phase separation via a transient wetting layer, this results in distinct stages of roughening as the film evolves between different phase equilibria.

Lateral phase separation in polymer-blend thin films: surface bifurcation.

We use simulations of a binary polymer blend confined between selectively attracting walls to identify and explain the mechanism of lateral phase separation via a transient wetting layer. We first show that equilibrium phases in the film are described by one-dimensional phase equilibria in the vertical (depth) dimension, and demonstrate that effective boundary conditions imposed by the film walls pin the film profile at the walls. We then show that, prior to lateral phase separation, distortion of the interface in a transient wetting layer is coupled to lateral phase separation at the walls. Using Hamiltonian phase portraits, we explain a "surface bifurcation mechanism" whereby the volume fraction at the walls evolves and controls the dynamics of the phase separation. We suggest how solvent evaporation may assist our mechanism.

Breakup of a transient wetting layer in polymer blend thin films: unification with 1D phase equilibria.

We show that lateral phase separation in polymer blend thin films can proceed via the formation of a transient wetting layer which breaks up to give a laterally segregated film. We show that the growth of lateral inhomogeneities at the walls in turn causes the distortion of the interface in the transient wetting layer. By addressing the 1D phase equilibria of a polymer blend thin film confined between selectively attracting walls, we show that the breakup of a transient wetting layer is due to wall-blend interactions; there are multiple values of the volume fraction at the walls which solve equilibrium boundary conditions. This mechanism of lateral phase separation should be general.

Surface roughening in polymer blend thin films by lateral phase separation: a thermodynamic mechanism.

We use a phase portrait approach to study the phase equilibria of a symmetric binary polymer blend confined between an attracting wall and a neutral wall. We find multiple solutions for this wall regime where only one solution exists for antisymmetric walls. We also argue that, when one wall is a free surface, roughening of the free surface upon lateral phase separation is expected, since thermodynamics alone implies that the coexisting phases should exist at different heights to minimise energy. We focus on the graphical insight provided by this phase portrait approach.