ABSTRACT
Quantitative susceptibility mapping (QSM) has been applied to map brain iron distribution after mild traumatic brain in-jury (mTBI), to understand properties of neural tissue which may be related to microstructural damage. However, mTBI is a heterogeneous injury associated with microstructural brain changes, and ‘traditional’ group-wise statistical approaches may lead to a loss of clinically relevant information, as subtle individual-level changes can be obscured by averages and confounded by within-group variability. More precise and individualised approaches are needed to characterise mTBI better and elucidate potential cellular mechanisms to improve intervention and rehabilitation. To address this issue, we build individualised profiles of regional positive (iron-related) magnetic susceptibility across 34 bilateral cortical regions of interest (ROIs) following mTBI. Healthy population templates were constructed for each cortical area using standardised z-scores derived from 25 age-matched male controls, serving as a reference against which z-scores of 35 males with acute (< 14 days) sports-related mTBI (sr-mTBI) were compared. Secondary analyses sensitive to cortical depth and curvature were also generated to approximate the location of iron accumulation in the cortical laminae and the effect of gyrification. Our primary analyses indicated that approximately one-third (11/35; 31%) of mTBI participants exhibited elevated positive sus-ceptibility indicative of abnormal iron profiles relative to the healthy control population, a finding that was mainly concentrated in ROIs within the temporal lobe. Injury severity was significantly higher (p < 0.01) for these mTBI participants than their iron-normal counterparts, suggesting a link between injury severity, symptom burden, and elevated cortical iron. Secondary analyses of cortical depth and curvature profiles revealed abnormal iron accumulation in 83% (29/35) of mTBI participants, enabling better localisation of mTBI-related changes in iron content to specific loci within each ROI and identifying effects that may be more subtle and lost in ROI-wise averaging. Our findings suggest that individualised approaches can further elucidate the clinical relevance of iron in mTBI. Differences in injury severity between iron-normal and iron-abnormal mTBI participants highlight not only why precise investigation is required to understand the link between objective changes in the brain and subjective symptomatology, but also identify iron as a candidate biomarker for tissue damage after mTBI.
INTRODUCTION
Exposure to mild traumatic brain injury (mTBI) is a significant public and personal health concern, accounting for approximately 90% of the 50-60 million annual cases of TBI worldwide.1 Global financial losses related to mTBI are estimated at ~USD $400 billion per year,1,2 however, beyond the economic impacts mTBI can increase the risk of neurodegeneration, dementia,3,4 and premature death.5 In the short term, mTBI can result in a range of symptoms with significant inter-individual variability, including cognitive, emotional, and physiological disturbances such as sleep disruption, light sensitivity, fatigue, headaches, vertigo, vestibular problems, depression, and anxiety, which significantly impact quality of life and participation in day-to-day activities for many.4 In some cases, these symptoms can persist even up to three decades post-injury.6,7 Numerous factors contribute to differences in injury severity, symptom burden, in-vivo brain tissue pathology, and even autopsy findings. These include individual differences prior to injury such as genetic predispositions, age, gender, IQ, psychiatric history, prior exposure to mTBI, and substance use history, as well as differences in the mechanisms and loci of injury.8 In sports-related mTBI (sr-mTBI), for example, variability in the sport and even player position can affect injury severity, lead to diverse effects on brain structure and function, and divergence in symptom burden and cluster.8
The heterogeneity of mTBI is apparent at even the cellular level. The rapid changes in inertia (acceleration/deceleration/rotation) or exogenous skull impact associated with mTBI cause the transmission of mechanical forces to the brain, resulting in a mechanisticallyspecific primary insult and microstructural tissue damage.8,9 This initiates a variable cascade of secondary cellular processes, including disruption of the blood-brain barrier (BBB), cerebrovascular dysfunction, oxidative stress, axonal degeneration, and neuroinflammation9,10 which can propagate for months after the initial impact.11 However, the pathophysiology of mTBI remains poorly understood, and specific biomarkers indicative of mTBI remain, to date, elusive. Unlike moderate-to-severe TBI (ms-TBI), where lesions, haemorrhages, or macroscopic morphological abnormalities can be detected, routine magnetic resonance imaging (MRI) methods are often insensitive to mTBI-related neuropathology.12,13 This limitation necessitates the use of advanced MRI techniques not typically employed in conventional medical settings to identify the subtle changes in brain structure characteristic of this ‘mild’ injury.14,15 Integrating these advanced imaging modalities into routine patient care requires further validation and the establishment of clinically and individually relevant biomarkers for mTBI diagnosis and treatment.
Iron accumulation is increasingly recognised as a component of neuropathology following mTBI, contributing not only to acute-phase secondary injury and later cell death,9 but also cognitive dysfunction after mTBI.16,17 Quantitative susceptibility mapping (QSM) is an advanced MRI technique that can be used to estimate the substrates of biological tissue by leveraging the inherent magnetic properties, such as paramagnetism exhibited by iron in response to an applied magnetic field.18–23 Non-heme iron (particularly ferritin-bound iron), is the main source of paramagnetism on QSM18,24–26 and widely recognised as the form of iron most involved in secondary injury after mTBI.9,27,28 Iron dyshomeostasis can trigger auto-toxic circuits that drive neurodegenerative processes,29 including the generation of reactive oxygen species (ROS), which at high levels can lead to cytotoxic oxidative stress,30 lipid damage, and increased permeability of the cell membrane,9 as well as iron-regulated cell death (ferroptosis).31 As such, elevated levels of iron in cortical regions would suggest localisation of injury-related pathological processes and changes in brain structure. These changes may be related, but not limited to, mTBI-induced permeability of the BBB32 and neuroinflammation,33 both of which are known to be involved in iron accumulation.9,30 Iron has also been implicated in the hyperphosphorylation of tau proteins (p-tau)28 observed in mTBI-related tauopathies; its co-localisation with p-tau thus identifies it as a promising early indicator of neurodegeneration.34,35
A limited number of studies have employed QSM to investigate the role of brain iron in microstructural tissue damage following mTBI, focusing mainly on subcortical nuclei or global grey and/or white matter,36–44 with only a few studies including cerebral regions of interest (ROIs)38 or investigating the relevance of cortical morphology.45 However, the diversity of mTBI effects may not be discernible at the group level, which currently constitutes the standard statistical approach. Individual-level investigations of injury-specific effects may better characterise mTBI-related neuropathology, and are increasingly recognised as providing more biologically informative data than group-level studies, especially in clinical populations where targeted interventions are both useful and necessary,46 such as mTBI. Personalised profiles can be generated by leveraging z-scores to compare the results of individual quantitative measures to the distribution of a healthy normative population. This approach allows for a clearer understanding of where the individual falls relative to normal ranges for selected metrics.47 Individual analytic approaches have been successfully applied in the context of mTBI using T2 relaxometry as a marker of neuroinflammation,48 and diffusion-weighted imaging to investigate white matter fiber tracts.49 Under the TBI umbrella more broadly, individual analyses have been applied to fixel-based analysis of diffusion MRI50 and diffusion tensor imaging51 to investigate white matter integrity, as well as structural connectomics,52 in ms-TBI. One study has used QSM to generate individualised profiles of iron deposition in ms-TBI47 however, to date, dedicated personalised investigations of iron deposition at the individual level following mild TBI are, to the authors knowledge, lacking.
To address these research gaps, we conducted the first dedicated individual-level investigation of iron-related mTBI effects. This study aimed to: 1) generate individual profiles of cortical iron deposition following sr-mTBI, and; 2) extend these findings by deriving profiles sensitive to cortical architectonics (depth and curvature) as a supplemental, secondary approach. Our prior research suggests that elevated magnetic susceptibility should be most evident in ROIs in the temporal lobe.45 However, given the preliminary and exploratory nature of this study, we did not have specific a priori hypotheses about the direction of effects in all cortical regions.
MATERIALS AND METHODS
Ethical approval for this research was obtained from the Health and Disabilities Ethics Committee (HDEC) (Date: 18/02/2022, Reference: 2022 EXP 11078) and institutional approval was also obtained from the Auckland University of Technology Ethics Committee (AUTEC) (Date: 18/02/2022, Reference: 22/12). In accordance with the Declaration of Helsinki, all participants provided written informed consent prior to data collection.
Participants
35 male contact sports players (M = 21.60 years [range: 16-33], SD = 4.98) with acute sr-mTBI (sustained within 14 days of MRI scanning (M = 10.40 days, SD = 3.01)) and 25 age-matched male controls (M = 21.10 years [range: 16-32], SD = 4.35) were recruited for this observational, case-control study (see Table 1). Ages were not significantly different between groups (t(58) = −0.44, p = .66). Clinical (sr-mTBI) participants were recruited through three Axis Sports Medicine clinics (Auckland, New Zealand), via print and social media advertisements, word-of-mouth, and through community-based pathways including referrals from healthcare professionals and sports team management. Each clinical participant was required to have a confirmed sr-mTBI diagnosis by a licensed physician as a prerequisite for study inclusion, and symptom severity was assessed using the Brain Injury Screening Tool (BIST)53 either upon presentation to Axis clinics or electronically following recruitment. Healthy controls (HC) were recruited through print and social media advertisements, and word-of-mouth. Exclusion criteria for all participants included a history of significant medical or neurological conditions unrelated to the study’s objectives and contraindications for MRI. Additionally, controls were excluded if they had any recent history of mTBI events (< 12 months) or were living with any long-term effects of previous mTBI. All participants completed a brief demographic questionnaire and attended a 1-hour MRI scan at The Centre for Advanced MRI (CAMRI), Auckland, New Zealand. All scans were reviewed by a certified neuro radiologist and consultant for clinically significant findings. There were no clinically significant diagnoses identified from MRI in either group. While some incidental findings were identified, none were considered to be clinically significant and so no further follow up action was needed (see Table 1).
Neuroimaging
Details on image acquisition and processing have been previously reported,45 and are summarised here for brevity.
Acquisition
MRI data were acquired on a 3T Siemens MAGNETOM Vida Fit scanner (Siemens Healthcare, Erlangen, Germany) equipped with a 20-channel head coil. A 3D flow-compensated GRE sequence was used to obtain magnitude and unfiltered phase images suitable for QSM reconstruction. Data were collected at 1 mm isotropic voxel size with matrix size = 180 x 224 x 160 mm, TR = 30 ms; TE = 20 ms; FA = 15°; FoV = 224 mm in a total acquisition time of ~ 3.43 minutes. For each participant, a high-resolution 3D T1-weighted (T1w) anatomical image volume was acquired for coregistration, parcellation and segmentation using a Magnetisation-Prepared Rapid Acquisition Gradient Echo (MPRAGE) sequence (TR = 1940.0 ms; TE = 2.49 ms, FA= 9°; slice thickness = 9 mm; FoV = 230 mm; matrix size = 192 x 512 x 512 mm; GRAPPA = 2; voxel size 0.45 x 0.45 x 0.90 mm) for a total acquisition time of ~ 4.31 minutes. DICOMs were converted to NIfTI files and transformed to brain imaging data structure (BIDS)54 for further processing using Dcm2Bids55 version 3.1.1, which is a wrapper for dcm2niix56 (v1.0.20230411).
Image processing
Bias field-corrected57,58 T1w images were processed in FreeSurfer59 to 1) delineate pial and grey matter/white matter (GM/WM) boundary meshes, and; 2) generate estimates of cortical thickness and curvature for each vertex.60 QSM images were reconstructed using a rapid open-source minimum spanning tree algorithm (ROMEO),61 background field removal with projection onto dipole fields (PDF),? and sparsity-based rapid two-step dipole inversion (RTS);62 a pipeline congruent with recent consensus statement recommendations for best-practice QSM reconstruction.63 All QSM reconstruction was carried out via QSMxT64 v6.4.2 (https://qsmxt.github.io/QSMxT/) and used a robust two-pass combination method for hole-filling and artefact reduction.65
Subsequent processing was performed using FSL.66–68 For each subject, the raw magnitude image was skull-stripped69 and binarised. These binary masks were used to erode non-brain signal around the brain perimeter using fslmaths. Magnitude images were linearly coregistered to the T1w image using FMRIB’s Linear Transformation Tool (FLIRT)70–72 with 12 degrees of freedom. Due to variability in acquisition type, FoV and matrix size between subjects’ QSM and T1w images, the 12 DoF linear registration provided more accurate alignment compared to the 6 DoF alternative, allowing for better compensation of non-rigid anatomical variations upon visual inspection. The resulting transformation matrix was used for spatial normalisation of the QSM images to T1w space. In line with prior research,60 QSM maps were then thresholded into separate inter-voxel sign (positive and negative) maps with fslmaths. Traditional QSM maps already represent average voxel-wise susceptibility values,73 this averaging can obscure individual susceptibility sources as well as introduce further confounding effects when further averaging across voxels during analysis. Thresholding may help address the latter limitation by isolating voxels containing paramagnetic substrates, such as iron, which could enable more targeted analyses of these specific susceptibility sources. Only positive sign maps were used in further processing and analyses to target cortical iron distribution.
Cortical column generation
To generate cortical columns and sample positive susceptibility values, we adapted a pipeline previously applied to DWI74 and QSM45,60 for depth- and curvature-specific cortical analysis. First, the T1w FreeSurfer59 recon served as an input into the easy_lausanne tool (https://github.com/mattcieslak/easy_lausanne.git) based on the open-source Connectome Mapper75 to separate the cortex into 34 ROIs per hemisphere according to the Lausanne multi-scale atlas (equivalent to the Desikan-Killiany atlas76 native to FreeSurfer).59
Depth
Cortical columns were created for each hemisphere in T1w space with write_mrtrix _tracks77 in MATLAB (version R2024a), which was used to connect vertex pairs between the pial and GM/WM boundary surface meshes. Each cortical column was segmented into 6 equidistant depths extending from the pial surface to the GM/WM boundary using MRtrix3 tckresample.77 This approach differs from previous studies utilising this technique, where 21 equidistant depths were employed.45,60,74 Here, we used 6 depths to decrease the number of depth-wise comparisons, and to better approximate the structure of the intracortical layers.78,79 It should be noted here that these represent equidistant segmentations rather than specific cellular laminae (L1 to L6) of the cortex as it is important to distinguish results produced using this approach from ultra-high field investigations of cyto- and myelo-architecture in the cerebral cortex; results described herein are related to cortical depth, rather than layer.
Curvature
The columns were also categorised based on cortical curvature, derived from FreeSurfer’s59 Gaussian curvature values at each GM/WM boundary vertex80 and quantified in units of 1/mm2. The categories included the gyral crown (curvature values: −0.6 to −0.1), sulcal bank (−0.1 to 0.1), sulcal fundus (0.1 to 0.6), or whole-ROI (−0.6 to 0.6).60 Positive curvature values indicated sulci, while negative values indicated gyri, with higher values corresponding to deeper curvatures.60 Only columns ranging from 0.5 mm to 6 mm in length were included in the analysis to capture plausible cortical morphology.81 Depth was measured in percentage of cortical thickness rather than absolute metrics (mm) to mitigate any variability between participants.
Personalised QSM profiles
We generated individual QSM profiles for each ROI at the bilateral level using MATLAB (2024a) (see Figure 1 for visualisation). Mean positive susceptibility values were extracted across the whole ROI (curvature and depth combined), as well as 3 curvature bins (gyral crown, sulcal bank, and sulcal fundus) and 6 cortical depths independently for all 34 ROIs. For the whole-ROI profiles, z-scores were calculated for all participants (both healthy controls and mTBI), by subtracting the HC group mean from each individual’s sus-ceptibility value and dividing by the HC group standard deviation; a method commonly used in prior research.47,48,51,52To bring the HC data closer to a normal distribution, outlier scores for the HC group were filtered47 if they fell outside 2 times the interquartile range (IQR), a more stringent criterion than the methods used to identify mild outliers 1.5 times the IQR, but less extreme than the more conservative filter of 3 times the IQR.82 As a result, data from n = 1 HC participants were excluded in 3 of the 34 ROIs, and data from n = 2 HC participants were excluded in 1 of the 34 ROIs. After filtering, the Shapiro-Wilk normality test yielded an average W value of M = 0.96 (SD = 0.02) across z-distributions for all ROIs, indicating that the data distribution within each ROI closely approximates normality. The final equation for deriving the whole-ROI z-scores for individual mTBI participants was as follows: where zmTBI represents the ROI-wise z-score for each individual mTBI participant; XmTBI is the ROI-wise mean QSM value for each mTBI participant; μHCnorm is the mean ROI-wise QSM value of the HC group after outlier filtering to normalise the distribution, and; σHCnorm is the ROI-wise standard deviation of the HC group QSM values after outlier filtering. This approach ensures that mTBI participants’ susceptibility values are directly comparable to the healthy range reflecting a normalised distribution, allowing the detection of deviations that might reflect underlying pathology.
STATISTICAL ANALYSES
To assess statistical significance for whole-ROI mTBI z-scores, two-tailed p-values were calculated from the z-scores using the cu-mulative distribution function of the standard normal distribution. A false discovery rate (FDR) correction 83 was applied to the p- values for 34 ROI-wise comparisons. The same process was repeated for each depth at each curvature bin, however, the IQR filter was omitted due to the number of comparison points.
To conduct secondary exploratory statistical tests, we divided the mTBI group into two subgroups: those whose z-scores significantly deviated from HC norms in the whole-ROI analysis (i.e., iron-abnormal) and those whose scores did not (i.e., iron-normal). Although there was no statistically significant difference in age between mTBI participants and controls, we performed a one-way analysis of variance (ANOVA) across these three groups to confirm that age was not driving the observed results at the individual level. Additionally, we used an independent samples t-test to assess whether injury severity (BIST53 scores) differed significantly between iron-abnormal and iron-normal mTBI participants, excluding mTBI-16 for this analysis only due to missing BIST data.
RESULTS
Regional individualised cortical iron profiles
We calculated personalised profiles of iron-related differences in positive susceptibility across 34 cortical ROIs for each mTBI participant, to understand the effects of mild brain trauma at the individual-level. Of the 35 mTBI participants, 11 (31%) exhibited significantly elevated positive susceptibility for at least one ROI relative to the healthy control population template (see Table 2), likely indicating elevated iron. No clinical participants’ z-scores were significantly lower than the healthy control population.
In these 11 clinical participants, injury-related elevated susceptibility was evident across all cortical lobes. Notably, a high density of affected ROIs was observed in the temporal lobe for 82% (9 out of 11) of participants with abnormal iron profiles (see Table 2). In contrast, 45% (5/11) had abnormal iron in occipital ROIs, 27% (3/11) in frontal ROIs, 18% (2/11) in the insula or cingulate, and only one participant (9%) had an abnormal profile inclusive of parietal ROIs. This variability reflects the diverse regions impacted by exposure to mild brain trauma (see Figure 2). For example, elevated cortical iron was localised to various temporal lobe ROIs in participants mTBI-03 (parahippocampal gyrus), mTBI-15 (middle temporal gyrus and entorhinal cortex), mTBI-16 (temporal pole), and mTBI-17 (temporal pole). In contrast, mTBI-30 showed extensive frontal lobe effects (caudal medial prefrontal cortex (mPFC), pars opercularis, and superior frontal gyrus), while the remaining participants exhibited multi-focal elevated positive susceptibility affecting two or more lobes (see Table 2 and Figure 2).
After subdividing the mTBI participants into ironnormal (24/35; 69%) and iron-abnormal (11/35; 31%) based on their individual ROI-wise profiles, we used a one-way ANOVA to test if age was significantly different between the three groups (iron-normal mTBI, iron-abnormal mTBI, and controls) as well as an independent samples t-test to understand whether injury severity (BIST53) was significantly different for iron-normal and iron-abnormal mTBI participants. The ANOVA yielded a non-significant effect of age, F(2, 26) = 2.00, p = .16. BIST53 scores were significantly higher for the iron-abnormal mTBI group (M = 59.6, SD = 40.5) than the iron-normal mTBI group (M = 29.2, SD = 23.3), t(32) = 2.77, p = < .01, suggesting a link between injury severity and elevated regional cortical iron deposition (see Figure 3).
Exploratory depth- and curvature-specific individualised iron profiles
Exploratory analyses sensitive to cortical depth (1 through 6, from the pial surface to the GM/WM interface) and curvature (crown, bank, and fundus) further highlighted the heterogeneity of iron deposition in mTBI. Only 17% (6/35) of participants retained normal iron profiles; the remaining 83% (29/35) showed elevated susceptibility in at least one ROI, for at least one depth, and for at least one curvature. Isolated instances of negative z-scores (indicating lower iron compared to HC) were observed in 7 of the 29 participants, but were typically limited to a single ROI/depth combination. Overall, abnormal iron accumulation was most pronounced in the sulcal fundus, followed by the sulcal bank, and was least evident in the gyral crown (see Figure 4 for reference). These distributions were consistent across both superficial (1-2), mid (3-4), and deeper (5-6) cortical depths, however, each unique combination of depth and curvature for each participant reflects individual variability in iron deposition across the cortex. For example, abnormal iron deposition for mTBI-25 was exclusive to superficial depths (1 and 2) in the gyral crown and sulcal fundus of the ironelevated ROIs. In contrast, mTBI-34 exhibited more widespread iron accumulation which became more pronounced at deeper depths, and was observed in the gyral crown, sulcal bank, sulcal fundus, or a combination thereof, depending on the ROI. Participant mTBI-29 showed iron accumulation at all depths but this was most pronounced at mid-to-deep depths (3-5) in the sulcal bank and fundus (see Figure 4).
For the 11 participants who demonstrated increased iron in the initial whole-ROI analyses, a similar ROI-wise distribution of elevated iron was observed in the depth- and curvature-specific analyses. For example, in whole-ROI analyses mTBI-30 showed elevated susceptibility in the caudal mPFC, pars opercularis, and superior frontal gyrus (see Table 2 and Figure 2). Depth and curvature analyses revealed that iron deposition in the caudal mPFC occurred across superficial-to-deep depths (2-5) in the sulcal bank and at mid-to-deep depths (3-6) in the sulcal fundus, with no deviation from the control template at the gyral crown (see Figure 4). In the pars opercularis, iron accumulation was apparent at mid-to-deeper depths (4-5) in the bank and at superficial-to-mid depths (2-4) in the sulcal fundus, again with no accumulation at the gyral crown. Similarly, the superior frontal region showed mid-to-deep (4-5) abnormalities in the bank. Regions that did not appear in the initial analyses, including the insula, lateral occipital cortex (LOC), pars triangularis, posterior cingulate, precentral, and superior parietal areas, exhibited elevated iron in depth- and curvature-specific analyses, mostly concentrated in the sulcal bank or fundus. Conversely, mTBI-25 exemplifies a case where no abnormal iron profile was detected in the whole-ROI analysis, but became evident with analyses sensitive to anatomical morphology (see Figure 4).
DISCUSSION
Previous studies using QSM to examine the role of brain iron following mTBI have primarily focused on subcortical brain areas or global grey and/or white matter.36–44 Only two investigations have included cortical ROIs and, of these, only our previous work45 accounted for anatomical variations in cortical depth and curvature. However, these studies relied on group-level statistical analyses, which can obscure individual brain changes due to the subtle nature of cell damage associated with mTBI. This approach may limit our understanding of this heterogeneous injury and inhibits the identification of personalised injury profiles, hindering the implementation of individualised rehabilitation strategies and treatments. We believe, therefore, that research that incorporates comparisons of individual clinical participant data to healthy normative ranges can play a key role in informing individualised neural and pharmacological interventions. We conducted the first investigation of individual differences in cortical magnetic susceptibility after sports-related mTBI across 34 cortical ROIs, using a healthy population template as a reference. Secondary exploratory analyses of individual differences in susceptibility distribution at different cortical depths and curvatures were included to better characterise injury profiles for each participant
Our findings revealed that a substantial subset of individuals with mTBI exhibit elevated levels of cortical magnetic susceptibility, indicating injury-related iron accumulation. Our primary investigation evidenced abnormal iron profiles in just under one-third (31%; 11/35) of participants relative to the healthy control population. In these 11 ‘iron-abnormal’ clinical participants, elevated susceptibility was evident across all cortical lobes, with a notably higher density of affected ROIs localised to the temporal lobe. At the individual level, elevated iron profiles were either focal to specific lobes (such as the temporal lobe in mTBIs −03, −15, −16, and −17, and widespread frontal anomalies in mTBI-30) or bi-focal (affecting two lobes). Only one participant (mTBI-29) exhibited multi-focal cortex-wide elevated susceptibility involving the frontal, temporal, parietal, occipital, and insular/cingulate cortices. These results suggest that whilst there is significant inter-subject variability, iron accumulation after mTBI preferentially affects temporal regions. Additionally, we found that mTBI participants with abnormal iron accumulation experienced more severe symptoms. Taken together, our findings support an iron-related mechanism of secondary injury that modulates symptom presentation.
Iron dyshomeostasis in the basal nuclei is known to impair cognitive function after mTBI,16,17 however, little is known about the effect of cortical iron aggregation on mTBI symptomatology. In our previous work,45 we reported few correlations between regions of cortical iron accumulation and BIST53 scores, a measure of injury severity and dominant symptom cluster (physical, vestibular, and cognitive). This lack of correlation may be accounted for by known variability in the accuracy of symptom reporting or purposeful underreporting of symptoms, a phenomenon notorious among sports players.84 However, group-level examinations may also obscure individual differences in addition to inhibiting the implementation of more appropriate statistical approaches. By assessing the effect of mTBI at the individual level, we were able to facilitate more precise between-group analyses of injury severity by differentiating between iron-normal and iron-abnormal mTBI participants. Our results revealed a significantly higher symptom burden for participants with mTBI when their iron profiles were also abnormal.
The field generally lacks reliable correlations between subjective assessments of injury severity and objective measures of brain injury and recovery,85 as well as alignment between cognitive or clinical findings and neuroimaging results.86 Identifying reliable, objective markers of structural changes that are related to subjective self-reported symptoms is crucial because individual variations in brain injury location and severity can lead to vastly different clinical presentations and recovery trajectories but may be missed in group-level analyses.8 Although research indicates that most individuals recover well from mTBI, between 15%11 and 30%8 of patients experience significant, and in some cases life-changing, long-term clinical sequelae. Understanding the underlying pathophysiological drivers of these poorer outcomes is essential for enabling precise, patient-specific clinical interventions. We found that 31% of participants exhibited abnormal iron profiles substantial enough to be detected in the primary ROI-wise analyses, which were linked to poorer outcomes. This aligns with findings from other individualised studies, which have reported that a similar percentage (28%) of participants with ms-TBI exhibit white matter anomalies in the subacute phase, associated with worse cognitive outcomes than participants without.51 In addition, standardised susceptibility values in the basal nuclei are reported to correlate with mTBI symptom duration, but only in a sub-group of participants with persistent symptomatology for at least a week.40 In keeping with these findings, our results reinforce the importance of individualised analyses in revealing associations between the extent of microstructural pathology and negative outcomes following mTBI as well as the ability of this approach to differentiate at-risk sub-cohorts. Here, the value of personalised approaches to understanding mild brain injury becomes strikingly apparent. Further longitudinal studies that track participants through recovery could help to determine whether elevated iron levels are associated with prolonged recovery, persistent post-concussive symptoms, or adverse outcomes later in life. Such research may offer insights into why a subset of individuals with mTBI fail to recover fully.
The promise of individualised assessments to identify biomarkers for mild brain injury is particularly salient given the current absence of objective markers for mTBI diagnosis. Diagnostic decisions are limited to subjective self-report and assessments of physiological function,13,85 as heterogeneity of mTBI complicates efforts to identify reliable biomarkers or imaging signatures that can be applied universally across patients. Whether cortical iron accumulation reflects inflammatory processes, BBB disruption,9 ferroptosis,31,87 or other cytotoxic processes is beyond the scope of the current research. However, we posit that the significantly higher symptom burden observed in the iron-abnormal mTBI cohort supports an iron-mediated mechanism of brain changes related to injury severity and functional impairment, and marks iron as a promising mTBI biomarker. Of particular note, many of the symptoms observed in mTBI resemble those seen in other iron-mediated processes, such as cognitive decline in normal aging,88 and the cognitive and motor dysfunctions characteristic of neurodegenerative diseases like Alzheimer’s, Huntington’s, Parkinson’s, and Friedrich’s ataxia, as well as multiple sclerosis, where iron dysregulation is a hallmark feature.9,89 Genetic disorders of iron over-load, such as neuroferritinopathy, also present with cognitive and motor symptoms.90
The susceptibility of the temporal lobe to iron accumulation post-mTBI, evident in this and previous research,45 aligns with the memory impairments characteristic of mTBI.91 Individual-level data from our mTBI sample emphasises the link between high temporal iron accumulation and memory deficits, as exemplified by BIST53 scores related to memory (i.e., mTBIs −01, −03, −04, −05, −15, −17, −29, and −34; see Table 2). The distributions of p-tau focal to temporal (and frontal) cortex are considered pathognomonic of CTE,92,93 an mTBI-related neurodegenerative disorder. Co-localisation of iron with p-tau in CTE has been highlighted in histological examinations.34 Whilst the precise relationship between iron overload and the downstream hyperphosphorylation of tau proteins remains an active area of research, further exploration of the interplay between acute cortical iron elevation, symptom burden, temporal recovery dynamics, and long-term brain health outcomes is warranted.
The diverse range of mechanistic antecedents, pathological mechanisms, and clinical outcomes associated with mTBI reflects the complexity of the underlying pathophysiology.8 mTBI is not a uniform injury, and structural indicators, such as those commonly observed in other forms of TBI, do not always correlate with clinical symptoms or outcomes.94 Here, we speculate that abnormal iron accumulation in specific cortical ROIs may be related to participant symptomatology. For instance, the superior temporal sulcus (STS) plays a key role in social cognition, empathy, mentalising about others’ emotional states, and ‘theory of mind.’95–97 Structural changes to this region may explain the severe irritability reported by mTBI-01 (see Table 2), along with complaints of severe tinnitus.98 As a hub for audiovisual integration,99 the superior temporal region could also be involved in phonophobia (sound sensitivity), visual disturbances, and vestibular dysfunction;100 symptoms experienced by both mTBI-01 and mTBI-29 whose profiles indicated local elevations in cortical iron. Similarly, the lingual gyrus is active during migraine episodes and responds to luminous stimuli, suggesting its involvement in photophobia, visual processing anomalies, and cephalalgia (headache).101 These symptoms were observed in all clinical participants (mTBIs −01, −04, −05, −18; see Table 2) with iron aggregation in this region. Our findings allude to specific areas of cortical iron accumulation that show a relationship to clinical sequelae, suggesting that regions with higher iron burden may be evidence of microstructural cell damage that disrupts normal function. Future research should integrate functional MRI to enhance the translation of findings and improve the mapping of structural changes to deficits in functional connectivity and brain networks. Without additional research, observations about regions of iron accumulation and symptom presentation are speculative and much less convincing than case-matching between brain lesions in gross TBI and neurobehavioural symptom presentation.102
Depth- and curvature-specific iron accumulation
Secondary analyses revealed significant inter-individual heterogeneity in depth- and curvature-specific cortical iron accumulation, however, there was a consistent trend for iron deposition in the sulci, particularly concentrated at the fundus, followed by the sulcal bank. This pattern may be attributed to the heightened vulnerability of the sulcal fundus to mTBI-related injury, which is sus-ceptible to mechanical deformation due to the ‘water hammer effect,’ where cerebrospinal fluid (CSF) is forced into the depths of the sulci, causing local damage.103 Supporting this, previous research has demonstrated that mTBI increases cortical curvature in the sulcus,104 as well as widens the sulci and causes focal vascular injury and microhemorrhages in the fundus, as evidenced by susceptibility-weighted imaging (SWI).103 Conversely, in more severe TBI, contusions are often concentrated at the gyri,91 suggesting that gyral iron accumulation may represent a less severe version of this type of injury. This highlights how crucial injury biomechanics may be to understanding antecedents to differences in the location of neuropathology. Research shows that different types of head impacts can result in varying brain deformations and injury patterns, with sulci being particularly vulnerable to mechanical strain, which is consistent with, and can be predicted by, patterns of tauopathy observed in neurodegenerative conditions.105 Personalised iron accumulation patterns may provide insight into the specific injury mechanisms and related cellular disruption experienced by each individual participant at the acute stage. Elucidating this link, as well as how this may be related to mTBI-induced neuropathology in later life, should be a focus for future research.
Histological studies, including Perls iron staining and ultra-high-field (7T) R2* mapping of tissue samples106 have localised iron deposition to specific cortical layers, reflecting distinct cyto- and myelo-architectural features with layer-specific distributions that show congruence with in-vivo QSM.107 In healthy populations, iron concentrations typically increase from the pial surface to-ward the GM/WM boundary; deviations away from baselines for each layer suggest an injury-specific model of cortical cellular trauma. For instance, Layer I primarily contains axons and dendrites, with the cell bodies of these processes residing in deeper layers,108 and iron accumulation in different layers may point to diverse pathologies affecting different parts of the cell. Depth-wise patterns may also be related to injury biomechanics: superficial iron accumulation may be a result of perivascular trauma, which is often linked to microhemorrhages and microglial activation after mTBI,27,28,91 whereas deeper iron deposition may reflect more severe shear forces, which are known to cause significant damage near the GM/WM interface in mTBI.109 This is supported by computational modelling showing shear forces are concentrated in this region,105 which are also a common site of microbleeds.110 Tying this in with cortical curvature, this would most likely be apparent at the depths of the sulci.103,105,111 However, contusions of the cerebral gyri damage often follow a layer-specific pattern where most damage occurs at the superficial crest but extends through the cortical mantle to the GM/WM boundary in a ‘wedge’ of haemorrhage and necrotic tissue.91 Iron-dependent cell death, ferroptosis,31,87 may account for some of the similarities observed between this and less severe forms of TBI.
It is important to note that both injuries at the pial surface111 and closer to the GM/WM border109 have been related to adverse outcomes after mTBI. However, subpial iron deposition near small blood vessels in the fundus is congruent with the pathognomonic distributions of sulcal tauopathy seen in CTE,111 which may be related to breaches in the BBB.105 Injury-induced microvascular dysfunction may increase BBB permeability112,113 which likely increases active transport of non-heme iron via vascular endothelial cells into the superficial layers of the cortex9,29,114which could account for perivascular accumulations.91,92 The potential long-term effects of iron deposition in this young cohort warrants further investigation, given the known relationship between iron accumulation and tau pathology.34
LIMITATIONS AND FUTURE RESEARCH
The ability of neuroimaging modalities to infer underlying biological processes is inherently limited, primarily by spatial constraints and inference about tissue composition from indirect markers. While the use of QSM can provide insight into potential iron accumulation, it remains a surrogate measure of iron-tissue content. Integrating complementary modalities, such as positron emission tomography (PET), would offer a more comprehensive understanding of the underlying neurobiology, for example by using tracers sensitive to inflammation or tau115 which may co-localise with iron overload. Additionally, the inclusion of blood-based biomarkers or other protein assays would enhance the interpretation of the neurobiological consequences of sr-mTBI, providing a clearer link between observed imaging abnormalities and molecular pathology. The cross-sectional nature of this study restricts our ability to establish a causal relationship between iron accumulation and symptom progression; here, longitudinal studies that track participants over time would be more effective in determining how iron deposition influences recovery and long-term outcomes, potentially also offering insights into the development of neurodegenerative conditions. This study does not address the precise mechanisms by which different types of mechanical strain or impact contribute to the observed depth- and curvature-specific patterns of iron accumulation. Whether these patterns are driven by depth-wise shear forces near the GM/WM interface105,109 or reflect superficial perivascular damage,27,28,91 or whether curvature-related deposition may be due to the ‘water hammer effect’ at the sulcal fundus103 or contusions of the gyri,91 or a combination of all of these factors, remains unclear. Future research should track athletes throughout the sporting season to better understand the specific biomechanics of their injuries at the time of insult. Longitudinal assessments here would also be useful for elucidating how patterns of iron deposition and the biomechanics of injury relate to the development of neuropathology across different cortical depths and layers.
It should be noted that this study did not control several potential confounding factors, including prior injuries, genetic predispositions, and environmental influences, all known to affect injury severity and recovery trajectories.8 Future studies may benefit from incorporating measures to account for these variables. Another limitation relates to the generation of a normal healthy control z-distribution, as there is no widely established best-practice method for this process. While this study filtered for outliers using robust techniques at two times the interquartile range, other individualised QSM research has employed less stringent criteria, such as filtering at three times the interquartile range to detect more severe TBI pathology.47 Future research would benefit from standardised practices for outlier detection. In addition, this research focused solely on net paramagnetic susceptibility voxel values, leaving net diamagnetic voxels, containing substances such as certain proteins or calcifications, unexplored. Previous studies in both mTBI45 and Alzheimer’s disease (AD)60 have suggested interactions between iron deposition and changes in diamagnetism on QSM. Future research should incorporate an analysis of both paramagnetic and diamagnetic susceptibility distributions to better characterise tissue changes after injury.
The current study also relied on a single-echo QSM sequence, which constrains susceptibility thresholding to a between voxels approach. Multi-echo sequences, which enable the separation of susceptibility sources within voxels,107,116,117 should be used in future research to further refine magnetic source separation. This approach would allow for more precise mapping of microstructural changes and improve the biological interpretation of the data. While this study approximated cortical architectonics by analysing six cortical depths, voxel resolution at a 3T field strength is limited and cannot delineate cortical laminae. Higher-field imaging (e.g., at 7T) would enable better resolution and alignment with cortical architectonics and should be considered for future research using a depth-wise approach to iron mapping. Finally, the focus on the male participants for this research limits the generalisability of the findings to females for whom hormonal variations,118 the use of oral contraceptives,119 and differences in neck musculature120 lend themselves to further heterogeneity in injury severity and outcome. Planned future research will extend these investigations to include female athletes to address this gap.
CONCLUSIONS
To investigate the mechanisms of iron dyshomeostasis in acute mild traumatic brain injury, we conducted the first QSM study to assess individual profiles of cortical iron deposition. Our findings revealed significant heterogeneity in iron accumulation, which may be influenced by cortical morphology, highlighting the importance of examining mTBI at an individual level rather than relying on group analyses. This variability likely complicates the search for universal biomarkers, further underscoring the need for a personalised approach that integrates advanced imaging and detailed symptom profiling. Our results emphasise the need for more targeted, individualised interventions to improve outcomes based on personal profiling, and suggest that iron-mediated cell damage plays a key role in mTBI pathology.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
AUTHOR CONTRIBUTIONS
Christi A. Essex (Conceptualization, Methodology, Project Administration, Validation, Software, Formal Analysis, Investigation, Re-sources, Data Curation, Writing - Original Draft, Writing - Review & Editing, Visualization); Devon K. Overson (Methodology, Software, Visualization, Validation, Writing - Review & Editing); Jenna L. Merenstein (Methodology, Writing - Review & Editing, Visual-ization); Trong-Kha Truong (Methodology, Software, Writing - Review & Editing); David J. Madden (Methodology, Software, Writing - Review & Editing); Mayan J. Bedggood (Writing - Review & Editing, Project administration, Investigation); Catherine Morgan (Methodology, Writing - Review & Editing); Helen Murray (Writing - Review & Editing); Samantha J. Holdsworth (Writing - Review & Editing); Ashley W. Stewart (Writing - Review & Editing); Richard L.M. Faull (Writing - Review & Editing); Patria Hume (Writing - Review & Editing); Alice Theadom (Conceptualization, Methodology, Writing - Review & Editing, Funding acquisition, Supervision); Mangor Pedersen (Conceptualization, Methodology, Writing - Review & Editing, Funding acquisition, Supervision)
FUNDING
This project was funded by a grant from the Health Research Council of New Zealand (HRC), grant #21/622.
CONFLICT OF INTEREST
None to declare.
DATA AVAILABILITY
De-identified MRI data and code used for image processing and statistical analysis can be made available upon request to the corresponding author. Parent codes for cortical column generation can be made available upon request from co-authors (JLM, TKT) based at the Brain Imaging and Analysis Center at Duke University Medical Center.
ACKNOWLEDGMENTS
We extend thanks to Amabelle Voice-Powell and Cassandra McGregor for their contribution to the data collection, and Tania Ka’ai for bringing her perspective to cultural considerations on this study. In addition, we thank Axis Sports Concussion Clinics, particularly Dr Stephen Kara, for their assistance with recruiting sr-mTBI participants and personnel at the Centre for Advanced Magnetic Resonance Imaging (CAMRI) for their assistance collecting MRI data. We also acknowledge Dr Tim Elliot for radiological reporting of all participants and Siemens Healthineers for the use of a work-in-progress (WIP) prototype sequence for the acquisition data used to perform QSM.
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].↵
- [115].↵
- [116].
- [117].
- [118].↵
- [119].↵
- [120].↵