Predicting Concussion Outcome by Integrating Finite Element Modeling and Network Analysis

The biomechanical signature of loss of consciousness: computational modelling of elite athlete head injuries

Sports related head injuries can cause transient neurological events including loss of consciousness and dystonic posturing. However, it is unknown why head impacts that appear similar produce distinct neurological effects. The biomechanical effect of impacts can be estimated using computational models of strain within the brain. Here, we investigate the strain and strain rates produced by professional American football impacts that led to loss of consciousness, posturing or no neurological signs. We reviewed 1280 National Football League American football games and selected cases where the team's medical personnel made a diagnosis of concussion. Videos were then analysed for signs of neurological events. We identified 20 head impacts that showed clear video signs of loss of consciousness and 21 showing clear abnormal posturing. Forty-one control impacts were selected where there was no observable evidence of neurological signs, resulting in 82 videos of impacts for analysis. Video analysis was used to guide physical reconstructions of these impacts, allowing us to estimate the impact kinematics. These were then used as input to a detailed 3D high-fidelity finite element model of brain injury biomechanics to estimate strain and strain rate within the brain. We tested the hypotheses that impacts producing loss of consciousness would be associated with the highest biomechanical forces, that loss of consciousness would be associated with high forces in brainstem nuclei involved in arousal and that dystonic posturing would be associated with high forces in motor regions. Impacts leading to loss of consciousness compared to controls produced higher head acceleration (linear acceleration; 81.5 g ± 39.8 versus 47.9 ± 21.4; P = 0.004, rotational acceleration; 5.9 krad/s2 ± 2.4 versus 3.5 ± 1.6; P < 0.001) and in voxel-wise analysis produced larger brain deformation in many brain regions, including parts of the brainstem and cerebellum. Dystonic posturing was also associated with higher deformation compared to controls, with brain deformation observed in cortical regions that included the motor cortex. Loss of consciousness was specifically associated with higher strain rates in brainstem regions implicated in maintenance of consciousness, including following correction for the overall severity of impact. These included brainstem nuclei including the locus coeruleus, dorsal raphé and parabrachial complex. The results show that in head impacts producing loss of consciousness, brain deformation is disproportionately seen in brainstem regions containing nuclei involved in arousal, suggesting that head impacts produce loss of consciousness through a biomechanical effect on key brainstem nuclei involved in the maintenance of consciousness.

Approximating subject-specific brain injury models via scaling based on head-brain morphological relationships

Most human head/brain models represent a generic adult male head/brain. They may suffer in accuracy when investigating traumatic brain injury (TBI) on a subject-specific basis. Subject-specific models can be developed from neuroimages; however, neuroimages are not typically available in practice. In this study, we establish simple and elegant regression models between brain outer surface morphology and head dimensions measured from neuroimages along with age and sex information (N = 191; 141 males and 50 females with age ranging 14-25 years). The regression models are then used to approximate subject-specific brain models by scaling a generic counterpart, without using neuroimages. Model geometrical accuracy is assessed using adjusted [Formula: see text] and absolute percentage error (e.g., 0.720 and 3.09 ± 2.38%, respectively, for brain volume when incorporating tragion-to-top). For a subset of 11 subjects (from smallest to largest in brain volume), impact-induced brain strains are compared with those from "morphed models" derived from neuroimage-based mesh warping. We find that regional peak strains from the scaled subject-specific models are comparable to those of the morphed counterparts but could be considerably different from those of the generic model (e.g., linear regression slope of 1.01-1.03 for gray and white matter regions versus 1.16-1.19, or up to ~ 20% overestimation for the smallest brain studied). These results highlight the importance of incorporating brain morphological variations in impact simulation and demonstrate the feasibility of approximating subject-specific brain models without neuroimages using age, sex, and easily measurable head dimensions. The scaled models may improve subject specificity for future TBI investigations.

Biomechanics of Traumatic Head and Neck Injuries on Women: A State-of-the-Art Review and Future Directions

The biomechanics of traumatic injuries of the human body as a consequence of road crashes, falling, contact sports, and military environments have been studied for decades. In particular, traumatic brain injury (TBI), the so-called "silent epidemic", is the traumatic insult responsible for the greatest percentage of death and disability, justifying the relevance of this research topic. Despite its great importance, only recently have research groups started to seriously consider the sex differences regarding the morphology and physiology of women, which differs from men and may result in a specific outcome for a given traumatic event. This work aims to provide a summary of the contributions given in this field so far, from clinical reports to numerical models, covering not only the direct injuries from inertial loading scenarios but also the role sex plays in the conditions that precede an accident, and post-traumatic events, with an emphasis on neuroendocrine dysfunctions and chronic traumatic encephalopathy. A review on finite element head models and finite element neck models for the study of specific traumatic events is also performed, discussing whether sex was a factor in validating them. Based on the information collected, improvement perspectives and future directions are discussed.

Use of Brain Biomechanical Models for Monitoring Impact Exposure in Contact Sports

Head acceleration measurement sensors are now widely deployed in the field to monitor head kinematic exposure in contact sports. The wealth of impact kinematics data provides valuable, yet challenging, opportunities to study the biomechanical basis of mild traumatic brain injury (mTBI) and subconcussive kinematic exposure. Head impact kinematics are translated into brain mechanical responses through physics-based computational simulations using validated brain models to study the mechanisms of injury. First, this article reviews representative legacy and contemporary brain biomechanical models primarily used for blunt impact simulation. Then, it summarizes perspectives regarding the development and validation of these models, and discusses how simulation results can be interpreted to facilitate injury risk assessment and head acceleration exposure monitoring in the context of contact sports. Recommendations and consensus statements are presented on the use of validated brain models in conjunction with kinematic sensor data to understand the biomechanics of mTBI and subconcussion. Mainly, there is general consensus that validated brain models have strong potential to improve injury prediction and interpretation of subconcussive kinematic exposure over global head kinematics alone. Nevertheless, a major roadblock to this capability is the lack of sufficient data encompassing different sports, sex, age and other factors. The authors recommend further integration of sensor data and simulations with modern data science techniques to generate large datasets of exposures and predicted brain responses along with associated clinical findings. These efforts are anticipated to help better understand the biomechanical basis of mTBI and improve the effectiveness in monitoring kinematic exposure in contact sports for risk and injury mitigation purposes.

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.

Brain architecture-based vulnerability to traumatic injury

The white matter tracts forming the intricate wiring of the brain are subject-specific; this heterogeneity can complicate studies of brain function and disease. Here we collapse tractography data from the Human Connectome Project (HCP) into structural connectivity (SC) matrices and identify groups of similarly wired brains from both sexes. To characterize the significance of these architectural groupings, we examined how similarly wired brains led to distinct groupings of neural activity dynamics estimated with Kuramoto oscillator models (KMs). We then lesioned our networks to simulate traumatic brain injury (TBI) and finally we tested whether these distinct architecture groups’ dynamics exhibited differing responses to simulated TBI. At each of these levels we found that brain structure, simulated dynamics, and injury susceptibility were all related to brain grouping. We found four primary brain architecture groupings (two male and two female), with similar architectures appearing across both sexes. Among these groupings of brain structure, two architecture types were significantly more vulnerable than the remaining two architecture types to lesions. These groups suggest that mesoscale brain architecture types exist, and these architectural differences may contribute to differential risks to TBI and clinical outcomes across the population.

American Football Helmet Effectiveness Against a Strain-Based Concussion Mechanism

Brain strain is increasingly being used in helmet design and safety performance evaluation as it is generally considered as the primary mechanism of concussion. In this study, we investigate whether different helmet designs can meaningfully alter brain strains using two commonly used metrics, peak maximum principal strain (MPS) of the whole brain and cumulative strain damage measure (CSDM). A convolutional neural network (CNN) that instantly produces detailed brain strains is first tested for accuracy for helmeted head impacts. Based on N = 144 impacts in 12 impact conditions from three random and representative helmet models, we conclude that the CNN is sufficiently accurate for helmet testing applications, for elementwise MPS (success rate of 98.6%), whole-brain peak MPS and CSDM (coefficient of determination of 0.977 and 0.980, with root mean squared error of 0.015 and 0.029, respectively). We then apply the technique to 23 football helmet models (N = 1104 impacts) to reproduce elementwise MPS. Assuming a concussion would occur when peak MPS or CSDM exceeds a threshold, we sweep their thresholds across the value ranges to evaluate the number of predicted hypothetical concussions that different helmets sustain across the impact conditions. Relative to the 12 impact conditions tested, we find that the "best" and "worst" helmets differ by an average of 22.5% in terms of predicted concussions, ranging from 0 to 42% (the latter achieved at the threshold value of 0.28 for peak MPS and 0.4 for CSDM, respectively). Such a large variation among helmets in strain-based concussion predictions demonstrate that helmet designs can still be optimized in a clinically meaningful way. The robustness and accuracy of the CNN tool also suggest its potential for routine use for helmet design and safety performance evaluation in the future. The CNN is freely available online at https://github.com/Jilab-biomechanics/CNN-brain-strains .

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.).

Individual-specific networks for prediction modelling – A scoping review of methods

Background

Recent advances in biotechnology enable the acquisition of high-dimensional data on individuals, posing challenges for prediction models which traditionally use covariates such as clinical patient characteristics. Alternative forms of covariate representations for the features derived from these modern data modalities should be considered that can utilize their intrinsic interconnection. The connectivity information between these features can be represented as an individual-specific network defined by a set of nodes and edges, the strength of which can vary from individual to individual. Global or local graph-theoretical features describing the network may constitute potential prognostic biomarkers instead of or in addition to traditional covariates and may replace the often unsuccessful search for individual biomarkers in a high-dimensional predictor space.

Methods

We conducted a scoping review to identify, collate and critically appraise the state-of-art in the use of individual-specific networks for prediction modelling in medicine and applied health research, published during 2000–2020 in the electronic databases PubMed, Scopus and Embase.

Results

Our scoping review revealed the main application areas namely neurology and pathopsychology, followed by cancer research, cardiology and pathology (N = 148). Network construction was mainly based on Pearson correlation coefficients of repeated measurements, but also alternative approaches (e.g. partial correlation, visibility graphs) were found. For covariates measured only once per individual, network construction was mostly based on quantifying an individual’s contribution to the overall group-level structure. Despite the multitude of identified methodological approaches for individual-specific network inference, the number of studies that were intended to enable the prediction of clinical outcomes for future individuals was quite limited, and most of the models served as proof of concept that network characteristics can in principle be useful for prediction.

Conclusion

The current body of research clearly demonstrates the value of individual-specific network analysis for prediction modelling, but it has not yet been considered as a general tool outside the current areas of application. More methodological research is still needed on well-founded strategies for network inference, especially on adequate network sparsification and outcome-guided graph-theoretical feature extraction and selection, and on how networks can be exploited efficiently for prediction modelling.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12874-022-01544-6.

An interdisciplinary computational model for predicting traumatic brain injury: Linking biomechanics and functional neural networks

The brain is a complex network consisting of neuron cell bodies in the gray matter and their axonal projections, forming the white matter tracts. These neurons are supported by an equally complex vascular network as well as glial cells. Traumatic brain injury (TBI) can lead to the disruption of the structural and functional brain networks due to disruption of both neuronal cell bodies in the gray matter as well as their projections and supporting cells. To explore how an impact can alter the function of brain networks, we integrated a finite element (FE) brain mechanics model with linked models of brain dynamics (Kuramoto oscillator) and vascular perfusion (Balloon-Windkessel) in this study. We used empirical resting-state functional magnetic resonance imaging (MRI) data to optimize the fit of our brain dynamics and perfusion models to clinical data. Results from the FE model were used to mimic injury in these optimized brain dynamics models: injury to the nodes (gray matter) led to a decrease in the nodal oscillation frequency, while damage to the edges (axonal connections/white matter) progressively decreased coupling among connected nodes. A total of 53 cases, including 33 non-injurious and 20 concussive head impacts experienced by professional American football players were simulated using this integrated model. We examined the correlation of injury outcomes with global measures of structural connectivity, neural dynamics, and functional connectivity of the brain networks when using different lesion methods. Results show that injurious head impacts cause significant alterations in global network topology regardless of lesion methods. Changes between the disrupted and healthy functional connectivity (measured by Pearson correlation) consistently correlated well with injury outcomes (AUC≥0.75), although the predictive performance is not significantly different (p>0.05) to that of traditional kinematic measures (angular acceleration). Intriguingly, our lesion model for gray matter damage predicted increases in global efficiency and clustering coefficient with increases in injury risk, while disrupting axonal connections led to lower network efficiency and clustering. When both injury mechanisms were combined into a single injury prediction model, the injury prediction performance depended on the thresholds used to determine neurodegeneration and mechanical tolerance for axonal injury. Together, these results point towards complex effects of mechanical trauma to the brain and provide a new framework for understanding brain injury at a causal mechanistic level and developing more effective diagnostic methods and therapeutic interventions.

A Machine Learning Approach to Investigate the Uncertainty of Tissue-Level Injury Metrics for Cerebral Contusion

Controlled cortical impact (CCI) on porcine brain is often utilized to investigate the pathophysiology and functional outcome of focal traumatic brain injury (TBI), such as cerebral contusion (CC). Using a finite element (FE) model of the porcine brain, the localized brain strain and strain rate resulting from CCI can be computed and compared to the experimentally assessed cortical lesion. This way, tissue-level injury metrics and corresponding thresholds specific for CC can be established. However, the variability and uncertainty associated with the CCI experimental parameters contribute to the uncertainty of the provoked cortical lesion and, in turn, of the predicted injury metrics. Uncertainty quantification via probabilistic methods (Monte Carlo simulation, MCS) requires a large number of FE simulations, which results in a time-consuming process. Following the recent success of machine learning (ML) in TBI biomechanical modeling, we developed an artificial neural network as surrogate of the FE porcine brain model to predict the brain strain and the strain rate in a computationally efficient way. We assessed the effect of several experimental and modeling parameters on four FE-derived CC injury metrics (maximum principal strain, maximum principal strain rate, product of maximum principal strain and strain rate, and maximum shear strain). Next, we compared the in silico brain mechanical response with cortical damage data from in vivo CCI experiments on pig brains to evaluate the predictive performance of the CC injury metrics. Our ML surrogate was capable of rapidly predicting the outcome of the FE porcine brain undergoing CCI. The now computationally efficient MCS showed that depth and velocity of indentation were the most influential parameters for the strain and the strain rate-based injury metrics, respectively. The sensitivity analysis and comparison with the cortical damage experimental data indicate a better performance of maximum principal strain and maximum shear strain as tissue-level injury metrics for CC. These results provide guidelines to optimize the design of CCI tests and bring new insights to the understanding of the mechanical response of brain tissue to focal traumatic brain injury. Our findings also highlight the potential of using ML for computationally efficient TBI biomechanics investigations.

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.

Influence of Anisotropic White Matter on Electroosmotic Flow Induced by Direct Current

Treatment of cerebral edema remains a major challenge in clinical practice and new innovative therapies are needed. This study presents a novel approach for mitigating cerebral edema by inducing bulk fluid transport utilizing the brain’s electroosmotic property using an anatomically detailed finite element head model incorporating anisotropy in the white matter (WM). Three representative anisotropic conductivity algorithms are employed for the WM and compared with isotropic WM. The key results are (1) the electroosmotic flow (EOF) is driven from the edema region to the subarachnoid space under an applied electric field with its magnitude linearly correlated to the electric field and direction following current flow pathways; (2) the extent of EOF distribution variation correlates highly with the degree of the anisotropic ratio of the WM regions; (3) the directions of the induced EOF in the anisotropic models deviate from its isotropically defined pathways and tend to move along the principal fiber direction. The results suggest WM anisotropy should be incorporated in head models for more reliable EOF evaluations for cerebral edema mitigation and demonstrate the promise of the electroosmosis based approach to be developed as a new therapy for edema treatment as evaluated with enhanced head models incorporating WM anisotropy.

Instantaneous Whole-Brain Strain Estimation in Dynamic Head Impact

Head injury models are notoriously time consuming and resource demanding in simulations, which prevents routine application. Here, we extend a convolutional neural network (CNN) to instantly estimate element-wise distribution of peak maximum principal strain (MPS) of the entire brain (>36 k speedup accomplished on a low-end computing platform). To achieve this, head impact rotational velocity and acceleration temporal profiles are combined into two-dimensional images to serve as CNN input for training and prediction of MPS. Compared with the directly simulated counterparts, the CNN-estimated responses (magnitude and distribution) are sufficiently accurate for 92.1% of the cases via 10-fold cross-validation using impacts drawn from the real world (n = 5661; range of peak rotational velocity in augmented data extended to 2-40 rad/sec). The success rate further improves to 97.1% for "in-range" impacts (n = 4298). When using the same CNN architecture to train (n = 3064) and test on an independent, reconstructed National Football League (NFL) impact dataset (n = 53; 20 concussions and 33 non-injuries), 51 out of 53, or 96.2% of the cases, are sufficiently accurate. The estimated responses also achieve virtually identical concussion prediction performances relative to the directly simulated counterparts, and they often outperform peak MPS of the whole brain (e.g., accuracy of 0.83 vs. 0.77 via leave-one-out cross-validation). These findings support the use of CNN for accurate and efficient estimation of spatially detailed brain strains across the vast majority of head impacts in contact sports. Our technique may hold the potential to transform traumatic brain injury (TBI) research and the design and testing standards of head protective gears by facilitating the transition from acceleration-based approximation to strain-based design and analysis. This would have broad implications in the TBI biomechanics field to accelerate new scientific discoveries. The pre-trained CNN is freely available online at https://github.com/Jilab-biomechanics/CNN-brain-strains.

A Review of Validation Methods for the Intracranial Response of FEHM to Blunt Impacts

The following is a review of the processes currently employed when validating the intracranial response of Finite Element Head Models (FEHM) against blunt impacts. The authors aim to collate existing validation tools, their applications and findings on their effectiveness to aid researchers in the validation of future FEHM and potential efforts in improving procedures. In this vain, publications providing experimental data on the intracranial pressure, relative brain displacement and brain strain responses to impacts in human subjects are surveyed and key data are summarised. This includes cases that have previously been used in FEHM validation and alternatives with similar potential uses. The processes employed to replicate impact conditions and the resulting head motion are reviewed, as are the analytical techniques used to judge the validity of the models. Finally, publications exploring the validation process and factors affecting it are critically discussed. Reviewing FEHM validation in this way highlights the lack of a single best practice, or an obvious solution to create one using the tools currently available. There is clear scope to improve the validation process of FEHM, and the data available to achieve this. By collecting information from existing publications, it is hoped this review can help guide such developments and provide a point of reference for researchers looking to validate or investigate FEHM in the future, enabling them to make informed choices about the simulation of impacts, how they are generated numerically and the factors considered during output assessment, whilst being aware of potential limitations in the process.