Diffuse myocardial fibrosis precedes subclinical functional myocardial impairment and provides prognostic information in systemic sclerosis

AIMS

Myocardial involvement is common in patients with systemic sclerosis (SSc) and causes myocardial fibrosis and subtle ventricular dysfunction. However, the temporal onset of myocardial involvement during the progression of the disease and its prognostic value are yet unknown. We used cardiovascular magnetic resonance (CMR) to investigate subclinical functional impairment and diffuse myocardial fibrosis in patients with very early diagnosis of SSc (VEDOSS) and established SSc and examined whether this was associated with mortality.

METHODS AND RESULTS

One hundred and ten SSc patients (86 established SSc, 24 VEDOSS) and 15 healthy controls were prospectively recruited. The patients were followed-up for a median duration of 7.0 years (interquartile range 6.0-7.3 years). Study subjects underwent CMR including assessment of myocardial fibrosis [native T1 and extracellular volume (ECV)] and measurement of global longitudinal (GLS) and circumferential (GCS) myocardial strain. Native T1 values and ECV were elevated in VEDOSS and SSc patients compared with controls (P < 0.001). GLS was similar in VEDOSS and controls but significantly impaired in patients with established SSc (P < 0.001). GCS was similar over all groups (P = 0.88). There were 12 deaths during follow-up. Elevated native T1 [hazard ratio (HR) 5.8, 95% confidence interval (CI): 1.7-20.4; P = 0.006] and reduced GLS (HR 6.1, 95% CI: 1.3-29.9; P = 0.038) identified subjects with increased risk of death. Only native T1 was predictive for cardiovascular mortality (P < 0.001).

CONCLUSION

Subclinical myocardial involvement first manifests as diffuse myocardial fibrosis identified by the expansion of ECV and increased native T1 in VEDOSS patients while subtle functional impairment only occurs in established SSc. Native T1 and GLS have prognostic value for all-cause mortality in SSc patients.

Characterization and correction of diffusion gradient-induced eddy currents in second-order motion-compensated echo-planar and spiral cardiac DTI

PURPOSE

Very high gradient amplitudes played out over extended time intervals as required for second-order motion-compensated cardiac DTI may violate the assumption of a linear time-invariant gradient system model. The aim of this work was to characterize diffusion gradient-related system nonlinearity and propose a correction approach for echo-planar and spiral spin-echo motion-compensated cardiac DTI.

METHODS

Diffusion gradient-induced eddy currents of 9 diffusion directions were characterized at b values of 150 s/mm² and 450 s/mm² for a 1.5 Tesla system and used to correct phantom, ex vivo, and in vivo motion-compensated cardiac DTI data acquired with echo-planar and spiral trajectories. Predicted trajectories were calculated using gradient impulse response function and diffusion gradient strength- and direction-dependent zeroth- and first-order eddy current responses. A reconstruction method was implemented using the predicted <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>k</mml:mi></mml:mrow> <mml:annotation>$$ k $$</mml:annotation></mml:semantics> </mml:math> -space trajectories to additionally include off-resonances and concomitant fields. Resulting images were compared to a reference reconstruction omitting diffusion gradient-induced eddy current correction.

RESULTS

Diffusion gradient-induced eddy currents exhibited nonlinear effects when scaling up the gradient amplitude and could not be described by a 3D basis alone. This indicates that a gradient impulse response function does not suffice to describe diffusion gradient-induced eddy currents. Zeroth- and first-order diffusion gradient-induced eddy current effects of up to -1.7 rad and -16 to +12 rad/m, respectively, were identified. Zeroth- and first-order diffusion gradient-induced eddy current correction yielded improved image quality upon image reconstruction.

CONCLUSION

The proposed approach offers correction of diffusion gradient-induced zeroth- and first-order eddy currents, reducing image distortions to promote improvements of second-order motion-compensated spin-echo cardiac DTI.

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.

Dos and don'ts in large animal models of aortic insufficiency

Aortic insufficiency caused by paravalvular leakage (PVL) is one of the most feared complications following transcatheter aortic valve replacement (TAVI) in patients. Domestic pigs (Sus scrofa domestica) are a popular large animal model to study such conditions and develop novel diagnostic and therapeutic techniques. However, the models based on prosthetic valve implantation are time intensive, costly, and often hamper further hemodynamic measurements such as PV loop and 4D MRI flow by causing implantation-related wall motion abnormalities and degradation of MR image quality. This study describes in detail, the establishment of a minimally invasive porcine model suitable to study the effects of mild-to-moderate “paravalvular“ aortic regurgitation on left ventricular (LV) performance and blood flow patterns, particularly under the influence of altered afterload, preload, inotropic state, and heart rate. Six domestic pigs (Swiss large white, female, 60–70 kg of body weight) were used to establish this model. The defects on the hinge point of aortic leaflets and annulus were created percutaneously by the pierce-and-dilate technique either in the right coronary cusp (RCC) or in the non-coronary cusp (NCC). The hemodynamic changes as well as LV performance were recorded by PV loop measurements, while blood flow patterns were assessed by 4D MRI. LV performance was additionally challenged by pharmaceutically altering cardiac inotropy, chronotropy, and afterload. The presented work aims to elaborate the dos and don'ts in porcine models of aortic insufficiency and intends to steepen the learning curve for researchers planning to use this or similar models by giving valuable insights ranging from animal selection to vascular access choices, placement of PV Loop catheter, improvement of PV loop data acquisition and post-processing and finally the induction of paravalvular regurgitation of the aortic valve by a standardized and reproducible balloon induced defect in a precisely targeted region of the aortic valve.

Quantification of left ventricular strain and torsion by joint analysis of 3D tagging and cine MR images

Cardiovascular magnetic resonance (CMR) imaging is the gold standard for the non-invasive assessment of left-ventricular (LV) function. Prognostic value of deformation metrics extracted directly from regular SSFP CMR images has been shown by numerous studies in the clinical setting, but with some limitations to detect torsion of the myocardium. Tagged CMR introduces trackable features in the myocardium that allow for the assessment of local myocardial deformation, including torsion; it is, however, limited in the quantification of radial strain, which is a decisive metric for assessing the contractility of the heart. In order to improve SSFP-only and tagged-only approaches, we propose to combine the advantages of both image types by fusing global shape motion obtained from SSFP images with the local deformation obtained from tagged images. To this end, tracking is first performed on SSFP images, and subsequently, the resulting motion is utilized to mask and track tagged data. Our implementation is based on a recent finite element-based motion tracking tool with mechanical regularization. Joint SSFP and tagged images registration performance is assessed based on deformation metrics including LV strain and twist using human and in-house porcine datasets. Results show that joint analysis of SSFP and 3DTAG images provides better quantification of LV strain and twist as either data source alone.

Validation of cardiac diffusion tensor imaging sequences: A multicentre test-retest phantom study

Cardiac diffusion tensor imaging (DTI) is an emerging technique for the in vivo characterisation of myocardial microstructure, and there is a growing need for its validation and standardisation. We sought to establish the accuracy, precision, repeatability and reproducibility of state-of-the-art pulse sequences for cardiac DTI among 10 centres internationally. Phantoms comprising 0%-20% polyvinylpyrrolidone (PVP) were scanned with DTI using a product pulsed gradient spin echo (PGSE; N = 10 sites) sequence, and a custom motion-compensated spin echo (SE; N = 5) or stimulated echo acquisition mode (STEAM; N = 5) sequence suitable for cardiac DTI in vivo. A second identical scan was performed 1-9 days later, and the data were analysed centrally. The average mean diffusivities (MDs) in 0% PVP were (1.124, 1.130, 1.113) x 10^-3  mm² /s for PGSE, SE and STEAM, respectively, and accurate to within 1.5% of reference data from the literature. The coefficients of variation in MDs across sites were 2.6%, 3.1% and 2.1% for PGSE, SE and STEAM, respectively, and were similar to previous studies using only PGSE. Reproducibility in MD was excellent, with mean differences in PGSE, SE and STEAM of (0.3 ± 2.3, 0.24 ± 0.95, 0.52 ± 0.58) x 10^-5  mm² /s (mean ± 1.96 SD). We show that custom sequences for cardiac DTI provide accurate, precise, repeatable and reproducible measurements. Further work in anisotropic and/or deforming phantoms is warranted.

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.

Detection of Intramyocardial Iron in Patients Following ST-Elevation Myocardial Infarction Using Cardiac Diffusion Tensor Imaging

Background

Intramyocardial hemorrhage (IMH) following ST‐elevation myocardial infarction (STEMI) is associated with poor prognosis. In cardiac magnetic resonance (MR), T2* mapping is the reference standard for detecting IMH while cardiac diffusion tensor imaging (cDTI) can characterize myocardial architecture via fractional anisotropy (FA) and mean diffusivity (MD) of water molecules. The value of cDTI in the detection of IMH is not currently known.

Hypothesis

cDTI can detect IMH post‐STEMI.

Study Type

Prospective.

Subjects

A total of 50 patients (20% female) scanned at 1‐week (V1) and 3‐month (V2) post‐STEMI.

Field Strength/Sequence

A 3.0 T; inversion‐recovery T1‐weighted‐imaging, multigradient‐echo T2* mapping, spin‐echo cDTI.

Assessment

T2* maps were analyzed to detect IMH (defined as areas with T2* < 20 msec within areas of infarction). cDTI images were co‐registered to produce averaged diffusion‐weighted‐images (DWIs), MD, and FA maps; hypointense areas were manually planimetered for IMH quantification.

Statistics

On averaged DWI, the presence of hypointense signal in areas matching IMH on T2* maps constituted to true‐positive detection of iron. Independent samples t‐tests were used to compare regional cDTI values. Results were considered statistically significant at P ≤ 0.05.

Results

At V1, 24 patients had IMH on T2*. On averaged DWI, all 24 patients had hypointense signal in matching areas. IMH size derived using averaged‐DWI was nonsignificantly greater than from T2* (2.0 ± 1.0 cm² vs 1.89 ± 0.96 cm², P = 0.69). Compared to surrounding infarcted myocardium, MD was significantly reduced (1.29 ± 0.20 × 10⁻³ mm²/sec vs 1.75 ± 0.16 × 10⁻³ mm²/sec) and FA was significantly increased (0.40 ± 0.07 vs 0.23 ± 0.03) within areas of IMH. By V2, all 24 patients with acute IMH continued to have hypointense signals on averaged‐DWI in the affected area. T2* detected IMH in 96% of these patients. Overall, averaged‐DWI had 100% sensitivity and 96% specificity for the detection of IMH.

Data Conclusion

This study demonstrates that the parameters MD and FA are susceptible to the paramagnetic properties of iron, enabling cDTI to detect IMH.

Evidence Level

1

Technical Efficacy

Stage 2

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.

Pathology and Advanced Imaging—Characterization of a Congenital Cardiac Defect and Complex Hemodynamics in a Pig: A Case Report

Domestic pigs are widely used in cardiovascular research as the porcine circulatory system bears a remarkable resemblance to that of humans. In order to reduce variability, only clinically healthy animals enter the study as their health status is assessed in entry examination. Like humans, pigs can also suffer from congenital heart disease, such as an atrial septal defect (ASD), which often remains undetected. Due to the malformation of the endocardial cushion during organ development, mitral valve defects (e.g., mitral clefts) are sometimes associated with ASDs, further contributing to hemodynamic instability. In this work, we report an incidental finding of a hemodynamically highly relevant ASD in the presence of incompetent mitral and tricuspid valves, in an asymptomatic, otherwise healthy juvenile pig. In-depth characterization of the cardiac blood flow by four-dimensional (4D) flow magnetic resonance imaging (MRI) revealed a prominent diastolic left-to-right and discrete systolic right-to-left shunt, resulting in a pulmonary-to-systemic flow ratio of 1.8. Severe mitral (15 mL/stroke) and tricuspid (22 mL/stroke) regurgitation further reduced cardiac output. Pathological examination confirmed the presence of an ostium primum ASD and found a serous cyst of lymphatic origin that was filled with clear fluid partially occluding the ASD. A large mitral cleft was identified as the most likely cause of severe regurgitation, and histology showed mild to moderate endocardiosis in the coaptation area of both atrio-ventricular valves. In summary, although not common, congenital heart defects could play a role as a cause of experimental variability or even intra-experimental mortality when working with apparently heathy, juvenile pigs.

In-silico study of accuracy and precision of left-ventricular strain quantification from 3D tagged MRI

Cardiac Magnetic Resonance Imaging (MRI) allows quantifying myocardial tissue deformation and strain based on the tagging principle. In this work, we investigate accuracy and precision of strain quantification from synthetic 3D tagged MRI using equilibrated warping. To this end, synthetic biomechanical left-ventricular tagged MRI data with varying tag distance, spatial resolution and signal-to-noise ratio (SNR) were generated and processed to quantify errors in radial, circumferential and longitudinal strains relative to ground truth. Results reveal that radial strain is more sensitive to image resolution and noise than the other strain components. The study also shows robustness of quantifying circumferential and longitudinal strain in the presence of geometrical inconsistencies of 3D tagged data. In conclusion, our study points to the need for higher-resolution 3D tagged MRI than currently available in practice in order to achieve sufficient accuracy of radial strain quantification.

Cardiovascular magnetic resonance imaging of functional and microstructural changes of the heart in a longitudinal pig model of acute to chronic myocardial infarction

BACKGROUND

We examined the dynamic response of the myocardium to infarction in a longitudinal porcine study using relaxometry, functional as well as diffusion cardiovascular magnetic resonance (CMR). We sought to compare non contrast CMR methods like relaxometry and in-vivo diffusion to contrast enhanced imaging and investigate the link of microstructural and functional changes in the acute and chronically infarcted heart.

METHODS

CMR was performed on five myocardial infarction pigs and four healthy controls. In the infarction group, measurements were obtained 2 weeks before 90 min occlusion of the left circumflex artery, 6 days after ischemia and at 5 as well as 9 weeks as chronic follow-up. The timing of measurements was replicated in the control cohort. Imaging consisted of functional cine imaging, 3D tagging, T2 mapping, native as well as gadolinium enhanced T1 mapping, cardiac diffusion tensor imaging, and late gadolinium enhancement imaging.

RESULTS

Native T1, extracellular volume (ECV) and mean diffusivity (MD) were significantly elevated in the infarcted region while fractional anisotropy (FA) was significantly reduced. During the transition from acute to chronic stages, native T1 presented minor changes (< 3%). ECV as well as MD increased from acute to the chronic stages compared to baseline: ECV: 125 ± 24% (day 6) 157 ± 24% (week 5) 146 ± 60% (week 9), MD: 17 ± 7% (day 6) 33 ± 14% (week 5) 29 ± 15% (week 9) and FA was further reduced: - 31 ± 10% (day 6) - 38 ± 8% (week 5) - 36 ± 14% (week 9). T2 as marker for myocardial edema was significantly increased in the ischemic area only during the acute stage (83 ± 3 ms infarction vs. 58 ± 2 ms control p < 0.001 and 61 ± 2 ms in the remote area p < 0.001). The analysis of functional imaging revealed reduced left ventricular ejection fraction, global longitudinal strain and torsion in the infarct group. At the same time the transmural helix angle (HA) gradient was steeper in the chronic follow-up and a correlation between longitudinal strain and transmural HA gradient was detected (r = 0.59 with p < 0.05). Comparing non-gadolinium enhanced data T2 mapping showed the largest relative change between infarct and remote during the acute stage (+ 33 ± 4% day 6, with p = 0.013 T2 vs. MD, p = 0.009 T2 vs. FA and p = 0.01 T2 vs. T1) while FA exhibited the largest relative change between infarct and remote during the chronic follow-up (+ 31 ± 2% week 5, with p = N.S. FA vs. MD, p = 0.03 FA vs. T2 and p = 0.003 FA vs. T1). Overall, diffusion parameters provided a higher contrast (> 23% for MD and > 27% for FA) during follow-up compared to relaxometry (T1 17-18%/T2 10-20%).

CONCLUSION

During chronic follow-up after myocardial infarction, cardiac diffusion tensor imaging provides a higher sensitivity for mapping microstructural alterations when compared to non-contrast enhanced relaxometry with the added benefit of providing directional tensor information to assess remodelling of myocyte aggregate orientations, which cannot be otherwise assessed.

Septaly Oriented Mild Aortic Regurgitant Jets Negatively Influence Left Ventricular Blood Flow-Insights From 4D Flow MRI Animal Study

Objectives: Paravalvular leakage (PVL) and eccentric aortic regurgitation remain a major clinical concern in patients receiving transcatheter aortic valve replacement (TAVR), and regurgitant volume remains the main readout parameter in clinical assessment. In this work we investigate the effect of jet origin and trajectory of mild aortic regurgitation on left ventricular hemodynamics in a porcine model.

Methods: A pig model of mild aortic regurgitation/PVL was established by transcatheter piercing and dilating the non-coronary (NCC) or right coronary cusp (RCC) of the aortic valve close to the valve annulus. The interaction between regurgitant blood and LV hemodynamics was assessed by 4D flow cardiovascular MRI.

Results: Six RCC, six NCC, and two control animals were included in the study and with one dropout in the NCC group, the success rate of model creation was 93%. Regurgitant jets originating from NCC were directed along the ventricular side of the anterior mitral leaflet and integrated well into the diastolic vortex forming in the left ventricular outflow tract. However, jets from the RCC were orientated along the septum colliding with flow within the vortex, and progressing down to the apex. As a consequence, the presence as well as the area of the vortex was reduced at the site of impact compared to the NCC group. Impairment of vortex formation was localized to the area of impact and not the entire vortex ring. Blood from the NCC jet was largely ejected during the following systole, whereas ejection of large portion of RCC blood was protracted.

Conclusions: Even for mild regurgitation, origin and trajectory of the regurgitant jet does cause a different effect on LV hemodynamics. Septaly oriented jets originating from RCC collide with the diastolic vortex, reduce its size, and reach the apical region of the left ventricle where blood resides extendedly. Hence, RCC jets display hemodynamic features which may have a potential negative impact on the long-term burden to the heart.

Emerging Techniques in Cardiac Magnetic Resonance Imaging

Cardiovascular disease is the leading cause of death and a significant contributor of health care costs. Noninvasive imaging plays an essential role in the management of patients with cardiovascular disease. Cardiac magnetic resonance (MR) can noninvasively assess heart and vascular abnormalities, including biventricular structure/function, blood hemodynamics, myocardial tissue composition, microstructure, perfusion, metabolism, coronary microvascular function, and aortic distensibility/stiffness. Its ability to characterize myocardial tissue composition is unique among alternative imaging modalities in cardiovascular disease. Significant growth in cardiac MR utilization, particularly in Europe in the last decade, has laid the necessary clinical groundwork to position cardiac MR as an important imaging modality in the workup of patients with cardiovascular disease. Although lack of availability, limited training, physician hesitation, and reimbursement issues have hampered widespread clinical adoption of cardiac MR in the United States, growing clinical evidence will ultimately overcome these challenges. Advances in cardiac MR techniques, particularly faster image acquisition, quantitative myocardial tissue characterization, and image analysis have been critical to its growth. In this review article, we discuss recent advances in established and emerging cardiac MR techniques that are expected to strengthen its capability in managing patients with cardiovascular disease. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

Diffuse myocardial fibrosis precedes impairment of myocardial strain in patients with systemic sclerosis

Funding Acknowledgements

Type of funding sources: None.

Background - Myocardial involvement is common in patients with systemic sclerosis (SSc) and causes myocardial fibrosis and subtle ventricular dysfunction. However, the temporal onset of myocardial involvement during the progression of the disease is yet unknown.

Purpose - To investigated the presence of subclinical functional impairment and diffuse myocardial fibrosis in patients with very early diagnosis of SSc (VEDOSS) and to compared the findings to patients with established SSc and healthy controls.

Methods - 110 SSc patients (86 with established SSc and 24 with VEDOSS) and 15 healthy controls were prospectively recruited. The study subjects underwent cardiovascular magnetic resonance on a clinical 1.5T system. Pre- and post-contrast T1 mapping was performed using a MOLLI (Modified Look-Locker Inversion Recovery) sequence. For extracellular volume (ECV) measurements, a single bolus protocol with image acquisition 15-20 min. post-contrast injection was used. For the assessment of subtle functional impairment, global longitudinal (GLS) and circumferential (GCS) myocardial strain were evaluated.

Results - Native T1 values and ECV were elevated in VEDOSS and in patients with established SSc compared to controls (p &lt; 0.001; Figure 1 A & B). GLS was similar in VEDOSS and controls but significantly reduced in patients with established SSc (p &lt; 0.001; Figure 1 C). GCS was similar over all groups (p = 0.88). Patients with clinical evidence of pulmonary or gastrointestinal involvement had higher ECV or T1 values, respectively. Patients with clinical signs of cardiac involvement had lower absolute GLS. SSc subtype, classification or disease duration were not associated with the extent of myocardial fibrosis or impaired strain.

Conclusion - Subclinical myocardial involvement first manifests as diffuse myocardial fibrosis identified by expansion of ECV and increased native T1 in VEDOSS patients while subtle functional impairment as measured by GLS only occurs in established SSc. No single clinical feature of SSc shows a strong association with subtle myocardial involvement.

The relationship between myocardial microstructure and strain in chronic infarcts, assessed using diffusion tensor imaging and feature tracking

Funding Acknowledgements

Type of funding sources: Foundation. Main funding source(s): British Heart Foundation

Background

Cardiac diffusion tensor imaging (cDTI) is a novel technique for the non-invasive assessment of myocardial microstructure.  It allows in-vivo characterisation of microstructural changes post myocardial infarction (MI). Previously published evidence shows significant loss of sheetlet orientation as derived by cDTI secondary eigenvector (E2A), and loss of subendocardial cardiomyocytes derived by reductions in the proportions of myocytes with right-handed orientation (RHM) on helix angle (HA) maps. The assessment of myocardial strain by feature tracking (FT) allows the measurement of radial strain (RS), thought to be driven by the dynamic reorientation of laminar sheetlets, and longitudinal strain (LS), which is thought to relate to subendocardial function. We sought to explore the relationship between the strain and cDTI parameters in patients at 3 months following ST-elevation MI (STEMI).

Methods

Twenty five STEMI patients (M:F = 18:7, mean age 58 ± 9) underwent 3T CMR scan (mean interval 106 ± 17 days) with the following protocol: second order motion compensated (M2), free-breathing spin echo DTI (3 slices, 18 diffusion directions at b-values 100s/mm2, 200s/mm2 and 500s/mm2, acquired resolution was 2.20*2.27*8mm3; cine gradient echo and Late Gadolinium Enhancement (LGE) imaging. HA maps were described by dividing values into left-handed HA (LHM, -90&lt; HA &lt; -30), circumferential HA (CM, -30° &lt; HA &lt; 30°), and right-handed HA (RHM, 30° &lt; HA &lt; 90°) and reported as relative proportions. Segmental analysis were undertaken to derive: HA proportions, E2A, longitudinal strain and LGE%. Segments positive for LGE were classed as infarct segments.

Results

cDTI acquisition was successful in all patients (acquisition time 13 ± 5mins). Mean ejection fraction was 47 ± 8% with mean LGE in the infarcted segment of 57 ± 27%. Mean radial strain was 21 (95% confidence interval, 15-26). The mean E2A was 44 (95% confidence interval 41-47). There was a significant correlation between segmental radial strain and segmental E2A in infarcted segments (p &lt; 0.001, figure 1). In addition, segmental longitudinal strain correlated with the proportion of RHM on HA maps (p &lt; 0.02, figure 2).

Conclusion

Through the combined use of cDTI and FT in patients with chronic infarcts, our results show that the loss of sheetlet orientation assessed using E2A, correlates with worsening radial strain. Segments with less subendocardial cardiomyocytes, evidenced by a lower proportion of myocytes with right-handed orientation on HA maps, correlated with worse longitudinal strain. While this could potentially elucidate the complex association between myocardial microstructure and regional function, further studies are needed to define the incremental clinical value of cDTI.

Investigating the reference domain influence in personalised models of cardiac mechanics : Effect of unloaded geometry on cardiac biomechanics

A major concern in personalised models of heart mechanics is the unknown zero-pressure domain, a prerequisite for accurately predicting cardiac biomechanics. As the reference configuration cannot be captured by clinical data, studies often employ in-vivo frames which are unlikely to correspond to unloaded geometries. Alternatively, zero-pressure domain is approximated through inverse methodologies, which, however, entail assumptions pertaining to boundary conditions and material parameters. Both approaches are likely to introduce biases in estimated biomechanical properties; nevertheless, quantification of these effects is unattainable without ground-truth data. In this work, we assess the unloaded state influence on model-derived biomechanics, by employing an in-silico modelling framework relying on experimental data on porcine hearts. In-vivo images are used for model personalisation, while in-situ experiments provide a reliable approximation of the reference domain, creating a unique opportunity for a validation study. Personalised whole-cycle cardiac models are developed which employ different reference domains (image-derived, inversely estimated) and are compared against ground-truth model outcomes. Simulations are conducted with varying boundary conditions, to investigate the effect of data-derived constraints on model accuracy. Attention is given to modelling the influence of the ribcage on the epicardium, due to its close proximity to the heart in the porcine anatomy. Our results find merit in both approaches for dealing with the unknown reference domain, but also demonstrate differences in estimated biomechanical quantities such as material parameters, strains and stresses. Notably, they highlight the importance of a boundary condition accounting for the constraining influence of the ribcage, in forward and inverse biomechanical models.

A 3D personalized cardiac myocyte aggregate orientation model using MRI data-driven low-rank basis functions

Cardiac myocyte aggregate orientation has a strong impact on cardiac electrophysiology and mechanics. Studying the link between structural characteristics, strain, and stresses over the cardiac cycle and cardiac function requires a full volumetric representation of the microstructure. In this work, we exploit the structural similarity across hearts to extract a low-rank representation of predominant myocyte orientation in the left ventricle from high-resolution magnetic resonance ex-vivo cardiac diffusion tensor imaging (cDTI) in porcine hearts. We compared two reduction methods, Proper Generalized Decomposition combined with Singular Value Decomposition and Proper Orthogonal Decomposition. We demonstrate the existence of a general set of basis functions of aggregated myocyte orientation which defines a data-driven, personalizable, parametric model featuring higher flexibility than existing atlas and rule-based approaches. A more detailed representation of microstructure matching the available patient data can improve the accuracy of personalized computational models. Additionally, we approximate the myocyte orientation of one ex-vivo human heart and demonstrate the feasibility of transferring the basis functions to humans.

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.

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.

Motion and eddy current-induced signal dephasing in in vivo cardiac DTI

PURPOSE

To address motion in cardiac DWI, stimulated-echo acquisition mode (STEAM) and second-order motion-compensated spin-echo (SE) sequences have been proposed. Despite applying motion-compensation strategies, residual motion can cause misleading signal attenuation. The purpose of this study is to estimate the motion-induced error in both sequences by analysis of image phase.

METHODS

Diffusion-weighted motion-compensated SE sequences and STEAM imaging was applied in vivo with diffusion encoding along 3 orthogonal directions. A b-value range of 100 to 600 s/mm² and trigger delays of 25%, 50%, and 75% of end systole and middiastole were used. Eddy-current contributions were obtained from phantom measurements. After computation of motion-induced phase maps, the amount of signal dephasing was computed from phase gradients, and the resulting errors in diffusion tensor parameters were calculated.

RESULTS

Motion-induced dephasing from the STEAM sequence showed less dependency on the b-value and no dependency on the heart phase, whereas SE imaging performed best at 75% end systole followed by 50% end systole and middiastole. For a typical experimental setting, errors of 3.3%/3.0% mean diffusivity, 4.9%/4.8% fractional anisotropy, 2.9º/3.2º helix angulation, 0.8º/0.7º transverse angulation, and 9.9º/10.0º sheet angulation (SE/STEAM) were calculated.

CONCLUSION

Image phase contains valuable information regarding uncompensated motion and eddy currents in cardiac DTI. Although the trigger delay window for SE is narrower compared with the STEAM-based approach, imaging in both systole and diastole is feasible and both sequences perform similarly if the trigger delays are selected carefully with SE.

Characterizing cardiac involvement in amyloidosis using cardiovascular magnetic resonance diffusion tensor imaging

BACKGROUND

In-vivo cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) allows imaging of alterations of cardiac fiber architecture in diseased hearts. Cardiac amyloidosis (CA) causes myocardial infiltration of misfolded proteins with unknown consequences for myocardial microstructure. This study applied CMR DTI in CA to assess microstructural alterations and their consequences for myocardial function compared to healthy controls.

METHODS

Ten patients with CA (8 AL, 2 ATTR) and ten healthy controls were studied using a diffusion-weighed second-order motion-compensated spin-echo sequence at 1.5 T. Additionally, left ventricular morphology, ejection fraction, strain and native T1 values were obtained in all subjects. In CA patients, T1 mapping was repeated after the administration of gadolinium for extracellular volume fraction (ECV) calculation. CMR DTI analysis was performed to yield the scalar diffusion metrics mean diffusivity (MD) and fractional anisotropy (FA) as well as the characteristics of myofiber orientation including helix, transverse and E2A sheet angle (HA, TA, E2A).

RESULTS

MD and FA were found to be significantly different between CA patients and healthy controls (MD 1.77 ± 0.17 10- 3 vs 1.41 ± 0.07 10- 3 mm2/s, p <  0.001; FA 0.25 ± 0.04 vs 0.35 ± 0.03, p <  0.001). MD demonstrated an excellent correlation with native T1 (r = 0.908, p <  0.001) while FA showed a significant correlation with ECV in the CA population (r = - 0.851, p <  0.002). HA exhibited a more circumferential orientation of myofibers in CA patients, in conjunction with a higher TA standard deviation and a higher absolute E2A sheet angle. The transmural HA slope was found to be strongly correlated with the global longitudinal strain (r = 0.921, p < 0.001).

CONCLUSION

CMR DTI reveals significant alterations of scalar diffusion metrics in CA patients versus healthy controls. Elevated MD and lower FA values indicate myocardial disarray with higher diffusion in CA that correlates well with native T1 and ECV measures. In CA patients, CMR DTI showed pronounced circumferential orientation of the myofibers, which may provide the rationale for the reduction of global longitudinal strain that occurs in amyloidosis patients. Accordingly, CMR DTI captures specific features of amyloid infiltration, which provides a deeper understanding of the microstructural consequences of CA.

Motion-Induced Signal Loss in In Vivo Cardiac Diffusion-Weighted Imaging

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:319-320.

Left ventricular blood flow patterns at rest and under dobutamine stress in healthy pigs

Intracardiac blood flow patterns are affected by the morphology of cardiac structures and are set up to support the heart's pump function. Exercise affects contractility and chamber size as well as pre- and afterload. The aim of this study was to test the feasibility of four-dimensional phase contrast cardiovascular MRI under pharmacological stress and to study left ventricular blood flow under stress. 4D flow data were successfully acquired and analysed in 12 animals. During dobutamine infusion, heart rate and ejection fraction increased (82 ± 5 bpm versus 124 ± 3 bpm/46 ± 9% versus 65 ± 7%; both p < 0.05). A decrease in left ventricular end-diastolic volume (72 ± 14 mL versus 55 ± 8 mL; p < 0.05) and end-systolic volume (40 ± 15 mL versus 19 ± 6 mL; p < 0.05) but no change in stroke volume were observed. Trans-mitral diastolic inflow velocity increased under dobutamine and the trajectory of inflowing blood was directed towards the anterior septum with increased inflow angle (26 ± 5°) when compared with controls (15 ± 2°). In 5/6 animals undergoing stress diastolic vortices developed later, and in 3/6 animals vortices collapsed earlier with significantly smaller cross-sectional area during diastole. The vorticity index was not affected. Under the stress condition direct flow (% ejection within the next heart beat) increased from 43 ± 6% to 53 ± 8%. 4D MRI blood flow acquisition and analysis are feasible in pig hearts under dobutamine-induced stress. Flow patterns characterized by high blood velocity and antero-septally oriented diastolic inflow as well as decreased ventricular volumes are unfavourable conditions for diastolic vortex development under pharmacological stress, and cardiac output is increased by a rise in heart rate and directly ejected left ventricular blood volume.

Enhancing intravoxel incoherent motion parameter mapping in the brain using k-b PCA

Intravoxel incoherent motion (IVIM) imaging of diffusion and perfusion parameters in the brain using parallel imaging suffers from local noise amplification. To address the issue, signal correlations in space and along the diffusion encoding dimension are exploited jointly using a constrained image reconstruction approach. IVIM imaging was performed on a clinical 3 T MR system with diffusion weighting along six gradient directions and 16 b-values encoded per direction across a range of 0-900 s/mm² . Data were collected in 11 subjects, retrospectively undersampled in k-space with net factors ranging from 2 to 6 and reconstructed using CG-SENSE and the proposed k-b PCA approach. Results of k-b PCA and CG-SENSE from retrospectively undersampled data were compared with those from the fully sampled reference. In addition, prospective single-shot k-b undersampling was implemented and data were acquired in five additional volunteers. IVIM parameter maps were derived using a segmented least-squares method. The proposed k-b PCA method outperformed CG-SENSE in terms of reconstruction errors for effective undersampling factors of 3 and beyond. Undersampling artifacts were effectively removed with k-b PCA up to sixfold undersampling. At net sixfold undersampling, relative errors (compared with the fully sampled reference) of image magnitude and IVIM parameters (D, f and D^* ) were (median ± interquartile range): 3.5 ± 3.7 versus 25.3 ± 25.8%, 2.7 ± 3.6 versus 14.2 ± 20.4%, 15.1 ± 26.1 versus 96.6 ± 67.4% and 14.8 ± 26.6 versus 100 ± 195.1% for k-b PCA versus CG-SENSE, respectively. Acquisition with sixfold prospective undersampling yielded average IVIM parameters in the brain of 0.79 ± 0.18 × 10^-3  mm² /s for D, 7.35 ± 7.27% for f and 7.11 ± 2.39 × 10^-3  mm² /s for D^* . Constrained reconstruction using k-b PCA improves IVIM parameter mapping from undersampled data when compared with CG-SENSE reconstruction. Prospectively undersampled single-shot echo planar imaging acquisition was successfully employed using k-b PCA, demonstrating a reduction of image artifacts and noise relative to parallel imaging.

Quantitative myocardial first-pass cardiovascular magnetic resonance perfusion imaging using hyperpolarized [1-13C] pyruvate

BACKGROUND

The feasibility of absolute myocardial blood flow quantification and suitability of hyperpolarized [1-13C] pyruvate as contrast agent for first-pass cardiovascular magnetic resonance (CMR) perfusion measurements are investigated with simulations and demonstrated in vivo in a swine model.

METHODS

A versatile simulation framework for hyperpolarized CMR subject to physical, physiological and technical constraints was developed and applied to investigate experimental conditions for accurate perfusion CMR with hyperpolarized [1-13C] pyruvate. Absolute and semi-quantitative perfusion indices were analyzed with respect to experimental parameter variations and different signal-to-noise ratio (SNR) levels. Absolute myocardial blood flow quantification was implemented with an iterative deconvolution approach based on Fermi functions. To demonstrate in vivo feasibility, velocity-selective excitation with an echo-planar imaging readout was used to acquire dynamic myocardial stress perfusion images in four healthy swine. Arterial input functions were extracted from an additional image slice with conventional excitation that was acquired within the same heartbeat.

RESULTS

Simulations suggest that obtainable SNR and B0 inhomogeneity in vivo are sufficient for the determination of absolute and semi-quantitative perfusion with ≤25% error. It is shown that for expected metabolic conversion rates, metabolic conversion of pyruvate can be neglected over the short duration of acquisition in first-pass perfusion CMR. In vivo measurements suggest that absolute myocardial blood flow quantification using hyperpolarized [1-13C] pyruvate is feasible with an intra-myocardial variability comparable to semi-quantitative perfusion indices.

CONCLUSION

The feasibility of quantitative hyperpolarized first-pass perfusion CMR using [1-13C] pyruvate has been investigated in simulations and demonstrated in swine. Using an approved and metabolically active compound is envisioned to increase the value of hyperpolarized perfusion CMR in patients.

Equilibrated warping: Finite element image registration with finite strain equilibrium gap regularization

In this paper, we propose a novel continuum finite strain formulation of the equilibrium gap regularization for image registration. The equilibrium gap regularization essentially penalizes any deviation from the solution of a hyperelastic body in equilibrium with arbitrary loads prescribed at the boundary. It thus represents a regularization with strong mechanical basis, especially suited for cardiac image analysis. We describe the consistent linearization and discretization of the regularized image registration problem, in the framework of the finite elements method. The method is implemented using FEniCS & VTK, and distributed as a freely available python library. We show that the equilibrated warping method is effective and robust: regularization strength and image noise have minimal impact on motion tracking, especially when compared to strain-based regularization methods such as hyperelastic warping. We also show that equilibrated warping is able to extract main deformation features on both tagged and untagged cardiac magnetic resonance images.

Myocardial Scar Delineation Using Diffusion Tensor Magnetic Resonance Tractography

Background

Late gadolinium enhancement (LGE) is the current standard for myocardial scar delineation. In this study, we introduce the tractographic propagation angle (PA), a metric of myofiber curvature (degrees/unit distance) derived from diffusion tensor imaging (DTI), and compare its use to LGE and invasive scar assessment by endocardial voltage mapping.

Methods and Results

DTI was performed on 7 healthy human volunteers, 5 patients with myocardial infarction, 6 normal mice, and 7 mice with myocardial infarction. LGE to delineate the infarct and border zones was performed with a 2‐dimensional inversion recovery gradient‐echo sequence. Ex vivo DTI was performed on 5 normal human and 5 normal sheep hearts. Endocardial electroanatomic mapping and subsequent ex vivo DTI was performed on 5 infarcted sheep hearts. PA in the normal human hearts varied smoothly and was generally <4. The mean PA in the infarct zone was significantly elevated (10.34±1.02 versus 4.05±0.45, P<0.05). Regions with a PA ≤4 consistently had a bipolar voltage ≥1.5 mV, whereas those with PA values between 4 and 10 had voltages between 0.5 and 1.5 mV. A PA threshold >4 was the most accurate DTI‐derived measure of infarct size and demonstrated the greatest correlation with LGE (r=0.95).

Conclusions

We found a strong correlation between infarct size by PA and LGE in both mice and humans. There was also an inverse relationship between PA values and endocardial voltage. The use of PA may enable myocardial scar delineation and characterization of arrhythmogenic substrate without the need for exogenous contrast agents.

Bayesian intravoxel incoherent motion parameter mapping in the human heart

BACKGROUND

Intravoxel incoherent motion (IVIM) imaging of diffusion and perfusion in the heart suffers from high parameter estimation error. The purpose of this work is to improve cardiac IVIM parameter mapping using Bayesian inference.

METHODS

A second-order motion-compensated diffusion weighted spin-echo sequence with navigator-based slice tracking was implemented to collect cardiac IVIM data in early systole in eight healthy subjects on a clinical 1.5 T CMR system. IVIM data were encoded along six gradient optimized directions with b-values of 0-300 s/mm2. Subjects were scanned twice in two scan sessions one week apart to assess intra-subject reproducibility. Bayesian shrinkage prior (BSP) inference was implemented to determine IVIM parameters (diffusion D, perfusion fraction F and pseudo-diffusion D*). Results were compared to least-squares (LSQ) parameter estimation. Signal-to-noise ratio (SNR) requirements for a given fitting error were assessed for the two methods using simulated data. Reproducibility analysis of parameter estimation in-vivo using BSP and LSQ was performed.

RESULTS

BSP resulted in reduced SNR requirements when compared to LSQ in simulations. In-vivo, BSP analysis yielded IVIM parameter maps with smaller intra-myocardial variability and higher estimation certainty relative to LSQ. Mean IVIM parameter estimates in eight healthy subjects were (LSQ/BSP): 1.63 ± 0.28/1.51 ± 0.14·10-3 mm2/s for D, 13.13 ± 19.81/13.11 ± 5.95% for F and 201.45 ± 313.23/13.11 ± 14.53·10-3 mm2/s for D ∗. Parameter variation across all volunteers and measurements was lower with BSP compared to LSQ (coefficient of variation BSP vs. LSQ: 9% vs. 17% for D, 45% vs. 151% for F and 111% vs. 155% for D ∗). In addition, reproducibility of the IVIM parameter estimates was higher with BSP compared to LSQ (Bland-Altman coefficients of repeatability BSP vs. LSQ: 0.21 vs. 0.26·10-3 mm2/s for D, 5.55 vs. 6.91% for F and 15.06 vs. 422.80·10-3 mm2/s for D*).

CONCLUSION

Robust free-breathing cardiac IVIM data acquisition in early systole is possible with the proposed method. BSP analysis yields improved IVIM parameter maps relative to conventional LSQ fitting with fewer outliers, improved estimation certainty and higher reproducibility. IVIM parameter mapping holds promise for myocardial perfusion measurements without the need for contrast agents.

Hyperpolarized 13C urea myocardial first-pass perfusion imaging using velocity-selective excitation

BACKGROUND

A velocity-selective binomial excitation scheme for myocardial first-pass perfusion measurements with hyperpolarized 13C substrates, which preserves bolus magnetization inside the blood pool, is presented. The proposed method is evaluated against gadolinium-enhanced 1H measurements in-vivo.

METHODS

The proposed excitation with an echo-planar imaging readout was implemented on a clinical CMR system. Dynamic myocardial stress perfusion images were acquired in six healthy pigs after bolus injection of hyperpolarized 13C urea with the velocity-selective vs. conventional excitation, as well as standard 1H gadolinium-enhanced images. Signal-to-noise, contrast-to-noise (CNR) and homogeneity of semi-quantitative perfusion measures were compared between methods based on first-pass signal-intensity time curves extracted from a mid-ventricular slice. Diagnostic feasibility is demonstrated in a case of septal infarction.

RESULTS

Velocity-selective excitation provides over three-fold reduction in blood pool signal with a two-fold increase in myocardial CNR. Extracted first-pass perfusion curves reveal a significantly reduced variability of semi-quantitative first-pass perfusion measures (12-20%) for velocity-selective excitation compared to conventional excitation (28-93%), comparable to that of reference 1H gadolinium data (9-15%). Overall image quality appears comparable between the velocity-selective hyperpolarized and gadolinium-enhanced imaging.

CONCLUSION

The feasibility of hyperpolarized 13C first-pass perfusion CMR has been demonstrated in swine. Comparison with reference 1H gadolinium data revealed sufficient data quality and indicates the potential of hyperpolarized perfusion imaging for human applications.

Studying Dynamic Myofiber Aggregate Reorientation in Dilated Cardiomyopathy Using In Vivo Magnetic Resonance Diffusion Tensor Imaging

Supplemental Digital Content is available in the text.

Background—

The objective of this study is to assess the dynamic alterations of myocardial microstructure and strain between diastole and systole in patients with dilated cardiomyopathy relative to healthy controls using the magnetic resonance diffusion tensor imaging, myocardial tagging, and biomechanical modeling.

Methods and Results—

Dual heart-phase diffusion tensor imaging was successfully performed in 9 patients and 9 controls. Tagging data were acquired for the diffusion tensor strain correction and cardiac motion analysis. Mean diffusivity, fractional anisotropy, and myocyte aggregate orientations were compared between both cohorts. Cardiac function was assessed by left ventricular ejection fraction, torsion, and strain. Computational modeling was used to study the impact of cardiac shape on fiber reorientation and how fiber orientations affect strain. In patients with dilated cardiomyopathy, a more longitudinal orientation of diastolic myofiber aggregates was measured compared with controls. Although a significant steepening of helix angles (HAs) during contraction was found in the controls, consistent change in HAs during contraction was absent in patients. Left ventricular ejection fraction, cardiac torsion, and strain were significantly lower in the patients compared with controls. Computational modeling revealed that the dilated heart results in reduced HA changes compared with a normal heart. Reduced torsion was found to be exacerbated by steeper HAs.

Conclusions—

Diffusion tensor imaging revealed reduced reorientation of myofiber aggregates during cardiac contraction in patients with dilated cardiomyopathy relative to controls. Left ventricular remodeling seems to be an important factor in the changes to myocyte orientation. Steeper HAs are coupled with a worsening in strain and torsion. Overall, the findings provide new insights into the structural alterations in patients with dilated cardiomyopathy.

Second-order motion compensated PRESS for cardiac spectroscopy

PURPOSE

Second-order motion compensation for point-resolved spectroscopy (PRESS) is proposed to allow for robust single-voxel cardiac spectroscopy throughout the entire cardiac cycle and at various heart rates.

METHODS

Bipolar FID spoiling gradient pairs compensating for first and second-order motion were designed and implemented into a cardiac-triggered PRESS sequence on a clinical MR system. A numerical three-dimensional model of cardiac motion was used to optimize and validate the gradient waveforms. In vivo measurements in healthy volunteers were obtained to assess the signal-to-noise ratio (SNR) and triglyceride-to-water ratio (TG/W). SNR gains and variability of TG/W of the proposed approach were evaluated against a conventional PRESS sequence with optimized gradients.

RESULTS

The proposed sequence increases the mean SNR by 32% (W) and 23% (TG) on average with significantly lower variability for different trigger delays. The variability of TG/W quantification over the cardiac cycle is significantly decreased with second-order motion compensated PRESS when compared with conventional PRESS with reduced-spoiler gradients (coefficient of variation: 0.1 ± 0.02 versus 0.37 ± 0.26).

CONCLUSION

Second-order motion compensated PRESS effectively reduces cardiac motion-induced signal degradation during FID spoiling, providing higher SNR and less variability for TG/W quantification. The sequence is considered promising to assess the TG/W modulation during various interventions including pharmacologically induced stress. Magn Reson Med 77:57-64, 2017. © 2016 Wiley Periodicals, Inc.

Spin echo versus stimulated echo diffusion tensor imaging of the in vivo human heart

Purpose

To compare signal‐to‐noise ratio (SNR) efficiency and diffusion tensor metrics of cardiac diffusion tensor mapping using acceleration‐compensated spin‐echo (SE) and stimulated echo acquisition mode (STEAM) imaging.

Methods

Diffusion weighted SE and STEAM sequences were implemented on a clinical 1.5 Tesla MR system. The SNR efficiency of SE and STEAM was measured (b = 50–450 s/mm²) in isotropic agar, anisotropic diffusion phantoms and the in vivo human heart. Diffusion tensor analysis was performed on mean diffusivity, fractional anisotropy, helix and transverse angles.

Results

In the isotropic phantom, the ratio of SNR efficiency for SE versus STEAM, SNR_t(SE/STEAM), was 2.84 ± 0.08 for all tested b‐values. In the anisotropic diffusion phantom the ratio decreased from 2.75 ± 0.05 to 2.20 ± 0.13 with increasing b‐value, similar to the in vivo decrease from 2.91 ± 0.43 to 2.30 ± 0.30. Diffusion tensor analysis revealed reduced deviation of helix angles from a linear transmural model and reduced transverse angle standard deviation for SE compared with STEAM. Mean diffusivity and fractional anisotropy were measured to be statistically different (P < 0.001) between SE and STEAM.

Conclusion

Cardiac DTI using motion‐compensated SE yields a 2.3–2.9× increase in SNR efficiency relative to STEAM and improved accuracy of tensor metrics. The SE method hence presents an attractive alternative to STEAM based approaches. Magn Reson Med 76:862–872, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Second-order motion-compensated spin echo diffusion tensor imaging of the human heart

PURPOSE

Myocardial microstructure has been challenging to probe in vivo. Spin echo-based diffusion-weighted sequences allow for single-shot acquisitions but are highly sensitive to cardiac motion. In this study, the use of second-order motion-compensated diffusion encoding was compared with first-order motion-compensated diffusion-weighted imaging during systolic contraction of the heart.

METHODS

First- and second-order motion-compensated diffusion encoding gradients were incorporated into a triggered single-shot spin echo sequence. The effect of contractile motion on the apparent diffusion coefficients and tensor orientations was investigated in vivo from basal to apical level of the heart.

RESULTS

Second-order motion compensation was found to increase the range of systolic trigger delays from 30%-55% to 15%-77% peak systole at the apex and from 25%-50% to 15%-79% peak systole at the base. Diffusion tensor analysis yielded more physiological transmural distributions when using second-order motion-compensated diffusion tensor imaging.

CONCLUSION

Higher-order motion-compensated diffusion encoding decreases the sensitivity to cardiac motion, thereby enabling cardiac DTI over a wider range of time points during systolic contraction of the heart.

High-resolution diffusion tensor imaging of the human kidneys using a free-breathing, multi-slice, targeted field of view approach

Fractional anisotropy (FA) obtained by diffusion tensor imaging (DTI) can be used to image the kidneys without any contrast media. FA of the medulla has been shown to correlate with kidney function. It is expected that higher spatial resolution would improve the depiction of small structures within the kidney. However, the achievement of high spatial resolution in renal DTI remains challenging as a result of respiratory motion and susceptibility to diffusion imaging artefacts. In this study, a targeted field of view (TFOV) method was used to obtain high-resolution FA maps and colour-coded diffusion tensor orientations, together with measures of the medullary and cortical FA, in 12 healthy subjects. Subjects were scanned with two implementations (dual and single kidney) of a TFOV DTI method. DTI scans were performed during free breathing with a navigator-triggered sequence. Results showed high consistency in the greyscale FA, colour-coded FA and diffusion tensors across subjects and between dual- and single-kidney scans, which have in-plane voxel sizes of 2 × 2 mm(2) and 1.2 × 1.2 mm(2) , respectively. The ability to acquire multiple contiguous slices allowed the medulla and cortical FA to be quantified over the entire kidney volume. The mean medulla and cortical FA values were 0.38 ± 0.017 and 0.21 ± 0.019, respectively, for the dual-kidney scan, and 0.35 ± 0.032 and 0.20 ± 0.014, respectively, for the single-kidney scan. The mean FA between the medulla and cortex was significantly different (p < 0.001) for both dual- and single-kidney implementations. High-spatial-resolution DTI shows promise for improving the characterization and non-invasive assessment of kidney function.

Dual-phase cardiac diffusion tensor imaging with strain correction

Purpose

In this work we present a dual-phase diffusion tensor imaging (DTI) technique that incorporates a correction scheme for the cardiac material strain, based on 3D myocardial tagging.

Methods

In vivo dual-phase cardiac DTI with a stimulated echo approach and 3D tagging was performed in 10 healthy volunteers. The time course of material strain was estimated from the tagging data and used to correct for strain effects in the diffusion weighted acquisition. Mean diffusivity, fractional anisotropy, helix, transverse and sheet angles were calculated and compared between systole and diastole, with and without strain correction. Data acquired at the systolic sweet spot, where the effects of strain are eliminated, served as a reference.

Results

The impact of strain correction on helix angle was small. However, large differences were observed in the transverse and sheet angle values, with and without strain correction. The standard deviation of systolic transverse angles was significantly reduced from 35.9±3.9° to 27.8°±3.5° (p<0.001) upon strain-correction indicating more coherent fiber tracks after correction. Myocyte aggregate structure was aligned more longitudinally in systole compared to diastole as reflected by an increased transmural range of helix angles (71.8°±3.9° systole vs. 55.6°±5.6°, p<0.001 diastole). While diastolic sheet angle histograms had dominant counts at high sheet angle values, systolic histograms showed lower sheet angle values indicating a reorientation of myocyte sheets during contraction.

Conclusion

An approach for dual-phase cardiac DTI with correction for material strain has been successfully implemented. This technique allows assessing dynamic changes in myofiber architecture between systole and diastole, and emphasizes the need for strain correction when sheet architecture in the heart is imaged with a stimulated echo approach.

Characterization and correction of eddy-current artifacts in unipolar and bipolar diffusion sequences using magnetic field monitoring

Diffusion tensor imaging (DTI) of moving organs is gaining increasing attention but robust performance requires sequence modifications and dedicated correction methods to account for system imperfections. In this study, eddy currents in the "unipolar" Stejskal-Tanner and the velocity-compensated "bipolar" spin-echo diffusion sequences were investigated and corrected for using a magnetic field monitoring approach in combination with higher-order image reconstruction. From the field-camera measurements, increased levels of second-order eddy currents were quantified in the unipolar sequence relative to the bipolar diffusion sequence while zeroth and linear orders were found to be similar between both sequences. Second-order image reconstruction based on field-monitoring data resulted in reduced spatial misalignment artifacts and residual displacements of less than 0.43 mm and 0.29 mm (in the unipolar and bipolar sequences, respectively) after second-order eddy-current correction. Results demonstrate the need for second-order correction in unipolar encoding schemes but also show that bipolar sequences benefit from second-order reconstruction to correct for incomplete intrinsic cancellation of eddy-currents.