- Open Access
Uniqueness of local myocardial strain patterns with respect to activation time and contractility of the failing heart: a computational study
© The Author(s) 2018
- Received: 26 May 2018
- Accepted: 27 November 2018
- Published: 5 December 2018
Myocardial deformation measured by strain is used to detect electro-mechanical abnormalities in cardiac tissue. Estimation of myocardial properties from regional strain patterns when multiple pathologies are present is therefore a promising application of computer modelling. However, if different tissue properties lead to indistinguishable strain patterns (‘degeneracy’), the applicability of any such method will be limited. We investigated whether estimation of local activation time (AT) and contractility from myocardial strain patterns is theoretically possible.
For four different global cardiac pathologies local myocardial strain patterns for 1025 combinations of AT and contractility were simulated with a computational model (CircAdapt). For each strain pattern, a cohort of similar patterns was found within estimated measurement error using the sum of least-squared differences. Cohort members came from (1) the same pathology only, and (2) all four pathologies. Uncertainty was calculated as accuracy and precision of cohort members in parameter space. Connectedness within the cohorts was also studied.
We found that cohorts drawn from one pathology had parameters with adjacent values although their distribution was neither constant nor symmetrical. In comparison cohorts drawn from four pathologies had disconnected components with drastically different parameter values and accuracy and precision values up to three times higher.
Global pathology must be known when extracting AT and contractility from strain patterns, otherwise degeneracy occurs causing unacceptable uncertainty in derived parameters.
- Activation time
- Cardiac strain
- Computer simulation
- Reverse engineering
- System degeneration
In the healthy heart, deformation of the ventricular myocardium is relatively uniform. In contrast, many heart failure patients exhibit non-uniform myocardial deformation that may arise from heterogeneities in tissue properties such as activation time (AT) and contractility [1–3]. Non-invasively recorded myocardial deformation patterns may therefore be considered a source of important diagnostic information . Determination of regional myocardial AT and contractility is of particular clinical relevance for heart failure patients with an indication for cardiac resynchronization therapy (CRT). Clinical studies have shown that a patients response to CRT depends on both the amount and type of electrical activation delay at baseline , and on the amount and location of scarred tissue present [6, 7]. Furthermore, positioning of the left ventricular (LV) lead in or near scarred myocardium is associated with reduced benefit from CRT .
Currently, peak strain, peak strain rate, and onset time of shortening are used in the clinic to measure regional contractility and activation time from strain recordings. These parameters have limitations as they can be influenced by both activation time and regional changes in contractility . Computational models of the heart can directly relate regional myocardial tissue properties to local myocardial deformation. It may therefore be possible to estimate local AT and contractility from clinical strain pattern by fitting a computational model to reproduce the observed myocardial deformation.
The most common use for fitting to date is patient-specific computer modeling [9–12]. Currently, the models are used for fitting contractility or AT but not for both together. Elaborate fitting techniques are used to reduce the number of required simulations because each simulation is time consuming. For fitting contractility a global approach is used based on LV ejection fraction or invasively measured LV cavity pressures and volumes [13–15]. Data assimilation approaches that fit regional tissue contractility using unscented Kalman filters [16, 17] or gradient-based minimization procedures [18, 19] through comparing simulated regional myocardial deformation with cine MRI have also been shown to be tractable, raising the possibility of extracting regional tissue properties via a model.
AT is not generally taken into account when fitting a model to regional mechanical deformation, and is normally determined based on invasive electrical mapping. Fitting of regional AT is known to carry a degree of uncertainty arising from measurement noise. Konukoglu et al.  and Wallman et al.  have addressed this uncertainty using Bayesian approaches for estimation of electrophysiological conductivity parameters from endocardial activation maps. Statistical learning has also been used to estimate myocardial diffusivities based on the 12-lead electrocardiogram in patients with dilated cardiomyopathy .
Cardiac wall mechanics are dependent on many tissue and geometry parameters with non-linear relationships. The result may be that radically different sets of parameters can lead to indistinguishably similar local myocardial deformation patterns, given a realistic level of measurement error. The uncertainty in any extracted tissue parameters therefore risks becoming impractically large when fitting to real data. This degeneracy in model behavior with respect to model parameters is well-known throughout systems biology [23, 24], and has also been demonstrated in a cardiac context in electrophysiological models [25, 26].
Our aim is to establish whether or not the estimation of local AT and contractility based on myocardial strain patterns is theoretically possible by investigating whether degeneracy in regional mechanical deformation patterns can occur when AT and contractility are varied. We map degeneracy in AT and contractility using parameter space exploration and simulations of regional mechanical deformation using Hill based computational model of cardiac mechanics and cardiovascular dynamics . The model has previously been shown to provide simulations of myocardial strain in the failing heart with regional differences in AT and contractility [28, 29], and is a fast alternative for simulation of cardiac mechanics and cardiovascular system hemodynamics as compared to finite element-based cardiac models [27, 30, 31].
In this study, we determine parameter values that generate strain patterns similar to those arising from simulations using any other parameter values within the parameter space. Similarity is determined based on direct comparison within a realistic amount of measurement noise. Using connected component analysis, we determine whether these parameters lie in similar regions of parameter space as well as the size and shape of these regions. Our analysis is performed for four common global cardiac pathologies: global heart failure (HF) alone and with a severely hypocontractile region (HF + HYPO), a left bundle branch block activation pattern (HF + LBBB), and the combination thereof (HF + HYPO + LBBB).
General description of the computational model
Changing local AT alters the onset time of rise of C and, hence, determines local myocardial pre- and afterload. A full description of the sarcomere model has previously been published by Walmsley et al. .
Relating wall tension to fiber strain
In the CircAdapt model, cavity volume and pressure are related to tissue deformation through wall tension and curvature as published by Lumens et al. . Each myocardial wall may be broken up into patches. Within each patch, tissue parameters and mechanics are homogeneous, while parameters can differ between patches. The patches making up a wall are considered to act in series and so they share a common wall tension and curvature as described previously by Walmsley et al. . Briefly, cavity volume allows calculation of wall area and curvature. The value of the intrinsic sarcomere length Lsi is used to estimate stiffness in each patch and the unloaded area of each segment. Summing segmental unloaded areas and stiffness gives the wall unloaded area and wall stiffness. Wall area, wall stiffness, and wall unloaded area allow the calculation of the current wall tension that is common to all segments. The wall tension and curvature then allow calculation of transmural pressure through Laplaces law. Pressures determine flows and consequently the changes in cavity volumes. Furthermore, wall tension, segmental stiffness, and unloaded segment area allow calculation of the current area of each segment, and therefore the current local strain and stress in each segment. These properties are used to update the intrinsic sarcomere length in each segment.
Baseline heart failure simulations
A simulation of normal cardiovascular mechanics and hemodynamics representing the healthy adult cardiovascular system under baseline resting conditions was obtained as described previously  and served as starting point for four different pathology simulations.
Firstly, a global heart failure (HF) simulation was obtained by decreasing myocardial contractility (f,act) of the three ventricular walls to 50% of its normal value, so that LV ejection fraction was 30%. All three ventricular walls were activated simultaneously, 200 ms after the right atrium. The resulting HF simulation was then used as starting point for the remaining three patient simulations. Secondly, a failing heart with left bundle-branch block (HF + LBBB) was simulated by delaying onset time of septal and LV free wall activation by 25 ms and 75 ms, respectively, relative to the RV free wall. Thirdly, a failing heart with regional LV hypocontractility (HF + HYPO) was simulated by subdividing the LV free wall in two equally sized patches  and decreasing contractility of one patch to 20% of its normal value. Finally, a failing heart with both LBBB and hypocontractility (HF + HYPO + LBBB) was simulated by combining both pathologies as described above.
In all simulations, systemic peripheral resistance and total blood volume were adjusted to represent physiological homeostatic control, so that mean arterial pressure and cardiac output were maintained at their resting values (92 mmHg and 4.2 L/min, respectively). Heart rate was kept constant at 80 beats/min.
Local myofiber strain
To evaluate the effects of local AT and contractility on myocardial deformation, a small myocardial tissue segment, representing 0.3% of total LV wall mass, was created (Fig. 1a). The size of this segment was chosen so that changes of its tissue properties did not significantly affect global cardiac pump mechanics and hemodynamics. Therefore, the small segments mechanical boundary conditions, represented by the time course of wall tension, remained similar for all simulations within each of the four different baseline HF simulations introduced in the previous section.
Parameter space exploration
For each baseline HF simulation, we performed 1025 different simulations, each with a different combination of AT and contractility in the small segment. AT was varied from − 100 to + 100 ms in 41 steps of 5 ms. The value of AT is given relatively to the time of activation of the large segment which means that when AT is negative the activation of small segment precedes the large segment and when AT is positive the activation of the small segment is delayed with respect to the large segment. In some patients with heart failure, for example those with LBBB, different regions of cardiac tissue may begin to contract considerably earlier or later than others. The contractility was scaled from 2 to 200% of the large segments value, in 25 equidistant steps of 8.25%.
Strain pattern cohort
Two different types of cohorts were assembled. In the first case, only strain patterns generated from the same baseline pathology situation (HF, LBBB, HYPO, or LBBB + HYPO) as the parent strain pattern were included in the cohort. In the second case, strain patterns generated from the other baseline pathologies were also included in the cohort. The first case reflects the situation when the global cardiac pathology is known and the second when the global cardiac pathology is unknown.
Evaluation of uncertainty and degeneracy
The uncertainty of strain patterns in the parameter space was evaluated from the distance between the cohort members and parent strain pattern. Greater distance represents a larger uncertainty in how well a given strain pattern uniquely represents one set of parameters. The uncertainty was evaluated using both accuracy and precision, which were calculated for each parent strain pattern for AT and contractility, separately. The accuracy was calculated as the absolute difference between the parameter value of the parent strain and the average parameter value of the corresponding cohort members The precision was calculated as the standard deviation of the parameter values within the cohort.
We used connected component analysis (Matlab, The MathWorks Inc., Natick, Massachusetts, USA) to determine whether all parameters in the cohort were adjacent, meaning that all cohort members differed by not more than one parameter step from at least one another member from the cohort, i.e. ± 5 ms AT and 8% contractility, respectively.
Uncertainty within one global cardiac pathology
Connected component analysis showed that the parameters in all cohorts form one connected component (island) in the parameter space (Fig. 4: top panel). The cohorts were not symmetrically distributed around the parent parameter set, and the shape of each cohort was different. The shape of cohorts 1 and 4 (low contractility) in the parameter space was elongated along the AT axis relative to the contractility axis, which indicates higher uncertainty in AT as compared to contractility. In contrast, the shape of cohorts 2 and 5 (high contractility) was elongated along the contractility axis relative to the AT axis, which indicates higher uncertainty in contractility. The cohort number 3, which is centered on the baseline simulation with standard AT and contractility has circular shape, thus uncertainty of AT and contractility is similar near the baseline parameters (AT = 0 ms, Contractility = 100%). Cohorts from different regions of the parameter space generally have distinctively different shapes. Consequently, the accuracy and precision were different for each location in the parameter space. Among baseline HF simulations (HF, HF + HYPO, HF + LBBB and HF + LBBB + HYPO) the shapes of cohort parameter components were similar. The shapes of all cohort components in the parameter space and the corresponding strain patterns are shown for all four baseline simulations in Additional files 1, 2, 3, 4.
Mean uncertainty in parameter space within one global cardiac pathology (C: contractility, *absolute value)
HF + HYPO
HF + LBBB
HF + LBBB + HYPO
Uncertainty among all four global cardiac pathologies
Mean uncertainty in parameter space within combined four global cardiac pathologies (C: contractility, *absolute value)
HF + HYPO
HF + LBBB
HF + LBBB + HYPO
In this computational study, we found that when the global cardiac pathology is known, local myocardial strain patterns uniquely represent realistic values of activation time and contractility. Thus no theoretical limitation was found for extracting local tissue properties from local strain measurements. When the global pathology is not known, however, the system becomes degenerated and strain patterns can be similar despite large differences in the tissues AT and contractility. Consequently the uncertainty of AT and contractility becomes large which makes extraction of tissue parameters impractical.
Our study confirms that any methodology used for extraction of mechanical and electrical myocardial tissue properties from local deformation patterns requires a priori setting of model boundary conditions representing the patients global pathology, such as width and morphology of the QRS complex or LV ejection fraction. Even when the global pathology is known and taken into account, large heterogeneities in accuracy and precision exist across the parameter space. For example, late-activated regions will be more difficult to detect when the tissue is also hypocontractile.
Scalability of proposed parameter space exploration technique
Our findings would be challenging to reproduce using a geometrically detailed finite-element modelling approach due to the large number of simulations required (> 1000 per global pathology). By using the CircAdapt model, with its simplified geometry and consequently faster simulation speed, we were able to perform the 1025 different simulations of cardiac mechanics and hemodynamics per global pathology in < 10 h on a regular PC (64 bit 2.7 GHz CPU and 8 GB RAM). In this study, only two tissue properties were varied that are known to be relevant for the cardiac pathology in patients receiving CRT for heart failure treatment. The number of simulations grows exponentially with the number of studied parameters and could thus become impractical when including additional parameters. Meta-modelling numerical techniques for parameter space exploration can be applied to circumvent this issue. Multivariate regression [35, 36], data-driven reduction  and Bayesian sensitivity analysis  can all reduce the number of required simulations while preserving parameter space exploration. Efficient sampling schemes such as Latin hypercube sampling and polynomial chaos expansions can also be used to reduce the number of simulations required when exploring model behavior across a high-dimensional parameter space [39, 40].
Another problem which emerges with increased number of studied parameters is difficult to visualize all of the different characteristics of the studied parameter space. While two parameters can be presented in a 2D graph, three present a challenge, any more than three parameters would need to be approached differently, for example using dimensional stacking . Manifold learning could help reduce the dimensionality of the output space of the model, by building a reduced-dimensional basis for it . Statistical learning could then allow the definition of a surrogate model, which would help to infer the properties of the input/output map. In this study, we provided videos as Additional files 1, 2, 3, 4 to visualize parameter uncertainty and the corresponding strain pattern cohorts for each simulated pathology.
Appropriateness of assumed measurement noise
The similarity of two simulated strain patterns was quantified as the squared difference between the two patterns, averaged over all time points within a predefined systolic range. All strain values in this range were therefore assumed equally important for determining similarity. In the future, a more selective analysis of specific strain characteristics could be considered to improve accuracy and precision for strain-based parameter estimation.
The value of the similarity threshold used in this study (Di,p≤ 0.02) can be considered an absolute strain error of ± 2%. A recent study in subjects with normal cardiac function that compared global longitudinal strain measurements obtained using nine commercial echocardiographic deformation imaging devices  reported intra-observer and inter-observer variabilities of global longitudinal strain ranging from 0.9 to 1.7% in absolute strain %points. Furthermore, another study using MRI-tagging to measure LV deformation reported inter-study variability in peak systolic circumferential strain of ± 1.9% in absolute strain %-points . Our 2% similarity threshold thus represents a realistic level of strain measurement error.
An additional potential error that we did not take into account in this study is signal drifting. Typically, the R-top or onset of the QRS complex on the ECG is used as the reference point of zero strain, which approximately corresponds to mitral valve closure. From the reference point onwards the strain measurement can accumulate residual error due to signal drifting that can be up to a relative error of 10% by the time of peak systolic strain .
Model of cardiovascular system dynamics
The simulation results presented in this study should be interpreted as the behavior of a single Hill-type myofiber model (the small element) embedded in a larger model of the heart and circulation. In fact, the mechanical boundary conditions (tension) felt by the myofiber model can be provided by other models of the heart. Therefore, our results are likely to apply to other, more complex, models of the heart and circulation provided they use a similar Hilltype model of myofiber mechanics [13–15]. Various studies have shown that the product of simulations of local myocardial tissue mechanics and global cardiac pump function generated by the used model can be found in measurements of dyssynchronous heart failure [3, 28, 29, 31, 46].
Our study reiterates that AT and contractility interact to determine regional myocardial strain. AT may be poorly determined in regions where contractility is low, and regions of relatively high contractility may also be hard to distinguish from one another. Our simulations support the combination of electrical and mechanical measurements in order to reduce the uncertainty in the parameters to be estimated. Furthermore, global mechanical interactions between ventricular walls and between regions within the ventricular walls should be considered when examining regional strain curves.
This study is also a step towards reliable extraction of mechanical and electrical myocardial tissue properties from noninvasive myocardial deformation measurements using a model of cardiac mechanics. Electromechanical simulation-based mapping of local ventricular tissue properties, such as AT and contractility, would enable the clinician to characterise the patients complex heterogeneity of pathological electro-mechanical tissue substrates and to optimise CRT delivery by personalising pacing therapy based upon the individual patients underlying myocardial pathology.
In this simulation study we demonstrated that similarity of local myocardial strain patterns implies similarity of underlying local activation time and contractility, provided that the global pathological condition is known. These results suggest that we found no theoretical limitation for simulation-based mapping of local myocardial activation time and contractility from myocardial deformation measurements however, heterogeneous distribution of parameter uncertainty throughout the evaluated parameter space should be considered. In contrast, when the global cardiac pathology is not known, similarity of local myocardial strain patterns does not guarantee similarity of underlying local tissue parameters. The relation can become degenerate, with different regions of parameter space producing similar strain patterns. The result is a theoretical constraint on parameter extraction from myocardial deformation due to large uncertainty in the derived parameters.
BK designed the study and interpreted the data. JL provided model simulations. JW was a major contributor in writing the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
This study was financed in part by research program P3-0019 ARRS, Slovenian Research Agency, and by a grant within the framework of the Dr. E. Dekker program of the Dutch Heart Foundation (Grant 2015T082) to J. Lumens.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Prinzen FW, et al. Mapping of regional myocardial strain and work during ventricular pacing: experimental study using magnetic resonance imaging tagging. J Am Coll Cardiol. 1999;33(6):1735–42.View ArticleGoogle Scholar
- Risum N, et al. Identification of typical left bundle branch block contraction by strain echocardiography is additive to electrocardiography in prediction of long-term outcome after cardiac resynchronization therapy. J Am Coll Cardiol. 2015;66(6):631–41.View ArticleGoogle Scholar
- Lumens J, et al. Differentiating electromechanical from non-electrical substrates of mechanical discoordination to identify responders to cardiac resynchronization therapy. Circ Cardiovasc Imaging. 2015;8(9):e003744.View ArticleGoogle Scholar
- Smiseth OA, et al. Myocardial strain imaging: how useful is it in clinical decision making? Eur Heart J. 2016;37(15):1196.View ArticleGoogle Scholar
- Cleland JG, et al. An individual patient meta-analysis of five randomized trials assessing the effects of cardiac resynchronization therapy on morbidity and mortality in patients with symptomatic heart failure. Eur Heart J. 2013;34(46):3547–56.View ArticleGoogle Scholar
- Delgado V, et al. Relative merits of left ventricular dyssynchrony, left ventricular lead position, and myocardial scar to predict long-term survival of ischemic heart failure patients undergoing cardiac resynchronization therapy. Circulation. 2011;123(1):70–8.View ArticleGoogle Scholar
- Hummel JP, et al. Extent of myocardial viability predicts response to biventricular pacing in ischemic cardiomyopathy. Heart Rhythm. 2005;2(11):1211–7.View ArticleGoogle Scholar
- Sade LE, et al. The association of left ventricular lead position related to regional scar by speckle-tracking echocardiography with clinical outcomes in patients receiving cardiac resynchronization therapy. J Am Soc Echocardiogr. 2014;27(6):648–56.View ArticleGoogle Scholar
- Trayanova NA. Whole-heart modeling: applications to cardiac electrophysiology and electromechanics. Circ Res. 2011;108(1):113–28.View ArticleGoogle Scholar
- Crozier A, et al. The relative role of patient physiology and device optimisation in cardiac resynchronisation therapy: a computational modelling study. J Mol Cell Cardiol. 2016;96:93–100.View ArticleGoogle Scholar
- Sermesant M, et al. Patient-specific electromechanical models of the heart for the prediction of pacing acute effects in CRT: a preliminary clinical validation. Med Image Anal. 2012;16(1):201–15.View ArticleGoogle Scholar
- Zhao S, et al. Robust segmentation of intima-media borders with different morphologies and dynamics during the cardiac cycle. IEEE J Biomed Health Inform. 2018;22(5):1571–82.View ArticleGoogle Scholar
- Aguado-Sierra J, et al. Patient-specific modeling of dyssynchronous heart failure: a case study. Prog Biophys Mol Biol. 2011;107(1):147–55.View ArticleGoogle Scholar
- Krishnamurthy A, et al. Patient-specific models of cardiac biomechanics. J Comput Phys. 2013;244:4–21.View ArticleGoogle Scholar
- Niederer SA, et al. Length-dependent tension in the failing heart and the efficacy of cardiac resynchronization therapy. Cardiovasc Res. 2011;89(2):336–43.View ArticleGoogle Scholar
- Chabiniok R, et al. Estimation of tissue contractility from cardiac cine-MRI using a biomechanical heart model. Biomech Model Mechanobiol. 2012;11(5):609–30.View ArticleGoogle Scholar
- Marchesseau S, et al. Fast parameter calibration of a cardiac electromechanical model from medical images based on the unscented transform. Biomech Model Mechanobiol. 2013;12(4):815–31.View ArticleGoogle Scholar
- Delingette H, et al. Personalization of cardiac motion and contractility from images using variational data assimilation. IEEE Trans Biomed Eng. 2012;59(1):20–4.View ArticleGoogle Scholar
- Sermesant M, et al. Cardiac function estimation from MRI using a heart model and data assimilation: advances and difficulties. Med Image Anal. 2006;10(4):642–56.View ArticleGoogle Scholar
- Konukoglu E, et al. Efficient probabilistic model personalization integrating uncertainty on data and parameters: application to eikonal-diffusion models in cardiac electrophysiology. Prog Biophys Mol Biol. 2011;107(1):134–46.View ArticleGoogle Scholar
- Wallman M, Smith NP, Rodriguez B. Computational methods to reduce uncertainty in the estimation of cardiac conduction properties from electroanatomical recordings. Med Image Anal. 2014;18(1):228–40.View ArticleGoogle Scholar
- Zettinig O, et al. Data-driven estimation of cardiac electrical diffusivity from 12-lead ECG signals. Med Image Anal. 2014;18(8):1361–76.View ArticleGoogle Scholar
- Edelman GM, Gally JA. Degeneracy and complexity in biological systems. Proc Natl Acad Sci USA. 2001;98(24):13763–8.View ArticleGoogle Scholar
- Gutenkunst RN, et al. Universally sloppy parameter sensitivities in systems biology models. PLoS Comput Biol. 2007;3(10):1871–8.MathSciNetView ArticleGoogle Scholar
- Zaniboni M, et al. How different two almost identical action potentials can be: a model study on cardiac repolarization. Math Biosci. 2010;228(1):56–70.MathSciNetMATHView ArticleGoogle Scholar
- Wang LW, et al. Noninvasive computational imaging of cardiac electrophysiology for 3-D infarct. IEEE Trans Biomed Eng. 2011;58(4):1033–43.View ArticleGoogle Scholar
- Arts T, et al. Adaptation to mechanical load determines shape and properties of heart and circulation: the CircAdapt model. Am J Physiol Heart Circ Physiol. 2005;288(4):H1943–54.View ArticleGoogle Scholar
- Lumens J, et al. Comparative electromechanical and hemodynamic effects of left ventricular and biventricular pacing in dyssynchronous heart failure: electrical resynchronization versus left–right ventricular interaction. J Am Coll Cardiol. 2013;62(25):2395–403.View ArticleGoogle Scholar
- Leenders GE, et al. Septal rebound stretch is a strong predictor of outcome after cardiac resynchronization therapy. J Cardiac Fail. 2012;18(5):404–12.View ArticleGoogle Scholar
- Lumens J, et al. Three-wall segment (TriSeg) model describing mechanics and hemodynamics of ventricular interaction. Ann Biomed Eng. 2009;37(11):2234–55.View ArticleGoogle Scholar
- Walmsley J, et al. Fast simulation of mechanical heterogeneity in the electrically asynchronous heart using the multipatch module. PLoS Comput Biol. 2015;11(7):e1004284.View ArticleGoogle Scholar
- Arts T, et al. Relation between left ventricular cavity pressure and volume and systolic fiber stress and strain in the wall. Biophys J. 1991;59(1):93–102.View ArticleGoogle Scholar
- de Tombe PP, ter Keurs HE. Force and velocity of sarcomere shortening in trabeculae from rat heart. Effects of temperature. Circ Res. 1990;66(5):1239–54.View ArticleGoogle Scholar
- ter Keurs HE, et al. Tension development and sarcomere length in rat cardiac trabeculae. Evidence of length-dependent activation. Circ Res. 1980;46(5):703–14.View ArticleGoogle Scholar
- Sobie EA. Parameter sensitivity analysis in electrophysiological models using multivariable regression. Biophys J. 2009;96(4):1264–74.View ArticleGoogle Scholar
- Tondel K, et al. Hierarchical multivariate regression-based sensitivity analysis reveals complex parameter interaction patterns in dynamic models. Chemometric Intell Lab Syst. 2013;120:25–41.View ArticleGoogle Scholar
- Mansi BG, Hussan J, Hunter PJ, Kamen A, Comaniciu D. Data driven reduction of a cardiac myofilament model. In: Ourselin S, Rueckert D, Smith N, eds. Functional imaging and modeling of the heart: 7th international conference, FIMH 2013, London, UK, June 20–22, 2013. Berlin, Heidelberg: Springer; 2013. p. 232–40.View ArticleGoogle Scholar
- Chang ET, Strong M, Clayton RH. Bayesian sensitivity analysis of a cardiac cell model using a gaussian process emulator. PLoS ONE. 2015;10(6):e0130252.View ArticleGoogle Scholar
- Eck VG, et al. A guide to uncertainty quantification and sensitivity analysis for cardiovascular applications. Int J Numer Method Biomed Eng. 2016;32:8.MathSciNetView ArticleGoogle Scholar
- Muszkiewicz A, et al. Variability in cardiac electrophysiology: using experimentally-calibrated populations of models to move beyond the single virtual physiological human paradigm. Prog Biophys Mol Biol. 2016;120(1–3):115–27.View ArticleGoogle Scholar
- Taylor AL, et al. Structure and visualization of high-dimensional conductance spaces. J Neurophysiol. 2006;96(2):891–905.MathSciNetView ArticleGoogle Scholar
- Eynard D, et al. Multimodal manifold analysis by simultaneous diagonalization of laplacians. IEEE Trans Pattern Anal Mach Intell. 2015;37(12):2505–17.View ArticleGoogle Scholar
- Farsalinos KE, et al. Head-to-head comparison of global longitudinal strain measurements among nine different vendors: the EACVI/ASE inter-vendor comparison study. J Am Soc Echocardiogr. 2015;28(10):1171–81.View ArticleGoogle Scholar
- Donekal S, et al. Inter-study reproducibility of cardiovascular magnetic resonance tagging. J Cardiovasc Magn Reson. 2013;15:37.View ArticleGoogle Scholar
- Chen X, Liu X, Sahn D, Kim K, Xie H, O’Donnell M. Temporal and spatial registration for cardiac strain rate imaging. Ultrasonics. 2003;2:2134–7.Google Scholar
- Zhang H, et al. A meshfree representation for cardiac medical image computing. IEEE J Transl Eng Health Med. 2018;6:1800212.Google Scholar