Abstract
Concussion is a public health crisis which results in a complex cascade of neurochemical changes in the brain that can have life changing consequences. Subconcussions are considered less serious and were overlooked until recently, but we now realise repetitive subconcussions, such as repetitive head impacts, can lead to serious neurological deficits. Subconcussions are common in contact sports, and the military where certain personnel are exposed to repetitive occupational blast overpressure. Postmortem studies in athletes reveal that cumulative duration of play and force from collisions are better predictors than concussion history for the presence and severity of chronic traumatic encephalopathy – a progressive and fatal neurodegenerative tauopathy, distinct from concussion, thought to be caused by repetitive head impacts, and only diagnosable postmortem – thus, an in vivo predictive biomarker would be game changing. Magnetoencephalography has exceptional temporal sampling for imaging the dynamics of neuronal electrochemical action, and functional MRI shows that functional connectivity is associated with tauopathy patterns. Therefore, both imaging modalities could provide a surrogate biomarker of tauopathy.
In this cross-sectional study, we examined the effects of repetitive subconcussion on neuronal activity and functional connectivity using magnetoencephalography and functional MRI, and on neurological symptoms and mental health in a military sample. For magnetoencephalography and outcome analyses, 81 participants were split into ‘high’ and ‘low’ blast exposure groups using the generalized blast exposure value: n=41 high blast (26.4–65.7 years; 4 females); n=40 low blast (28.0–63.3 years; 8 females). For fMRI, two high blast male participants without data were excluded: n=39 (29.6–65.7 years).
Magnetoencephalography revealed disrupted neuronal activity in participants with a greater history of repetitive subconcussions, including: neural slowing (excessive delta activity) in the right frontal and temporal lobes and subcortical regions (hippocampus, amygdala, caudate, pallidum, thalamus), and functional dysconnectivity in the posterior default mode network (low and high gamma). Notably, these irregularities were independent of concussion or traumatic stress history, and magnetoencephalography revealed functional dysconnectivity not detected with functional MRI. Besides regional slowing and functional disconnection in crucial brain hubs, those with greater blast exposure had poorer neurological outcomes in somatic and cognitive domains, with no blast-related differences in mental health and no associations between neurological symptoms and neuronal activity.
This study suggests that repetitive subconcussions have insidious effects on the brain and that magnetoencephalography provides a potential avenue for both treatment targets by identifying affected brain regions and in prevention by identifying those at risk of cumulative subconcussive neurotrauma.
1. Introduction
Concussion, or mild traumatic brain injury (mTBI), involves a complex cascade of neurochemical changes and neurological symptoms that can have life changing consequences.1 Repetitive subconcussions were thought to be less serious, and have been largely overlooked up until recently2–4 but are now considered insidious, and occur from repetitive head impacts (RHI) common in contact sports5,6 and in military settings with occupational repetitive blast overpressure (ReBOP) exposure (e.g., explosions, small arms, heavy weapons).7–11 Postmortem studies in athletes reveal that the accumulated duration of play5,6 and total force6 are better predictors than concussion history for the presence and severity of chronic traumatic encephalopathy (CTE), a progressive and fatal neurodegenerative tauopathy, which is distinct from concussions and persistent postconcussion symptoms (PPCS). Yet, the pathogenic mechanisms of CTE are not well understood,12 although histopathology suggests that expression of glial fibrillary acidic proteins at the cortical grey-white matter interface are a neuropathological hallmark of cumulative mild blast or impact neurotrauma.13 However, postmortem examinations have shown that CTE was not often detected in the brains of military members, as evidenced by modest neuropathologic changes, and CTE risk ratios were higher for decedents who experienced other traumatic brain injuries (TBI) outside of the military (e.g., contact sports) than those who only had blast exposure or military-related TBI.14
Given that CTE and proteinopathic dementias can only be diagnosed postmortem, an in vivo predictive biomarker would be game changing. Positron emission tomography (PET) radiotracer binding to tau with Flortaucipir offers a potential avenue,15 with uptake correlating with blast exposure in frontal, temporal, occipital, and cerebellar regions in warfighters and Veterans.11 However, a limitation of PET includes its high cost and concern of non-specific binding.16 A PET-tau study in Alzheimer’s dementia found that tauopathy colocalized with abnormal synaptic function quantified by magnetoencephalography (MEG),17 and another investigation in primary 4-repeat tauopathies showed that tauopathy patterns were associated with functional connectivity quantified by functional magnetic resonance imaging (fMRI).18 Together this research suggests that MEG and fMRI can quantify/identify dysfunction associated with tauopathy at a fraction of the cost of PET.
Electrophysiology with MEG offers a highly sensitive, non-invasive means to image the dynamics of neuronal electrochemical action by measuring magnetic fields generated by primary neuronal currents, allowing for micro- to macro-level examinations of neuronal circuits. Neuronal oscillations are important surrogate measures of cell physiology, indexing the rhythmic activity of brain excitability and inhibition, with abnormalities in the temporal pattern of oscillatory activity indicating pathology.19,20 MEG can measure neuronal oscillations because of its exceptional temporal sampling rate, whereas fMRI relies on the blood oxygen level dependent (BOLD) signal, a slower haemodynamic proxy for neural activity.
Pathological rhythmic neuronal activity – or ‘oscillopathies’ – and functional connectivity are important pathogenic markers of psychiatric and neurological disorders, including concussion,21–26 post-traumatic stress disorder (PTSD),21 Parkinson’s disease (PD),27 and Alzheimer’s disease (AD).28 “Pathological slowing” of rhythmic activity and changes in functional connectivity are common features of brain injury and neurodegeneration in military personnel and Veterans with concussion.24 In neurodegenerative disease specifically, neuronal slowing in PD has been observed in parieto-occipital cortices27 and neuronal slowing in AD predicts functional impairment and amyloid burden.28
Subconcussions have been studied in the context of RHI, but occupational blast exposure provides a model for the effects of cumulative, repetitive subconcussions common in certain military roles.29 Military personnel and Veterans exhibit neurobehavioural symptoms that can mimic persistent post-concussion symptoms, dementia, and PTSD. These include somatic, cognitive, and emotional symptoms and mental health issues (e.g., anxiety, depression).7–10 Also, higher lifetime subconcussive blast exposure is associated with poorer brain-related outcomes, including worse neurobehavioural symptoms (e.g., sleep, fatigue).30
Although there are functional neuroimaging studies of military personnel who have sustained an mTBI from blast exposure,31–34 there are no multimodal studies on the impact of repetitive subconcussions, or those dissociating the effects of concussion on brain function.35 There is one MEG study in combat Veterans which identified reduced functional connectivity in blast-related concussion with PTSD34; however, they did not disentangle repetitive subconcussive blast exposure from concussion. A magnetic resonance imaging (MRI) study showed that blast exposure resulted in worse cognitive, emotional, and somatic symptoms and smaller cortical volumes in the right superior frontal gyrus.10 An fMRI study revealed the effect that a single day of breacher training with acute, multiple blast exposures has on the brain, wherein the researchers gave half of the participants a jugular vein compression collar – those without the collars presented with blast-related impairment in working memory and auditory network connectivity, suggesting that the effects of repetitive subconcussion may be managed by simple protective devices.36 In sports, a longitudinal fMRI study showed that high school football players exhibited dysfunctional connectivity across a single season of play despite not showing any symptoms, suggesting pre-symptomatic brain changes.37 These studies highlight the value of functional neuroimaging to identify in vivo biomarkers of subconcussive neurotrauma before debilitating symptoms emerge.
The pathological mechanism of repetitive subconcussions – including blast – that might lead to CTE remains largely unknown, but repetitive blast exposure provides a model to investigate pathogenesis and novel markers.38 Such a model can inform basic science and provide translatable paths to refine training procedures and improve overall health outcomes for operational readiness.7–10,30,36,37 There is evidence that repetitive blast damages the central nervous system and/or the central or peripheral vestibular system due to blast waves entering the ear canal and subsequently causing mechanical reverberations throughout the brain.39,40 However, other studies argue that the mechanism extends beyond acceleration–deceleration injuries, due to the effect of blast waves on both air-filled organs and/or organs surrounded by fluid-filled cavities within the body.41,42
The purpose of this study was to understand the interaction between lifetime blast exposure estimates, neurobehavioural symptoms, mental health, and brain function outcomes in Canadian Armed Forces (CAF) and Royal Canadian Mounted Police (RCMP) personnel and Veterans exposed to the continuum of repetitive subconcussions from blast overpressure, controlling for concussion and traumatic stress history. The generalized blast exposure value (GBEV) tool29 was used to capture blast history – we performed a median split of the cohort into two groups of ‘high’ and ‘low’ blast exposure as described in this other protocol,43 to understand the impact and relationships between blast, health outcomes, and neural activity. Multimodal studies suggest that proteinopathic burden and neurodegeneration correlates with neuronal slowing and changes in brain communication17 – therefore, we predict that participants with higher blast load will show neuronal slowing (MEG) and abnormal functional connectivity (MEG, fMRI) alongside worse neurobehavioural and mental health symptoms.
2. Materials and methods
2.1 Participants
Participants were recruited as part of a large, multidimensional, cross-sectional, and longitudinal study with a transdiagnostic design to examine military brain and mental health. Methods of recruitment included dissemination of study information directly to military units, as well as health service channels and philanthropic organizations dealing with military and Veteran brain health (e.g., Project Enlist, Veterans Affairs Canada, and other Veteran’s groups). Participants were eligible for this study if they were active CAF/RCMP personnel or a Veteran of the CAF/RCMP at the time of study enrolment, had sufficient English language skills and cognitive capacity to participate in the study, and had no contradictions to MRI or MEG. Participant data collection occurred from September 2020 to May 2023, and along with all analyses, was conducted at the Hospital for Sick Children in Toronto, ON, Canada.
Initially, 84 total participants were recruited; 3 participants did not complete imaging (MEG or fMRI), and of the 81 remaining, 2 completed the MEG scan but not fMRI. For all MEG and outcome analyses, 81 CAF/RCMP personnel and Veterans from various backgrounds, including those from combat and emergency response roles involving weapons use, armoured, engineering, operators, explosive disposal, and artillery, were divided into two groups based on blast exposure levels assessed by the GBEV questionnaire: n = 41 high blast exposure group (aged 26.4 to 65.7 years; mean 48.1 ± 9.0 years; 4 females) and n = 40 low blast exposure group (aged 28.0 to 63.3 years; mean 47.8 ± 9.5 years; 8 females). For fMRI analyses, two male participants without fMRI data were removed from the high blast group described above so that: n = 39 high blast group (aged 29.6 to 65.7 years; mean 48.6 ± 8.5 years; 4 females). The GBEV is a tool for standardizing a lifetime of blast exposure.29 It quantifies units of blast exposure by querying participants about exposure over a lifetime across five categories of blast type/severity. This study was approved by the Hospital for Sick Children Ethics Board. Written informed consent was obtained from all participants in accordance with the Declaration of Helsinki.
2.2 Neurobehavioural and psychiatric outcomes
All participants completed neurobehavioural and mental health assessments which included: the Generalized Anxiety Disorder 7 screener (GAD7)44; Patient Health Questionnaire (PHQ9)45; PTSD Checklist Military Version (PCL-M)46; and the neurobehavioural symptom evaluation from the Sports Concussion Assessment Tool 5 (SCAT5).47 Participants also completed brief screeners related to diagnosed lifetime concussion history, using the Acute Concussion Evaluation (ACE)48 and defined by Department of Defence/Veterans Affairs guidelines for mTBI as injury resulting in loss of consciousness <30 minutes, post traumatic amnesia <24 hours, and a Glasgow Coma Score of 13 or more. Traumatic stress history was also estimated with the Brief Trauma Questionnaire (BTQ).49
2.3 Magnetoencephalography acquisition and analysis
Resting state MEG data were acquired on a 151 channel CTF system at the Hospital for Sick Children in Toronto. Participants were scanned in the supine position with fiducial coils attached at the nasion and bilateral preauricular pits for head motion recording. Data sampling rate was 600 Hz, bandpass filtered offline (high-pass, 1 Hz; low-pass, 150 Hz; notch, 60 Hz), and analyzed in the Fieldtrip Toolbox.50 Sensor level time series were mean centred. Cardiac and ocular artefacts were removed using independent component analysis (ICA). Time-series data were divided into 10 s epochs and the maximum number were retained from the 5 minute recording based on the following exclusion criteria: (a) head position did deviate more than 5 mm during a given epoch; (b) epochs that contained superconducting quantum interference device resets or exceeded a ±2 pT threshold following ICA component rejection.
MEG co-registration was performed using a single-shell head model per participant based on their anatomical T1-weighted MRI data normalized to Montreal Neurological Institute (MNI) space. A beamformer was used to resolve time series data from 90 regions of the Automated Anatomical Labelling (AAL) atlas51 with the parcel centroid of each AAL region used as the node location for beamformer reconstruction for which the Fieldtrip50 linearly constrained minimum variance vector beamformer52 was used, with 5% Tikhonov regularization. At each AAL node, lead fields were calculated from the template single shell head model for a unit current dipole across three dimensions. Beamformer weights of nodal neural activity were calculated through projection of the sensor weights along the dimensional axis with the highest singular value decomposition variance which produces a neural activity time series for each node.
The retained time series data from virtual electrodes were z-scored (i.e., mean centred, variance normalised). These broadband regional time courses per node were then filtered into delta (1-3 Hz), theta (3-7 Hz), alpha (8-14 Hz), beta (15-25 Hz), low gamma 1 (30-55 Hz), low gamma 2 (65-80 Hz), and high gamma (80-150 Hz) ranges. To reduce artificially inflated coupling that can result from beamformer leakage, a symmetric orthogonalization correction procedure was used.53 Welch’s Method was used to obtain the power spectrum density (PSD) at each AAL node. Lobe-wise average power was calculated per hemisphere for each participant by averaging across each AAL node within a given lobe.54 The Hilbert transform was applied to the band-limited time course to generate instantaneous estimates of the amplitude envelope, which was then down-sampled to 1 Hz.55 Pearson correlations were calculated between all AAL node pairs to quantify functional coupling. Amplitude envelope coupling (AEC) was chosen over other measures of neural communication (for example, phase synchronisation) as it is the most reliable measure of connectivity across sessions and over individuals,21,56,57 pointing to the greatest replicability as well as the lowest susceptibility to co-registration related errors.58
2.4 MRI acquisition and analysis
MRI data were acquired on a 3T Siemens PrismaFit scanner with a 20-channel head and neck coil at the Hospital for Sick Children in Toronto. For source localization of MEG and fMRI data, 3D magnetization-prepared rapid gradient echo T1-weighted images were acquired with 0.8 mm isotropic voxel size in 5:01 minutes (repetition time, TR, 1870 ms; echo time, TE, 3.1 ms; inversion time, TI, 945 ms; field-of-view, FOV, 240×256 mm with 240 0.8 mm slices).
Resting state fMRI data were acquired with 3 mm isotropic voxel size in 5 minutes (TR, 1500 ms; TE, 30 ms) and standard preprocessing was completed with the Data Processing Assistant for Resting-State fMRI toolbox59 and the SPM8 toolbox60 (http://www.fil.ion.ucl.ac.uk/spm). First, slice time correction was applied to the interleaved sequence acquisition, then realignment of the volumes was accomplished by calculating a six-parameter rigid body spatial transformation.61 None of the participants exhibited >3 degrees rotation or >3 mm displacement in any direction across the acquisition; therefore, no volumes were discarded. Next, normalization to MNI space was performed through unified T1-weighted image segmentation. The global signal was filtered out along with nuisance covariates including head motion parameters, white matter (WM) signal, and cerebrospinal fluid (CSF) signal to minimize the effects of non-neuronal oscillations and motion.62–65 Next, images were band-pass filtered at 0.01-0.1 Hz to retain low frequency oscillations in the resting state fMRI which are believed to represent neuronal activity,66 while suppressing low frequency drifts and high frequency noise resulting from other physiological activity.67,68 BOLD time domain signals were extracted from every voxel of the same 90 AAL atlas regions as MEG.51 For each individual AAL region, the mean timeseries were calculated across a parcel, then used to estimate functional connectivity with Pearson correlations.69
2.5 Statistical analysis
Statistical analyses were performed using MATLAB version R2022b, JASP version 0.18.0 (JASP, 2023), and the Network-Based Statistic (NBS) toolbox.70 False positives due to multiple comparisons were controlled using the Benjamini-Hochberg false discovery rate (FDR-BH) method at a significance threshold of p < 0.05. There are known effects of sex, age, concussion, and psychological trauma in these brain and outcome measures; thus, these confounders were controlled for as described below.
2.5.1 Participant demographics and outcomes
Analyses of covariance (ANCOVA) were used to test for blast-related group differences in outcome measures (GAD7, PHD9, PCL-M, SCAT5) while controlling for age, sex, traumatic stress history (BTQ), and number of diagnosed concussions (ACE). ANCOVA was also used to test for group differences in outcomes, stratified by concussion/no concussion history, while controlling for age and sex.
2.5.2 Power spectrum density analysis
ANCOVAs were used to test for blast-related group differences in power across all seven frequency bands across hemispheric lobes by calculating the average power for each AAL atlas node in the left and right frontal, temporal, parietal, and occipital lobes, and subcortical regions (grouping hippocampus, amygdala, caudate, pallidum, thalamus) while controlling for age, sex, psychological trauma, and number of concussions.
2.5.3 Network based analysis
The NBS toolbox was used to test group differences in functional connectivity within eight networks based on AAL atlas nodes (anterior DMN, aDMN; posterior DMN, pDMN; sensorimotor network, SMN; salience network, SN; memory network, MN; visual network, VN; central executive network, CEN; attention network, AN) in the seven MEG frequency bands and single fMRI connectomes, while controlling for age, sex, psychological trauma, and number of concussions (statistical model: threshold = 3.1; F-test; 5000 permutations; NBS method; Extent component size).
2.5.4 Concussion group analysis
Conveniently, of the 81 total participants, n = 38 (aged 26.4 to 65.7 years; mean 48.6 ± 9.5 years; 5 females) had a history of a clinical and confirmed diagnosis of at least one concussion, and the remaining n = 43 (aged 28.0 to 61.9 years; mean 47.3 ± 8.9 years; 7 females) did not; therefore, we also conducted an exploratory between groups analysis based on concussion history, covarying for age, sex, psychological trauma, and blast exposure as quantified by the GBEV, to examine concussion specific changes in brain function with MEG. For the fMRI analyses, the two male participants without fMRI data were removed: one from the concussion group (n = 37; aged 29.6 to 65.7 years; mean 49.2 ± 8.9 years; 5 females) and another from the non-concussion group (n = 42; aged 28.0 to 61.9 years; mean 47.3 ± 9.0 years; 7 females).
2.5.5 Brain function and outcome association analyses
As a post-hoc analysis, regressions were conducted between any brain measures (MEG power or connectivity, fMRI connectivity) and outcomes (mental health, neurological) that showed significant blast-level group differences, while controlling for age, sex, psychological trauma, and number of concussions.
3. Results
3.1 Demographics
The age of the high blast group (n = 41; range: 26.4 to 65.7 years; mean: 48.1 ± 9.0 years) did not differ from the low blast group (n = 40; range: 28.0 to 63.3 years; mean: 47.8 ± 9.5 years), nor did sex ratios (high blast group: 37 males, 4 females; low blast group: 32 males, 8 females) (Figure 1A), and there were no age or sex differences when the two high blast male participants were removed for fMRI analyses due to lack of data (n = 39; range: 29.6 to 65.7 years; mean 48.6 ± 8.5 years; 35 males, 4 females). Additionally, when re-stratifying the cohort by concussion history, the age of the concussion group (n = 38; range: 26.4 to 65.7 years; mean: 48.6 ± 9.5 years) did not differ from the non-concussion group (n = 43; range: 28.0 to 61.9 years; mean: 47.3 ± 8.9 years), nor did sex ratios (concussion group: 33 males, 5 females; non-concussion group: 36 males, 7 females), and there were no age or sex differences when the two male participants without fMRI data were removed, including one from the concussion group (n = 37; range: 29.6 to 65.7 years; mean: 49.2 ± 8.9 years; 32 males, 5 females) and another from the non-concussion group (n = 42; range: 28.0 to 61.9 years; mean: 47.3 ± 9.0 years; 35 males, 7 females) (Figure 1B).
3.2 Neurological symptoms due to repetitive subconcussive overpressure, not concussion history
The high blast exposure group showed significantly worse overall neurobehavioural symptom severity (Figure 2A), driven by poorer outcomes on somatic (Figure 2B) and cognitive subscales (Figure 2C), with no blast-related differences in emotional (Figure 2D) or sleep symptoms (Figure 2E) when controlling for age, sex, psychological trauma history, and number of diagnosed concussions.
There were no group differences between high and low blast in the number of diagnosed concussions (Figure 3A), psychological trauma severity (Figure 3B), and mental health outcomes, including depression (Figure 3C), anxiety (Figure 3D), or PTSD (Figure 3E) symptom severity. Moreover, there were no significant group differences for neurobehavioural, anxiety, depression, or PTSD symptoms in the concussion stratified analysis when controlling for age, sex, psychological trauma, and blast exposure history. These results suggest that the elevated neurological symptoms experienced by individuals with greater blast exposure is due to repetitive subconcussions and not concussion or traumatic stress history.
3.3 Neural slowing is evident in fronto-temporal and subcortical regions after high blast exposure
Comparisons of lobe-wise power spectrum reveal larger alpha amplitude in the high blast group, but otherwise comparable and expected 1/f power spectrum curves for both groups (Figure 4).
Delta activity (1-3 Hz) was significantly higher in the right frontal and right temporal lobes, as well as in right subcortical regions (hippocampus, amygdala, caudate, pallidum, thalamus) in the high blast compared to low blast group when controlling for age, sex, traumatic stress, and number of diagnosed concussions (Figure 5), but not for any other lobes at the delta frequency. There were no group differences for any other of the six frequency bands.
Moreover, there were no group differences when comparing diagnosed concussion versus no concussion history in regional power when controlling for age, sex, trauma, and blast exposure. Together, these results suggest that neural slowing in the right fronto-temporal lobes and subcortical regions is due to repetitive subconcussive neurotrauma and not concussion history.
3.4 Decreased pDMN connectivity in high blast exposure as measured by MEG but not fMRI
MEG revealed significant reductions in pDMN functional connectivity for the high blast group at low gamma (30-55 Hz; 7 nodes and 7 edges; Figure 6A) and high gamma (80-150 Hz; 14 nodes and 17 edges; Figure 6B).
For low gamma connectivity, the high blast group exhibited dysconnectivity in the pDMN including the right parahippocampus, left fusiform gyrus, bilateral superior parietal cortex, right inferior parietal cortex, right precuneus, and left middle temporal cortex (marked with asterisks, Figure 6A) – with no differences in the other seven networks (Supplemental Figures S1, S2). For high gamma connectivity, the high blast group exhibited dysconnectivity across pDMN areas including the bilateral posterior cingulate cortex (PCC), left hippocampus, left fusiform gyrus, bilateral superior and inferior parietal cortex, bilateral angular gyrus, bilateral precuneus, and bilateral middle temporal cortex (marked with asterisks, Figures 6B) – with no differences in the other seven networks (Supplemental Figures S3, S4).
The regions with the greatest difference (edges scaled by the F test statistic, Figures 6A & 6B) were mainly interhemispheric connections across both gamma ranges in the pDMN. There were no blast exposure group differences in any of the other individual networks at any frequency band, including the AN, CEN, aDMN, MN, SMN, SN, and VN. In contrast to the MEG results, fMRI did not reveal any significant blast-related differences (Supplemental Figures S5, S6). Additionally, there were no significant group differences in either MEG or fMRI connectivity when stratified by concussion history. These results suggest that (1) the dysconnectivity of the pDMN in the high blast group is due to repetitive subconcussions, and not concussion history; (2) MEG is more sensitive to network dysregulation from repetitive subconcussions than fMRI.
3.5 MEG delta power and functional dysconnectivity does not correlate with neurobehavioural outcomes
Despite group differences in delta activity (right frontal and temporal lobes and subcortical regions), pDMN functional connectivity, and neurological symptom outcomes (overall neurological symptom severity, and subscales of cognitive and somatic symptom severity), there were no significant associations between MEG delta power or functional connectivity and neurological symptom severity when controlling for age, sex, psychological trauma, and number of concussions (Figure 7A-I). For functional connectivity in the pDMN, node strength per participant was calculated for the node with the highest degree of connections with blast-related differences – namely, the right inferior parietal hub at low gamma and the left posterior cingulate hub at high gamma. There were also no significant associations between MEG functional connectivity and neurological symptom severity when controlling for age, sex, psychological trauma, and number of concussions (Figure 7J-O).
4. Discussion
4.1 Summary
We present the first multimodal functional neuroimaging evidence of dysregulated neurophysiological functioning in a sample of military members and Veterans with a history of repetitive subconcussions, resulting from varying degrees of exposure to blast overpressure – importantly, these results were observed independently of concussion or traumatic stress history. Along with regional slowing of activity, indexed by increased delta activity, and functional dysconnectivity in key brain hubs, those with higher blast exposure also exhibited worse neurological symptoms, especially in the somatic and cognitive domains, with no blast-related differences in mental health outcomes. Despite regional slowing and dysconnectivity and worse neurological symptoms in higher blast exposure, there were no significant associations between MEG power and network functional connectivity and neurological symptom severity.
4.2 Neural slowing in fronto-temporal and subcortical regions
Those with greater blast exposure exhibited neural slowing – elevated delta activity – in fronto-temporal and subcortical regions that was independent of concussion history. Moreover, we observed functional dysconnectivity in the pDMN – an effect that was absent in our fMRI measurements using the exact same cohort. Despite group blast-related differences in neurological outcomes and neuronal activity and functional connectivity, there were no associations between neurological outcomes and functional measures.
Animal models show a dose-response curve of blast exposure frequency with tauopathy71 that results in abnormal neuronal activity and neurobehavioural deficits.71–78 This work is consistent with our findings which show worse cognitive, somatic, and physical symptoms, neural slowing in the right fronto-temporal and subcortical regions, and lower functional connectivity in the pDMN in those with greater subconcussive load. Rodent studies have shown that even a single blast (i.e., mild-moderate shockwave) can have insidious effects, resulting in poorer memory function and elevated markers of neurodegeneration, microglia activation, and gliosis with decreased mature neurons in the prefrontal cortex and limbic areas, including the hippocampus, at three months post-exposure – suggesting these areas are particularly susceptible to the effects of blast.78
Another rodent study reported brain-wide and hippocampal CTE-linked neuropathology, including myelinated axonopathy, chronic neuroinflammation, and neurodegeneration two weeks following a single blast, and at one month, cognitive impairment, reduced synaptic plasticity, and slowed axonal conduction velocity72 – mechanistically, this would explain our observations of neural slowing. Importantly, the neurobiological mechanism of neuronal slowing differs between acute brain injury and neurodegenerative disease processes: slowing in brain injury is putatively caused by axonal shearing at the grey-white matter boundary,79 whereas in neurodegeneration, slowing is caused by proteinopathy, such as microtubule degradation, hyperphosphorylation of tau, and resultant neurofibrillary tangles (NFT).80,81 For example, in both civilian and military mTBI, multimodal neuroimaging with MEG detected neuronal slowing and diffusion tensor imaging (DTI) found proximal axonal / white matter damage.82
Our data suggest that delta oscillations, a stable in vivo biomarker of neural slowing, could have clinical utility in a broad range of neurodegenerative diseases, including Parkinson’s disease (PD),27 Alzheimer’s disease (AD),28 and other dementias such as fronto-temporal and vascular dementia.83 A study by Goldstein et al.72 reported CTE-linked cortical and hippocampal neuropathology; it included a complementary ex vivo analysis of brains from blast-exposed military veterans, American football athletes, and typical controls which revealed CTE-linked tauopathy as identified by tau-immunoreactive NFTs in the frontal, parietal, and temporal cortices.72 Together, these findings suggest that blast exposure causes neurodegeneration which leads to abnormal neuronal activity, providing a potentially viable non-invasive biomarker of disease states.
4.3 High-frequency mediated DMN dysconnectivity
We observed high frequency pDMN dysconnectivity, converging on the PCC. The PCC is well-established as a multimodal, polysensory, rich club node, with a metabolic rate on average ∼40% higher than the rest of the brain.84 Prior studies have shown PCC dysfunction in numerous acquired brain injury and neurodegenerative states, with PCC lesions severely impacting dynamical functional network configurations associated with general cognitive performance.85,86 Patients with mTBI have shown reduced DMN connectivity involving the posterior cingulate,87 and interhemispheric posterior DMN functional connectivity is predictive of the level of consciousness in TBI.88 Moreover, a longitudinal pre-/post-season fMRI study by DeSimone et al.89 reported that the number of subconcussive RHIs were key in causing reduced DMN functional connectivity (which was not present pre-season) in a group of football players in the absence of concussion after only a single season of play relative to control athletes from non-contact sports, particularly including dysconnectivity with the left PCC. In our study, we found reduced posterior DMN functional connectivity in the gamma band, gamma being a mechanism thought to selectively and flexibly couple proximal or distant cortical regions to temporally organize neuronal activity through the action of interneurons.90
Animal models of tauopathy suggest that tau protein propagation in the brain is promoted by neuronal activity.91–93 A recent study of patients with participants along the AD spectrum (preclinical, mild cognitive impairment, and dementia) indicated that tau (identified with PET) spread along functionally connected (identified by MEG) brain regions.94 Another AD study in amyloid-positive patients showed lower functional connectivity (identified by fMRI) in the DMN for those with elevated tau (identified with PET) but higher functional connectivity for those with lower tau.95 Along this line of work, a study by García-Colomo et al.96 reported longitudinal increases in functional connectivity (identified by MEG) in correlation with tau (identified with PET) in the precuneus in individuals with a family history of AD relative to those without – a brain region which was also identified here as part of the disconnected pDMN. This work suggests that changes in functional connectivity as identified by MEG can be indicative of tauopathy, in turn allowing us to not only track disease progression but also to pre-emptively predict where tauopathy will spread based on functional connectivity. MEG could be utilized in the treatment and diagnosis of brain injury by identifying regions of the brain in affected individuals to deliver precision treatment to reduce further spreading and predict individuals at risk and intervene before tauopathy spreads.
4.4 Association with symptoms
Although worse self-reported neurological symptoms were present with greater blast exposure, there were no significant direct associations between neuronal slowing or dysconnectivity and symptoms. However, neuronal slowing in fronto-temporal and subcortical regions is likely pathogenic in terms of neurobehavioural deficits. The frontal lobe is involved in executive functioning including mental flexibility, goal-directed behaviours, language, learning, and memory, and related processes like attention.97 The hippocampus – part of our subcortical analysis – is involved in multiple aspects of memory, including visual and auditory working and episodic memory,98–101 as well as a breadth of other processes including executive functioning, attention, social behaviours, spatial navigation, and language.101,102 Together, the self-reported cognitive difficulties align with the regional dysfunction we observed. We believe that additional research is needed to explore more fully the association between MEG network functional connectivity and neurological symptom severity, perhaps based on a more comprehensive set of outcome measures known to be impacted by ReBOP.
4.5 Translational use
Notably, MEG but not fMRI revealed pDMN functional dysconnectivity in greater blast exposure, suggesting increased sensitivity to detect network-level effects – these effects were in the gamma range, high-frequency oscillatory activity to which fMRI is blind. Gamma activity reflects the coordinated action of neural excitability and inhibition, driven by local microcircuits involving excitatory pyramidal neurons and inhibitory interneurons, described as the pyramidal-interneuron network gamma (PING) or interneuron network gamma (ING) models.103–105
The application of MEG as a potential surrogate marker for neurodegenerative states would inform translatable paths to procedural and institutional change to reduce the insidious effects of blast exposure in the military and in contact sports. Biomarkers that can identify individuals at risk of neurodegeneration would be critical in preventing further brain injury, and can play a role in preventative medicine. Along with other safety measures – translated from rodent models that suggest head immobilization during blast exposure can reduce blast-related brain dysfunction72 and that jugular compression during breacher training provides preventative utility36 – timely integration of longitudinal monitoring with brain-based biomarkers can improve the operational readiness of warfighters and extend the healthy working lifespan of military members and Veterans.
4.6 Limitations & future directions
Limitations of this study include the use of a subjective self-reported probe of lifetime blast exposure with the GBEV scale; of course, objective blast exposure data – such as lifetime pressure gauge use – would be advantageous, but due to practical limitations, this data is impossible to capture in our current cohort. Second, multimodal MEG and PET imaging data with Flortaucipir could confirm the presence of tau deposition and spatially colocalized neuronal slowing. Confirmatory tau-slowing associations in a single cohort with suspected CTE would be a step change in our understanding of the relationship between tau-induced neurodegeneration and neural slowing, and confirm the latter as a robust surrogate marker for tauopathy. Third, there were no significant relationships between brain measures and self-reported outcomes – a longitudinal approach, with baseline capture before explosives weapons training, would be crucial in establishing slowing as a marker of neurological functioning, given the between-subject heterogeneity and within-subject stability of brain markers, as well as using objective clinical measures (e.g. neuropsychological testing).54
In future work, we plan to examine white matter microstructure in this same cohort, and establish whether white matter damage leads to neural desynchrony, as white matter damage has been reported in contact sport players (football) in the absence of concussions after only a single season of play,106,107 as well as MEG tasks in cognitive control, mental flexibility, and memory. Furthermore, there has been major, recent interest in expanding MEG analyses beyond canonical frequency band, periodic analyses as applied here and quantifying the aperiodic component of the neurophysiological signal recorded by MEG.108 Studies have shown age-related changes in the aperiodic component,109 and aperiodic alterations in AD110 and Parkinson’s disease and dementia with Lewy bodies.111 Finally, future research should also aim to pair MEG with PET-tau and PET-amyloid as seen in the tauopathy studies mentioned above (e.g., Schoonhoven et al.94; Schultz et al.95; García-Colomo et al.96), as well as other various biomarkers,35 including cerebrovascular112 and genetic113 serum and CSF biomarkers.114
4.7 Conclusions
Electrophysiological imaging using MEG produced the first evidence of dysregulated neuronal functioning in a cohort with a history of repetitive subconcussions, including military members and Veterans with varying degrees of exposure to repetitive blast overpressure. Notably, these aberrations were observed independently of concussion or traumatic stress history. In addition to regional slowing of activity as indicated by increased delta activity and functional dysconnectivity in key brain hubs and networks, those with worse blast exposure also exhibited worse neurological symptoms, especially in the somatic and cognitive domains, with no blast-related differences in mental health outcomes. MEG has potential utility as a biomarker for neurodegeneration which can translate to procedural and institutional improvements to minimize harm from subconcussions in military blast exposure and in contact sports.
Data Availability
Defence considerations related to the confidential nature of the human data collected means it is not publicly available.
Funding
This research was funded by awards to BTD from the Department of National Defence, through the Innovation for Defence Excellence and Security (IDEaS) program, Defence Research and Development Canada (DRDC), and industry research funding from MYndspan Ltd.
Competing interests
BTD is Chief Science Officer at MYndspan Ltd. The remaining authors report no competing interests.
Acknowledgements
This research was funded by awards to BTD from the Department of National Defence, and their Innovation for Defence Excellence and Security (IDEaS) program, Defence Research and Development Canada (DRDC), and industry research funding from MYndspan Ltd.
IDEaS is a defence innovation program which invests in research and technology aimed at meeting the demands of the current and future global defence and security environment. Since launching in April 2018, IDEaS has been working with Canadian innovators to develop the S&T landscape, and is helping turn innovative thinking into tangible solutions for the Department of National Defence and the Canadian Armed Forces (DND/CAF), as well as Canadians. Visit the IDEaS web site to get updates on the latest challenges issued by the program.
Footnotes
References
- 1.↵
- 2.↵
- 3.
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.
- 9.
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- 19.↵
- 20.↵
- 21.↵
- 22.
- 23.
- 24.↵
- 25.
- 26.↵
- 27.↵
- 28.↵
- 29.↵
- 30.↵
- 31.↵
- 32.
- 33.
- 34.↵
- 35.↵
- 36.↵
- 37.↵
- 38.↵
- 39.↵
- 40.↵
- 41.↵
- 42.↵
- 43.↵
- 44.↵
- 45.↵
- 46.↵
- 47.↵
- 48.↵
- 49.↵
- 50.↵
- 51.↵
- 52.↵
- 53.↵
- 54.↵
- 55.↵
- 56.↵
- 57.↵
- 58.↵
- 59.↵
- 60.↵
- 61.↵
- 62.↵
- 63.
- 64.
- 65.↵
- 66.↵
- 67.↵
- 68.↵
- 69.↵
- 70.↵
- 71.↵
- 72.↵
- 73.
- 74.
- 75.
- 76.
- 77.
- 78.↵
- 79.↵
- 80.↵
- 81.↵
- 82.↵
- 83.↵
- 84.↵
- 85.↵
- 86.↵
- 87.↵
- 88.↵
- 89.↵
- 90.↵
- 91.↵
- 92.
- 93.↵
- 94.↵
- 95.↵
- 96.↵
- 97.↵
- 98.↵
- 99.
- 100.
- 101.↵
- 102.↵
- 103.↵
- 104.
- 105.↵
- 106.↵
- 107.↵
- 108.↵
- 109.↵
- 110.↵
- 111.↵
- 112.↵
- 113.↵
- 114.↵