Abstract
Background
A significant number of patients with mild traumatic brain injury (mTBI) would experience cognitive deficit.
Purpose
To investigate the brain structural changes in sub-acute mTBI by diffusion kurtosis imaging (DKI) and volumetric analysis, and to assess the relationship between brain structural changes and cognitive functions.
Material and Methods
A total of 23 patients with sub-acute mTBI and 24 control participants were recruited. All the participants underwent examinations of neuropsychological tests, DKI, and magnetic resonance imaging (MRI)-based morphological scans. Images were investigated using whole brain-based analysis and further regions of interest-based analysis for subcortical nuclei. The neuropsychological tests were compared between the mTBI and the control group. Correlation analysis was performed to examine the relationship between gray matter (GM) volume, DKI parameters, and cognitive functions.
Results
Compared with control participants, mTBI patients performed worse in the domains of verbal memory, attention and executive function (P < 0.05). No regional GM volume differences were observed between the mTBI and control groups (P > 0.05). Using DKI, patients with mTBI showed lower mean kurtosis (MK) in widespread white matter (WM) regions and several subcortical nuclei (P < 0.05), and higher mean diffusivity (MD) in the right pallidum (P < 0.05). Lower MK value of multiple WM regions and several subcortical nuclei correlated with cognitive impairment (P < 0.05).
Conclusion
DKI was sensitive in detecting brain microstructural changes in patients with sub-acute mTBI showing lower MK value in widespread WM regions and several subcortical nuclei, which were statistically associated with cognitive deficits.
Keywords
Introduction
Mild traumatic brain injury (mTBI) is a major public health problem, affecting approximately 42 million individuals worldwide, annually (1). Despite its name, mTBI could lead to various neurological and non-neurological disorders, among which cognitive impairment is the most common one (2). Previous studies revealed that a significant number of patients with mTBI experienced cognitive deficit during the first two weeks, and the cognitive deficit would persist for about three months after the injury (3,4). The possible mechanisms for cognitive impairment include both gray matter (GM) and/or white matter (WM) damage. Finding a reliable biomarker for accurate assessment of the brain pathological changes underlying cognitive deficit after mTBI is of great significance.
Magnetic resonance imaging (MRI), as a non-invasive tool, is playing an increasingly important role in the assessment of TBI from its acute to chronic phases (5). Notably, diffusion tensor imaging (DTI), which describes water diffusion in biological tissues, was frequently used to quantify the microstructural changes of WM (6). Previous studies have proved the value of DTI to assess WM microstructural changes in sub-acute mTBI using a region of interest (ROI) or tract-based spatial statistics (TBSS) method, and the DTI parameters changes in certain WM tracts were associated with cognitive impairment (7–10).
However, DTI is insufficient to reflect the actual non-Gaussian diffusion of water molecules in brain tissues. In particular, it cannot accurately estimate the circumstance of WM with multiple fibers in different directions. Diffusion kurtosis imaging (DKI), in contrast, is a method that uses the non-Gaussian model of water diffusion that could provide complementary information such as mean kurtosis (MK) to depict the complexities of the brain microstructure, especially for GM (11). Theoretically, DKI is a more accurate method than DTI for investigating the pathological alterations in brain tissues. Presently, only a few studies have proved that DKI could provide more information regarding the mTBI-induced brain microstructural changes (12–15).
In addition to WM tracts, GM damage was also observed in patients with mTBI. Loss of GM volume has been shown in patients with chronic mTBI by using a voxel-based morphometry (VBM) method (16–18). However, whether the volume loss in mTBI begins to appear in an earlier sub-acute stage and its relationship with cognitive function is unknown. Furthermore, the study by Zagorchev et al. (19) observed volume loss in subcortical nuclei including hippocampus, amygdala, caudate, putamen, and thalamus in patients with chronic mTBI. As for the microstructural changes, only a few studies proved that DKI parameters in the thalamus might be useful for early prediction of persistent brain damage and cognitive outcome after mTBI (12,13,20). Other subcortical nuclei have not been studied using the DKI technique in cognitive impairment after mTBI.
Determining brain structural changes that underlie cognitive impairment will help reveal the pathophysiology of cognitive deficits due to mTBI. The aim of the present study was to combine DKI and volumetric analysis to assess the GM and WM microscopic structure changes in patients with sub-acute mTBI and, if present, to correlate these changes with the neuropsychological performance. We hypothesized that: (i) brain volume and DKI parameter alterations could be found in sub-acute mTBI and would be associated with cognitive impairment; and (ii) compared with DTI, DKI could detect subtle microstructural changes in sub-acute mTBI. These findings would enhance our knowledge regarding cognitive impairment due to mTBI.
Material and Methods
Participants
The present study was performed in accordance with the recommendations of the local Ethics Committee with written informed consent from all the participants in accordance with the Declaration of Helsinki. The study included 23 consecutive patients with mTBI in the sub-acute stage (from two weeks to two months) recruited from the neurosurgery clinic. mTBI was diagnosed based on the Mayo classification system for injury severity (21). In addition, 24 healthy participants without any history of mTBI were also included in this study. All the participants were assessed using a battery of neuropsychological tests followed by an MRI examination. The exclusion criteria were as follows: a history of drug abuse; cerebrovascular disease; psychiatric or neurological illness before the head injury; contraindications to MR scanning; and patients with obvious brain abnormality on conventional MRI.
Magnetic resonance imaging
The MRI scans were performed on a 3.0-T MRI scanner (Simens Verio, Erlangen, Germany) with a 32-channel head coil at our institution. A three-dimensional MPRAGE T1 sequence was used for measuring the GM volume with the following parameters: repetition time (TR)/echo time (TE)/inversion time = 1500/2.96/900 ms; flip angle = 9°; matrix = 256 × 256; field of view (FOV) = 230 mm; 120 axial slices; slice thickness = 1 mm. The DKI sequence was performed using a spin-echo echo-planar imaging diffusion sequence with 30 different diffusion-encoding directions. For each direction, six b-values (b = 0, 500, 1000, 1500, 2000, 2500 s/mm2) were acquired. Other imaging parameters were as follows: TR/TE = 5400/109 ms; matrix = 128 × 128; FOV = 230 mm; 35 axial slices; slice thickness = 3 mm. The scan duration for DKI was 14 min 4 s.
Image analysis
The whole-brain volumetric analysis was done using FSL-VBM, which is a voxel-based morphometry analysis method (22) and is a part of the FMRIB Software Library (FSL, version 5.0.1, Oxford, UK, http://www.fmrib.ox.ac.uk/fsl/) (23). As for the subcortical nuclei volumetric analysis, the “FIRST” toolbox (24) in the FSL software package was used to extract the bilateral volumes of the following six subcortical nuclei: hippocampus; thalamus; amygdala; caudate nucleus; putamen; and globus pallidus.
For DKI image analysis, motion artifacts and eddy-current–induced distortion were first corrected using the affine alignment of each diffusion-weighted image to the b0 image by employing the FSL software toolbox. Subsequently, the Diffusion Kurtosis Estimator software (DKE) (https://www.nitrc.org/projects/dke/) was used to calculate the DKI parametric maps including MK, as well as mean diffusivity (MD) and fractional anisotropy (FA) (25). The whole-brain DKI parameters were analyzed using TBSS (26) in the FSL. As for the subcortical nuclei DKI analysis, subcortical nuclei masks were aligned to the diffusion space using Flirt tool with nine degrees of freedom in FSL software package (27,28). All subcortical nuclei masks in the diffusion space were checked visually for registration errors and no errors were observed. All the parameters were manifested in the ITK-SNAP software (www.itksnap.org) (29). All the image analyses were done by an experienced board-certified neuroradiologist (M.L.W.). Supplemental Figure 1 depicts an example of subcortical nuclei segmentation and masks in the diffusion space.
Neuropsychological assessment
All participants completed a battery of neuropsychological tests that are sensitive to detect TBI-induced cognitive impairment on the same day, at which the participants underwent MRI. This study focused on three cognitive domains: verbal memory; attention; and executive function. Verbal memory was assessed using the Huashan Auditory Verbal Learning Test (HAVLT) (30). The total learning score for each participant was measured on the basis of the total score of three learning free-recall trials of 12 words. Furthermore, the attention abilities were examined using Digit Span Forward. The executive function was assessed using Digit Span Backward. Another experienced board-certified neuroradiologist (M.M.Y) performed the neuropsychological assessments.
Statistical analysis
Statistical tests were performed by the Statistical Package for Social Sciences (IBM Corp., Armonk, NY, USA) software for Windows version 20.0. The demographics, neuropsychological tests results, and subcortical nuclei MRI parameters between the mTBI group and the control group were compared using the chi-square test for qualitative variables and the Student’s t-test for normal distribution quantitative variables or the Mann–Whitney U-test for non-normal distribution quantitative variables. Pearson or Spearman correlation analysis was performed to examine the relationship between subcortical nuclei volume, DKI parameters and cognitive function. P values < 0.05 indicated statistical significance. The whole brain-based analysis was performed using a voxel-wise general linear model (GLM) for between-group comparisons and correlation analysis (age and gender included as covariates), with permutation-based non-parametric testing (5000 permutations). Threshold-free cluster enhancement (TFCE) (31) was used to form significant clusters, corrected for the multiple comparisons by family-wise error (FWE) rate. The corrected threshold for significance was P < 0.05.
Results
Patient demographics and neuropsychological assessment
Patient demographics and specific mTBI characteristics are shown in Table 1. The major reason for mTBI is traffic accident (15/23, 65.22%). Nineteen patients had a Glasgow Coma Scale (GCS) score of 15, three had a score of 14, and one had a score of 13. MRI and neuropsychological tests were performed at a mean of 39 days after the injury. Five patients with mTBI had cerebral microbleeds and no patients with mTBI were diagnosed with diffuse axonal injury. Compared with the control participants, the patients with mTBI had no significant differences in age, sex, and education. However, patients with mTBI performed worse in the domains of verbal memory, attention, and executive function (P < 0.05).
Patients demographics and neuropsychological results.
Values are given as n (%), mean ± SD (using Student’s t-test), or median (IQR).
*P < 0.05.
†Chi-square test.
‡Mann–Whitney U-test.
HAVLT, Huashan Auditory Verbal Learning Test; IQR, interquartile range; SD, standard deviation; TBI, traumatic brain injury.
Whole brain-based analysis in brain DKI parameters and volume
In the whole brain-based volume analysis, no regional GM volume differences were observed between the mTBI and control groups (P > 0.05, FWE corrected). As for the whole brain-based DKI parameters analysis, patients with mTBI showed lower MK value in widespread WM regions (P < 0.05, FWE corrected) than control participants. Fig. 1 depicts the specific significant regions mainly including the intrahemispheric association fibers of bilateral superior longitudinal fasciculus (SLF), cingulum, right inferior longitudinal fasciculus (ILF), inferior fronto-occipital fasciculus (IFF) and uncinate fasciculus (UF), inter-hemispheric fibers of body of corpus callosum, projection fibers of corticospinal tract (CST), anterior thalamic radiation (ATR), and bilateral internal capsule and external capsule. However, no statistical significance of FA and MD was observed between the mTBI and control groups.

Axial images derived from TBSS results and rendered on T1-weighted MR images from the ch2bet atlas showing lower MK in widespread WM regions in the mTBI group than that in the control group (P < 0.05 FWE corrected for multiple comparisons). Red = decrease in patients with mTBI. FWE, family-wise error; MK, mean kurtosis; MR, magnetic resonance; mTBI, mild traumatic brain injury; TBSS, tract-based spatial statistics; WM, white matter.
Further correlational study revealed that attention deficit revealed by lower digit span forward score was significantly associated with lower MK value in multiple WM regions (P < 0.05, FWE corrected) in the mTBI group. Fig. 2 depicts the specific significant regions, mainly including: intrahemispheric association fibers of bilateral SLF, ILF, IFF, UF, and cingulum; inter-hemispheric fibers of body of corpus callosum; projection fibers of CST, ATR, bilateral internal capsule, external capsule, and middle cerebellar peduncle.

Axial images correlational analysis rendered on T1-weighted MR images from the ch2bet atlas showing positive correlations between MK values and digit span test forward scores in patients with mTBI (P < 0.05 FWE corrected for multiple comparisons). Red = positive correlation in patients with mTBI. FWE, family-wise error; MK, mean kurtosis; MR, magnetic resonance; mTBI, mild traumatic brain injury.
Abnormalities of DKI parameters and volume in subcortical nuclei
Similar to the whole brain-based volume analysis, no significant volume differences were observed in any subcortical nuclei between the mTBI and control groups (P > 0.05). Table 2 shows the DKI parameters of subcortical nulei in the mTBI group and the control group. The mTBI group had a lower MK value in the right hippocampus, left thalamus, left caudate, right putamen, and right pallidum than those in the control group (P < 0.05). The MD value of the right pallidum in the mTBI group was also higher than that in the control group (P < 0.05). In addition, we observed no statistical difference in the FA value of subcortical nuclei between the mTBI group and the control group.
Statistical analysis of subcortical nuclei DKI parameters between mTBI and control group.
Values are given as mean ± SD (Student’s t-test).
*P < 0.05.
L, left; MD, mean diffusivity; MK, mean kurtosis; R, right; TBI, traumatic brain injury.
Further correlational study revealed that the MK values of the left amygdala (r = 0.625, P = 0.010), right thalamus (r = 0.634, P = 0.008), right putamen (r = 0.582, P = 0.018), right hippocampus (r= 0.576, P = 0.020), and right amygdala (r = 0.503, P = 0.047) were positively correlated with the attention ability. The MD value of the left thalamus was negatively correlated with the verbal memory (r = –0.678, P = 0.022). The FA value of the left amygdala was positively correlated with the verbal memory (r = 0.611, P = 0.016) (Fig. 3).

Correlational analysis between DKI parameters and neuropsychological scale scores. (a–e) Positive correlation between the MK of several subcortical nuclei with the attention ability revealed by the Digit Span Forward Score (P < 0.05). (f) Negative correlation between MD of the left thalamus with verbal memory revealed by Auditory Verbal Learning Test Score (P < 0.05). (g) Positive correlation between FA of the left amygdala with verbal memory revealed by Auditory Verbal Learning Test Score (P < 0.05). DKI, diffusion kurtosis imaging; FA, fractional anisotropy; MD, mean diffusivity; MK, mean kurtosis.
Discussion
The present study aimed to combine DKI and volumetric analysis to quantify the brain GM and WM microscopic structural changes in patients with sub-acute mTBI, and, if present, to correlate these changes with the neuropsychological test performance. Consistent with our hypothesis, DKI was sensitive to find brain microstructural changes revealing lower MK value in widespread WM regions and several subcortical nuclei in sub-acute mTBI, and were statistically associated with cognitive function changes. However, the whole brain-based analysis revealed no significant difference in the GM volume, FA, and MD value between the mTBI group and the control group.
Previous studies have revealed that the loss of GM volume was associated with cognitive impairment in patients with mTBI (16–18). However, all the present studies mainly focused on chronic mTBI. Whether the loss of GM volume begins to appear in an earlier sub-acute stage is unknown. In our study, we observed no significant changes in GM volume including subcortical nuclei in patients with sub-acute mTBI. We speculated that the neuron loss might not be obvious in the sub-acute stage. To date, the “atrophy curve,” that is, the relationship between brain volume and injury time is still unknown (32). Our study indicates that the loss of GM volume would not appear at the sub-acute stage (from two weeks to two months) after mTBI.
As for the DKI parameters, our study found lower MK value in widespread WM regions and several subcortical nuclei in patients with sub-acute mTBI. The study by Grossman et al. (13) also revealed lower MK values in individuals with mTBI. However, the study by Davenport et al. (33) reported higher MK values in high school football players. In a recent study of individuals with sport-related concussion, the concussed group showed significantly higher axial kurtosis than the control group (34). The contradictory results are most likely caused by the difference in injury mechanism, as football may only cause subtle injury to the brain. Furthermore, there was no statistical significance in FA and MD between the mTBI group and control group in our study. As most of the individuals with mTBI (19/23) had a GCS of 15 in our study, the impact of mTBI on the brain could be relatively weak, and thus, the differences both in FA and MD would be insignificant. The study by Hart et al. (35) also found no statistical differences in DTI-derived FA and MD values using TBSS between control and university amateur boxers. Thus, the present study suggested the MK parameter was more sensitive than DTI parameters including FA and MD parameters in detecting brain microstructural changes after mTBI.
The underlying brain pathological changes of DKI parameters are not yet completely understood. TBI is a heterogeneous and complex pathological process (36), including the following: (i) cell death and inflammation; (ii) cell proliferation for tissue replacement; and (iii) tissue remodeling (37). Reasonably, the DKI parameter would reflect an integration of multiple events after the injury. Our previous animal study revealed that low MK value in the TBI rat brain was associated with neuron loss (38). The lower MK value might result at least partially from the contribution of degenerative processes like neuron loss and changes in axonal and myelin density (39), leading to the decrement of diffusional heterogeneity.
Furthermore, the present study showed lower MK values in multiple WM regions and several subcortical nuclei and correlated with the attention deficit. This was reasonable as the lower MK value in the mTBI group might represent neuron loss and demyelinating, which could further cause the disruption of the brain attention circuit including the basal ganglia, amygdala, SLF, and so on (40). Therefore, our study suggested the MK value could be used to assess cognitive impairment in mTBI.
The present study has some limitations. First, the sample size in the study is relatively small, and future studies with larger sample sizes are needed to replicate our results. Second, as the baseline volume and DKI data of subcortical nuclei before mTBI were not obtained, we cannot conclude that the DKI parameter differences were attributable to the injury. The DKI parameter differences might exist before mTBI occurred and even predispose an individual to have an mTBI. Third, this study was a cross-sectional study, and thus we cannot determine the exact order of causal effects between neuroimaging parameters changes and cognitive impairment. To determine their order of causal effect, future large longitudinal studies will be required to better assess the relationships between neuroimaging changes and cognitive impairment after mTBI. Fourth, the MD and FA values were derived from the DKI model rather than the DTI model. Comparison with previous studies using DTI-derived MD and FA values could not be made as this could cause confusion in interpretation.
In conclusion, the present study found that DKI was sensitive to detect brain microstructural changes in patients with sub-acute mTBI, showing lower MK value in widespread WM regions and several subcortical nuclei, which were statistically associated with cognitive deficits. No loss of GM volume was observed in patients with sub-acute mTBI.
Supplemental Material
sj-pdf-1-acr-10.1177_0284185121998317 - Supplemental material for Cognitive impairment in mild traumatic brain injury: a diffusion kurtosis imaging and volumetric study
Supplemental material, sj-pdf-1-acr-10.1177_0284185121998317 for Cognitive impairment in mild traumatic brain injury: a diffusion kurtosis imaging and volumetric study by Ming-Liang Wang, Xiao-Er Wei, Meng-Meng Yu and Wen-Bin Li in Acta Radiologica
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received the following financial support for the research, authorship, and/or publication of this article: This study is supported by National Natural Science Foundation of China (Nos. 81901727, 81271540, 81301213), Natural Science Foundation of Xinjiang Province (No. 2016D01C083), 2018 Youth Medical Talents – Medical Imaging Practitioner Program, and Shanghai key discipline of medical imaging (No. 2017ZZ02005).
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
