Cardiac MR elastography: comparison with left ventricular pressure measurement

Multiple motion encoding in phase-contrast MRI: A general theory and application to elastography imaging

While MRI allows to encode the motion of tissue in the magnetization's phase, it remains yet a challenge to obtain high fidelity motion images due to wraps in the phase for high encoding efficiencies. Therefore, we propose an optimal multiple motion encoding method (OMME) and exemplify it in Magnetic Resonance Elastography (MRE) data. OMME is formulated as a non-convex least-squares problem for the motion using an arbitrary number of phase-contrast measurements with different motion encoding gradients (MEGs). The mathematical properties of OMME are proved in terms of standard deviation and dynamic range of the motion's estimate for arbitrary MEGs combination which are confirmed using synthetically generated data. OMME's performance is assessed on MRE data from in vivo human brain experiments and compared to dual encoding strategies. The unwrapped images are further used to reconstruct stiffness maps and compared to the ones obtained using conventional unwrapping methods. OMME allowed to successfully combine several MRE phase images with different MEGs, outperforming dual encoding strategies in either motion-to-noise ratio (MNR) or number of successfully reconstructed voxels with good noise stability. This lead to stiffness maps with greater resolution of details than obtained with conventional unwrapping methods. The proposed OMME method allows for a flexible and noise robust increase in the dynamic range and thus provides wrap-free phase images with high MNR. In MRE, the method may be especially suitable when high resolution images with high MNR are needed.

Cardiovascular magnetic resonance elastography: A review

Cardiovascular diseases are the leading cause of death worldwide. These cardiovascular diseases are associated with mechanical changes in the myocardium and aorta. It is known that stiffness is altered in many diseases, including the spectrum of ischemia, diastolic dysfunction, hypertension and hypertrophic cardiomyopathy. In addition, the stiffness of the aortic wall is altered in multiple diseases, including hypertension, coronary artery disease and aortic aneurysm formation. For example, in diastolic dysfunction in which the ejection fraction is preserved, stiffness can potentially be an important biomarker. Similarly, in aortic aneurysms, stiffness can provide valuable information with regard to rupture potential. A number of studies have addressed invasive and non-invasive approaches to test and measure the mechanical properties of the myocardium and aorta. One of the non-invasive approaches is magnetic resonance elastography (MRE). MRE is a phase-contrast magnetic resonance imaging technique that measures tissue stiffness non-invasively. This review article highlights the technical details and application of MRE in the quantification of myocardial and aortic stiffness in different disease states.

Relative identifiability of anisotropic properties from magnetic resonance elastography

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

Stiffness reconstruction methods for MR elastography

Assessment of tissue stiffness is desirable for clinicians and researchers, as it is well established that pathophysiological mechanisms often alter the structural properties of tissue. Magnetic resonance elastography (MRE) provides an avenue for measuring tissue stiffness and has a long history of clinical application, including staging liver fibrosis and stratifying breast cancer malignancy. A vital component of MRE consists of the reconstruction algorithms used to derive stiffness from wave-motion images by solving inverse problems. A large range of reconstruction methods have been presented in the literature, with differing computational expense, required user input, underlying physical assumptions, and techniques for numerical evaluation. These differences, in turn, have led to varying accuracy, robustness, and ease of use. While most reconstruction techniques have been validated against in silico or in vitro phantoms, performance with real data is often more challenging, stressing the robustness and assumptions of these algorithms. This article reviews many current MRE reconstruction methods and discusses the aforementioned differences. The material assumptions underlying the methods are developed and various approaches for noise reduction, regularization, and numerical discretization are discussed. Reconstruction methods are categorized by inversion type, underlying assumptions, and their use in human and animal studies. Future directions, such as alternative material assumptions, are also discussed.

Estimation of transversely isotropic material properties from magnetic resonance elastography using the optimised virtual fields method

Magnetic resonance elastography (MRE) has been used to estimate isotropic myocardial stiffness. However, anisotropic stiffness estimates may give insight into structural changes that occur in the myocardium as a result of pathologies such as diastolic heart failure. The virtual fields method (VFM) has been proposed for estimating material stiffness from image data. This study applied the optimised VFM to identify transversely isotropic material properties from both simulated harmonic displacements in a left ventricular (LV) model with a fibre field measured from histology as well as isotropic phantom MRE data. Two material model formulations were implemented, estimating either 3 or 5 material properties. The 3-parameter formulation writes the transversely isotropic constitutive relation in a way that dissociates the bulk modulus from other parameters. Accurate identification of transversely isotropic material properties in the LV model was shown to be dependent on the loading condition applied, amount of Gaussian noise in the signal, and frequency of excitation. Parameter sensitivity values showed that shear moduli are less sensitive to noise than the other parameters. This preliminary investigation showed the feasibility and limitations of using the VFM to identify transversely isotropic material properties from MRE images of a phantom as well as simulated harmonic displacements in an LV geometry.

Time-Harmonic Ultrasound elastography of the Descending Abdominal Aorta: Initial Results

Stiffening of central large vessels is considered a key pathophysiologic factor within the cardiovascular system. Current diagnostic parameters such as pulse wave velocity (PWV) indirectly measure aortic stiffness, a hallmark of coronary diseases. The aim of the present study was to perform elastography of the proximal abdominal aorta based on externally induced time-harmonic shear waves. Experiments were performed in 30 healthy volunteers (25 young, 5 old, >50 y) and 5 patients with longstanding hypertension (PWV >10 m/s). B-Mode-guided sonographic time-harmonic elastography was used for measurement of externally induced shear waves at 30-Hz vibration frequency. Thirty-hertz shear wave amplitudes (SWAs) within the abdominal aorta were measured and displayed in real time and processed offline for differences in SWA between systole and diastole (ΔSWA). Data were analyzed using the Kruskal-Wallis test and receiver operating characteristic curve analysis. The change in SWA over the cardiac cycle was reduced significantly in all patients as assessed with ΔSWA (volunteers: mean = 10 ± 5 μm, patients: mean = 4 ± 1 μm; p < 0.001). The best separation of healthy volunteers from patients was obtained with a ΔSWA threshold of 4.7 μm, resulting in a sensitivity of 0.9 and a specificity of 1.0, with an overall area under the curve of 0.96. Time harmonic elastography of the abdominal aorta is feasible and shows promise for the exploitation of time-varying shear wave amplitudes as a diagnostic marker for aortic wall stiffening. Patients with elevated PWVs suggesting increased aortic wall stiffness were best identified by ΔSWA-a parameter that could be related to the ability of the vessel walls to distend on passages of the pulse wave.

Cardiac MR elastography for quantitative assessment of elevated myocardial stiffness in cardiac amyloidosis

Purpose

To evaluate if cardiac magnetic resonance elastography (MRE) can measure increased stiffness in patients with cardiac amyloidosis. Myocardial tissue stiffness plays an important role in cardiac function. A noninvasive quantitative imaging technique capable of measuring myocardial stiffness could aid in disease diagnosis, therapy monitoring, and disease prognostic strategies. We recently developed a high‐frequency cardiac MRE technique capable of making noninvasive stiffness measurements.

Materials and Methods

In all, 16 volunteers and 22 patients with cardiac amyloidosis were enrolled in this study after Institutional Review Board approval and obtaining formal written consent. All subjects were imaged head‐first in the supine position in a 1.5T closed‐bore MR imager. 3D MRE was performed using 5 mm isotropic resolution oblique short‐axis slices and a vibration frequency of 140 Hz to obtain global quantitative in vivo left ventricular stiffness measurements. The median stiffness was compared between the two cohorts. An octahedral shear strain signal‐to‐noise ratio (OSS‐SNR) threshold of 1.17 was used to exclude exams with insufficient motion amplitude.

Results

Five volunteers and six patients had to be excluded from the study because they fell below the 1.17 OSS‐SNR threshold. The myocardial stiffness of cardiac amyloid patients (median: 11.4 kPa, min: 9.2, max: 15.7) was significantly higher (P = 0.0008) than normal controls (median: 8.2 kPa, min: 7.2, max: 11.8).

Conclusion

This study demonstrates the feasibility of 3D high‐frequency cardiac MRE as a contrast‐agent‐free diagnostic imaging technique for cardiac amyloidosis.

Level of Evidence: 2

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2017;46:1361–1367.

Investigating Shear Wave Physics in a Generic Pediatric Left Ventricular Model via In Vitro Experiments and Finite Element Simulations

Shear wave elastography (SWE) is a potentially valuable tool to noninvasively assess ventricular function in children with cardiac disorders, which could help in the early detection of abnormalities in muscle characteristics. Initial experiments demonstrated the potential of this technique in measuring ventricular stiffness; however, its performance remains to be validated as complicated shear wave (SW) propagation characteristics are expected to arise due to the complex non-homogenous structure of the myocardium. In this work, we investigated the (i) accuracy of different shear modulus estimation techniques (time-of-flight (TOF) method and phase velocity analysis) across myocardial thickness and (ii) effect of the ventricular geometry, surroundings, acoustic loading, and material viscoelasticity on SW physics. A generic pediatric (10-15-year old) left ventricular model was studied numerically and experimentally. For the SWE experiments, a polyvinylalcohol replicate of the cardiac geometry was fabricated and SW acquisitions were performed on different ventricular areas using varying probe orientations. Additionally, the phantom's stiffness was obtained via mechanical tests. The results of the SWE experiments revealed the following trends for stiffness estimation across the phantom's thickness: a slight stiffness overestimation for phase speed analysis and a clear stiffness underestimation for the TOF method for all acquisitions. The computational model provided valuable 3-D insights in the physical factors influencing SW patterns, especially the surroundings (water), interface force, and viscoelasticity. In conclusion, this paper presents a validation study of two commonly used shear modulus estimators for different ventricular locations and the essential role of SW modeling in understanding SW physics in the pediatric myocardium.

Interplay of Aging and Hypertension in Cardiac Remodeling: A Mathematical Geometric Model

Hypertensive disorder can cause cardiac deformities. Elastic characteristic parameters, like Young’s modulus of elasticity (E) derived from a traditional cylindrical model, increase significantly with aging. However, the geometric and component changes of aging hearts because of chronic hypertension remain unknown. To better describe the effects, we propose an elliptical elastic and mathematical model to evaluate myocardial stiffness. Ninety-six hypertensive patients (HTN^Pos) (men: 59.3%; age ≥ 65 years: 20.8%) were enrolled and compared with normotensive controls (HTN^Neg) (n = 47, 48.9%). HTN^Pos patients had a thicker interventricular septum in diastole (IVSd) (HTN^Pos: 0.96 ± 0.21 cm vs. HTN^Neg: 0.77 ± 0.15; p = 0.005) and higher intracardiac pressure (e/e′: 9.06 ± 4.85 cm vs. 7.76 ± 3.41; p = 0.01), especially the elderly (> 65 years) (IVSd: 1.03 ± 0.19 cm, e/e′: 11.39 ± 1.99; p = 0.006 and 0.01, respectively). Nevertheless, the internal dimension decreased more significantly in the HTN^Pos rather than in the HTN^Neg elderly (5.23 ± 0.46 vs. 4.74 ± 0.69 cm; p = 0.02). We found different directions of cardiac remodeling with normotensive and hypertensive loads. Different from the longitudinal and circumferential strain, E and Poisson’s ratio (υ) are values that directly present the rigidity of myocardium. E was significantly higher in the elderly (8011.92 ± 2431.85 vs. 6052.43 ± 3121.50; p = 0.02), whereas υ was significantly higher in all HTN^Pos patients (0.73 ± 0.12 vs. 0.61 ± 0.07; p < 0.001). Because E and υ reflected the material changes of myocardium in the HTN^Pos elderly, the proposed elliptical mathematical heart model better describes the geometric deformity induced by aging and hypertension.

In vivo quantification of myocardial stiffness in hypertensive porcine hearts using MR elastography

PURPOSE

To determine alteration in left ventricular (LV) myocardial stiffness (MS) with hypertension (HTN). Cardiac MR elastography (MRE) was used to estimate MS in HTN induced pigs and MRE-derived MS measurements were compared against LV pressure, thickness and circumferential strain.

MATERIALS AND METHODS

Renal-wrapping surgery was performed to induce HTN in eight pigs. LV catheterization (to measure pressure) and cardiac MRI (1.5 Tesla; gradient echo-MRE and tagging) was performed pre-surgery at baseline (Bx), and post-surgery at month 1 (M1) and month 2 (M2). Images were analyzed to estimate LV-MS, thickness, and circumferential strain across the cardiac cycle. The associations between end-diastolic (ED) and end-systolic (ES) MS and (i) mean LV pressure; (ii) ED and ES thickness, respectively; and (iii) circumferential strain were evaluated using Spearman's correlation method.

RESULTS

From Bx to M2, mean pressure, MRE-derived stiffness, and thickness increased while circumferential strain decreased significantly (slope test, P ≤ 0.05). Both ED and ES MS had significant positive correlation with (i) mean pressure (ED MS: ρ = 0.56; P = 0.005 and ES MS: ρ = 0.45; P = 0.03); (ii) ED thickness ( ρ = 0.73; P < 0.0001) and ES thickness ( ρ = 0.84; P < 0.0001), respectively; but demonstrated a negative trend with circumferential strain (ED MS: ρ = 0.31 and ES MS: ρ = 0.37).

CONCLUSION

This study demonstrated that, in a HTN porcine model, MRE-derived MS increased with increase in pressure and thickness.

LEVEL OF EVIDENCE

1 J. Magn. Reson. Imaging 2017;45:813-820.

General review of magnetic resonance elastography

Magnetic resonance elastography (MRE) is an innovative imaging technique for the non-invasive quantification of the biomechanical properties of soft tissues via the direct visualization of propagating shear waves in vivo using a modified phase-contrast magnetic resonance imaging (MRI) sequence. Fundamentally, MRE employs the same physical property that physicians utilize when performing manual palpation - that healthy and diseased tissues can be differentiated on the basis of widely differing mechanical stiffness. By performing “virtual palpation”, MRE is able to provide information that is beyond the capabilities of conventional morphologic imaging modalities. In an era of increasing adoption of multi-parametric imaging approaches for solving complex problems, MRE can be seamlessly incorporated into a standard MRI examination to provide a rapid, reliable and comprehensive imaging evaluation at a single patient appointment. Originally described by the Mayo Clinic in 1995, the technique represents the most accurate non-invasive method for the detection and staging of liver fibrosis and is currently performed in more than 100 centers worldwide. In this general review, the mechanical properties of soft tissues, principles of MRE, clinical applications of MRE in the liver and beyond, and limitations and future directions of this discipline -are discussed. Selected diagrams and images are provided for illustration.

In vivo, high-frequency three-dimensional cardiac MR elastography: Feasibility in normal volunteers

PURPOSE

Noninvasive stiffness imaging techniques (elastography) can image myocardial tissue biomechanics in vivo. For cardiac MR elastography (MRE) techniques, the optimal vibration frequency for in vivo experiments is unknown. Furthermore, the accuracy of cardiac MRE has never been evaluated in a geometrically accurate phantom. Therefore, the purpose of this study was to determine the necessary driving frequency to obtain accurate three-dimensional (3D) cardiac MRE stiffness estimates in a geometrically accurate diastolic cardiac phantom and to determine the optimal vibration frequency that can be introduced in healthy volunteers.

METHODS

The 3D cardiac MRE was performed on eight healthy volunteers using 80 Hz, 100 Hz, 140 Hz, 180 Hz, and 220 Hz vibration frequencies. These frequencies were tested in a geometrically accurate diastolic heart phantom and compared with dynamic mechanical analysis (DMA).

RESULTS

The 3D Cardiac MRE was shown to be feasible in volunteers at frequencies as high as 180 Hz. MRE and DMA agreed within 5% at frequencies greater than 180 Hz in the cardiac phantom. However, octahedral shear strain signal to noise ratios and myocardial coverage was shown to be highest at a frequency of 140 Hz across all subjects.

CONCLUSION

This study motivates future evaluation of high-frequency 3D MRE in patient populations. Magn Reson Med 77:351-360, 2017. © 2016 Wiley Periodicals, Inc.

Cardiac MR elastography of the mouse: Initial results

PURPOSE

Many cardiovascular diseases are associated with abnormal function of myocardial contractility or dilatability, which is related to elasticity changes of the myocardium over the cardiac cycle. The mouse is a common animal model in studies of the progression of various cardiomyopathies. We introduce a novel noninvasive approach using microscopic scale MR elastography (MRE) to measure the myocardium stiffness change during the cardiac cycle on a mouse model.

METHODS

A harmonic mechanical wave of 400 Hz was introduced into the mouse body. An electrocardiograph-gated and respiratory-gated fractional encoding cine-MRE pulse sequence was applied to encode the resulting oscillatory motion on a short-axis slice of the heart. Five healthy mice (age range, 3-13.5 mo) were examined. The weighted summation effective stiffness of the left ventricle wall during the cardiac cycle was estimated.

RESULTS

The ratio of stiffness at end diastole and end systole was 0.5-0.67. Additionally, variation in shear wave amplitude in the left ventricle wall throughout the cardiac cycle was measured and found to correlate with estimates of stiffness variation.

CONCLUSION

This study demonstrates the feasibility of implementing cardiac MRE on a mouse model. Magn Reson Med 76:1879-1886, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Investigation of the effects of myocardial anisotropy for shear wave elastography using impulsive force and harmonic vibration

The myocardium is known to be an anisotropic medium where the muscle fiber orientation changes through the thickness of the wall. Shear wave elastography methods use propagating waves which are measured by ultrasound or magnetic resonance imaging (MRI) techniques to characterize the mechanical properties of various tissues. Ultrasound- or MR-based methods have been used and the excitation frequency ranges for these various methods cover a large range from 24-500 Hz. Some of the ultrasound-based methods have been shown to be able to estimate the fiber direction. We constructed a model with layers of elastic, transversely isotropic materials that were oriented at different angles to simulate the heart wall in systole and diastole. We investigated the effect of frequency on the wave propagation and the estimation of fiber direction and wave speeds in the different layers of the assembled models. We found that waves propagating at low frequencies such as 30 or 50 Hz showed low sensitivity to the fiber direction but also had substantial bias in estimating the wave speeds in the layers. Using waves with higher frequency content (>200 Hz) allowed for more accurate fiber direction and wave speed estimation. These results have particular relevance for MR- and ultrasound-based elastography applications in the heart.

Quantification of regional aortic stiffness using MR elastography: A phantom and ex-vivo porcine aorta study

MR Elastography (MRE) is a noninvasive technique for measuring tissue stiffness that has been used to assess the average stiffness of the abdominal aorta. The utility of aortic MRE would be improved if it could provide information about local variations in aortic stiffness. We hypothesize that regional variations in aortic stiffness can also be measured with MRE and the purpose of this work was to demonstrate that MRE can measure regional stiffness variations in a vascular phantom and in ex vivo porcine aortas. A vascular phantom was fabricated, containing two silicone tubes embedded in gel. A segment of one of the tubes was modified to increase its stiffness. MRE was performed on the phantom with a continuous flow of water through the tubes. The stiffness distribution along the modified tube was measured and compared to the reference tube. MRE was also performed in porcine aortas embedded in gel with segments treated with saline or formalin for 4 days. The stiffness difference between saline- and formalin-treated aortic segments was measured by MRE and mechanical tests. A positive correlation was found between the regional stiffnesses measured by MRE and mechanical tests. The results indicate that MRE can be used to evaluate the local stiffness distribution in silicone tubes and ex vivo porcine aortas. It may therefore be possible to apply MRE to measure regional stiffness variations of the aorta in vivo.

Measuring age-dependent myocardial stiffness across the cardiac cycle using MR elastography: A reproducibility study

PURPOSE

To assess reproducibility in measuring left ventricular (LV) myocardial stiffness in volunteers throughout the cardiac cycle using MR elastography (MRE) and to determine its correlation with age.

METHODS

Cardiac MRE (CMRE) was performed on 29 normal volunteers, with ages ranging from 21 to 73 years. For assessing reproducibility of CMRE-derived stiffness measurements, scans were repeated per volunteer. Wave images were acquired throughout the LV myocardium, and were analyzed to obtain mean stiffness during the cardiac cycle. CMRE-derived stiffness values were correlated to age.

RESULTS

Concordance correlation coefficient revealed good interscan agreement with rc of 0.77, with P-value < 0.0001. Significantly higher myocardial stiffness was observed during end-systole (ES) compared with end-diastole (ED) across all subjects. Additionally, increased deviation between ES and ED stiffness was observed with increased age.

CONCLUSION

CMRE-derived stiffness is reproducible, with myocardial stiffness changing cyclically across the cardiac cycle. Stiffness is significantly higher during ES compared with ED. With age, ES myocardial stiffness increases more than ED, giving rise to an increased deviation between the two.

Quantitative elasticity measurement of urinary bladder wall using laser-induced surface acoustic waves

The maintenance of urinary bladder elasticity is essential to its functions, including the storage and voiding phases of the micturition cycle. The bladder stiffness can be changed by various pathophysiological conditions. Quantitative measurement of bladder elasticity is an essential step toward understanding various urinary bladder disease processes and improving patient care. As a nondestructive, and noncontact method, laser-induced surface acoustic waves (SAWs) can accurately characterize the elastic properties of different layers of organs such as the urinary bladder. This initial investigation evaluates the feasibility of a noncontact, all-optical method of generating and measuring the elasticity of the urinary bladder. Quantitative elasticity measurements of ex vivo porcine urinary bladder were made using the laser-induced SAW technique. A pulsed laser was used to excite SAWs that propagated on the bladder wall surface. A dedicated phase-sensitive optical coherence tomography (PhS-OCT) system remotely recorded the SAWs, from which the elasticity properties of different layers of the bladder were estimated. During the experiments, series of measurements were performed under five precisely controlled bladder volumes using water to estimate changes in the elasticity in relation to various urinary bladder contents. The results, validated by optical coherence elastography, show that the laser-induced SAW technique combined with PhS-OCT can be a feasible method of quantitative estimation of biomechanical properties.

Detection and characterisation of biopsy tissue using quantitative optical coherence elastography (OCE) in men with suspected prostate cancer

We present first quantitative three-dimensional (3D) data sets recorded using optical coherence elastography (OCE) for the diagnosis and detection of prostate cancer (PCa). 120 transrectal ultrasound guided prostate biopsy specimens from 10 men suspected with prostate cancer were imaged using OCE. 3D quantitative mechanical assessment of biopsy specimens obtained in kilopascals (kPa) at an interval of 40 µm was compared with histopathology. Sensitivity, specificity, and positive and negative predictive values were calculated for OCE in comparison to histopathology. The results show OCE imaging could reliably differentiate between benign prostate tissue, acinar atypical hyperplasia, prostatic intraepithelial neoplasia and malignant PCa. The sensitivity and specificity of OCE for the detection of prostate cancer was 0.98 and 0.91 with AUC > 0.99. Quantitative 3D OCE based on the assessment of mechanical properties of tissues can reliably differentiate prostate tissue specimen in an ex-vivo setting. This is a promising imaging modality for characterising different grades of cancers.

In vivo high-resolution magnetic resonance elastography of the uterine corpus and cervix

OBJECTIVES

To apply 3D multifrequency MR elastography (3DMMRE) to the uterus and analyse the viscoelasticity of the uterine tissue in healthy volunteers considering individual variations and variations over the menstrual cycle.

METHODS

Sixteen healthy volunteers participated in the study, one of whom was examined 12 times over two menstrual cycles. Pelvic 3DMMRE was performed on a 1.5-T scanner with seven vibration frequencies (30-60 Hz) using a piezoelectric driver. Two mechanical parameter maps were obtained corresponding to the magnitude (|G (*) |) and the phase angle (φ) of the complex shear modulus.

RESULTS

On average, the uterine corpus had higher elasticity, but similar viscosity compared with the cervix, reflected by |G (*) |uterine corpus = 2.58 ± 0.52 kPa vs. |G (*) |cervix = 2.00 ± 0.34 kPa (p < 0.0001) and φ uterine corpus = 0.54 ± 0.08, φ cervix = 0.57 ± 0.12 (p = 0.428). With 2.23 ± 0.26 kPa, |G (*) | of the myometrium was lower in the secretory phase (SP) compared with that of the proliferative phase (PP, |G (*) | = 3.01 ± 0.26 kPa). For the endometrium, the value of |G (*) | in SP was 68% lower than during PP (PP, |G (*) | = 3.34 ± 0.42 kPa; SP, |G (*) | = 1.97 ± 0.34 kPa; p = 0.0061).

CONCLUSION

3DMMRE produces high-resolution mechanical parameter maps of the uterus and cervix and shows sensitivity to structural and functional changes of the endometrium and myometrium during the menstrual cycle.

KEY POINTS

MR elastography provided for the first time spatially resolved viscoelasticity maps of uterus. Uterine corpus had a higher elasticity, but similar viscosity compared with cervix. The stiffness of both endometrium and myometrium decreases during the menstrual cycle.

Spatially-resolved hydraulic conductivity estimation via poroelastic magnetic resonance elastography

Poroelastic magnetic resonance elastography is an imaging technique that could recover mechanical and hydrodynamical material properties of in vivo tissue. To date, mechanical properties have been estimated while hydrodynamical parameters have been assumed homogeneous with literature-based values. Estimating spatially-varying hydraulic conductivity would likely improve model accuracy and provide new image information related to a tissue's interstitial fluid compartment. A poroelastic model was reformulated to recover hydraulic conductivity with more appropriate fluid-flow boundary conditions. Simulated and physical experiments were conducted to evaluate the accuracy and stability of the inversion algorithm. Simulations were accurate (property errors were < 2%) even in the presence of Gaussian measurement noise up to 3%. The reformulated model significantly decreased variation in the shear modulus estimate (p << 0.001) and eliminated the homogeneity assumption and the need to assign hydraulic conductivity values from literature. Material property contrast was recovered experimentally in three different tofu phantoms and the accuracy was improved through soft-prior regularization. A frequency-dependence in hydraulic conductivity contrast was observed suggesting that fluid-solid interactions may be more prominent at low frequency. In vivo recovery of both structural and hydrodynamical characteristics of tissue could improve detection and diagnosis of neurological disorders such as hydrocephalus and brain tumors.

Shear-wave amplitudes measured with cardiac MR elastography for diagnosis of diastolic dysfunction

PURPOSE

To test whether shear-wave amplitudes (SWAs) in the myocardium measured with cardiac magnetic resonance (MR) elastography enable diagnosis of myocardial relaxation abnormalities in patients with diastolic dysfunction.

MATERIALS AND METHODS

Each subject gave written informed consent to participate in this institutional review board-approved prospective study. Electrocardiographically triggered SWA-based cardiac MR elastography with 24.13-Hz external vibration frequency was performed in 50 subjects grouped into asymptomatic young (n = 10, 18-39 years) and asymptomatic old (n = 10, 40-68 years) subjects and patients with echocardiographically proved mild, moderate, or severe diastolic dysfunction (n = 30, 44-73 years). SWA images were analyzed in the left ventricular (LV) region and were normalized against reference SWA of the thoracic wall. Analysis of variance with Bonferroni-corrected pairwise comparison and Pearson correlation were used for statistical evaluation.

RESULTS

Young and old control subjects had normalized mean LV SWA of 0.67 ± 0.04 (standard error of the mean) and 0.56 ± 0.04 (P = .18, F test), respectively. Compared with the control groups, patients with mild, moderate, and severe diastolic dysfunction displayed significantly reduced normalized mean LV SWA of 0.37 ± 0.04, 0.34 ± 0.04, and 0.29 ± 0.04 (P < .001, F test), respectively, which was inversely correlated to the severity of diastolic dysfunction (R = -0.61, P < .001). The best cutoff value to differentiate between asymptomatic volunteers and patients was 0.43, yielding an area under the receiver operating characteristic curve of 0.92, with 90% sensitivity and 89.7% specificity.

CONCLUSION

LV SWA measured with cardiac MR elastography provides image contrast sensitive to myocardial relaxation abnormalities and shows significantly lower values in patients with diastolic dysfunction.

Isovolumetric elasticity alteration in the human heart detected by in vivo time-harmonic elastography

Time harmonic elastography (THE) has recently been introduced for measurement of the periodic alteration in myocardial shear modulus based on externally induced low-frequency acoustic vibrations produced by a loudspeaker. In this study, we propose further developments of cardiac THE toward a clinical modality including integration of the vibration source into the patient bed and automated parameter extraction from harmonic shear wave amplitudes, wall motion profiles and synchronized electrocardiographic records. This method has enabled us to evaluate the delay between wall motion and wave amplitude alteration for the measurement of isovolumetric times of elasticity alteration during contraction (τ(C)) and relaxation (τ(R)) in a group of 32 healthy volunteers. On average, the wave amplitudes changed between systole and diastole by a factor of 1.7 ± 0.3, with a τ(C) of 137 ± 61 ms and a τ(R) of 68 ± 73 ms, which agrees with results obtained with the more time-consuming and expensive cardiac magnetic resonance elastography. Furthermore, because of the high sampling rate, elasto-morphometric parameters such as transition times and the area of wave amplitude-cardiac motion cycles can be processed in an automated way for the future clinical detection of myocardial relaxation abnormalities.

Multi-parametric MRI as an indirect evaluation tool of the mechanical properties of in-vitro cardiac tissues

Background

Early detection of heart failure is essential to effectively reduce related mortality. The quantification of the mechanical properties of the myocardium, a primordial indicator of the viability of the cardiac tissue, is a key element in patient’s care. Despite an incremental utilization of multi-parametric magnetic resonance imaging (MRI) for cardiac tissue characteristics and function, the link between multi-parametric MRI and the mechanical properties of the heart has not been established. We sought to determine the parametric relationship between the myocardial mechanical properties and the MR parameters. The specific aim was to develop a reproducible evaluative quantitative tool of the mechanical properties of cardiac tissue using multi-parametric MRI associated to principal component analysis.

Methods

Samples from porcine hearts were submitted to a multi-parametric MRI acquisition followed by a uniaxial tensile test. Multi linear regressions were performed between dependent (Young’s modulus E) and independent (relaxation times T1, T2 and T2*, magnetization transfer ratio MTR, apparent diffusion coefficient ADC and fractional anisotropy FA) variables. A principal component analysis was used to convert the set of possibly correlated variables into a set of linearly uncorrelated variables.

Results

Values of 46.1±12.7 MPa for E, 729±21 ms for T1, 61±6 ms for T2, 26±7 for T2*, 35±5% for MTRx100, 33.8±4.7 for FAx10^-2, and 5.85±0.21 mm²/s for ADCx10^-4 were measured. Multi linear regressions showed that only 45% of E can be explained by the MRI parameters. The principal component analysis reduced our seven variables to two principal components with a cumulative variability of 63%, which increased to 80% when considering the third principal component.

Conclusions

The proposed multi-parametric MRI protocol associated to principal component analysis is a promising tool for the evaluation of mechanical properties within the left ventricle in the in vitro porcine model. Our in vitro experiments will now allow us focused in vivo testing on healthy and infracted hearts in order to determine useful quantitative MR-based biomarkers.

Review of MR elastography applications and recent developments

The technique of MR elastography (MRE) has emerged as a useful modality for quantitatively imaging the mechanical properties of soft tissues in vivo. Recently, MRE has been introduced as a clinical tool for evaluating chronic liver disease, but many other potential applications are being explored. These applications include measuring tissue changes associated with diseases of the liver, breast, brain, heart, and skeletal muscle including both focal lesions (e.g., hepatic, breast, and brain tumors) and diffuse diseases (e.g., fibrosis and multiple sclerosis). The purpose of this review article is to summarize some of the recent developments of MRE and to highlight some emerging applications.

Magnetic resonance elastography as a method to estimate myocardial contractility

PURPOSE

To determine whether increasing epinephrine infusion in an in vivo pig model is associated with an increase in end-systolic magnetic resonance elastography (MRE)-derived effective stiffness.

MATERIALS AND METHODS

Finite element modeling (FEM) was performed to determine the range of myocardial wall thicknesses that could be used for analysis. Then MRE was performed on five pigs to measure the end-systolic effective stiffness with epinephrine infusion. Epinephrine was continuously infused intravenously in each pig to increase the heart rate in increments of 20%. For each such increase end-systolic effective stiffness was measured using MRE. In each pig, Student's t-test was used to compare effective end-systolic stiffness at baseline and at initial infusion of epinephrine. Least-square linear regression was performed to determine the correlation between normalized end-systolic effective stiffness and increase in heart rate with epinephrine infusion.

RESULTS

FEM showed that phase gradient inversion could be performed on wall thickness ≈≥1.5 cm. In pigs, effective end-systolic stiffness significantly increased from baseline to the first infusion in all pigs (P = 0.047). A linear correlation was found between normalized effective end-systolic stiffness and percent increase in heart rate by epinephrine infusion with R(2) ranging from 0.86-0.99 in four pigs. In one of the pigs the R(2) value was 0.1. A linear correlation with R(2) = 0.58 was found between normalized effective end-systolic stiffness and percent increase in heart rate when pooling data points from all pigs.

CONCLUSION

Noninvasive MRE-derived end-systolic effective myocardial stiffness may be a surrogate for myocardial contractility.

Evaluating elastic properties of heterogeneous soft tissue by surface acoustic waves detected by phase-sensitive optical coherence tomography

The combined use of surface acoustic wave (SAW) and phase-sensitive optical coherence tomography (PhS-OCT) is useful to evaluate the elasticity of layered biological tissues, such as normal skin. However, the pathological tissue is often originated locally, leading to the alternation of mechanical properties along both axial and lateral directions. We present a feasibility study on whether the SAW technique is sensitive to detect the alternation of mechanical property along the lateral direction within tissue, which is important for clinical utility of this technique to localize diseased tissue. Experiments are carried out on purposely designed tissue phantoms and ex vivo chicken breast samples, simulating the localized change of elasticity. A PhS-OCT system is employed not only to provide the ultra-high sensitive measurement of the generated surface waves on the tissue surface, but also to provide the real time imaging of the tissue to assist the elasticity evaluation of the heterogeneous tissue. The experimental results demonstrate that with PhS-OCT used as a pressure sensor, the SAW is highly sensitive to the elasticity change of the specimen in both vertical and lateral directions with a sensing depth of ∼5 mm with our current system setup, thus promising its useful clinical applications where the quantitative elasticity of localized skin diseases is needed to aid in diagnosis and treatment.

Quantitative elastography provided by surface acoustic waves measured by phase-sensitive optical coherence tomography

We report on a quantitative elastography technique achieved by combining phase-sensitive optical coherence tomography (PhS-OCT) with the surface acoustic wave (SAW) method. Different from traditional optical coherence elastography, the elastography is achieved by impulse-stimulated SAW, rather than by shear waves. PhS-OCT serves not only as a detector to measure SAW signals but also as a means to provide a cross-sectional image of the sample. The experimental results indicate that the combination of PhS-OCT with SAW is feasible to provide quantitative elastography of heterogeneous tissue samples.

In Vivo time harmonic elastography of the human heart

Time harmonic elastography is introduced as a modality for assessing myocardial elasticity changes during the cardiac cycle. It is based on external stimulation and real-time analysis of 30-Hz harmonic shear waves in axial direction of a parasternal line of sight through the lateral heart wall. In 20 healthy volunteers, the externally induced waves showed smaller amplitudes during systole (76.0 ± 30.8 μm) and higher amplitudes during diastole (126.7 ± 52.1 μm). This periodic wave amplitude alteration preceded ventricular contraction and dilation by about 100 ms. The amplitude ratio of 1.75 ± 0.49 indicates a relative change in myocardial shear elasticity on the order of 14 ± 11. These results well agree with observations made by cardiac magnetic resonance elastography for a similar displacement component and region of the heart. The proposed method provides reproducible elastodynamic information on the heart in real-time and may help in diagnosing myocardial relaxation abnormalities in the future.

On Lamb and Rayleigh wave convergence in viscoelastic tissues

Characterization of the viscoelastic material properties of soft tissue has become an important area of research over the last two decades. Our group has been investigating the feasibility of using a shear wave dispersion ultrasound vibrometry (SDUV) method to excite Lamb waves in organs with plate-like geometry to estimate the viscoelasticity of the medium of interest. The use of Lamb wave dispersion ultrasound vibrometry to quantify the mechanical properties of viscoelastic solids has previously been reported. Two organs, the heart wall and the spleen, can be readily modeled using plate-like geometries. The elasticity of these two organs is important because they change in pathological conditions. Diastolic dysfunction is the inability of the left ventricle (LV) of the heart to supply sufficient stroke volumes into the systemic circulation and is accompanied by the loss of compliance and stiffening of the LV myocardium. It has been shown that there is a correlation between high splenic stiffness in patients with chronic liver disease and strong correlation between spleen and liver stiffness. Here, we investigate the use of the SDUV method to quantify the viscoelasticity of the LV free-wall myocardium and spleen by exciting Rayleigh waves on the organ's surface and measuring the wave dispersion (change of wave velocity as a function of frequency) in the frequency range 40–500 Hz. An equation for Rayleigh wave dispersion due to cylindrical excitation was derived by modeling the excised myocardium and spleen with a homogenous Voigt material plate immersed in a nonviscous fluid. Boundary conditions and wave potential functions were solved for the surface wave velocity. Analytical and experimental convergence between the Lamb and Rayleigh waves is reported in a finite element model of a plate in a fluid of similar density, gelatin plate and excised porcine spleen and left-ventricular free-wall myocardium.

AN OVERVIEW OF ELASTOGRAPHY - AN EMERGING BRANCH OF MEDICAL IMAGING

From times immemorial manual palpation served as a source of information on the state of soft tissues and allowed detection of various diseases accompanied by changes in tissue elasticity. During the last two decades, the ancient art of palpation gained new life due to numerous emerging elasticity imaging (EI) methods. Areas of applications of EI in medical diagnostics and treatment monitoring are steadily expanding. Elasticity imaging methods are emerging as commercial applications, a true testament to the progress and importance of the field.In this paper we present a brief history and theoretical basis of EI, describe various techniques of EI and, analyze their advantages and limitations, and overview main clinical applications. We present a classification of elasticity measurement and imaging techniques based on the methods used for generating a stress in the tissue (external mechanical force, internal ultrasound radiation force, or an internal endogenous force), and measurement of the tissue response. The measurement method can be performed using differing physical principles including magnetic resonance imaging (MRI), ultrasound imaging, X-ray imaging, optical and acoustic signals.Until recently, EI was largely a research method used by a few select institutions having the special equipment needed to perform the studies. Since 2005 however, increasing numbers of mainstream manufacturers have added EI to their ultrasound systems so that today the majority of manufacturers offer some sort of Elastography or tissue stiffness imaging on their clinical systems. Now it is safe to say that some sort of elasticity imaging may be performed on virtually all types of focal and diffuse disease. Most of the new applications are still in the early stages of research, but a few are becoming common applications in clinical practice.

Cardiac magnetic resonance elastography: toward the diagnosis of abnormal myocardial relaxation

AIM

To assess the potential of cardiac magnetic resonance elastography (MRE) for elasticity-based detection of abnormal left ventricular (LV) relaxation.

MATERIALS AND METHODS

Cardiac MRE was performed in 3 groups: young volunteers (n = 11; mean age, 31.7 years), older volunteers (n = 5; mean age, 54.8 years), and a group with relaxation abnormalities (n = 11; mean age, 58 years) identified by transthoracic echocardiography. Cine MR imaging served to measure LV volumes and global LV systolic function. Wave-amplitude-sensitive electrocardiograph-gated steady-state MRE was performed using an extended piston driver attached to the anterior chest wall. Phase contrast shear wave images were acquired in all 3 Cartesian components and combined to generate amplitude maps. This was done using the time-gradient operator for linear high-pass filtering and phase unwrapping followed by temporal Fourier transformation for extracting externally induced 24.13-Hz shear oscillations from intrinsic motion and blood flow. Amplitudes were evaluated in the left ventricle and normalized by wave amplitudes outside the heart, adjacent to the right ventricle.

RESULTS

One patient and 1 young volunteer had to be excluded from final analysis because of considerable body movement during the acquisition of the MRE scans. Mean wave amplitudes in the remaining subjects were 0.22 ± 0.05 mm in young volunteers, 0.23 ± 0.09 in older volunteers, and 0.14 ± 0.03 mm in patients. The mean ratio of amplitudes inside the ventricle to the anterior chest wall was 0.62 ± 0.15 for young volunteers, 0.50 ± 0.09 for older volunteers, and 0.33 ± 0.08 for patients.

CONCLUSION

MRE identifies significantly reduced LV shear wave amplitudes in patients with mild relaxation abnormality. Thus, cardiac MRE provides a promising modality for an elasticity-based diagnosis of dysfunctional myocardial relaxation.

Review of journal of cardiovascular magnetic resonance 2010

There were 75 articles published in the Journal of Cardiovascular Magnetic Resonance (JCMR) in 2010, which is a 34% increase in the number of articles since 2009. The quality of the submissions continues to increase, and the editors were delighted with the recent announcement of the JCMR Impact Factor of 4.33 which showed a 90% increase since last year. Our acceptance rate is approximately 30%, but has been falling as the number of articles being submitted has been increasing. In accordance with Open-Access publishing, the JCMR articles go on-line as they are accepted with no collating of the articles into sections or special thematic issues. Last year for the first time, the Editors summarized the papers for the readership into broad areas of interest or theme, which we felt would be useful to practitioners of cardiovascular magnetic resonance (CMR) so that you could review areas of interest from the previous year in a single article in relation to each other and other recent JCMR articles 1. This experiment proved very popular with a very high rate of downloading, and therefore we intend to continue this review annually. The papers are presented in themes and comparison is drawn with previously published JCMR papers to identify the continuity of thought and publication in the journal. We hope that you find the open-access system increases wider reading and citation of your papers, and that you will continue to send your quality manuscripts to JCMR for publication.

In vivo assessment of MR elastography-derived effective end-diastolic myocardial stiffness under different loading conditions

PURPOSE

To compare magnetic resonance elastography (MRE) effective stiffness to end-diastolic pressure at different loading conditions to demonstrate a relationship between myocardial MRE effective stiffness and end-diastolic left ventricular (LV) pressure.

MATERIALS AND METHODS

MRE was performed on four pigs to measure the end-diastolic effective stiffness under different loading conditions. End-diastolic pressure was increased by infusing Dextran-40 (20% of blood volume). For each infusion of Dextran-40, end-diastolic pressure was recorded and end-diastolic effective stiffness was measured using MRE. In each pig, least-square linear regression was performed to determine the correlation between end-diastolic effective stiffness and end-diastolic LV pressure.

RESULTS

A linear correlation was found between end-diastolic LV pressure and end-diastolic effective stiffness with R(2) ranging from 0.73-0.9. A linear correlation with R(2) = 0.26 was found between end-diastolic LV pressure and end-diastolic effective stiffness when pooling data points from all pigs.

CONCLUSION

End-diastolic effective myocardial stiffness increases linearly with end-diastolic LV pressure.

Elasticity-based determination of isovolumetric phases in the human heart

BACKGROUND/MOTIVATION: To directly determine isovolumetric cardiac time intervals by magnetic resonance elastography (MRE) using the magnitude of the complex signal for deducing morphological information combined with the phase of the complex signal for tension-relaxation measurements.

METHODS

Thirty-five healthy volunteers and 11 patients with relaxation abnormalities were subjected to transthoracic wave stimulation using vibrations of approximately 25 Hz. A k-space-segmented, ECG-gated gradient-recalled echo steady-state sequence with a 500-Hz bipolar motion-encoding gradient was used for acquiring a series of 360 complex images of a short-axis view of the heart at a frame rate of less than 5.2 ms. Magnitude images were employed for measuring the cross-sectional area of the left ventricle, while phase images were used for analyzing the amplitudes of the externally induced waves. The delay between the decrease in amplitude and onset of ventricular contraction was determined in all subjects and assigned to the time of isovolumetric tension. Conversely, the delay between the increase in wave amplitude and ventricular dilatation was used for measuring the time of isovolumetric elasticity relaxation.

RESULTS

Wave amplitudes decreased during systole and increased during diastole. The variation in wave amplitude occurred ahead of morphological changes. In healthy volunteers the time of isovolumetric elasticity relaxation was 75 ± 31 ms, which is significantly shorter than the time of isovolumetric tension of 136 ± 36 ms (P < 0.01). In patients with relaxation abnormalities (mild diastolic dysfunction, n = 11) isovolumetric elasticity relaxation was significantly prolonged, with 133 ± 57 ms (P < 0.01), whereas isovolumetric tension time was in the range of healthy controls (161 ± 45 ms; P = 0.053).

CONCLUSION

The complex MRE signal conveys complementary information on cardiac morphology and elasticity, which can be combined for directly measuring isovolumetric tension and elasticity relaxation in the human heart.

Review of Journal of Cardiovascular Magnetic Resonance 2009

There were 56 articles published in the Journal of Cardiovascular Magnetic Resonance in 2009. The editors were impressed with the high quality of the submissions, of which our acceptance rate was about 40%. In accordance with open-access publishing, the articles go on-line as they are accepted with no collating of the articles into sections or special thematic issues. We have therefore chosen to briefly summarise the papers in this article for quick reference for our readers in broad areas of interest, which we feel will be useful to practitioners of cardiovascular magnetic resonance (CMR). In some cases where it is considered useful, the articles are also put into the wider context with a short narrative and recent CMR references. It has been a privilege to serve as the Editor of the JCMR this past year. I hope that you find the open-access system increases wider reading and citation of your papers, and that you will continue to send your quality manuscripts to JCMR for publication.