Thigh muscle stiffness assessed with magnetic resonance elastography in hyperthyroid patients before and after medical treatment

Analysis of thigh muscle stiffness from childhood to adulthood using magnetic resonance elastography (MRE) technique

BACKGROUND

Magnetic resonance elastography has been performed in healthy and pathological muscles in order to provide clinicians with quantitative muscle stiffness data. However, there is a lack of data on pediatric muscle. Therefore, the present work studies age-related changes of the mechanical properties.

METHODS

26 healthy subjects composed of 7 children (8-12 years), 9 young adults (24-29 years) and 10 middle-aged adults (53-58 years) underwent a magnetic resonance elastography test. Shear modulus (μ) and its spatial distribution, as well as the attenuation coefficient (α) were measured on the vastus medialis muscle at rest and at contracted conditions (10% and 20% of the maximum voluntary contraction) for each group.

FINDINGS

The shear modulus linearly increases with the degree of contraction for young adults while it is maximum at 10% of the maximum voluntary contraction for children (μ_(children_10%)=14.9kPa (SD 2.18)) and middle-aged adults (μ_(middle-aged_10%)=10.42kPa (SD 1.38)). Mapping of shear modulus revealed a diffuse distribution of colors reflecting differences in muscle physiological activity as a function of age. The attenuation coefficient showed a similar behavior for all groups, i.e. a decrease from the relaxed to the contracted states.

INTERPRETATION

This study demonstrates that the magnetic resonance elastography technique is sensitive enough to detect changes in muscle mechanical properties for children, middle-aged and young adults and could provide clinicians with a muscle reference data base as a function of age, improving the diagnosis of muscular dystrophy.

Calculation of shear stiffness in noise dominated magnetic resonance elastography data based on principal frequency estimation

Magnetic resonance elastography (MRE) is a non-invasive phase-contrast-based method for quantifying the shear stiffness of biological tissues. Synchronous application of a shear wave source and motion encoding gradient waveforms within the MRE pulse sequence enable visualization of the propagating shear wave throughout the medium under investigation. Encoded shear wave-induced displacements are then processed to calculate the local shear stiffness of each voxel. An important consideration in local shear stiffness estimates is that the algorithms employed typically calculate shear stiffness using relatively high signal-to-noise ratio (SNR) MRE images and have difficulties at an extremely low SNR. A new method of estimating shear stiffness based on the principal spatial frequency of the shear wave displacement map is presented. Finite element simulations were performed to assess the relative insensitivity of this approach to decreases in SNR. Additionally, ex vivo experiments were conducted on normal rat lungs to assess the robustness of this approach in low SNR biological tissue. Simulation and experimental results indicate that calculation of shear stiffness by the principal frequency method is less sensitive to extremely low SNR than previously reported MRE inversion methods but at the expense of loss of spatial information within the region of interest from which the principal frequency estimate is derived.

A mechano-acoustic indentor system for in vivo measurement of nonlinear elastic properties of soft tissue

OBJECTIVES

Soft tissue exhibits nonlinear stress-strain behavior under compression. Characterizing its nonlinear elasticity may aid detection, diagnosis, and treatment of soft tissue abnormality. The purposes of this study were to develop a rate-controlled Mechano-Acoustic Indentor System and a corresponding finite element optimization method to extract nonlinear elastic parameters of soft tissue and evaluate its test-retest reliability.

METHODS

An indentor system using a linear actuator to drive a force-sensitive probe with a tip-mounted ultrasound transducer was developed. Twenty independent sites at the upper lateral quadrant of the buttock from 11 asymptomatic subjects (7 men and 4 women from a chiropractic college) were indented at 6% per second for 3 sessions, each consisting of 5 trials. Tissue thickness, force at 25% deformation, and area under the load-deformation curve from 0% to 25% deformation were calculated. Optimized hyperelastic parameters of the soft tissue were calculated with a finite element model using a first-order Ogden material model. Load-deformation response on a standardized block was then simulated, and the corresponding area and force parameters were calculated. Between-trials repeatability and test-retest reliability of each parameter were evaluated using coefficients of variation and intraclass correlation coefficients, respectively.

RESULTS

Load-deformation responses were highly reproducible under repeated measurements. Coefficients of variation of tissue thickness, area under the load-deformation curve from 0% to 25% deformation, and force at 25% deformation averaged 0.51%, 2.31%, and 2.23%, respectively. Intraclass correlation coefficients ranged between 0.959 and 0.999, indicating excellent test-retest reliability.

CONCLUSIONS

The automated Mechano-Acoustic Indentor System and its corresponding optimization technique offers a viable technology to make in vivo measurement of the nonlinear elastic properties of soft tissue. This technology showed excellent between-trials repeatability and test-retest reliability with potential to quantify the effects of a wide variety of manual therapy techniques on the soft tissue elastic properties.

A method for assessing the fit of a constitutive material model to experimental stress-strain data

Higher-order polynomial functions can be used as a constitutive model to represent the mechanical behaviour of biological materials. The goal of this study was to present a method for assessing the fit of a given constitutive three-dimensional material model. Goodness of fit was assessed using multiple parameters including the root mean square error and Hotelling's T 2-test. Specifically, a polynomial model was used to characterise the stress-strain data, varying the number of model terms used (45 combinations of between 3 and 11 terms) and the manner of optimisation used to establish model coefficients (i.e. determining coefficients either by parameterisation of all data simultaneously or averaging coefficients obtained by parameterising individual data trials). This framework for model fitting helps to ensure that a given constitutive formulation provides the best characterisation of biological material mechanics.

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.

Elastography: modality-specific approaches, clinical applications, and research horizons

Manual palpation has been used for centuries to provide a relative indication of tissue health and disease. Engineers have sought to make these assessments increasingly quantitative and accessible within daily clinical practice. Since many of the developed techniques involve image-based quantification of tissue deformation in response to an applied force (i.e., "elastography"), such approaches fall squarely within the domain of the radiologist. While commercial elastography analysis software is becoming increasingly available for clinical use, the internal workings of these packages often remain a "black box," with limited guidance on how to usefully apply the methods toward a meaningful diagnosis. The purpose of the present review article is to introduce some important approaches to elastography that have been developed for the most widely used clinical imaging modalities (e.g., ultrasound, MRI), to provide a basic sense of the underlying physical principles, and to discuss both current and potential (musculoskeletal) applications. The article also seeks to provide a perspective on emerging approaches that are rapidly developing in the research laboratory (e.g., optical coherence tomography, fibered confocal microscopy), and which may eventually gain a clinical foothold.

Evaluation of muscles affected by myositis using magnetic resonance elastography

INTRODUCTION

Idiopathic inflammatory myopathies (IIMs, or myositis) represent a group of autoimmune diseases that result in decreased muscle strength and/or endurance. Non-invasive tools to assess muscle may improve our understanding of the clinical and functional consequences of myopathies and their response to treatment. In this study we examine magnetic resonance elastography (MRE), a non-invasive technique that assesses the shear modulus (stiffness) of muscle, in IIM subjects.

METHODS

Nine subjects with active myositis completed the MRE protocol. Participants lay in a positioning device, and scans of the vastus medialis (VM) were taken in the relaxed state and at two contraction levels. Manual inversion was used to estimate the stiffness.

RESULTS

A significant reduction in muscle stiffness was seen in myositis subjects compared with healthy controls during the "relaxed" condition.

DISCUSSION

The use of non-invasive technologies such as MRE may provide greater understanding of the pathophysiology of IIM and improve assessment of treatment efficacy.

Characterization of muscle architecture in children and adults using magnetic resonance elastography and ultrasound techniques

The purpose of this study is to characterize the muscle architecture of children and adults using magnetic resonance elastography and ultrasound techniques. Five children (8-12 yr) and seven adults (24-58 yr) underwent both tests on the vastus medialis muscle at relaxed and contracted (10% and 20% of MVC) states. Longitudinal ultrasonic images were performed in the same area as the phase image showing the shear wave's propagation. Two geometrical parameters were defined: the wave angle (α(_MRE)) corresponding to the shear wave propagation and the fascicule angle (α(_US)) tracking the path of fascicles. Moreover, shear modulus was measured at different localizations within the muscle and in the subcutaneous adipose tissue. The association of both techniques demonstrates that the shear wave propagation follows the muscle fascicles path, reflecting the internal muscle architecture. At rest, ultrasound images revealed waves propagating parallel to the children fascicle while adults showed oblique waves corresponding to already oriented (α(_US)=15.4±2.54°) muscle fascicles. In contraction, the waves' propagation were in an oblique direction for children (α(_US_10%MVC)=10.6±2.27°, α(_US_20%MVC)=10.2±2.29°) as well as adults (α(_US_10%MVC)=15.4±2.54°, α(_US_20%MVC)=17.2±2.44°). A stiffness variation (1 kPa) was found between the upper and lower parts of the adult VM muscle and a lower stiffness (1.85±0.17 kPa) was measured in the subcutaneous adipose tissue. This study demonstrates the feasibility of the MRE technique to provide geometrical insights from the children and adults muscles and to characterize different physiological media.

Viscoelasticity-based MR elastography of skeletal muscle

An in vivo multifrequency magnetic resonance elastography (MRE) protocol was developed for studying the viscoelastic properties of human skeletal muscle in different states of contraction. Low-frequency shear vibrations in the range of 25-62.5 Hz were synchronously induced into the femoral muscles of seven volunteers and measured in a cross-sectional view by encoding the fast-transverse shear wave component parallel to the muscle fibers. The so-called springpot model was used for deriving two viscoelastic constants, μ and α, from the dispersion functions of the complex shear modulus in relaxed and in loaded muscle. Representing the shear elasticity parallel to the muscle fibers, μ increased in all volunteers upon contraction from 2.68 ± 0.23 kPa to 3.87 ± 0.50 kPa. Also α varied with load, indicating a change in the geometry of the mechanical network of muscle from relaxation (α = 0.253 ± 0.009) to contraction (α = 0.270 ± 0.009). These results provide a reference for a future assessment of muscular dysfunction using rheological parameters.

Magnetic resonance elastography as a method for the assessment of effective myocardial stiffness throughout the cardiac cycle

MR elastography (MRE) is a noninvasive technique in which images of externally generated waves propagating in tissue are used to measure stiffness. The first aim is to determine, from a range of driver configurations, the optimal driver for the purpose of generating waves within the heart in vivo. The second aim is to quantify the shear stiffness of normal myocardium throughout the cardiac cycle using MRE and to compare MRE stiffness to left ventricular chamber pressure in an in vivo pig model. MRE was performed in six pigs with six different driver setups, including no motion, three noninvasive drivers, and two invasive drivers. MRE wave displacement amplitudes were calculated for each driver. During the same MRI examination, left ventricular pressure and MRI-measured left ventricular volume were obtained, and MRE myocardial stiffness was calculated for 20 phases of the cardiac cycle. No discernible waves were imaged when no external motion was applied, and a single pneumatic drum driver produced higher amplitude waves than the other noninvasive drivers (P < 0.05). Pressure-volume loops overlaid onto stiffness-volume loops showed good visual agreement. Pressure and MRE-measured effective stiffness showed good correlation (R(2) = 0.84). MRE shows potential as a noninvasive method for estimating effective myocardial stiffness throughout the cardiac cycle.

Endosonographic elastography of the anal sphincter in patients with fecal incontinence

OBJECTIVE

In fecal incontinence the role of elastography has not yet been evaluated. We performed a trial to further characterize the internal and external anal sphincter in patients with fecal incontinence and compared a visual assessment scale with a computerized program for quantifying elastic properties of the anal sphincter.

MATERIAL AND METHODS

Fifty consecutive patients with fecal incontinence were studied (n = 31 following lower anterior resection, n = 8 with Crohn's disease, n = 9 following colon surgery, n = 2 others). Elastogram color distribution within the sphincter representing elastic properties was quantified using a visual analog scale and an off-line computerized area calculation program.

RESULTS

The main finding was that the inner anal sphincter (IAS) differed significantly from the external anal sphincter (EAS) with regard to elastogram color distribution. There were no significant correlations with clinical and functional parameters. There was, however, a non-significant increase in the percentage of blue (hard) areas in the IAS in patients neoadjuvantly irradiated for rectal or cervical cancer compared to non-irradiated patients, which was accompanied by a significant decrease in the resting sphincter pressure (p < 0.009).

CONCLUSIONS

The IAS, a smooth muscle, and the EAS, a striated muscle, have different elastogram color distributions, probably reflecting their different elastic properties. The absence of significant correlations with the major clinical and functional parameters suggests that in routine clinical practice ultrasound real-time elastography may not yield additional information in patients with fecal incontinence. There may be exceptions, particularly in irradiated patients.

Myositis associated with the decline of thyroid hormone levels in thyrotoxicosis: a syndrome?

BACKGROUND

Musculoskeletal complaints are common in patients with thyroid dysfunction. Both thyrotoxic and hypothyroid myopathy have been well described, and there are distinct presentations, laboratory findings, and clinical outcomes between the two groups. Myopathy has also been reported in hyperthyroid patients only after beginning treatment, suggesting that relative hypothyroidism may also contribute to musculoskeletal disease. A confounding factor in these cases was that these patients were on antithyroid drugs that may also have direct effects on the muscle, irrespective of the rate of decline in thyroid hormone levels.

SUMMARY

We report a patient with Graves' disease who developed myalgias with elevated creatine kinase levels after total thyroidectomy. Addition of triiodothyronine quickly resolved her symptoms and creatine kinase levels, whereas discontinuation of triiodothyronine, despite having normal to elevated total thyroxine levels, led to a relapse.

CONCLUSION

Myositis after correction of thyrotoxicosis may constitute a syndrome that should be assessed for in hyperthyroid patients complaining of myalgias after starting treatment.

Development and application of magnetic resonance elastography of the normal and pathological thyroid gland in vivo

PURPOSE

To noninvasively assess the shear stiffness of the thyroid gland in vivo in order to determine whether magnetic resonance elastography (MRE) might hold clinical utility in the diagnosis of thyroid disease.

MATERIALS AND METHODS

Quantitative parametric images of thyroid stiffness in normal volunteers and patients were produced and quantitative stiffness values measured. Average gland stiffness was determined by region of interest analysis of the parametric images. This technique was used to assess stiffness of the thyroid in normal individuals (n = 12), patients with Hashimoto's thyroiditis (HT; n = 5), and patients with a solitary benign (n = 8) or malignant (n = 2) thyroid nodule.

RESULTS

Mean shear modulus of normal thyroid glands was 1.9 +/- 0.6 kPa at 100 Hz and 1.3 +/- 0.5 kPa at 80 Hz, while that of HT glands was 2.8 +/- 0.6 kPa and 1.8 +/- 0.6 kPa at 80 Hz, respectively (P = 0.004 at 100 Hz). Elastographic parameters could not differentiate benign from malignant thyroid nodules in these small sample sizes.

CONCLUSION

We developed a method for the application of MRE to the study of thyroid gland pathology. The results show that the HT gland can be differentiated from normal thyroid. The clinical utility of this imaging modality in the diagnosis and management of thyroid disease awaits further study.

Evaluation of a rapid, multiphase MRE sequence in a heart-simulating phantom

The aims of this study were to validate stiffness estimates of a phantom undergoing cyclic deformation obtained using a multiphase magnetic resonance elastography (MRE) imaging sequence by comparison with those obtained using a single-phase MRE sequence and to quantify the stability of the multiphase-derived stiffness estimates as a function of deformation frequency and imaging parameters. A spherical rubber shell of 10 cm diameter and 1 cm thickness was connected to a computerized flow pump to produce cyclic pressure variations within the phantom. The phantom was imaged at cyclic pressures between 18-72 bpm using single-phase and multiphase MRE acquisitions. The shear stiffness of the phantom was resolved using a spherical shell wave inversion algorithm. Shear stiffness was averaged over the slice of interest and plotted against pressure within the phantom. A linear correlation was observed between stiffness and pressure. Good correlation (R(2) = 0.98) was observed between the stiffness estimates obtained using the standard single-phase and the multiphase pulse sequences. Stiffness estimates obtained using multiphase MRE were stable when the fraction of the deformation period required for acquisition of a single image was not greater than 42%. The results demonstrate the potential of multiphase MRE technique for imaging dynamic organs, such as the heart.

Wave attenuation as a measure of muscle quality as measured by magnetic resonance elastography: initial results

Advances in imaging technologies such as magnetic resonance elastography (MRE) have allowed researchers to gain insights into muscle function in vivo. MRE has been used to examine healthy and diseased muscle by calculating shear modulus. However, additional information can be measured from visualizing a mechanical wave as it passes through a tissue. One such measurable quantity is wave attenuation. The purpose of this study was to determine if a simple measure of wave attenuation could be used to distinguish between healthy and diseased muscle. Twenty seven subjects (14 healthy controls, 7 hyperthyroid myopathy patients, 6 myositis patients) participated in this study. Wave amplitude was determined along a linear profile through the center of the muscle, and an exponential decay curve was fit to the data. This measure was able to find significant differences in attenuation between healthy and diseased muscle. Furthermore, four hyperthyroid myopathy subjects who were tested following treatment all showed improvement by this measure. A likely reason for patients with hyperthyroid myopathy and myositis behaving similarly is that this measurement may reflect similar changes in the muscle extracellular matrix. In addition to modulus, attenuation seems to be an important parameter to measure in skeletal muscle. Further research is needed to investigate other potential measures of attenuation as well as examining other potential measures that can be found from visualizing wave propagation. Future studies should also include muscle biopsies to confirm that the changes seen are as a result of changes in extracellular matrix structure.

Measurement of liver stiffness with two imaging techniques: magnetic resonance elastography and ultrasound elastometry

PURPOSE

To cross-validate the magnetic resonance elastography (MRE) technique with a clinical device, based on an ultrasound elastometry system called Fibroscan.

MATERIALS AND METHODS

Ten healthy subjects underwent an MRE and a Fibroscan test. The MRE technique used a round pneumatic driver at 60 Hz to generate shear waves inside the liver. An elastogram representing a map of the liver stiffness was generated allowing for the measurement of the average liver stiffness inside a region of interest. The Fibroscan technique used an ultrasound probe (3.5 MHz) composed of a vibrator that sent low-frequency (50 Hz) shear waves inside the right liver lobe. The probe acts as an emitter-receptor that measures the velocity of the waves propagated inside the liver tissue.

RESULTS

The mean shear stiffness measured with the MRE and Fibroscan techniques were 1.95+/-0.06 kPa and 1.79+/-0.30 kPa, respectively. A higher standard deviation was found for the same subject with Fibroscan.

CONCLUSION

This study shows why MRE should be investigated beyond the Fibroscan. The MRE technique provided elasticity of the entire liver, meanwhile the Fibroscan provided values of elasticity locally.

Longitudinal strain estimation in incompressible cylindrical tissues from magnetic resonance imaging

In this paper, we present a simple approach for estimating the average longitudinal strains from models reconstructed from medical images. It can be used for many incompressible generalized cylindrical tissues, such as tendons, ligaments, and fusiform muscles; the major deformation directions of these soft tissues are along the longitudinal axes. The method is especially useful when pre- and post-deformation tissue correspondences are difficult to establish directly from images for various reasons, such as insufficient image resolution, homogenous image intensity, and noise. Incompressibility, which is accepted as a good approximation for soft tissues, is exploited as a constraint on the tissue deformation. Experiments with Magnetic Resonance Imaging (MRI) of tissue phantoms and computer simulations show that the method is accurate and practical even in the presence of noise. Finally, we demonstrate the usefulness of our approach on studying extraocular muscle deformation.

Quantitative sonoelastography for the in vivo assessment of skeletal muscle viscoelasticity

A novel quantitative sonoelastography technique for assessing the viscoelastic properties of skeletal muscle tissue was developed. Slowly propagating shear wave interference patterns (termed crawling waves) were generated using a two-source configuration vibrating normal to the surface. Theoretical models predict crawling wave displacement fields, which were validated through phantom studies. In experiments, a viscoelastic model was fit to dispersive shear wave speed sonoelastographic data using nonlinear least-squares techniques to determine frequency-independent shear modulus and viscosity estimates. Shear modulus estimates derived using the viscoelastic model were in agreement with that obtained by mechanical testing on phantom samples. Preliminary sonoelastographic data acquired in healthy human skeletal muscles confirm that high-quality quantitative elasticity data can be acquired in vivo. Studies on relaxed muscle indicate discernible differences in both shear modulus and viscosity estimates between different skeletal muscle groups. Investigations into the dynamic viscoelastic properties of (healthy) human skeletal muscles revealed that voluntarily contracted muscles exhibit considerable increases in both shear modulus and viscosity estimates as compared to the relaxed state. Overall, preliminary results are encouraging and quantitative sonoelastography may prove clinically feasible for in vivo characterization of the dynamic viscoelastic properties of human skeletal muscle.

Ability of magnetic resonance elastography to assess taut bands

BACKGROUND

Myofascial taut bands are central to diagnosis of myofascial pain. Despite their importance, we still lack either a laboratory test or imaging technique capable of objectively confirming either their nature or location. This study explores the ability of magnetic resonance elastography to localize and investigate the mechanical properties of myofascial taut bands on the basis of their effects on shear wave propagation.

METHODS

This study was conducted in three phases. The first involved the imaging of taut bands in gel phantoms, the second a finite element modeling of the phantom experiment, and the third a preliminary evaluation involving eight human subjects-four of whom had, and four of whom did not have myofascial pain. Experiments were performed with a 1.5 T magnetic resonance imaging scanner. Shear wave propagation was imaged and shear stiffness was reconstructed using matched filtering stiffness inversion algorithms.

FINDINGS

The gel phantom imaging and finite element calculation experiments supported our hypothesis that taut bands can be imaged based on its outstanding shear stiffness. The preliminary human study showed a statistically significant 50-100% (P=0.01) increase of shear stiffness in the taut band regions of the involved subjects relative to that of the controls or in nearby uninvolved muscle.

INTERPRETATION

This study suggests that magnetic resonance elastography may have a potential for objectively characterizing myofascial taut bands that have been up to now detectable only by the clinician's fingers.