Effect of flow-encoding strength on intravoxel incoherent motion in the liver

PURPOSE

To study the impact of variable flow-encoding strength on intravoxel incoherent motion (IVIM) liver imaging of diffusion and perfusion.

THEORY

Signal attenuation in DWI arises from (1) intravoxel microvascular blood flow, which depends on the flow-encoding strength α (first gradient moment) of the diffusion-encoding waveform, and (2) intravoxel spin diffusion, which depends on the b-value of the diffusion-encoding gradient waveforms α and b-value. Both are linked to the diffusion-encoding gradient waveform and conventionally are not independently controlled.

METHODS

In this work a convex optimization framework was used to generate gradient waveforms with independent α and b-value. Thirty-six unique α and b-value sample points from 5 different gradient waveforms were used to reconstruct perfusion fraction (f), coefficient of diffusion (D), and blood velocity standard deviation (V_b ) maps using a recently proposed IVIM model. Faster acquisition strategies were evaluated with 1000 random subsampling strategies of 16, 8, and 4 α and b-value. Among the subsampled reconstructions, the sampling schemes that minimized the difference with the fully sampled reconstruction were reported.

RESULTS

Healthy volunteers (N = 9) were imaged on a 3T scanner. Liver perfusion and diffusion estimates using the fully sampled IVIM method were f = 0.19 ± 0.06, D = 1.15 ± 0.15 × 10^-3 mm² /s, and V_b = 5.22 ± 3.86 mm/s. No statistical differences were found between the fully sampled and 2-times undersampled reconstruction (f = 0.2 ± 0.07, D = 1.19 ± 0.15 × 10^-3 mm² /s, V_b = 5.79 ± 3.43 mm/s); 4-times undersampled (f = 0.2 ± 0.06, D = 1.15 ± 0.17 × 10^-3 mm² /s, V_b = 4.66 ± 3.61 mm/s), or 8-times undersampled ( f = 0.2 ± 0.06, D = 1.23 ± 0.22 × 10^-3 mm² /s, V_b = 4.99 ± 3.82 mm/s) approaches.

CONCLUSION

We demonstrate the IVIM signal's dependence on the b-value, the diffusion-encoding time and the flow-encoding strength and observe in vivo the ballistic regime signature of microperfusion in the liver. This work also demonstrates that using an IVIM model and sampling scheme matched to the ballistic regime, pixel-wise IVIM parameter maps are possible when sampling as few as 4 IVIM signals.

Microstructural Infarct Border Zone Remodeling in the Post-infarct Swine Heart Measured by Diffusion Tensor MRI

Introduction: Computational models of the heart increasingly require detailed microstructural information to capture the impact of tissue remodeling on cardiac electromechanics in, for example, hearts with myocardial infarctions. Myocardial infarctions are surrounded by the infarct border zone (BZ), which is a site of electromechanical property transition. Magnetic resonance imaging (MRI) is an emerging method for characterizing microstructural remodeling and focal myocardial infarcts and the BZ can be identified with late gadolinium enhanced (LGE) MRI. Microstructural remodeling within the BZ, however, remains poorly characterized by MRI due, in part, to the fact that LGE and DT-MRI are not always available for the same heart. Diffusion tensor MRI (DT-MRI) can evaluate microstructural remodeling by quantifying the DT apparent diffusion coefficient (ADC, increased with decreased cellularity), fractional anisotropy (FA, decreased with increased fibrosis), and tissue mode (decreased with increased fiber disarray). The purpose of this work was to use LGE MRI in post-infarct porcine hearts (N = 7) to segment remote, BZ, and infarcted myocardium, thereby providing a basis to quantify microstructural remodeling in the BZ and infarcted regions using co-registered DT-MRI. Methods: Chronic porcine infarcts were created by balloon occlusion of the LCx. 6-8 weeks post-infarction, MRI contrast was administered, and the heart was potassium arrested, excised, and imaged with LGE MRI (0.33 × 0.33 × 0.33 mm) and co-registered DT-MRI (1 × 1 × 3 mm). Myocardium was segmented as remote, BZ, or infarct by LGE signal intensity thresholds. DT invariants were used to evaluate microstructural remodeling by quantifying ADC, FA, and tissue mode. Results: The BZ significantly remodeled compared to both infarct and remote myocardium. BZ demonstrated a significant decrease in cellularity (increased ADC), significant decrease in tissue organization (decreased FA), and a significant increase in fiber disarray (decreased tissue mode) relative to remote myocardium (all p < 0.05). Microstructural remodeling in the infarct was similar, but significantly larger in magnitude (all p < 0.05). Conclusion: DT-MRI can identify regions of significant microstructural remodeling in the BZ that are distinct from both remote and infarcted myocardium.

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

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

Construction and Validation of Subject-Specific Biventricular Finite-Element Models of Healthy and Failing Swine Hearts From High-Resolution DT-MRI

Predictive computational modeling has revolutionized classical engineering disciplines and is in the process of transforming cardiovascular research. This is particularly relevant for investigating emergent therapies for heart failure, which remains a leading cause of death globally. The creation of subject-specific biventricular computational cardiac models has been a long-term endeavor within the biomedical engineering community. Using high resolution (0.3 × 0.3 × 0.8 mm) ex vivo data, we constructed a precise fully subject-specific biventricular finite-element model of healthy and failing swine hearts. Each model includes fully subject-specific geometries, myofiber architecture and, in the case of the failing heart, fibrotic tissue distribution. Passive and active material properties are prescribed using hyperelastic strain energy functions that define a nearly incompressible, orthotropic material capable of contractile function. These materials were calibrated using a sophisticated multistep approach to match orthotropic tri-axial shear data as well as subject-specific hemodynamic ventricular targets for pressure and volume to ensure realistic cardiac function. Each mechanically beating heart is coupled with a lumped-parameter representation of the circulatory system, allowing for a closed-loop definition of cardiovascular flow. The circulatory model incorporates unidirectional fluid exchanges driven by pressure gradients of the model, which in turn are driven by the mechanically beating heart. This creates a computationally meaningful representation of the dynamic beating of the heart coupled with the circulatory system. Each model was calibrated using subject-specific experimental data and compared with independent in vivo strain data obtained from echocardiography. Our methods produced highly detailed representations of swine hearts that function mechanically in a remarkably similar manner to the in vivo subject-specific strains on a global and regional comparison. The degree of subject-specificity included in the models represents a milestone for modeling efforts that captures realism of the whole heart. This study establishes a foundation for future computational studies that can apply these validated methods to advance cardiac mechanics research.

TIME RESOLVED DISPLACEMENT-BASED REGISTRATION OF IN VIVO CDTI CARDIOMYOCYTE ORIENTATIONS

In vivo cardiac microstructure acquired using cardiac diffusion tensor imaging (cDTI) is a critical component of patient-specific models of cardiac electrophysiology and mechanics. In order to limit bulk motion artifacts and acquisition time, cDTI microstructural data is acquired at a single cardiac phase necessitating registration to the reference configuration on which the patient-specific computational models are based. Herein, we propose a method to register subject-specific microstructural data to an arbitrary cardiac phase using measured cardiac displacements. We validate our approach using a subject-specific computational phantom based on data from human subjects. Compared to a geometry-based non-rigid registration method, the displacement-based registration leads to improved accuracy (less than 1° versus 10° average median error in cardiomyocyte angular differences) and tighter confidence interval (3° versus 65° average upper limit of the 95% confidence interval).

Velocity reconstruction with nonconvex optimization for low-velocity-encoding phase-contrast MRI

PURPOSE

To introduce and demonstrate a nonconvex optimization method for reconstructing velocity data from low-velocity-encoding (V_enc ) phase-contrast MRI data.

THEORY AND METHODS

Solving for velocity values from phase-contrast MRI data was formulated as a nonconvex optimization problem. Weighting was added to account for intravoxel dephasing, and a Laplacian-based regularization was used to account for residual velocity aliasing. The reconstruction was tested with two low-V_enc schemes: dual-V_enc and a multidirectional high-moment encoding. The reconstruction method was tested in a digital simulation, in flow phantoms, and in healthy volunteers (N = 5).

RESULTS

The nonconvex-optimization reconstruction velocity error was lower than the conventional reconstruction in simulations (4.6 versus 3.0 cm/s for multidirectional high moment, 8.3 versus 3.8 cm/s for dual-V_enc ) and in flow phantoms (23.9 versus 5.9 cm/s for multidirectional high moment, 15.2 versus 6.4 cm/s for dual-V_enc ). Qualitative assessment of velocity fields in all experiments, including healthy volunteers, showed decreased apparent noise in the velocity fields and fewer phase wraps. No additional velocity bias in measured velocities was seen in volunteers with the proposed method.

CONCLUSIONS

The proposed nonconvex-optimization reconstruction method incorporates additional information to solve for velocities when using any type of low-V_enc (high-moment) acquisition. The method reduces the amount of residual phase aliasing, and decreases velocity errors that result from intravoxel dephasing. These improvements allow for more robust acquisitions, and for V_enc to be lowered 2 to 4 times relative to conventional acquisitions, thereby increasing the velocity-to-noise ratio. Magn Reson Med, 2017. © 2017 International Society for Magnetic Resonance in Medicine. Magn Reson Med 80:42-52, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Method for the unique identification of hyperelastic material properties using full-field measures. Application to the passive myocardium material response

Quantitative measurement of the material properties (eg, stiffness) of biological tissues is poised to become a powerful diagnostic tool. There are currently several methods in the literature to estimating material stiffness, and we extend this work by formulating a framework that leads to uniquely identified material properties. We design an approach to work with full-field displacement data-ie, we assume the displacement field due to the applied forces is known both on the boundaries and also within the interior of the body of interest-and seek stiffness parameters that lead to balanced internal and external forces in a model. For in vivo applications, the displacement data can be acquired clinically using magnetic resonance imaging while the forces may be computed from pressure measurements, eg, through catheterization. We outline a set of conditions under which the least-square force error objective function is convex, yielding uniquely identified material properties. An important component of our framework is a new numerical strategy to formulate polyconvex material energy laws that are linear in the material properties and provide one optimal description of the available experimental data. An outcome of our approach is the analysis of the reliability of the identified material properties, even for material laws that do not admit unique property identification. Lastly, we evaluate our approach using passive myocardium experimental data at the material point and show its application to identifying myocardial stiffness with an in silico experiment modeling the passive filling of the left ventricle.

Cardiac MRI: a Translational Imaging Tool for Characterizing Anthracycline-Induced Myocardial Remodeling

Cardiovascular side effects of cancer therapeutics are the leading causes of morbidity and mortality in cancer survivors. Anthracyclines (AC) serve as the backbone of many anti-cancer treatment strategies, but dose-dependent myocardial injury limits their use. Cumulative AC exposure can disrupt the dynamic equilibrium of the myocardial microarchitecture while repeated injury and repair leads to myocyte loss, interstitial myocardial fibrosis, and impaired contractility. Although children are assumed to have greater myocardial plasticity, AC exposure at a younger age portends worse prognosis. In older patients, there is lower overall survival once they develop cardiovascular disease. Because aberrations in the myocardial architecture predispose the heart to a decline in function, early detection with sensitive imaging tools is crucial and the implications for resource utilization are substantial. As a comprehensive imaging modality, cardiac magnetic resonance (CMR) imaging is able to go beyond quantification of ejection fraction and myocardial deformation to characterize adaptive microstructural and microvascular changes that are important to myocardial tissue health. Herein, we describe CMR as an established translational imaging tool that can be used clinically to characterize AC-associated myocardial remodeling.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Microstructurally Anchored Cardiac Kinematics by Combining In Vivo DENSE MRI and cDTI

Metrics of regional myocardial function can detect the onset of cardiovascular disease, evaluate the response to therapy, and provide mechanistic insight into cardiac dysfunction. Knowledge of local myocardial microstructure is necessary to distinguish between isotropic and anisotropic contributions of local deformation and to quantify myofiber kinematics, a microstructurally anchored measure of cardiac function. Using a computational model we combine in vivo cardiac displacement and diffusion tensor data to evaluate pointwise the deformation gradient tensor and isotropic and anisotropic deformation invariants. In discussing the imaging methods and the model construction, we identify potential improvements to increase measurement accuracy. We conclude by demonstrating the applicability of our method to compute myofiber strain in five healthy volunteers.

Simultaneous measurement of T2 and apparent diffusion coefficient (T2 +ADC) in the heart with motion-compensated spin echo diffusion-weighted imaging

PURPOSE

To evaluate a technique for simultaneous quantitative T₂ and apparent diffusion coefficient (ADC) mapping in the heart (T₂ +ADC) using spin echo (SE) diffusion-weighted imaging (DWI).

THEORY AND METHODS

T₂ maps from T₂ +ADC were compared with single-echo SE in phantoms and with T₂ -prepared (T₂ -prep) balanced steady-state free precession (bSSFP) in healthy volunteers. ADC maps from T₂ +ADC were compared with conventional DWI in phantoms and in vivo. T₂ +ADC was also demonstrated in a patient with acute myocardial infarction (MI).

RESULTS

Phantom T₂ values from T₂ +ADC were closer to a single-echo SE reference than T₂ -prep bSSFP (-2.3 ± 6.0% vs 22.2 ± 16.3%; P < 0.01), and ADC values were in excellent agreement with DWI (0.28 ± 0.4%). In volunteers, myocardial T₂ values from T₂ +ADC were significantly shorter than T₂ -prep bSSFP (35.8 ± 3.1 vs 46.8 ± 3.8 ms; P < 0.01); myocardial ADC was not significantly (N.S.) different between T₂ +ADC and conventional motion-compensated DWI (1.39 ± 0.18 vs 1.38 ± 0.18 mm² /ms; P = N.S.). In the patient, T₂ and ADC were both significantly elevated in the infarct compared with remote myocardium (T₂ : 40.4 ± 7.6 vs 56.8 ± 22.0; P < 0.01; ADC: 1.47 ± 0.59 vs 1.65 ± 0.65 mm² /ms; P < 0.01).

CONCLUSION

T₂ +ADC generated coregistered, free-breathing T₂ and ADC maps in healthy volunteers and a patient with acute MI with no cost in accuracy, precision, or scan time compared with DWI. Magn Reson Med 79:654-662, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

A viscoactive constitutive modeling framework with variational updates for the myocardium

We present a constitutive modeling framework for contractile cardiac mechanics by formulating a single variational principle from which incremental stress-strain relations and kinetic rate equations for active contraction and relaxation can all be derived. The variational framework seamlessly incorporates the hyperelastic behavior of the relaxed and contracted tissue along with the rate - and length - dependent generation of contractile force. We describe a three-element, Hill-type model that unifies the active tension and active deformation approaches. As in the latter approach, we multiplicatively decompose the total deformation gradient into active and elastic parts, with the active deformation parametrizing the contractile Hill element. We adopt as internal variables the fiber, cross-fiber, and sheet normal stretch ratios. The kinetics of these internal variables are modeled via definition of a kinetic potential function derived from experimental force-velocity relations. Additionally, we account for dissipation during tissue deformation by adding a Newtonian viscous potential. To model the force activation, the kinetic equations are coupled with the calcium transient obtained from a cardiomyocyte electrophysiology model. We first analyze our model at the material point level using stress and strain versus time curves for different viscosity values. Subsequently, we couple our constitutive framework with the finite element method (FEM) and study the deformation of three-dimensional tissue slabs with varying cardiac myocyte orientation. Finally, we simulate the contraction and relaxation of an ellipsoidal left ventricular model and record common kinematic measures, such as ejection fraction, and myocardial tissue volume changes.

Terahertz Imaging of Cutaneous Edema: Correlation With Magnetic Resonance Imaging in Burn Wounds

OBJECTIVE

In vivo visualization and quantification of edema, or 'tissue swelling' following injury, remains a clinical challenge. Herein, we investigate the ability of reflective terahertz (THz) imaging to track changes in tissue water content (TWC)-the direct indicator of edema-by comparison to depth-resolved magnetic resonance imaging (MRI) in a burn-induced model of edema.

METHODS

A partial thickness and full thickness burns were induced in an in vivo rat model to elicit unique TWC perturbations corresponding to burn severity. Concomitant THz surface maps and MRI images of both burn models were acquired with a previously reported THz imaging system and T₂-weighted MRI, respectively, over 270 min. Reflectivity was analyzed for the burn contact area in THz images, while proton density (i.e., mobile TWC) was analyzed for the same region at incrementally increasing tissue depths in companion, transverse MRI images. A normalized cross correlation of THz and depth-dependent MRI measurements was performed as a function of time in histologically verified burn wounds.

RESULTS

For both burn types, strong positive correlations were evident between THz reflectivity and MRI data analyzed at greater tissue depths (>258 μm). MRI and THz results also revealed biphasic trends consistent with burn edema pathogenesis.

CONCLUSION

This paper offers the first in vivo correlative assessment of mobile TWC-based contrast and the sensing depth of THz imaging.

SIGNIFICANCE

The ability to implement THz imaging immediately following injury, combined with TWC sensing capabilities that compare to MRI, further support THz sensing as an emerging tool to track fluid in tissue.

Sympathetic modulation of electrical activation in normal and infarcted myocardium: implications for arrhythmogenesis

The influence of cardiac sympathetic innervation on electrical activation in normal and chronically infarcted ventricular myocardium is not understood. Yorkshire pigs with normal hearts (NL, n = 12) or anterior myocardial infarction (MI, n = 9) underwent high-resolution mapping of the anteroapical left ventricle at baseline and during left and right stellate ganglion stimulation (LSGS and RSGS, respectively). Conduction velocity (CV), activation times (ATs), and directionality of propagation were measured. Myocardial fiber orientation was determined using diffusion tensor imaging and histology. Longitudinal CV (CV_L) was increased by RSGS (0.98 ± 0.11 vs. 1.2 ± 0.14m/s, P < 0.001) but not transverse CV (CV_T). This increase was abrogated by β-adrenergic receptor and gap junction (GJ) blockade. Neither CV_L nor CV_T was increased by LSGS. In the peri-infarct region, both RSGS and LSGS shortened ARIs in sinus rhythm (423 ± 37 vs. 322 ± 30 ms, P < 0.001, and 423 ± 36 vs. 398 ± 36 ms, P = 0.035, respectively) and altered activation patterns in all animals. CV, as estimated by mean ATs, increased in a directionally dependent manner by RSGS (14.6 ± 1.2 vs. 17.3 ± 1.6 ms, P = 0.015), associated with GJ lateralization. RSGS and LSGS inhomogeneously modulated AT and induced relative or absolute functional activation delay in parts of the mapped regions in 75 and 67%, respectively, in MI animals, and in 0 and 15%, respectively, in control animals (P < 0.001 for both). In conclusion, sympathoexcitation increases CV in normal myocardium and modulates activation propagation in peri-infarcted ventricular myocardium. These data demonstrate functional control of arrhythmogenic peri-infarct substrates by sympathetic nerves and in part explain the temporal nature of arrhythmogenesis.NEW & NOTEWORTHY This study demonstrates regional control of conduction velocity in normal hearts by sympathetic nerves. In infarcted hearts, however, not only is modulation of propagation heterogeneous, some regions showed paradoxical conduction slowing. Sympathoexcitation altered propagation in all infarcted hearts studied, and we describe the temporal arrhythmogenic potential of these findings.Listen to this article's corresponding podcast at http://ajpheart.podbean.com/e/sympathetic-nerves-and-cardiac-propagation/.

John MA, Grundfest WS, Taylor ZD. Non-invasive terahertz imaging of tissue water content for flap viability assessment

Accurate and early prediction of tissue viability is the most significant determinant of tissue flap survival in reconstructive surgery. Perturbation in tissue water content (TWC) is a generic component of the tissue response to such surgeries, and, therefore, may be an important diagnostic target for assessing the extent of flap viability in vivo. We have previously shown that reflective terahertz (THz) imaging, a non-ionizing technique, can generate spatially resolved maps of TWC in superficial soft tissues, such as cornea and wounds, on the order of minutes. Herein, we report the first in vivo pilot study to investigate the utility of reflective THz TWC imaging for early assessment of skin flap viability. We obtained longitudinal visible and reflective THz imagery comparing 3 bipedicled flaps (i.e. survival model) and 3 fully excised flaps (i.e. failure model) in the dorsal skin of rats over a postoperative period of 7 days. While visual differences between both models manifested 48 hr after surgery, statistically significant (p < 0.05, independent t-test) local differences in TWC contrast were evident in THz flap image sets as early as 24 hr. Excised flaps, histologically confirmed as necrotic, demonstrated a significant, yet localized, reduction in TWC in the flap region compared to non-traumatized skin. In contrast, bipedicled flaps, histologically verified as viable, displayed mostly uniform, unperturbed TWC across the flap tissue. These results indicate the practical potential of THz TWC sensing to accurately predict flap failure 24 hours earlier than clinical examination.

Phase-contrast MRI with hybrid one and two-sided flow-encoding and velocity spectrum separation

PURPOSE

To develop and evaluate a phase-contrast MRI (PC-MRI) technique with hybrid one and two-sided flow-encoding and velocity spectrum separation (HOTSPA) for accelerated blood flow and velocity measurement.

METHODS

In the HOTSPA technique, the two-sided flow encoding (FE) is used for two FE directions and one-sided is used for the remaining FE direction. Such a temporal modulation of the FE strategy allows for separations of the Fourier velocity spectrum into components for the flow-compensated and the three-directional velocity waveforms, accelerating PC-MRI by encoding three-directional velocities using only two repetition times (TRs) instead of four TRs as in standard PC-MRI. The HOTSPA was evaluated and compared with standard PC-MRI in the common carotid arteries of six healthy volunteers.

RESULTS

Total volumetric flow and peak velocity measurements based on HOTSPA and the conventional PC-MRI were in good agreement with a bias of -0.005 mL (-0.1% relative bias error) for total volumetric flow and 1.21 cm/s (1.1% relative bias error) for peak velocity, although the total acquisition time was 50% of the conventional PC-MRI.

CONCLUSION

The proposed HOTSPA technique achieved nearly two-fold acceleration of PC-MRI while maintaining accuracy for total volumetric flow and peak velocity quantification by separating the paired acquisitions in the Fourier velocity spectrum domain. Magn Reson Med 78:182-192, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Scar voltage threshold determination using ex vivo magnetic resonance imaging integration in a porcine infarct model: Influence of interelectrode distances and three-dimensional spatial effects of scar

BACKGROUND

Studies analyzing optimal voltage thresholds for scar detection with electroanatomic mapping frequently lack a gold standard for comparison.

OBJECTIVE

The purpose of this study was to use a porcine infarct model with ex vivo magnetic resonance imaging (MRI) integration to characterize the relationship between interelectrode spacing and bipolar voltage thresholds and examine the influence of 3-dimensional scar on unipolar voltages.

METHODS

Thirty-two combined endocardial-epicardial electroanatomic maps were created in 8 postinfarct porcine subjects (bipolar 2-mm, 5-mm, and 8-mm interelectrode spacing and unipolar) for comparison with ex vivo MRI. Two thresholds were compared: (1) 95% normal distribution and (2) best fit to MRI. Direct electrogram analysis was performed in regions across from MRI-defined scar and adjacent to scar border zone.

RESULTS

A linear increase in optimal thresholds was observed with wider bipole spacing. The 95% thresholds for scar were lower than MRI-matched thresholds with moderate sensitivity for nontransmural scar (54% endo, 63% epi). Unipolar endocardial scar area exceeded MRI-defined scar, resulting in mismatched false scar in 5 of 8 (63%). Endocardial and epicardial unipolar voltages were lower than normal in regions adjacent and across from scar.

CONCLUSION

Variations in interelectrode spacing necessitate tailored bipolar voltage thresholds to optimize scar detection. Statistical 95% thresholds appear to be conservative and not fully sensitive for the detection of scar defined by high-resolution ex vivo MRI. In the presence of endocardial scar, unipolar mapping to quantitatively characterize epicardial scar may be overly sensitive due to 3-dimensional spatial averaging.

Convex optimized diffusion encoding (CODE) gradient waveforms for minimum echo time and bulk motion-compensated diffusion-weighted MRI

PURPOSE

To evaluate convex optimized diffusion encoding (CODE) gradient waveforms for minimum echo time and bulk motion-compensated diffusion-weighted imaging (DWI).

METHODS

Diffusion-encoding gradient waveforms were designed for a range of b-values and spatial resolutions with and without motion compensation using the CODE framework. CODE, first moment (M₁ ) nulled CODE-M₁ , and first and second moment (M₂ ) nulled CODE-M₁ M₂ were used to acquire neuro, liver, and cardiac ADC maps in healthy subjects (n=10) that were compared respectively to monopolar (MONO), BIPOLAR (M₁  = 0), and motion-compensated (MOCO, M₁  + M₂  = 0) diffusion encoding.

RESULTS

CODE significantly improved the SNR of neuro ADC maps compared with MONO (19.5 ± 2.5 versus 14.5 ± 1.9). CODE-M₁ liver ADCs were significantly lower (1.3 ± 0.1 versus 1.8 ± 0.3 × 10^-3 mm² /s, ie, less motion corrupted) and more spatially uniform (6% versus 55% ROI difference) than MONO and had higher SNR than BIPOLAR (SNR = 14.9 ± 5.3 versus 8.0 ± 3.1). CODE-M₁ M₂ cardiac ADCs were significantly lower than MONO (1.9 ± 0.6 versus 3.8 ± 0.3 x10^-3 mm² /s) throughout the cardiac cycle and had higher SNR than MOCO at systole (9.1 ± 3.9 versus 7.0 ± 2.6) while reporting similar ADCs (1.5 ± 0.2 versus 1.4 ± 0.6 × 10^-3 mm² /s).

CONCLUSIONS

CODE significantly improved SNR for ADC mapping in the brain, liver and heart, and significantly improved DWI bulk motion robustness in the liver and heart. Magn Reson Med 77:717-729, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Effect of free-breathing on left ventricular rotational mechanics in healthy subjects and patients with duchenne muscular dystrophy

PURPOSE

Cardiovascular magnetic resonance imaging exams can be performed during free-breathing. This may be especially important for boys with Duchenne muscular dystrophy (DMD) given their frequently limited breath-hold abilities. The impact of the respiratory compensation method on quantitative measurements of left ventricular (LV) rotational mechanics is incompletely understood. The purpose of this study was to evaluate differences in LV rotational mechanics acquired during breath-holding (BH), free-breathing with averaging (AVG), and free-breathing with respiratory bellows gating (BEL).

METHODS

LV short-axis tagged images from healthy subjects (N = 16) and DMD patients (N = 5) were acquired with BH, AVG, and BEL. LV twist and circumferential-longitudinal shear (CL-shear) angle were measured using the Fourier Analysis of STimulated echoes (FAST) method.

RESULTS

Peak LV twist estimates using BEL were significantly lower compared with BH in both healthy subjects (10.2 ± 3.6 ° versus 12.9 ± 2.3 °, P = 0.003) and patients with DMD (8.6 ± 3.6 ° versus 10.5 ± 3.6 °, P = 0.004). AVG results were in between BEL and BH. No significant differences in CL-shear were detected between BEL and BH.

CONCLUSION

Breath-holding directly affects estimates of peak LV twist, but not CL-shear. Using a free-breathing strategy for the evaluation of cardiac function is important for intrasubject longitudinal studies, intersubject comparisons, and multicenter trials for patients with DMD. Magn Reson Med 77:864-869, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Free-breathing variable flip angle balanced SSFP cardiac cine imaging with reduced SAR at 3T

PURPOSE

To develop a free-breathing variable flip angle (VFA) balanced steady-state free precession (bSSFP) cardiac cine imaging technique with reduced specific absorption rate (SAR) at 3 Tesla.

METHODS

Free-breathing VFA (FB-VFA) images in the short-axis and four-chamber views were acquired using an optimal VFA scheme, then compared with conventional breath-hold constant flip angle (BH-CFA) acquisitions. Two cardiac MRI experts used a 5-point scale to score images from healthy subjects (N = 10). The left ventricular ejection fraction, end diastolic volume (LVEDV), end systolic volume, stroke volume (LVSV), and end diastolic myocardial mass (LVEDM) were determined by manual contour analysis for BH-CFA and FB-VFA. A pilot evaluation of FB-VFA was performed in one patient with Duchenne muscular dystrophy.

RESULTS

FB-VFA SAR was 25% lower than BH-CFA with similar blood-myocardium contrast. The qualitative FB-VFA score was lower than the BH-CFA for the short-axis (3.1 ± 0.5 versus 4.3 ± 0.8; P < 0.05) and the four-chamber view (3.4 ± 0.4 versus 4.6 ± 0.6; P < 0.05). The LVEDV and the LVSV were 5% and 12% larger (P < 0.05) for FB-VFA compared with BH-CFA. There was no difference in LVEDM.

CONCLUSION

FB-VFA bSSFP cardiac cine imaging decreased the SAR at 3T with image quality sufficient to perform cardiac functional analysis. Magn Reson Med 76:1210-1216, 2016. © 2015 Wiley Periodicals, Inc.

Variable flip angle balanced steady-state free precession for lower SAR or higher contrast cardiac cine imaging

PURPOSE

Cardiac cine balanced steady-state free precession (bSSFP) imaging uses a high flip angle (FA) to obtain high blood-myocardium signal-to-noise and contrast-to-noise ratios (CNR). Use of high FAs, however, results in substantially increased SAR. Our objective was to develop a variable FA bSSFP cardiac cine imaging technique with: (1) low SAR and blood-myocardium CNR similar to conventional constant FA bSSFP (CFA-bSSFP) or (2) increased blood-myocardium CNR compared to CFA-bSSFP with similar SAR.

METHODS

Variable FA bSSFP cardiac cine imaging was achieved using an asynchronous k-space acquisition, which is asynchronous to the cardiac cycle (aVFA-bSSFP). Bloch simulations and phantom experiments were performed to compare the signal, resolution, and frequency response of the variable FA bSSFP and CFA-bSSFP schemes. Ten volunteers were imaged with different aVFA-bSSFP and asynchronous CFA-bSSFP schemes and compared to conventional segmented CFA-bSSFP.

RESULTS

The SAR of aVFA-bSSFP is significantly decreased by 36% compared to asynchronous CFA-bSSFP (1.9 ± 0.2 vs. 3.0 ± 0.2 W/kg, P <  10(-10)) for similar blood-myocardium CNR (34 ± 6 vs. 35 ± 9, P = 0.5). Alternately, the CNR of the aVFA-bSSFP is improved by 28% compared to asynchronous CFA-bSSFP (49 ± 9 vs. 38 ± 8, P <  10(-4)) with similar SAR (3.2 ± 0.5 vs. 3.3 ± 0.5 W/kg, P = 0.6).

CONCLUSION

aVFA-bSSFP can be used for lower SAR or higher contrast cardiac cine imaging compared to the conventional segmented CFA-bSSFP imaging.

Testing Foundations of Biological Scaling Theory Using Automated Measurements of Vascular Networks

Scientists have long sought to understand how vascular networks supply blood and oxygen to cells throughout the body. Recent work focuses on principles that constrain how vessel size changes through branching generations from the aorta to capillaries and uses scaling exponents to quantify these changes. Prominent scaling theories predict that combinations of these exponents explain how metabolic, growth, and other biological rates vary with body size. Nevertheless, direct measurements of individual vessel segments have been limited because existing techniques for measuring vasculature are invasive, time consuming, and technically difficult. We developed software that extracts the length, radius, and connectivity of in vivo vessels from contrast-enhanced 3D Magnetic Resonance Angiography. Using data from 20 human subjects, we calculated scaling exponents by four methods—two derived from local properties of branching junctions and two from whole-network properties. Although these methods are often used interchangeably in the literature, we do not find general agreement between these methods, particularly for vessel lengths. Measurements for length of vessels also diverge from theoretical values, but those for radius show stronger agreement. Our results demonstrate that vascular network models cannot ignore certain complexities of real vascular systems and indicate the need to discover new principles regarding vessel lengths.

Vascular networks distribute resources and constrain metabolic rate. Founded on a few key principles, biological scaling theories predict characteristic patterns for vascular networks as they branch from large to small vessels. These theories also predict seemingly unrelated phenomena, such as size limits on mammals. However, vascular networks are difficult to measure because there are billions of vessels that range in size from meters to micrometers. To test the foundations of biological scaling theories, we developed software that quickly measures thousands of in vivo vessels based on MRI. Data for vessel radii match predicted patterns but lengths do not. Our work suggests the need for new theoretical principles and should facilitate comparisons across organisms, spatial scales, and healthy and diseased tissue.

Abstract 328: Effects of Sympatho-excitation on Directional Anisotropy in Ventricular Myocardium

Introduction: Control of electrical propagation exerted by the sympathetic nervous system is has not been quantified in-depth.

Methods: High-resolution mapping (64-electrode plaque, 8x8, 1.96cm 2 ) of the anterior LV myocardium in porcine model (n=6) was performed before & during left stellate ganglion stimulation (LSGS). Activation times (AT), activation recovery intervals (ARIs), conduction velocities (CV), and CV anisotropy were obtained during pacing. Ex vivo diffusion tensor MRI and histology were performed to define myocardial fiber orientation in the mapped regions.

Results: LSGS shortened ARI (314.3±7.8ms vs. 287.8±6.6ms, p<0.001). At baseline, longitudinal CV (CV L ) was greater than transverse CV (CV T ) (1.2±0.2m/s vs 0.7±0.1m/s, p<0.001). LSGS did not increase CV L (1.16±0.16m/s vs 1.17±0.15m/s, p=0.2), or CV T (0.67±0.07m/s vs 0.7±0.04m/s, p=0.2). However, CV in the retrograde direction along fiber orientation was significantly increased by LSGS (0.7±0.05m/s vs 0.9±0.06m/s, p<0.01). This resulted in a significant reduction in directional anisotropy (CV antegrade /CV retrograde ) along fiber direction (1.6±0.3 vs 1.3±0.3, p<0.001), but not CV L /CV T (1.78±0.2 vs. 1.66±0.2 for BL and LSGS, respectively, p>0.2). Heterogeneity of activation (AT disp ) was greater in the transverse than longitudinal direction (20±2ms 2 /cm 2 vs. 13±1.5ms 2 /cm 2 , p= 0.022). LSGS decreased transverse AT disp to 18±1.8 ms 2 /cm 2 (p=0.015), but not longitudinal AT disp .

Conclusion: Cardiac sympatho-excitation modulates CV and AT disp in a fiber orientation-dependent manner. These uneven changes may contribute to arrhythmogenic mechanisms seen during sympatho-excitation.

Simulation Methods and Validation Criteria for Modeling Cardiac Ventricular Electrophysiology

We describe a sequence of methods to produce a partial differential equation model of the electrical activation of the ventricles. In our framework, we incorporate the anatomy and cardiac microstructure obtained from magnetic resonance imaging and diffusion tensor imaging of a New Zealand White rabbit, the Purkinje structure and the Purkinje-muscle junctions, and an electrophysiologically accurate model of the ventricular myocytes and tissue, which includes transmural and apex-to-base gradients of action potential characteristics. We solve the electrophysiology governing equations using the finite element method and compute both a 6-lead precordial electrocardiogram (ECG) and the activation wavefronts over time. We are particularly concerned with the validation of the various methods used in our model and, in this regard, propose a series of validation criteria that we consider essential. These include producing a physiologically accurate ECG, a correct ventricular activation sequence, and the inducibility of ventricular fibrillation. Among other components, we conclude that a Purkinje geometry with a high density of Purkinje muscle junctions covering the right and left ventricular endocardial surfaces as well as transmural and apex-to-base gradients in action potential characteristics are necessary to produce ECGs and time activation plots that agree with physiological observations.

Left ventricular twist and shear in patients with primary mitral regurgitation

PURPOSE

To evaluate the relationship between left ventricular (LV) twist, shear, and twist-per-volume and the severity of mitral regurgitation (MR). Primary MR is a valvular disorder that induces LV dysfunction. There exist several measures of LV rotational mechanics, but it remains unclear which measure of LV dysfunction best accords with the severity of MR. We hypothesized that LV systolic twist-per-volume slope would decrease with increasing severity of MR because of both decreases in rotational mechanics and increases in stroke volumes.

MATERIALS AND METHODS

Normal subjects (n = 54), moderate MR patients (n = 29), and severe MR patients (n = 54) were studied. Magnetic resonance imaging (MRI) was performed on a 1.5T scanner and grid-tagged LV images were collected at the LV base and LV apex. Measures of LV rotational mechanics were derived from tagged images using Fourier Analysis of STimulated echoes (FAST).

RESULTS

Peak systolic twist-per-volume slope was significantly different for all pairwise comparisons (P < 0.0001) and compared to normal subjects (-0.14 ± 0.05°/mL) was decreased in moderate MR (-0.12 ± 0.04°/mL) and further decreased in severe MR (-0.07 ± 0.03°/mL).

CONCLUSION

Peak systolic twist-per-volume slope significantly decreased with increasing severity of MR and is therefore a suitable quantitative imaging biomarker for LV dysfunction in patients with MR.

Probing dynamic myocardial microstructure with cardiac magnetic resonance diffusion tensor imaging

This article is an invited editorial comment on the paper entitled "In vivo cardiovascular magnetic resonance diffusion tensor imaging shows evidence of abnormal myocardial laminar orientations and mobility in hypertrophic cardiomyopathy" by Ferreira et al., and published as Journal of Cardiovascular Magnetic Resonance 2014; 16:87.

Fast 3D T2 -weighted imaging using variable flip angle transition into driven equilibrium (3D T2 -TIDE) balanced SSFP for prostate imaging at 3T

PURPOSE

Three-dimensional (3D) T2 -weighted fast spin echo (FSE) imaging of the prostate currently requires long acquisition times. Our objective was to develop a fast 3D T2 -weighted sequence for prostate imaging at 3T using a variable flip angle transition into driven equilibrium (T2 -TIDE) scheme.

METHODS

3D T2 -TIDE uses interleaved spiral-out phase encode ordering to efficiently sample the ky -kz phase encodes and also uses the transient balanced steady-state free precession signal to acquire the center of k-space for T2 -weighted imaging. Bloch simulations and images from 10 healthy subjects were acquired to evaluate the performance of 3D T2 -TIDE compared to 3D FSE.

RESULTS

3D T2 -TIDE images were acquired in 2:54 minutes compared to 7:02 minutes for 3D FSE with identical imaging parameters. The signal-to-noise ratio (SNR) efficiency was significantly higher for 3D T2 -TIDE compared to 3D FSE in nearly all tissues, including periprostatic fat (45 ± 12 vs. 31 ± 7, P < 0.01), gluteal fat (48 ± 8 vs. 41 ± 10, P = 0.12), right peripheral zone (20 ± 4 vs. 16 ± 8, P = 0.12), left peripheral zone (17 ± 2 vs. 12 ± 3, P < 0.01), and anterior fibromuscular stroma (12 ± 4 vs. 4 ± 2, P < 0.01).

CONCLUSION

3D T2 -TIDE images of the prostate can be acquired quickly with SNR efficiency that exceeds that of 3D FSE.

Complementary radial tagging for improved myocardial tagging contrast

PURPOSE

To develop and evaluate complementary radial tagging (CRT) for improved myocardial tagging contrast.

METHODS

We sought to develop and evaluate CRT, which aims to preserve the radial tag contrast throughout the cardiac cycle. Similar to complementary spatial modulation of magnetization, CRT acquires two sets of images with a phase shift in the tag pattern. The combination of a ramped imaging flip angle and image subtraction enhances tag contrast throughout the cardiac cycle. The proposed CRT technique uses a small table shift away from the isocenter to improve the uniformity of the radial tag pattern. We provide a mathematical solution for the optimal table shift and validate the solution in using a retrospective analysis of images from 500 patients in the Cardiac Atlas Project database.

RESULTS

CRT simulations, phantom experiments, and in vivo images all demonstrate the improved tag contrast of CRT compared to RT. The retrospective evaluation demonstrated that acceptable CRT images could be acquired in over 98% of the clinical exams.

CONCLUSION

The CRT technique improves radial tag contrast throughout the cardiac cycle and should produce high quality tag patterns in nearly all patients.

Phase contrast MRI with flow compensation view sharing

PURPOSE

To develop and evaluate a technique for accelerating phase contrast MRI (PC-MRI) acquisitions without significant compromise in flow quantification accuracy.

METHODS

PC-MRI is commonly acquired using interleaved flow-compensated (FC) and flow-encoded (FE) echoes. We hypothesized that FC data, which represent background phase, do not change significantly over time. Therefore, we proposed to undersample the FC data and use an FC view sharing (FCVS) approach to synthesize a composite FC frame for each corresponding FE frame. FCVS was evaluated in a flow phantom and healthy volunteers and compared with a standard FC/FE PC-MRI.

RESULTS

The FCVS sequence resulted in an error of 0.0% for forward flow and 2.0% for reverse flow volume when compared with FC/FE PC-MRI in a flow phantom. Measurements in the common carotid arteries showed that the FCVS method had -1.16 cm/s bias for maximum peak velocity and -0.019 mL bias in total flow, when compared with FC/FE with the same temporal resolution, but double the total acquisition time. These results represent ≤1.3% bias error in velocity and volumetric flow quantification.

CONCLUSION

FCVS can accelerate PC-MRI acquisitions while maintaining flow and velocity measurement accuracy when there is limited temporal variation in the FC data.

Device artifact reduction for magnetic resonance imaging of patients with implantable cardioverter-defibrillators and ventricular tachycardia: late gadolinium enhancement correlation with electroanatomic mapping

BACKGROUND

Late gadolinium enhancement (LGE) magnetic resonance imaging (MRI) of ventricular scar has been shown to be accurate for detection and characterization of arrhythmia substrates. However, the majority of patients referred for ventricular tachycardia (VT) ablation have an implantable cardioverter-defibrillator (ICD), which obscures image integrity and the clinical utility of MRI.

OBJECTIVE

The purpose of this study was to develop and validate a wideband LGE MRI technique for device artifact removal.

METHODS

A novel wideband LGE MRI technique was developed to allow for improved scar evaluation on patients with ICDs. The wideband technique and the standard LGE MRI were tested on 18 patients with ICDs. VT ablation was performed in 13 of 18 patients with either endocardial and/or epicardial approach and the correlation between the scar identified on MRI and electroanatomic mapping (EAM) was analyzed.

RESULTS

Hyperintensity artifact was present in 16 of 18 of patients using standard MRI, which was eliminated using the wideband LGE and allowed for MRI interpretation in 15 of 16 patients. All patients had ICD lead characteristics confirmed as unchanged post-MRI and had no adverse events. LGE scar was seen in 11 of 18 patients. Among the 15 patients in whom wideband LGE allowed visualization of myocardium, 10 had LGE scar and 5 had normal myocardium in the regions with image artifacts when using the standard LGE. The left ventricular scar size measurements using wideband MRI and EAM were correlated with R(2) = 0.83 and P = .00003.

CONCLUSION

Wideband LGE MRI improves the ability to visualize myocardium for clinical interpretation, which correlated well with EAM findings during VT ablation.

Noncontrast enhanced four-dimensional dynamic MRA with golden angle radial acquisition and K-space weighted image contrast (KWIC) reconstruction

PURPOSE

To explore the feasibility of 2D and 3D golden-angle radial acquisition strategies in conjunction with k-space weighted image contrast (KWIC) temporal filtering to achieve noncontrast enhanced dynamic MRA (dMRA) with high spatial resolution, low streaking artifacts and high temporal fidelity.

METHODS

Simulations and in vivo examinations in eight normal volunteers and an arteriovenous malformation patient were carried out. Both 2D and 3D golden angle radial sequences, preceded by spin tagging, were used for dMRA of the brain. The radial dMRA data were temporally filtered using the KWIC strategy and compared with matched standard Cartesian techniques.

RESULTS

The 2D and 3D dynamic MRA image series acquired with the proposed radial techniques demonstrated excellent image quality without discernible temporal blurring compared with standard Cartesian based approaches. The image quality of radial dMRA was equivalent to or higher than that of Cartesian dMRA by visual inspection. A reduction factor of up to 10 and 3 in scan time was achieved for 2D and 3D radial dMRA compared with the Cartesian-based counterparts.

CONCLUSION

The proposed 2D and 3D radial dMRA techniques demonstrated image quality comparable or even superior to those obtained with standard Cartesian methods, but within a fraction of the scan time.

Intra- and interscan reproducibility using Fourier Analysis of STimulated Echoes (FAST) for the rapid and robust quantification of left ventricular twist

PURPOSE

To assess the intra- and interscan reproducibility of LV twist using FAST. Assessing the reproducibility of the measurement of new MRI biomarkers is an important part of validation. Fourier Analysis of STimulated Echoes (FAST) is a new MRI tissue tagging method that has recently been shown to compare favorably with conventional estimates of left ventricular (LV) twist from cardiac tagged images, but with significantly reduced user interaction time.

MATERIALS AND METHODS

Healthy volunteers (N = 10) were scanned twice using FAST over 1 week. On day 1, two measurements of LV twist were collected for intrascan comparisons. Measurements for LV twist were again collected on day 8 for interscan assessment. LV short-axis tagged images were acquired on a 3 Tesla (T) scanner to ensure detectability of tags during early and mid-diastole. Peak LV twist is reported as mean ± SD. Reproducibility was assessed using the concordance correlation coefficient (CCC) and the repeatability coefficient (RC) (95% confidence interval [CI] range).

RESULTS

Mean peak twist measurements were 13.4 ± 4.3° (day 1, scan 1), 13.6 ± 3.7° (day 1, scan 2), and 13.0 ± 2.7° (day 8). Bland-Altman analysis resulted in intra- and interscan bias and 95% CI of -0.6° [-1.0°, 1.6°] and 1.4° (-1.0°, 3.0°), respectively. The Bland-Altman RC for peak LV twist was 2.6° and 4.0° for intra- and interscan, respectively. The CCC was 0.9 and 0.6 for peak LV twist for intra- and interscan, respectively.

CONCLUSION

FAST is a semi-automated method that provides a quick and quantitative assessment of LV systolic and diastolic twist that demonstrates high intrascan and moderate interscan reproducibility in preliminary studies.

Off-resonance insensitive complementary SPAtial Modulation of Magnetization (ORI-CSPAMM) for quantification of left ventricular twist

PURPOSE

To evaluate Off Resonance Insensitive Complementary SPAtial Modulation of Magnetization (ORI-CSPAMM) and Fourier Analysis of STimulated echoes (FAST) for the quantification of left ventricular (LV) systolic and diastolic function and compare it with the previously validated FAST+SPAMM technique.

MATERIALS AND METHODS

LV short-axis tagged images were acquired with ORI-CSPAMM and SPAMM in healthy volunteers (n = 13). The FAST method was used to automatically estimate LV systolic and diastolic twist parameters from rotation of the stimulated echo and stimulated anti-echo about the middle of k-space subsequent to ∼3 min of user interaction.

RESULTS

There was no significant difference between measures obtained for FAST+ORI-CSPAMM and FAST+SPAMM for mean peak twist (12.9 ± 3.4° versus 11.9 ± 4.0°; P = 0.4), torsion (3.3 ± 0.9°/cm versus 2.9 ± 1.0°/cm, P = 0.3), circumferential-longitudinal shear angle (9.1 ± 3.0° versus 8.2 ± 3.4°, P = 0.3), twisting rate (79.6 ± 20.2°/s versus 68.2 ± 23.4°/s, P = 0.1), untwisting rate (-117.5 ± 31.4°/s versus -106.6 ± 32.4°/s, P = 0.3), normalized untwisting rate (-9.3 ± 2.0/s versus -9.9 ± 4.4/s, P = 0.7), and time of peak twist (281 ± 18 ms versus 293 ± 25 ms, P = 0.04). FAST+ORI-CSPAMM also provided measures of duration of untwisting (148 ± 21 ms) and the ratio of rapid untwisting to peak twist (0.9 ± 0.3). Bland-Altman analysis of FAST+ORI-CSPAMM and FAST+SPAMM twist data demonstrates excellent agreement with a bias of -0.1° and 95% confidence intervals of (-1.0°, 3.2°).

CONCLUSION

FAST+ORI-CSPAMM is a semi-automated method that provides a quick and quantitative assessment of LV systolic and diastolic twist and torsion. ORI-CSPAMM corrects off-resonance accrued during tagging preparation and readout and visibly removes chemical shift from the tagging pattern, which confers greater robustness to the derived quantitative measures.