Probing cardiomyocyte mobility with multi-phase cardiac diffusion tensor MRI

The effects of field strength on stimulated echo and motion compensated spin echo diffusion tensor cardiovascular magnetic resonance sequences

BACKGROUND

In-vivo diffusion tensor CMR (DT-CMR) is an emerging technique for microstructural tissue characterisation in the myocardium. Most studies are performed at 3T, where higher signal to noise ratio (SNR) should benefit this signal starved method. However, a few studies have suggested that DT-CMR is possible at 1.5T, where EPI artefacts may be less severe and 1.5T hardware is more widely available.

METHODS

We recruited 20 healthy volunteers and performed mid-ventricular short axis DT-CMR at 1.5 T and 3 T. Acquisitions were performed at peak systole and end-diastole using both stimulated echo acquisition mode (STEAM) and motion compensated spin-echo (MCSE) sequences at matched spatial resolutions. DT-CMR parameters were averaged over the LV and compared between 1.5 T and 3 T sequences using both datasets with and without the blow reference data included.

RESULTS

Eleven (1.5T) and 12 (3T) diastolic MCSE acquisitions were rejected as the helix angle (HA) demonstrated <50% normal appearance circumferentially or the acquisition was abandoned due to poor image quality; a maximum of one acquisition was rejected for other datasets. Subjective HA map quality was significantly better at 3T than 1.5T for STEAM (p<0.05), but not for MCSE and other DT-CMR quality measures were consistent with improvements in STEAM at 3T over 1.5T. When blow data was excluded, no significant differences in mean diffusivity were observed between field strengths, but fractional anisotropy was significantly higher at 1.5T than 3T for STEAM systole (p<0.05). Absolute second eigenvector orientation (E2A, sheetlet angle) was significantly higher at 1.5T than 3T for MCSE systole and STEAM diastole, but significantly lower for STEAM systole (all p<0.05). Transmural HA distribution was less steep at 1.5T than 3T for STEAM diastole data (p<0.05). SNR in the blow images was higher at 3T than 1.5T for all acquisitions (p<0.05).

CONCLUSION

While 3T provides benefits in terms of SNR, both STEAM and MCSE can be performed at 1.5T. However, MCSE is unreliable in diastole at both field strengths and STEAM benefits from the improved SNR at 3T over 1.5T. Future clinical research studies may be able to leverage the wider availability of 1.5T CMR hardware where MCSE acquisitions are desirable.

Cardiac DTI using short-axis PROPELLER: A feasibility study

PURPOSE

We aimed to develop a free-breathing (FB) cardiac DTI (cDTI) method based on short-axis PROPELLER (SAP) and M2 motion compensated spin-echo EPI (SAP-M2-EPI) to mitigate geometric distortion and eliminate aliasing in acquired diffusion-weighted (DW) images, particularly in patients with a higher body mass index (BMI).

THEORY AND METHODS

The study involved 10 healthy volunteers whose BMI values fell into specific categories: BMI <25 (4 volunteers), 25< BMI <28 (5 volunteers), and BMI >30 (1 volunteer). We compared DTI parameters, including fractional anisotropy (FA), mean diffusivity (MD), and helix angle transmurality (HAT), between SAP-M2-EPI and M2-ssEPI. To evaluate the performance of SAP-M2-EPI in reducing geometric distortions in the left ventricle (LV) compared to CINE and M2-ssEPI, we utilized the DICE similarity coefficient (DSC) and assessed misregistration area.

RESULTS

In all volunteers, SAP-M2-EPI yielded high-quality LV DWIs without aliasing, demonstrating significantly reduced geometric distortion (with an average DSC of 0.92 and average misregistration area of 90 mm2) and diminished signal loss due to bulk motion when compared to M2-ssEPI. DTI parameter maps exhibited consistent patterns across slices without motion related artifacts.

CONCLUSION

SAP-M2-EPI facilitates free-breathing cDTI of the entire LV, effectively eliminating aliasing and minimizing geometric distortion compared to M2-ssEPI. Furthermore, it preserves accurate quantification of myocardial microstructure.

Diffusion Tensor Phenomapping of the Healthy and Pressure-Overloaded Human Heart

Current techniques to image the microstructure of the heart with diffusion tensor MRI (DTI) are highly under-resolved. We present a technique to improve the spatial resolution of cardiac DTI by almost 10-fold and leverage this to measure local gradients in cardiomyocyte alignment or helix angle (HA). We further introduce a phenomapping approach based on voxel-wise hierarchical clustering of these gradients to identify distinct microstructural microenvironments in the heart. Initial development was performed in healthy volunteers (n=8). Thereader, subjects with severe but well-compensated aortic stenosis (AS, n=10) were compared to age-matched controls (CTL, n=10). Radial HA gradient was significantly reduced in AS (8.0±0.8°/mm vs. 10.2±1.8°/mm, p=0.001) but the other HA gradients did not change significantly. Four distinct microstructural clusters could be idenJfied in both the CTL and AS subjects and did not differ significantly in their properties or distribution. Despite marked hypertrophy, our data suggest that the myocardium in well-compensated AS can maintain its microstructural coherence. The described phenomapping approach can be used to characterize microstructural plasticity and perturbation in any organ system and disease.

Characterizing the Myoarchitecture of the Supraspinatus and Infraspinatus Muscles With MRI Using Diffusion Tensor Imaging

BACKGROUND

The societal cost of shoulder disabilities in our aging society keeps rising. Providing biomarkers of early changes in the microstructure of rotator cuff (RC) muscles might improve surgical planning. Elevation angle (E1A) and pennation angle (PA) assessed by ultrasound change with RC tears. Furthermore, ultrasounds lack repeatability.

PURPOSE

To propose a repeatable framework to quantify the myocyte angulation in RC muscles.

STUDY TYPE

Prospective.

SUBJECTS

Six asymptomatic healthy volunteers (1 female aged 30 years; 5 males, mean age 35 years, range 25-49 years), who underwent three repositioned scanning sessions (10 minutes apart) of the right infraspinatus muscle (ISPM) and supraspinatus muscle (SSPM).

FIELD STRENGTH/SEQUENCE

3-T, T1-weighted and diffusion tensor imaging (DTI; 12 gradient encoding directions, b-values of 500 and 800 s/mm2 ).

ASSESSMENT

Each voxel was binned in percentage of depth defined by the shortest distance in the antero-posterior direction (manual delineation), i.e. the radial axis. A second order polynomial fit for PA across the muscle depth was used, while E1A described a sigmoid across depth: E 1 A sig = E 1 A range × sigmf 1 : 100 % depth , - EA 1 grad   ,   E 1 A asym + E 1 A shift .

STATISTICAL TESTS

Repeatability was assessed with the nonparametric Wilcoxon's rank-sum test for paired comparisons across repeated scans in each volunteer for each anatomical muscle region and across repeated measures of the radial axis. A P-value <0.05 was considered statistically significant.

RESULTS

In the ISPM, E1A was constantly negative, became helicoidal, then mainly positive across the antero-posterior depth, respective at the caudal, central and cranial regions. In the SSPM, posterior myocytes ran more parallel to the intramuscular tendon ( PA ≈ 0 ° ), while anterior myocytes inserted with a pennation angle ( PA ≈ - 20 ° ). E1A and PA were repeatable in each volunteer (error < 10%). Intra-repeatability of the radial axis was achieved (error < 5%).

DATA CONCLUSION

ElA and PA in the proposed framework of the ISPM and SSPM are repeatable with DTI. Variations of myocyte angulation in the ISPM and SSPM can be quantified across volunteers.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Cardiac structure discontinuities revealed by ex-vivo microstructural characterization. A focus on the basal inferoseptal left ventricle region

BACKGROUND

While the microstructure of the left ventricle (LV) has been largely described, only a few studies investigated the right ventricular insertion point (RVIP). It was accepted that the aggregate cardiomyocytes organization was much more complex due to the intersection of the ventricular cavities but a precise structural characterization in the human heart was lacking even if clinical phenotypes related to right ventricular wall stress or arrhythmia were observed in this region.

METHODS

MRI-derived anatomical imaging (150 µm3) and diffusion tensor imaging (600 µm3) were performed in large mammalian whole hearts (human: N = 5, sheep: N = 5). Fractional anisotropy, aggregate cardiomyocytes orientations and tractography were compared within both species. Aggregate cardiomyocytes orientation on one ex-vivo sheep whole heart was then computed using structure tensor imaging (STI) from 21 µm isotropic acquisition acquired with micro computed tomography (MicroCT) imaging. Macroscopic and histological examination were performed. Lastly, experimental cardiomyocytes orientation distribution was then compared to the usual rule-based model using electrophysiological (EP) modeling. Electrical activity was modeled with the monodomain formulation.

RESULTS

The RVIP at the level of the inferior ventricular septum presented a unique arrangement of aggregate cardiomyocytes. An abrupt, mid-myocardial change in cardiomyocytes orientation was observed, delimiting a triangle-shaped region, present in both sheep and human hearts. FA's histogram distribution (mean ± std: 0.29 ± 0.06) of the identified region as well as the main dimension (22.2 mm ± 5.6 mm) was found homogeneous across samples and species. Averaged volume is 0.34 cm3 ± 0.15 cm3. Both local activation time (LAT) and morphology of pseudo-ECGs were strongly impacted with delayed LAT and change in peak-to-peak amplitude in the simulated wedge model.

CONCLUSION

The study was the first to describe the 3D cardiomyocytes architecture of the basal inferoseptal left ventricle region in human hearts and identify the presence of a well-organized aggregate cardiomyocytes arrangement and cardiac structural discontinuities. The results might offer a better appreciation of clinical phenotypes like RVIP-late gadolinium enhancement or uncommon idiopathic ventricular arrhythmias (VA) originating from this region.

A tomographic microscopy-compatible Langendorff system for the dynamic structural characterization of the cardiac cycle

Introduction

Cardiac architecture has been extensively investigated ex vivo using a broad spectrum of imaging techniques. Nevertheless, the heart is a dynamic system and the structural mechanisms governing the cardiac cycle can only be unveiled when investigating it as such.

Methods

This work presents the customization of an isolated, perfused heart system compatible with synchrotron-based X-ray phase contrast imaging (X-PCI).

Results

Thanks to the capabilities of the developed setup, it was possible to visualize a beating isolated, perfused rat heart for the very first time in 4D at an unprecedented 2.75 μm pixel size (10.6 μm spatial resolution), and 1 ms temporal resolution.

Discussion

The customized setup allows high-spatial resolution studies of heart architecture along the cardiac cycle and has thus the potential to serve as a tool for the characterization of the structural dynamics of the heart, including the effects of drugs and other substances able to modify the cardiac cycle.

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.

Comparison of interpolation methods of predominant cardiomyocyte orientation from in vivo and ex vivo cardiac diffusion tensor imaging data

Cardiac electrophysiology and cardiac mechanics both depend on the average cardiomyocyte long-axis orientation. In the realm of personalized medicine, knowledge of the patient-specific changes in cardiac microstructure plays a crucial role. Patient-specific computational modelling has emerged as a tool to better understand disease progression. In vivo cardiac diffusion tensor imaging (cDTI) is a vital tool to non-destructively measure the average cardiomyocyte long-axis orientation in the heart. However, cDTI suffers from long scan times, rendering volumetric, high-resolution acquisitions challenging. Consequently, interpolation techniques are needed to populate bio-mechanical models with patient-specific average cardiomyocyte long-axis orientations. In this work, we compare five interpolation techniques applied to in vivo and ex vivo porcine input data. We compare two tensor interpolation approaches, one rule-based approximation, and two data-driven, low-rank models. We demonstrate the advantage of tensor interpolation techniques, resulting in lower interpolation errors than do low-rank models and rule-based methods adapted to cDTI data. In an ex vivo comparison, we study the influence of three imaging parameters that can be traded off against acquisition time: in-plane resolution, signal to noise ratio, and number of acquired short-axis imaging slices.

Cardiac Diffusion Tensor Biomarkers of Chronic Infarction Based on In Vivo Data

In vivo cardiac diffusion tensor imaging (cDTI) data were acquired in swine subjects six to ten weeks post-myocardial infarction (MI) to identify microstructural-based biomarkers of MI. Diffusion tensor invariants, diffusion tensor eigenvalues, and radial diffusivity (RD) are evaluated in the infarct, border, and remote myocardium, and compared with extracellular volume fraction (ECV) and native T1 values. Additionally, to aid the interpretation of the experimental results, the diffusion of water molecules was numerically simulated as a function of ECV. Finally, findings based on in vivo measures were confirmed using higher-resolution and higher signal-to-noise data acquired ex vivo in the same subjects. Mean diffusivity, diffusion tensor eigenvalues, and RD increased in the infarct and border regions compared to remote myocardium, while fractional anisotropy decreased. Secondary (e₂) and tertiary (e₃) eigenvalues increased more significantly than the primary eigenvalue in the infarct and border regions. These findings were confirmed by the diffusion simulations. Although ECV presented the largest increase in infarct and border regions, e₂, e₃, and RD increased the most among non-contrast-based biomarkers. RD is of special interest as it summarizes the changes occurring in the radial direction and may be more robust than e₂ or e₃ alone.

Impact of Intraventricular Septal Fiber Orientation on Cardiac Electromechanical Function

Cardiac fiber direction is an important factor determining the propagation of electrical activity, as well as the development of mechanical force. In this article, we imaged the ventricles of several species with special attention to the intraventricular septum to determine the functional consequences of septal fiber organization. First, we identified a dual-layer organization of the fiber orientation in the intraventricular septum of ex vivo sheep hearts using diffusion tensor imaging at high field MRI. To expand the scope of the results, we investigated the presence of a similar fiber organization in five mammalian species (rat, canine, pig, sheep, and human) and highlighted the continuity of the layer with the moderator band in large mammalian species. We implemented the measured septal fiber fields in three-dimensional electromechanical computer models to assess the impact of the fiber orientation. The downward fibers produced a diamond activation pattern superficially in the right ventricle. Electromechanically, there was very little change in pressure volume loops although the stress distribution was altered. In conclusion, we clarified that the right ventricular septum has a downwardly directed superficial layer in larger mammalian species, which can have modest effects on stress distribution.

NEW & NOTEWORTHY A dual-layer organization of the fiber orientation in the intraventricular septum was identified in ex vivo hearts of large mammals. The RV septum has a downwardly directed superficial layer that is continuous with the moderator band. Electrically, it produced a diamond activation pattern. Electromechanically, little change in pressure volume loops were noticed but stress distribution was altered. Fiber distribution derived from diffusion tensor imaging should be considered for an accurate strain and stress analysis.

Personalization of biomechanical simulations of the left ventricle by in-vivo cardiac DTI data: Impact of fiber interpolation methods

Simulations of cardiac electrophysiology and mechanics have been reported to be sensitive to the microstructural anisotropy of the myocardium. Consequently, a personalized representation of cardiac microstructure is a crucial component of accurate, personalized cardiac biomechanical models. In-vivo cardiac Diffusion Tensor Imaging (cDTI) is a non-invasive magnetic resonance imaging technique capable of probing the heart's microstructure. Being a rather novel technique, issues such as low resolution, signal-to noise ratio, and spatial coverage are currently limiting factors. We outline four interpolation techniques with varying degrees of data fidelity, different amounts of smoothing strength, and varying representation error to bridge the gap between the sparse in-vivo data and the model, requiring a 3D representation of microstructure across the myocardium. We provide a workflow to incorporate in-vivo myofiber orientation into a left ventricular model and demonstrate that personalized modelling based on fiber orientations from in-vivo cDTI data is feasible. The interpolation error is correlated with a trend in personalized parameters and simulated physiological parameters, strains, and ventricular twist. This trend in simulation results is consistent across material parameter settings and therefore corresponds to a bias introduced by the interpolation method. This study suggests that using a tensor interpolation approach to personalize microstructure with in-vivo cDTI data, reduces the fiber uncertainty and thereby the bias in the simulation results.