Publications by Author: Curtis L. Johnson

Introduction

Computational head injury models are promising tools for understanding and predicting traumatic brain injuries. However, most available head injury models are “average” models that employ a single set of head geometry (e.g., 50th-percentile U.S. male) without considering variability in these parameters across the human population. A significant variability of head shapes exists in U.S. Army soldiers, evident from the Anthropometric Survey of U.S. Army Personnel (ANSUR II). The objective of this study is to elucidate the effects of head shape on the predicted risk of traumatic brain injury from computational head injury models.

Materials and Methods

Magnetic resonance imaging scans of 25 human subjects are collected. These images are registered to the standard MNI152 brain atlas, and the resulting transformation matrix components (called head shape parameters) are used to quantify head shapes of the subjects. A generative machine learning model is used to generate 25 additional head shape parameter datasets to augment our database. Head injury models are developed for these head shapes, and a rapid injurious head rotation event is simulated to obtain several brain injury predictor variables (BIPVs): Peak cumulative maximum principal strain (CMPS), average CMPS, and the volume fraction of brain exceeding an injurious CMPS threshold. A Gaussian process regression model is trained between head shape parameters and BIPVs, which is then used to study the relative sensitivity of the various BIPVs on individual head shape parameters. We distinguish head shape parameters into 2 types: Scaling components Txx⁠, Tyy⁠, and Tzz that capture the breadth, length, and height of the head, respectively, and shearing components (⁠Txy,Txz,Tyx,Tyz,Tzx⁠, and Tzy⁠) that capture the relative skewness of the head shape.

Results

An overall positive correlation is evident between scaling components and BIPVs. Notably, a very high, positive correlation is seen between the BIPVs and the head volume. As an example, a 57% increase in peak CMPS was noted between the smallest and the largest investigated head volume parameters. The variation in shearing components Txy,Txz,Tyx,Tyz,Tzx⁠, and Tzy on average does not cause notable changes in the BIPVs. From the Gaussian process regression model, all 3 BIPVs showed an increasing trend with each of the 3 scaling components, but the BIPVs are found to be most sensitive to the height dimension of the head. From the Sobol sensitivity analysis, the Tzz scaling parameter contributes nearly 60% to the total variance in peak and average CMPS; Tyy contributes approximately 20%, whereas Txx contributes less than 5%. The remaining contribution is from the 6 shearing components. Unlike peak and average CMPS, the VF-CMPS BIPV is associated with relatively evenly distributed Sobol indices across the 3 scaling parameters. Furthermore, the contribution of shearing components on the total variance in this case is negligible.

Conclusions

Head shape has a considerable influence on the injury predictions of computational head injury models. Available “average” head injury models based on a 50th-percentile U.S. male are likely associated with considerable uncertainty. In general, larger head sizes correspond to greater BIPV magnitudes, which point to potentially a greater injury risk under rapid neck rotation for people with larger heads.

Strain energy and kinetic energy in the human brain were estimated by magnetic resonance elastography (MRE) during harmonic excitation of the head, and compared to characterize the effect of loading direction and frequency on brain deformation. In brain MRE, shear waves are induced by external vibration of the skull and imaged by a modified MR imaging sequence; the resulting harmonic displacement fields are typically “inverted” to estimate mechanical properties, like stiffness or damping. However, measurements of tissue motion from MRE also illuminate key features of the response of the brain to skull loading. In this study, harmonic excitation was applied in two different directions and at five different frequencies from 20 to 90 Hz. Lateral loading induced primarily left-right head motion and rotation in the axial plane; occipital loading induced anterior-posterior head motion and rotation in the sagittal plane. The ratio of strain energy to kinetic energy (SE/KE) depended strongly on both direction and frequency. The ratio of SE/KE was approximately four times larger for lateral excitation than for occipital excitation and was largest at the lowest excitation frequencies studied. These results are consistent with clinical observations that suggest lateral impacts are more likely to cause injury than occipital or frontal impacts, and also with observations that the brain has low-frequency (∼10 Hz) natural modes of oscillation. The SE/KE ratio from brain MRE is potentially a simple and powerful dimensionless metric of brain vulnerability to deformation and injury.

Magnetic Resonance Elastography (MRE) is a noninvasive method for quantitatively assessing the viscoelastic properties of tissues, such as the brain. MRE has been successfully used to measure the material properties and diagnose diseases based on the difference in mechanical properties between diseased and normal tissue. However, MRE is still an emerging technology that is not part of routine clinical imaging like structural Magnetic Resonance Imaging (MRI), and the acquisition equipment is not widely available. Thus, it is challenging to collect MRE, but there is an increasing interest in it. In this study, we explore using structural MRI images to synthesize the MRE-derived material properties of the human brain. We use deep networks that employ both MRI and Diffusion Tensor Imaging (DTI) to explore the best input images for MRE image synthesis. This work is the first study to report on the feasibility of MRE synthesis from structural MRI and DTI.

Purpose

To determine how the biomechanical vulnerability of the human brain is affected by features of individual anatomy and loading.

Methods

To identify the features that contribute most to brain vulnerability, we imparted mild harmonic acceleration to the head and measured the resulting brain motion and deformation using magnetic resonance elastography (MRE). Oscillatory motion was imparted to the heads of adult participants using a lateral actuator (n = 24) or occipital actuator (n = 24) at 20 Hz, 30 Hz, and 50 Hz. Displacement vector fields and strain tensor fields in the brain were obtained from MRE measurements. Anatomical images, as well as displacement and strain fields from each participant were rigidly and deformably aligned to a common atlas (MNI-152). Vulnerability of the brain to deformation was quantified by the ratio of strain energy (SE) to kinetic energy (KE) for each participant. Similarity of deformation patterns between participants was quantified using strain field correlation (C_V). Linear regression models were used to identify the effect of similarity of brain size, shape, and age, as well as similarity of loading, on C_V.

Results

The SE/KE ratio decreased with frequency and was larger for participants undergoing lateral, rather than occipital, actuation. Head rotation about the inferior–superior axis was correlated with larger SE/KE ratio. Strain field correlations were primarily affected by the similarity of rigid-body motion.

Conclusion

The motion applied to the skull is the most important factor in determining both the vulnerability of the brain to deformation and the similarity between strain fields in different individuals.

Noninvasive measurements of brain deformation in human participants in vivo are needed to develop models of brain biomechanics and understand traumatic brain injury (TBI). Tagged magnetic resonance imaging (tagged MRI) and magnetic resonance elastography (MRE) are two techniques to study human brain deformation; these techniques differ in the type of motion and difficulty of implementation. In this study, oscillatory strain fields in the human brain caused by impulsive head acceleration and measured by tagged MRI were compared quantitatively to strain fields measured by MRE during harmonic head motion at 10 and 50 Hz. Strain fields were compared by registering to a common anatomical template, then computing correlations between the registered strain fields. Correlations were computed between tagged MRI strain fields in six participants and MRE strain fields at 10 Hz and 50 Hz in six different participants. Correlations among strain fields within the same experiment type were compared statistically to correlations from different experiment types. Strain fields from harmonic head motion at 10 Hz imaged by MRE were qualitatively and quantitatively similar to modes excited by impulsive head motion, imaged by tagged MRI. Notably, correlations between strain fields from 10 Hz MRE and tagged MRI did not differ significantly from correlations between strain fields from tagged MRI. These results suggest that low-frequency modes of oscillation dominate the response of the brain during impact. Thus, low-frequency MRE, which is simpler and more widely available than tagged MRI, can be used to illuminate the brain's response to head impact.

Traumatic brain injuries (TBIs) occur from rapid head motion that results in brain deformation. Computational models are typically used to estimate brain deformation to predict risk of injury and evaluate the effectiveness of safety countermeasures. The accuracy of these models relies on validation to experimental brain deformation data. In this study, we create the first group-average biomechanical responses of the brain, including structure, material properties, and deformation response, by age and sex from 157 subjects. Subjects were sorted intro three age groups—young, mid-age, and older—and by sex to create group-average neuroanatomy, material properties, and brain deformation response to non-injurious loading using structural and specialized magnetic resonance imaging data. Computational models were also built using the group-average geometry and material properties for each of the six groups. The material properties did not depend on sex, but showed a decrease in shear stiffness in the older adult group. The brain deformation response also showed differences in the distribution of strain and a decrease in the magnitude of maximum strain in the older adult group. The computational models were simulated using the same non-injurious loading conditions as the subject data. While the models’ strain response showed differences among the models, there were no clear relationships with age. Further studies, both modeling and experimental, with more data from subjects in each age group, are needed to clarify the mechanisms underlying the observed changes in strain response with age, and for computational models to better match the trends observed across the group-average responses.