Advances in computational modelling for personalised medicine after myocardial infarction

Myocardial infarction (MI) is a leading cause of premature morbidity and mortality worldwide. Determining which patients will experience heart failure and sudden cardiac death after an acute MI is notoriously difficult for clinicians. The extent of heart damage after an acute MI is informed by cardiac imaging, typically using echocardiography or sometimes, cardiac magnetic resonance (CMR). These scans provide complex data sets that are only partially exploited by clinicians in daily practice, implying potential for improved risk assessment. Computational modelling of left ventricular (LV) function can bridge the gap towards personalised medicine using cardiac imaging in patients with post-MI. Several novel biomechanical parameters have theoretical prognostic value and may be useful to reflect the biomechanical effects of novel preventive therapy for adverse remodelling post-MI. These parameters include myocardial contractility (regional and global), stiffness and stress. Further, the parameters can be delineated spatially to correspond with infarct pathology and the remote zone. While these parameters hold promise, there are challenges for translating MI modelling into clinical practice, including model uncertainty, validation and verification, as well as time-efficient processing. More research is needed to (1) simplify imaging with CMR in patients with post-MI, while preserving diagnostic accuracy and patient tolerance (2) to assess and validate novel biomechanical parameters against established prognostic biomarkers, such as LV ejection fraction and infarct size. Accessible software packages with minimal user interaction are also needed. Translating benefits to patients will be achieved through a multidisciplinary approach including clinicians, mathematicians, statisticians and industry partners

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

Myocardial tagging by Cardiovascular Magnetic Resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications

Cardiovascular magnetic resonance (CMR) tagging has been established as an essential technique for measuring regional myocardial function. It allows quantification of local intramyocardial motion measures, e.g. strain and strain rate. The invention of CMR tagging came in the late eighties, where the technique allowed for the first time for visualizing transmural myocardial movement without having to implant physical markers. This new idea opened the door for a series of developments and improvements that continue up to the present time. Different tagging techniques are currently available that are more extensive, improved, and sophisticated than they were twenty years ago. Each of these techniques has different versions for improved resolution, signal-to-noise ratio (SNR), scan time, anatomical coverage, three-dimensional capability, and image quality. The tagging techniques covered in this article can be broadly divided into two main categories: 1) Basic techniques, which include magnetization saturation, spatial modulation of magnetization (SPAMM), delay alternating with nutations for tailored excitation (DANTE), and complementary SPAMM (CSPAMM); and 2) Advanced techniques, which include harmonic phase (HARP), displacement encoding with stimulated echoes (DENSE), and strain encoding (SENC). Although most of these techniques were developed by separate groups and evolved from different backgrounds, they are in fact closely related to each other, and they can be interpreted from more than one perspective. Some of these techniques even followed parallel paths of developments, as illustrated in the article. As each technique has its own advantages, some efforts have been made to combine different techniques together for improved image quality or composite information acquisition. In this review, different developments in pulse sequences and related image processing techniques are described along with the necessities that led to their invention, which makes this article easy to read and the covered techniques easy to follow. Major studies that applied CMR tagging for studying myocardial mechanics are also summarized. Finally, the current article includes a plethora of ideas and techniques with over 300 references that motivate the reader to think about the future of CMR tagging

Flow pattern analysis for magnetic resonance velocity imaging

Blood flow in the heart is highly complex. Although blood flow patterns have been investigated by both computational modelling and invasive/non-invasive imaging techniques, their evolution and intrinsic connection with cardiovascular disease has yet to be explored. Magnetic resonance (MR) velocity imaging provides a comprehensive distribution of multi-directional in vivo flow distribution so that detailed quantitative analysis of flow patterns is now possible. However, direct visualisation or quantification of vector fields is of little clinical use, especially for inter-subject or serial comparison of changes in flow patterns due to the progression of the disease or in response to therapeutic measures. In order to achieve a comprehensive and integrated description of flow in health and disease, it is necessary to characterise and model both normal and abnormal flows and their effects. To accommodate the diversity of flow patterns in relation to morphological and functional changes, we have described in this thesis an approach of detecting salient topological features prior to analytical assessment of dynamical indices of the flow patterns. To improve the accuracy of quantitative analysis of the evolution of topological flow features, it is essential to restore the original flow fields so that critical points associated with salient flow features can be more reliably detected. We propose a novel framework for the restoration, abstraction, extraction and tracking of flow features such that their dynamic indices can be accurately tracked and quantified. The restoration method is formulated as a constrained optimisation problem to remove the effects of noise and to improve the consistency of the MR velocity data. A computational scheme is derived from the First Order Lagrangian Method for solving the optimisation problem. After restoration, flow abstraction is applied to partition the entire flow field into clusters, each of which is represented by a local linear expansion of its velocity components. This process not only greatly reduces the amount of data required to encode the velocity distribution but also permits an analytical representation of the flow field from which critical points associated with salient flow features can be accurately extracted. After the critical points are extracted, phase portrait theory can be applied to separate them into attracting/repelling focuses, attracting/repelling nodes, planar vortex, or saddle. In this thesis, we have focused on vortical flow features formed in diastole. To track the movement of the vortices within a cardiac cycle, a tracking algorithm based on relaxation labelling is employed. The constraints and parameters used in the tracking algorithm are designed using the characteristics of the vortices. The proposed framework is validated with both simulated and in vivo data acquired from patients with sequential MR examination following myocardial infarction. The main contribution of the thesis is in the new vector field restoration and flow feature abstraction method proposed. They allow the accurate tracking and quantification of dynamic indices associated with salient features so that inter- and intra-subject comparisons can be more easily made. This provides further insight into the evolution of blood flow patterns and permits the establishment of links between blood flow patterns and localised genesis and progression of cardiovascular disease.Open acces

Imaging the right heart: the use of integrated multimodality imaging

During recent years, right ventricular (RV) structure and function have been found to be an important determinant of outcome in different cardiovascular and also pulmonary diseases. Currently, echocardiography and cardiac magnetic resonance (CMR) imaging are the two imaging modalities most commonly used to visualize the RV. Most structural abnormalities of the RV can be reliably described by echocardiography but due its complex geometrical shape, echocardiographic assessment of RV function is more challenging. Newer promising echocardiographic techniques are emerging but lack of validation and limited normal reference data influence their routine clinical application. Cardiac magnetic resonance is generally considered the clinical reference technique due to its unlimited imaging planes, superior image resolution, and three-dimensional volumetric rendering. The accuracy and reliability of CMR measurements make it the ideal tool for serial examinations of RV function. Multidetector computed tomography (MDCT) plays an important role in the diagnosis of pulmonary emboli but can also be used for assessing RV ischaemic disease or as an alternative for CMR if contra-indicated. Radionuclide techniques have become more obsolete in the current era. The different imaging modalities should be considered complimentary and each plays a role for different indication

Analytical method to measure three-dimensional strain patterns in the left ventricle from single slice displacement data

Background: Displacement encoded Cardiovascular MR (CMR) can provide high spatial resolution measurements of three-dimensional (3D) Lagrangian displacement. Spatial gradients of the Lagrangian displacement field are used to measure regional myocardial strain. In general, adjacent parallel slices are needed in order to calculate the spatial gradient in the through-slice direction. This necessitates the acquisition of additional data and prolongs the scan time. The goal of this study is to define an analytic solution that supports the reconstruction of the out-of-plane components of the Lagrangian strain tensor in addition to the in-plane components from a single-slice displacement CMR dataset with high spatio-temporal resolution. The technique assumes incompressibility of the myocardium as a physical constraint. Results: The feasibility of the method is demonstrated in a healthy human subject and the results are compared to those of other studies. The proposed method was validated with simulated data and strain estimates from experimentally measured DENSE data, which were compared to the strain calculation from a conventional two-slice acquisition. Conclusion: This analytical method reduces the need to acquire data from adjacent slices when calculating regional Lagrangian strains and can effectively reduce the long scan time by a factor of two

Left ventricular remodeling and function in ischemic heart disease and aortic valve disease

Background: Cardiac remodeling is a broad term that refers to structural and functional alterations of the heart in response to chronic changes in loading conditions or left ventricular (LV) contractile performance. Different loading conditions will affect the heart in different ways, some leading to impaired heart function, symptoms of heart failure, or even death. However, the process of remodeling may not be permanent. If the heart is relieved of the underlying cause of the remodeling, the heart function and structure may normalize in a process referred to as reverse remodeling. The complex interplay of factors that determine the process of reverse remodeling is not fully elucidated. Cardiac remodeling can be evaluated by many different diagnostic modalities, but the most widely used diagnostic tool is two-dimensional echocardiography (2DE). In recent years, three-dimensional echocardiography (3DE) has emerged with possible advantages in the assessment of LV volume and function. The thesis aimed to evaluate 3DE in the assessment of LV function and remodeling, and to study different aspects of remodeling in response to pressure and volume overload in patients with aortic stenosis (AS) and aortic regurgitation (AR), respectively. Methods: Studies I and II investigated patients with ischemic heart disease (n = 15 and n = 32, respectively). In Study I, the assessments of LV volume and ejection fraction (EF) were compared using 3DE, cardiac magnetic resonance (CMR), and single-photon emission computer tomography (SPECT). Study II compared the performance of 2DE, contrast-enhanced 2DE, 3DE, and contrast-enhanced 3DE in the assessment LV volumes and EF, using CMR as a reference standard. In Studies III and IV, 65 patients with severe AR and 120 patients with severe AS, respectively, were examined using 2DE and 3DE before and at one year after aortic valve replacement (AVR). In Study III, LV volumes, systolic and diastolic LV function, and left atrial strain (LAS) were analyzed to identify predictors of impaired LV reverse remodeling in AR. Study IV assessed LV functional indices, including 2D global longitudinal strain (GLS) and 3D strain, to assess predictors of incomplete reverse remodeling in AS. Results and conclusions: There were significant differences among 3DE, SPECT and CMR regarding the measurement of LV volumes. However, the estimation of EF showed good agreement. 3DE was more accurate and showed more favorable reproducibility than 2DE for the assessment of EF and LV volumes. Contrast enhancement improved accuracy and reproducibility for both 2DE and 3DE. One-third of patients with AR had signs of impaired LV diastolic function. After AVR, diastolic LV functional indices improved, LV and left atrial (LA) volumes decreased, and indices of LA function increased. LA conduit strain had an incremental prognostic value for the prediction of impaired LV functional and structural recovery. In patients with AS, AVR was associated with a decrease in LV mass, an improvement in 2D GLS, and a decrease in LV twist. 2D GLS and left ventricular mass index were predictive of incomplete reverse remodeling during the follow-up period. 3D GLS did not add discriminatory or predictive information over 2D GLS