Comparison of diffusion tensor imaging by cardiovascular magnetic resonance and gadolinium enhanced 3D image intensity approaches to investigation of structural anisotropy in explanted rat hearts

Monte Carlo simulation of water diffusion through cardiac tissue models

Monte Carlo diffusion simulations are commonly used to establish a reliable ground truth of tissue microstructure, including for the validation of diffusion-weighted MRI. However, selecting simulation parameters is challenging and affects validity and reproducibility. We conducted experiments to investigate critical conditions in Monte Carlo simulations, such as tissue representation complexity, simulated molecules, update duration, and compartment size. Results show significant changes in microstructure characteristics when parameters are altered, emphasizing the importance of careful control for a reliable ground truth.

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.

Quantification of Myocyte Disarray in Human Cardiac Tissue

Proper three-dimensional (3D)-cardiomyocyte orientation is important for an effective tension production in cardiac muscle. Cardiac diseases can cause severe remodeling processes in the heart, such as cellular misalignment, that can affect both the electrical and mechanical functions of the organ. To date, a proven methodology to map and quantify myocytes disarray in massive samples is missing. In this study, we present an experimental pipeline to reconstruct and analyze the 3D cardiomyocyte architecture in massive samples. We employed tissue clearing, staining, and advanced microscopy techniques to detect sarcomeres in relatively large human myocardial strips with micrometric resolution. Z-bands periodicity was exploited in a frequency analysis approach to extract the 3D myofilament orientation, providing an orientation map used to characterize the tissue organization at different spatial scales. As a proof-of-principle, we applied the proposed method to healthy and pathologically remodeled human cardiac tissue strips. Preliminary results suggest the reliability of the method: strips from a healthy donor are characterized by a well-organized tissue, where the local disarray is log-normally distributed and slightly depends on the spatial scale of analysis; on the contrary, pathological strips show pronounced tissue disorganization, characterized by local disarray significantly dependent on the spatial scale of analysis. A virtual sample generator is developed to link this multi-scale disarray analysis with the underlying cellular architecture. This approach allowed us to quantitatively assess tissue organization in terms of 3D myocyte angular dispersion and may pave the way for developing novel predictive models based on structural data at cellular resolution.

3D magnetization transfer (MT) for the visualization of cardiac free-running Purkinje fibers: an ex vivo proof of concept

OBJECTIVES

We investigate the possibility to exploit high-field MRI to acquire 3D images of Purkinje network which plays a crucial role in cardiac function. Since Purkinje fibers (PF) have a distinct cellular structure and are surrounded by connective tissue, we investigated conventional contrast mechanisms along with the magnetization transfer (MT) imaging technique to improve image contrast between ventricular structures of differing macromolecular content.

METHODS

Three fixed porcine ventricular samples were used with free-running PFs on the endocardium. T1, T2*, T2, and M0 were evaluated on 2D slices for each sample at 9.4 T. MT parameters were optimized using hard pulses with different amplitudes, offset frequencies and durations. The cardiac structure was assessed through 2D and 3D T1w images with isotropic resolutions of 150 µm. Histology, immunofluorescence, and qPCR were performed to analyze collagen contents of cardiac tissue and PF.

RESULTS

An MT preparation module of 350 ms duration inserted into the sequence with a B1 = 10 µT and frequency offset = 3000 Hz showed the best contrast, approximately 0.4 between PFs and myocardium. Magnetization transfer ratio (MTR) appeared higher in the cardiac tissue (MTR = 44.7 ± 3.5%) than in the PFs (MTR = 25.2 ± 6.3%).

DISCUSSION

MT significantly improves contrast between PFs and ventricular myocardium and appears promising for imaging the 3D architecture of the Purkinje network.

Comprehensive assessment of myocardial remodeling in ischemic heart disease by synchrotron propagation based X-ray phase contrast imaging

Cardiovascular research is in an ongoing quest for a superior imaging method to integrate gross-anatomical information with microanatomy, combined with quantifiable parameters of cardiac structure. In recent years, synchrotron radiation-based X-ray Phase Contrast Imaging (X-PCI) has been extensively used to characterize soft tissue in detail. The objective was to use X-PCI to comprehensively quantify ischemic remodeling of different myocardial structures, from cell to organ level, in a rat model of myocardial infarction. Myocardial infarction-induced remodeling was recreated in a well-established rodent model. Ex vivo rodent hearts were imaged by propagation based X-PCI using two configurations resulting in 5.8 µm and 0.65 µm effective pixel size images. The acquired datasets were used for a comprehensive assessment of macrostructural changes including the whole heart and vascular tree morphology, and quantification of left ventricular myocardial thickness, mass, volume, and organization. On the meso-scale, tissue characteristics were explored and compared with histopathological methods, while microstructural changes were quantified by segmentation of cardiomyocytes and calculation of cross-sectional areas. Propagation based X-PCI provides detailed visualization and quantification of morphological changes on whole organ, tissue, vascular as well as individual cellular level of the ex vivo heart, with a single, non-destructive 3D imaging modality.

A 3D high resolution MRI method for the visualization of cardiac fibro-fatty infiltrations

Modifications of the myocardial architecture can cause abnormal electrical activity of the heart. Fibro-fatty infiltrations have been implicated in various cardiac pathologies associated with arrhythmias and sudden cardiac death, such as arrhythmogenic right ventricular cardiomyopathy (ARVC). Here, we report the development of an MRI protocol to observe these modifications at 9.4 T. Two fixed ex vivo human hearts, one healthy and one ARVC, were imaged with an Iterative decomposition with echo asymmetry and least-square estimations (IDEAL) and a magnetization transfer (MT) 3D sequences. The resulting fat fraction and MT ratio (MTR) were analyzed and compared to histological analysis of the three regions (“ARVC triangle”) primarily involved in ARVC structural remodeling. In the ARVC heart, high fat content was observed in the “ARVC triangle” and the superimposition of the MTR and fat fraction allowed the identification of fibrotic regions in areas without the presence of fat. The healthy heart exhibited twice less fat than the ARVC heart (31.9%, 28.7% and 1.3% of fat in the same regions, respectively). Localization of fat and fibrosis were confirmed by means of histology. This non-destructive approach allows the investigation of structural remodeling in human pathologies where fibrosis and/or fatty tissue infiltrations are expected to occur.

A Review of Ex Vivo X-ray Microfocus Computed Tomography-Based Characterization of the Cardiovascular System

Cardiovascular malformations and diseases are common but complex and often not yet fully understood. To better understand the effects of structural and microstructural changes of the heart and the vasculature on their proper functioning, a detailed characterization of the microstructure is crucial. In vivo imaging approaches are noninvasive and allow visualizing the heart and the vasculature in 3D. However, their spatial image resolution is often too limited for microstructural analyses, and hence, ex vivo imaging is preferred for this purpose. Ex vivo X-ray microfocus computed tomography (microCT) is a rapidly emerging high-resolution 3D structural imaging technique often used for the assessment of calcified tissues. Contrast-enhanced microCT (CE-CT) or phase-contrast microCT (PC-CT) improve this technique by additionally allowing the distinction of different low X-ray-absorbing soft tissues. In this review, we present the strengths of ex vivo microCT, CE-CT and PC-CT for quantitative 3D imaging of the structure and/or microstructure of the heart, the vasculature and their substructures in healthy and diseased state. We also discuss their current limitations, mainly with regard to the contrasting methods and the tissue preparation.

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.

Cardiac multi-scale investigation of the right and left ventricle ex vivo: a review

The heart is a complex multi-scale system composed of components integrated at the subcellular, cellular, tissue and organ levels. The myocytes, the contractile elements of the heart, form a complex three-dimensional (3D) network which enables propagation of the electrical signal that triggers the contraction to efficiently pump blood towards the whole body. Cardiovascular diseases (CVDs), a major cause of mortality in developed countries, often lead to cardiovascular remodeling affecting cardiac structure and function at all scales, from myocytes and their surrounding collagen matrix to the 3D organization of the whole heart. As yet, there is no consensus as to how the myocytes are arranged and packed within their connective tissue matrix, nor how best to image them at multiple scales. Cardiovascular imaging is routinely used to investigate cardiac structure and function as well as for the evaluation of cardiac remodeling in CVDs. For a complete understanding of the relationship between structural remodeling and cardiac dysfunction in CVDs, multi-scale imaging approaches are necessary to achieve a detailed description of ventricular architecture along with cardiac function. In this context, ventricular architecture has been extensively studied using a wide variety of imaging techniques: ultrasound (US), optical coherence tomography (OCT), microscopy (confocal, episcopic, light sheet, polarized light), magnetic resonance imaging (MRI), micro-computed tomography (micro-CT) and, more recently, synchrotron X-ray phase contrast imaging (SR X-PCI). Each of these techniques have their own set of strengths and weaknesses, relating to sample size, preparation, resolution, 2D/3D capabilities, use of contrast agents and possibility of performing together with in vivo studies. Therefore, the combination of different imaging techniques to investigate the same sample, thus taking advantage of the strengths of each method, could help us to extract the maximum information about ventricular architecture and function. In this review, we provide an overview of available and emerging cardiovascular imaging techniques for assessing myocardial architecture ex vivo and discuss their utility in being able to quantify cardiac remodeling, in CVDs, from myocyte to whole organ.

Short microsecond pulses achieve homogeneous electroporation of elongated biological cells irrespective of their orientation in electric field

In gene electrotransfer and cardiac ablation with irreversible electroporation, treated muscle cells are typically of elongated shape and their orientation may vary. Orientation of cells in electric field has been reported to affect electroporation, and hence electrodes placement and pulse parameters choice in treatments for achieving homogeneous effect in tissue is important. We investigated how cell orientation influences electroporation with respect to different pulse durations (ns to ms range), both experimentally and numerically. Experimentally detected electroporation (evaluated separately for cells parallel and perpendicular to electric field) via Ca²⁺ uptake in H9c2 and AC16 cardiomyocytes was numerically modeled using the asymptotic pore equation. Results showed that cell orientation affects electroporation extent: using short, nanosecond pulses, cells perpendicular to electric field are significantly more electroporated than parallel (up to 100-times more pores formed), and with long, millisecond pulses, cells parallel to electric field are more electroporated than perpendicular (up to 1000-times more pores formed). In the range of a few microseconds, cells of both orientations were electroporated to the same extent. Using pulses of a few microseconds lends itself as a new possible strategy in achieving homogeneous electroporation in tissue with elongated cells of different orientation (e.g. electroporation-based cardiac ablation).

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.

Transmural Remodeling of Cardiac Microstructure in Aged Spontaneously Hypertensive Rats by Diffusion Tensor MRI

The long-standing high blood pressure (also known as hypertension) overworks the heart. Microstructural remodeling is a key factor of hypertensive heart disease progression. Diffusion tensor magnetic resonance imaging (DT-MRI) is a powerful tool for the rapid noninvasive nondestructive delineation of the cardiomyocyte organization. The spontaneously hypertensive rat (SHR) is a well-established model of genetic hypertension. The goal of this study was to employ high-resolution DT-MRI and the SHR animal model to assess the transmural layer-specific remodeling of myocardial microstructure associated with hypertension. Ex vivo experiments were performed on excised formalin-fixed hearts of aged SHRs (n = 4) and age-matched controls (n = 4). The DT-MRI-derived fractional anisotropy (FA), longitudinal diffusivity (λ_L), transversal diffusivity (λ_T), and mean diffusivity (MD) served as the readout parameters investigated at three transmural zones (i.e., endocardium, mesocardium, and epicardium). The helix angles (HAs) of the aggregated cardiomyocytes and the orientation of laminar sheetlets were also studied. Compared with controls, the SHRs exhibited decreased epicardial FA, while FA changes in the other two transmural regions were insignificant. No substantial differences were observed in the diffusivity parameters and the transmural course of HAs between the two groups. A consistent distribution pattern of laminar sheetlet orientation was not identified for either group. Our findings are in line with the known cellular microstructure from early painstaking histological studies. Biophysical explanations of the study outcomes are provided. In conclusion, our experimental findings indicate that the epicardial microstructure is more vulnerable to high blood pressure leading to more pronounced changes in this region during remodeling. DT-MRI is well-suited for elucidating these alterations. The revealed transmural nonuniformity of myocardial reorganization may shed light on the mechanisms of the microstructure-function relationship in hypertension progression. Our results provide insights into the management of patients with systemic arterial hypertension, thus prevent the progression toward heart failure. The findings of this study should be acknowledged by electromechanical models of the heart that simulate the specific cardiac pathology.

Microstructurally Motivated Constitutive Modeling of Heart Failure Mechanics

Heart failure (HF) is one of the leading causes of death worldwide. HF is associated with substantial microstructural remodeling, which is linked to changes in left ventricular geometry and impaired cardiac function. The role of myocardial remodeling in altering the mechanics of failing hearts remains unclear. Structurally based constitutive modeling provides an approach to improve understanding of the relationship between biomechanical function and tissue organization in cardiac muscle during HF. In this study, we used cardiac magnetic resonance imaging and extended-volume confocal microscopy to quantify the remodeling of left ventricular geometry and myocardial microstructure of healthy and spontaneously hypertensive rat hearts at the ages of 12 and 24 months. Passive cardiac mechanical function was characterized using left ventricular pressure-volume compliance measurements. We have developed a, to our knowledge, new structurally based biomechanical constitutive equation built on parameters quantified directly from collagen distributions observed in confocal images of the myocardium. Three-dimensional left ventricular finite element models were constructed from subject-specific in vivo magnetic resonance imaging data. The structurally based constitutive equation was integrated into geometrically subject-specific finite element models of the hearts and used to investigate the underlying mechanisms of ventricular dysfunction during HF. Using a single pair of material parameters for all hearts, we were able to produce compliance curves that reproduced all of the experimental compliance measurements. The value of this study is not limited to reproducing the mechanical behavior of healthy and diseased hearts, but it also provides important insights into the structure-function relationship of diseased myocardium that will help pave the way toward more effective treatments for HF.

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.

Quantitative descriptions of extended volume cardiac tissue architecture from multiple large hearts

The arrangement of cardiac cells into strand and sheet-like structures within the heart wall, confers important electrical properties onto heart tissue. Unraveling cardiomyocyte architecture in both healthy and diseased hearts is fundamental to understanding the mechanisms generating normal rhythm and arrhythmia. We analyzed five extended volume serial image stacks of normal pig left ventricular tissue. Analysis included: (1) reconstruction of original tissue volume and shape with non-linear correction maps; (2) segmentation and higher-order descriptions, areas and orientations of laminar structures through the heart wall; (3)computation of fiber directions; (4) computation of tissue connectivity using a shell filter. These measures contributed to a deeper and more objective understanding of cardiac tissue structures and their spatial variation than previously possible.

Relative identifiability of anisotropic properties from magnetic resonance elastography

Although magnetic resonance elastography (MRE) has been used to estimate isotropic stiffness in the heart, myocardium is known to have anisotropic properties. This study investigated the determinability of global transversely isotropic material parameters using MRE and finite-element modeling (FEM). A FEM-based material parameter identification method, using a displacement-matching objective function, was evaluated in a gel phantom and simulations of a left ventricular (LV) geometry with a histology-derived fiber field. Material parameter estimation was performed in the presence of Gaussian noise. Parameter sweeps were analyzed and characteristics of the Hessian matrix at the optimal solution were used to evaluate the determinability of each constitutive parameter. Four out of five material stiffness parameters (Young's modulii E1 and E3 , shear modulus G13 and damping coefficient s), which describe a transversely isotropic linear elastic material, were well determined from the MRE displacement field using an iterative FEM inversion method. However, the remaining parameter, Poisson's ratio, was less identifiable. In conclusion, Young's modulii, shear modulii and damping can theoretically be well determined from MRE data, but Poisson's ratio is not as well determined and could be set to a reasonable value for biological tissue (close to 0.5).

[The antagonistic function of the heart muscle sustains the autoregulation according to Frank and Starling : Part I: Structure and function of heart muscle]

In the tradition of Harvey and according to Otto Frank the heart muscle structure is arranged in a strictly tangential fashion hence all contractile forces act in the direction of ventricular ejection. In contrast, morphology confirms that the heart consists of a 3-dimensional network of muscle fibers with up to two fifths of the chains of aggregated myocytes deviating from a tangential alignment at variable angles. Accordingly, the myocardial systolic forces contain, in addition to a constrictive also a (albeit smaller) radially acting component. Using needle force probes we have correspondingly measured an unloading type of force in a tangential direction and an auxotonic type in dilatative transversal direction of the ventricular walls to show that the myocardial body contracts actively in a 3-dimensional pattern. This antagonism supports the autoregulation of heart muscle function according to Frank and Starling, preserving ventricular shape, enhances late systolic fast dilation and attenuates systolic constriction of the ventricle wall. Auxotonic dilating forces are particularly sensitive to inotropic medication. Low dose beta-blocker is able to attenuate the antagonistic activity. All myocardial components act against four components of afterload, the hemodynamic, the myostructural, the stromatogenic and the hydraulic component. This complex interplay critically complicates clinical diagnostics. Clinical implications are far-reaching (see Part II, https://doi.org/10.1007/s00059-018-4735-x).

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.

Review of Journal of Cardiovascular Magnetic Resonance (JCMR) 2015-2016 and transition of the JCMR office to Boston

The Journal of Cardiovascular Magnetic Resonance (JCMR) is the official publication of the Society for Cardiovascular Magnetic Resonance (SCMR). In 2016, the JCMR published 93 manuscripts, including 80 research papers, 6 reviews, 5 technical notes, 1 protocol, and 1 case report. The number of manuscripts published was similar to 2015 though with a 12% increase in manuscript submissions to an all-time high of 369. This reflects a decrease in the overall acceptance rate to <25% (excluding solicited reviews). The quality of submissions to JCMR continues to be high. The 2016 JCMR Impact Factor (which is published in June 2016 by Thomson Reuters) was steady at 5.601 (vs. 5.71 for 2015; as published in June 2016), which is the second highest impact factor ever recorded for JCMR. The 2016 impact factor means that the JCMR papers that were published in 2014 and 2015 were on-average cited 5.71 times in 2016.In accordance with Open-Access publishing of Biomed Central, the JCMR articles are published on-line in the order that they are accepted with no collating of the articles into sections or special thematic issues. For this reason, over the years, the Editors have felt that it is useful to annually summarize the publications into broad areas of interest or themes, so that readers can view areas of interest in a single article in relation to each other and other recent JCMR articles. The papers are presented in broad themes with previously published JCMR papers to guide continuity of thought in the journal. In addition, I have elected to open this publication with information for the readership regarding the transition of the JCMR editorial office to the Beth Israel Deaconess Medical Center, Boston and the editorial process.Though there is an author publication charge (APC) associated with open-access to cover the publisher's expenses, this format provides a much wider distribution/availability of the author's work and greater manuscript citation. For SCMR members, there is a substantial discount in the APC. I hope that you will continue to send your high quality manuscripts to JCMR for consideration. Importantly, I also ask that you consider referencing recent JCMR publications in your submissions to the JCMR and elsewhere as these contribute to our impact factor. I also thank our dedicated Associate Editors, Guest Editors, and reviewers for their many efforts to ensure that the review process occurs in a timely and responsible manner and that the JCMR continues to be recognized as the leading publication in our field.

The impact of signal-to-noise ratio, diffusion-weighted directions and image resolution in cardiac diffusion tensor imaging - insights from the ex-vivo rat heart

BACKGROUND

Cardiac diffusion tensor imaging (DTI) is limited by scan time and signal-to-noise (SNR) restrictions. This invariably leads to a trade-off between the number of averages, diffusion-weighted directions (ND), and image resolution. Systematic evaluation of these parameters is therefore important for adoption of cardiac DTI in clinical routine where time is a key constraint.

METHODS

High quality reference DTI data were acquired in five ex-vivo rat hearts. We then retrospectively set 2 ≤ SNR ≤ 97, 7 ≤ ND ≤ 61, varied the voxel volume by up to 192-fold and investigated the impact on the accuracy and precision of commonly derived parameters.

RESULTS

For maximal scan efficiency, the accuracy and precision of the mean diffusivity is optimised when SNR is maximised at the expense of ND. With typical parameter settings used clinically, we estimate that fractional anisotropy may be overestimated by up to 13% with an uncertainty of ±30%, while the precision of the sheetlet angles may be as poor as ±31°. Although the helix angle has better precision of ±14°, the transmural range of helix angles may be under-estimated by up to 30° in apical and basal slices, due to partial volume and tapering myocardial geometry.

CONCLUSIONS

These findings inform a baseline of understanding upon which further issues inherent to in-vivo cardiac DTI, such as motion, strain and perfusion, can be considered. Furthermore, the reported bias and reproducibility provides a context in which to assess cardiac DTI biomarkers.

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.

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.

Extraction of the 3D local orientation of myocytes in human cardiac tissue using X-ray phase-contrast micro-tomography and multi-scale analysis

This paper presents a methodology to access the 3D local myocyte arrangements in fresh human post-mortem heart samples. We investigated the cardiac micro-structure at a high and isotropic resolution of 3.5 µm in three dimensions using X-ray phase micro-tomography at the European Synchrotron Radiation Facility. We then processed the reconstructed volumes to extract the 3D local orientation of the myocytes using a multi-scale approach with no segmentation. We created a simplified 3D model of tissue sample made of simulated myocytes with known size and orientations, to evaluate our orientation extraction method. Afterwards, we applied it to 2D histological cuts and to eight 3D left ventricular (LV) cardiac tissue samples. Then, the variation of the helix angles, from the endocardium to the epicardium, was computed at several spatial resolutions ranging from 3.6³ mm³ to 112³ µm³. We measure an increased range of 20° to 30° from the coarsest resolution level to the finest level in the experimental samples. This result is in line with the higher values measured from histology. The displayed tractography demonstrates a rather smooth evolution of the transmural helix angle in six LV samples and a sudden discontinuity of the helix angle in two septum samples. These measurements bring a new vision of the human heart architecture from macro- to micro-scale.

Review of Journal of Cardiovascular Magnetic Resonance 2015

There were 116 articles published in the Journal of Cardiovascular Magnetic Resonance (JCMR) in 2015, which is a 14 % increase on the 102 articles published in 2014. The quality of the submissions continues to increase. The 2015 JCMR Impact Factor (which is published in June 2016) rose to 5.75 from 4.72 for 2014 (as published in June 2015), which is the highest impact factor ever recorded for JCMR. The 2015 impact factor means that the JCMR papers that were published in 2013 and 2014 were cited on average 5.75 times in 2015. The impact factor undergoes natural variation according to citation rates of papers in the 2 years following publication, and is significantly influenced by highly cited papers such as official reports. However, the progress of the journal's impact over the last 5 years has been impressive. Our acceptance rate is <25 % and has been falling because the number of articles being submitted has been increasing. In accordance with Open-Access publishing, the JCMR articles go on-line as they are accepted with no collating of the articles into sections or special thematic issues. For this reason, the Editors have felt that it is useful once per calendar year to summarize the papers for the readership into broad areas of interest or theme, so that areas of interest can be reviewed in a single article in relation to each other and other recent JCMR articles. The papers are presented in broad themes and set in context with related literature and previously published JCMR papers to guide continuity of thought in the journal. We hope that you find the open-access system increases wider reading and citation of your papers, and that you will continue to send your quality papers to JCMR for publication.

On the averaging of cardiac diffusion tensor MRI data: the effect of distance function selection

Diffusion tensor magnetic resonance imaging (DT-MRI) allows a unique insight into the microstructure of highly-directional tissues. The selection of the most proper distance function for the space of diffusion tensors is crucial in enhancing the clinical application of this imaging modality. Both linear and nonlinear metrics have been proposed in the literature over the years. The debate on the most appropriate DT-MRI distance function is still ongoing. In this paper, we presented a framework to compare the Euclidean, affine-invariant Riemannian and log-Euclidean metrics using actual high-resolution DT-MRI rat heart data. We employed temporal averaging at the diffusion tensor level of three consecutive and identically-acquired DT-MRI datasets from each of five rat hearts as a means to rectify the background noise-induced loss of myocyte directional regularity. This procedure is applied here for the first time in the context of tensor distance function selection. When compared with previous studies that used a different concrete application to juxtapose the various DT-MRI distance functions, this work is unique in that it combined the following: (i) metrics were judged by quantitative-rather than qualitative-criteria, (ii) the comparison tools were non-biased, (iii) a longitudinal comparison operation was used on a same-voxel basis. The statistical analyses of the comparison showed that the three DT-MRI distance functions tend to provide equivalent results. Hence, we came to the conclusion that the tensor manifold for cardiac DT-MRI studies is a curved space of almost zero curvature. The signal to noise ratio dependence of the operations was investigated through simulations. Finally, the 'swelling effect' occurrence following Euclidean averaging was found to be too unimportant to be worth consideration.

Diffusion tensor imaging and histology of developing hearts

Diffusion tensor imaging (DTI) has emerged as a promising method for noninvasive quantification of myocardial microstructure. However, the origin and behavior of DTI measurements during myocardial normal development and remodeling remain poorly understood. In this work, conventional and bicompartmental DTI in addition to three-dimensional histological correlation were performed in a sheep model of myocardial development from third trimester to postnatal 5 months of age. Comparing the earliest time points in the third trimester with the postnatal 5 month group, the scalar transverse diffusivities preferentially increased in both left ventricle (LV) and right ventricle (RV): secondary eigenvalues D2 increased by 54% (LV) and 36% (RV), whereas tertiary eigenvalues D3 increased by 85% (LV) and 67% (RV). The longitudinal diffusivity D1 changes were small, which led to a decrease in fractional anisotropy by 41% (LV) and 33% (RV) in 5 month versus fetal hearts. Histological analysis suggested that myocardial development is associated with hyperplasia in the early stages of the third trimester followed by myocyte growth in the later stages up to 5 months of age (increased average myocyte width by 198%, myocyte length by 128%, and decreased nucleus density by 70% between preterm and postnatal 5 month hearts.) In a few histological samples (N = 6), correlations were observed between DTI longitudinal diffusivity and myocyte length (r = 0.86, P < 0.05), and transverse diffusivity and myocyte width (r = 0.96, P < 0.01). Linear regression analysis showed that transverse diffusivities are more affected by changes in myocyte size and nucleus density changes than longitudinal diffusivities, which is consistent with predictions of classical models of diffusion in porous media. Furthermore, primary and secondary DTI eigenvectors during development changed significantly. Collectively, the findings demonstrate a role for DTI to monitor and quantify myocardial development, and potentially cardiac disease. Copyright © 2016 John Wiley & Sons, Ltd.

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.

Myofiber Architecture of the Human Atria as Revealed by Submillimeter Diffusion Tensor Imaging

BACKGROUND

Accurate knowledge of the human atrial fibrous structure is paramount in understanding the mechanisms of atrial electric function in health and disease. Thus far, such knowledge has been acquired from destructive sectioning, and there is a paucity of data about atrial fiber architecture variability in the human population.

METHODS AND RESULTS

In this study, we have developed a customized 3-dimensional diffusion tensor magnetic resonance imaging sequence on a clinical scanner that makes it possible to image an entire intact human heart specimen ex vivo at submillimeter resolution. The data from 8 human atrial specimens obtained with this technique present complete maps of the fibrous organization of the human atria. The findings demonstrate that the main features of atrial anatomy are mostly preserved across subjects although the exact location and orientation of atrial bundles vary. Using the full tractography data, we were able to cluster, visualize, and characterize the distinct major bundles in the human atria. Furthermore, quantitative characterization of the fiber angles across the atrial wall revealed that the transmural fiber angle distribution is heterogeneous throughout different regions of the atria.

CONCLUSIONS

The application of submillimeter diffusion tensor magnetic resonance imaging provides an unprecedented level of information on both human atrial structure, as well as its intersubject variability. The high resolution and fidelity of this data could enhance our understanding of structural contributions to atrial rhythm and pump disorders and lead to improvements in their targeted treatment.

Resolving Fine Cardiac Structures in Rats with High-Resolution Diffusion Tensor Imaging

Cardiac architecture is fundamental to cardiac function and can be assessed non-invasively with diffusion tensor imaging (DTI). Here, we aimed to overcome technical challenges in ex vivo DTI in order to extract fine anatomical details and to provide novel insights in the 3D structure of the heart. An integrated set of methods was implemented in ex vivo rat hearts, including dynamic receiver gain adjustment, gradient system scaling calibration, prospective adjustment of diffusion gradients, and interleaving of diffusion-weighted and non-diffusion-weighted scans. Together, these methods enhanced SNR and spatial resolution, minimised orientation bias in diffusion-weighting, and reduced temperature variation, enabling detection of tissue structures such as cell alignment in atria, valves and vessels at an unprecedented level of detail. Improved confidence in eigenvector reproducibility enabled tracking of myolaminar structures as a basis for segmentation of functional groups of cardiomyocytes. Ex vivo DTI facilitates acquisition of high quality structural data that complements readily available in vivo cardiac functional and anatomical MRI. The improvements presented here will facilitate next generation virtual models integrating micro-structural and electro-mechanical properties of the heart.

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.

Cardiac MR imaging: current status and future direction

Coronary artery disease is currently a worldwide epidemic with increasing impact on healthcare systems. Magnetic resonance imaging (MRI) sequences give complementary information on LV function, regional perfusion, angiogenesis, myocardial viability and orientations of myocytes. T2-weighted short-tau inversion recovery (T2-STIR), fat suppression and black blood sequences have been frequently used for detecting edematous area at risk (AAR) of infarction. T2 mapping, however, indicated that the edematous reaction in acute myocardial infarct (AMI) is not stable and warranted the use of edematous area in evaluating therapies. On the other hand, cine MRI demonstrated reproducible data on LV function in healthy volunteers and LV remodeling in patients. Noninvasive first pass perfusion, using exogenous tracer (gadolinium-based contrast media) and arterial spin labeling MRI, using endogenous tracer (water), are sensitive and useful techniques for evaluating myocardial perfusion and angiogenesis. Recently, new strategies have been developed to quantify myocardial viability using T1-mapping and equilibrium contrast enhanced MR techniques because existing delayed contrast enhancement MRI (DE-MRI) sequences are limited in detecting patchy microinfarct and diffuse fibrosis. These new techniques were successfully used for characterizing diffuse myocardial fibrosis associated with myocarditis, amyloidosis, sarcoidosis heart failure, aortic hypertrophic cardiomyopathy, congenital heart disease, restrictive cardiomyopathy, arrhythmogenic right ventricular dysplasia and hypertension). Diffusion MRI provides information regarding microscopic tissue structure, while diffusion tensor imaging (DTI) helps to characterize the myocardium and monitor the process of LV remodeling after AMI. Novel trends in hybrid imaging, such as cardiac positron emission tomography (PET)/MRI and optical imaging/MRI, are recently under intensive investigation. With the promise of higher spatial-temporal resolution and 3D coverage in the near future, cardiac MRI will be an indispensible tool in the diagnosis of cardiac diseases, coronary intervention and myocardial therapeutic delivery.