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

Use of the intrusion angle to describe the radial orientation of local cardiomyocytes in the left ventricle

The projected transverse angle and the nonprojected intrusion angle can be used to describe the radial orientation of local cardiomyocytes in the left ventricle wall, although to date their descriptive relevance has not been demonstrated. This paper compares the evolution of the transverse angle and the intrusion angle in five left ventricle wall samples, and investigates in more detail their respective behaviors when the nonprojected helical angle varies. We show that the intrusion angle avoids the "projection" effect, and contrary to the transverse angle, it remains stable whatever the values taken by the nonprojected helical angle, even when this approaches 90°. The intrusion angle is the better choice, rather than the transverse angle, in describing the radial orientation of local cardiomyocytes. Furthermore, the oscillation of the intrusion angle in the samples is assessed, whose results indicate that the intrusion angle's oscillation amplitude and period are regional and related to the local tissue architecture.

Measuring cardiomyocyte cellular characteristics in cardiac hypertrophy using diffusion-weighted MRI

PURPOSE

This paper presents a hierarchical modeling approach for estimating cardiomyocyte major and minor diameters and intracellular volume fraction (ICV) using diffusion-weighted MRI (DWI) data in ex vivo mouse hearts.

METHODS

DWI data were acquired on two healthy controls and two hearts 3 weeks post transverse aortic constriction (TAC) using a bespoke diffusion scheme with multiple diffusion times (Δ $$ \Delta $$ ), q-shells and diffusion encoding directions. Firstly, a bi-exponential tensor model was fitted separately at each diffusion time to disentangle the dependence on diffusion times from diffusion weightings, that is, b-values. The slow-diffusing component was attributed to the restricted diffusion inside cardiomyocytes. ICV was then extrapolated atΔ = 0 $$ \Delta =0 $$ using linear regression. Secondly, given the secondary and the tertiary diffusion eigenvalue measurements for the slow-diffusing component obtained at different diffusion times, major and minor diameters were estimated assuming a cylinder model with an elliptical cross-section (ECS). High-resolution three-dimensional synchrotron X-ray imaging (SRI) data from the same specimen was utilized to evaluate the biophysical parameters.

RESULTS

Estimated parameters using DWI data were (control 1/control 2 vs. TAC 1/TAC 2): major diameter-17.4μ $$ \mu $$ m/18.0μ $$ \mu $$ m versus 19.2μ $$ \mu $$ m/19.0μ $$ \mu $$ m; minor diameter-10.2μ $$ \mu $$ m/9.4μ $$ \mu $$ m versus 12.8μ $$ \mu $$ m/13.4μ $$ \mu $$ m; and ICV-62%/62% versus 68%/47%. These findings were consistent with SRI measurements.

CONCLUSION

The proposed method allowed for accurate estimation of biophysical parameters suggesting cardiomyocyte diameters as sensitive biomarkers of hypertrophy in the heart.

Microstructural and Microvascular Phenotype of Sarcomere Mutation Carriers and Overt Hypertrophic Cardiomyopathy

BACKGROUND

In hypertrophic cardiomyopathy (HCM), myocyte disarray and microvascular disease (MVD) have been implicated in adverse events, and recent evidence suggests that these may occur early. As novel therapy provides promise for disease modification, detection of phenotype development is an emerging priority. To evaluate their utility as early and disease-specific biomarkers, we measured myocardial microstructure and MVD in 3 HCM groups-overt, either genotype-positive (G+LVH+) or genotype-negative (G-LVH+), and subclinical (G+LVH-) HCM-exploring relationships with electrical changes and genetic substrate.

METHODS

This was a multicenter collaboration to study 206 subjects: 101 patients with overt HCM (51 G+LVH+ and 50 G-LVH+), 77 patients with G+LVH-, and 28 matched healthy volunteers. All underwent 12-lead ECG, quantitative perfusion cardiac magnetic resonance imaging (measuring myocardial blood flow, myocardial perfusion reserve, and perfusion defects), and cardiac diffusion tensor imaging measuring fractional anisotropy (lower values expected with more disarray), mean diffusivity (reflecting myocyte packing/interstitial expansion), and second eigenvector angle (measuring sheetlet orientation).

RESULTS

Compared with healthy volunteers, patients with overt HCM had evidence of altered microstructure (lower fractional anisotropy, higher mean diffusivity, and higher second eigenvector angle; all P<0.001) and MVD (lower stress myocardial blood flow and myocardial perfusion reserve; both P<0.001). Patients with G-LVH+ were similar to those with G+LVH+ but had elevated second eigenvector angle (P<0.001 after adjustment for left ventricular hypertrophy and fibrosis). In overt disease, perfusion defects were found in all G+ but not all G- patients (100% [51/51] versus 82% [41/50]; P=0.001). Patients with G+LVH- compared with healthy volunteers similarly had altered microstructure, although to a lesser extent (all diffusion tensor imaging parameters; P<0.001), and MVD (reduced stress myocardial blood flow [P=0.015] with perfusion defects in 28% versus 0 healthy volunteers [P=0.002]). Disarray and MVD were independently associated with pathological electrocardiographic abnormalities in both overt and subclinical disease after adjustment for fibrosis and left ventricular hypertrophy (overt: fractional anisotropy: odds ratio for an abnormal ECG, 3.3, P=0.01; stress myocardial blood flow: odds ratio, 2.8, P=0.015; subclinical: fractional anisotropy odds ratio, 4.0, P=0.001; myocardial perfusion reserve odds ratio, 2.2, P=0.049).

CONCLUSIONS

Microstructural alteration and MVD occur in overt HCM and are different in G+ and G- patients. Both also occur in the absence of hypertrophy in sarcomeric mutation carriers, in whom changes are associated with electrocardiographic abnormalities. Measurable changes in myocardial microstructure and microvascular function are early-phenotype biomarkers in the emerging era of disease-modifying therapy.

Cardiac q-space trajectory imaging by motion-compensated tensor-valued diffusion encoding in human heart in vivo

PURPOSE

Tensor-valued diffusion encoding can probe more specific features of tissue microstructure than what is available by conventional diffusion weighting. In this work, we investigate the technical feasibility of tensor-valued diffusion encoding at high b-values with q-space trajectory imaging (QTI) analysis, in the human heart in vivo.

METHODS

Ten healthy volunteers were scanned on a 3T scanner. We designed time-optimal gradient waveforms for tensor-valued diffusion encoding (linear and planar) with second-order motion compensation. Data were analyzed with QTI. Normal values and repeatability were investigated for the mean diffusivity (MD), fractional anisotropy (FA), microscopic FA (μFA), isotropic, anisotropic and total mean kurtosis (MKi, MKa, and MKt), and orientation coherence (Cc ). A phantom, consisting of two fiber blocks at adjustable angles, was used to evaluate sensitivity of parameters to orientation dispersion and diffusion time.

RESULTS

QTI data in the left ventricular myocardium were MD = 1.62 ± 0.07 μm2 /ms, FA = 0.31 ± 0.03, μFA = 0.43 ± 0.07, MKa = 0.20 ± 0.07, MKi = 0.13 ± 0.03, MKt = 0.33 ± 0.09, and Cc  = 0.56 ± 0.22 (mean ± SD across subjects). Phantom experiments showed that FA depends on orientation dispersion, whereas μFA was insensitive to this effect.

CONCLUSION

We demonstrated the first tensor-valued diffusion encoding and QTI analysis in the heart in vivo, along with first measurements of myocardial μFA, MKi, MKa, and Cc . The methodology is technically feasible and provides promising novel biomarkers for myocardial tissue characterization.

Slice-specific tracking for free-breathing diffusion tensor cardiac MRI

BACKGROUND

Diffusion tensor cardiac magnetic resonance (DT-CMR) imaging has great potential to characterize myocardial microarchitecture. However, its accuracy is limited by respiratory and cardiac motion and long scan times. Here, we develop and evaluate a slice-specific tracking method to improve the efficiency and accuracy of DT-CMR acquisition during free breathing.

METHODS

Coronal images were obtained along with signals from a diaphragmatic navigator. Respiratory and slice displacements were obtained from the navigator signals and coronal images, respectively, and these displacements were fitted with a linear model to obtain the slice-specific tracking factors. This method was evaluated in DT-CMR examinations of 17 healthy subjects, and the results were compared with those obtained using a fixed tracking factor of 0.6. DT-CMR with breath-holding was used for reference. Quantitative and qualitative evaluation methods were used to analyze the performance of the slice-specific tracking method and the consistency between the obtained diffusion parameters.

RESULTS

In the study, the slice-specific tracking factors showed an upward trend from the basal to the apical slice. Residual in-plane movements were lower in slice-specific tracking than in fixed-factor tracking (RMSE: 2.748 ± 1.171 versus 5.983 ± 2.623, P < 0.001). The diffusion parameters obtained using slice-specific tracking were not significantly different from those obtained from breath-holding acquisition (P > 0.05).

CONCLUSION

In free-breathing DT-CMR imaging, the slice-specific tracking method reduced misalignment of the acquired slices. The diffusion parameters obtained using this approach were consistent with those obtained with the breath-holding technique.

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.

Detailed quantification of cardiac ventricular myocardial architecture in the embryonic and fetal mouse heart by application of structure tensor analysis to high resolution episcopic microscopic data

The mammalian heart, which is one of the first organs to form and function during embryogenesis, develops from a simple tube into a complex organ able to efficiently pump blood towards the rest of the body. The progressive growth of the compact myocardium during embryonic development is accompanied by changes in its structural complexity and organisation. However, how myocardial myoarchitecture develops during embryogenesis remain poorly understood. To date, analysis of heart development has focused mainly on qualitative descriptions using selected 2D histological sections. High resolution episcopic microscopy (HREM) is a novel microscopic imaging technique that enables to obtain high-resolution three-dimensional images of the heart and perform detailed quantitative analyses of heart development. In this work, we performed a detailed characterization of the development of myocardial architecture in wildtype mice, from E14.5 to E18.5, by means of structure tensor analysis applied to HREM images of the heart. Our results shows that even at E14.5, myocytes are already aligned, showing a gradual change in their helical angle from positive angulation in the endocardium towards negative angulation in the epicardium. Moreover, there is gradual increase in the degree of myocardial organisation concomitant with myocardial growth. However, the development of the myoarchitecture is heterogeneous showing regional differences between ventricles, ventricular walls as well as between myocardial layers, with different growth patterning between the endocardium and epicardium. We also found that the percentage of circumferentially arranged myocytes within the LV significantly increases with gestational age. Finally, we found that fractional anisotropy (FA) within the LV gradually increases with gestational age, while the FA within RV remains unchanged.

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.

Biomedical Imaging in Experimental Models of Cardiovascular Disease

Major advances in biomedical imaging have occurred over the last 2 decades and now allow many physiological, cellular, and molecular processes to be imaged noninvasively in small animal models of cardiovascular disease. Many of these techniques can be also used in humans, providing pathophysiological context and helping to define the clinical relevance of the model. Ultrasound remains the most widely used approach, and dedicated high-frequency systems can obtain extremely detailed images in mice. Likewise, dedicated small animal tomographic systems have been developed for magnetic resonance, positron emission tomography, fluorescence imaging, and computed tomography in mice. In this article, we review the use of ultrasound and positron emission tomography in small animal models, as well as emerging contrast mechanisms in magnetic resonance such as diffusion tensor imaging, hyperpolarized magnetic resonance, chemical exchange saturation transfer imaging, magnetic resonance elastography and strain, arterial spin labeling, and molecular imaging.

Assessing Myocardial Microstructure With Biophysical Models of Diffusion MRI

Biophysical models are a promising means for interpreting diffusion weighted magnetic resonance imaging (DW-MRI) data, as they can provide estimates of physiologically relevant parameters of microstructure including cell size, volume fraction, or dispersion. However, their application in cardiac microstructure mapping (CMM) has been limited. This study proposes seven new two-compartment models with combination of restricted cylinder models and a diffusion tensor to represent intra- and extracellular spaces, respectively. Three extended versions of the cylinder model are studied here: cylinder with elliptical cross section (ECS), cylinder with Gamma distributed radii (GDR), and cylinder with Bingham distributed axes (BDA). The proposed models were applied to data in two fixed mouse hearts, acquired with multiple diffusion times, q-shells and diffusion encoding directions. The cylinderGDR-pancake model provided the best performance in terms of root mean squared error (RMSE) reducing it by 25% compared to diffusion tensor imaging (DTI). The cylinderBDA-pancake model represented anatomical findings closest as it also allows for modelling dispersion. High-resolution 3D synchrotron X-ray imaging (SRI) data from the same specimen was utilized to evaluate the biophysical models. A novel tensor-based registration method is proposed to align SRI structure tensors to the MR diffusion tensors. The consistency between SRI and DW-MRI parameters demonstrates the potential of compartment models in assessing physiologically relevant parameters.

Exploring arterial tissue microstructural organization using non-Gaussian diffusion magnetic resonance schemes

The purpose of this study was to characterize the alterations in microstructural organization of arterial tissue using higher-order diffusion magnetic resonance schemes. Three porcine carotid artery models namely; native, collagenase treated and decellularized, were used to estimate the contribution of collagen and smooth muscle cells (SMC) on diffusion signal attenuation using gaussian and non-gaussian schemes. The samples were imaged in a 7 T preclinical scanner. High spatial and angular resolution diffusion weighted images (DWIs) were acquired using two multi-shell (max b-value = 3000 s/mm²) acquisition protocols. The processed DWIs were fitted using monoexponential, stretched-exponential, kurtosis and bi-exponential schemes. Directionally variant and invariant microstructural parametric maps of the three artery models were obtained from the diffusion schemes. The parametric maps were used to assess the sensitivity of each diffusion scheme to collagen and SMC composition in arterial microstructural environment. The inter-model comparison showed significant differences across the considered models. The bi-exponential scheme based slow diffusion compartment (Ds) was highest in the absence of collagen, compared to native and decellularized microenvironments. In intra-model comparison, kurtosis along the radial direction was the highest. Overall, the results of this study demonstrate the efficacy of higher order dMRI schemes in mapping constituent specific alterations in arterial microstructure.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

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.

Phenotyping hypertrophic cardiomyopathy using cardiac diffusion magnetic resonance imaging: the relationship between microvascular dysfunction and microstructural changes

Aims

Microvascular dysfunction in hypertrophic cardiomyopathy (HCM) is predictive of clinical decline, however underlying mechanisms remain unclear. Cardiac diffusion tensor imaging (cDTI) allows in vivo characterization of myocardial microstructure by quantifying mean diffusivity (MD), fractional anisotropy (FA) of diffusion, and secondary eigenvector angle (E2A). In this cardiac magnetic resonance (CMR) study, we examine associations between perfusion and cDTI parameters to understand the sequence of pathophysiology and the interrelation between vascular function and underlying microstructure.

Methods and results

Twenty HCM patients underwent 3.0T CMR which included: spin-echo cDTI, adenosine stress and rest perfusion mapping, cine-imaging, and late gadolinium enhancement (LGE). Ten controls underwent cDTI. Myocardial perfusion reserve (MPR), MD, FA, E2A, and wall thickness were calculated per segment and further divided into subendocardial (inner 50%) and subepicardial (outer 50%) regions. Segments with wall thickness ≤11 mm, MPR ≥2.2, and no visual LGE were classified as ‘normal’. Compared to controls, ‘normal’ HCM segments had increased MD (1.61 ± 0.09 vs. 1.46 ± 0.07 × 10⁻³ mm²/s, P = 0.02), increased E2A (60 ± 9° vs. 38 ± 12°, P < 0.001), and decreased FA (0.29 ± 0.04 vs. 0.35 ± 0.02, P = 0.002). Across all HCM segments, subendocardial regions had higher MD and lower MPR than subepicardial (MD_endo 1.61 ± 0.08 × 10⁻³ mm²/s vs. MD_epi 1.56 ± 0.18 × 10⁻³ mm²/s, P = 0.003, MPR_endo 1.85 ± 0.83, MPR_epi 2.28 ± 0.87, P < 0.0001).

Conclusion

In HCM patients, even in segments with normal wall thickness, normal perfusion, and no scar, diffusion is more isotropic than in controls, suggesting the presence of underlying cardiomyocyte disarray. Increased E2A suggests the myocardial sheetlets adopt hypercontracted angulation in systole. Increased MD, most notably in the subendocardium, is suggestive of regional remodelling which may explain the reduced subendocardial blood flow.

A Multi-Scale Computational Model for the Rat Ventricle: Construction, Parallelization, and Applications

BACKGROUND

Cardiovascular diseases are the top killer of human beings. The ventricular arrhythmia, as a type of malignant cardiac arrhythmias, typically leads to death if not treated within minutes. The multi-scale virtual heart provides an idealized tool for exploring the underlying mechanisms, by means of incorporating abundant experimental data at the level of ion channels and analyzing the subsequent pathological changes at organ levels. However, there are few studies on building a virtual heart model for rats-a species most widely used in experiments.

OBJECTIVE

To build a multi-scale computational model for rats, with detailed methodology for the model construction, computational optimization, and its applications.

METHODS

First, approaches for building multi-scale models ranging from cellular to 3-D organ levels are introduced, with detailed descriptions of handling the ventricular myocardium heterogeneity, geometry processing, and boundary conditions, etc. Next, for dealing with the expensive computational costs of 3-D models, optimization approaches including an optimized representation and a GPU-based parallelization method are introduced. Finally, methods for reproducing of some key phenomenon (e.g., electrocardiograph, spiral/scroll waves) are demonstrated.

RESULTS

Three types of heterogeneity, including the transmural heterogeneity, the interventricular heterogeneity, and the base-apex heterogeneity are incorporated into the model. The normal and reentrant excitation waves, as well as the corresponding pseudo-ECGs are reproduced by the constructed ventricle model. In addition, the temporal and spatial vulnerability to reentry arrhythmias are quantified based on the evaluation experiments of vulnerable window and the critical length.

CONCLUSIONS

The constructed multi-scale rat ventricle model is able to reproduce both the physiological and the pathological phenomenon in different scales. Evaluation experiments suggest that the apex is the most susceptible area to arrhythmias. The model can be a promising tool for the investigation of arrhythmogenesis and the screening of anti-arrhythmic drugs.

Passive mechanical properties in healthy and infarcted rat left ventricle characterised via a mixture model

During the cardiac cycle, electrical excitation is coupled with mechanical response of the myocardium. Besides the active contraction, passive mechanics plays an important role, and its behaviour differs in healthy and diseased hearts as well as among different animal species. The aim of this study is the characterisation of passive mechanical properties in healthy and infarcted rat myocardium by means of mechanical testing and subsequent parameter fitting. Elasticity assessments via uniaxial extension tests are performed on healthy and infarcted tissue samples from left ventricular rat myocardium. In order to fully characterise the orthotropic cardiac tissue, our experimental data are combined with other previously published tests in rats - shear tests on healthy myocardium and equibiaxial tests on infarcted tissue. In a first step, we calibrate the Holzapfel-Ogden strain energy function in the healthy case. Sa far, this orthotropic constitutive law for the passive myocardium has been fitted to experimental data in several species, however there is a lack of an appropriate parameter set for the rat. With our determined parameters, a finite element simulation of the end-diastolic filling is performed. In a second step, we propose a model for the infarcted tissue. It is represented as a mixture of intact myocardium and a transversely isotropic scar structure. In our mechanical experiments, the tissue after myocardial infarction shows significantly stiffer behaviour than in the healthy case, and the stiffness correlates with the amount of fibrosis. A similar relationship is observed in the computational simulation of the end-diastolic filling. We conclude that our new proposed material model can capture the behaviour of two kinds of tissues - healthy and infarcted rat myocardium, and its calibration with the fitted parameters represents the experimental data well.

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.

Abnormal skeletal muscle blood flow, contractile mechanics and fibre morphology in a rat model of obese-HFpEF

KEY POINTS

Heart failure is characterised by limb and respiratory muscle impairments that limit functional capacity and quality of life. However, compared with heart failure with reduced ejection fraction (HFrEF), skeletal muscle alterations induced by heart failure with preserved ejection fraction (HFpEF) remain poorly explored. Here we report that obese-HFpEF induces multiple skeletal muscle alterations in the rat hindlimb, including impaired muscle mechanics related to shortening velocity, fibre atrophy, capillary loss, and an impaired blood flow response to contractions that implies a perfusive oxygen delivery limitation. We also demonstrate that obese-HFpEF is characterised by diaphragmatic alterations similar to those caused by denervation - atrophy in Type IIb/IIx (fast/glycolytic) fibres and hypertrophy in Type I (slow/oxidative) fibres. These findings extend current knowledge in HFpEF skeletal muscle physiology, potentially underlying exercise intolerance, which may facilitate future therapeutic approaches.

ABSTRACT

Peripheral skeletal muscle and vascular alterations induced by heart failure with preserved ejection fraction (HFpEF) remain poorly identified, with limited therapeutic targets. This study used a cardiometabolic obese-HFpEF rat model to comprehensively phenotype skeletal muscle mechanics, blood flow, microvasculature and fibre atrophy. Lean (n = 8) and obese-HFpEF (n = 8) ZSF1 rats were compared. Skeletal muscles (soleus and diaphragm) were assessed for in vitro contractility (isometric and isotonic properties) alongside indices of fibre-type cross-sectional area, myosin isoform, and capillarity, and estimated muscle PO₂ . In situ extensor digitorum longus (EDL) contractility and femoral blood flow were assessed. HFpEF soleus demonstrated lower absolute maximal force by 22%, fibre atrophy by 24%, a fibre-type shift from I to IIa, and a 17% lower capillary-to-fibre ratio despite increased capillary density (all P < 0.05) with preserved muscle PO₂ (P = 0.115) and isometric specific force (P > 0.05). Soleus isotonic properties (shortening velocity and power) were impaired by up to 17 and 22%, respectively (P < 0.05), while the magnitude of the exercise hyperaemia was attenuated by 73% (P = 0.012) in line with higher muscle fatigue by 26% (P = 0.079). Diaphragm alterations (P < 0.05) included Type IIx fibre atrophy despite Type I/IIa fibre hypertrophy, with increased indices of capillarity alongside preserved contractile properties during isometric, isotonic, and cyclical contractions. In conclusion, obese-HFpEF rats demonstrated blunted skeletal muscle blood flow during contractions in parallel to microvascular structural remodelling, fibre atrophy, and isotonic contractile dysfunction in the locomotor muscles. In contrast, diaphragm phenotype remained well preserved. This study identifies numerous muscle-specific impairments that could exacerbate exercise intolerance in obese-HFpEF.

Insight Into Myocardial Microstructure of Athletes and Hypertrophic Cardiomyopathy Patients Using Diffusion Tensor Imaging

BACKGROUND

Hypertrophic cardiomyopathy (HCM) remains the commonest cause of sudden cardiac death among young athletes. Differentiating between physiologically adaptive left ventricular (LV) hypertrophy observed in athletes' hearts and pathological HCM remains challenging. By quantifying the diffusion of water molecules, diffusion tensor imaging (DTI) MRI allows voxelwise characterization of myocardial microstructure.

PURPOSE

To explore microstructural differences between healthy volunteers, athletes, and HCM patients using DTI.

STUDY TYPE

Prospective cohort.

POPULATION

Twenty healthy volunteers, 20 athletes, and 20 HCM patients.

FIELD STRENGTH/SEQUENCE

3T/DTI spin echo.

ASSESSMENT

In-house MatLab software was used to derive mean diffusivity (MD) and fractional anisotropy (FA) as markers of amplitude and anisotropy of the diffusion of water molecules, and secondary eigenvector angles (E2A)-reflecting the orientations of laminar sheetlets.

STATISTICAL TESTS

Independent samples t-tests were used to detect statistical significance between any two cohorts. Analysis of variance was utilized for detecting the statistical difference between the three cohorts. Statistical tests were two-tailed. A result was considered statistically significant at P ≤ 0.05.

RESULTS

DTI markers were significantly different between HCM, athletes, and volunteers. HCM patients had significantly higher global MD and E2A, and significantly lower FA than athletes and volunteers. (MDHCM = 1.52 ± 0.06 × 10-3 mm2 /s, MDAthletes = 1.49 ± 0.03 × 10-3 mm2 /s, MDvolunteers = 1.47 ± 0.02 × 10-3 mm2 /s, P < 0.05; E2AHCM = 58.8 ± 4°, E2Aathletes = 47 ± 5°, E2Avolunteers = 38.5 ± 7°, P < 0.05; FAHCM = 0.30 ± 0.02, FAAthletes = 0.35 ± 0.02, FAvolunteers = 0.36 ± 0.03, P < 0.05). HCM patients had significantly higher E2A in their thickest segments compared to the remote (E2Athickest = 66.8 ± 7, E2Aremote = 51.2 ± 9, P < 0.05).

DATA CONCLUSION

DTI depicts an increase in amplitude and isotropy of diffusion in the myocardium of HCM compared to athletes and volunteers as reflected by increased MD and decreased FA values. While significantly higher E2A values in HCM and athletes reflect steeper configurations of the myocardial sheetlets than in volunteers, HCM patients demonstrated an eccentric rise in E2A in their thickest segments, while athletes demonstrated a concentric rise. Further studies are required to determine the diagnostic capabilities of DTI.

EVIDENCE LEVEL

1 TECHNICAL EFFICACY STAGE: 2.

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.

Microarchitecture of the hearts in term and former-preterm lambs using diffusion tensor imaging

Diffusion tensor imaging (DTI) is an MRI technique that can be used to map cardiomyocyte tracts and estimate local cardiomyocyte and sheetlet orientation within the heart. DTI measures diffusion distances of water molecules within the myocardium, where water diffusion generally occurs more freely along the long axis of cardiomyocytes and within the extracellular matrix, but is restricted by cell membranes such that transverse diffusion is limited. DTI can be undertaken in fixed hearts and it allows the three-dimensional mapping of the cardiac microarchitecture, including cardiomyocyte organization, within the whole heart. The objective of this study was to use DTI to compare the cardiac microarchitecture and cardiomyocyte organization in archived fixed left ventricles of lambs that were born either preterm (n = 5) or at term (n = 7), at a postnatal timepoint equivalent to about 6 years of age in children. Although the findings support the feasibility of retrospective DTI scanning of fixed hearts, several hearts were excluded from DTI analysis because of poor scan quality, such as ghosting artifacts. The preliminary findings from viable DTI scans (n = 3/group) suggest that the extracellular compartment is altered and that there is an immature microstructural phenotype early in postnatal life in the LV of lambs born preterm. Our findings support a potential time-efficient imaging role for DTI in detecting abnormal changes in the microstructure of fixed hearts of former-preterm neonates, although further investigation into factors that affect scan quality is required.

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.

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.

Improved compressed sensing and super-resolution of cardiac diffusion MRI with structure-guided total variation

Purpose

Structure‐guided total variation is a recently introduced prior that allows reconstruction of images using knowledge of the location and orientation of edges in a reference image. In this work, we demonstrate the advantages of a variant of structure‐guided total variation known as directional total variation (DTV), over traditional total variation (TV), in the context of compressed‐sensing reconstruction and super‐resolution.

Methods

We compared TV and DTV in retrospectively undersampled ex vivo diffusion tensor imaging and diffusion spectrum imaging data from healthy, sham, and hypertrophic rat hearts.

Results

In compressed sensing at an undersampling factor of 8, the RMS error of mean diffusivity and fractional anisotropy relative to the fully sampled ground truth were 44% and 20% lower in DTV compared with TV. In super‐resolution, these values were 29% and 14%, respectively. Similarly, we observed improvements in helix angle, transverse angle, sheetlet elevation, and sheetlet azimuth. The RMS error of the diffusion kurtosis in the undersampled data relative to the ground truth was uniformly lower (22% on average) with DTV compared to TV.

Conclusion

Acquiring one fully sampled non‐diffusion‐weighted image and 10 diffusion‐weighted images at 8× undersampling would result in an 80% net reduction in data needed. We demonstrate efficacy of the DTV algorithm over TV in reducing data sampling requirements, which can be translated into higher apparent resolution and potentially shorter scan times. This method would be equally applicable in diffusion MRI applications outside the heart.

Investigation of the Role of Myocyte Orientations in Cardiac Arrhythmia Using Image-Based Models

Cardiac electrical excitation-propagation is influenced by myocyte orientations (cellular organization). Quantitatively understanding this relationship presents a significant research challenge, especially during arrhythmias in which excitation patterns become complex. Tissue-scale simulations of cardiac electrophysiology, incorporating both dynamic action potential behavior and image-based myocardial architecture, provide an approach to investigate three-dimensional (3D) propagation of excitation waves in the heart. In this study, we aimed to assess the importance of natural variation in myocyte orientations on cardiac arrhythmogenesis using 3D tissue electrophysiology simulations. Three anatomical models (i.e., describing myocyte orientations) of healthy rat ventricles—obtained using diffusion tensor imaging at 100 μm resolution—were registered to a single biventricular geometry (i.e., a single cardiac shape), in which the myocyte orientations could be represented by each of the diffusion tensor imaging data sets or by an idealized rule-based description. The Fenton-Karma cellular excitation model was modified to reproduce rat ventricular action potential duration restitution to create reaction-diffusion cardiac electrophysiology models. Over 250 3D simulations were performed to investigate the effects of myocyte orientations on the following: 1) ventricular activation, 2) location-dependent arrhythmia induction via rapid pacing, and 3) dynamics of re-entry averaged over multiple episodes. It was shown that 1) myocyte orientation differences manifested themselves in local activation times, but the influence on total activation time was small; 2) differences in myocyte orientations could critically affect the inducibility and persistence of arrhythmias for specific stimulus-location/cycle-length combinations; and 3) myocyte orientations alone could be an important determinant of scroll wave break, although no significant differences were observed in averaged arrhythmia dynamics between the four myocyte orientation scenarios considered. Our results show that myocyte orientations are an important determinant of arrhythmia inducibility, persistence, and scroll wave break. These findings suggest that where specificity is desired (for example, when predicting location-dependent, patient-specific arrhythmia inducibility), subject-specific myocyte orientations may be important.

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.

A comparison of two quasi-static computational models for assessment of intra-myocardial injection as a therapeutic strategy for heart failure

Myocardial infarction, or heart attack, is the leading cause of mortality globally. Although the treatment of myocardial infarct has improved significantly, scar tissue that persists can often lead to increased stress and adverse remodeling of surrounding tissue and ultimately to heart failure. Intra-myocardial injection of biomaterials represents a potential treatment to attenuate remodeling, mitigate degeneration, and reverse the disease process in the tissue. In vivo experiments on animal models have shown functional benefits of this therapeutic strategy. However, a poor understanding of the optimal injection pattern, volume, and material properties has acted as a barrier to its widespread clinical adoption. In this study, we developed two quasistatic finite element simulations of the left ventricle to investigate the mechanical effect of intra-myocardial injection. The first model employed an idealized left ventricular geometry with rule-based cardiomyocyte orientation. The second model employed a subject-specific left ventricular geometry with cardiomyocyte orientation from diffusion tensor magnetic resonance imaging. Both models predicted cardiac parameters including ejection fraction, systolic wall thickening, and ventricular twist that matched experimentally reported values. All injection simulations showed cardiomyocyte stress attenuation, offering an explanation for the mechanical reinforcement benefit associated with injection. The study also enabled a comparison of injection location and the corresponding effect on cardiac performance at different stages of the cardiac cycle. While the idealized model has lower fidelity, it predicts cardiac function and differentiates the effects of injection location. Both models represent versatile in silico tools to guide optimal strategy in terms of injection number, volume, site, and material properties.

Myoarchitectural disarray of hypertrophic cardiomyopathy begins pre-birth

Myoarchitectural disarray – the multiscalar disorganisation of myocytes, is a recognised histopathological hallmark of adult human hypertrophic cardiomyopathy (HCM). It occurs before the establishment of left ventricular hypertrophy (LVH) but its early origins and evolution around the time of birth are unknown. Our aim is to investigate whether myoarchitectural abnormalities in HCM are present in the fetal heart. We used wild‐type, heterozygous and homozygous hearts (n = 56) from a Mybpc3‐targeted knock‐out HCM mouse model and imaged the 3D micro‐structure by high‐resolution episcopic microscopy. We developed a novel structure tensor approach to extract, display and quantify myocyte orientation and its local angular uniformity by helical angle, angle of intrusion and myoarchitectural disarray index, respectively, immediately before and after birth. In wild‐type, we demonstrate uniformity of orientation of cardiomyocytes with smooth transitions of helical angle transmurally both before and after birth but with traces of disarray at the septal insertion points of the right ventricle. In comparison, heterozygous mice free of LVH, and homozygous mice showed not only loss of the normal linear helical angulation transmural profiles observed in wild‐type but also fewer circumferentially arranged myocytes at birth. Heterozygous and homozygous showed more disarray with a wider distribution than in wild‐type before birth. In heterozygous mice, disarray was seen in the anterior, septal and inferior walls irrespective of stage, whereas in homozygous mice it extended to the whole LV circumference including the lateral wall. In conclusion, myoarchitectural disarray is detectable in the fetal heart of an HCM mouse model before the development of LVH.

High-Resolution Ex Vivo Microstructural MRI After Restoring Ventricular Geometry via 3D Printing

Computational modeling of the heart requires accurately incorporating both gross anatomical detail and local microstructural information. Together, these provide the necessary data to build 3D meshes for simulation of cardiac mechanics and electrophysiology. Recent MRI advances make it possible to measure detailed heart motion in vivo, but in vivo microstructural imaging of the heart remains challenging. Consequently, the most detailed measurements of microstructural organization and microanatomical infarct details are obtained ex vivo. The objective of this work was to develop and evaluate a new method for restoring ex vivo ventricular geometry to match the in vivo configuration. This approach aids the integration of high-resolution ex vivo microstructural information with in vivo motion measurements. The method uses in vivo cine imaging to generate surface meshes, then creates a 3D printed left ventricular (LV) blood pool cast and a pericardial mold to restore the ex vivo cardiac geometry to a mid-diastasis reference configuration. The method was evaluated in healthy (N = 7) and infarcted (N = 3) swine. Dice similarity coefficients were calculated between in vivo and ex vivo images for the LV cavity (0.93 ± 0.01), right ventricle (RV) cavity (0.80 ± 0.05), and the myocardium (0.72 ± 0.04). The R ² coefficient between in vivo and ex vivo LV and RV cavity volumes were 0.95 and 0.91, respectively. These results suggest that this method adequately restores ex vivo geometry to match in vivo geometry. This approach permits a more precise incorporation of high-resolution ex vivo anatomical and microstructural data into computational models that use in vivo data for simulation of cardiac mechanics and electrophysiology.

Comprehensive Analysis of Animal Models of Cardiovascular Disease using Multiscale X-Ray Phase Contrast Tomography

Cardiovascular diseases (CVDs) affect the myocardium and vasculature, inducing remodelling of the heart from cellular to whole organ level. To assess their impact at micro and macroscopic level, multi-resolution imaging techniques that provide high quality images without sample alteration and in 3D are necessary: requirements not fulfilled by most of current methods. In this paper, we take advantage of the non-destructive time-efficient 3D multiscale capabilities of synchrotron Propagation-based X-Ray Phase Contrast Imaging (PB-X-PCI) to study a wide range of cardiac tissue characteristics in one healthy and three different diseased rat models. With a dedicated image processing pipeline, PB-X-PCI images are analysed in order to show its capability to assess different cardiac tissue components at both macroscopic and microscopic levels. The presented technique evaluates in detail the overall cardiac morphology, myocyte aggregate orientation, vasculature changes, fibrosis formation and nearly single cell arrangement. Our results agree with conventional histology and literature. This study demonstrates that synchrotron PB-X-PCI, combined with image processing tools, is a powerful technique for multi-resolution structural investigation of the heart ex-vivo. Therefore, the proposed approach can improve the understanding of the multiscale remodelling processes occurring in CVDs, and the comprehensive and fast assessment of future interventional approaches.

A rule-based method to model myocardial fiber orientation in cardiac biventricular geometries with outflow tracts

Rule-based methods are often used for assigning fiber orientation to cardiac anatomical models. However, existing methods have been developed using data mostly from the left ventricle. As a consequence, fiber information obtained from rule-based methods often does not match histological data in other areas of the heart such as the right ventricle, having a negative impact in cardiac simulations beyond the left ventricle. In this work, we present a rule-based method where fiber orientation is separately modeled in each ventricle following observations from histology. This allows to create detailed fiber orientation in specific regions such as the endocardium of the right ventricle, the interventricular septum, and the outflow tracts. We also carried out electrophysiological simulations involving these structures and with different fiber configurations. In particular, we built a modeling pipeline for creating patient-specific volumetric meshes of biventricular geometries, including the outflow tracts, and subsequently simulate the electrical wavefront propagation in outflow tract ventricular arrhythmias with different origins for the ectopic focus. The resulting simulations with the proposed rule-based method showed a very good agreement with clinical parameters such as the 10 ms isochrone ratio in a cohort of nine patients suffering from this type of arrhythmia. The developed modeling pipeline confirms its potential for an in silico identification of the site of origin in outflow tract ventricular arrhythmias before clinical intervention.

Whole heart detailed and quantitative anatomy, myofibre structure and vasculature from X-ray phase-contrast synchrotron radiation-based micro computed tomography

Background

While individual cardiac myocytes only have a limited ability to shorten, the heart efficiently pumps a large volume-fraction thanks to a cell organization in a complex 3D fibre structure. Subclinical subtle cardiac structural remodelling is often present before symptoms arise. Understanding and early detection of these subtle changes is crucial for diagnosis and prevention. Additionally, personalized computational modelling requires knowledge on the multi-scale structure of the whole heart and vessels.

Methods and Results

We developed a rapid acquisition together with visualization and quantification methods of the integrated microstructure of whole in-vitro rodents hearts using synchrotron based X-ray phase-contrast tomography. These images are formed not only by X-ray absorption by the tissue but also by wave propagation phenomena, enhancing structural information, thus allowing to raise tissue contrast to an unprecedented level. We used a (ex-vivo) normal rat heart and fetal rabbit hearts suffering intrauterine growth restriction as a model of subclinical cardiac remodelling to illustrate the strengths and potential of the technique. For comparison, histology and diffusion tensor magnetic resonance imaging was performed.

Conclusions

We have developed a novel, high resolution, image acquisition, and quantification approach to study a whole in-vitro heart at myofibre resolution, providing integrated 3D structural information at microscopic level without any need of tissue slicing and processing. This superior imaging approach opens up new possibilities for a systems approach towards analysing cardiac structure and function, providing rapid acquisition of quantitative microstructure of the heart in a near native state.

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.

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.

Monte Carlo Simulations of Diffusion Weighted MRI in Myocardium: Validation and Sensitivity Analysis

A model of cardiac microstructure and diffusion MRI is presented, and compared with experimental data from ex vivo rat hearts. The model includes a simplified representation of individual cells, with physiologically correct cell size and orientation, as well as intra- to extracellular volume ratio. Diffusion MRI is simulated using a Monte Carlo model and realistic MRI sequences. The results show good correspondence between the simulated and experimental MRI signals. Similar patterns are observed in the eigenvalues of the diffusion tensor, the mean diffusivity (MD), and the fractional anisotropy (FA). A sensitivity analysis shows that the diffusivity is the dominant influence on all three eigenvalues of the diffusion tensor, the MD, and the FA. The area and aspect ratio of the cell cross-section affect the secondary and tertiary eigenvalues, and hence the FA. Within biological norms, the cell length, volume fraction of cells, and rate of change of helix angle play a relatively small role in influencing tissue diffusion. Results suggest that the model could be used to improve understanding of the relationship between cardiac microstructure and diffusion MRI measurements, as well as in testing and refinement of cardiac diffusion MRI protocols.

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.

Evaluation of non-Gaussian diffusion in cardiac MRI

PURPOSE

The diffusion tensor model assumes Gaussian diffusion and is widely applied in cardiac diffusion MRI. However, diffusion in biological tissue deviates from a Gaussian profile as a result of hindrance and restriction from cell and tissue microstructure, and may be quantified better by non-Gaussian modeling. The aim of this study was to investigate non-Gaussian diffusion in healthy and hypertrophic hearts.

METHODS

Thirteen rat hearts (five healthy, four sham, four hypertrophic) were imaged ex vivo. Diffusion-weighted images were acquired at b-values up to 10,000 s/mm² . Models of diffusion were fit to the data and ranked based on the Akaike information criterion.

RESULTS

The diffusion tensor was ranked best at b-values up to 2000 s/mm² but reflected the signal poorly in the high b-value regime, in which the best model was a non-Gaussian "beta distribution" model. Although there was considerable overlap in apparent diffusivities between the healthy, sham, and hypertrophic hearts, diffusion kurtosis and skewness in the hypertrophic hearts were more than 20% higher in the sheetlet and sheetlet-normal directions.

CONCLUSION

Non-Gaussian diffusion models have a higher sensitivity for the detection of hypertrophy compared with the Gaussian model. In particular, diffusion kurtosis may serve as a useful biomarker for characterization of disease and remodeling in the heart. Magn Reson Med 78:1174-1186, 2017. © 2016 International Society for Magnetic Resonance in Medicine.