Abnormal Tissue Zone Detection and Average Active Stress Estimation in Patients with LV Dysfunction

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

Insights from echocardiography, magnetic resonance imaging, and microcomputed tomography relative to the mid-myocardial left ventricular echogenic zone.

BACKGROUND: The anatomical substrate for the mid-mural ventricular hyperechogenic zone remains uncertain, but it may represent no more than ultrasound reflected from cardiomyocytes orientated orthogonally to the ultrasonic beam. We sought to ascertain the relationship between the echogenic zone and the orientation of the cardiomyocytes. METHODS: We used 3D echocardiography, diffusion tensor imaging, and microcomputed tomography to analyze the location and orientation of cardiomyocytes within the echogenic zone. RESULTS: We demonstrated that visualization of the echogenic zone is dependent on the position of the transducer and is most clearly seen from the apical window. Diffusion tensor imaging and microcomputed tomography show that the echogenic zone seen from the apical window corresponds to the position of the circumferentially orientated cardiomyocytes. An oblique band seen in the parasternal view relates to cardiomyocytes orientated orthogonally to the ultrasonic beam. CONCLUSIONS: The mid-mural ventricular hyperechogenic zone represents reflected ultrasound from cardiomyocytes aligned orthogonal to the ultrasonic beam. The echogenic zone does not represent a space, a connective tissue sheet, a boundary between ascending and descending limbs of a hypothetical helical ventricular myocardial band, nor an abrupt change in cardiomyocyte orientation

Towards automating cine DENSE MRI image analysis : segmentation, tissue tracking and strain computation

Includes bibliographical references (p. 192-206).Over the past two decades, magnetic resonance imaging (MRI) has developed into a powerful imaging tool for the heart. Imaging cardiac morphology is now commonplace in clinical practice, and a plethora of quantitative techniques have also arisen on the research front. Myocardial tagging is an established quantitative cardiac MRI method that involves magnetically tagging the heart with a set of saturated bands, and monitoring the deformation of these bands as the heart contracts

Computational biomechanics of acute myocardial infarction and its treatment

The intramyocardial injection of biomaterials is an emerging therapy for myocardial infarction. Computational methods can help to study the mechanical effect s of biomaterial injectates on the infarcted heart s and can contribute to advance and optimise the concept of this therapy. The distribution of polyethylene glycol hydrogel injectate delivered immediately after the infarct induction was studied using rat infarct model. A micro-structural three-dimensional geometrical model of the entire injectate was reconstructed from histological micro graphs. The model provides a realistic representation of biomaterial injectates in computational models at macroscopic and microscopic level. Biaxial and compression mechanical testing was conducted for healing rat myocardial infarcted tissue at immediate (0 day), 7, 14 and 28 days after infarction onset. Infarcts were found to be mechanically anisotropic with the tissue being stiffer in circumferential direction than in longitudinal direction. The 0, 7, 14 and 28 days infarcts showed 443, 670, 857 and 1218 kPa circumferential tensile moduli. The 28 day infarct group showed a significantly higher compressive modulus compared to the other infarct groups (p= 0.0055, 0.028, and 0.018 for 0, 7 and 14 days groups). The biaxial mechanical data were utilized to establish material constitutive models of rat healing infarcts. Finite element model s and genetic algorithms were employed to identify the parameters of Fung orthotropic hyperelastic strain energy function for the healing infarcts. The provided infarct mechanical data and the identified constitutive parameters offer a platform for investigations of mechanical aspects of myocardial infarction and therapies in the rat, an experimental model extensively used in the development of infarct therapies. Micro-structurally detailed finite element model of a hydrogel injectate in an infarct was developed to provide an insight into the micromechanics of a hydrogel injectate and infarct during the diastolic filling. The injectate caused the end-diastolic fibre stresses in the infarct zone to decrease from 22.1 to 7.7 kPa in the 7 day infarct and from 35.7 to 9.7 kPa in the 28 day infarct. This stress reduction effect declined as the stiffness of the biomaterial increased. It is suggested that the gel works as a force attenuating system through micromechanical mechanisms reducing the force acting on tissue layers during the passive diastolic dilation of the left ventricle and thus reducing the stress induced in these tissue layers

Novel cardiovascular magnetic resonance phenotyping of the myocardium

INTRODUCTION Left ventricular (LV) microstructure is unique, composed of a winding helical pattern of myocytes and rotating aggregations of myocytes called sheetlets. Hypertrophic cardiomyopathy (HCM) is a cardiovascular disease characterised by left ventricular hypertrophy (LVH), however the link between LVH and underlying microstructural aberration is poorly understood. In vivo cardiovascular diffusion tensor imaging (cDTI) is a novel cardiovascular MRI (CMR) technique, capable of characterising LV microstructural dynamics non-invasively. In vivo cDTI may therefore improve our understanding microstructural-functional relationships in health and disease. METHODS AND RESULTS The monopolar diffusion weighted stimulated echo acquisition mode (DW-STEAM) sequence was evaluated for in vivo cDTI acquisitions at 3Tesla, in healthy volunteers (HV), patients with hypertensive LVH, and HCM patients. Results were contextualised in relation to extensively explored technical limitations. cDTI parameters demonstrated good intra-centre reproducibility in HCM, and good inter-centre reproducibility in HV. In all subjects, cDTI was able to depict the winding helical pattern of myocyte orientation known from histology, and the transmural rate of change in myocyte orientation was dependent on LV size and thickness. In HV, comparison of cDTI parameters between systole and diastole revealed an increase in transmural gradient, combined with a significant re-orientation of sheetlet angle. In contrast, in HCM, myocyte gradient increased between phases, however sheetlet angulation retained a systolic-like orientation in both phases. Combined analysis with hypertensive patients revealed a proportional decrease in sheetlet mobility with increasing LVH. CONCLUSION In vivo DW-STEAM cDTI can characterise LV microstructural dynamics non-invasively. The transmural rate of change in myocyte angulation is dependent on LV size and wall thickness, however inter phase changes in myocyte orientation are unaffected by LVH. In contrast, sheetlet dynamics demonstrate increasing dysfunction, in proportion to the degree of LVH. Resolving technical limitations is key to advancing this technique, and improving the understanding of the role of microstructural abnormalities in cardiovascular disease expression.Open Acces

Considering New Regularization Parameter-Choice Techniques for the Tikhonov Method to Improve the Accuracy of Electrocardiographic Imaging

The electrocardiographic imaging (ECGI) inverse problem highly relies on adding constraints, a process called regularization, as the problem is ill-posed. When there are no prior information provided about the unknown epicardial potentials, the Tikhonov regularization method seems to be the most commonly used technique. In the Tikhonov approach the weight of the constraints is determined by the regularization parameter. However, the regularization parameter is problem and data dependent, meaning that different numerical models or different clinical data may require different regularization parameters. Then, we need to have as many regularization parameter-choice methods as techniques to validate them. In this work, we addressed this issue by showing that the Discrete Picard Condition (DPC) can guide a good regularization parameter choice for the two-norm Tikhonov method. We also studied the feasibility of two techniques: The U-curve method (not yet used in the cardiac field) and a novel automatic method, called ADPC due its basis on the DPC. Both techniques were tested with simulated and experimental data when using the method of fundamental solutions as a numerical model. Their efficacy was compared with the efficacy of two widely used techniques in the literature, the L-curve and the CRESO methods. These solutions showed the feasibility of the new techniques in the cardiac setting, an improvement of the morphology of the reconstructed epicardial potentials, and in most of the cases of their amplitude