In vivo estimates of axonal stretch and 3D brain deformation during mild head impact

The rapid deformation of brain tissue in response to head impact can lead to traumatic brain injury. In vivo measurements of brain deformation during non-injurious head impacts are necessary to understand the underlying mechanisms of traumatic brain injury and compare to computational models of brain biomechanics. Using tagged magnetic resonance imaging (MRI), we obtained measurements of three-dimensional strain tensors that resulted from a mild head impact after neck rotation or neck extension. Measurements of maximum principal strain (MPS) peaked shortly after impact, with maximal values of 0.019-0.053 that correlated strongly with peak angular velocity. Subject-specific patterns of MPS were spatially heterogeneous and consistent across subjects for the same motion, though regions of high deformation differed between motions. The largest MPS values were seen in the cortical gray matter and cerebral white matter for neck rotation and the brainstem and cerebellum for neck extension. Axonal fiber strain (Ef) was estimated by combining the strain tensor with diffusion tensor imaging data. As with MPS, patterns of Ef varied spatially within subjects, were similar across subjects within each motion, and showed group differences between motions. Values were highest and most strongly correlated with peak angular velocity in the corpus callosum for neck rotation and in the brainstem for neck extension. The different patterns of brain deformation between head motions highlight potential areas of greater risk of injury between motions at higher loading conditions. Additionally, these experimental measurements can be directly compared to predictions of generic or subject-specific computational models of traumatic brain injury.

Multi-Excitation Magnetic Resonance Elastography of the Brain: Wave Propagation in Anisotropic White Matter

Magnetic resonance elastography (MRE) has emerged as a sensitive imaging technique capable of providing a quantitative understanding of neural microstructural integrity. However, a reliable method for the quantification of the anisotropic mechanical properties of human white matter is currently lacking, despite the potential to illuminate the pathophysiology behind neurological disorders and traumatic brain injury. In this study, we examine the use of multiple excitations in MRE to generate wave displacement data sufficient for anisotropic inversion in white matter. We show the presence of multiple unique waves from each excitation which we combine to solve for parameters of an incompressible, transversely isotropic (ITI) material: shear modulus, μ, shear anisotropy, ϕ, and tensile anisotropy, ζ. We calculate these anisotropic parameters in the corpus callosum body and find the mean values as μ = 3.78 kPa, ϕ = 0.151, and ζ = 0.099 (at 50 Hz vibration frequency). This study demonstrates that multi-excitation MRE provides displacement data sufficient for the evaluation of the anisotropic properties of white matter.

Estimation of Anisotropic Material Properties of Soft Tissue by MRI of Ultrasound-Induced Shear Waves

This paper describes a new method for estimating anisotropic mechanical properties of fibrous soft tissue by imaging shear waves induced by focused ultrasound (FUS) and analyzing their direction-dependent speeds. Fibrous materials with a single, dominant fiber direction may exhibit anisotropy in both shear and tensile moduli, reflecting differences in the response of the material when loads are applied in different directions. The speeds of shear waves in such materials depend on the propagation and polarization directions of the waves relative to the dominant fiber direction. In this study, shear waves were induced in muscle tissue (chicken breast) ex vivo by harmonically oscillating the amplitude of an ultrasound beam focused in a cylindrical tissue sample. The orientation of the fiber direction relative to the excitation direction was varied by rotating the sample. Magnetic resonance elastography (MRE) was used to visualize and measure the full 3D displacement field due to the ultrasound-induced shear waves. The phase gradient (PG) of radially propagating "slow" and "fast" shear waves provided local estimates of their respective wave speeds and directions. The equations for the speeds of these waves in an incompressible, transversely isotropic (TI), linear elastic material were fitted to measurements to estimate the shear and tensile moduli of the material. The combination of focused ultrasound and MR imaging allows noninvasive, but comprehensive, characterization of anisotropic soft tissue.

Shear Wave Propagation and Estimation of Material Parameters in a Nonlinear, Fibrous Material

This paper describes the propagation of shear waves in a Holzapfel-Gasser-Ogden (HGO) material and investigates the potential of magnetic resonance elastography (MRE) for estimating parameters of the HGO material model from experimental data. In most MRE studies the behavior of the material is assumed to be governed by linear, isotropic elasticity or viscoelasticity. In contrast, biological tissue is often nonlinear and anisotropic with a fibrous structure. In such materials, application of a quasi-static deformation (predeformation) plays an important role in shear wave propagation. Closed form expressions for shear wave speeds in an HGO material with a single family of fibers were found in a reference (undeformed) configuration and after imposed predeformations. These analytical expressions show that shear wave speeds are affected by the parameters (μ0, k1, k2, κ) of the HGO model and by the direction and amplitude of the predeformations. Simulations of corresponding finite element (FE) models confirm the predicted influence of HGO model parameters on speeds of shear waves with specific polarization and propagation directions. Importantly, the dependence of wave speeds on the parameters of the HGO model and imposed deformations could ultimately allow the noninvasive estimation of material parameters in vivo from experimental shear wave image data.

Quantitative Validation of MRI-Based Motion Estimation for Brain Impact Biomechanics

Head impact can cause traumatic brain injury (TBI) through axonal overstretch or subsequent inflammation and understanding the biomechanics of the impact event is useful for TBI prevention research. Tagged magnetic resonance imaging (MRI) acquired during a mild-acceleration impact has enabled measurement and visualization of brain deformation in vivo. However, measurements using MRI are subject to error, and having independent validation while imaging in vivo is very difficult. Thus, characterizing the accuracy of these measurements needs to be done in a separate experiment using a phantom where a gold standard is available. This study describes a method for error quantification using a calibration phantom compatible with MRI and high-speed video (the gold standard). During linear acceleration, the maximum shear strain (MSS) in the phantom ranged from 0 to 12%, which is similar to in vivo brain deformation at a similar acceleration. The mean displacement error against video was 0.3±0.3 mm, and the MSS error was 1.4±0.3%. To match resolutions, video data was filtered temporally using an averaging filter. Compared to the unfiltered results, resolution matching improved the agreement between MRI and video results by 15%. In conclusion, tagged MRI analysis compares well to video data provided that resolutions are matched-a finding that is also applicable when using MRI to validate simulations.

Statistical Characterization of Human Brain Deformation During Mild Angular Acceleration Measured In Vivo by Tagged Magnetic Resonance Imaging

Understanding of in vivo brain biomechanical behavior is critical in the study of traumatic brain injury (TBI) mechanisms and prevention. Using tagged magnetic resonance imaging, we measured spatiotemporal brain deformations in 34 healthy human volunteers under mild angular accelerations of the head. Two-dimensional (2D) Lagrangian strains were examined throughout the brain in each subject. Strain metrics peaked shortly after contact with a padded stop, corresponding to the inertial response of the brain after head deceleration. Maximum shear strain of at least 3% was experienced at peak deformation by an area fraction (median±standard error) of 23.5±1.8% of cortical gray matter, 15.9±1.4% of white matter, and 4.0±1.5% of deep gray matter. Cortical gray matter strains were greater in the temporal cortex on the side of the initial contact with the padded stop and also in the contralateral temporal, frontal, and parietal cortex. These tissue-level deformations from a population of healthy volunteers provide the first in vivo measurements of full-volume brain deformation in response to known kinematics. Although strains differed in different tissue type and cortical lobes, no significant differences between male and female head accelerations or strain metrics were found. These cumulative results highlight important kinematic features of the brain's mechanical response and can be used to facilitate the evaluation of computational simulations of TBI.

Acoustic trap-and-release for rapid assessment of cell motility

Functional cilia and flagella are crucial to the propulsion of physiological fluids, motile cells, and microorganisms. Motility assessment of individual cells allows discrimination of normal from dysfunctional behavior, but cell-scale analysis of individual trajectories to represent a population is laborious and impractical for clinical, industrial, and even research applications. We introduce an assay that quantifies swimming capability as a function of the variation in polar moment of inertia of cells released from an acoustic trap. Acoustic confinement eliminates the need to trace discrete trajectories and enables automated analysis of hundreds of cells in minutes. The approach closely approximates the average speed estimated from the mean squared displacement of individual cells for wild-type Chlamydomonas reinhardtii and two mutants (ida3 and oda5) that display aberrant swimming behaviors. Large-population acoustic trap-and-release rapidly differentiates these cell types based on intrinsic motility, which provides a highly sensitive and efficient alternative to conventional particle tracing.

Insights Into Traumatic Brain Injury From MRI of Harmonic Brain Motion

Measurements of dynamic deformation of the human brain, induced by external harmonic vibration of the skull, were analyzed to illuminate the mechanics of mild traumatic brain injury (TBI). Shear wave propagation velocity vector fields were obtained to illustrate the role of the skull and stiff internal membranes in transmitting motion to the brain. Relative motion between the cerebrum and cerebellum was quantified to assess the vulnerability of connecting structures. Mechanical deformation was quantified throughout the brain to investigate spatial patterns of strain and axonal stretch. Strain magnitude was generally attenuated as shear waves propagated into interior structures of the brain; this attenuation was greater at higher frequencies. Analysis of shear wave propagation direction indicates that the stiff membranes (falx and tentorium) greatly affect brain deformation during imposed skull motion as they serve as sites for both initiation and reflection of shear waves. Relative motion between the cerebellum and cerebrum was small in comparison with the overall motion of both structures, which suggests that such relative motion might play only a minor role in TBI mechanics. Strain magnitudes and the amount of axonal stretch near the bases of sulci were similar to those in other areas of the cortex, and local strain concentrations at the gray-white matter boundary were not observed. We tentatively conclude that observed differences in neuropathological response in these areas might be due to heterogeneity in the response to mechanical deformation rather than heterogeneity of the deformation itself.

How Does Cilium Length Affect Beating?

The effects of cilium length on the dynamics of cilia motion were investigated by high-speed video microscopy of uniciliated mutants of the swimming alga, Chlamydomonas reinhardtii. Cells with short cilia were obtained by deciliating cells via pH shock and allowing cilia to reassemble for limited times. The frequency of cilia beating was estimated from the motion of the cell body and of the cilium. Key features of the ciliary waveform were quantified from polynomial curves fitted to the cilium in each image frame. Most notably, periodic beating did not emerge until the cilium reached a critical length between 2 and 4 μm. Surprisingly, in cells that exhibited periodic beating, the frequency of beating was similar for all lengths with only a slight decrease in frequency as length increased from 4 μm to the normal length of 10-12 μm. The waveform average curvature (rad/μm) was also conserved as the cilium grew. The mechanical metrics of ciliary propulsion (force, torque, and power) all increased in proportion to length. The mechanical efficiency of beating appeared to be maximal at the normal wild-type length of 10-12 μm. These quantitative features of ciliary behavior illuminate the biophysics of cilia motion and, in future studies, may help distinguish competing hypotheses of the underlying mechanism of oscillation.

A 3D Computational Head Model Under Dynamic Head Rotation and Head Extension Validated Using Live Human Brain Data, Including the Falx and the Tentorium

We employ an advanced 3D computational model of the head with high anatomical fidelity, together with measured tissue properties, to assess the consequences of dynamic loading to the head in two distinct modes: head rotation and head extension. We use a subject-specific computational head model, using the material point method, built from T1 magnetic resonance images, and considering the anisotropic properties of the white matter which can predict strains in the brain under large rotational accelerations. The material model now includes the shear anisotropy of the white matter. We validate the model under head rotation and head extension motions using live human data, and advance a prior version of the model to include biofidelic falx and tentorium. We then examine the consequences of incorporating the falx and tentorium in terms of the predictions from the computational head model.

Frequency-dependent shear properties of annulus fibrosus and nucleus pulposus by magnetic resonance elastography

Aging and degeneration are associated with changes in mechanical properties in the intervertebral disc, generating interest in the establishment of mechanical properties as early biomarkers for the degenerative cascade. Magnetic resonance elastography (MRE) of the intervertebral disc is usually limited to the nucleus pulposus, as the annulus fibrosus is stiffer and less hydrated. The objective of this work was to adapt high-frequency needle MRE to the characterization of the shear modulus of both the nucleus pulposus and annulus fibrosus. Bovine intervertebral discs were removed from fresh oxtails and characterized by needle MRE. The needle was inserted in the center of the disc and vibrations were generated by an amplified piezoelectric actuator. MRE acquisitions were performed on a 4.7-T small-animal MR scanner using a spin echo sequence with sinusoidal motion encoding gradients. Acquisitions were repeated over a frequency range of 1000-1800 Hz. The local frequency estimation inversion algorithm was used to compute the shear modulus. Stiffness maps allowed the visualization of the soft nucleus pulposus surrounded by the stiffer annulus fibrosus surrounded by the homogeneous gel. A significant difference in shear modulus between the nucleus pulposus and annulus fibrosus, and an increase in the shear modulus with excitation frequency, were observed, in agreement with the literature. This study demonstrates that global characterization of both the nucleus pulposus and annulus fibrosus of the intervertebral disc is possible with needle MRE using a preclinical magnetic resonance imaging (MRI) scanner. MRE can be a powerful method for the mapping of the complex properties of the intervertebral disc. The developed method could be adapted for in situ use by preserving adjacent vertebrae and puncturing the side of the intervertebral disc, thereby allowing an assessment of the contribution of osmotic pressure to the mechanical behavior of the intervertebral disc.

Quantifying Ciliary Dynamics during Assembly Reveals Stepwise Waveform Maturation in Airway Cells

Motile cilia are essential for clearance of particulates and pathogens from airways. For effective transport, ciliary motor proteins and axonemal structures interact to generate the rhythmic, propulsive bending, but the mechanisms that produce a dynamic waveform remain incompletely understood. Biomechanical measures of human ciliary motion and their relationships to ciliary assembly are needed to illuminate the biophysics of normal ciliary function and to quantify dysfunction in ciliopathies. To these ends, we analyzed ciliary motion by high-speed video microscopy of ciliated cells sampled from human lung airways compared with primary culture cells that undergo ciliogenesis in vitro. Quantitative assessment of waveform parameters showed variations in waveform shape between individual cilia; however, general trends in waveform parameters emerged, associated with progression of cilia length and stage of differentiation. When cilia emerged from cultured cells, beat frequency was initially elevated, then fell and remained stable as cilia lengthened. In contrast, the average bending amplitude and the ability to generate force gradually increased and eventually approached values observed in ex vivo samples. Dynein arm motor proteins DNAH5, DNAH9, DNAH11, and DNAH6 were localized within specific regions of the axoneme in the ex vivo cells; however, distinct stages of in vitro waveform development identified by biomechanical features were associated with the progressive movement of dyneins to the appropriate proximal or distal sections of the cilium. These observations suggest that the stepwise variation in waveform development during ciliogenesis is dependent on cilia length and potentially on outer dynein arm assembly.

Mechanics of cortical folding: stress, growth and stability

Cortical folding, or gyrification, coincides with several important developmental processes. The folded shape of the human brain allows the cerebral cortex, the thin outer layer of neurons and their associated projections, to attain a large surface area relative to brain volume. Abnormal cortical folding has been associated with severe neurological, cognitive and behavioural disorders, such as epilepsy, autism and schizophrenia. However, despite decades of study, the mechanical forces that lead to cortical folding remain incompletely understood. Leading hypotheses have focused on the roles of (i) tangential growth of the outer cortex, (ii) spatio-temporal patterns in the birth and migration of neurons, and (iii) internal tension in axons. Recent experimental studies have illuminated not only the fundamental cellular and molecular processes underlying cortical development, but also the stress state, mechanical properties and spatio-temporal patterns of growth in the developing brain. The combination of mathematical modelling and physical measurements has allowed researchers to evaluate hypothesized mechanisms of folding, to determine whether each is consistent with physical laws. This review summarizes what physical scientists have learned from models and recent experimental observations, in the context of recent neurobiological discoveries regarding cortical development. Here, we highlight evidence of a combined mechanism, in which spatio-temporal patterns bias the locations of primary folds (i), but tangential growth of the cortical plate induces mechanical instability (ii) to propagate primary and higher-order folds.This article is part of the Theo Murphy meeting issue 'Mechanics of development'.

Quantifying Tensor Field Similarity With Global Distributions and Optimal Transport

Strain tensor fields quantify tissue deformation and are important for functional analysis of moving organs such as the heart and the tongue. Strain data can be readily obtained using medical imaging. However, quantification of similarity between different data sets is difficult. Strain patterns vary in space and time, and are inherently multidimensional. Also, the same type of mechanical deformation can be applied to different shapes; hence, automatic quantification of similarity should be unaffected by the geometry of the objects being deformed. This work introduces the application of global distributions used to classify shapes and vector fields in the pattern recognition literature, in the context of tensorial strain data. In particular, the distribution of mechanical properties of a field are approximated using a 3D histogram, and the Wasserstein distance from optimal transport theory is used to measure the similarity between histograms. To measure the method's consistency in matching deformations across different objects, the proposed approach was evaluated by sorting strain fields according to their similarity. Performance was compared to sorting via maximum shear distribution (a 1D histogram) and tensor residual magnitude (in perfectly registered objects). The technique was also applied to correlate muscle activation to muscular contraction observed via tagged MRI. The results show that the proposed approach accurately matches deformation regardless of the shape of the object being deformed. Sorting accuracy surpassed 1D shear distribution and was on par with residual magnitude, but without the need for registration between objects.

Pre-clinical MR elastography: Principles, techniques, and applications

Magnetic resonance elastography (MRE) is a method for measuring the mechanical properties of soft tissue in vivo, non-invasively, by imaging propagating shear waves in the tissue. The speed and attenuation of waves depends on the elastic and dissipative properties of the underlying material. Tissue mechanical properties are essential for biomechanical models and simulations, and may serve as markers of disease, injury, development, or recovery. MRE is already established as a clinical technique for detecting and characterizing liver disease. The potential of MRE for diagnosing or characterizing disease in other organs, including brain, breast, and heart is an active research area. Studies involving MRE in the pre-clinical setting, in phantoms and artificial biomaterials, in the mouse, and in other mammals, are critical to the development of MRE as a robust, reliable, and useful modality.

Dynamic patterns of cortical expansion during folding of the preterm human brain

During the third trimester of human brain development, the cerebral cortex undergoes dramatic surface expansion and folding. Physical models suggest that relatively rapid growth of the cortical gray matter helps drive this folding, and structural data suggest that growth may vary in both space (by region on the cortical surface) and time. In this study, we propose a unique method to estimate local growth from sequential cortical reconstructions. Using anatomically constrained multimodal surface matching (aMSM), we obtain accurate, physically guided point correspondence between younger and older cortical reconstructions of the same individual. From each pair of surfaces, we calculate continuous, smooth maps of cortical expansion with unprecedented precision. By considering 30 preterm infants scanned two to four times during the period of rapid cortical expansion (28-38 wk postmenstrual age), we observe significant regional differences in growth across the cortical surface that are consistent with the emergence of new folds. Furthermore, these growth patterns shift over the course of development, with noninjured subjects following a highly consistent trajectory. This information provides a detailed picture of dynamic changes in cortical growth, connecting what is known about patterns of development at the microscopic (cellular) and macroscopic (folding) scales. Since our method provides specific growth maps for individual brains, we are also able to detect alterations due to injury. This fully automated surface analysis, based on tools freely available to the brain-mapping community, may also serve as a useful approach for future studies of abnormal growth due to genetic disorders, injury, or other environmental variables.

Relationships between scalp, brain, and skull motion estimated using magnetic resonance elastography

The objective of this study was to characterize the relationships between motion in the scalp, skull, and brain. In vivo estimates of motion transmission from the skull to the brain may illuminate the mechanics of traumatic brain injury. Because of challenges in directly sensing skull motion, it is useful to know how well motion of soft tissue of the head, i.e., the scalp, can approximate skull motion or predict brain tissue deformation. In this study, motion of the scalp and brain were measured using magnetic resonance elastography (MRE) and separated into components due to rigid-body displacement and dynamic deformation. Displacement estimates in the scalp were calculated using low motion-encoding gradient strength in order to reduce "phase wrapping" (an ambiguity in displacement estimates caused by the 2 π-periodicity of MRE phase contrast). MRE estimates of scalp and brain motion were compared to skull motion estimated from three tri-axial accelerometers. Comparison of the relative amplitudes and phases of harmonic motion in the scalp, skull, and brain of six human subjects indicate that data from scalp-based sensors should be used with caution to estimate skull kinematics, but that fairly consistent relationships exist between scalp, skull, and brain motion. In addition, the measured amplitude and phase relationships of scalp, skull, and brain can be used to evaluate and improve mathematical models of head biomechanics.

Mechanical properties of porcine brain tissue in vivo and ex vivo estimated by MR elastography

The mechanical properties of brain tissue in vivo determine the response of the brain to rapid skull acceleration. These properties are thus of great interest to the developers of mathematical models of traumatic brain injury (TBI) or neurosurgical simulations. Animal models provide valuable insight that can improve TBI modeling. In this study we compare estimates of mechanical properties of the Yucatan mini-pig brain in vivo and ex vivo using magnetic resonance elastography (MRE) at multiple frequencies. MRE allows estimations of properties in soft tissue, either in vivo or ex vivo, by imaging harmonic shear wave propagation. Most direct measurements of brain mechanical properties have been performed using samples of brain tissue ex vivo. It has been observed that direct estimates of brain mechanical properties depend on the frequency and amplitude of loading, as well as the time post-mortem and condition of the sample. Using MRE in the same animals at overlapping frequencies, we observe that porcine brain tissue in vivo appears stiffer than porcine brain tissue samples ex vivo at frequencies of 100 Hz and 125 Hz, but measurements show closer agreement at lower frequencies.

Finite element models of flagella with sliding radial spokes and interdoublet links exhibit propagating waves under steady dynein loading

It remains unclear how flagella generate propulsive, oscillatory waveforms. While it is well known that dynein motors, in combination with passive cytoskeletal elements, drive the bending of the axoneme by applying shearing forces and bending moments to microtubule doublets, the origin of rhythmicity is still mysterious. Most conceptual models of flagellar oscillation involve dynein regulation or switching, so that dynein activity first on one side of the axoneme, then the other, drives bending. In contrast, a "viscoelastic flutter" mechanism has recently been proposed, based on a dynamic structural instability. Simple mathematical models of coupled elastic beams in viscous fluid, subjected to steady, axially distributed, dynein forces of sufficient magnitude, can exhibit oscillatory motion without any switching or dynamic regulation. Here we introduce more realistic finite element (FE) models of 6-doublet and 9-doublet flagella, with radial spokes and interdoublet links that slide along the central pair or corresponding doublet. These models demonstrate the viscoelastic flutter mechanism. Above a critical force threshold, these models exhibit an abrupt onset of propulsive, wavelike oscillations typical of flutter instability. Changes in the magnitude and spatial distribution of steady dynein force, or to viscous resistance, lead to behavior qualitatively consistent with experimental observations. This study demonstrates the ability of FE models to simulate nonlinear interactions between axonemal components during flagellar beating, and supports the plausibility of viscoelastic flutter as a mechanism of flagellar oscillation.

Measurement of anisotropic mechanical properties in porcine brain white matter ex vivo using magnetic resonance elastography

The mechanical properties of brain tissue, particularly those of white matter (WM), need to be characterized accurately for use in finite element (FE) models of brain biomechanics and traumatic brain injury (TBI). Magnetic resonance elastography (MRE) is a powerful tool for non-invasive estimation of the mechanical properties of soft tissues. While several studies involving direct mechanical tests of brain tissue have shown mechanical anisotropy, most MRE studies of brain tissue assume an isotropic model. In this study, an incompressible transversely isotropic (TI) material model parameterized by minimum shear modulus (μ₂), shear anisotropy parameter (ϕ), and tensile anisotropy parameter (ζ) is applied to analyze MRE measurements of ex vivo porcine white matter (WM) brain tissue. To characterize shear anisotropy, "slow" (pure transverse) shear waves were propagated at 100, 200 and 300Hz through sections of ex vivo brain tissue including both WM and gray matter (GM). Shear waves were found to propagate with elliptical fronts, consistent with TI material behavior. Shear wave fields were also analyzed within regions of interest (ROI) to find local shear wavelengths parallel and perpendicular to fiber orientation. FE simulations of a TI material with a range of plausible shear modulus (μ₂) and shear anisotropy parameters (ϕ) were run and the results were analyzed in the same fashion as the experimental case. Parameters of the FE simulations which most closely matched each experiment were taken to represent the mechanical properties of that particular sample. Using this approach, WM in the ex vivo porcine brain was found to be mildly anisotropic in shear with estimates of minimum shear modulus (actuation frequencies listed in parenthesis): μ₂= 1.04 ± 0.12 kPa (at 100Hz), μ₂= 1.94 ± 0.29 kPa (at 200Hz), and μ₂= 2.88 ± 0.34 kPa (at 300Hz) and corresponding shear anisotropy factors of ϕ= 0.27 ± 0.09 (at 100Hz), ϕ= 0.29 ± 0.14 (at 200Hz) and ϕ= 0.34 ± 0.13 (at 300Hz). Future MRE studies will focus on tensile anisotropy, which will require both slow and fast shear waves for accurate estimation.

The Relationship of Three-Dimensional Human Skull Motion to Brain Tissue Deformation in Magnetic Resonance Elastography Studies

In traumatic brain injury (TBI), membranes such as the dura mater, arachnoid mater, and pia mater play a vital role in transmitting motion from the skull to brain tissue. Magnetic resonance elastography (MRE) is an imaging technique developed for noninvasive estimation of soft tissue material parameters. In MRE, dynamic deformation of brain tissue is induced by skull vibrations during magnetic resonance imaging (MRI); however, skull motion and its mode of transmission to the brain remain largely uncharacterized. In this study, displacements of points in the skull, reconstructed using data from an array of MRI-safe accelerometers, were compared to displacements of neighboring material points in brain tissue, estimated from MRE measurements. Comparison of the relative amplitudes, directions, and temporal phases of harmonic motion in the skulls and brains of six human subjects shows that the skull-brain interface significantly attenuates and delays transmission of motion from skull to brain. In contrast, in a cylindrical gelatin "phantom," displacements of the rigid case (reconstructed from accelerometer data) were transmitted to the gelatin inside (estimated from MRE data) with little attenuation or phase lag. This quantitative characterization of the skull-brain interface will be valuable in the parameterization and validation of computer models of TBI.

Multimodal surface matching with higher-order smoothness constraints

In brain imaging, accurate alignment of cortical surfaces is fundamental to the statistical sensitivity and spatial localisation of group studies, and cortical surface-based alignment has generally been accepted to be superior to volume-based approaches at aligning cortical areas. However, human subjects have considerable variation in cortical folding, and in the location of functional areas relative to these folds. This makes alignment of cortical areas a challenging problem. The Multimodal Surface Matching (MSM) tool is a flexible, spherical registration approach that enables accurate registration of surfaces based on a variety of different features. Using MSM, we have previously shown that driving cross-subject surface alignment, using areal features, such as resting state-networks and myelin maps, improves group task fMRI statistics and map sharpness. However, the initial implementation of MSM's regularisation function did not penalize all forms of surface distortion evenly. In some cases, this allowed peak distortions to exceed neurobiologically plausible limits, unless regularisation strength was increased to a level which prevented the algorithm from fully maximizing surface alignment. Here we propose and implement a new regularisation penalty, derived from physically relevant equations of strain (deformation) energy, and demonstrate that its use leads to improved and more robust alignment of multimodal imaging data. In addition, since spherical warps incorporate projection distortions that are unavoidable when mapping from a convoluted cortical surface to the sphere, we also propose constraints that enforce smooth deformation of cortical anatomies. We test the impact of this approach for longitudinal modelling of cortical development for neonates (born between 31 and 43 weeks of post-menstrual age) and demonstrate that the proposed method increases the biological interpretability of the distortion fields and improves the statistical significance of population-based analysis relative to other spherical methods.