Measurement of strain in physical models of brain injury: a method based on HARP analysis of tagged magnetic resonance images (MRI)

Two-dimensional (2-D) strain fields were estimated non-invasively in two simple experimental models of closed-head brain injury. In the first experimental model, shear deformation of a gel was induced by angular acceleration of its spherical container In the second model the brain of a euthanized rat pup was deformed by indentation of its skull. Tagged magnetic resonance images (MRI) were obtained by gated image acquisition during repeated motion. Harmonic phase (HARP) images corresponding to the spectral peaks of the original tagged MRI were obtained, following procedures proposed by Osman, McVeigh and Prince. Two methods of HARP strain analysis were applied, one based on the displacement of tag line intersections, and the other based on the gradient of harmonic phase. Strain analysis procedures were also validated on simulated images of deformed grids. Results show that it is possible to visualize deformation and to quantify strain efficiently in animal models of closed head injury.     

Deformation of the human brain induced by mild angular head acceleration

Deformation of the human brain was measured in tagged magnetic resonance images (MRI) obtained dynamically during angular acceleration of the head. This study was undertaken to provide quantitative experimental data to illuminate the mechanics of traumatic brain injury (TBI). Mild angular acceleration was imparted to the skull of a human volunteer inside an MR scanner, using a custom MR-compatible device to constrain motion. A grid of MR "tag" lines was applied to the MR images via spatial modulation of magnetization (SPAMM) in a fast gradient echo imaging sequence. Images of the moving brain were obtained dynamically by synchronizing the imaging process with the motion of the head. Deformation of the brain was characterized quantitatively via Lagrangian strain. Consistent patterns of radial-circumferential shear strain occur in the brain, similar to those observed in models of a viscoelastic gel cylinder subjected to angular acceleration. Strain fields in the brain, however, are clearly mediated by the effects of heterogeneity, divisions between regions of the brain (such as the central fissure and central sulcus) and the brain's tethering and suspension system, including the dura mater, falx cerebri, and tentorium membranes.     

Imaging two-dimensional displacements and strains in skeletal muscle during joint motion by cine DENSE MR

The objective of this study was to apply cine magnetic resonance imaging (MRI) using displacement encoding with stimulated echoes (DENSE) to measure the dynamic two-dimensional (2D) displacement and Lagrangian strain fields in the biceps brachii muscle. Six healthy volunteers underwent cine DENSE MRI during repeated elbow flexion against the load of gravity. Displacement encoded dynamic images of the upper arm were acquired with spatial and temporal resolutions of 1.9 x 1.9 mm(2) and 30 ms, respectively. Pixel-wise Lagrangian displacement and strain fields were calculated from the measured images. We extracted the first and second principal strains (E1 and E2) along the centerline and anterior regions of the muscle. E1 and E2 were relatively uniform along the anterior region. However, E1 and E2 were both non-uniform along the centerline region-normalized values for E1 and E2 varied over the ranges of 0.27-1.35, and 0.45-2.36, respectively. The directions of the first and second principal strains varied throughout the muscle and showed that the direction of principal shortening is not necessarily aligned with fascicle direction. This study demonstrates the utility of cine DENSE MRI for analyzing skeletal muscle mechanics and provides data describing the in vivo mechanics of muscle tissue to a level of detail that has not been previously possible.     

Measurement of the dynamic shear modulus of mouse brain tissue in vivo by magnetic resonance elastography

In this study, the magnetic resonance (MR) elastography technique was used to estimate the dynamic shear modulus of mouse brain tissue in vivo. The technique allows visualization and measurement of mechanical shear waves excited by lateral vibration of the skull. Quantitative measurements of displacement in three dimensions during vibration at 1200 Hz were obtained by applying oscillatory magnetic field gradients at the same frequency during a MR imaging sequence. Contrast in the resulting phase images of the mouse brain is proportional to displacement. To obtain estimates of shear modulus, measured displacement fields were fitted to the shear wave equation. Validation of the procedure was performed on gel characterized by independent rheometry tests and on data from finite element simulations. Brain tissue is, in reality, viscoelastic and nonlinear. The current estimates of dynamic shear modulus are strictly relevant only to small oscillations at a specific frequency, but these estimates may be obtained at high frequencies (and thus high deformation rates), noninvasively throughout the brain. These data complement measurements of nonlinear viscoelastic properties obtained by others at slower rates, either ex vivo or invasively.     

Measurement of in vivo anterior cruciate ligament strain during dynamic jump landing

Despite recent attention in the literature, anterior cruciate ligament (ACL) injury mechanisms are controversial and incidence rates remain high. One explanation is limited data on in vivo ACL strain during high-risk, dynamic movements. The objective of this study was to quantify ACL strain during jump landing. Marker-based motion analysis techniques were integrated with fluoroscopic and magnetic resonance (MR) imaging techniques to measure dynamic ACL strain non-invasively. First, eight subjects' knees were imaged using MR. From these images, the cortical bone and ACL attachment sites of the tibia and femur were outlined to create 3D models. Subjects underwent motion analysis while jump landing using reflective markers placed directly on the skin around the knee. Next, biplanar fluoroscopic images were taken with the markers in place so that the relative positions of each marker to the underlying bone could be quantified. Numerical optimization allowed jumping kinematics to be superimposed on the knee model, thus reproducing the dynamic in vivo joint motion. ACL length, knee flexion, and ground reaction force were measured. During jump landing, average ACL strain peaked 55±14 ms (mean and 95% confidence interval) prior to ground impact, when knee flexion angles were lowest. The peak ACL strain, measured relative to its length during MR imaging, was 12±7%. The observed trends were consistent with previously described neuromuscular patterns. Unrestricted by field of view or low sampling rate, this novel approach provides a means to measure kinematic patterns that elevate ACL strains and that provide new insights into ACL injury mechanisms.     

Optimization of Tagged MRI for Quantification of Liver Stiffness Using Computer Simulated Data

The heartbeat has been proposed as an intrinsic source of motion that can be used in combination with tagged Magnetic Resonance Imaging (MRI) to measure displacements induced in the liver as an index of liver stiffness. Optimizing a tagged MRI acquisition protocol in terms of sensitivity to these displacements, which are in the order of pixel size, is necessary to develop the method as a quantification tool for staging fibrosis. We reproduced a study of cardiac-induced strain in the liver at 3T and simulated tagged MR images with different grid tag patterns to evaluate the performance of the Harmonic Phase (HARP) image analysis method and its dependence on the parameters of tag spacing and grid angle. The Partial Volume Effect (PVE), T1 relaxation, and different levels of noise were taken into account. Four displacement fields of increasing intensity were created and applied to the tagged MR images of the liver. These fields simulated the deformation at different liver stiffnesses. An Error Index (EI) was calculated to evaluate the estimation accuracy for various parameter values. In the absence of noise, the estimation accuracy of the displacement fields increased as tag spacings decreased. EIs for each of the four displacement fields were lower at 0° and the local minima of the EI were found to correspond to multiples of pixel size. The accuracy of the estimation decreased for increasing levels of added noise; as the level increased, the improved estimation caused by decreasing the tag spacing tended to zero. The optimal tag spacing turned out to be a compromise between the smallest tag period that is a multiple of the pixel size and is achievable in a real acquisition and the tag spacing that guarantees an accurate liver displacement measure in the presence of realistic levels of noise.     

In vivo tibiofemoral cartilage deformation during the stance phase of gait

The knowledge of articular cartilage contact biomechanics in the knee joint is important for understanding the joint function and cartilage pathology. However, the in vivo tibiofemoral articular cartilage contact biomechanics during gait remains unknown. The objective of this study was to determine the in vivo tibiofemoral cartilage contact biomechanics during the stance phase of treadmill gait. Eight healthy knees were magnetic resonance (MR) scanned and imaged with a dual fluoroscopic system during gait on a treadmill. The tibia, femur and associated cartilage were constructed from the MR images and combined with the dual fluoroscopic images to determine in vivo cartilage contact deformation during the stance phase of gait. Throughout the stance phase of gait, the magnitude of peak compartmental contact deformation ranged between 7% and 23% of the resting cartilage thickness and occurred at regions with thicker cartilage. Its excursions in the anteroposterior direction were greater in the medial tibiofemoral compartment as compared to those in the lateral compartment. The contact areas throughout the stance phase were greater in the medial compartment than in the lateral compartment. The information on in vivo tibiofemoral cartilage contact biomechanics during gait could be used to provide physiological boundaries for in vitro testing of cartilage. Also, the data on location and magnitude of deformation among non-diseased knees during gait could identify where loading and later injury might occur in diseased knees.     

In vivo diffusion-perfusion magnetic resonance imaging of acute cerebral ischemia.

We compared the anatomic extent and severity of ischemic brain injury shown on diffusion-weighted magnetic resonance (MR) images, with cerebral tissue perfusion deficits demonstrated by a nonionic intravascular T2*-shortening magnetic susceptibility contrast agent used in conjunction with standard T2-weighted spin-echo and gradient-echo echo-planar images. Diffusion-weighted images displayed increased signal intensity in the vascular territory of the middle cerebral artery 25-40 min after permanent occlusion, whereas T2-weighted images without contrast were negative or equivocal for at least 2-3 h after stroke was induced. Contrast-enhanced T2-weighted and echo-planar images revealed perfusion deficits that were spatially closely related to the anatomic regions of ischemic tissue injury. These data indicate that diffusion-weighted MR images are very sensitive to early onset pathophysiologic changes induced by acute cerebral ischemia. Combined sequential diffusion-perfusion imaging enables noninvasive in vivo examination of the relationship between hypoperfusion and evolving ischemic brain injury.     

Contrast-Enhanced FLAIR (Fluid-Attenuated Inversion Recovery) for Evaluating Mild Traumatic Brain Injury

Purpose

To evaluate whether adding a contrast-enhanced fluid-attenuated inversion recovery (FLAIR) sequence to routine magnetic resonance imaging (MRI) can detect additional abnormalities in the brains of symptomatic patients with mild traumatic brain injury.

Materials and Methods

Fifty-four patients with persistent symptoms following mild closed head injury were included in our retrospective study (M∶F = 32∶22, mean age: 59.8±16.4, age range: 26–84 years). All MRI examinations were obtained within 14 days after head trauma (mean: 3.2±4.1 days, range: 0.2–14 days). Two neuroradiologists recorded (1) the presence of traumatic brain lesions on MR images with and without contrast-enhanced FLAIR images and (2) the pattern and location of meningeal enhancement depicted on contrast-enhanced FLAIR images. The number of additional traumatic brain lesions diagnosed with contrast-enhanced FLAIR was recorded. Correlations between meningeal enhancement and clinical findings were also evaluated.

Results

Traumatic brain lesions were detected on routine image sequences in 25 patients. Three additional cases of brain abnormality were detected with the contrast-enhanced FLAIR images. Meningeal enhancement was identified on contrast-enhanced FLAIR images in 9 cases while the other routine image sequences showed no findings of traumatic brain injury. Overall, the additional contrast-enhanced FLAIR images revealed more extensive abnormalities than routine imaging in 37 cases (p<0.001). In multivariate logistic regression analysis, subdural hematoma and posttraumatic loss of consciousness showed a significant association with meningeal enhancement on contrast-enhanced FLAIR images, with odds ratios 13.068 (95% confidence interval 2.037 to 83.852), and 15.487 (95% confidence interval 2.545 to 94.228), respectively.

Conclusion

Meningeal enhancement on contrast-enhanced FLAIR images can help detect traumatic brain lesions as well as additional abnormalities not identified on routine unenhanced MRI. Therefore contrast-enhanced FLAIR MR imaging is recommended when a contrast MR study is indicated in a patient with a symptomatic prior closed mild head injury.     

The effect of the elongation of the proximal aorta on the estimation of the aortic wall distensibility

The compliance of the proximal aortic wall is a major determinant of cardiac afterload. Aortic compliance is often estimated based on cross-sectional area changes over the pulse pressure, under the assumption of a negligible longitudinal stretch during the pulse. However, the proximal aorta is subjected to significant axial stretch during cardiac contraction. In the present study, we sought to evaluate the importance of axial stretch on compliance estimation by undertaking both an in silico and an in vivo approach. In the computational analysis, we developed a 3-D finite element model of the proximal aorta and investigated the discrepancy between the actual wall compliance to the value estimated after neglecting the longitudinal stretch of the aorta. A parameter sensitivity analysis was further conducted to show how increased material stiffness and increased aortic root motion might amplify the estimation errors (discrepancies between actual and estimated distensibility ranging from − 20 to − 62%). Axial and circumferential aortic deformation during ventricular contraction was also evaluated in vivo based on MR images of the aorta of 3 healthy young volunteers. The in vivo results were in good qualitative agreement with the computational analysis (underestimation errors ranging from − 26 to − 44%, with increased errors reflecting higher aortic root displacement). Both the in silico and in vivo findings suggest that neglecting the longitudinal strain during contraction might lead to severe underestimation of local aortic compliance, particularly in the case of women who tend to have higher aortic root motion or in subjects with stiff aortas.

     

Noninvasive high resolution mechanical strain maps of the spine intervertebral disc using nonrigid registration of magnetic resonance images

High resolution strain measurements are of particular interest in load bearing tissues such as the intervertebral disc (IVD), permitting characterization of biomechanical conditions which could lead to injury and degenerative outcomes. Magnetic resonance (MR) imaging produces excellent image contrast in cartilaginous tissues, allowing for image-based strain determination. Nonrigid registration (NRR) of MR images has previously demonstrated sub-voxel registration accuracy although its accuracy and precision in determining strain has not been evaluated. Accuracy and precision of NRR-derived strain measurements were evaluated using computer generated deformations applied to both computer generated images and MR images. Two different measures of registration similarity--the cost function which drives the registration algorithm--were compared: Mutual Information (MI) and Least Squares (LS). Strain error was evaluated with respect to signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and strain heterogeneity. Additionally, the creep strain response from an in vitro loaded porcine IVD is shown and comparisons between similarity measures are presented. MI showed a decrease in strain precision with increasing CNR and decreasing SNR while LS was insensitive to both. Both similarity measures showed a decrease in strain precision with increasing strain heterogeneity. When computer generated heterogeneous strains were applied to MR images of the IVD, LS showed substantially lower strain error in comparison to MI. Results suggest that LS-driven NRR provides a more accurate image-based method for mapping large and heterogeneous strain fields and this method can be applied to studies of the IVD and, potentially, other soft tissues which present sufficient image texture.     

Patient-specific non-linear finite element modelling for predicting soft organ deformation in real-time; Application to non-rigid neuroimage registration

Long computation times of non-linear (i.e. accounting for geometric and material non-linearity) biomechanical models have been regarded as one of the key factors preventing application of such models in predicting organ deformation for image-guided surgery. This contribution presents real-time patient-specific computation of the deformation field within the brain for six cases of brain shift induced by craniotomy (i.e. surgical opening of the skull) using specialised non-linear finite element procedures implemented on a graphics processing unit (GPU). In contrast to commercial finite element codes that rely on an updated Lagrangian formulation and implicit integration in time domain for steady state solutions, our procedures utilise the total Lagrangian formulation with explicit time stepping and dynamic relaxation. We used patient-specific finite element meshes consisting of hexahedral and non-locking tetrahedral elements, together with realistic material properties for the brain tissue and appropriate contact conditions at the boundaries. The loading was defined by prescribing deformations on the brain surface under the craniotomy. Application of the computed deformation fields to register (i.e. align) the preoperative and intraoperative images indicated that the models very accurately predict the intraoperative deformations within the brain. For each case, computing the brain deformation field took less than 4 s using an NVIDIA Tesla C870 GPU, which is two orders of magnitude reduction in computation time in comparison to our previous study in which the brain deformation was predicted using a commercial finite element solver executed on a personal computer.     

Heparin affin regulatory peptide in milk: its involvement in mammary gland homeostasis

HARP (heparin affin regulatory peptide) is a heparin binding growth factor implicated in cellular growth and differentiation. Previously, HARP had been localized in the human mammary, in both alveolar epithelial and myoepithelial cells although HARP mRNAs were only expressed by myoepithelial cells [J. Histochem. Cytochem. 45 (1997) 1]. In the present study, we demonstrate that HARP is secreted in human mature milk with concentrations ranging from 17.68+/-6.4ng/ml in mature milk to 59.9+/-11.22ng/ml in colostrum. In vitro, HARP was found to be mitogenic on human mammary epithelial and myoepithelial cell lines and correlated with the expression of its high affinity receptor tyrosine kinase ALK (anaplastic lymphoma kinase). In vivo, ALK is expressed in both mammary epithelial and myoepithelial cells, suggesting that HARP could act in vivo as a paracrine and autocrine growth factor in the regulation of the mammary gland development and its homeostatic maintenance during pregnancy and lactation.     

Improved measurement of brain deformation during mild head acceleration using a novel tagged MRI sequence

In vivo measurements of human brain deformation during mild acceleration are needed to help validate computational models of traumatic brain injury and to understand the factors that govern the mechanical response of the brain. Tagged magnetic resonance imaging is a powerful, noninvasive technique to track tissue motion in vivo which has been used to quantify brain deformation in live human subjects. However, these prior studies required from 72 to 144 head rotations to generate deformation data for a single image slice, precluding its use to investigate the entire brain in a single subject. Here, a novel method is introduced that significantly reduces temporal variability in the acquisition and improves the accuracy of displacement estimates. Optimization of the acquisition parameters in a gelatin phantom and three human subjects leads to a reduction in the number of rotations from 72 to 144 to as few as 8 for a single image slice. The ability to estimate accurate, well-resolved, fields of displacement and strain in far fewer repetitions will enable comprehensive studies of acceleration-induced deformation throughout the human brain in vivo.     

Inertia effects on characterization of dynamic response of brain tissue

Modeling and simulation of traumatic brain injury (TBI) resulted from collision or blast loading requires characterization of mechanical response over a wide range of loading rates under valid testing conditions. In this study, mechanical response of fresh bovine brain tissue was studied using the two modified Kolsky bar techniques. Radial deformation behavior of annular specimens, which are typically used to characterize the dynamic uniaxial compressive response of biological tissues, was examined using a modified Kolsky bar and a high speed camera to collect images while the specimen deforms at an axial strain rate of 2000s(-1). The high-speed images revealed inhomogeneous specimen deformation possibly brought about by radial inertia and causing a multi-axial stress state. To acquire valid stress-strain results that can be used to produce constitutive behavior of the soft materials, a novel torsion technique was developed to obtain pure shear response at dynamic loading rates. Experimental results show clear differences in the material response using the two methods. These results indicate that the previously demonstrated annular specimen geometry aimed at reducing inertia induced stress components for high rate soft materials uniaxial-compressive testing may still possess a significant component of radial inertia induced radial stress which consequently caused the observed inhomogeneous deformation in brain tissue test samples.     

磁共振肿瘤矩特征信息研究

目的:采用MR脑肿瘤图像分割与矩方法进行结合,以获取特定器官及组织的轮廓。方法:对MR脑肿瘤图像进行分割,并对分割的结果进行矩描述。通过分析当前常用的医学图像分割方法,采用了一种基于形变模型的医学图像分割方法,并按照相应的理论算法模型和实现步骤对医学图像进行了处理,最后用Visual C 6.0编程,并对MR脑肿瘤图像进行分割实验。结果:从切割的图形中可以看出,本分割方法分割边界清晰,总体不确定性较小,利用矩技术所提取的图像特征在基于内容的图像检索中是有效的。结论:本分割方法切实可行,分割效果较好,为进一步的MR脑肿瘤图像分析和研究提供了一种有效工具。     

Tissue-level thresholds for axonal damage in an experimental model of central nervous system white matter injury

In vivo, tissue-level, mechanical thresholds for axonal injury were determined by comparing morphological injury and electrophysiological impairment to estimated tissue strain in an in vivo model of axonal injury. Axonal injury was produced by dynamically stretching the right optic nerve of an adult male guinea pig to one of seven levels of ocular displacement (Nlevel = 10; Ntotal = 70). Morphological injury was detected with neurofilament immunohistochemical staining (NF68, SM132). Simultaneously, functional injury was determined by the magnitude of the latency shift of the N35 peak of the visual evoked potentials (VEPs) recorded before and after stretch. A companion set of in situ experiments (Nlevel = 5) was used to determine the empirical relationship between the applied ocular displacement and the magnitude of optic nerve stretch. Logistic regression analysis, combined with sensitivity and specificity measures and receiver operating characteristic (ROC) curves were used to predict strain thresholds for axonal injury. From this analysis, we determined three Lagrangian strain-based thresholds for morphological damage to white matter. The liberal threshold, intended to minimize the detection of false positives, was a strain of 0.34, and the conservative threshold strain that minimized the false negative rate was 0.14. The optimal threshold strain criterion that balanced the specificity and sensitivity measures was 0.21. Similar comparisons for electrophysiological impairment produced liberal, conservative, and optimal strain thresholds of 0.28, 0.13, and 0.18, respectively. With these threshold data, it is now possible to predict more accurately the conditions that cause axonal injury in human white matter.     

In vivo biomechanical behavior of the human heel pad during the stance phase of gait

A technique is introduced for simultaneous measurements of the heel pad tissue deformation and the heel–ground contact stresses developing during the stance phase of gait. Subjects walked upon a gait platform integrating the contact pressure display optical method for plantar pressure measurements and a digital radiographic fluoroscopy system for skeletal and soft tissue motion recording. Clear images of the posterior-plantar aspect of the calcaneus and enveloping soft tissues were obtained simultaneously with the pressure distribution under the heel region throughout the stance phase of gait. The heel pad was shown to undergo a rapid compression during initial contact and heel strike, reaching a strain of 0.39±0.05 in about 150 ms. The stress–strain relation of the heel pad was shown to be highly non-linear, with a compression modulus of 105±11 kPa initially and 306±16 kPa at 30% strain. The energy dissipation during heel strike was evaluated to be 17.8±0.8%. The present technique is useful for biomechanical as well as clinical evaluation of the stress–strain and energy absorption characteristics of the heel pad in vivo, during natural gait.     

The influence of deceleration forces on ACL strain during single-leg landing: a simulation study

Anterior cruciate ligament (ACL) injury commonly occurs during single limb landing or stopping from a run, yet the conditions that influence ACL strain are not well understood. The purpose of this study was to develop, test and apply a 3D specimen-specific dynamic simulation model of the knee designed to evaluate the influence of deceleration forces during running to a stop (single-leg landing) on ACL strain. This work tested the conceptual development of the model by simulating a physical experiment that provided direct measurements of ACL strain during vertical impact loading (peak value 1294N) with the leg near full extension. The properties of the soft tissue structures were estimated by simulating previous experiments described in the literature. A key element of the model was obtaining precise anatomy from segmented MR images of the soft tissue structures and articular geometry for the tibiofemoral and patellofemoral joints of the knee used in the cadaver experiment. The model predictions were correlated (Pearson correlation coefficient 0.889) to the temporal and amplitude characteristic of the experimental strains. The simulation model was then used to test the balance between ACL strain produced by quadriceps contraction and the reductions in ACL strain associated with the posterior braking force. When posterior forces that replicated in vivo conditions were applied, the peak ACL strain was reduced. These results suggest that the typical deceleration force that occurs during running to a single limb landing can substantially reduce the strain in the ACL relative to conditions associated with an isolated single limb landing from a vertical jump.     

In phantom evaluation of targeting accuracy in MRI-based brain radiosurgery

PurposeTo determine the targeting accuracy of brain radiosurgery when planning procedures employing different MRI and MRI + CT combinations are adopted.Materials and methodA new phantom, the BrainTool, has been designed and realized to test image co-registration and targeting accuracy in a realistic anatomical situation. The phantom was created with a 3D printer and materials that mimic realistic brain MRI and CT contrast using a model extracted from a synthetic MRI study of a human brain. Eight markers distributed within the BrainTool provide for assessment of the accuracy of image registrations while two cavities that host an ionization chamber are used to perform targeting accuracy measurements with an iterative cross-scan method. Two procedures employing 1.5 T MRI-only or a combination of MRI (taken with 1.5 T or 3 T scanners) and CT to carry out Gamma Knife treatments were investigated. As distortions can impact targeting accuracy, MR images were preliminary evaluated to assess image deformation extent using GammaTool phantom.ResultsMR images taken with both scanners showed average and maximum distortion of 0.3 mm and 1 mm respectively. The marker distances in co-registered images resulted below 0.5 mm for both MRI scans. The targeting mismatches obtained were 0.8 mm, 1.0 mm and 1.2 mm for MRI-only and MRI + CT (1,5T and 3 T), respectively.ConclusionsProcedures using a combination of MR and CT images provide targeting accuracies comparable to those of MRI-only procedures. The BrainTool proved to be a suitable tool for carrying out co-registration and targeting accuracy of Gamma Knife brain radiosurgery treatments.