The ES pressure values for LBBB patients measured through aortic valves.
Detection of regional ventricular dysfunction is a challenging problem. This study presents an efficient method based on ultrasound (US) imaging and finite element (FE) analysis, for detecting akinetic and dyskinetic regions in the left ventricle (LV). The underlying hypothesis is that the contraction of a healthy LV is approximately homogeneous. Therefore, any deviations between the image-based measured deformation and a homogeneous contraction FE model should correspond to a pathological region. The method was first successfully applied to synthetic data simulating an acute ischemia; it demonstrated that the pathological areas were revealed with a higher contrast than those observed directly in the deformation maps. The technique was then applied to a cohort of eight left bundle branch block (LBBB) patients. For this group, the heterogeneities were significantly less pronounced than those revealed for the synthetic cases but the method was still able to identify the abnormal regions of the LV. This study indicated the potential clinical utility of the method by its simplicity in a patient-specific context and its ability to quickly identify various heterogeneities in LV function. Further studies are required to determine the model accuracy in other pathologies and to investigate its robustness to noise and image artifacts.
- acute ischemia
- left ventricle
- cardiac mechanics
- GE Healthcare US
Realistic finite element (FE) simulation of regional defects caused by pathologies such as acute ischemia remains a significant problem for physiologically based FE analysis of the heart. From the moment that the cell death starts to the thinning and the expansion of the scar surface of myocardium tissue, it takes approximately 21 days before the remodeling phase starts . Therefore, an accessible method to detect and characterize the influenced region can provide a better plan for diagnosis of these patients with medications and myocardial injection to avoid tissue progressive remodeling [2, 3].
Tracking the deformation of the LV has been actively pursued for the last decades, but important information is lacking about in-vivo wall stress reconstruction in medical imaging modalities for computational purpose . Accurate FE models based on high-resolution medical images can be applied for stress reconstruction, but these models are still limited to medical research and not generally available in the clinic.
Generally, the primary tool for ischemia detection is the analysis of standard electrocardiogram (ECG) signals. The majority of research focus regarding mathematical and computational tools of ischemic zone detection were addressed by electro-physiology models [5, 6, 7, 8]. An alternative method for ischemia detection is to assess cardiac function with magnetic resonance imaging (MRI) modalities. However, the abnormal zone detection in MRI modalities is limited and remains a challenging problem [9, 10, 11]. In this context, the speckle tracking technique which employs the non-invasive ultrasound (US)  has become a commonplace technique in clinical daily practice  due to the accessibility of this system.
In this chapter, we study the possibilities of computing a region of abnormal function from 4D strain acquired with the EchoPac® system. The strain estimates are based on speckle tracking, and capture 3D strain through time for 17 LV regions based on the standardized segments of the American Heart Association (AHA) .
Coupling US data with computations based on the FE method may provide more detailed information about the mechanical state of the patient’s heart. The purpose of the present chapter is to apply FE simulations to detect regions of abnormal function in the LV, based on the hypothesis that a homogeneous contraction is a regional assumption for a healthy LV, and large deviations between measured deformations and FE simulations based on homogeneous contraction will correspond to an abnormal region. We aim to investigate the validity of this hypothesis using the 4D strain data acquired with EchoPac®.
The chapter is organized as follows. First, we present an overview of our pipeline and of the equations for the systolic phase of the heart cycle. Next, we test the accuracy of the method using synthetic data, simulating zones of acute ischemia in the LV. Finally, the methodology is applied to a cohort of patients with left bundle branch block (LBBB), and we discuss the results and the strengths and limitations of the proposed method.
2. Material and method
Akinetic or dyskinetic regions such as ischemic regions or infarcts may be detected in medical images as regions of abnormal deformation. However, healthy tissue surrounding a pathological region may also show abnormal patterns of deformation, which makes it challenging to locate the pathological region. The purpose of the present paper is to investigate an alternative method for detecting abnormal regions, by comparing measured strains with strains resulting from a homogeneous contraction. As noted earlier, the underlying hypothesis is that the healthy contraction is approximately homogeneous [15, 16], and the largest deviation will therefore be associated with the pathological region. To test our hypothesis, we propose to:
Build a patient-specific mesh from the images (using an in-house mesh morphing method ),
Apply a uniform contraction and tune this to give the closest match to the patient’s strain data,
Evaluate the spatial deviations, and compare the region of largest deviation to the (known) regions of abnormalities.
These main steps of the methodology are illustrated in Figure 1. When applying the algorithm to the patient data, it is necessary to generate a patient-specific FE mesh for each individual case. Besides, the accuracy of the geometry associated with the generic LV model is greater than the patient-specific geometries, hence, the need to adapt the generic model to each patient. This task is accomplished by a previously developed mesh morphing method, which takes a generic LV mesh and adapts it to the patients’ image data . The rest of the methodology consists of a minimization loop on the active tension
The main reason why we focus on the systolic phase is that in acute ischemia, the passive behavior of the abnormal zone during diastolic phase stays intact, whereas the first appearance of the defected region is during the active contraction . As a result, remodeling of abnormal tissue induces a different deformation due to the average stress in regional LV motion.
2.1. Reference FE simulation model for systolic phase
To examine the abnormal tissue size effects on global and regional deformation, we need to explain the model which was used for FE systolic phase. We employed a healthy volunteer LV extracted from 4D US  at end diastole (ED). A 3D mesh with eight-noded linear hexahedral elements was created using Abaqus® mesh generation software. The reference mesh had 155,172 nodes and 141,405 elements.
To incorporate the known material anisotropy of the LV during the systolic phase, we defined a local curvilinear coordinate system aligned with the local fibre orientation, to model myocyte contraction [19, 20]. We assumed a helical fibre distribution, with a helical angle varying from −70° at the epicardial surface to 60° at the endocardial surface (0° is the circumferential direction). At each point, we defined the local orthonormal basis as (), where () is aligned with the local fibre .
We were interested in computing the mechanical state of the tissue at ES. The tissue deformation in the systolic phase results from a combination of active and passive mechanical properties. We modeled this interaction using a so-called active stress approach, where the total tissue stress is decomposed into an active () and a passive () part.
Following our assumption of homogeneous contraction, the value of
In the literature, different strain energy functions have been proposed for the cardiac tissue [24, 25, 26]. To reduce the complexity and number of parameters in the model, and based on the assumption that at ES the tissue stress is dominated by the active component, we used a simple Mooney-Rivlin strain energy model, as previously used by Marchesseau et al. :
We performed a reference FE simulation using the following MPs:
We used these boundary conditions, MPs and contraction model for the
2.2. Application on synthetic acute ischemia
Two reference meshes have been used for application of our developed pipeline.
2.2.1. Ellipsoid reference mesh
A truncated ellipsoid mesh has been generated using Abaqus® software. The height of the mesh was 9 cm, and the wall thickness was 1 cm. The fibre directions were assigned as outlined earlier. We defined a small zone at the equatorial level of the septal wall (see Figure 2) as acutely ischemic, i.e. non-contractile. The active stress in this region was set to zero, while the rest of the reference FE model,
To provide a first test of the hypothesis underlying this research, we performed an FE simulation with an AAS, and tuned the active stress parameter to match the strains of the synthetic reference model as closely as possible. The cost function was defined as the sum of the squared differences (LSQ) between the AAS strains and the reference strains, and a simple minimization was performed by gradually increasing
2.2.2. Healthy reference mesh
The healthy mesh obtained from US imaging in Section 2.1 was used in this study case. The abnormal zone for this mesh was defined on the mid-cavity anterior wall (Figure 2). The boundary condition, systolic blood pressure, contraction model and the MPs were attributed as in the ellipsoid study case. The midwall regions were also generated in order to follow the regional deformations through the systolic phase. Reorientation of the fibre directions in abnormal LV tissue is an open question [31, 32, 33, 34]. Therefore, we leave the fibres unchanged in the abnormal zone as defined for a healthy LV case in Section 2.1.
2.3. Prediction of abnormal tissue zone for LBBB patients
A cohort of patients suffering from LBBB was available for studying our proposed method for abnormal zone detection (with chronic infarcts). Left bundle branch block (LBBB) is a cardiac conduction abnormality which causes the LV to contract later than the right ventricle. We presented, previously, our FE method for a healthy LV contraction. In order to introduce the patient data to the FE software, we developed an in-house code to morph the reference FE mesh to the patient geometries. We demonstrated the steps for estimating homogeneous active contraction value data of acute ischemia in Section 2.2.
Eight LBBB patients were recruited after informed consent at Rikshospitalet, Oslo, Norway. LV triangulated mesh geometries were acquired by a standard GE Healthcare echocardiography examination . The patient’s LV cavity pressures were measured via an arterial catheter inside the left ventricle’s cavity. The 17 regional strains were measured using 4D strain tracking EchoPac® US system manufactured by GE Healthcare . EchoPac® algorithm calculates the classical longitudinal, circumferential and area strains for each segment at the midwall between endocardial and epicardial surfaces from the segment’s dimensions  (Figure 4).
From area strain (
We morphed the reference mesh (refer to Section 2.1) onto the geometry of each patient at ED time step for which US geometries were acquired. For this, endocardial and epicardial nodes of the deformable mesh were projected onto the patient’s triangulated surfaces from a defined LV
For each patient, after application of the mesh morphing method, we employed the pressure values at ES from pressure curves obtained during the medical intervention (Table 1). This cavity pressure was applied as the boundary condition on the endocardial surface as previously mentioned in Section 2.1. For the sake of simulation costs, for each patient, we considered a set of
|ES pressure (kPa)||16.8||14||14.27||15.47||12.67||16.53||14.53||8666|
In this section, we describe the results of applying the proposed method to synthetic reference data and patient data from a group of LBBB patients.
3.1. Application to synthetic data
For the ellipsoid mesh, the closest match was achieved for
Figure 6 shows the results of applying the method on the healthy LV geometry, with synthetically generated acute ischemic data. In this case, a minimum cost function value of 0.0387 was obtained for
3.2. Application on LBBB patients’ data
We applied the mesh morphing method explained in Section 2.3 to construct patient-specific FE meshes for 8 LBBBB patients, shown in the figure by Behdadfar et al. . Table 2 shows the minimum cost value and the AAS values for each patient. In addition, the measured and simulated ED and ES cavity volumes are shown, as well as the EFs for both cases.
|EDV-S (ml)||ESV-S (ml)||EF-S%||EDV-P (ml)||ESV-P (ml)||EF-P%|
The deformations obtained from US are in circumferential, longitudinal and area strain of the midwall mesh. Then, the RS was obtained by incompressibility hypothesis of the LV segments as explained by Heimdal . The radial, circumferential and longitudinal strains are shown in Figure 7 for patients 4 and 7 where the pipeline detected the abnormal zone. The FE strain results were subtracted from the patient’s respective data and are shown in this figure as well.
Figure 8 shows the transversal and longitudinal cuts of LBBB patients’ geometries superimposed by the FE simulation results for the optimal AAS. These cuts are tuned to show the maximum error zone (from 17 segments), obtained from LSQ application of regional deformation. Patient 5 had convergence issues and was not considered for further analysis.
In this study, a new approach for detecting potentially infarcted areas was first validated on synthetic data and then applied to a cohort of 8 LBBB patients before resynchronization. The results revealed zones of moderate heterogeneities which were often of smaller thickness but the heterogeneities were significantly less pronounced than that revealed for a synthetic case of acute ischemia. The method also revealed that the tissue tends to be stiffer in the lateral wall of LBBB patients. It is very interesting to notice that the obtained results confirm what was already observed by other researchers using standard medical imaging. More specifically, Veress et al.  showed the typical septum-related abnormal wall motion and impairment of wall thickening during systole, which caused LV remodeling. Often, technologies such as MRI or PET (positron emission tomography) were used to characterize the infarct heterogeneities as well as histological images [36, 37, 38]. We examined if such ability to detect abnormal tissue can stand for novel 4D strain EchoPac® system.
In order to detect the abnormal tissue, we identified an AAS for all FE meshes which was the closest match to the patient’s midwall mesh deformation (Figure 1). This AAS represents a homogeneous contraction for each patient [15, 16]. In this condition, it is expected that for a healthy and homogeneous contraction, the FE deformation responses will be different in abnormal tissue than the rest of the LV model. In addition, this value represents the average movement of remote (uninfarcted) and injured tissue.
Often, patients suffering from abnormal contractility are evaluated by an optimized gradient of activation stress map. Wenk et al.  have developed an animal-specific FE model for such reconstruction by minimizing the regional deformation error between experimental data and the FE model. In this study, a border zone was considered for early stage of tissue damaging process (due to calcium concentrations) which relates the injured to the remote tissue. They confirm that the infarct zone has no contractile action (zero active value).
4.1. Impact of MPs on the FE LV deformation
In Table 2, the AASs are mentioned for optimal cost values. These identified values are strongly coupled to the wall thickness, the blood pressure and the cavity volume. The FE results are based on the midwall strain deformation which is not representative data for the complete wall motion and might be affected by artifacts. Some meshes are inflated such as in cases 4 and 7 in Table 2, even if, these cases are the most successful results of the abnormal tissue detection. The reason for this inflation might be the compliance MPs for the studied cases. To study this hypothesis, we increased the elastic parameters (
For this test, patient 7 has been selected as the posterior and neighboring walls were deformed dramatically in FE simulation while these walls had not been moved during ED to ES movement. The cost value for this test from 0.7901 with 175 kPa of AAS was reduced to [0.1893, 0.0775, 0.0626, 0.0592, 0.0591, 0.0564], respectively, and it is shown in Figure 4 for factor 12 (
Figure 8 shows several successful zone detections wherein the maximum error happens to be at the abnormal region. In case 1, the maximum deformation was observed at the Inf-Sept-AntSept regions where this neighborhood wall slightly inflates (blue lines at Max-error in longitudinal cut) and the method detected this region successfully. However, in case 6, the inhomogeneous strain was observed at the septal wall wherein the method detected the posterior wall as the maximum error. These results show that the strain data in some cases can be difficult to be judged on the detected zone and is subjective to the visual detection of the patient’s strain and geometry data. However, a FE simulation of pathology cases permit us to narrow down the potentially defeated zone to less regions, especially, in case 4 and 7 in Figure 7 and to study the mechanical state of the cardiac tissue by analyzing the wall stress as a post-processing step.
The results showed a good agreement with several patient deformations in RS, which has been shown to be a better indicator than circumferential or longitudinal deformations due to the homogeneous distribution of these deformation values, for discrimination of possible heterogeneities in the active stress and abnormal tissue movement. We observed that the identified active stress is highly coupled to the patient geometry and wall thickness . This is the advantage of the new generations of US system in extracting the RS from measuring the surface of a given segment that is not possible for modalities such as tagging MRI .
4.2. Remodeling in LBBB patients and its impact on the pipeline results
The LBBB, itself, causes tissue remodeling due to the redistribution of LV workload in the long run, especially, in circumferential shortening, cardiac mass, septal hypoperfusion and myocardial blood flow [44, 45, 46, 47]. However, this pathology is, often, the result of other cardiac diseases. Statistically, 40% of patients with cardiomyopathy (CMP), weakening of the heart muscle and congestive heart failure have unsynchronized ventricular contraction. The resynchronization procedure helps in reverse remodeling and hospitalization-free survival rate of 70–90% in these patients [48, 49, 50].
The observed regional differences did not exist in the case of ellipsoid and healthy meshes illustrated previously, where the wall thickness is nearly homogeneous. Thus, the change in MPs or wall thickness in LBBB patients with previous infarct history should reduce the efficiency of this method for some patients with high heterogeneity in the wall deformation due to tissue remodeling. In Figure 8 and Table 2, we can observe that the FE simulation responses in the maximal discrepancies with the patient’s data are mostly over or under estimates of the tissue deformation. These maximal discrepancies have been observed to be at the high heterogeneities in wall thickness.
Patient 5 had convergence issues due to large cavity volume (319 ml at ED) and thin wall thickness (remodeled tissue). The large cavity volume is related to an aneurysmal bulging which made contraction simulation fail to reach a solution.
Several studies showed that LBBB patients respond poorly (30–50%) to resynchronization procedure for multiple reasons and cases with developed scarred tissue tend to not respond at all [51, 52, 53]. This study is a tool for patient cases which are not yet affected by this remodeling phase and are in the early stages of LBBB pathology. Information on areas with pathology within the left ventricular wall may be helpful in planning therapy with a resynchronization pacemaker .
One major limitation in this study is the choice of material and contraction models to simplify the complexity and computational costs. In addition, we defined an
Another limitation is the lack of patient’s history information. In this study, the patients were known to have LBBB pathology with infarct history but there was no other information on co-morbidities that may cause depressed myocardial contraction. The large cavity volume, low EF and low wall thickness can be related to various pathologies such as idiopathic dilated CMP, and ischemic cardiomyopathy. So in this work, it is not possible to categorize the patients with their pathologies (such as patient 5) and the identified AAS. We considered a stress-free configuration to simulate the systolic phase for our cohort so this issue should be further addressed.
In this study, a new approach for detecting potentially infarcted areas was introduced, validated and then applied onto a cohort of eight LBBB patients before their cardiac resynchronization therapy. The results revealed zones of moderate heterogeneities which were often of regions of thinner wall. However, the heterogeneities were significantly less pronounced than what was revealed for a synthetic case of acute ischemia, which we interpreted as an effect of remodeling induced by the therapy. The method also revealed that the tissue tends to be stiffer in the lateral wall of LBBB patients, as the deformations are not as pronounced as they are in the simulation. This study is promising for the assessment of LBBB and its quantification using FE simulations. Again, this illustrates the importance of patient-specific FE simulations in the domain of cardiac biomechanics.
Further work is required to transfer the promising synthetic results to real acute ischemic patient cases. A possible improvement is to reconstruct the patient fiber orientations from MRI imaging modalities such as diffusion tensor and tagging MRI to improve the deformation trajectories during contraction . The possibility to compare this study results with cardiac MRI for identifying the abnormal tissue is also an interesting future step. Finally, it would be interesting to compare the reconstructions made with our methodology on LBBB patients at different stages of remodeling after the cardiac resynchronization treatment.
The authors gratefully acknowledge the support and data provided by the Center for Cardiological Innovation, funded by the Research Council of Norway. SA is grateful to the European Research Council for grant ERC-2014-CoG BIOLOCHANICS.
Fishbein MC, Maclean D, Maroko PR. The histopathologic evolution of myocardial infarction. Chest. 1978; 73(6):843-849
Sabbah HN, Wang M, Gupta RC, Rastogi S, Ilsar I, Sabbah MS, Kohli S, Helgerson S, Lee RJ. Augmentation of left ventricular wall thickness with alginate hydrogel implants improves left ventricular function and prevents progressive remodeling in dogs with chronic heart failure. JACC: Heart Failure. 2013; 1(3):252-258
Lee LC, Wall ST, Klepach D, Ge L, Zhang Z, Lee RJ, Hinson A, Gorman JH, Gorman RC, Guccione JM. Algisyl-lvr-tm with coronary artery bypass grafting reduces left ventricular wall stress and improves function in the failing human heart. International Journal of Cardiology. 2013; 168(3):2022-2028
Veress AI, Paul Segars W, Tsui BMW, Gullberg GT. Incorporation of a left ventricle finite element model defining infarction into the xcat imaging phantom. IEEE Transactions on Medical Imaging. 2011; 30(4):915-927
Sundnes J, Lines GT, Cai X, Nielsen BF, Mardal K-A, Tveito A. In: Computing the electrical activity in the heart. Vol. 1. Springer Science & Business Media; 2007
Nejib Z. Theoretical and numerical study of the electric activity of the heart. Modeling and Numerical Simulation of Electrocardiograms [PhD thesis]. Université Paris Sud-Paris XI; 2009
Shi WV, Chang TN, Zhou MC. Method to detect cardiac abnormalities based on electrocardiography and sinoatrial pacemaker model. In: 2010 International Conference on Mechatronics and Automation (ICMA), IEEE; 2010. pp. 566-571
Oduneye SO, Pop M, Shurrab M, Biswas L, Ramanan V, Barry J, Crystal E, Wright GA. Distribution of abnormal potentials in chronic myocardial infarction using a real time magnetic resonance guided electrophysiology system. Journal of Cardiovascular Magnetic Resonance. 2015; 17(1):1
Underwood SR, Rees RS, Savage PE, Klipstein RH, Firmin DN, Fox KM, Poole-Wilson PA, Longmore DB. Assessment of regional left ventricular function by magnetic resonance. British Heart Journal. 1986; 56(4):334-340
Wael M, Ibrahim E-SH, Fahmy AS. Detection of cardiac function abnormality from mri images using normalized wall thickness temporal patterns. International Journal of Biomedical Imaging. 2016; 2016
Ibrahim E-SH. Myocardial tagging by cardiovascular magnetic resonance: Evolution of techniques–pulse sequences, analysis algorithms, and applications. Journal of Cardiovascular Magnetic Resonance. 2011; 13(1):1
Okrasinski SJ, Ramachandran B, Konofagou EE. Assessment of myocardial elastography performance in phantoms under combined physiologic motion configurations with preliminary in vivo feasibility. Physics in Medicine and Biology. 2012; 57(17):5633
Dandel M, Hetzer R. Echocardiographic strain and strain rate imaging clinical applications. International Journal of Cardiology. 2009; 132(1):11-24
Cerqueira MD, Weissman NJ, Dilsizian V, Jacobs AK, Kaul S, Laskey WK, Pennell DJ, Rumberger JA, Ryan T, Verani MS, et al. Standardized myocardial segmentation and nomenclature for tomographic imaging of the heart a statement for healthcare professionals from the cardiac imaging committee of the council on clinical cardiology of the american heart association. Circulation. 2002; 105(4):539-542
McVeigh ER, Prinzen FW, Wyman BT, Tsitlik JE, Halperin HR, Hunter WC. Imaging asynchronous mechanical activation of the paced heart with tagged MRI. Magnetic Resonance in Medicine. 1998; 39(4):507-513
Wyman BT, Hunter WC, Prinzen FW, McVeigh ER. Mapping propagation of mechanical activation in the paced heart with mri tagging. American Journal of Physiology-Heart and Circulatory Physiology. 1999; 276(3):H881-H891
Behdadfar S, Navarro L, Sundnes J, Maleckar M, Ross S, Odland HH, Avril S. A centerline based model morphing algorithm for patient-specific finite element modeling of the left ventricle. IEEE Transactions on Biomedical Engineering. 2017
Hansegård J, Urheim S, Lunde K, Malm S, Rabben SI. Semiautomated quantification of left ventricular volumes and ejection fraction by real-time three-dimensional echocardiography. Cardiovascular Ultrasound. 2009; 7(1):1
Huyghe JM, Arts T, van Campen DH, Reneman RS. Porous medium finite element model of the beating left ventricle. American Journal of Physiology-Heart and Circulatory Physiology. 1992; 262(4):H1256-H1267
Bovendeerd PHM, Huyghe JM, Arts T, Van Campen DH, Reneman RS. Influence of endocardial-epicardial crossover of muscle fibers on left ventricular wall mechanics. Journal of Biomechanics. 1994; 27(7):941-951
Behdadfar S, Navarro L, Sundnes J, Maleckar MM, Avril S. Importance of material parameters and strain energy function on the wall stresses in the left ventricle. Computer Methods in Biomechanics and Biomedical Engineering. 2017; 20(11):1223-1232
Walker JC, Ratcliffe MB, Zhang P, Wallace AW, Fata B, Hsu EW, Saloner D, Guccione JM. Mri-based finite-element analysis of left ventricular aneurysm. American Journal of Physiology-Heart and Circulatory Physiology. 2005; 289(2):H692-H700
Dorri F, Niederer PF, Lunkenheimer PP. A finite element model of the human left ventricular systole. Computer Methods in Biomechanics and Biomedical Engineering. 2006; 9(5):319-341
Holzapfel GA, Ogden RW. Constitutive modelling of passive myocardium: A structurally based framework for material characterization. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences. 2009; 367(1902):3445-3475
Guccione JM, Costa KD, McCulloch AD. Finite element stress analysis of left ventricular mechanics in the beating dog heart. Journal of Biomechanics. 1995; 28(10):1167-1177
Usyk TP, Mazhari R, McCulloch AD. Effect of laminar orthotropic myofiber architecture on regional stress and strain in the canine left ventricle. Journal of Elasticity and the Physical Science of Solids. 2000; 61(1–3):143-164
Marchesseau S, Delingette H, Sermesant M, Ayache N. Fast parameter calibration of a cardiac electromechanical model from medical images based on the unscented transform. Biomechanics and Modeling in Mechanobiology. 2013; 12(4):815-831
McCulloch A, Waldman L, Rogers J, Guccione J. Large-scale finite element analysis of the beating heart. Critical Reviews in Biomedical Engineering. 1991; 20(5–6):427-449
Guccione JM, McCulloch AD. Mechanics of active contraction in cardiac muscle: Part I. Constitutive relations for fiber stress that describe deactivation. Journal of Biomechanical Engineering. 1993; 115(1):72-81
Hunter PJ, McCulloch AD, Ter Keurs HEDJ. Modelling the mechanical properties of cardiac muscle. Progress in Biophysics and Molecular Biology. 1998; 69(2):289-331
Holmes JW, Nunez JA, Covell JW. Functional implications of myocardial scar structure. American Journal of Physiology-Heart and Circulatory Physiology. 1997; 272(5):H2123-H2130
Holmes JW, Borg TK, Covell JW. Structure and mechanics of healing myocardial infarcts. Annu. Rev. Biomed. Eng. 2005; 7:223-253
Walker JC, Guccione JM, Jiang Y, Zhang P, Wallace AW, Hsu EW, Ratcliffe MB. Helical myofiber orientation after myocardial infarction and left ventricular surgical restoration in sheep. The Journal of Thoracic and Cardiovascular Surgery. 2005; 129(2):382-390
Chen J, Song S-K, Liu W, McLean M, Stacy Allen J, Tan J, Wickline SA, Xin Y. Remodeling of cardiac fiber structure after infarction in rats quantified with diffusion tensor MRI. American Journal of Physiology-Heart and Circulatory Physiology. 2003; 285(3):H946-H954
Heimdal A. 4d Strain: Advanced Research Application for Quantitative Echocardiography. GE Healthcare: United Kingdom; 2011
Walker JC, Ratcliffe MB, Zhang P, Wallace AW, Hsu EW, Saloner DA, Guccione JM. Magnetic resonance imaging-based finite element stress analysis after linear repair of left ventricular aneurysm. The Journal of Thoracic and Cardiovascular Surgery. 2008; 135(5):1094-1102
Arevalo HJ, Vadakkumpadan F, Guallar E, Jebb A, Malamas P, Wu KC, Trayanova NA. Arrhythmia risk stratification of patients after myocardial infarction using personalized heart models. Nature Communications. 2016; 7
Pop M, Sermesant M, Mansi T, Crystal E, Detsky J, Yang Y, Fefer P, McVeigh ER, Dick A, Ayache N, et al. Characterization of post-infarct scars in a porcine model—A combined experimental and theoretical study. In: International Conference on Functional Imaging and Modeling of the Heart. Springer; 2009. pp. 1-10
Wenk JF, Sun K, Zhang Z, Soleimani M, Ge L, Saloner D, Wallace AW, Ratcliffe MB, Guccione JM. Regional left ventricular myocardial contractility and stress in a finite element model of posterobasal myocardial infarction. Journal of Biomechanical Engineering. 2011; 133(4):044501
Novak VP, Yin FCP, Humphrey JD. Regional mechanical properties of passive myocardium. Journal of Biomechanics. 1994; 27(4):403-412
Kelley ST, Malekan R, Gorman JH, Jackson BM, Gorman RC, Suzuki Y, Plappert T, Bogen DK, Martin G, Sutton SJ, Henry Edmunds L. Restraining infarct expansion preserves left ventricular geometry and function after acute anteroapical infarction. Circulation. 1999; 99(1):135-142
Götte MJW, van Rossum AC, Twisk JWR, Kuijer JPA, Tim Marcus J, Visser CA. Quantification of regional contractile function after infarction: Strain analysis superior to wall thickening analysis in discriminating infarct from remote myocardium. Journal of the American College of Cardiology. 2001; 37(3):808-817
Denney TS, Gerber BL, Yan L. Unsupervised reconstruction of a three-dimensional left ventricular strain from parallel tagged cardiac images. Magnetic Resonance in Medicine. 2003; 49(4):743-754
Kasai T, Gordon DePuey E, Shah AA. Decreased septal wall thickening in patients with left bundle branch block. Journal of Nuclear Cardiology. 2004; 11(1):32-37
Mahrholdt H, Zhydkov A, Hager S, Meinhardt G, Vogelsberg H, Wagner A, Sechtem U. Left ventricular wall motion abnormalities as well as reduced wall thickness can cause false positive results of routine spect perfusion imaging for detection of myocardial infarction. European Heart Journal. 2005; 26(20):2127-2135
Sugihara H, Tamaki N, Nozawa M, Ohmura T, et al. Septal perfusion and wall thickening in patients with left bundle branch block assessed by technetium-99m-sestamibi gated tomography. The. Journal of Nuclear Medicine. 1997; 38(4):545
Vernooy K, Verbeek XAAM, Peschar M, Crijns HJGM, Arts T, Cornelussen RNM, Prinzen FW. Left bundle branch block induces ventricular remodelling and functional septalhypoperfusion. European Heart Journal. 2005; 26(1):91-98
Prinzen FW, Cheriex EC, Delhaas T, van Oosterhout MFM, Arts T, Wellens HJJ, Reneman RS. Asymmetric thickness of the left ventricular wall resulting from asynchronous electric activation: A study in dogs with ventricular pacing and in patients with left bundle branch block. American Heart Journal. 1995; 130(5):1045-1053
Breithardt O-A, Stellbrink C, Herbots L, Claus P, Sinha AM, Bijnens B, Hanrath P, Sutherland GR. Cardiac resynchronization therapy can reverse abnormal myocardial strain distribution in patients with heart failure and left bundle branch block. Journal of the American College of Cardiology. 2003; 42(3):486-494
Ypenburg C, van Bommel RJ, Jan Willem Borleffs C, Bleeker GB, Boersma E, Schalij MJ, Bax JJ. Long-term prognosis after cardiac resynchronization therapy is related to the extent of left ventricular reverse remodeling at midterm follow-up. Journal of the American College of Cardiology. 2009; 53(6):483-490
Auricchio A, Stellbrink C, Sack S, Block M, Vogt J, Bakker P, Huth C, Schöndube F, Wolfhard U, Böcker D, et al. Long-term clinical effect of hemodynamically optimized cardiac resynchronization therapy in patients with heart failure and ventricular conduction delay. Journal of the American College of Cardiology. 2002; 39(12):2026-2033
Abraham WT, Fisher WG, Smith AL, Delurgio DB, Leon AR, Loh E, Kocovic DZ, Packer M, Clavell AL, Hayes DL, et al. Cardiac resynchronization in chronic heart failure. New England Journal of Medicine. 2002; 346(24):1845-1853
Kerckhoffs RCP, McCulloch AD, Omens JH, Mulligan LJ. Effects of biventricular pacing and scar size in a computational model of the failing heart with left bundle branch block. Medical Image Analysis. 2009; 13(2):362-369
Constantino J, Hu Y, Trayanova NA. A computational approach to understanding the cardiac electromechanical activation sequence in the normal and failing heart, with translation to the clinical practice of crt. Progress in Biophysics and Molecular Biology. 2012; 110(2):372-379
Ming-Ting W, Tseng W-YI, Su M-YM, Liu C-P, Chiou K-R, Wedeen VJ, Reese TG, Yang C-F. Diffusion tensor magnetic resonance imaging mapping the fiber architecture remodeling in human myocardium after infarction correlation with viability and wall motion. Circulation. 2006; 114(10):1036-1045