Direct histological validation of diffusion tensor MRI in formaldehyde-fixed myocardium

Investigating the three-dimensional myocardial micro-architecture in the laminar structure using X-ray phase-contrast microtomography

A comprehensive grasp of the myocardial micro-architecture is essential for understanding diverse heart functions. This study aimed to investigate three-dimensional (3D) cardiomyocyte arrangement in the laminar structure using X-ray phase-contrast microtomography. Using the ID-19 beamline at the European Synchrotron Radiation Facility, we imaged human left ventricular (LV) wall transparietal samples and reconstructed them with an isotropic voxel edge length of 3.5 μm. From the reconstructed volumes, we extracted different regions to analyze the orientation distribution of local cardiomyocyte aggregates, presenting findings in terms of helix and intrusion angles. In regions containing one sheetlet population, we observed cardiomyocyte aggregates running along the local LV wall's radial direction at the border of sheetlets, branching and merging into a complex network around connecting points of different sheetlets, and bending to accommodate vessel passages. In regions with two sheetlet populations, the helix angle of local cardiomyocyte aggregates experiences a nonmonotonic change, and some cardiomyocyte aggregates run along the local radial direction. X-ray phase-contrast microtomography is a valuable technique for investigating the 3D local myocardial architecture at microscopic level. The arrangement of local cardiomyocyte aggregates in the LV wall proves to be both regional and complex, intricately linked to the local laminar structure.

Using diffusion tensor imaging to depict myocardial changes after matured pluripotent stem cell-derived cardiomyocyte transplantation

BACKGROUND

Novel treatment strategies are needed to improve the structure and function of the myocardium post-infarction. In vitro-matured pluripotent stem cell-derived cardiomyocytes (PSC-CMs) have been shown to be a promising regenerative strategy. We hypothesized that mature PSC-CMs will have anisotropic structure and improved cell alignment when compared to immature PSC-CMs using cardiovascular magnetic resonance (CMR) in a guinea pig model of cardiac injury.

METHODS

Guinea pigs (n = 16) were cryoinjured on day -10, followed by transplantation of either 108 polydimethylsiloxane (PDMS)-matured PSC-CMs (n = 6) or 108 immature tissue culture plastic (TCP)-generated PSC-CMs (n = 6) on day 0. Vehicle (sham-treated) subjects were injected with a pro-survival cocktail devoid of cells (n = 4), while healthy controls (n = 4) did not undergo cryoinjury or treatment. Animals were sacrificed on either day +14 or day +28 post-transplantation. Animals were imaged ex vivo on a 7T Bruker MRI. A 3D diffusion tensor imaging (DTI) sequence was used to quantify structure via fractional anisotropy (FA), mean diffusivity (MD), and myocyte alignment measured by the standard deviation of the transverse angle (TA).

RESULTS

MD and FA of mature PDMS grafts demonstrated anisotropy was not significantly different than the healthy control hearts (MD = 1.1 ± 0.12 × 10-3 mm2/s vs 0.93 ± 0.01 × 10-3 mm2/s, p = 0.4 and FA = 0.22 ± 0.05 vs 0.26 ± 0.001, p = 0.5). Immature TCP grafts exhibited significantly higher MD than the healthy control (1.3 ± 0.08 × 10-3 mm2/s, p < 0.05) and significantly lower FA than the control (0.12 ± 0.02, p < 0.05) but were not different from mature PDMS grafts in this small cohort. TA of healthy controls showed low variability and was not significantly different than mature PDMS grafts (p = 0.4) while immature TCP grafts were significantly different (p < 0.001). DTI parameters of mature graft tissue trended toward that of the healthy myocardium, indicating the grafted cardiomyocytes may have a similar phenotype to healthy tissue. Contrast-enhanced magnetic resonance images corresponded well to histological staining, demonstrating a non-invasive method of localizing the repopulated cardiomyocytes within the scar.

CONCLUSIONS

The DTI measures within graft tissue were indicative of anisotropic structure and showed greater myocyte organization compared to the scarred territory. These findings show that MRI is a valuable tool to assess the structural impacts of regenerative therapies.

A computational study of right ventricular mechanics in a rat model of pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) presents a significant challenge to right ventricular (RV) function due to progressive pressure overload, necessitating adaptive remodeling in the form of increased wall thickness, enhanced myocardial contractility and stiffness to maintain cardiac performance. However, the impact of these remodeling mechanisms on RV mechanics in not clearly understood. In addition, there is a lack of quantitative understanding of how each mechanism individually influences RV mechanics. Utilizing experimental data from a rat model of PAH at three distinct time points, we developed biventricular finite element models to investigate how RV stress and strain evolved with PAH progression. The finite element models were fitted to hemodynamic and morphological data to represent different disease stages and used to analyze the impact of RV remodeling as well as the altered RV pressure. Furthermore, we performed a number of theoretical simulation studies with different combinations of morphological and physiological remodeling, to assess and quantify their individual impact on overall RV load and function. Our findings revealed a substantial 4-fold increase in RV stiffness and a transient 2-fold rise in contractility, which returned to baseline by week 12. These changes in RV material properties in addition to the 2-fold increase in wall thickness significantly mitigated the increase in wall stress and strain caused by the progressive increase in RV afterload. Despite the PAH-induced cases showing increased wall stress and strain at end-diastole and end-systole compared to the control, our simulations suggest that without the observed remodeling mechanisms, the increase in stress and strain would have been much more pronounced. Our model analysis also indicated that while changes in the RV's material properties-particularly increased RV stiffness - have a notable effect on its mechanics, the primary compensatory factor limiting the stress and strain increase in the early stages of PAH was the significant increase in wall thickness. These findings underscore the importance of RV remodeling in managing the mechanical burden on the right ventricle due to pressure overload.

A simple approach for image-based modelling of the heart that enables robust simulation of highly heterogeneous electrical excitation

Remodelling of cardiac tissue structure, including intercellular electrical coupling, is a major determinant of the complex and heterogeneous excitation patterns associated with cardiac arrhythmias. Evaluation of the precise mechanisms by which local tissue structure determines global arrhythmic excitation patterns is a major challenge that may be critically important for the development of effective treatment strategies. Computational modelling is a key tool in the study of cardiac arrhythmias, yet the established approaches for organ-scale modelling are unsuitable to capture the impact of local conduction heterogeneities; a novel approach is required to provide this multi-scale mechanistic insight. We present a fundamentally simple yet powerful approach to simulate electrical excitation in highly heterogeneous whole-heart models that exploits the underlying discreteness of the myocardium. Preliminary simulations demonstrate that this approach can capture lower conduction velocities and reproduce wave breakdown and the development of re-entry in a range of conditions.

The relationship between myocardial microstructure and strain in chronic infarction using cardiovascular magnetic resonance diffusion tensor imaging and feature tracking

BACKGROUND

Cardiac diffusion tensor imaging (cDTI) using cardiovascular magnetic resonance (CMR) is a novel technique for the non-invasive assessment of myocardial microstructure. Previous studies have shown myocardial infarction to result in loss of sheetlet angularity, derived by reduced secondary eigenvector (E2A) and reduction in subendocardial cardiomyocytes, evidenced by loss of myocytes with right-handed orientation (RHM) on helix angle (HA) maps. Myocardial strain assessed using feature tracking-CMR (FT-CMR) is a sensitive marker of sub-clinical myocardial dysfunction. We sought to explore the relationship between these two techniques (strain and cDTI) in patients at 3 months following ST-elevation MI (STEMI).

METHODS

32 patients (F = 28, 60 ± 10 years) underwent 3T CMR three months after STEMI (mean interval 105 ± 17 days) with second order motion compensated (M2), free-breathing spin echo cDTI, cine gradient echo and late gadolinium enhancement (LGE) imaging. HA maps divided into left-handed HA (LHM, - 90 < HA < - 30), circumferential HA (CM, - 30° < HA < 30°), and right-handed HA (RHM, 30° < HA < 90°) were reported as relative proportions. Global and segmental analysis was undertaken.

RESULTS

Mean left ventricular ejection fraction (LVEF) was 44 ± 10% with a mean infarct size of 18 ± 12 g and a mean infarct segment LGE enhancement of 66 ± 21%. Mean global radial strain was 19 ± 6, mean global circumferential strain was - 13 ± - 3 and mean global longitudinal strain was - 10 ± - 3. Global and segmental radial strain correlated significantly with E2A in infarcted segments (p = 0.002, p = 0.011). Both global and segmental longitudinal strain correlated with RHM of infarcted segments on HA maps (p < 0.001, p = 0.003). Mean Diffusivity (MD) correlated significantly with the global infarct size (p < 0.008). When patients were categorised according to LVEF (reduced, mid-range and preserved), all cDTI parameters differed significantly between the three groups.

CONCLUSION

Change in sheetlet orientation assessed using E2A from cDTI correlates with impaired radial strain. Segments with fewer subendocardial cardiomyocytes, evidenced by a lower proportion of myocytes with right-handed orientation on HA maps, show impaired longitudinal strain. Infarct segment enhancement correlates significantly with E2A and RHM. Our data has demonstrated a link between myocardial microstructure and contractility following myocardial infarction, suggesting a potential role for CMR cDTI to clinically relevant functional impact.

Three-dimensional micro-structurally informed in silico myocardium-Towards virtual imaging trials in cardiac diffusion weighted MRI

In silico tissue models (viz. numerical phantoms) provide a mechanism for evaluating quantitative models of magnetic resonance imaging. This includes the validation and sensitivity analysis of imaging biomarkers and tissue microstructure parameters. This study proposes a novel method to generate a realistic numerical phantom of myocardial microstructure. The proposed method extends previous studies by accounting for the variability of the cardiomyocyte shape, water exchange between the cardiomyocytes (intercalated discs), disorder class of myocardial microstructure, and four sheetlet orientations. In the first stage of the method, cardiomyocytes and sheetlets are generated by considering the shape variability and intercalated discs in cardiomyocyte-cardiomyocyte connections. Sheetlets are then aggregated and oriented in the directions of interest. The morphometric study demonstrates no significant difference (p>0.01) between the distribution of volume, length, and primary and secondary axes of the numerical and real (literature) cardiomyocyte data. Moreover, structural correlation analysis validates that the in-silico tissue is in the same class of disorderliness as the real tissue. Additionally, the absolute angle differences between the simulated helical angle (HA) and input HA (reference value) of the cardiomyocytes (4.3°±3.1°) demonstrate a good agreement with the absolute angle difference between the measured HA using experimental cardiac diffusion tensor imaging (cDTI) and histology (reference value) reported by (Holmes et al., 2000) (3.7°±6.4°) and (Scollan et al. 1998) (4.9°±14.6°). Furthermore, the angular distance between eigenvectors and sheetlet angles of the input and simulated cDTI is much smaller than those between measured angles using structural tensor imaging (as a gold standard) and experimental cDTI. Combined with the qualitative results, these results confirm that the proposed method can generate richer numerical phantoms for the myocardium than previous studies.

Myocardial mesostructure and mesofunction

The complex and highly organized structural arrangement of some five billion cardiomyocytes directs the coordinated electrical activity and mechanical contraction of the human heart. The characteristic transmural change in cardiomyocyte orientation underlies base-to-apex shortening, circumferential shortening, and left ventricular torsion during contraction. Individual cardiomyocytes shorten ∼15% and increase in diameter ∼8%. Remarkably, however, the left ventricular wall thickens by up to 30-40%. To accommodate this, the myocardium must undergo significant structural rearrangement during contraction. At the mesoscale, collections of cardiomyocytes are organized into sheetlets, and sheetlet shear is the fundamental mechanism of rearrangement that produces wall thickening. Herein, we review the histological and physiological studies of myocardial mesostructure that have established the sheetlet shear model of wall thickening. Recent developments in tissue clearing techniques allow for imaging of whole hearts at the cellular scale, whereas magnetic resonance imaging (MRI) and computed tomography (CT) can image the myocardium at the mesoscale (100 µm to 1 mm) to resolve cardiomyocyte orientation and organization. Through histology, cardiac diffusion tensor imaging (DTI), and other modalities, mesostructural sheetlets have been confirmed in both animal and human hearts. Recent in vivo cardiac DTI methods have measured reorientation of sheetlets during the cardiac cycle. We also examine the role of pathological cardiac remodeling on sheetlet organization and reorientation, and the impact this has on ventricular function and dysfunction. We also review the unresolved mesostructural questions and challenges that may direct future work in the field.

A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: An application to the ventricular junctions

Background

Knowledge of the normal myocardial–myocyte orientation could theoretically allow the definition of relevant quantitative biomarkers in clinical routine to diagnose heart pathologies. A whole heart diffusion tensor template representative of the global myofiber organization over species is therefore crucial for comparisons across populations. In this study, we developed a groupwise registration and tractography framework to resolve the global myofiber arrangement of large mammalian sheep hearts. To demonstrate the potential application of the proposed method, a novel description of sub-regions in the intraventricular septum is presented.

Methods

Three explanted sheep (ovine) hearts (size ~12×8×6 cm3, heart weight ~ 150 g) were perfused with contrast agent and fixative and imaged in a 9.4T magnet. A group-wise registration of high-resolution anatomical and diffusion-weighted images were performed to generate anatomical and diffusion tensor templates. Diffusion tensor metrics (eigenvalues, eigenvectors, fractional anisotropy …) were computed to provide a quantitative and spatially-resolved analysis of cardiac microstructure. Then tractography was performed using deterministic and probabilistic algorithms and used for different purposes: i) Visualization of myofiber architecture, ii) Segmentation of sub-area depicting the same fiber organization, iii) Seeding and Tract Editing. Finally, dissection was performed to confirm the existence of macroscopic structures identified in the diffusion tensor template.

Results

The template creation takes advantage of high-resolution anatomical and diffusion-weighted images obtained at an isotropic resolution of 150 μm and 600 μm respectively, covering ventricles and atria and providing information on the normal myocardial architecture. The diffusion metric distributions from the template were found close to the one of the individual samples validating the registration procedure. Small new sub-regions exhibiting spatially sharp variations in fiber orientation close to the junctions of the septum and ventricles were identified. Each substructure was defined and represented using streamlines. The existence of a fiber-bundles in the posterior junction was validated by anatomical dissection. A complex structural organization of the anterior junction in comparison to the posterior junction was evidenced by the high-resolution acquisition.

Conclusions

A new framework combining cardiac template generation and tractography was applied on the whole sheep heart. The framework can be used for anatomical investigation, characterization of microstructure and visualization of myofiber orientation across samples. Finally, a novel description of the ventricular junction in large mammalian sheep hearts was proposed.

A Transmural Path Model Improves the Definition of the Orthotropic Tissue Structure in Heart Simulations

In the past decades, the structure of the heart, human as well as other species, has been explored in a detailed way, e.g., via histological studies or diffusion tensor magnetic resonance imaging. Nevertheless, the assignment of the characteristic orthotropic structure in a patient-specific finite element model remains a challenging task. Various types of rule-based models, which define the local fiber and sheet orientation depending on the transmural depth, have been developed. However, the correct assessment of the transmural depth is not trivial. Its accuracy has a substantial influence on the overall mechanical and electrical properties in rule-based models. The main purpose of this study is the development of a finite element-based approach to accurately determine the transmural depth on a general unstructured grid. Instead of directly using the solution of the Laplace problem as the transmural depth, we make use of a well-established model for the assessment of the transmural thickness. It is based on two hyperbolic first-order partial differential equations for the definition of a transmural path, whereby the transmural thickness is defined as the arc length of this path. Subsequently, the transmural depth is determined based on the position on the transmural path. Originally, the partial differential equations were solved via finite differences on structured grids. In order to circumvent the need of two grids and mapping between the structured (to determine the transmural depth) and unstructured (electromechanical heart simulation) grids, we solve the equations directly on the same unstructured tetrahedral mesh. We propose a finite-element-based discontinuous Galerkin approach. Based on the accurate transmural depth, we assign the local material orientation of the orthotropic tissue structure in a usual fashion. We show that this approach leads to a more accurate definition of the transmural depth. Furthermore, for the left ventricle, we propose functions for the transmural fiber and sheet orientation by fitting them to literature-based diffusion tensor magnetic resonance imaging data. The proposed functions provide a distinct improvement compared to existing rules from the literature.

Measurement of local orientation of cardiomyocyte aggregates in human left ventricle free wall samples using X-ray phase-contrast microtomography

Most cardiomyocytes in the left ventricle wall are grouped in aggregates of four to five units that are quasi-parallel to each other. When one or more "cardiomyocyte aggregates" are delimited by two cleavage planes, this defines a "sheetlet" that can be considered as a "work unit" that contributes to the thickening of the wall during the cardiac cycle. In this paper, we introduce the skeleton method to measure the local three-dimensional (3D) orientation of cardiomyocyte aggregates in the sheetlets in three steps: data segmentation; extraction of the skeleton of the sheetlets; and calculation of the local orientation of the cardiomyocyte aggregates inside the sheetlets. These data include a series of virtual tissue volumes and five transmural human left ventricle free wall samples, imaged with 3D synchrotron radiation phase-contrast microtomography, and reconstructed with a 3.5×3.5×3.5μm³ voxel size. We computed the local orientation of the cardiomyocyte aggregates inside the sheetlets with a working window of 112×112×112μm³ in size. These data demonstrate that the skeleton method can provide accurate 3D measurements and reliable screening of the 3D evolution of the orientation of cardiomyocyte aggregates within the sheetlets. We showed that in regions that contain one population of quasi-parallel sheetlets, the orientation of the cardiomyocyte aggregates undergo "oscillations" along the perpendicular direction of the sheetlets. In regions that contain two populations of sheetlets with a different angular range, we demonstrate some discontinuity of the helix angle of the cardiomyocyte aggregates at the interface between the two populations.

Arrhythmogenic potential of myocardial disarray in hypertrophic cardiomyopathy: genetic basis, functional consequences and relation to sudden cardiac death

Myocardial disarray is defined as disorganized cardiomyocyte spatial distribution, with loss of physiological fibre alignment and orientation. Since the first pathological descriptions of hypertrophic cardiomyopathy (HCM), disarray appeared as a typical feature of this condition and sparked vivid debate regarding its specificity to the disease and clinical significance as a diagnostic marker and a risk factor for sudden death. Although much of the controversy surrounding its diagnostic value in HCM persists, it is increasingly recognized that myocardial disarray may be found in physiological contexts and in cardiac conditions different from HCM, raising the possibility that central focus should be placed on its quantity and distribution, rather than a mere presence. While further studies are needed to establish what amount of disarray should be considered as a hallmark of the disease, novel experimental approaches and emerging imaging techniques for the first time allow ex vivo and in vivo characterization of the myocardium to a molecular level. Such advances hold the promise of filling major gaps in our understanding of the functional consequences of myocardial disarray in HCM and specifically on arrhythmogenic propensity and as a risk factor for sudden death. Ultimately, these studies will clarify whether disarray represents a major determinant of the HCM clinical profile, and a potential therapeutic target, as opposed to an intriguing but largely innocent bystander.

Multiphysics computational modelling of the cardiac ventricles

Development of cardiac multiphysics models has progressed significantly over the decades and simulations combining multiple physics interactions have become increasingly common. In this review, we summarise the progress in this field focusing on various approaches of integrating ventricular structures. electrophysiological properties, myocardial mechanics, as well as incorporating blood hemodynamics and the circulatory system. Common coupling approaches are discussed and compared, including the advantages and shortcomings of each. Currently used strategies for patient-specific implementations are highlighted and potential future improvements considered.

3D MRI of explanted sheep hearts with submillimeter isotropic spatial resolution: comparison between diffusion tensor and structure tensor imaging

Objective

The aim of the study is to compare structure tensor imaging (STI) with diffusion tensor imaging (DTI) of the sheep heart (approximately the same size as the human heart).

Materials and methods

MRI acquisition on three sheep ex vivo hearts was performed at 9.4 T/30 cm with a seven-element RF coil. 3D FLASH with an isotropic resolution of 150 µm and 3D spin-echo DTI at 600 µm were performed. Tensor analysis, angles extraction and segments divisions were performed on both volumes.

Results

A 3D FLASH allows for visualization of the detailed structure of the left and right ventricles. The helix angle determined using DTI and STI exhibited a smooth transmural change from the endocardium to the epicardium. Both the helix and transverse angles were similar between techniques. Sheetlet organization exhibited the same pattern in both acquisitions, but local angle differences were seen and identified in 17 segments representation.

Discussion

This study demonstrated the feasibility of high-resolution MRI for studying the myocyte and myolaminar architecture of sheep hearts. We presented the results of STI on three whole sheep ex vivo hearts and demonstrated a good correspondence between DTI and STI.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10334-021-00913-4.

Acute Microstructural Changes after ST-Segment Elevation Myocardial Infarction Assessed with Diffusion Tensor Imaging

Background Cardiac diffusion tensor imaging (cDTI) allows for in vivo characterization of myocardial microstructure. In cDTI, mean diffusivity and fractional anisotropy (FA)-markers of magnitude and anisotropy of diffusion of water molecules-are known to change after myocardial infarction. However, little is known about regional changes in helix angle (HA) and secondary eigenvector angle (E2A), which reflects orientations of laminar sheetlets, and their association with long-term recovery of left ventricular ejection fraction (LVEF). Purpose To assess serial changes in cDTI biomarkers in participants following ST-segment elevation myocardial infarction (STEMI) and to determine their associations with long-term left ventricular remodeling. Materials and Methods In this prospective study, 30 participants underwent cardiac MRI (3 T) after STEMI at 5 days and 3 months after reperfusion (National Institute of Health Research study no. 33963 and Research Ethics no. REC17/YH/0062). Spin-echo cDTI with second-order motion-compensation (approximate duration, 13 minutes; three sections; 18 noncollinear diffusion-weighted scans with b values of 100 sec/mm² [three acquisitions], 200 sec/mm² [three acquisitions], and 500 sec/mm² [12 acquisitions]), functional images, and late gadolinium enhancement images were obtained. Multiple regression analysis was used to assess associations between acute cDTI parameters and 3-month LVEF. Results Acutely infarcted myocardium had reduced FA, E2A, and myocytes with right-handed orientation (RHM) on HA maps compared with remote myocardium (mean remote FA = 0.36 ± 0.02 [standard deviation], mean infarcted FA = 0.25 ± 0.03, P < .001; mean remote E2A = 55° ± 9, mean infarcted E2A = 49° ± 10, P < .001; mean remote RHM = 16% ± 6, mean infarcted RHM = 9% ± 5, P < .001). All three parameters (FA, E2A, and RHM) correlated with 3-month LVEF (r = 0.68, r = 0.59, and r = 0.53, respectively), with acute FA being independently predictive of 3-month LVEF (standardized β = 0.56, P = .008) after multivariable analysis adjusting for factors, including acute LVEF and infarct size. Conclusion After ST-segment elevation myocardial infarction, diffusion becomes more isotropic in acutely infarcted myocardium as reflected by decreased fractional anisotropy. Reductions in secondary eigenvector angle suggest that the myocardial sheetlets are unable to adopt their usual steep orientations in systole, whereas reductions in myocytes with right-handed orientation on helix angle maps are likely reflective of a loss of organization among subendocardial myocytes. Correlations between these parameters and 3-month left ventricular ejection fraction highlight the potential clinical use of cardiac diffusion tensor imaging after myocardial infarction in predicting long-term remodeling. © RSNA, 2021 Online supplemental material is available for this article.

Assessing Myocardial Architecture: The Challenges and Controversies

In recent decades, investigators have strived to describe and quantify the orientation of the cardiac myocytes in an attempt to classify their arrangement in healthy and diseased hearts. There are, however, striking differences between the investigations from both a technical and methodological standpoint, thus limiting their comparability and impeding the drawing of appropriate physiological conclusions from the structural assessments. This review aims to elucidate these differences, and to propose guidance to establish methodological consensus in the field. The review outlines the theory behind myocyte orientation analysis, and importantly has identified pronounced differences in the definitions of otherwise widely accepted concepts of myocytic orientation. Based on the findings, recommendations are made for the future design of studies in the field of myocardial morphology. It is emphasised that projection of myocyte orientations, before quantification of their angulation, introduces considerable bias, and that angles should be assessed relative to the epicardial curvature. The transmural orientation of the cardiomyocytes should also not be neglected, as it is an important determinant of cardiac function. Finally, there is considerable disagreement in the literature as to how the orientation of myocardial aggregates should be assessed, but to do so in a mathematically meaningful way, the normal vector of the aggregate plane should be utilised.

An Introductory Overview of Image-Based Computational Modeling in Personalized Cardiovascular Medicine

Cardiovascular diseases account for the number one cause of deaths in the world. Part of the reason for such grim statistics is our limited understanding of the underlying mechanisms causing these devastating pathologies, which is made difficult by the invasiveness of the procedures associated with their diagnosis (e.g., inserting catheters into the coronal artery to measure blood flow to the heart). Likewise, it is also difficult to design and test assistive devices without implanting them in vivo. However, with the recent advancements made in biomedical scanning technologies and computer simulations, image-based modeling (IBM) has arisen as the next logical step in the evolution of non-invasive patient-specific cardiovascular medicine. Yet, due to its novelty, it is still relatively unknown outside of the niche field. Therefore, the goal of this manuscript is to review the current state-of-the-art and the limitations of the methods used in this area of research, as well as their applications to personalized cardiovascular investigations and treatments. Specifically, the modeling of three different physics – electrophysiology, biomechanics and hemodynamics – used in the cardiovascular IBM is discussed in the context of the physiology that each one of them describes and the mechanisms of the underlying cardiac diseases that they can provide insight into. Only the “bare-bones” of the modeling approaches are discussed in order to make this introductory material more accessible to an outside observer. Additionally, the imaging methods, the aspects of the unique cardiac anatomy derived from them, and their relation to the modeling algorithms are reviewed. Finally, conclusions are drawn about the future evolution of these methods and their potential toward revolutionizing the non-invasive diagnosis, virtual design of treatments/assistive devices, and increasing our understanding of these lethal cardiovascular diseases.

Ex vivo cardiovascular magnetic resonance diffusion weighted imaging in congenital heart disease, an insight into the microstructures of tetralogy of Fallot, biventricular and univentricular systemic right ventricle

PURPOSE

Common types of congenital heart disease exhibit a variety of structural and functional variations which may be accompanied by changes in the myocardial microstructure. We aimed to compare myocardial architecture from magnetic resonance diffusion tensor imaging (DTI) in preserved pathology specimens.

MATERIALS AND METHODS

Pathology specimens (n = 24) formalin-fixed for 40.8 ± 7.9 years comprised tetralogy of Fallot (TOF, n = 10), dextro-transposition of great arteries (D-TGA, n = 8) five with ventricular septal defect (VSD), systemic right ventricle (n = 4), situs inversus totalis (SIT, n = 1) and levo-TGA (L-TGA, n = 1). Specimens were imaged using a custom spin-echo sequence and segmented automatically according to tissue volume fraction. In each specimen T1, T2, fractional anisotropy, mean diffusivity, helix angle (HA) and sheet angle (E2A) were quantified. Pathologies were compared according to their HA gradient, HA asymmetry and E2A mean value in each myocardial segment (anterior, posterior, septal and lateral walls).

RESULTS

TOF and D-TGA with VSD had decreased helix angle gradient by - 0.34°/% and remained symmetric in the septum in comparison to D-TGA without VSD. Helix angle range was decreased by 45°. It was associated with a decreased HA gradient in the right ventricular (RV) wall, i.e. predominant circumferential myocytes. The sheet angle in the septum of TOF was opposing those of the left ventricular (LV) free wall. Univentricular systemic RV had the lowest HA gradient (- 0.43°/%) and the highest HA asymmetry (75%). HA in SIT was linear, asymmetric, and reversed with a sign change at about 70% of the depth at mid-ventricle. In L-TGA with VSD, HA was asymmetric (90%) and its gradients were decreased in the septum, anterior and lateral wall.

CONCLUSION

The organization of the myocytes as determined by DTI differs between TOF, D-TGA, L-TGA, systemic RV and SIT specimens. These differences in cardiac structure may further enlighten our understanding of cardiac function in these diverse congenital heart diseases.

Regional variations in ex-vivo diffusion tensor anisotropy are associated with cardiomyocyte remodeling in rats after left ventricular pressure overload

BACKGROUND

Pressure overload left ventricular (LV) hypertrophy is characterized by increased cardiomyocyte width and ventricle wall thickness, however the regional variation of this remodeling is unclear. Cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) may provide a non-invasive, comprehensive, and geometrically accurate method to detect regional differences in structural remodeling in hypertrophy. We hypothesized that DTI parameters, such as fractional and planar anisotropy, would reflect myocyte remodeling due to pressure overload in a regionally-dependent manner.

METHODS

We investigated the regional distributions of myocyte remodeling in rats with or without transverse aortic constriction (TAC) via direct measurement of myocyte dimensions with confocal imaging of thick tissue sections, and correlated myocyte cross-sectional area and other geometric features with parameters of diffusivity from ex-vivo DTI in the same regions of the same hearts.

RESULTS

We observed regional differences in several parameters from DTI between TAC hearts and SHAM controls. Consistent with previous studies, helix angles from DTI correlated strongly with those measured directly from histological sections (p < 0.001, R2 = 0.71). There was a transmural gradient in myocyte cross-sectional area in SHAM hearts that was diminished in the TAC group. We also found several regions of significantly altered DTI parameters in TAC LV compared to SHAM, especially in myocyte sheet angle dispersion and planar anisotropy. Among others, these parameters correlated significantly with directly measured myocyte aspect ratios.

CONCLUSIONS

These results show that structural remodeling in pressure overload LV hypertrophy is regionally heterogeneous, especially transmurally, with a greater degree of remodeling in the sub-endocardium compared to the sub-epicardium. Additionally, several parameters derived from DTI correlated significantly with measurements of myocyte geometry from direct measurement in histological sections. We suggest that DTI may provide a non-invasive, comprehensive method to detect regional structural myocyte LV remodeling during disease.

Identification of Myocardial Disarray in Patients With Hypertrophic Cardiomyopathy and Ventricular Arrhythmias

Background

Myocardial disarray is a likely focus for fatal arrhythmia in hypertrophic cardiomyopathy (HCM). This microstructural abnormality can be inferred by mapping the preferential diffusion of water along cardiac muscle fibers using diffusion tensor cardiac magnetic resonance (DT-CMR) imaging. Fractional anisotropy (FA) quantifies directionality of diffusion in 3 dimensions. The authors hypothesized that FA would be reduced in HCM due to disarray and fibrosis that may represent the anatomic substrate for ventricular arrhythmia.

Objectives

This study sought to assess FA as a noninvasive in vivo biomarker of HCM myoarchitecture and its association with ventricular arrhythmia.

Methods

A total of 50 HCM patients (47 ± 15 years of age, 77% male) and 30 healthy control subjects (46 ± 16 years of age, 70% male) underwent DT-CMR in diastole, cine, late gadolinium enhancement (LGE), and extracellular volume (ECV) imaging at 3-T.

Results

Diastolic FA was reduced in HCM compared with control subjects (0.49 ± 0.05 vs. 0.52 ± 0.03; p = 0.0005). Control subjects had a mid-wall ring of high FA. In HCM, this ring was disrupted by reduced FA, consistent with published histology demonstrating that disarray and fibrosis invade circumferentially aligned mid-wall myocytes. LGE and ECV were significant predictors of FA, in line with fibrosis contributing to low FA. Yet FA adjusted for LGE and ECV remained reduced in HCM (p = 0.028). FA in the hypertrophied segment was reduced in HCM patients with ventricular arrhythmia compared to patients without (n = 15; 0.41 ± 0.03 vs. 0.46 ± 0.06; p = 0.007). A decrease in FA of 0.05 increased odds of ventricular arrhythmia by 2.5 (95% confidence interval: 1.2 to 5.3; p = 0.015) in HCM and remained significant even after correcting for LGE, ECV, and wall thickness (p = 0.036).

Conclusions

DT-CMR assessment of left ventricular myoarchitecture matched patterns reported previously on histology. Low diastolic FA in HCM was associated with ventricular arrhythmia and is likely to represent disarray after accounting for fibrosis. The authors propose that diastolic FA could be the first in vivo marker of disarray in HCM and a potential independent risk factor.

Anatomically correct assessment of the orientation of the cardiomyocytes using diffusion tensor imaging

Diffusion tensor imaging has been used for assessing the orientation of cardiac myocytes for decades. Striking methodological differences exist between studies when quantifying these orientations. This limits the comparability between studies, and impedes collaboration and the drawing of appropriate physiological conclusions. We have sought to elucidate these differences, permitting us to propose a standardised "tool set" that might better establish consensus in future studies. We fixed hearts from seven 25 kg pigs in formalin, and scanned them using diffusion tensor imaging. Using various angle definitions as found in literature, we assessed the orientations of cardiomyocytes, comparing them in terms of helical and intrusion angles, along with the orientation of their aggregations. The difference between assessment of the helical angle with and without relation to the epicardial curvature was 25.2° (SD: 7.9) at the base, 5.8° (1.9) at the equatorial level, and 28.0° (7.0) at the apex, ANOVA P = 0.001. In comparable fashion, the intrusion angle differed by 25.9° (12.9), 7.6° (0.98) and 17.5° (4.7), P = 0.01, and the angle of the aggregates (E3-angle) differed by 25.0° (13.5) at the base, 9.4° (1.7) at the equator, and 23.1° (6.2) apically, P = 0.003. When assessing 14 definitions used in literature to calculate the orientation of aggregates, only 4 rendered identical results. The findings show that any attempt to use projection of eigenvectors introduces considerable bias. The epicardial curvature of the ventricular cone needs to be taken into account when seeking to provide accurate quantification of the orientation of the aggregated cardiomyocytes, especially in the apical and basal regions. This means that projection of eigenvectors should be avoided prior to quantifying myocyte orientation, especially when assessing radial orientation. Based on our results, we suggest appropriate methods for valid assessment of myocyte orientation using diffusion tensor imaging.

Multi-scale approaches for the simulation of cardiac electrophysiology: II - Tissue-level structure and function

Computational models of the heart, from cell-level models, through one-, two- and three-dimensional tissue-level simplifications, to biophysically-detailed three-dimensional models of the ventricles, atria or whole heart, allow the simulation of excitation and propagation of this excitation, and have provided remarkable insight into the normal and pathological functioning of the heart. In this article we present equations for modelling cellular excitation (i.e. the cell action potential) from both a phenomenological and a biophysical perspective. Hodgkin-Huxley formalism is discussed, along with the current generation of biophysically-detailed cardiac cell models. Alternative Markovian formulations for modelling ionic currents are also presented. Equations describing propagation of this cellular excitation, through one-, two- and three-dimensional idealised or realistic tissues, are then presented. For all types of model, from cell to tissue, methods for discretisation and integration of the underlying equations are discussed. The article finishes with a discussion of two tissue-level experimental imaging techniques - diffusion tensor magnetic resonance imaging and optical imaging - that can be used to provide data for parameterisation and validation of cell- and tissue-level cardiac models.

Cardiac Diffusion: Technique and Practical Applications

The 3D microarchitecture of the cardiac muscle underlies the mechanical and electrical properties of the heart. Cardiomyocytes are arranged helically through the depth of the wall, and their shortening leads to macroscopic torsion, twist, and shortening during cardiac contraction. Furthermore, cardiomyocytes are organized in sheetlets separated by shear layers, which reorientate, slip, and shear during macroscopic left ventricle (LV) wall thickening. Cardiac diffusion provides a means for noninvasive interrogation of the 3D microarchitecture of the myocardium. The fundamental principle of MR diffusion is that an MRI signal is attenuated by the self-diffusion of water in the presence of large diffusion-encoding gradients. Since water molecules are constrained by the boundaries in biological tissue (cell membranes, collagen layers, etc.), depicting their diffusion behavior elucidates the shape of the myocardial microarchitecture they are embedded in. Cardiac diffusion therefore provides a noninvasive means to understand not only the dynamic changes in cardiac microstructure of healthy myocardium during cardiac contraction but also the pathophysiological changes in the presence of disease. This unique and innovative technology offers tremendous potential to enable improved clinical diagnosis through novel microstructural and functional assessment. in vivo cardiac diffusion methods are immediately translatable to patients, opening new avenues for diagnostic investigation and treatment evaluation in a range of clinically important cardiac pathologies. This review article describes the 3D microstructure of the LV, explains in vivo and ex vivo cardiac MR diffusion acquisition and postprocessing techniques, as well as clinical applications to date. Level of Evidence: 1 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2019. J. Magn. Reson. Imaging 2020;52:348-368.

Accelerating cardiac diffusion tensor imaging combining local low-rank and 3D TV constraint

OBJECTIVE

Diffusion tensor magnetic resonance imaging (DT-MRI, or DTI) is a promising technique for invasively probing biological tissue structures. However, DTI is known to suffer from much longer acquisition time with respect to conventional MRI and the problem is worsened when dealing with in vivo acquisitions. Therefore, faster DTI for both ex vivo and in vivo scans is highly desired.

MATERIALS AND METHODS

This paper proposes a new compressed sensing (CS) reconstruction method that employs local low-rank (LLR) model and three-dimensional (3D) total variation (TV) constraint to reconstruct cardiac diffusion-weighted (DW) images from highly undersampled k-space data. The LLR model takes the set of DW images corresponding to different diffusion gradient directions as a 3D image volume and decomposes the latter into overlapping 3D blocks. Then, the 3D blocks are stacked as two-dimensional (2D) matrix. Finally, low-rank property is applied to each block matrix and the 3D TV constraint to the 3D image volume. The underlying constrained optimization problem is finally solved using the first-order fast method. The proposed method is evaluated on real ex vivo cardiac DTI data as a prerequisite to in vivo cardiac DTI applications.

RESULTS

The results on real human ex vivo cardiac DTI images demonstrate that the proposed method exhibits lower reconstruction errors for DTI indices, including fractional anisotropy (FA), mean diffusivities (MD), transverse angle (TA), and helix angle (HA), compared to existing CS-based DTI image reconstruction techniques.

CONCLUSION

The proposed method provides better reconstruction quality and more accurate DTI indices in comparison with the state-of-the-art CS-based DW image reconstruction methods.

A rule-based method for predicting the electrical activation of the heart with cardiac resynchronization therapy from non-invasive clinical data

Background

Cardiac Resynchronization Therapy (CRT) is one of the few effective treatments for heart failure patients with ventricular dyssynchrony. The pacing location of the left ventricle is indicated as a determinant of CRT outcome.

Objective

Patient specific computational models allow the activation pattern following CRT implant to be predicted and this may be used to optimize CRT lead placement.

Methods

In this study, the effects of heterogeneous cardiac substrate (scar, fast endocardial conduction, slow septal conduction, functional block) on accurately predicting the electrical activation of the LV epicardium were tested to determine the minimal detail required to create a rule based model of cardiac electrophysiology. Non-invasive clinical data (CT or CMR images and 12 lead ECG) from eighteen patients from two centers were used to investigate the models.

Results

Validation with invasive electro-anatomical mapping data identified that computer models with fast endocardial conduction were able to predict the electrical activation with a mean distance errors of 9.2 ± 0.5 mm (CMR data) or (CT data) 7.5 ± 0.7 mm.

Conclusion

This study identified a simple rule-based fast endocardial conduction model, built using non-invasive clinical data that can be used to rapidly and robustly predict the electrical activation of the heart. Pre-procedural prediction of the latest electrically activating region to identify the optimal LV pacing site could potentially be a useful clinical planning tool for CRT procedures.

Resolving the natural myocardial remodelling brought upon by cardiac contraction; a porcine ex-vivo cardiovascular magnetic resonance study of the left and right ventricle

BACKGROUND

The three-dimensional rearrangement of the right ventricular (RV) myocardium during cardiac deformation is unknown. Previous in-vivo studies have shown that myocardial left ventricular (LV) deformation is driven by rearrangement of aggregations of cardiomyocytes that can be characterised by changes in the so-called E3-angle. Ex-vivo imaging offers superior spatial resolution compared with in-vivo measurements, and can thus provide novel insight into the deformation of the myocardial microstructure in both ventricles. This study sought to describe the dynamic changes of the orientations of the cardiomyocytes in both ventricles brought upon by cardiac contraction, with particular interest in the thin-walled RV, which has not previously been described in terms of its micro-architecture.

METHODS

The hearts of 14 healthy 20 kg swine were excised and preserved in either a relaxed state or a contracted state. Myocardial architecture was assessed and compared between the two contractional states by quantification of the helical, transmural and E3-angles of the cardiomyocytes using high-resolution diffusion tensor imaging.

RESULTS

The differences between the two states of contraction were most pronounced in the endocardium where the E3-angle decreased from 78.6° to 24.8° in the LV and from 82.6° to 68.6° in the RV. No significant change in neither the helical nor the transmural angle was found in the cardiomyocytes of the RV. In the endocardium of the LV, however, the helical angle increased from 35.4° to 47.8° and the transmural angle increased from 3.1° to 10.4°.

CONCLUSION

The entire myocardium rearranges through the cardiac cycle with the change in the orientation of the aggregations of cardiomyocytes being the predominant mediator of myocardial wall thickening. Interestingly, differences also exist between the RV and LV, which helps in the explanation of the different physiological capabilities of the ventricles.

Diffusion tensor cardiovascular magnetic resonance

Cardiac structure and function are complex and inter-related. Current in vivo techniques assess the heart on a macroscopic scale, but a novel technique called diffusion tensor cardiovascular magnetic resonance (DT-CMR) can now assess the cardiac microstructure non-invasively. It provides information on the helical arrangement of cardiomyocytes that drives torsion and offers dynamic assessment of the sheetlets (aggregated cardiomyocytes) that rotate through the cardiac cycle to facilitate wall thickening. Through diffusion biomarkers, the expansion and organisation of the underlying myocardium can be described. DT-CMR has already identified novel microstructural abnormalities in cardiomyopathy, and ischaemic and congenital heart disease. This new knowledge supports the potential of DT-CMR to improve diagnostics and prognostication in various cardiac diseases.

A framework for biomechanics simulations using four-chamber cardiac models

Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly-fewer elements are needed to build the cardiac geometry-compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.

Personalized Cardiac Computational Models: From Clinical Data to Simulation of Infarct-Related Ventricular Tachycardia

In the chronic stage of myocardial infarction, a significant number of patients develop life-threatening ventricular tachycardias (VT) due to the arrhythmogenic nature of the remodeled myocardium. Radiofrequency ablation (RFA) is a common procedure to isolate reentry pathways across the infarct scar that are responsible for VT. Unfortunately, this strategy show relatively low success rates; up to 50% of patients experience recurrent VT after the procedure. In the last decade, intensive research in the field of computational cardiac electrophysiology (EP) has demonstrated the ability of three-dimensional (3D) cardiac computational models to perform in-silico EP studies. However, the personalization and modeling of certain key components remain challenging, particularly in the case of the infarct border zone (BZ). In this study, we used a clinical dataset from a patient with a history of infarct-related VT to build an image-based 3D ventricular model aimed at computational simulation of cardiac EP, including detailed patient-specific cardiac anatomy and infarct scar geometry. We modeled the BZ in eight different ways by combining the presence or absence of electrical remodeling with four different levels of image-based patchy fibrosis (0, 10, 20, and 30%). A 3D torso model was also constructed to compute the ECG. Patient-specific sinus activation patterns were simulated and validated against the patient's ECG. Subsequently, the pacing protocol used to induce reentrant VTs in the EP laboratory was reproduced in-silico. The clinical VT was induced with different versions of the model and from different pacing points, thus identifying the slow conducting channel responsible for such VT. Finally, the real patient's ECG recorded during VT episodes was used to validate our simulation results and to assess different strategies to model the BZ. Our study showed that reduced conduction velocities and heterogeneity in action potential duration in the BZ are the main factors in promoting reentrant activity. Either electrical remodeling or fibrosis in a degree of at least 30% in the BZ were required to initiate VT. Moreover, this proof-of-concept study confirms the feasibility of developing 3D computational models for cardiac EP able to reproduce cardiac activation in sinus rhythm and during VT, using exclusively non-invasive clinical data.

Postmortem cardiac magnetic resonance in sudden cardiac death

Postmortem imaging is increasingly used in forensic practice as good complementary tool to conventional autopsy investigations. Over the last decade, postmortem cardiac magnetic resonance (PMCMR) imaging was introduced in forensic investigations of natural deaths related to cardiovascular diseases, which represent the most common causes of death in developed countries. Postmortem CMR application has yielded interesting results in ischemic myocardium injury investigations and in visualizing other pathological findings in the heart. This review presents the actual state of postmortem imaging for cardiovascular pathologies in cases of sudden cardiac death (SCD), taking into consideration both the advantages and limitations of PMCMR application.

Effect of intra-myocardial Algisyl-LVR™ injectates on fibre structure in porcine heart failure

Recent preclinical trials have shown that alginate injections are a promising treatment for ischemic heart disease. Although improvements in heart function and global structure have been reported following alginate implants, the underlying structure is poorly understood. Using high resolution ex vivo MRI and DT-MRI of the hearts of normal control swine (n = 8), swine with induced heart failure (n = 5), and swine with heart failure and alginate injection treatment (n = 6), we visualized and quantified the fibre distribution and implant material geometry. Our findings show that the alginate injectates form solid ellipsoids with a retention rate of 68.7% ± 21.3% (mean ± SD) and a sphericity index of 0.37 ± 0.03. These ellipsoidal shapes solidified predominantly at the mid-wall position with an inclination of -4.9° ± 31.4° relative to the local circumferential direction. Overall, the change to left ventricular wall thickness and myofiber orientation was minor and was associated with heart failure and not the presence of injectates. These results show that alginate injectates conform to the pre-existing tissue structure, likely expanding along directions of least resistance as mass is added to the injection sites. The alginate displaces the myocardial tissue predominantly in the longitudinal direction, causing minimal disruption to the surrounding myofiber orientations.

Will the real ventricular architecture please stand up?

Ventricular twisting, essential for cardiac function, is attributed to the contraction of myocardial helical fibers. The exact relationship between ventricular anatomy and function remains to be determined, but one commonly used explanatory model is the helical ventricular myocardial band (HVMB) model of Torrent‐Guasp. This model has been successful in explaining many aspects of ventricular function, (Torrent‐Guasp et al. Eur. J. Cardiothorac. Surg., 25, 376, 2004; Buckberg et al. Eur. J. Cardiothorac. Surg., 47, 587, 2015; Buckberg et al. Eur. J. Cardiothorac. Surg. 47, 778, 2015) but the model ignores important aspects of ventricular anatomy and should probably be replaced. The purpose of this review is to compare the HVMB model with a different model (nested layers). A complication when interpreting experimental observations that relate anatomy to function is that, in the myocardium, shortening does not always imply activation and lengthening does not always imply inactivation.

Computational Modeling for Cardiac Resynchronization Therapy

Cardiac resynchronization therapy (CRT) is an effective treatment for heart failure (HF) patients with an electrical substrate pathology causing ventricular dyssynchrony. However 40-50% of patients do not respond to treatment. Cardiac modeling of the electrophysiology, electromechanics, and hemodynamics of the heart has been used to study mechanisms behind HF pathology and CRT response. Recently, multi-scale dyssynchronous HF models have been used to study optimal device settings and optimal lead locations, investigate the underlying cardiac pathophysiology, as well as investigate emerging technologies proposed to treat cardiac dyssynchrony. However the breadth of patient and experimental data required to create and parameterize these models and the computational resources required currently limits the use of these models to small patient numbers. In the future, once these technical challenges are overcome, biophysically based models of the heart have the potential to become a clinical tool to aid in the diagnosis and treatment of HF.

Right Ventricular Fiber Structure as a Compensatory Mechanism in Pressure Overload: A Computational Study

Right ventricular failure (RVF) is a lethal condition in diverse pathologies. Pressure overload is the most common etiology of RVF, but our understanding of the tissue structure remodeling and other biomechanical factors involved in RVF is limited. Some remodeling patterns are interpreted as compensatory mechanisms including myocyte hypertrophy, extracellular fibrosis, and changes in fiber orientation. However, the specific implications of these changes, especially in relation to clinically observable measurements, are difficult to investigate experimentally. In this computational study, we hypothesized that, with other variables constant, fiber orientation alteration provides a quantifiable and distinct compensatory mechanism during RV pressure overload (RVPO). Numerical models were constructed using a rabbit model of chronic pressure overload RVF based on intraventricular pressure measurements, CINE magnetic resonance imaging (MRI), and diffusion tensor MRI (DT-MRI). Biventricular simulations were conducted under normotensive and hypertensive boundary conditions using variations in RV wall thickness, tissue stiffness, and fiber orientation to investigate their effect on RV pump function. Our results show that a longitudinally aligned myocardial fiber orientation contributed to an increase in RV ejection fraction (RVEF). This effect was more pronounced in response to pressure overload. Likewise, models with longitudinally aligned fiber orientation required a lesser contractility for maintaining a target RVEF against elevated pressures. In addition to increased wall thickness and material stiffness (diastolic compensation), systolic mechanisms in the forms of myocardial fiber realignment and changes in contractility are likely involved in the overall compensatory responses to pressure overload.

Changes in overall ventricular myocardial architecture in the setting of a porcine animal model of right ventricular dilation

BACKGROUND

Chronic pulmonary regurgitation often leads to myocardial dysfunction and heart failure. It is not fully known why secondary hypertrophy cannot fully protect against the increase in wall stress brought about by the increased end-diastolic volume in ventricular dilation. It has been assumed that mural architecture is not deranged in this situation, but we hypothesised that there might be a change in the pattern of orientation of the aggregations of cardiomyocytes, which would contribute to contractile impairment.

METHODS

We created pulmonary valvular regurgitation by open chest, surgical suturing of its leaflets in seven piglets, performing sham operations in seven control animals. Using cardiovascular magnetic resonance imaging after 12 weeks of recovery, we demonstrated significantly increased right ventricular volumes in the test group. After sacrifice, diffusion tensor imaging of their hearts permitted measurement of the orientation of the cardiomyocytes.

RESULTS

The helical angles in the right ventricle approached a more circumferential orientation in the setting of right ventricular RV dilation (p = 0.007), with an increased proportion of surface-parallel cardiomyocytes. In contrast, this proportion decreased in the left ventricle. Also in the left ventricle a higher proportion of E3 angles with a value around zero was found, and conversely a lower proportion of angles was found with a numerical higher value. In the dilated right ventricle the proportion of E3 angles around -90° is increased, while the proportion around 90° is decreased.

CONCLUSION

Contrary to traditional views, there is a change in the orientation of both the left ventricular and right ventricular cardiomyocytes subsequent to right ventricular dilation. This will change their direction of contraction and hinder the achievement of normalisation of cardiomyocytic strain, affecting overall contractility. We suggest that the aetiology of the cardiac failure induced by right vetricular dilation may be partly explained by morphological changes in the myocardium itself.

Ventricular myocardium development and the role of connexins in the human fetal heart

The developmental timeline of the human heart remains elusive. The heart takes on its characteristic four chambered appearance by ~56 days gestational age (DGA). However, owing to the complexities (both technical and logistical) of exploring development in utero, we understand little of how the ventricular walls develop. To address this, we employed diffusion tensor magnetic resonance imaging to explore the architecture and tissue organization of the developing heart aged 95-143 DGA. We show that fractional anisotropy increases (from ~0.1 to ~0.5), diffusion coefficients decrease (from ~1 × 10-3mm2/sec to ~0.4 × 10-3mm2/sec), and fiber paths, extracted by tractography, increase linearly with gestation, indicative of the increasing organization of the ventricular myocytes. By 143 DGA, the developing heart has the classical helical organization observed in mature mammalian tissue. This was accompanied by an increase in connexin 43 and connexin 40 expression levels, suggesting their role in the development of the ventricular conduction system and that electrical propagation across the heart is facilitated in later gestation. Our findings highlight a key developmental window for the structural organization of the fetal heart.

Evaluation of the impact of strain correction on the orientation of cardiac diffusion tensors with in vivo and ex vivo porcine hearts

PURPOSE

To evaluate the importance of strain-correcting stimulated echo acquisition mode echo-planar imaging cardiac diffusion tensor imaging.

METHODS

Healthy pigs (n = 11) were successfully scanned with a 3D cine displacement-encoded imaging with stimulated echoes and a monopolar-stimulated echo-planar imaging diffusion tensor imaging sequence at 3 T during diastasis, peak systole, and strain sweet spots in a midventricular short-axis slice. The same diffusion tensor imaging sequence was repeated ex vivo after arresting the hearts in either a relaxed (KCl-induced) or contracted (BaCl₂ -induced) state. The displacement-encoded imaging with stimulated echoes data were used to strain-correct the in vivo cardiac diffusion tensor imaging in diastole and systole. The orientation of the primary (helix angles) and secondary (E2A) diffusion eigenvectors was compared with and without strain correction and to the strain-free ex vivo data.

RESULTS

Strain correction reduces systolic E2A significantly when compared without strain correction and ex vivo (median absolute E2A = 34.3° versus E2A = 57.1° (P = 0.01), E2A = 60.5° (P = 0.006), respectively). The systolic distribution of E2A without strain correction is closer to the contracted ex vivo distribution than with strain correction, root mean square deviation of 0.027 versus 0.038.

CONCLUSIONS

The current strain-correction model amplifies the contribution of microscopic strain to diffusion resulting in an overcorrection of E2A. Results show that a new model that considers cellular rearrangement is required. Magn Reson Med 79:2205-2215, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Validation of diffusion tensor MRI measurements of cardiac microstructure with structure tensor synchrotron radiation imaging

BACKGROUND

Diffusion tensor imaging (DTI) is widely used to assess tissue microstructure non-invasively. Cardiac DTI enables inference of cell and sheetlet orientations, which are altered under pathological conditions. However, DTI is affected by many factors, therefore robust validation is critical. Existing histological validation is intrinsically flawed, since it requires further tissue processing leading to sample distortion, is routinely limited in field-of-view and requires reconstruction of three-dimensional volumes from two-dimensional images. In contrast, synchrotron radiation imaging (SRI) data enables imaging of the heart in 3D without further preparation following DTI. The objective of the study was to validate DTI measurements based on structure tensor analysis of SRI data.

METHODS

One isolated, fixed rat heart was imaged ex vivo with DTI and X-ray phase contrast SRI, and reconstructed at 100 μm and 3.6 μm isotropic resolution respectively. Structure tensors were determined from the SRI data and registered to the DTI data.

RESULTS

Excellent agreement in helix angles (HA) and transverse angles (TA) was observed between the DTI and structure tensor synchrotron radiation imaging (STSRI) data, where HADTI-STSRI = -1.4° ± 23.2° and TADTI-STSRI = -1.4° ± 35.0° (mean ± 1.96 standard deviation across all voxels in the left ventricle). STSRI confirmed that the primary eigenvector of the diffusion tensor corresponds with the cardiomyocyte long-axis across the whole myocardium.

CONCLUSIONS

We have used STSRI as a novel and high-resolution gold standard for the validation of DTI, allowing like-with-like comparison of three-dimensional tissue structures in the same intact heart free of distortion. This represents a critical step forward in independently verifying the structural basis and informing the interpretation of cardiac DTI data, thereby supporting the further development and adoption of DTI in structure-based electro-mechanical modelling and routine clinical applications.

Diffusion MRI in the heart

Diffusion MRI provides unique information on the structure, organization, and integrity of the myocardium without the need for exogenous contrast agents. Diffusion MRI in the heart, however, has proven technically challenging because of the intrinsic non-rigid deformation during the cardiac cycle, displacement of the myocardium due to respiratory motion, signal inhomogeneity within the thorax, and short transverse relaxation times. Recently developed accelerated diffusion-weighted MR acquisition sequences combined with advanced post-processing techniques have improved the accuracy and efficiency of diffusion MRI in the myocardium. In this review, we describe the solutions and approaches that have been developed to enable diffusion MRI of the heart in vivo, including a dual-gated stimulated echo approach, a velocity- (M1 ) or an acceleration- (M2 ) compensated pulsed gradient spin echo approach, and the use of principal component analysis filtering. The structure of the myocardium and the application of these techniques in ischemic heart disease are also briefly reviewed. The advent of clinical MR systems with stronger gradients will likely facilitate the translation of cardiac diffusion MRI into clinical use. The addition of diffusion MRI to the well-established set of cardiovascular imaging techniques should lead to new and complementary approaches for the diagnosis and evaluation of patients with heart disease. © 2015 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Submillimeter diffusion tensor imaging and late gadolinium enhancement cardiovascular magnetic resonance of chronic myocardial infarction

BACKGROUND

Knowledge of the three-dimensional (3D) infarct structure and fiber orientation remodeling is essential for complete understanding of infarct pathophysiology and post-infarction electromechanical functioning of the heart. Accurate imaging of infarct microstructure necessitates imaging techniques that produce high image spatial resolution and high signal-to-noise ratio (SNR). The aim of this study is to provide detailed reconstruction of 3D chronic infarcts in order to characterize the infarct microstructural remodeling in porcine and human hearts.

METHODS

We employed a customized diffusion tensor imaging (DTI) technique in conjunction with late gadolinium enhancement (LGE) cardiovascular magnetic resonance (CMR) on a 3T clinical scanner to image, at submillimeter resolution, myofiber orientation and scar structure in eight chronically infarcted porcine hearts ex vivo. Systematic quantification of local microstructure was performed and the chronic infarct remodeling was characterized at different levels of wall thickness and scar transmurality. Further, a human heart with myocardial infarction was imaged using the same DTI sequence.

RESULTS

The SNR of non-diffusion-weighted images was >100 in the infarcted and control hearts. Mean diffusivity and fractional anisotropy (FA) demonstrated a 43% increase, and a 35% decrease respectively, inside the scar tissue. Despite this, the majority of the scar showed anisotropic structure with FA higher than an isotropic liquid. The analysis revealed that the primary eigenvector orientation at the infarcted wall on average followed the pattern of original fiber orientation (imbrication angle mean: 1.96 ± 11.03° vs. 0.84 ± 1.47°, p = 0.61, and inclination angle range: 111.0 ± 10.7° vs. 112.5 ± 6.8°, p = 0.61, infarcted/control wall), but at a higher transmural gradient of inclination angle that increased with scar transmurality (r = 0.36) and the inverse of wall thickness (r = 0.59). Further, the infarcted wall exhibited a significant increase in both the proportion of left-handed epicardial eigenvectors, and in the angle incoherency. The infarcted human heart demonstrated preservation of primary eigenvector orientation at the thinned region of infarct, consistent with the findings in the porcine hearts.

CONCLUSIONS

The application of high-resolution DTI and LGE-CMR revealed the detailed organization of anisotropic infarct structure at a chronic state. This information enhances our understanding of chronic post-infarction remodeling in large animal and human hearts.

Transmural gradients of myocardial structure and mechanics: Implications for fiber stress and strain in pressure overload

Although a truly complete understanding of whole heart activation, contraction, and deformation is well beyond our current reach, a significant amount of effort has been devoted to discovering and understanding the mechanisms by which myocardial structure determines cardiac function to better treat patients with cardiac disease. Several experimental studies have shown that transmural fiber strain is relatively uniform in both diastole and systole, in contrast to predictions from traditional mechanical theory. Similarly, mathematical models have largely predicted uniform fiber stress across the wall. The development of this uniform pattern of fiber stress and strain during filling and ejection is due to heterogeneous transmural distributions of several myocardial structures. This review summarizes these transmural gradients, their contributions to fiber mechanics, and the potential functional effects of their remodeling during pressure overload hypertrophy.

Anomalous Diffusion in Cardiac Tissue as an Index of Myocardial Microstructure

Diffusion in biological tissues is known to be hindered by the structural complexity of the underlying medium. In the heart, improved characterisation on how this complexity influences acquired diffusion weighted signals is key to advancing our interpretation of diffusion magnetic resonance imaging, as well as to propose novel biomarkers to further characterise myocardial microstructure. In this work, we propose stretched Mittag-Leffler signal decay models for the quantification of the anomalous decay observed in acquired diffusion weighted signals. Our results, analysed in ex vivo healthy, fixed rat ventricles, indicate that such a representation suffices to capture the anomalous signal decay observed in the myocardial syncytium. The subdiffusive order of signal decay is shown to encode independent information to that encapsulated by standard diffusion tensor metrics, and thus may provide additional information on tissue microstructure. Moreover, subdiffusion gradients are shown to be indicative of the total structural heterogeneity spanning the left ventricular wall, which includes progressive myolaminae branching and spatially varying densities of perimysial collagen, microvasculature and adipose tissue. The proposed approach may therefore have important implications for the characterisation of tissue microstructure, both in cardiac and other tissue types.

Evaluation of the Differences of Myocardial Fibers between Acute and Chronic Myocardial Infarction: Application of Diffusion Tensor Magnetic Resonance Imaging in a Rhesus Monkey Model

Objective

To understand microstructural changes after myocardial infarction (MI), we evaluated myocardial fibers of rhesus monkeys during acute or chronic MI, and identified the differences of myocardial fibers between acute and chronic MI.

Materials and Methods

Six fixed hearts of rhesus monkeys with left anterior descending coronary artery ligation for 1 hour or 84 days were scanned by diffusion tensor magnetic resonance imaging (MRI) to measure apparent diffusion coefficient (ADC), fractional anisotropy (FA) and helix angle (HA).

Results

Comparing with acute MI monkeys (FA: 0.59 ± 0.02; ADC: 5.0 ± 0.6 × 10^-4 mm²/s; HA: 94.5 ± 4.4°), chronic MI monkeys showed remarkably decreased FA value (0.26 ± 0.03), increased ADC value (7.8 ± 0.8 × 10^-4mm²/s), decreased HA transmural range (49.5 ± 4.6°) and serious defects on endocardium in infarcted regions. The HA in infarcted regions shifted to more components of negative left-handed helix in chronic MI monkeys (-38.3 ± 5.0°–11.2 ± 4.3°) than in acute MI monkeys (-41.4 ± 5.1°–53.1 ± 3.7°), but the HA in remote regions shifted to more components of positive right-handed helix in chronic MI monkeys (-43.8 ± 2.7°–66.5 ± 4.9°) than in acute MI monkeys (-59.5 ± 3.4°–64.9 ± 4.3°).

Conclusion

Diffusion tensor MRI method helps to quantify differences of mechanical microstructure and water diffusion of myocardial fibers between acute and chronic MI monkey's models.

Mapping cardiac microstructure of rabbit heart in different mechanical states by high resolution diffusion tensor imaging: A proof-of-principle study

Myocardial microstructure and its macroscopic materialisation are fundamental to the function of the heart. Despite this importance, characterisation of cellular features at the organ level remains challenging, and a unifying description of the structure of the heart is still outstanding. Here, we optimised diffusion tensor imaging data to acquire high quality data in ex vivo rabbit hearts in slack and contractured states, approximating diastolic and systolic conditions. The data were analysed with a suite of methods that focused on different aspects of the myocardium. In the slack heart, we observed a similar transmural gradient in helix angle of the primary eigenvector of up to 23.6°/mm in the left ventricle and 24.2°/mm in the right ventricle. In the contractured heart, the same transmural gradient remained largely linear, but was offset by up to +49.9° in the left ventricle. In the right ventricle, there was an increase in the transmural gradient to 31.2°/mm and an offset of up to +39.0°. The application of tractography based on each eigenvector enabled visualisation of streamlines that depict cardiomyocyte and sheetlet organisation over large distances. We observed multiple V- and N-shaped sheetlet arrangements throughout the myocardium, and insertion of sheetlets at the intersection of the left and right ventricle. This study integrates several complementary techniques to visualise and quantify the heart's microstructure, projecting parameter representations across different length scales. This represents a step towards a more comprehensive characterisation of myocardial microstructure at the whole organ level.

Parametric Modeling of the Mouse Left Ventricular Myocardial Fiber Structure

Magnetic resonance diffusion tensor imaging (DTI) has greatly facilitated detailed quantifications of myocardial structures. However, structural patterns, such as the distinctive transmural rotation of the fibers, remain incompletely described. To investigate the validity and practicality of pattern-based analysis, 3D DTI was performed on 13 fixed mouse hearts and fiber angles in the left ventricle were transformed and fitted to parametric expressions constructed from elementary functions of the prolate spheroidal spatial variables. It was found that, on average, the myocardial fiber helix angle could be represented to 6.5° accuracy by the equivalence of a product of 10th-order polynomials of the radial and longitudinal variables, and 17th-order Fourier series of the circumferential variable. Similarly, the fiber imbrication angle could be described by 10th-order polynomials and 24th-order Fourier series, to 5.6° accuracy. The representations, while relatively concise, did not adversely affect the information commonly derived from DTI datasets including the whole-ventricle mean fiber helix angle transmural span and atlases constructed for the group. The unique ability of parametric models for predicting the 3D myocardial fiber structure from finite number of 2D slices was also demonstrated. These findings strongly support the principle of parametric modeling for characterizing myocardial structures in the mouse and beyond.

Diffusion Tensor Imaging of Healthy and Infarcted Porcine Hearts: Study on the Impact of Formalin Fixation

BACKGROUND

Due to complexities of in-vivo cardiac diffusion tensor imaging (DTI), ex-vivo formalin-fixed specimens are used to investigate cardiac remodeling in diseases, and reported results have shown conflicting trends. This study investigates the impact of formalin-fixation on diffusion properties and optimizes tracking parameters based on controls to understand remodeling in myocardial-infarction (MI).

METHODS

DTI was performed on 4 healthy (controls) and 4 MI induced formalin-fixed (PoMI) ex-vivo porcine hearts. Controls were scanned pre-fixation (PrCtrl) and re-scanned (PoCtrl) after formalin-fixation. Fractional anisotropy (FA) and apparent diffusion coefficient (ADC) were estimated in all hearts. Tracking parameters (FA, tract termination angle (TTA), fiber-length) were optimized in controls and then used to investigate structural remodeling in PoMI hearts.

RESULTS

Fixation increased ADC and decreased FA. PoMI showed increased ADC but decreased FA in infarcted zone compared to remote zone. TTA showed sharp increase in slope from 5°-10°, which flattened after 25° in all groups. Mean fiber-length for different tracking length range showed that PoCtrl had shorter fibers compared to PrCtrl. Fibers around infarction were shorter in length and disarrayed compared to PoCtrl group.

CONCLUSION

Formalin-fixation affects diffusion properties and hence DTI parametric trends observed in pathology may be influenced by the fixation process which can cause contradictory findings.

Biomechanics Simulations Using Cubic Hermite Meshes with Extraordinary Nodes for Isogeometric Cardiac Modeling

Cubic Hermite hexahedral finite element meshes have some well-known advantages over linear tetrahedral finite element meshes in biomechanical and anatomic modeling using isogeometric analysis. These include faster convergence rates as well as the ability to easily model rule-based anatomic features such as cardiac fiber directions. However, it is not possible to create closed complex objects with only regular nodes; these objects require the presence of extraordinary nodes (nodes with 3 or >= 5 adjacent elements in 2D) in the mesh. The presence of extraordinary nodes requires new constraints on the derivatives of adjacent elements to maintain continuity. We have developed a new method that uses an ensemble coordinate frame at the nodes and a local-to-global mapping to maintain continuity. In this paper, we make use of this mapping to create cubic Hermite models of the human ventricles and a four-chamber heart. We also extend the methods to the finite element equations to perform biomechanics simulations using these meshes. The new methods are validated using simple test models and applied to anatomically accurate ventricular meshes with valve annuli to simulate complete cardiac cycle simulations.