Abstract
We explored the prognostic value of diffusion tensor imaging (DTI) parameters of selected white matter (WM) tracts in predicting neuropsychological outcome, both at baseline and 6 months later, among well-characterized patients diagnosed with mild traumatic brain injury (mTBI). Sixty-one patients with mTBI (mean age=27.08; standard deviation [SD], 8.55) underwent scanning at an average of 10 h (SD, 4.26) post-trauma along with assessment of their neuropsychological performance at an average of 4.35 h (SD, 7.08) upon full Glasgow Coma Scale recovery. Results were then compared to 19 healthy control participants (mean age=29.05; SD, 5.84), both in the acute stage and 6 months post-trauma. DTI and neuropsychological measures between acute and chronic phases were compared, and significant differences emerged. Specifically, chronic-phase fractional anisotropy and radial diffusivity values showed significant group differences in the corona radiata, anterior limb of internal capsule, cingulum, superior longitudinal fasciculus, optic radiation, and genu of corpus callosum. Findings also demonstrated associations between DTI indices and neuropsychological outcome across two time points. Our results provide new evidence for the use of DTI as an imaging biomarker and indicator of WM damage occurring in the context of mTBI, and they underscore the dynamic nature of brain injury and possible biological basis of chronic neurocognitive alterations.
Introduction
M
Mild neurotrauma is associated with traumatic axonal injury (TAI), which is described as a progressive event gradually evolving from focal axonal alteration to delayed axonal disconnection. 3 Importantly, TAI is thought to represent one of the more common injuries observed in the aftermath of mTBI. 2 These subtle alterations of brain tracts or fiber pathways have been visualized using diffusion tensor imaging (DTI), which enables better visualization of the extent of early microstructural changes post-mTBI. 2,4 –6 A variety of metrics can be generated through DTI scans, including fractional anisotrophy (FA), which is a per-voxel indication of the directionality of underlying water diffusion. FA values range from 0 to 1, where FA=0 would indicate nondirectional diffusion (completely isotropic) and FA=1 would indicate a single direction of diffusion where the water molecules are restricted to diffusion only along a single axis (completely anisotropic). 2,7 Reduced FA in the white matter (WM) is believed to reflect a loss of integrity, indicating possible damage to myelin or the axon membrane, reduced axonal packing, and/or decreased axonal coherence. 2,8 Mean diffusivity (MD), on the other hand, describes the per-voxel average magnitude of water diffusion, regardless of diffusion direction. Differences in MD are thought to reflect overall restrictions to the movement of water diffusion, examples being the variations within the intra- and extracellular space. 9 Radial diffusivity (RD) is defined as the diffusion of water perpendicular to WM fibers, 10 which increases in response to demyelination 11 and dysmyelination. 12 Changes in axonal diameter or density can also influence changes observed in RD. 13 FA and MD values are usually inversely correlated with one another, given that myelination, which enforces directionality (thus increasing FA), also represents a restriction to overall movement (thus lowering MD). 14,15 The demyelinating changes, as evinced by the changing RD, however, are not expected to occur in the first week post-mTBI, despite the presence of axonal swelling and synaptic disruption. 16
In mTBI, widespread changes in FA are frequently observed, especially in the frontal, mid-line, and temporal regions. Studies have shown that these changes can be detected as early as a few days to weeks after neurotrauma, as well as months or even a year after the initial insult. 14,15,17 The shearing forces of trauma can breach the vascular permeability of vessels (hence rupturing them), sever fibers, and lyse cells, which deregulates the normal homeostasis of the blood–brain barrier, and usually manifests as cerebral vasogenic and cytotoxic edema. 18 Evaluation of DTI indices enables the differentiation of these edemas, 19 which is crucial to predicting long-term neurological outcomes, including neuropsychological performance (NP). Vasogenic edema commonly observed in mTBI is characterized by reduced FA, increased MD, and RD and is considered reversible. In contrast, cytotoxic edema, characterized by an increased FA, reduced MD, and RD, is considered irreversible and therefore confers a poor prognosis. 2,19,20
To date, there are few studies that have been conducted longitudinally with acute mTBI samples to help elucidate the evolution of these DTI-based changes in mTBI over time. 21 –25 Unfortunately, existing studies have yielded generally equivocal findings. 2,15,26 –29 For example, some studies have reported decreased integrity of several tracts at different time intervals (acute and chronic), 30 –35 although others report elevated FA and reduced MD in the acute stage. 16,21,36 Given that these findings may be inconsistent owing to various methodological differences, including patient recruitment, imaging protocol differences, varying intervals studied, sample-size differences, and heterogeneity of injury severities, we aimed to clarify these longitudinal DTI changes using a whole-brain WM measurement strategy with tract-based spatial statistics (TBSS) 37 in mTBI and control groups. From the TBSS WM skeleton of comparison across subjects, we identified significant tract changes and correlated these regions with neuropsychological performances, both at admission and 6 months postinjury, in patients with mTBI. We also examined the relationship between anatomical correlates of tracts and cognition in an effort to improve prognostic values of DTI parameters in mTBI care.
Methods
Participants
Sixty-one patients with mTBI who presented to the emergency department (ED) of University of Malaya Medical Center (Kuala Lumpur, Malaysia) for a consecutive 11-month period between April 1, 2013 and March 1, 2014 were prospectively recruited for this study. Patients were selected based on the inclusion and exclusion criteria as presented in Figure 1. For the purposes of this study, mTBI was defined as acute TBI, consisting of nonpenetrating head trauma resulting in one or more of the following: confusion/disorientation; loss of consciousness (LOC) less than 30 min; post-traumatic amnesia (PTA; less than 24 h in duration) and/or transient focal neurological signs or seizures; and Glasgow Coma Scale (GCS) of 13–15 upon acute clinical evaluation. The flow of the study is presented in Figure 2. Nineteen healthy age-matched control participants were also recruited for this study.

Inclusion and exclusion criteria for patients with mTBI. mTBI, mild traumatic brain injury; GCS, Glasgow Coma Scale; CNS, central nervous system.

Complete flowchart of the study protocols, patient and healthy control recruitment, and neuroimaging and neuropsychological assessments. TBI, traumatic brain injury; CT, computed tomography; mTBI, mild traumatic brain injury; MRI, magnetic resonance imaging; FLAIR, fluid-attentuated inversion recovery; DTI, diffusion tensor imaging; DWI, diffusion-weighted imaging; FSPGR, fast spoiled gradient echo; S-NAB, Neuropsychological Assessment Battery-Screening.
All subjects meeting criteria for the study underwent computed tomography (CT) scans of the brain in the ED using a Siemens Somatom Sensation 16 CT scanner (Siemens AG, Berlin, Germany). Cross-sectional images of the brain were obtained craniocaudally from the base of skull to vertex. Scan parameters used were 120 kVp, 300 mAs, and collimation of 16×0.75 mm with standard brain and bone windowing. A neuroradiologist (N.R.) and a neurosurgeon (V.N.) who were blinded to the clinical diagnosis independently evaluated the images for each patient, and only patients who were deemed not requiring surgical intervention were included in this study. All subjects gave informed consent as required by the institutional research ethics committee and the hospital ethics committee (UM/EC Ref: 947.15)
Study protocols
Magnetic resonance imaging (MRI) and neuropsychological assessments were performed at admission and repeated again at 6 months post-trauma (Fig. 2). Healthy control participants were subjected to the same protocols as patients upon admission (i.e., MRI and neuropsychological assessment).
Magnetic resonance imaging data acquisition
All consented subjects were imaged on a 3T MRI scanner (Signa HDx; General Electric, Harvey, IL) using an eight-channel head coil. The imaging protocol included: 1) axial T1-weighted three-dimensional fast spoiled gradient echo (FSPGR; repetition time [TR]=minimum 6.7 ms; echo time [TE]=minimum 1.9 ms; field of view [FOV]=31 mm; matrix=256×256; slice thickness=1.2 mm; slice overlap=0.6 mm with image scan time of 3 min 48 sec); 2) axial T2-weighted fast spin echo; TR=4240 ms; TE=102 ms; FOV=24 mm; matrix=512×384; thickness=5 mm; spacing=1.5 mm; and image scan time of 2 min 30 sec; and 3) coronal gradient echo (TR=655 ms; TE=20 ms; flip angle=15 degrees; bandwidth=31.25; FOV=24 cm; matrix=320×256; thickness=5.0 m; spacing=1.5 mm; and image scan time of 2 min 7 sec. The DTI sequence was obtained using these parameters: TR=13,000 ms; TE=81.2 ms; FOV=24 mm; matrix=128×128; slice thickness=3.0 mm; 32 directions; diffusion-weighted factor, b=700 s/mm2; and image scan time of 7 min 22 sec.
Magnetic resonance imaging analysis
Tract-based spatial statistics
Voxel-wise statistical analysis of the diffusion-weighted data was carried out using TBSS,
37
part of the FSL (v5.0.6; University of Oxford, Oxford UK) software package. Initial preprocessing involved corrections for head movement and eddy currents, brain tissue extraction, and fitting of the diffusion tensor model. These were carried out using the FSL eddy_correct, bet, and dtifit tools, respectively. The standard TBSS analysis workflow was followed (V 1.2;
The voxel-wise statistical analysis was carried out using the FSL randomize tool. A two-group unpaired t-test design was used for comparing admission scans of control and mTBI subjects. A paired t-test design was used for comparing admission and follow-up scans of mTBI subjects. These statistical analyses were carried out separately for FA, MD, and RD values. In all cases, cluster-based thresholding was used, and 0.05 was adopted as the threshold for significant clusters.
Region of interest
In addition, we obtained mean FA, MD, and RD for all tracts identified on TBSS using region of interest (ROI) analysis. The image-processing pipeline consisted of preprocessing, image registration, and analysis, utilizing the FSL (v5.0.6; University of Oxford) and AFNI (v2011_12_21_1014; National Institute of Mental Health, Bethesda, MD) software packages. Initial preprocessing involved corrections for head movement and eddy currents, brain tissue extraction, and fitting of the diffusion tensor model. These were carried out using the FSL eddy_correct, bet, and dtifit tools, respectively. For image registration, the FSL tool fnirt was used to carry out nonlinear spatial registration of each subject to the FMRIB58_FA standard-space image, using the built-in FA_2_FMRIB58_1mm config file. This was assumed to be spatially compatible with the International Consortium of Brain Mapping (ICBM) DTI-81 atlas. Postregistration, composite axial slice images of the underlying FA and ICBM atlas tract outlines were generated at z=42 and z=82 for each subject. An experienced neuroradiologist (N.R.) verified the overall registration quality of these images. Finally, the AFNI 3dROIstats tool was used to map the predefined ROIs to each individual subject and calculate the median FA, MD, and RD values for each tract.
Neuropsychological assessment
All subjects underwent cognitive assessment using the Screening Module of the Neuropsychological Assessment Battery (S-NAB Form 1), which was performed by the neuropsychologist once the patient had recovered to a GCS score of 15, which occurred with an average turnaround time of 4.35 h (SD, 7.08) between time of trauma and full GCS recovery. The S-NAB comprises a comprehensive set of neuropsychological tests (refer to Table 1), with demographically corrected norms for adults between the ages of 18 and 97 years, assessing orientation and five cognitive domains (i.e., attention, memory, language, and visuospatial and executive functions). This battery consists of 12 individual tests across the five domains aforementioned. From these 12 tests, a total of 16 T scores are derived, 14 of which contribute toward five separate Screening Index (domain-specific) scores and one Total Screening Index score. 38 The same subtests were repeated at 6 months by the same neuropsychologist using the S-NAB Form 2 in order to minimize practice effects.
S-NAB, Neuropsychological Assessment Battery-Screening.
Statistical analysis
An independent-samples t-test was used to establish whether mean values of FA, MD, and RD of selected WM tracts were significantly different between healthy control and mTBI groups during the acute phase (or baseline exam). The same test was also used to investigate whether patients (at admission) performed differently from healthy control participants on the neuropsychological assessment. A paired t-test was used thereafter to ascertain how the WM tracts had changed over time. TBSS skelotonized image of the significant changes observed over time were processed to better visualize any signifcant changes. Spearman's rho correlation coefficient was used to examine the association between WM ROIs (nine selected tracts) and NP over the different phases. Last, a simple frequency analysis of neuropsychological performance at 6 months post-trauma was performed to determine the types of changes observed in neuropsychological status longitudinally (improved, unchanged, or worsened).
Results
Demographic and clinical data
The demographic characteristics of study patients and healthy controls are presented in Table 2. The TBI group ranged in age between 18 and 53 years (mean, 27.08; standard deviation [SD], 8.55), was predominantly male (88.5%), and had a mean age of 29.05 (SD, 5.84). There was no significant difference in the mean age of healthy controls versus the TBI patient group (t (78)=0.937; p=0.256). However, the mean years of education differed significantly between groups (TBI, 11.52 [SD, 1.94] vs. controls: 15.95 [SD, 2.01] years). Admission GCS scores of patients were within the range of 13–15 (mean=14.44; SD, 0.74). Seventy-seven percent of patients had experienced transient LOC whereas 73.8% reported PTA. Average time to scan of patients during the acute phase was 10.01 h (SD, 4.26), whereas the chronic phase scan was performed at 6 months postinjury (mean=6.05; SD, 0.12). Of the 61 patients, 33 underwent the chronic phase repeat scans (3 patients were later dropped during the postscan processing because their scans had major artifacts in echo planar imaging). All healthy controls had no significant neurological findings.
GCS, Glasgow Coma Scale; LOC, loss of consciousness; PTA, post-traumatic amnesia; GOSE, Glasgow Outcome Score Extended; MVAs, motor vehicle accidents; TBI, traumatic brain injury; SD, standard deviation; n/a, not applicable.
Diffusion metrics, intergroup differences, and intragroup changes over time
Table 3 presents mean FA, MD, and RD values of both patients and healthy control participants during the acute phase. At baseline, the mTBI group showed significantly lower splenium FA (t (78)=2.196; p=0.04), when compared to the control group, whereas the MDa value was significantly higher in the mTBI versus control group in the posterior limb of the internal capsule (PLIC; (t (78)=−2.179; p=0.03); cingulum (CG; t (78)=−2.29; p=0.02); optic radiation (OR; t (78)=−3.176; p=0.002); and splenium (SCC; t (78)=−2.514; p=0.02). Group differences for the RDa values were observed only in the OR (t (78)=2.637; p=0.02) and SCC (t (78)=2.519; p=0.02). See Figure 3 for the TBSS skeletonized image reflecting the significant area of differences between the groups.

Tract-based spatial statistics (TBSS) showing composite of fractional anisotropy/mean diffusivity
Variances in the group were similar for all comparisons.
DTI, diffusion tensor imaging; FA, fractional anisotropy; MD, medial diffusivity; RD, radial diffusivity; mTBI, mild traumatic brain injury; SD, standard deviation; CTRL, control; TBI, traumatic brain injury.
The paired t-test results (FA, MD, and RD) of patients' WM pathway changes, as observed at 6 months (chronic phase) against the baseline (acute phase) values, are presented in Table 4a, 4b, and 4c, along with their effect sizes (Cohen's d with <0.2=small effect, >0.2 <0.8=moderate effect, and >0.8=large effect). Various differences were observed in all DTI parameters, with greater alterations noted in FA values across regions. Specifically, FA of the corona radiata (CR; t (28)=3.497; p=0.002; d=0.661), anterior limb of internal capsule (ALIC; t (28)=2.582; p=0.016; d=0.488), CG (t(28)=2.973; p=0.006; d=0.562), superior longitudinal fasciculus (SLF; t (28)=2.404; p=0.024; d=0.454), OR (t (28)=2.643; p=0.014; d=0.499), and the genu of corpus callosum (GCC; t (28)=2.732; p=0.011; d=0.516) were significantly lower in the TBI group. Almost all of the MD values across the phases showed no significant changes, except the CG (t (28)=3.189; p=0.004; d=0.603), which was significantly higher in the TBI group. Although the changes in PLIC across the phases were found to be statistically nonsignificant (t (28)=1.494; p=0.15), the effect size, however, suggests a moderate level of change (d=0.282). Finally, there were no significant differences in RD values across time points. However, Cohen's d effect-size calculation indicated moderate effect sizes for changes in the RD of CR (t (28)=1.582; p=0.126; d=0.299), CG (t (28)=1.54; p=0.136; d=0.291), OR (t (28)=1.87; p=0.073; d=0.353), and GCC (t(28)=1.975; p=0.059; d=0.373).
FA, fractional anisotropy; SD, standard deviation.
MD, medial diffusivity; SD, standard deviation.
RD, radial diffusivity; SD, standard deviation.
All the above-stated significant differences were reflected in the TBSS images of the skeletonized tracts as well (see Fig. 4).

Tract-based spatial statistics (TBSS) showing composite of fractional anisotropy/mean diffusivity
Neuropsychological performance
Table 5 presents the mean interpretive categories score comparison for the domain-specific NP among mTBI and healthy control groups. During the acute phase, patients with mTBI performed poorly across all domains, in comparison to the healthy control group. The independent-samples t-tests of both groups and their NPs indicated that the mTBI group was significantly poorer (all p values, <0.001) on all but one of the neuropsychological domains (visuospatial functions, t (78)=0.055; p=0.956). Meanwhile, during the chronic phase, patients with mTBI continued to perform poorly across most domains, in exception of visuospatial functions, where the patients outperformed healthy controls (t (48)=−2.373; p=0.021).
Interpretive category score legend: 1=severely impaired; 2=severe to moderately impaired; 3=moderately impaired; 4=mildly to moderately impaired; 5=mildly impaired; 6=below average; 7=average; 8=above average; 9=superior; and 10=very superior.
S-NAB, Neuropsychological Assessment Battery-Screening; SD, standard deviation; CTRL, control; TBI, traumatic brain injury.
Associations between diffusion tensor imaging parameters and neuropsychological performance
Longitudinal analysis of WM tract changes against the NP at different intervals among a subset of study patients (n=30) are presented in Table 6. Most of the observed FA associations with neurocognitive status (acute and chronic) represented negative associations. FAa negatively correlated with the following acute NPa: attention versus CR (r=−0.429; p<0.05); language versus CR (r=−0.375; p<0.05); language versus SLF (r=−0.557; p<0.01); and language versus GCC (r=−0.443; p<0.05). Acute FAs (FAa) were also negatively correlated with NPc, as follows: language versus middle cerebelar peduncle (MCP; rs =−0.440; p<0.05); attention versus CR (rs =−0.441; p<0.05); language versus CR (rs =−0.415; p<0.05); attention versus SLF (rs =−0.417; p<0.05); language versus SLF (rs =−0.409; p<0.05); spatial versus CG (rs =−0.489; p<0.05). The two positive correlations observed in the chronic phase were the spatial versus CG (rs =0.489; p<0.05) and spatial versus SCC (r=0.402; p<0.05). The significantly reduced FAc had only two associations with the NPc: language versus CR (r=−0.400; p<0.05) and spatial versus SCC (r=0.402; p<0.05)
Correlation is significant at p<0.05 (two-tailed).
Correlation is significant at p<0.01 (two-tailed).
DTI, diffusion tensor imaging; FA, fractional anisotropy; MD, medial diffusivity; RD, radial diffusivity.
MDa showed limited associations with the NPc. The following were the only MDa values associated with the NPc: attention versus SLF (rs =0.404; p<0.05); spatial versus CG (rs =−0.390; p<0.05); and spatial versus GCC (rs =–0.404; p<0.05). No associations between the MDc and NPc were observed.
The RD of seven of nine tracts studied was significantly associated with attention, language, spatial, and executive function in both phases of the study. Specifically, RDa was associated with three NPa scores: attention versus CR (r=0.485; p<0.05), attention versus SLF (r=0.487; p<0.05); and spatial versus OR (r=−0.378; p<0.05). Additional associations were subsequently observed between RDa and NPc, which were as follows: language versus MCP (r=0.398; p<0.05); language versus CR (rs =0.529; p<0.01); attention versus SLF (r=0.450; p<0.05); language versus SLF (r=0.491; p<0.05); language versus CG (r=0.423; p<0.05); spatial versus CG (r=0.626; p<0.01); and attention versus GCC (r=−0.378; p<0.05). Increased RDc values were also associated with some domains of the NPc, including language versus MCP (rs =0.438; p<0.05), executive function versus CR (rs =0.389; p<0.05); executive function versus CG (rs =−0.404; p<0.05); and executive function versus SCC (rs =−0.391; p<0.05). No association was found between the domains of memory and any of the WM tracts investigated in this study. We also found no associations between the DTI parameters (FA, MD, and RD) and any domains of the NP in the healthy control group (kindly refer to Table 6b).
DTI, diffusion tensor imaging; FA, fractional anisotropy; MD, medial diffusivity; RD, radial diffusivity.
Neuropsychological outcomes
Table 7 presents the neuropsychological profiles of the subset of patients (n=30) who had completed neuropsychological evaluations, as well as imaging, both at admission and at month 6 follow-up. Within the domain of attention, 56.7% of the patients remained impaired at 6 months post-trauma, 3.33% worsened, and the remaining improved or remained unaffected. Within the domain of language function, 63.3% of patients remained impaired, 26.7% improved, 6.67% remained unaffected, and 3.33% worsened. Approximately 33.3% of patients remained impaired within the domain of memory 6 months post-trauma with essentially equivocal changes (26.7% worsened, 26.7% remained unaffected, and 13.3% with improved memory status). The majority of patients (53.3%) remained unaffected for spatial function, although 23.3% showed signs of delayed impairments in the chronic phase. A total of 70.0% of the subset patient group remained impaired with respect to executive functioning, with only 20.0% of these patients showing signs of improvement or recovery after the 6-month period post-trauma.
Interpretive categories: impaired=1; not impaired=2.
Discussion
We examined the relationship between microstructural changes and neuropsychological functioning that takes place in the immediate aftermath of an mTBI as well as 6 months post-TBI. Specifically, the initial neuropsychological assessment was completed, on average, 4.35 h after full GCS recovery, and the neuroimaging procedure was completed within an average of 10 h post-trauma in order to identify the structural changes at a very early stage before they were subject to confounds, such as changes specific to environmental and recovery parameters. Subsequently, both a repeat imaging and neuropsychological assessment were performed at an average of 6 months post-trauma to characterize any pertinent changes over time.
Results showed associations between DTI and neuropsychological indices at both acute and chronic phases. Specifically, in the acute phase, the TBI group showed significantly poorer WM integrity across several tracts, including the splenium (SCC), PLIC, CG, and OR. Several other WM tracts beyond those initially affected also showed reduced FA at 6 months post-injury (i.e., CR, ALIC, SLF, and GCC). Interestingly, moderate effect sizes for changes in RD were observed from baseline to follow-up in the following regions: CR, CG, and OR. Finally, with respect to neuropsychological functioning, patients with mTBI performed more poorly across all domains in the acute stage, when compared to healthy control participants. The majority of patients remained impaired with respect to executive functioning, and a sizable portion of the sample worsened over time on tasks of attention and language. Moreover, longitudinal analysis of WM tract changes, as they relate to cognition, showed several associations, indicating that reduced WM integrity was associated with specific neuropsychological decrements. Indeed, DTI indices in both acute and chronic stages related to cognition at both time points, and RD values in seven of the nine tracts studied was significantly associated with attention, language, spatial, and executive function in both phases of the study. Interestingly, no association was found between the domains of memory and any of the WM tracts investigated in this study.
Our findings of reduced FA, coupled with increased MD and RD in the acute phase, is considered indicative of vasogenic brain edema, which refers to the release of intracellular proteins into the brain parenchyma, also known as extracellular edema. 18 An increased FAa, reduced MDa, and RDa in the acute phase, on the other hand, may be indicative of cytotoxic edema, 2,36 also known as intracellular edema. This type of cerebral edema in the early course of injury has been implicated in poor outcome, 2,18 which could explain frequently observed neuropsychological impairment months or even years after the initial trauma event. The increased FA, MD, and rather an unchanged RD in the acute phase of this study is most likely indicative of reactive astrogliosis, which is a process where fibrous astrocytes migrate to the site of injury, locally increasing the density of the cells 39 and diffusivity of the affected tissue. 39,40
Vasogenic edema occurs in the acute phase, whereas demyelination occurs later. 16 As such, even though the DTI markers for these pathogenic processes are the same, the timeframe of imaging changes rather indicate an ongoing vasogenic edema. These changes in DTI parameters were observed in the SCC, PLIC, CG, and OR (with statistically significant changes in most parameters). Other tracts, including the ALIC, SLF, and GCC, showed similar trends, although they did not reach statistical significance. Results showing minimally reduced FAa in the context of significantly increased MD of the SLF and GCC in the patient group, and its negative association with language function, are similar to those reported by Ingelese and colleagues 33 and Arfanakis and colleagues. 24 Cognitive impairment observed acutely in this group of patients coupled with the changing DTI metrics are likely immediate signs reflecting an ongoing edematous process, which is not observed by conventional CT or MRI.
In contrast, changes observed in CRa and MCPa were more suggestive of reactive astrogliosis and dovetails with findings recently reported by Croall and colleagues. 41 Although some of these changes did not reach statistical significance, the effect size of 0.299 for the RD changes in CR, coupled with the effect size of 0.661 for the FA changes of CR, for instance, does provide some moderate-to-large evidence of significant change in this tract. Higher FA value of the CR also negatively correlated with attention and language functions acutely. Its corresponding positive association with RD at baseline best explains the influence of reactive astrogliosis on cognition in the acute stage, as previously stated. 40 –42 The importance of the above finding cannot be understated given that specific types of cerebral edema and gliosis in the acute phase have been implicated to negatively influence long-term neuropsychological performance in patients with mTBI. 43 –46
We also found that acute-phase DTI parameters in selected WM tracts were significantly associated with chronic domain-specific cognitive deficits. This includes those of the MCP, CR, SLF, CG, GCC, and SCC and their association with specific neurocognitive functions chronically, similar to findings of other studies. 39 –42,47,48 Miles and colleagues, for instance, noted that decreased FA values and increased MD values of the GCC, SCC, and PLIC during the acute phase were significantly related to executive dysfunction at 6 months follow-up. 2,47 Kumar and colleagues and Matshushita and colleagues reported similar findings involving the GCC, SCC, and NPc in their studies of mild and moderate head injury. 48,49 The mostly positive association between the unchanged or raised RDa of the MCP, CR, SLF, and CG with specific chronic deficits in attention, language, and spatial function likely reflect the long-term effects of acute vasogenic edema and reactive astrogliosis, a well-known immune response in the immediate aftermath of central nervous system (CNS) injury. 41 It is important to note that cerebral edema and cascading gliosis usually occur concurrently and are not necessarily independent of each other.
Chronically, the continued alteration of DTI parameters, especially the FAc and RDc, implies significant changes in WM integrity, where the myelin sheath or the axonal membrane or both may have been permanently damaged. 19,20,36 This disruption may be irreversible, especially when a reduced FAc and elevated MDc and RDc are noted. 2,19 Conversely, a reduced FAc, unchanged MD, and subtly altered RDc would mean that the microstructural changes to the WM are independent of gross tissue loss. 43 –46 The subtly altered RDc in this context would imply minor fiber damage without gross tissue loss. 44 Such changes have also been implicated as a part of the dual effects of gliosis. 46,50,51
FA measured at follow-up was significantly reduced, in comparison to the acute phase, with mostly unchanged MDc and subtly increased RDc across six tracts (i.e., CR, ALIC, CG, SLF, OR, and GCC). Whereas the changes observed in the ALIC and OR are indicative of possible minor fiber damage and thus vasogenic edema, the change pattern observed in CR suggests irreversible consequences of astrogliosis. The mixed association between certain neurocognitive performances in the chronic stage (e.g., executive, spatial, and language function) and FA and RD of several tracts (i.e., MCP, CR, CG, and SCC), though possibly paradoxical, is not completely unexpected. Similar observations were found in a recent study with attribution to the multifaceted chronic phase findings. 41 This included neuronal network reorganisation, 52 spurting of axons with smaller calibres, 53 glial scarring, 46,50,51,54 effects of accumulated phosypholarated neurofilaments, 55 –59 and disruptive neurofibrillary tangles. 60 Taken together, our findings further strengthen evidence for physiogenic influences on the prolonged functional and neuropsychological sequelae in patients with mTBI. 61 Clearly, disruption of the connectomes (i.e., the connectivities between specific cortical and subcortical structures) by various neuropathological 62,63 events post-trauma and the recovery process over time requires further study.
Interestingly, a common finding in our study was that the more severely affected fibers were generally the long tract fibers (SLF, CR, and CG) and commisural fibers (SCC and GCC). These fibers have an increased susceptibility to injury 64 given their relatively long length and high membrane-to-cytoplasmic ratios. 65,66 The corpus callosum, for example, as noted by Aoki and colleagues, is the major fiber bundle that enables communication between the hemispheres and is topographically organized. 27,67 The CC in general is divided into the genu (GCC), the body, and the splenium (SCC). 27 The CC has been long recognized as a frequently injured region owing to shear strain, 68 and topographical location and external accelerational forces can seriously injure the fiber. 69,70 The posterior region of the CC, namely, the SCC, often noted to be more vulnerable to injury than the anterior part (GCC), 67,71 –74 was replicated in this study as well; and the SCC was the only tract that showed statistically significant changes in all three DTI parameters at the acute stage.
Our findings lend support to the theory that microstructural changes occur within hours after the initial insult and affect the integrity of the WM, therefore leading to the varied manifestation of neuropsychological impairments, in agreement with some of the earlier studies reported in the literature. 31,34,41,75,76 Strengths of our study include examination within a short timeframe from time of injury to imaging and neuropsychological testing within the acute phase, as well as consistent test-retest interval at 6-month follow-up. Although FA and MD have been shown to be very sensitive in detecting subtle WM changes that correlate with neuropsychological findings, one other parameter that has been included in this study, namely, RD, 77 provides further insight into the presumed nature of the microstructural changes 75,78,79 given that it is sensitive to myelin integrity 23,76,80 and also other CNS immune response, such as gliosis. 41 These parameters enabled us to detect the early deep WM microstructural changes, which might otherwise be confounded by ongoing secondary damage or recovery mechanisms. Additionally, the current study is unique in its longitudinal follow-up of DTI indices and cognitive functioning at 6 months post-trauma in a large cohort of well-characterized patients with relatively homogenous types of injuries. That said, findings may be limited owing to the significant heterogenity of the trauma observed in our sample as well as possible obscuring of certain differences owing to the nature of brain morphing and averaging, which is inherent in template-based analyses. Specifically, any deviation present in a small number of voxels within a region may be hidden by the averaging and therefore lead to low sensitivity to detect FA and MD changes. 81,82
Conclusion
DTI is a useful technique used to assess the integrity of WM tracts in patients who have sustained mTBI. This study found significant correlations between neuropsychological deficits and WM tract integrity both within hours of neurotrauma and again 6 months later. These associations are likely attributable to alterations in the integrity of connectomes between specific cortical areas and subcortical structures, which is not evident in conventional MRI or CT procedures. Specific WM changes, as observed through the DTI parameters, especially in the acute period, and the corresponding neuropsychological impairments among the patients in the chronic phase, highlights the need to introduce appropriate imaging techniques early in patient management protocols for early prognostication and rehabilitative intervention. Further, our findings underscore the dynamic nature of brain injury and possible biological basis of cognitive dysfunction in the context of the postconcussive syndrome.
Footnotes
Acknowledgments
This research is partially funded by a University Malaya Research Grant (UMRG; RG008C-13HTM) and a High Impact Research Grant of University of Malaya (HIR-UM.C/625/1/HIR-MOHE/12).
Author Disclosure Statement
No competing financial interests exist.
