Lateral impacts correlate with falx cerebri displacement and corpus callosum trauma in sports-related concussions

Role of player-specific white matter parcellation and scaling in impact-induced strain responses

Head finite element models (hFEMs) are valuable in understanding injury mechanisms in head impacts. Personalizing hFEMs is important for capturing individualized brain responses, with brain volume scaling proving effective. However, the role of refined white matter (WM) parcellation in hFEMs for evaluating brain strain responses, particularly important in the context of subconcussive head impacts (SHIs) often assessed through changes in WM integrity, remains relatively underexplored. This study evaluated the effect of refined subject-specific WM parcellation in 34 WM segments on responses variability due to brain volume variations, using peak maximum principal strain (95MPS) and strain rate (95MPSr) as injury predictive metrics. Data from diffusion-weighted imaging of 21 Canadian varsity football players were utilized to personalize 21 hFEMs. Simulating four different head impacts, representing 50th and 99th percentile resultant accelerations in frontal and angled-top-right directions, refined player-specific WM parcellation better captured variability of strain responses compared to baseline parcellation. Up to 75.71% of 95MPS and 77.14% of 95MPSr responses were deemed different across refined WM segments for players, compared to a maximum of 16.19% of responses with baseline parcellation. These results suggest that player-specific refined WM parcellation improves the ability to capture player-specific responses. Both impact direction and intensity influenced variations in strain response, with angled-top head impacts combined with high intensity showing greater player-specificity compared to lower intensity and frontal head impacts. These findings highlight the potential benefit of model scaling along with player-specific refined WM parcellation in hFEMs for comprehensively evaluating strain responses. Detailed WM parcellation in hFEMs is important for comprehensive injury assessment, enhancing the alignment of hFEMs with imaging studies evaluating changes in WM integrity across segments. The simple and straightforward method presented herein to achieve player-specific strain response is promising for future SHI studies.

Differences between two maximal principal strain rate calculation schemes in traumatic brain analysis with in-vivo and in-silico datasets

Brain deformation caused by a head impact leads to traumatic brain injury (TBI). The maximum principal strain (MPS) was used to measure the extent of brain deformation and predict injury, and the recent evidence has indicated that incorporating the maximum principal strain rate (MPSR) and the product of MPS and MPSR, denoted as MPS × SR, enhances the accuracy of TBI prediction. However, ambiguities have arisen about the calculation of MPSR. Two schemes have been utilized: one is to use the time derivative of MPS (MPSR1), and another is to use the first eigenvalue of the strain rate tensor (MPSR2). Both MPSR1 and MPSR2 have been applied in previous studies to predict TBI. To quantify the discrepancies between these two methodologies, we compared them across eight in-vivo and one in-silico head impact datasets and found that 95MPSR1 was slightly larger than 95MPSR2 and 95MPS × SR1 was 4.85 % larger than 95MPS × SR2 in average. Across every element in all head impacts, the average MPSR1 was 12.73 % smaller than MPSR2, and MPS × SR1 was 11.95 % smaller than MPS × SR2. Furthermore, logistic regression models were trained to predict TBI using MPSR (or MPS × SR), and no significant difference was observed in the predictability. The consequence of misuse of MPSR and MPS × SR thresholds (i.e. compare threshold of 95MPSR1 with value from 95MPSR2 to determine if the impact is injurious) was investigated, and the resulting false rates were found to be around 1 %. The evidence suggested that these two methodologies were not significantly different in detecting TBI.

The White Matter Fiber Tract Deforms Most in the Perpendicular Direction During In Vivo Volunteer Impacts

White matter (WM) tract-related strains are increasingly used to quantify brain mechanical responses, but their dynamics in live human brains during in vivo impact conditions remain largely unknown. Existing research primarily looked into the normal strain along the WM fiber tracts (i.e., tract-oriented normal strain), but it is rarely the case that the fiber tract only endures tract-oriented normal strain during impacts. In this study, we aim to extend the in vivo measurement of WM fiber deformation by quantifying the normal strain perpendicular to the fiber tract (i.e., tract-perpendicular normal strain) and the shear strain along and perpendicular to the fiber tract (i.e., tract-oriented shear strain and tract-perpendicular shear strain, respectively). To achieve this, we combine the three-dimensional strain tensor from the tagged magnetic resonance imaging with the diffusion tensor imaging (DTI) from an open-access dataset, including 44 volunteer impacts under two head loading modes, i.e., neck rotations (N = 30) and neck extensions (N = 14). The strain tensor is rotated to the coordinate system with one axis aligned with DTI-revealed fiber orientation, and then four tract-related strain measures are calculated. The results show that tract-perpendicular normal strain peaks are the largest among the four strain types (p < 0.05, Friedman's test). The distribution of tract-related strains is affected by the head loading mode, of which laterally symmetric patterns with respect to the midsagittal plane are noted under neck extensions, but not under neck rotations. Our study presents a comprehensive in vivo strain quantification toward a multifaceted understanding of WM dynamics. We find that the WM fiber tract deforms most in the perpendicular direction, illuminating new fundamentals of brain mechanics. The reported strain images can be used to evaluate the fidelity of computational head models, especially those intended to predict fiber deformation under noninjurious conditions.

Head Impact Kinematics and Brain Tissue Strains in High School Lacrosse

Male lacrosse and female lacrosse have differences in history, rules, and equipment. There is current debate regarding the need for enhanced protective headwear in female lacrosse like that worn by male lacrosse players. To inform this discussion, 17 high school lacrosse players (6 female and 11 male) wore the Stanford Instrumented Mouthguard during 26 competitive games over the 2021 season. Time-windowing and video review were used to remove false-positive recordings and verify head acceleration events (HAEs). The HAE rate in high school female lacrosse (0.21 per athlete exposure and 0.24 per player hour) was approximately 35% lower than the HAE rate in high school male lacrosse (0.33 per athlete exposure and 0.36 per player hour). Previously collected kinematics data from the 2019 high school male and female lacrosse season were combined with the newly collected 2021 kinematics data, which were used to drive a finite element head model and simulate 42 HAEs. Peak linear acceleration (PLA), peak angular velocity (PAV), and 95th percentile maximum principal strain (MPS95) of brain tissue were compared between HAEs in high school female and male lacrosse. Median values for peak kinematics and MPS95 of HAEs in high school female lacrosse (PLA, 22.3 g; PAV, 10.4 rad/s; MPS95, 0.05) were lower than for high school male lacrosse (PLA, 24.2 g; PAV, 15.4 rad/s; MPS95, 0.07), but the differences were not statistically significant. Quantifying a lower HAE rate in high school female lacrosse compared to high school male lacrosse, but similar HAE magnitudes, provides insight into the debate regarding helmets in female lacrosse. However, due to the small sample size, additional video-verified data from instrumented mouthguards are required.

Biomechanical Response of Head Surrogate With and Without the Helmet

Measurements of brain deformations under injurious loading scenarios are actively sought. In this work, we report experimentally measured head kinematics and corresponding dynamic, two-dimensional brain simulant deformations in head surrogates under a blunt impact, with and without a helmet. Head surrogates used in this work consisted of skin, skull, dura, falx, tentorium, and brain stimulants. The head surrogate geometry was based on the global human body models consortium's head model. A base head surrogate consisting of skin-skull-brain was considered. In addition, the response of two other head surrogates, skin-skull-dura-brain, and skin-skull-dura-brain-falx-tentorium, was investigated. Head surrogate response was studied for sagittal and coronal plane rotations for impactor velocities of 1 and 3 m/s. Response of head surrogates was compared against strain measurements in PMHS. The strain pattern in the brain simulant was heterogenous, and peak strains were established within ∼30 ms. The choice of head surrogate affect the spatiotemporal evolution of strain. For no helmet case, peak MPS of ∼50-60% and peak MSS of ∼35-50% were seen in brain simulant corresponding to peak rotational accelerations of ∼5000-7000 rad/s2. Peak head kinematics and peak MPS have been reduced by up to 75% and 45%, respectively, with the conventional helmet and by up to 90% and 85%, respectively, with the helmet with antirotational pads. Overall, these results provide important, new data on brain simulant strains under a variety of loading scenarios-with and without the helmets.

Head injury: Importance of the deep brain nuclei in force transmission to the brain

Finite element modeling provides a digital representation of the human body. It is currently the most pertinent method to study the mechanisms of head injury, and is becoming a scientific reference in forensic expert reports. Improved biofidelity is a recurrent aim of research studies in biomechanics in order to improve earlier models whose mechanical properties conformed to simplified elastic behavior and mechanic laws. We aimed to study force transmission to the brain following impacts to the head, using a finite element head model with increased biofidelity. To the model developed by the Laboratory of Applied Biomechanics of Marseille, we added new brain structures (thalamus, central gray nuclei and ventricular systems) as well as three tracts involved in the symptoms of head injury: the corpus callosum, uncinate tracts and corticospinal tracts. Three head impact scenarios were simulated: an uppercut with the prior model and an uppercut with the improved model in order to compare the two models, and a lateral impact with an impact velocity of 6.5 m/s in the improved model. In these conditions, in uppercuts the maximum stress values did not exceed the injury risk threshold. On the other hand, the deep gray matter (thalamus and central gray nuclei) was the region at highest risk of injury during lateral impacts. Even if injury to the deep gray matter is not immediately life-threatening, it could explain the chronic disabling symptoms of even low-intensity head injury.

Regional brain strain dependance on direction of head rotation

Brain injuries in automated vehicles during crash events are likely to include mechanisms of head impact in non-standard positions and postures (i.e., occupants not facing forward in an upright position). Federal regulations currently focus on impact conditions in primary planes of motion, such as frontal or rear impacts (sagittal plane of motion) or side impact (coronal plane of motion) and do not account for out of position occupants or non-standard postures. The objective of the present study was to develop and use the anatomically accurate brain finite element model to parametrically determine the injury metrics under different vectors with head rotation. A custom developed brain finite element model with anatomical accuracy and several anatomical regions defined was used to evaluate whole-brain strain as well as regional brain strain. Cumulative Strain Damage Measure (CSDM) at a threshold of 20% strain and the 95th percentile of the maximum principal strain (MPS95) were calculated for the whole brain and each brain region under multiple rotational directions. The model was exposed to a sinusoidal angular acceleration pulse of 5000 rad per second squared (rad/s2-) over 12.5 ms. The same pulse was used in the primary axes of motion and (lateral bending, flexion, extension, axial rotation) and combined axes representing oblique flexion and oblique extension. Whole brain CSDM20 was highest for lateral bending. Whole brain MPS95 was highest for axial rotation. The rCSDM20 was more susceptible to impact direction, with several brain regions having substantial accumulation of strain for oblique flexion and lateral bending. Comparatively, rMPS95 was more consistent across all rotation directions. The present study quantified the regional brain strain response under multiple rotational vectors identifying a high amount of variability in the accumulation of strain (i.e., CSDM20) in the hypothalamus, hippocampus, and midbrain specifically. While there was a high amount of variability in the accumulation of strain for multiple regions, the maximum strain measured (i.e., MPS95) in the regions was more consistent.

Regional Strain Response of an Anatomically Accurate Human Finite Element Head Model Under Frontal Versus Lateral Loading

INTRODUCTION

Because brain regions are responsible for specific functions, regional damage may cause specific, predictable symptoms. However, the existing brain injury criteria focus on whole brain response. This study developed and validated a detailed human brain computational model with sufficient fidelity to include regional components and demonstrate its feasibility to obtain region-specific brain strains under selected loading.

METHODS

Model development used the Simulated Injury Monitor (SIMon) model as a baseline. Each SIMon solid element was split into 8, with each shell element split into 4. Anatomical regions were identified from FreeSurfer fsaverage neuroimaging template. Material properties were obtained from literature. The model was validated against experimental intracranial pressure, brain-skull displacement, and brain strain data. Model simulations used data from laboratory experiments with a rigid arm pendulum striking a helmeted head-neck system. Data from impact tests (6 m/s) at 2 helmet sites (front and left) were used.

RESULTS

Model validation showed good agreement with intracranial pressure response, fair to good agreement with brain-skull displacement, and good agreement for brain strain. CORrelation Analysis scores were between 0.72 and 0.93 for both maximum principal strain (MPS) and shear strain. For frontal impacts, regional MPS was between 0.14 and 0.36 (average of left and right hemispheres). For lateral impacts, MPS was between 0.20 and 0.48 (left hemisphere) and between 0.22 and 0.51 (right hemisphere). For frontal impacts, regional cumulative strain damage measure (CSDM20) was between 0.01 and 0.87. For lateral impacts, CSDM20 was between 0.36 and 0.99 (left hemisphere) and between 0.09 and 0.93 (right hemisphere).

CONCLUSIONS

Recognizing that neural functions are related to anatomical structures and most model-based injury metrics focus on whole brain response, this study developed an anatomically accurate human brain model to capture regional responses. Model validation was comparable with current models. The model provided sufficient anatomical detail to describe brain regional responses under different impact conditions.

Microstructural Alterations in Tract Development in College Football and Volleyball Players: A Longitudinal Diffusion MRI Study

BACKGROUND AND OBJECTIVES

Repeated impacts in high-contact sports such as American football can affect the brain's microstructure, which can be studied using diffusion MRI. Most imaging studies are cross-sectional, do not include low-contact players as controls, or lack advanced tract-specific microstructural metrics. We aimed to investigate longitudinal changes in high-contact collegiate athletes compared with low-contact controls using advanced diffusion MRI and automated fiber quantification.

METHODS

We examined brain microstructure in high-contact (football) and low-contact (volleyball) collegiate athletes with up to 4 years of follow-up. Inclusion criteria included university and team enrollment. Exclusion criteria included history of neurosurgery, severe brain injury, and major neurologic or substance abuse disorder. We investigated diffusion metrics along the length of tracts using nested linear mixed-effects models to ascertain the acute and chronic effects of subconcussive and concussive impacts, and associations between diffusion changes with clinical, behavioral, and sports-related measures.

RESULTS

Forty-nine football and 24 volleyball players (271 total scans) were included. Football players had significantly divergent trajectories in multiple microstructural metrics and tracts. Longitudinal increases in fractional anisotropy and axonal water fraction, and decreases in radial/mean diffusivity and orientation dispersion index, were present in volleyball but absent in football players (all findings |T-statistic|> 3.5, p value <0.0001). This pattern was present in the callosum forceps minor, superior longitudinal fasciculus, thalamic radiation, and cingulum hippocampus. Longitudinal differences were more prominent and observed in more tracts in concussed football players (n = 24, |T|> 3.6, p < 0.0001). An analysis of immediate postconcussion scans (n = 12) demonstrated a transient localized increase in axial diffusivity and mean/radial kurtosis in the uncinate and cingulum hippocampus (|T| > 3.7, p < 0.0001). Finally, within football players, those with high position-based impact risk demonstrated increased intracellular volume fraction longitudinally (T = 3.6, p < 0.0001).

DISCUSSION

The observed longitudinal changes seen in football, and especially concussed athletes, could reveal diminished myelination, altered axonal calibers, or depressed pruning processes leading to a static, nondecreasing axonal dispersion. This prospective longitudinal study demonstrates divergent tract-specific trajectories of brain microstructure, possibly reflecting a concussive and repeated subconcussive impact-related alteration of white matter development in football athletes.

Current State of Diffusion-Weighted Imaging and Diffusion Tensor Imaging for Traumatic Brain Injury Prognostication

Advanced imaging techniques are needed to assist in providing a prognosis for patients with traumatic brain injury (TBI), particularly mild TBI (mTBI). Diffusion tensor imaging (DTI) is one promising advanced imaging technique, but has shown variable results in patients with TBI and is not without limitations, especially when considering individual patients. Efforts to resolve these limitations are being explored and include developing advanced diffusion techniques, creating a normative database, improving study design, and testing machine learning algorithms. This article will review the fundamentals of DTI, providing an overview of the current state of its utility in evaluating and providing prognosis in patients with TBI.

Interhemispheric transfer time and concussion in adolescents: A longitudinal study using response time and event-related potential measures

Introduction

Concussion in children and adolescents is a public health concern with higher concussion incidence than adults and increased susceptibility to axonal injury. The corpus callosum is a vulnerable location of concussion-related white matter damage that can be associated with short- and long-term effects of concussion. Interhemispheric transfer time (IHTT) of visual information across the corpus callosum can be used as a direct measure of corpus callosum functioning that may be impacted by adolescent concussion with slower IHTT relative to matched controls. Longitudinal studies and studies testing physiological measures of IHTT following concussion in adolescents are lacking.

Methods

We used the N1 and P1 components of the scalp-recorded brain event-related potential (ERP) to measure IHTT in 20 adolescents (ages 12–19 years old) with confirmed concussion and 16 neurologically-healthy control participants within 3 weeks of concussion (subacute stage) and approximately 10 months after injury (longitudinal).

Results

Separate two-group (concussion, control) by two-time (3 weeks, 10 months) repeated measures ANOVAs on difference response times and IHTT latencies of the P1 and N1 components showed no significant differences by group (ps ≥ 0.25) nor by time (ps ≥ 0.64), with no significant interactions (ps ≥ 0.15).

Discussion

Results from the current sample suggest that measures of IHTT may not be strongly influenced at 3 weeks or longitudinally following adolescent concussion using the current IHTT paradigm.

Volumetric MRI Findings in Mild Traumatic Brain Injury (mTBI) and Neuropsychological Outcome

Region of interest (ROI) volumetric assessment has become a standard technique in quantitative neuroimaging. ROI volume is thought to represent a coarse proxy for making inferences about the structural integrity of a brain region when compared to normative values representative of a healthy sample, adjusted for age and various demographic factors. This review focuses on structural volumetric analyses that have been performed in the study of neuropathological effects from mild traumatic brain injury (mTBI) in relation to neuropsychological outcome. From a ROI perspective, the probable candidate structures that are most likely affected in mTBI represent the target regions covered in this review. These include the corpus callosum, cingulate, thalamus, pituitary-hypothalamic area, basal ganglia, amygdala, and hippocampus and associated structures including the fornix and mammillary bodies, as well as whole brain and cerebral cortex along with the cerebellum. Ventricular volumetrics are also reviewed as an indirect assessment of parenchymal change in response to injury. This review demonstrates the potential role and limitations of examining structural changes in the ROIs mentioned above in relation to neuropsychological outcome. There is also discussion and review of the role that post-traumatic stress disorder (PTSD) may play in structural outcome in mTBI. As emphasized in the conclusions, structural volumetric findings in mTBI are likely just a single facet of what should be a multimodality approach to image analysis in mTBI, with an emphasis on how the injury damages or disrupts neural network integrity. The review provides an historical context to quantitative neuroimaging in neuropsychology along with commentary about future directions for volumetric neuroimaging research in mTBI.

Finding the Spatial Co-Variation of Brain Deformation With Principal Component Analysis

OBJECTIVE

Strain and strain rate are effective traumatic brain injury metrics. In finite element (FE) head model, thousands of elements were used to represent the spatial distribution of these metrics. Owing that these metrics are resulted from brain inertia, their spatial distribution can be represented in more concise pattern. Since head kinematic features and brain deformation vary largely across head impact types (Zhan et al., 2021), we applied principal component analysis (PCA) to find the spatial co-variation of injury metrics (maximum principal strain (MPS), MPS rate (MPSR) and MPS × MPSR) in four impact types: simulation, football, mixed martial arts and car crashes, and used the PCA to find patterns in these metrics and improve the machine learning head model (MLHM).

METHODS

We applied PCA to decompose the injury metrics for all impacts in each impact type, and investigate the spatial co-variation using the first principal component (PC1). Furthermore, we developed a MLHM to predict PC1 and then inverse-transform to predict for all brain elements. The accuracy, the model complexity and the size of training dataset of PCA-MLHM are compared with previous MLHM (Zhan et al., 2021).

RESULTS

PC1 explained variance on the datasets. Based on PC1 coefficients, the corpus callosum and midbrain exhibit high variance on all datasets. Finally, the PCA-MLHM reduced model parameters by 74% with a similar MPS estimation accuracy.

CONCLUSION

The brain injury metric in a dataset can be decomposed into mean components and PC1 with high explained variance.

SIGNIFICANCE

The spatial co-variation analysis enables better interpretation of the patterns in brain injury metrics. It also improves the efficiency of MLHM.

Concussion Prone Scenarios: A Multi-Dimensional Exploration in Impact Directions, Brain Morphology, and Network Architectures Using Computational Models

While individual susceptibility to traumatic brain injury (TBI) has been speculated, past work does not provide an analysis considering how physical features of an individual's brain (e.g., brain size, shape), impact direction, and brain network features can holistically contribute to the risk of suffering a TBI from an impact. This work investigated each of these features simultaneously using computational modeling and analyses of simulated functional connectivity. Unlike the past studies that assess the severity of TBI based on the quantification of brain tissue damage (e.g., principal strain), we approached the brain as a complex network in which neuronal oscillations orchestrate to produce normal brain function (estimated by functional connectivity) and, to this end, both the anatomical damage location and its topological characteristics within the brain network contribute to the severity of brain function disruption and injury. To represent the variations in the population, we analyzed a publicly available database of brain imaging data and selected five distinct network architectures, seven different brain sizes, and three uniaxial head rotational conditions to study the consequences of 74 virtual impact scenarios. Results show impact direction produces the most significant change in connections across brain areas (structural connectome) and the functional coupling of activity across these brain areas (functional connectivity). Axial rotations were more injurious than those with sagittal and coronal rotations when the head kinematics were the same for each condition. When the impact direction was held constant, brain network architecture showed a significantly different vulnerability across axial and sagittal, but not coronal rotations. As expected, brain size significantly affected the expected change in structural and functional connectivity after impact. Together, these results provided groupings of predicted vulnerability to impact-a subgroup of male brain architectures exposed to axial impacts were most vulnerable, while a subgroup of female brain architectures was the most tolerant to the sagittal impacts studied. These findings lay essential groundwork for subject-specific analyses of concussion and provide invaluable guidance for designing personalized protection equipment.

Assessment of brain injury characterization and influence of modeling approaches

In this study, using computational biomechanics models, we investigated influence of the skull-brain interface modeling approach and the material property of cerebrum on the kinetic, kinematic and injury outputs. Live animal head impact tests of different severities were reconstructed in finite element simulations and DAI and ASDH injury results were compared. We used the head/brain models of Total HUman Model for Safety (THUMS) and Global Human Body Models Consortium (GHBMC), which had been validated under several loading conditions. Four modeling approaches of the skull-brain interface in the head/brain models were evaluated. They were the original models from THUMS and GHBMC, the THUMS model with skull-brain interface changed to sliding contact, and the THUMS model with increased shear modulus of cerebrum, respectively. The results have shown that the definition of skull-brain interface would significantly influence the magnitude and distribution of the load transmitted to the brain. With sliding brain-skull interface, the brain had lower maximum principal stress compared to that with strong connected interface, while the maximum principal strain slightly increased. In addition, greater shear modulus resulted in slightly higher the maximum principal stress and significantly lower the maximum principal strain. This study has revealed that using models with different modeling approaches, the same value of injury metric may correspond to different injury severity.

Multi-tract multi-symptom relationships in pediatric concussion

Background

The heterogeneity of white matter damage and symptoms in concussion has been identified as a major obstacle to therapeutic innovation. In contrast, most diffusion MRI (dMRI) studies on concussion have traditionally relied on group-comparison approaches that average out heterogeneity. To leverage, rather than average out, concussion heterogeneity, we combined dMRI and multivariate statistics to characterize multi-tract multi-symptom relationships.

Methods

Using cross-sectional data from 306 previously concussed children aged 9-10 from the Adolescent Brain Cognitive Development Study, we built connectomes weighted by classical and emerging diffusion measures. These measures were combined into two informative indices, the first representing microstructural complexity, the second representing axonal density. We deployed pattern-learning algorithms to jointly decompose these connectivity features and 19 symptom measures.

Results

Early multi-tract multi-symptom pairs explained the most covariance and represented broad symptom categories, such as a general problems pair, or a pair representing all cognitive symptoms, and implicated more distributed networks of white matter tracts. Further pairs represented more specific symptom combinations, such as a pair representing attention problems exclusively, and were associated with more localized white matter abnormalities. Symptom representation was not systematically related to tract representation across pairs. Sleep problems were implicated across most pairs, but were related to different connections across these pairs. Expression of multi-tract features was not driven by sociodemographic and injury-related variables, as well as by clinical subgroups defined by the presence of ADHD. Analyses performed on a replication dataset showed consistent results.

Conclusions

Using a double-multivariate approach, we identified clinically-informative, cross-demographic multi-tract multi-symptom relationships. These results suggest that rather than clear one-to-one symptom-connectivity disturbances, concussions may be characterized by subtypes of symptom/connectivity relationships. The symptom/connectivity relationships identified in multi-tract multi-symptom pairs were not apparent in single-tract/single-symptom analyses. Future studies aiming to better understand connectivity/symptom relationships should take into account multi-tract multi-symptom heterogeneity.

Funding

Financial support for this work came from a Vanier Canada Graduate Scholarship from the Canadian Institutes of Health Research (G.I.G.), an Ontario Graduate Scholarship (S.S.), a Restracomp Research Fellowship provided by the Hospital for Sick Children (S.S.), an Institutional Research Chair in Neuroinformatics (M.D.), as well as a Natural Sciences and Engineering Research Council CREATE grant (M.D.).

Head Impact Kinematics and Brain Deformation in Paired Opposing Youth Football Players

Head impact exposure is often quantified using peak resultant kinematics. While kinematics describes the inertial response of the brain to impact, they do not fully capture the dynamic brain response. Strain, a measure of the tissue-level response of the brain, may be a better predictor of injury. In this study, kinematic and strain metrics were compared to contact characteristics in youth football. Players on 2 opposing teams were instrumented with head impact sensors to record impact kinematics. Video was collected to identify contact scenarios involving opposing instrumented players (ie, paired contact scenarios) and code contact characteristics (eg, player role, impact location). A previously validated, high-resolution brain finite element model, the atlas-based brain model, was used to simulate head impacts and calculate strain metrics. Fifty-two paired contact scenarios (n = 105 impacts) were evaluated. Lighter players tended to have greater biomechanical metrics compared to heavier players. Impacts to the top of the helmet were associated with lower strain metrics. Overall, strain was better correlated with rotational kinematics, suggesting these metrics may be better predictors of the tissue-level brain response than linear kinematics. Understanding the effect of contact characteristics on brain strain will inform future efforts to improve sport safety.

Mechanical characterisation of the human dura mater, falx cerebri and superior sagittal sinus

The cranial meninges have been shown to play a pivotal role in traumatic brain injury mechanopathology. However, while the mechanical response of the brain and its many subregions have been studied extensively, the meninges have conventionally been overlooked. This paper presents the first comparative mechanical analysis of human dura mater, falx cerebri and superior sagittal sinus tissues. Biaxial tensile analysis identified that these tissues are mechanically heterogeneous, in contrast to the assumption that the tissues are mechanically homogeneous which is typically employed in FE model design. A thickness of 0.91 ± 0.05 (standard error) mm for the falx cerebri was also identified. This data can aid in improving the biofidelity of the influential falx structure in FE models. Additionally, the use of a collagen hybridizing peptide on the superior sagittal sinus suggests this structure is particularly susceptible to the effects of circumferential stretch, which may have important implications for clinical treatment of dural venous sinus pathologies. Collectively, this research progresses understanding of meningeal mechanical and structural characteristics and may aid in elucidating the behaviour of these tissues in healthy and diseased conditions. STATEMENT OF SIGNIFICANCE: This study presents the first evaluation of human falx cerebri and superior sagittal sinus mechanical, geometrical and structural properties, along with a comparison to cranial dura mater. To mechanically characterise the tissues, biaxial tensile testing is conducted on the tissues. This analysis identifies, for the first time, mechanical stiffness differences between these tissues. Additionally, geometrical analysis identifies that there are thickness differences between the tissues. The evaluation of human meningeal tissues allows for direct implementation of the novel data to finite element head injury models to enable improved biofidelity of these influential structures in traumatic brain injury simulations. This work also identifies that the superior sagittal sinus may be easily damaged during clinical angioplasty procedures, which may inform the treatment of dural sinus pathologies.

Toward a Comprehensive Delineation of White Matter Tract-Related Deformation

Finite element (FE) models of the human head are valuable instruments to explore the mechanobiological pathway from external loading, localized brain response, and resultant injury risks. The injury predictability of these models depends on the use of effective criteria as injury predictors. The FE-derived normal deformation along white matter (WM) fiber tracts (i.e., tract-oriented strain) recently has been suggested as an appropriate predictor for axonal injury. However, the tract-oriented strain only represents a partial depiction of the WM fiber tract deformation. A comprehensive delineation of tract-related deformation may improve the injury predictability of the FE head model by delivering new tract-related criteria as injury predictors. Thus, the present study performed a theoretical strain analysis to comprehensively characterize the WM fiber tract deformation by relating the strain tensor of the WM element to its embedded fiber tract. Three new tract-related strains with exact analytical solutions were proposed, measuring the normal deformation perpendicular to the fiber tracts (i.e., tract-perpendicular strain), and shear deformation along and perpendicular to the fiber tracts (i.e., axial-shear strain and lateral-shear strain, respectively). The injury predictability of these three newly proposed strain peaks along with the previously used tract-oriented strain peak and maximum principal strain (MPS) were evaluated by simulating 151 impacts with known outcome (concussion or non-concussion). The results preliminarily showed that four tract-related strain peaks exhibited superior performance than MPS in discriminating concussion and non-concussion cases. This study presents a comprehensive quantification of WM tract-related deformation and advocates the use of orientation-dependent strains as criteria for injury prediction, which may ultimately contribute to an advanced mechanobiological understanding and enhanced computational predictability of brain injury.

Predictive Factors of Kinematics in Traumatic Brain Injury from Head Impacts Based on Statistical Interpretation

Brain tissue deformation resulting from head impacts is primarily caused by rotation and can lead to traumatic brain injury. To quantify brain injury risk based on measurements of kinematics on the head, finite element (FE) models and various brain injury criteria based on different factors of these kinematics have been developed, but the contribution of different kinematic factors has not been comprehensively analyzed across different types of head impacts in a data-driven manner. To better design brain injury criteria, the predictive power of rotational kinematics factors, which are different in (1) the derivative order (angular velocity, angular acceleration, angular jerk), (2) the direction and (3) the power (e.g., square-rooted, squared, cubic) of the angular velocity, were analyzed based on different datasets including laboratory impacts, American football, mixed martial arts (MMA), NHTSA automobile crashworthiness tests and NASCAR crash events. Ordinary least squares regressions were built from kinematics factors to the 95% maximum principal strain (MPS95), and we compared zero-order correlation coefficients, structure coefficients, commonality analysis, and dominance analysis. The angular acceleration, the magnitude and the first power factors showed the highest predictive power for the majority of impacts including laboratory impacts, American football impacts, with few exceptions (angular velocity for MMA and NASCAR impacts). The predictive power of rotational kinematics about three directions (x: posterior-to-anterior, y: left-to-right, z: superior-to-inferior) of kinematics varied with different sports and types of head impacts.

Regional White Matter Diffusion Changes Associated with the Cumulative Tensile Strain and Strain Rate in Nonconcussed Youth Football Players

The purpose of this study is to assess the relationship between regional white matter diffusion imaging changes and finite element strain measures in nonconcussed youth football players. Pre- and post-season diffusion-weighted imaging was performed in 102 youth football subject-seasons, in which no concussions were diagnosed. The diffusion data were normalized to the IXI template. Percent change in fractional anisotropy (%ΔFA) images were generated. Using data from the head impact telemetry system, the cumulative maximum principal strain one times strain rate (CMPS1 × SR), a measure of the cumulative tensile brain strain and strain rate for one season, was calculated for each subject. Two linear regression analyses were performed to identify significant positive or inverse relationships between CMPS1 × SR and %ΔFA within the international consortium for brain mapping white matter mask. Age, body mass index, days between pre- and post-season imaging, previous brain injury, attention disorder diagnosis, and imaging protocol were included as covariates. False discovery rate correction was used with corrected alphas of 0.025 and voxel thresholds of zero. Controlling for all covariates, a significant, positive linear relationship between %ΔFA and CMPS1 × SR was identified in the bilateral cingulum, fornix, internal capsule, external capsule, corpus callosum, corona radiata, corticospinal tract, cerebral and middle cerebellar peduncle, superior longitudinal fasciculus, and right superior fronto-occipital fasciculus. Post hoc analyses further demonstrated significant %ΔFA differences between high-strain football subjects and noncollision control athletes, no significant %ΔFA differences between low-strain subjects and noncollision control athletes, and that CMPS1 × SR significantly explained more %ΔFA variance than number of head impacts alone.

Intracranial Displacement Measurements Within Targeted Anatomical Regions of a Postmortem Human Surrogate Brain Subjected to Impact

The dynamic response of the human brain subjected to impulsive loading conditions is of fundamental importance to the understanding of traumatic brain injuries. Due to the complexity of such measurements, the existing experimental datasets available to researchers are sparse. However, these measurements are used extensively in the validation of complex finite element models used in the design of protective equipment and the development of injury mitigation strategies. The primary objective of this study was to develop a comprehensive methodology to measure displacement in specific anatomical regions of the brain. A state-of-the-art high-speed cineradiography system was used to capture brain motion in post-mortem human surrogate specimens at a rate of 7500 fps. This paper describes the methodology used to capture these data and presents measurements from these tests. Two-dimensional displacement fields are presented and analyzed based on anatomical regions of the brain. These data demonstrated a multi-modal displacement response in several regions of the brain. The full response of the brain consisted of an elastic superposition of a series of bulk rotations of the brain about its centre of gravity. The displacement field could be linked directly to specific anatomical regions. The methods presented mark an improvement in temporal and spatial resolution of data collection, which has implications for our developing understanding of brain trauma.

Measurement of brain simulant strains in head surrogate under impact loading

Impact-induced traumatic brain injury (TBI) is a major source of disability and mortality. Knowledge of brain strains during impact (accelerative) loading is critical for the overall management of TBI, including the development of injury thresholds, personal protective equipment, and validation of computational models. Despite these needs, the current understanding of brain strains in humans or humanlike surrogates is limited, especially for injury causing loading magnitudes. Toward this end, we measured full-field, in-plane (2D) strains in a brain simulant using the hemispherical head surrogate. The hemispherical head was mounted on the Hybrid-III neck and subjected to impact loading using a linear impactor system. The resulting head kinematics was measured using a triaxial accelerometer and angular rate sensors. Dynamic, 2D strains in a brain simulant were obtained using high-speed imaging and digital image correlation. Concurrent finite element (FE) simulations of the experiment were also performed to gain additional insights. The role of stiff membranes of the head was also studied using experiments. Our results suggest that rotational modes dominate the response of the brain simulant. The wave propagation in the brain simulant as a result of impact has a timescale of ~100 ms. We obtain peak strains of ~20%, ~40%, ~60% for peak rotational accelerations of ~838, ~5170, ~11,860 rad/s², respectively. Further, peak strains in cortical regions are higher than subcortical regions by up to ~70%. The agreement between the experiments and FE simulations is reasonable in terms of spatiotemporal evolution of strain pattern and peak strain magnitudes. Experiments with the addition of falx and tentorium indicate significant strain concentration (up to 115%) in the brain simulant near the interface of falx or tentorium and brain simulant. Overall, this work provides important insights into the biomechanics of strain in the brain simulant during impact loading.

White Matter Disruption in Pediatric Traumatic Brain Injury: Results From ENIGMA Pediatric Moderate to Severe Traumatic Brain Injury

OBJECTIVE

Our study addressed aims (1) to test the hypothesis that moderate-severe traumatic brain injury (TBI) in pediatric patients is associated with widespread white matter (WM) disruption, (2) to test the hypothesis that age and sex affect WM organization after injury, and (3) to examine associations between WM organization and neurobehavioral outcomes.

METHODS

Data from 10 previously enrolled, existing cohorts recruited from local hospitals and clinics were shared with the Enhancing NeuroImaging Genetics Through Meta-Analysis (ENIGMA) Pediatric Moderate/Severe TBI (msTBI) working group. We conducted a coordinated analysis of diffusion MRI (dMRI) data using the ENIGMA dMRI processing pipeline.

RESULTS

Five hundred seven children and adolescents (244 with complicated msTBI and 263 controls) were included. Patients were clustered into 3 postinjury intervals: acute/subacute, <2 months; postacute, 2 to 6 months; and chronic, ≥6 months. Outcomes were dMRI metrics and postinjury behavioral problems as indexed by the Child Behavior Checklist. Our analyses revealed altered WM diffusion metrics across multiple tracts and all postinjury intervals (effect sizes range d = -0.5 to -1.3). Injury severity is a significant contributor to the extent of WM alterations but explained less variance in dMRI measures with increasing time after injury. We observed a sex-by-group interaction: female patients with TBI had significantly lower fractional anisotropy in the uncinate fasciculus than controls (β = 0.043), which coincided with more parent-reported behavioral problems (β = -0.0027).

CONCLUSIONS

WM disruption after msTBI is widespread, persistent, and influenced by demographic and clinical variables. Future work will test techniques for harmonizing neurocognitive data, enabling more advanced analyses to identify symptom clusters and clinically meaningful patient subtypes.

A new open-access platform for measuring and sharing mTBI data

Despite numerous research efforts, the precise mechanisms of concussion have yet to be fully uncovered. Clinical studies on high-risk populations, such as contact sports athletes, have become more common and give insight on the link between impact severity and brain injury risk through the use of wearable sensors and neurological testing. However, as the number of institutions operating these studies grows, there is a growing need for a platform to share these data to facilitate our understanding of concussion mechanisms and aid in the development of suitable diagnostic tools. To that end, this paper puts forth two contributions: (1) a centralized, open-access platform for storing and sharing head impact data, in collaboration with the Federal Interagency Traumatic Brain Injury Research informatics system (FITBIR), and (2) a deep learning impact detection algorithm (MiGNet) to differentiate between true head impacts and false positives for the previously biomechanically validated instrumented mouthguard sensor (MiG2.0), all of which easily interfaces with FITBIR. We report 96% accuracy using MiGNet, based on a neural network model, improving on previous work based on Support Vector Machines achieving 91% accuracy, on an out of sample dataset of high school and collegiate football head impacts. The integrated MiG2.0 and FITBIR system serve as a collaborative research tool to be disseminated across multiple institutions towards creating a standardized dataset for furthering the knowledge of concussion biomechanics.

Challenges and opportunities for neuroimaging in young patients with traumatic brain injury: a coordinated effort towards advancing discovery from the ENIGMA pediatric moderate/severe TBI group

Traumatic brain injury (TBI) is a major cause of death and disability in children in both developed and developing nations. Children and adolescents suffer from TBI at a higher rate than the general population, and specific developmental issues require a unique context since findings from adult research do not necessarily directly translate to children. Findings in pediatric cohorts tend to lag behind those in adult samples. This may be due, in part, both to the smaller number of investigators engaged in research with this population and may also be related to changes in safety laws and clinical practice that have altered length of hospital stays, treatment, and access to this population. The ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) Pediatric Moderate/Severe TBI (msTBI) group aims to advance research in this area through global collaborative meta-analysis of neuroimaging data. In this paper, we discuss important challenges in pediatric TBI research and opportunities that we believe the ENIGMA Pediatric msTBI group can provide to address them. With the paucity of research studies examining neuroimaging biomarkers in pediatric patients with TBI and the challenges of recruiting large numbers of participants, collaborating to improve statistical power and to address technical challenges like lesions will significantly advance the field. We conclude with recommendations for future research in this field of study.

Molecular and Diffusion Tensor Imaging Biomarkers of Traumatic Brain Injury: Principles for Investigation and Integration

The last 20 years have seen the advent of new technologies that enhance the diagnosis and prognosis of traumatic brain injury (TBI). There is recognition that TBI affects the brain beyond initial injury, in some cases inciting a progressive neuropathology that leads to chronic impairments. Medical researchers are now searching for biomarkers to detect and monitor this condition. Perhaps the most promising developments are in the biomolecular and neuroimaging domains. Molecular assays can identify proteins indicative of neuronal injury and/or degeneration. Diffusion imaging now allows sensitive evaluations of the brain's cellular microstructure. As the pace of discovery accelerates, it is important to survey the research landscape and identify promising avenues of investigation. In this review, we discuss the potential of molecular and diffusion tensor imaging (DTI) biomarkers in TBI research. Integration of these technologies could advance models of disease prognosis, ultimately improving care. To date, however, few studies have explored relationships between molecular and DTI variables in patients with TBI. Here, we provide a short primer on each technology, review the latest research, and discuss how these biomarkers may be incorporated in future studies.

Surface reconstruction from routine CT-scan shows large anatomical variations of falx cerebri and tentorium cerebelli

BACKGROUND

Finite element modeling of the human head offers an alternative to experimental methods in understanding the biomechanical response of the head in trauma brain injuries. Falx, tentorium, and their notches are important structures surrounding the brain, and data about their anatomical variations are sparse.

OBJECTIVE

To describe and quantify anatomical variations of falx cerebri, tentorium cerebelli, and their notches.

METHODS

3D reconstruction of falx and tentorium was performed by points identification on 40 brain CT-scans in a tailored Matlab program. A scatter plot was obtained for each subject, and 8 anatomical landmarks were selected. A reference frame was defined to determine the coordinates of landmarks. Segments and areas were computed. A reproducibility study was done.

RESULTS

The height of falx was 34.9 ± 3.9 mm and its surface area 56.5 ± 7.7 cm². The width of tentorium was 99.64 ± 4.79 mm and its surface area 57.6 ± 5.8 cm². The mean length, height, and surface area of falx notch were respectively 96.9 ± 8 mm, 41.8 ± 5.9 mm, and 28.8 ± 5.8 cm² (range 15.8-40.5 cm²). The anterior and maximal widths of tentorial notch were 25.5 ± 3.5 mm and 30.9 ± 2.5 mm; its length 54.9 ± 5.2 mm and its surface area 13.26 ± 1.6 cm². The length of falx notch correlated with the length of tentorial notch (r = 0.62, P < 0.05).

CONCLUSION

We observe large anatomical variations of falx, tentorium, and notches, crucial to better understand the biomechanics of brain injury, in personalized finite element models.

Analysis of HIC and Hydrostatic Pressure in the Human Head during NOCSAE Tests of American Football Helmets

Brain damage is a serious economic and social burden. Contact sports such as American football, are one of the most common sources of concussions. The biomechanical response of the head-helmet system caused by dynamic loading plays a major role. The literature has focused on measuring the resultant kinematics that act on the head and helmet during tackles. However, few studies have focused on helmet validation tests, supported by recent findings and emerging numerical approaches. The future of helmet standards could benefit from insights at the level of injury mechanisms, using numerical tools to assess the helmets. Therefore, in this work, a numerical approach is employed to investigate the influence of intracranial pressure (ICP) on brain pathophysiology during and after helmeted impacts, which are common in American football. The helmeted impacts were performed at several impact locations according to the NOCSAE standard (configurations A, AP, B, C, D, F, R, UT). In order to evaluate the ICP levels, the αHEAD finite element head and brain model was combined with a Hybrid III-neck structure and then coupled with an American football helmet to simulate the NOCSAE impacts. In addition, the ICP level was analyzed together with the resulting HIC value, since the latter is commonly used, in this application and others, as the injury criterion. The obtained results indicate that ICP values exceed the common threshold of head injury criteria and do not correlate with HIC values. Thus, this work raises concern about applying the HIC to predict brain injury in American football direct head impacts, since it does not correlate with ICP predicted with the FE head model.

White Matter Tract-Oriented Deformation Is Dependent on Real-Time Axonal Fiber Orientation

Traumatic axonal injury (TAI) is a critical public health issue with its pathogenesis remaining largely elusive. Finite element (FE) head models are promising tools to bridge the gap between mechanical insult, localized brain response, and resultant injury. In particular, the FE-derived deformation along the direction of white matter (WM) tracts (i.e., tract-oriented strain) has been shown to be an appropriate predictor for TAI. The evolution of fiber orientation in time during the impact and its potential influence on the tract-oriented strain remains unknown, however. To address this question, the present study leveraged an embedded element approach to track real-time fiber orientation during impacts. A new scheme to calculate the tract-oriented strain was proposed by projecting the strain tensors from pre-computed simulations along the temporal fiber direction instead of its static counterpart directly obtained from diffuse tensor imaging. The results revealed that incorporating the real-time fiber orientation not only altered the direction but also amplified the magnitude of the tract-oriented strain, resulting in a generally more extended distribution and a larger volume ratio of WM exposed to high deformation along fiber tracts. These effects were exacerbated with the impact severities characterized by the acceleration magnitudes. Results of this study provide insights into how best to incorporate fiber orientation in head injury models and derive the WM tract-oriented deformation from computational simulations, which is important for furthering our understanding of the underlying mechanisms of TAI.

Mechanical Properties of the Cranial Meninges: A Systematic Review

The meninges are membranous tissues that are pivotal in maintaining homeostasis of the central nervous system. Despite the importance of the cranial meninges in nervous system physiology and in head injury mechanics, our knowledge of the tissues' mechanical behavior and structural composition is limited. This systematic review analyzes the existing literature on the mechanical properties of the meningeal tissues. Publications were identified from a search of Scopus, Academic Search Complete, and Web of Science and screened for eligibility according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. The review details the wide range of testing techniques employed to date and the significant variability in the observed experimental findings. Our findings identify many gaps in the current literature that can serve as a guide for future work for meningeal mechanics investigators. The review identifies no peer-reviewed mechanical data on the falx and tentorium tissues, both of which have been identified as key structures in influencing brain injury mechanics. A dearth of mechanical data for the pia-arachnoid complex also was identified (no experimental mechanics studies on the human pia-arachnoid complex were identified), which is desirable for biofidelic modeling of human head injuries. Finally, this review provides recommendations on how experiments can be conducted to allow for standardization of test methodologies, enabling simplified comparisons and conclusions on meningeal mechanics.

Brain response of a computational head model for prescribed skull kinematics and simulated football helmet impact boundary conditions

Computational human body models (HBM) present a novel approach to predict brain response in football impact scenarios, with prescribed kinematic boundary conditions for the HBM skull typically used at present. However, computational optimization of helmets requires simulation of the coupled helmet and HBM model; which is much more complex and has not been assessed in the context of brain deformation and existing simplified approaches. In the current study, two boundary conditions and the resulting brain deformations were compared using a HBM head model: (1) a prescribed skull kinematics (PK) boundary condition using measured head kinematics from experimental impacts; and (2) a novel detailed simulation of a HBM head and neck, helmet and linear impactor (HBM-S). While lateral and rear impacts exhibited similar levels of maximum principal strain (MPS) in the brain tissue using both boundary conditions, differences were noted in the frontal orientation (at 9.3 m/s, MPS was 0.39 for PK, 0.54 for HBM-S). Importantly, both PK and HBM-S boundary conditions produced a similar distribution of MPS throughout the brain for each impact orientation considered. Within the corpus callosum and thalamus, high MPS was associated with lateral impacts and lower values with frontal and rear impacts. The good correspondence of both boundary conditions is encouraging for future optimization of helmet designs. A limitation of the PK approach is the need for experimental head kinematics data, while the HBM-S can predict brain response for varying impact conditions and helmet configurations, with potential as a tool to improve helmet protection performance.

Mechanical and structural characterisation of the dural venous sinuses

The dural venous sinuses play an integral role in draining venous blood from the cranial cavity. As a result of the sinuses anatomical location, they are of significant importance when evaluating the mechanopathology of traumatic brain injury (TBI). Despite the importance of the dural venous sinuses in normal neurophysiology, no mechanical analyses have been conducted on the tissues. In this study, we conduct mechanical and structural analysis on porcine dural venous sinus tissue to help elucidate the tissues’ function in healthy and diseased conditions. With longitudinal elastic moduli values ranging from 33 to 58 MPa, we demonstrate that the sinuses exhibit higher mechanical stiffness than that of native dural tissue, which may be of interest to the field of TBI modelling. Furthermore, by employing histological staining and a colour deconvolution protocol, we show that the sinuses have a collagen-dominant extracellular matrix, with collagen area fractions ranging from 84 to 94%, which likely explains the tissue’s large mechanical stiffness. In summary, we provide the first investigation of the dural venous sinus mechanical behaviour with accompanying structural analysis, which may aid in understanding TBI mechanopathology.

Finite element simulation of head impacts in mixed martial arts

Thirteen MMA athletes were fitted with the MiG2.0 Stanford instrumented mouthguard. 451 video confirmed impacts were recoded during sparring sessions and competitive events. The competitive events resulted in five concussions. The impact with the highest angular acceleration from each event was simulated using the GHBMC head model. Average strain in the corpus callosum of concussed fighters was 0.27, which was 87.9% higher than uninjured fighters and was the best strain indicator of concussion. The best overall predictor of concussion found in this study was shear stress in the corpus callosum which differed by 111.4% between concussed and uninjured athletes.

Mechanisms and variances of rotation-induced brain injury: a parametric investigation between head kinematics and brain strain

There lacks a comprehensive understanding of the correlation between head kinematics and brain strain especially deep-brain strain, partially resulting the deficiency of understanding brain injury mechanisms and the difficulty of choosing appropriate brain injury metrics. Hence, we simulated 76 impacts that were focused on concussion-relevant rotational kinematics and evaluated cumulative strain damage measure (CSDM) and average strain that could represent brain strain distribution. For the whole brain, axial rotation induced the highest CSDM, while lateral bending produced the lowest CSDM. However, for the deep-brain components, lateral bending produced the highest CSDM to the corpus callosum and thalamus. We further confirmed that brain strain was mainly produced by rotational kinematics, for which the effect of rotational deceleration could not be ignored with the deceleration influencing CSDM20 up to 27%. Our data supported that peak rotational velocity correlated to brain strain with an average R² of 0.77 across various impact directions and different shapes of loading curves. The correlation between peak rotational velocity and brain strain reached to an average R² of 0.99 for each specific impact direction. Our results supported using direction-specific peak rotation velocity for predicting strain-related brain injury. Additionally, we highlighted the importance of investigating whole-brain and deep-brain strain, as well as considering rotational deceleration.

Three-dimensional finite element analysis of the dural folds and the human skull under head acceleration

Bone and collagen fiber architecture adapt to external mechanical loads. In humans, due to the low insertion of the temporal muscle, mastication does not lead to a physiological loading of the calvaria. Forces applied to the skull by the dural folds can lead to compressive stresses in the calvaria. To investigate the relationship between mechanical loads and form in the skull and its membranes, in a finite element three-dimensional model of the human skull, loads due to head acceleration in daily activities are applied to the falx cerebri and the tentorium cerebelli. The dural folds are modeled as membranes. The stress paths in the dural folds correlate with anatomical fiber direction. Head accelerations of 9 g lead to compressive stress in the calvaria. Finite element analysis of the falx cerebri and the tentorium cerebelli can be used to study the influence of mechanical stresses on the ossification of the dural folds and their impact on calvarial growth. This study presents an example of functional loading of bone by fibrous membranes and describes a possible mechanism by which Wolff's law works on the bone of the calvaria creating evolutionarily beneficial lightweight constructions.

Incorporation of vasculature in a head injury model lowers local mechanical strains in dynamic impact

Cerebral vasculature is several orders of magnitude stiffer than the brain tissue. However, only a handful of studies have investigated its potential stiffening effect on dynamic brain strains; yet, they report contradictory findings. Here, we reanalyze the cerebrovascular stiffening effect by incorporating vasculature derived from the latest neuroimaging atlases into a re-meshed Worcester Head Injury Model using an embedded element method. Regional brain strains with and without vasculature were simulated using a reconstructed, predominantly sagittal head impact. Using the two previously adopted linear or non-linear vessel material models, we reproduced the earlier conflicting results (~40% vs. ~1-6% in regional strain reductions). Nevertheless, with refitted non-linear material models chosen to represent the average dynamic tension behaviors of arteries and veins, respectively, inclusion of vasculature reduced regional brain strains by ~13-36% relative to the baselines without vasculature. Compared to the whole brain baseline response, inclusion of vasculature led to an element-wise linear regression slope of 0.8 and a Pearson correlation coefficient of 0.8. The vascular stiffening effect appears mild for the whole brain but more significant locally, which should not be ignored in head injury models. Nevertheless, more work is necessary to investigate the cerebrovascular mechanical behaviors and loading environment to allow for more biofidelic modeling of the brain in the future.

A network-based response feature matrix as a brain injury metric

Conventional brain injury metrics are scalars that treat the whole head/brain as a single unit but do not characterize the distribution of brain responses. Here, we establish a network-based "response feature matrix" to characterize the magnitude and distribution of impact-induced brain strains. The network nodes and edges encode injury risks to the gray matter regions and their white matter interconnections, respectively. The utility of the metric is illustrated in injury prediction using three independent, real-world datasets: two reconstructed impact datasets from the National Football League (NFL) and Virginia Tech, respectively, and measured concussive and non-injury impacts from Stanford University. Injury predictions with leave-one-out cross-validation are conducted using the two reconstructed datasets separately, and then by combining all datasets into one. Using support vector machine, the network-based injury predictor consistently outperforms four baseline scalar metrics including peak maximum principal strain of the whole brain (MPS), peak linear/rotational acceleration, and peak rotational velocity across all five selected performance measures (e.g., maximized accuracy of 0.887 vs. 0.774 and 0.849 for MPS and rotational acceleration with corresponding positive predictive values of 0.938, 0.772, and 0.800, respectively, using the reconstructed NFL dataset). With sufficient training data, real-world injury prediction is similar to leave-one-out in-sample evaluation, suggesting the potential advantage of the network-based injury metric over conventional scalar metrics. The network-based response feature matrix significantly extends scalar metrics by sampling the brain strains more completely, which may serve as a useful framework potentially allowing for other applications such as characterizing injury patterns or facilitating targeted multi-scale modeling in the future.

Convolutional neural network for efficient estimation of regional brain strains

Head injury models are important tools to study concussion biomechanics but are impractical for real-world use because they are too slow. Here, we develop a convolutional neural network (CNN) to estimate regional brain strains instantly and accurately by conceptualizing head rotational velocity profiles as two-dimensional images for input. We use two impact datasets with augmentation to investigate the CNN prediction performances with a variety of training-testing configurations. Three strain measures are considered, including maximum principal strain (MPS) of the whole brain, MPS of the corpus callosum, and fiber strain of the corpus callosum. The CNN is further tested using an independent impact dataset (N = 314) measured in American football. Based on 2592 training samples, it achieves a testing R² of 0.916 and root mean squared error (RMSE) of 0.014 for MPS of the whole brain. Combining all impact-strain response data available (N = 3069), the CNN achieves an R² of 0.966 and RMSE of 0.013 in a 10-fold cross-validation. This technique may enable a clinical diagnostic capability to a sophisticated head injury model, such as facilitating head impact sensors in concussion detection via a mobile device. In addition, it may transform current acceleration-based injury studies into focusing on regional brain strains. The trained CNN is publicly available along with associated code and examples at https://github.com/Jilab-biomechanics/CNN-brain-strains. They will be updated as needed in the future.

Novel strain analysis informs about injury susceptibility of the corpus callosum to repeated impacts

Increasing evidence for the cumulative effects of head trauma on structural integrity of the brain has emphasized the need to understand the relationship between tissue mechanic properties and injury susceptibility. Here, diffusion tensor imaging, helmet accelerometers and amplified magnetic resonance imaging were combined to gather insight about the region-specific vulnerability of the corpus callosum to microstructural changes in white-matter integrity upon exposure to sub-concussive impacts. A total of 33 male Canadian football players (mean_age = 20.3 ± 1.4 years) were assessed at three time points during a football season (baseline pre-season, mid-season and post-season). The athletes were split into a LOW (N = 16) and HIGH (N = 17) exposure group based on the frequency of sub-concussive impacts sustained on a per-session basis, measured using the helmet-mounted accelerometers. Longitudinal decreases in fractional anisotropy were observed in anterior and posterior regions of the corpus callosum (average cluster size = 40.0 ± 4.4 voxels; P < 0.05, corrected) for athletes from the HIGH exposure group. These results suggest that the white-matter tract may be vulnerable to repetitive sub-concussive collisions sustained over the course of a football season. Using these findings as a basis for further investigation, a novel exploratory analysis of strain derived from sub-voxel motion of brain tissues in response to cardiac impulses was developed using amplified magnetic resonance imaging. This approach revealed specific differences in strain (and thus possibly stiffness) along the white-matter tract (P < 0.0001) suggesting a possible signature relationship between changes in white-matter integrity and tissue mechanical properties. In light of these findings, additional information about the viscoelastic behaviour of white-matter tissues may be imperative in elucidating the mechanisms responsible for region-specific differences in injury susceptibility observed, for instance, through changes in microstructural integrity following exposure to sub-concussive head impacts.

Vulnerable locations on the head to brain injury and implications for helmet design

In studying traumatic brain injury (TBI), it has been long hypothesized that the head is more vulnerable to injury from impacts in certain directions or locations, as the relationship between impact force and the resulting neurological outcome is complex and can vary significantly between individual cases. Many studies have identified head angular acceleration to be the putative cause of brain trauma, but it is not well understood how impact location can affect the resulting head kinematics and tissue strain. Here, we identify the susceptibility of the head to accelerations and brain strain from normal forces at contact points across the surface of the skull and jaw using a 3-dimensional, 20 degree-of-freedom rigid-body head and cervical spine model. We find that head angular acceleration and brain tissue strain resulting from an input force can vary by orders of magnitude based on impact location on the skull, with the mandible as the most vulnerable region. Conversely, head linear acceleration is not sensitive to contact location. Using these analyses, we present an optimization scheme to distribute helmet padding thickness to minimize angular acceleration, resulting in a reduction of angular acceleration by an estimated 25% at the most vulnerable contact point compared to uniform thickness padding. This work gives intuition behind the relationship between input force and resulting brain injury risk, and presents a framework for developing and evaluating novel head protection gear.