Search

Introduction: Congenital hypoplasia of the distal half of the tunica albuginea has not been previously described. Aim: To review a patient presenting with erectile dysfunction secondary to congenital penile hypoplasia. Methods: History, physical exam and penile Doppler ultrasound of the patient, followed by a discussed of treatment options and definitive management. Results: Successful operative treatment of our patient with insertion of an inflatable penile prosthesis. Conclusion: We present a case of congenital hypoplasia of the distal tunica albuginea and a successful treatment strategy. We highlight the need for further study of penile embryology.Punjani N, McGarry P, Brock G. Successful Use of an Inflatable Penile Prosthesis for the Treatment of Distal Deficiency of the Tunica Albuginea and Cavernous Tissue. Sex Med 2018;6:356–359. Key Words: Penile Hypoplasia, Erectile Dysfunction, Penile Embryology

Magnetic resonance elastography (MRE) is a phase contrast MRI technique which uses external palpation to create maps of brain mechanical properties noninvasively and in vivo. These mechanical properties are sensitive to tissue microstructure and reflect tissue integrity. MRE has been used extensively to study aging and neurodegeneration, and to assess individual cognitive differences in adults, but little is known about mechanical properties of the pediatric brain. Here we use high-resolution MRE imaging in participants of ages ranging from childhood to adulthood to understand brain mechanical properties across brain maturation. We find that brain mechanical properties differ considerably between childhood and adulthood, and that neuroanatomical subregions have differing maturational trajectories. Overall, we observe lower brain stiffness and greater brain damping ratio with increasing age from 5 to 35 years. Gray and white matter change differently during maturation, with larger changes occurring in gray matter for both stiffness and damping ratio. We also found that subregions of cortical and subcortical gray matter change differently, with the caudate and thalamus changing the most with age in both stiffness and damping ratio, while cortical subregions have different relationships with age, even between neighboring regions. Understanding how brain mechanical properties mature using high-resolution MRE will allow for a deeper understanding of the neural substrates supporting brain function at this age and can inform future studies of atypical maturation., Competing Interests: Conflict of Interest There are no known conflicts of interest associated with this publication and there has been no significant financial support for this work that could have influenced its outcome., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Intrinsic activation MR elastography (IA-MRE) is a novel technique which seeks to estimate brain mechanical properties non-invasively and without external mechanical drivers. The method eliminates actuation hardware and patient discomfort while capitalizing on the brain's intrinsic low frequency motion. This study explores low frequency actuation (1 Hz) MR elastography in phantoms and analyzes performance of non-linear inversion (NLI) of viscoelastic and poroelastic mechanical models as a framework for assessing clinical results from IA-MRE. We present results from four gelatin phantoms and report stiffness resolution of 6 mm (two measurement voxels) with a stiffness contrast ratio of 4.21 relative to the background and 9 mm (three measurement voxels) with a lower stiffness contrast ratio of near 1.77. Stiffness edge resolution was also evaluated using edge spread and line spread functions and yielded a stiffness edge response distance of 9 mm. The intraclass correlation coefficient was high (0.93) between mechanical testing and poroelastic estimates, although quantitative agreement was affected by model-data mismatch. Viscoelastic MRE at low frequencies has issues with non-uniqueness due to small inertial forces, and performed worse than poroelastic MRE in terms of inclusion detection and consistency with mechanical testing. These results present the first evaluation of MR elastography using displacement measurements from an actuation frequency less than 5 Hz and support the validity of brain IA-MRE to recover spatially resolved stiffness changes. They provide a baseline of performance in terms of standard metrics for future animal and human brain stiffness studies and analyses based on intrinsic motion.

Cardiovascular diseases (CVD) are the most prevalent cause of death in the Western World, and their prevalence is only expected to rise. Several screening modalities aim at detecting CVD at the early stages. A common target for early screening is common carotid artery (CCA) stiffness, as reflected in the pulse wave velocity (PWV). For assessing the CCA stiffness using ultrasound (US), one-dimensional (1D) measurements along the CCA axis are typically used, ignoring possible boundary conditions of neck anatomy and the US probe itself. In this study, the effect of stresses and deformations induced by the US probe, and the effect of anatomy surrounding CCA on a simulated 1D stiffness measurement (PWVus) is compared with the ground truth stiffness (PWVgt) in 60 finite-element models (FEM) derived from anatomical computed tomography (CT) scans of ten healthy male volunteers. Based on prior knowledge from the literature, and from results in this study, we conclude that it is safe to approximate arterial stiffness using 1D measurements of compliance or pulse wave velocity, regardless of boundary conditions emerging from the anatomy or from the measurement procedure.

This study evaluated non-linear inversion MRE (NLI-MRE) based on viscoelastic governing equations to determine its sensitivity to small, low contrast inclusions and interface changes in shear storage modulus and damping ratio. Reconstruction parameters identical to those used in recent in vivo MRE studies of mechanical property variations in small brain structures were applied. NLI-MRE was evaluated on four phantoms with contrast in stiffness and damping ratio. Image contrast to noise ratio was assessed as a function of inclusion diameter and property contrast, and edge and line spread functions were calculated as measures of imaging resolution. Phantoms were constructed from silicone, agar, and tofu materials. Reconstructed property estimates were compared with independent mechanical testing using dynamic mechanical analysis (DMA). The NLI-MRE technique detected inclusions as small as 8 mm with a stiffness contrast as low as 14%. Storage modulus images also showed an interface edge response distance of 11 mm. Damping ratio images distinguished inclusions with a diameter as small as 8 mm, and yielded an interface edge response distance of 10 mm. Property differences relative to DMA tests were in the 15%-20% range in most cases. In this study, NLI-MRE storage modulus estimates resolved the smallest inclusion with the lowest stiffness contrast, and spatial resolution of attenuation parameter images was quantified for the first time. These experiments and image quality metrics establish quantitative guidelines for the accuracy expected in vivo for MRE images of small brain structures, and provide a baseline for evaluating future improvements to the NLI-MRE pipeline.

Pulse wave imaging (PWI) is a noninvasive technique for tracking the propagation of the pulse wave along the arterial wall. The 3-D ultrasound imaging would aid in objectively estimating the pulse wave velocity (PWV) vector. This paper aims to introduce a novel PWV estimation method along the propagation direction, validate it in phantoms, and test its feasibility in vivo. A silicone vessel phantom consisting of a stiff and a soft segment along the longitudinal axis and a silicone vessel with a plaque were constructed. A 2-D array with a center frequency of 2.5 MHz was used. Propagation was successfully visualized in 3-D in each phantom and in vivo in six healthy subjects. In three of the healthy subjects, results were compared against conventional PWI using a linear array. PWVs were estimated in the stiff (3.42 ± 0.23 m [Formula: see text]) and soft (2.41 ± 0.07 m [Formula: see text]) phantom segments. Good agreement was found with the corresponding static testing values (stiff: 3.41 m [Formula: see text] and soft: 2.48 m [Formula: see text]). PWI-derived vessel compliance values were validated with dynamic testing. Comprehensive views of pulse propagation in the plaque phantom were generated and compared against conventional PWI acquisitions. Good agreement was found in vivo between the results of 4-D PWI (4.80 ± 1.32 m [Formula: see text]) and conventional PWI (4.28±1.20 m [Formula: see text]) ( n=3 ). PWVs derived for all of the healthy subjects ( n = 6 ) were within the physiological range. Thus, the 4-D PWI was successfully validated in phantoms and used to image the pulse wave propagation in normal human subjects in vivo.

Pulse wave velocity (PWV) is a surrogate marker of arterial stiffness linked to cardiovascular morbidity. Pulse wave imaging (PWI) is a technique developed by our group for imaging the pulse wave propagation in vivo. PWI requires high temporal and spatial resolution, which conventional ultrasonic imaging is unable to simultaneously provide. Coherent compounding is known to address this tradeoff and provides full aperture images at high frame rates. This study aims to implement PWI using coherent compounding within a GPU-accelerated framework. The results of the implemented method were validated using a silicone phantom against static mechanical testing. Reproducibility of the measured PWVs was assessed in the right common carotid of six healthy subjects (n  =  6) approximately 10-15 mm before the bifurcation during two cardiac cycles over the course of 1-3 d. Good agreement of the measured PWVs (3.97  ±  1.21 m s -1 , 4.08  ±  1.15 m s -1 , p  =  0.74) was obtained. The effects of frame rate, transmission angle and number of compounded plane waves on PWI performance were investigated in the six healthy volunteers. Performance metrics such as the reproducibility of the PWVs, the coefficient of determination (r 2 ), the SNR of the PWI axial wall velocities ([Formula: see text]) and the percentage of lateral positions where the pulse wave appears to arrive at the same time-point, indicating inadequacy of the temporal resolution (i.e. temporal resolution misses) were used to evaluate the effect of each parameter. Compounding plane waves transmitted at 1° increments with a linear array yielded optimal performance, generating significantly higher r 2 and [Formula: see text] values (p  ⩽  0.05). Higher frame rates (⩾1667 Hz) produced improvements with significant gains in the r 2 coefficient (p  ⩽  0.05) and significant increase in both r 2 and [Formula: see text] from single plane wave imaging to 3-plane wave compounding (p  ⩽  0.05). Optimal performance was established at 2778 Hz with 3 plane waves and at 1667 Hz with 5 plane waves.

We describe an efficient gradient computation for solving inverse problems arising in magnetic resonance elastography (MRE). The algorithm can be considered as a generalized 'adjoint method' based on a Lagrangian formulation. One requirement for the classic adjoint method is assurance of the self-adjoint property of the stiffness matrix in the elasticity problem. In this paper, we show this property is no longer a necessary condition in our algorithm, but the computational performance can be as efficient as the classic method, which involves only two forward solutions and is independent of the number of parameters to be estimated. The algorithm is developed and implemented in material property reconstructions using poroelastic and viscoelastic modeling. Various gradient- and Hessian-based optimization techniques have been tested on simulation, phantom and in vivo brain data. The numerical results show the feasibility and the efficiency of the proposed scheme for gradient calculation.

Neurological disorders are one of the most important public health concerns in developed countries. Established brain imaging techniques such as magnetic resonance imaging (MRI) and x-ray computerised tomography (CT) have been essential in the identification and diagnosis of a wide range of disorders, although usually are insufficient in sensitivity for detecting subtle pathological alterations to the brain prior to the onset of clinical symptoms-at a time when prognosis for treatment is more favourable. The mechanical properties of biological tissue provide information related to the strength and integrity of the cellular microstructure. In recent years, mechanical properties of the brain have been visualised and measured non-invasively with magnetic resonance elastography (MRE), a particularly sensitive medical imaging technique that may increase the potential for early diagnosis. This review begins with an introduction to the various methods used for the acquisition and analysis of MRE data. A systematic literature search is then conducted to identify studies that have specifically utilised MRE to investigate the human brain. Through the conversion of MRE-derived measurements to shear stiffness (kPa) and, where possible, the loss tangent (rad), a summary of results for global brain tissue and grey and white matter across studies is provided for healthy participants, as potential baseline values to be used in future clinical investigations. In addition, the extent to which MRE has revealed significant alterations to the brain in patients with neurological disorders is assessed and discussed in terms of known pathophysiology. The review concludes by predicting the trends for future MRE research and applications in neuroscience.

The aim of this study was to introduce remote wave excitation for high-resolution cerebral multifrequency MR elastography (mMRE). mMRE of 25-45-Hz drive frequencies by head rocker stimulation was compared with mMRE by remote wave excitation based on a thorax mat in 12 healthy volunteers. Maps of the magnitude |G*| and phase φ of the complex shear modulus were reconstructed using multifrequency dual elasto-visco (MDEV) inversion. After the scan, the subjects and three operators assessed the comfort and convenience of cerebral mMRE using two methods of stimulating the brain. Images were acquired in a coronal view in order to identify anatomical regions along the spinothalamic pathway. In mMRE by remote actuation, all subjects and operators appreciated an increased comfort and simplified procedural set-up. The resulting strain amplitudes in the brain were sufficiently large to analyze using MDEV inversion, and yielded high-resolution viscoelasticity maps which revealed specific anatomical details of brain mechanical properties: |G*| was lowest in the pons (0.97 ± 0.08 kPa) and decreased within the corticospinal tract in the caudal-cranial direction from the crus cerebri (1.64 ± 0.26 kPa) to the capsula interna (1.29 ± 0.14 kPa). By avoiding onerous mechanical stimulation of the head, remote excitation of intracranial shear waves can be used to measure viscoelastic parameters of the brain with high spatial resolution. Therewith, the new mMRE method is suitable for neuroradiological examinations in the clinic., (Copyright © 2015 John Wiley & Sons, Ltd.)

Purpose: Descriptions of the structure of brain tissue as a porous cellular matrix support application of a poroelastic (PE) mechanical model which includes both solid and fluid phases. However, the majority of brain magnetic resonance elastography (MRE) studies use a single phase viscoelastic (VE) model to describe brain tissue behavior, in part due to availability of relatively simple direct inversion strategies for mechanical property estimation. A notable exception is low frequency intrinsic actuation MRE, where PE mechanical properties are imaged with a nonlinear inversion algorithm., Methods: This paper investigates the effect of model choice at each end of the spectrum of in vivo human brain actuation frequencies. Repeat MRE examinations of the brains of healthy volunteers were used to compare image quality and repeatability for each inversion model for both 50 Hz externally produced motion and ≈1 Hz intrinsic motions. Additionally, realistic simulated MRE data were generated with both VE and PE finite element solvers to investigate the effect of inappropriate model choice for ideal VE and PE materials., Results: In vivo, MRE data revealed that VE inversions appear more representative of anatomical structure and quantitatively repeatable for 50 Hz induced motions, whereas PE inversion produces better results at 1 Hz. Reasonable VE approximations of PE materials can be derived by equating the equivalent wave velocities for the two models, provided that the timescale of fluid equilibration is not similar to the period of actuation. An approximation of the equilibration time for human brain reveals that this condition is violated at 1 Hz but not at 50 Hz. Additionally, simulation experiments when using the "wrong" model for the inversion demonstrated reasonable shear modulus reconstructions at 50 Hz, whereas cross-model inversions at 1 Hz were poor quality. Attenuation parameters were sensitive to changes in the forward model at both frequencies, however, no spatial information was recovered because the mechanisms of VE and PE attenuation are different., Conclusions: VE inversions are simpler with fewer unknown properties and may be sufficient to capture the mechanical behavior of PE brain tissue at higher actuation frequencies. However, accurate modeling of the fluid phase is required to produce useful mechanical property images at the lower frequencies of intrinsic brain motions.

Purpose: To develop an acquisition scheme for generating MR elastography (MRE) displacement data with whole-brain coverage, high spatial resolution, and adequate signal-to-noise ratio (SNR) in a short scan time., Theory and Methods: A 3D multislab, multishot acquisition for whole-brain MRE with 2.0 mm isotropic spatial resolution is proposed. The multislab approach allowed for the use of short repetition time to achieve very high SNR efficiency. High SNR efficiency allowed for a reduced acquisition time of only 6 min while the minimum SNR needed for inversion was maintained., Results: The mechanical property maps estimated from whole-brain displacement data with nonlinear inversion (NLI) demonstrated excellent agreement with neuroanatomical features, including the cerebellum and brainstem. A comparison with an equivalent 2D acquisition illustrated the improvement in SNR efficiency of the 3D multislab acquisition. The flexibility afforded by the high SNR efficiency allowed for higher resolution with a 1.6 mm isotropic voxel size, which generated higher estimates of brainstem stiffness compared with the 2.0 mm isotropic acquisition., Conclusion: The acquisition presented allows for the capture of whole-brain MRE displacement data in a short scan time, and may be used to generate local mechanical property estimates of neuroanatomical features throughout the brain., (Copyright © 2013 Wiley Periodicals, Inc.)

The noninvasive measurement of the mechanical properties of brain tissue using magnetic resonance elastography (MRE) has emerged as a promising method for investigating neurological disorders. To date, brain MRE investigations have been limited to reporting global mechanical properties, though quantification of the stiffness of specific structures in the white matter architecture may be valuable in assessing the localized effects of disease. This paper reports the mechanical properties of the corpus callosum and corona radiata measured in healthy volunteers using MRE and atlas-based segmentation. Both structures were found to be significantly stiffer than overall white matter, with the corpus callosum exhibiting greater stiffness and less viscous damping than the corona radiata. Reliability of both local and global measures was assessed through repeated experiments, and the coefficient of variation for each measure was less than 10%. Mechanical properties within the corpus callosum and corona radiata demonstrated correlations with measures from diffusion tensor imaging pertaining to axonal microstructure., (Copyright © 2013 Elsevier Inc. All rights reserved.)

Magnetic resonance elastography (MRE) has been introduced in clinical practice as a possible surrogate for mechanical palpation, but its application to study the human brain in vivo has been limited by low spatial resolution and the complexity of the inverse problem associated with biomechanical property estimation. Here, we report significant improvements in brain MRE data acquisition by reporting images with high spatial resolution and signal-to-noise ratio as quantified by octahedral shear strain metrics. Specifically, we have developed a sequence for brain MRE based on multishot, variable-density spiral imaging, and three-dimensional displacement acquisition and implemented a correction scheme for any resulting phase errors. A Rayleigh damped model of brain tissue mechanics was adopted to represent the parenchyma and was integrated via a finite element-based iterative inversion algorithm. A multiresolution phantom study demonstrates the need for obtaining high-resolution MRE data when estimating focal mechanical properties. Measurements on three healthy volunteers demonstrate satisfactory resolution of gray and white matter, and mechanical heterogeneities correspond well with white matter histoarchitecture. Together, these advances enable MRE scans that result in high-fidelity, spatially resolved estimates of in vivo brain tissue mechanical properties, improving upon lower resolution MRE brain studies that only report volume averaged stiffness values., (© 2012 Wiley Periodicals, Inc.)

Purpose: Breast cancer is a major public health issue for women, and early detection significantly increases survival rate. Currently, there is increased research interest in elastographic soft-tissue imaging techniques based on the correlation between pathology and mechanical stiffness. Anthropomorphic breast phantoms are critical for ex vivo validation of emerging elastographic technologies. This research develops heterogeneous breast phantoms for use in testing elastographic imaging modalities., Methods: Mechanical property estimation of eight different elastomers is performed to determine storage moduli (E') and damping ratios (ζ) using a dynamic mechanical analyzer. Dynamic compression testing was carried out isothermally at room temperature over a range of 4-50 Hz. Silicone compositions with physiologically realistic storage modulus were chosen for mimicking skin adipose, cancerous tumors, and pectoral muscles and 13 anthropomorphic breast phantoms were constructed for ex vivo trials of digital image elastotomography (DIET) breast cancer screening system. A simpler fabrication was used to assess the possibility of multiple tumor detection using magnetic resonance elastography (MRE)., Results: Silicone materials with ranges of storage moduli (E') from 2 to 570 kPa and damping ratios (ζ) from 0.03 to 0.56 were identified. The resulting phantoms were tested in two different elastographic breast cancer diagnostic modalities. A significant contrast was successfully identified between healthy tissues and cancerous tumors both in MRE and DIET., Conclusions: The phantoms presented promise aid to researchers in elastographic imaging modalities for breast cancer detection and provide a foundation for silicone based phantom materials for mimicking soft tissues of other human organs.

Many pathologies alter the mechanical properties of tissue. Magnetic resonance elastography (MRE) has been developed to noninvasively characterize these quantities in vivo. Typically, small vibrations are induced in the tissue of interest with an external mechanical actuator. The resulting displacements are measured with phase contrast sequences and are then used to estimate the underlying mechanical property distribution. Several MRE studies have quantified brain tissue properties. However, the cranium and meninges, especially the dura, are very effective at damping externally applied vibrations from penetrating deeply into the brain. Here, we report a method, termed 'intrinsic activation', that eliminates the requirement for external vibrations by measuring the motion generated by natural blood vessel pulsation. A retrospectively gated phase contrast MR angiography sequence was used to record the tissue velocity at eight phases of the cardiac cycle. The velocities were numerically integrated via the Fourier transform to produce the harmonic displacements at each position within the brain. The displacements were then reconstructed into images of the shear modulus based on both linear elastic and poroelastic models. The mechanical properties produced fall within the range of brain tissue estimates reported in the literature and, equally important, the technique yielded highly reproducible results. The mean shear modulus was 8.1 kPa for linear elastic reconstructions and 2.4 kPa for poroelastic reconstructions where fluid pressure carries a portion of the stress. Gross structures of the brain were visualized, particularly in the poroelastic reconstructions. Intra-subject variability was significantly less than the inter-subject variability in a study of six asymptomatic individuals. Further, larger changes in mechanical properties were observed in individuals when examined over time than when the MRE procedures were repeated on the same day. Cardiac pulsation, termed intrinsic activation, produces sufficient motion to allow mechanical properties to be recovered. The poroelastic model is more consistent with the measured data from brain at low frequencies than the linear elastic model. Intrinsic activation allows MRE to be performed without a device shaking the head so the patient notices no differences between it and the other sequences in an MR examination.

Purpose: Nonlinear inversion (NLI) in MR elastography requires discretization of the displacement field for a finite element (FE) solution of the "forward problem", and discretization of the unknown mechanical property field for the iterative solution of the "inverse problem". The resolution requirements for these two discretizations are different: the forward problem requires sufficient resolution of the displacement FE mesh to ensure convergence, whereas lowering the mechanical property resolution in the inverse problem stabilizes the mechanical property estimates in the presence of measurement noise. Previous NLI implementations use the same FE mesh to support the displacement and property fields, requiring a trade-off between the competing resolution requirements., Methods: This work implements and evaluates multiresolution FE meshes for NLI elastography, allowing independent discretizations of the displacements and each mechanical property parameter to be estimated. The displacement resolution can then be selected to ensure mesh convergence, and the resolution of the property meshes can be independently manipulated to control the stability of the inversion., Results: Phantom experiments indicate that eight nodes per wavelength (NPW) are sufficient for accurate mechanical property recovery, whereas mechanical property estimation from 50 Hz in vivo brain data stabilizes once the displacement resolution reaches 1.7 mm (approximately 19 NPW). Viscoelastic mechanical property estimates of in vivo brain tissue show that subsampling the loss modulus while holding the storage modulus resolution constant does not substantially alter the storage modulus images. Controlling the ratio of the number of measurements to unknown mechanical properties by subsampling the mechanical property distributions (relative to the data resolution) improves the repeatability of the property estimates, at a cost of modestly decreased spatial resolution., Conclusions: Multiresolution NLI elastography provides a more flexible framework for mechanical property estimation compared to previous single mesh implementations.

A signal-to-noise ratio (SNR) measure based on the octahedral shear strain (the maximum shear strain in any plane for a 3D state of strain) is presented for magnetic resonance elastography (MRE), where motion-based SNR measures are commonly used. The shear strain, γ, is directly related to the shear modulus, μ, through the definition of shear stress, τ = μγ. Therefore, noise in the strain is the important factor in determining the quality of motion data, rather than the noise in the motion. Motion and strain SNR measures were found to be correlated for MRE of gelatin phantoms and the human breast. Analysis of the stiffness distributions of phantoms reconstructed from the measured motion data revealed a threshold for both strain and motion SNR where MRE stiffness estimates match independent mechanical testing. MRE of the feline brain showed significantly less correlation between the two SNR measures. The strain SNR measure had a threshold above which the reconstructed stiffness values were consistent between cases, whereas the motion SNR measure did not provide a useful threshold, primarily due to rigid body motion effects.

Purpose: Recently, the attenuating behavior of soft tissue has been addressed in magnetic resonance elastography by the inclusion of a damping mechanism in the methods used to reconstruct the resulting mechanical property image. To date, this mechanism has been based on a viscoelastic model for material behavior. Rayleigh, or proportional, damping provides a more generalized model for elastic energy attenuation that uses two parameters to characterize contributions proportional to elastic and inertial forces. In the case of time-harmonic vibration, these two parameters lead to both the elastic modulus and the density being complex valued (as opposed to the case of pure viscoelasticity, where only the elastic modulus is complex valued)., Methods: This article presents a description of Rayleigh damping in the time-harmonic case, discussing the differences between this model and the viscoelastic damping models. In addition, the results from a subzone based Rayleigh damped elastography study of gelatin and tofu phantoms are discussed, along with preliminary results from in vivo breast data., Results: Both the phantom and the tissue studies presented here indicate a change in the Rayleigh damping structure, described as Rayleigh composition, between different material types, with tofu and healthy tissue showing lower Rayleigh composition values than gelatin or cancerous tissue., Conclusions: It is possible that Rayleigh damping elastography and the concomitant Rayleigh composition images provide a mechanism for differentiating tissue structure in addition to measuring elastic stiffness and attenuation. Such information could be valuable in the use of Rayleigh damped magnetic resonance elastography as a diagnostic imaging tool.

The mechanical model commonly used in magnetic resonance elastography (MRE) is linear elasticity. However, soft tissue may exhibit frequency- and direction-dependent (FDD) shear moduli in response to an induced excitation causing a purely linear elastic model to provide an inaccurate image reconstruction of its mechanical properties. The goal of this study was to characterize the effects of reconstructing FDD data using a linear elastic inversion (LEI) algorithm. Linear and FDD phantoms were manufactured and LEI images were obtained from time-harmonic MRE acquisitions with variations in frequency and driving signal amplitude. LEI responses to artificially imposed uniform phase shifts in the displacement data from both purely linear elastic and FDD phantoms were also evaluated. Of the variety of FDD phantoms considered, LEI appeared to tolerate viscoelastic data-model mismatch better than deviations caused by poroelastic and anisotropic mechanical properties in terms of visual image contrast. However, the estimated shear modulus values were substantially incorrect relative to independent mechanical measurements even in the successful viscoelastic cases and the variations in mean values with changes in experimental conditions associated with uniform phase shifts, driving signal frequency and amplitude were unpredictable. Overall, use of LEI to reconstruct data acquired in phantoms with FDD material properties provided biased results under the best conditions and significant artifacts in the worst cases. These findings suggest that the success with which LEI is applied to MRE data in tissue will depend on the underlying mechanical characteristics of the tissues and/or organs systems of clinical interest.

Imaging of the mechanical properties of in vivo brain tissue could eventually lead to non-invasive diagnosis of hydrocephalus, Alzheimer's disease and other pathologies known to alter the intracranial environment. The purpose of this work is to (1) use time-harmonic magnetic resonance elastography (MRE) to estimate the mechanical property distribution of cerebral tissue in the normal feline brain and (2) compare the recovered properties of grey and white matter. Various in vivo and ex vivo brain tissue property measurement strategies have led to the highly variable results that have been reported in the literature. MR elastography is an imaging technique that can estimate mechanical properties of tissue non-invasively and in vivo. Data was acquired in 14 felines and elastic parameters were estimated using a globo-regional nonlinear image reconstruction algorithm. Results fell within the range of values reported in the literature and showed a mean shear modulus across the subject group of 7-8 kPa with all but one animal falling within 5-15 kPa. White matter was statistically stiffer (p<0.01) than grey matter by about 1 kPa on a per subject basis. To the best of our knowledge, the results reported represent the most extensive set of estimates in the in vivo brain which have been based on MRE acquisition of the three-dimensional displacement field coupled to volumetric shear modulus image reconstruction achieved through nonlinear parameter estimation. However, the inter-subject variation in mean shear modulus indicates the need for further study, including the possibility of applying more advanced models to estimate the relevant tissue mechanical properties from the data., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

A Rayleigh damping model applied to magnetic resonance elastography incorporates attenuation behavior proportionally related to both elastic and inertial forces, and allows two damping parameters to be extracted from an MRI motion dataset. Under time-harmonic conditions, the model can be implemented by the use of complex shear modulus and density, whereas viscoelastic damping models commonly used in elastography consist of only a complex shear modulus, and model only a single damping effect. Simulation studies reveal that the differences between damped elastic behavior resulting from a purely complex shear modulus (CSM damping) and from a purely complex density (CD damping) become larger as the overall level of damping present (indicated by the damping ratio) increases. A plot of results generated from the finite element (FE) model indicate the relative motion differences estimated for a range of damping ratios and CSM/CD damping combinations increase with damping ratio, and can be up to 15% at a damping ratio of 50% and therefore using the correct model for a Rayleigh damped material becomes increasingly important as damping levels increase. Resonance-related effects cause values from this plot to vary by as much as 3% as parameters such as wave speed, frequency, and problem size are altered. These motion differences can be compared to expected noise levels to estimate the parameter resolution achievable by a reconstruction algorithm. An optimization-based global property reconstruction algorithm was developed, and used for testing Rayleigh damping parameter reconstructions with gaussian noise added to the simulated motion input data. The coherent motion errors resulting from altering the combination of the two damping parameters are large enough to allow accurate determination of both of the Rayleigh damping parameters with incoherent noise levels comparable to MR measurements. The accuracy achieved by the global reconstructions was significantly better than would be predicted by examining the motion differences for differing CSM/CD damping combinations, which is likely to be due to the low ratio between number of reconstructed parameters and number of noisy measurements.