Original Research ARTICLE
Investigating Microstructural Abnormalities and Neurocognition in Sub-Acute and Chronic Traumatic Brain Injury Patients with Normal-Appearing White Matter: A Preliminary Diffusion Tensor Imaging Study
- 1Department of Neurosurgery, St. Michael’s Hospital, Toronto, ON, Canada
- 2Department of Neurology, University of California at San Francisco, San Francisco, CA, USA
- 3Department of Brain and Behavioral Sciences, University of Pavia, Italy
- 4Department of Radiology, St. Michael’s Hospital, Toronto, ON, Canada
- 5Sunnybrook Research Institute, University of Toronto, Toronto, ON, Canada
- 6Department of Medical Imaging, St. Michael’s Hospital, University of Toronto, Toronto, ON, Canada
- 7Department of Medical Imaging and Neurosurgery at the University of Toronto, Toronto, ON, Canada
- 8Faculty of Medicine, Department of Surgery, Division of Neurosurgery, University of Toronto, Toronto, ON, Canada
For a significant percentage of subjects, with chronic traumatic brain injury (TBI), who report persisting cognitive impairment and functional loss, the diagnosis is often impeded by the fact that routine neuroimaging often does not reveal any abnormalities. In this paper, we used diffusion tensor imaging (DTI) to investigate the apparently normal white matter (as assessed by routine magnetic resonance imaging) in the brains of 19 subjects with sub-acute (9) and chronic (10) TBI. We also assessed memory, executive function, and visual-motor coordination in these subjects. Using a voxel-wise approach, we investigated if parameters of diffusion were significantly different between TBI subjects and 17 healthy controls (HC), who were demographically matched to the TBI group. We also investigated if changes in DTI parameters were associated with neuropsychological performance in either group. Our results indicate significantly increased mean and axial diffusivity (MD and AD, respectively) values in widespread brain locations in TBI subjects, while controlling for age, sex, and time since injury. HC performed significantly better than the TBI subjects on tests of memory and executive function, indicating the persisting functional loss in chronic TBI. We found no correlation between diffusion parameters and performance on test of executive function in either group. We found negative correlation between FA and composite memory scores, and positive correlation between RD and visuomotor coordination test scores, in various tracts in both groups. Our study suggests that changes in MD and AD can indicate persisting micro-structure abnormalities in normal-appearing white matter in the brains of subjects with chronic TBI. Our results also suggest that FA in major white matter tracts is correlated with memory in health and in disease, alike; larger and longitudinal studies are needed to discern potential differences in these correlations in the two groups.
Traumatic brain injury (TBI) often results in cognitive and functional deficits, impairs daily life functioning, and degrades the quality of life of the injured person significantly in the acute stage (1). Furthermore, the cognitive and functional deficits may persist and yet remain untreated in the chronic stage, particularly when routine diagnostic imaging [e.g., structural magnetic resonance imaging (MRI), and computed tomography (CT)] does not reveal any gross structural abnormalities, indicative of the cognitive or functional loss, in the brain of the injured patients. A large literature indicates that these persisting deficits are linked to microstructural damage in the brain (2). Diffusion tensor imaging (DTI) (3), which uses the diffusion of water molecules to assess microstructure in both pathological and normal brain tissue, has frequently been used in the past decade to investigate microstructural damage in patients with TBI (4–6) [also see the review by Hulkower et al. (7)]. The DTI studies of TBI often investigate the change in one or more of the classic diffusion parameters—axial diffusivity (AD), radial diffusivity (RD), mean diffusivity (MD), and fractional anisotropy (FA)—using either region of interest (ROI) (8, 9) or whole brain voxel-wise (10, 11) approaches. Most of these studies have focused on FA alone (10) or FA and MD (12, 13), the two parameters that describe the overall or average diffusion trends. These studies suggest a general trend of decrease in directionality of diffusion as measured by FA (7, 10, 14, 15) and an increase in the overall diffusion as measured by MD (7, 16, 17) in subjects with chronic TBI and associate these changes with diffuse axonal injury (DAI).
The DAI, an example of microstructural disruption resulting from TBI, is thought to be one of the major causes of neurocognitive deficits after TBI (18, 19). However, given the complexity of the various cognitive and executive function networks in the brain, it is important to investigate the role of DAI in neurocognitive deficits after TBI in situations where no other confounding structural anomalies (e.g., focal lesions indicative of the cognitive or functional loss) are present and by adopting a whole-brain analysis approach. The DTI is an excellent technique to investigate the DAI in TBI, however, with notable exceptions (20, 21), most of the DTI studies investigating the correlation between neurocognitive function and DTI parameters did not exclude subjects with white matter abnormalities as assessed with routine diagnostic imaging (22–25). Furthermore, many of these DTI studies investigating neurocognitive function after TBI are restricted only to FA and MD and only to TBI patients, excluding healthy controls (HC) from the correlation analysis (9, 22, 23, 26–28). A few studies investigated this correlation in all subjects including HC using an ROI approach (22, 23, 29, 30) and only a fraction of these extended the correlation analysis to whole brain white matter tracts using the voxel-wise approach (31–33). While the aforementioned studies present similar conditions, these studies either did not investigate the presence of gross or focal white matter abnormalities (31, 33) or did not exclude the subjects with such abnormalities from the extended investigation of neurocognitive deficits. The present study builds on these latter studies of the whole brain white matter tracts using a voxel-wise analysis approach to elucidate the relationship between neuro-cognition and white matter microstructure in sub-acute and chronic TBI patients with normal-appearing white matter (NAWM).
In this pilot study, we investigate if the classic DTI parameters (FA, MD, RD and AD) exhibit any changes in patients with sub-acute and chronic TBI of all severities but with no gross abnormalities on routine MRI. We also aim to explore if the relationship between the classic DTI parameters and neurocognitive function is altered when compared with HC. This knowledge is essential as part of considering whether classic DTI parameters are useful biomarkers of the degree of neurocognitive deficit after TBI.
Materials and Methods
This study was performed at St. Michael’s Hospital in Toronto, ON, Canada with approval from the hospital’s Research Ethics Board. Patients with a positive history of TBI were recruited from the Trauma/Neurosurgery Ward, Neurosurgery and Head Injury Clinics, and Trauma Registry Centre at the hospital. Inclusion criteria for TBI patients required that they were within 16–70 years of age and had sustained an isolated, closed head TBI with Glasgow Coma Score [GCS; (34)] in the range 5–15 as assessed within 12 h of injury. The HC were recruited through the word-of-mouth and through online advertisement. For HC, the inclusion criterion was to be within 16–70 years of age, while the exclusion criterion included having a history of TBI, motor vehicle collisions, or falls. The exclusion criteria, as applicable to both groups, included the presence of non-hemorrhagic or micro-hemorrhagic brain lesions as assessed by fluid-attenuated inversion recovery (FLAIR) and gradient recalled echo (GRE) MRI. The presence of small superficial contusions and superficial hemosiderin staining was not considered to be an exclusion criterion, as these lesions are not expected to affect the central white matter tracts that this study aimed to analyze. The exclusion criterion for both groups also included incomplete assessment of neuro-cognition as tested with standardized neuropsychological (NP) tests (see below for more details). Subjects with existing substance abuse were excluded from the study, as this dependency can potentially alter cognition (35), and is known to be associated with microstructural abnormalities with altered FA values (36). Both groups also were subject to standard MRI exclusion criteria (e.g., claustrophobia, ferromagnetic implants) and were required to possessed adequate verbal English language skills. All participants provided informed consent to participate in the study, and received monetary compensation for the time spent on the study.
Cognitive function was assessed in all subjects in the three categories of memory, executive function, and visual-motor coordination. These assessments were performed using the following standardized NP tests. Memory was assessed with the Digit Span subtest (37) of the Wechsler Adult Intelligence Scale (WAIS-III) (38) and the Hopkins Verbal Learning Test (HVLT) (39). Executive function was assessed using the Token Test (40), the Phonemic Fluency subtest, which is a part of the Neurosensory Center Comprehensive Examination for Aphasia (41), and the Digit Symbol subtest of the Wechsler Adult Intelligence Scale (WAIS-III) (38). The visuomotor coordination skill was assessed using the Grooved Pegboard test (42). The National Adult Reading Test (NART) (43) was used as a measure of verbal IQ. Except for the Grooved Pegboard, all tests were administered using pen and paper.
In the Digit Span test, the subjects were asked to recall numbers in forward and backward order and the number of correct responses was recorded in each case. The HVLT consisted of three trials; the total correct responses, incorrect responses, and intrusions were recorded for each trial and the number of correct delayed recall responses at the end of the three trials was also recorded. For the Token test, the total correct responses in each of the two trials performed were recorded. The Phonemic Fluency test required subjects to list words starting with each of the letters F, A, and S; total correct responses, incorrect/repeated responses, and intrusions were recorded in each case. For the Digit Symbol test, the total number of correctly completed, incorrectly completed, skipped, and incomplete boxes were recorded. For the Grooved Pegboard test, the time (in seconds) taken to put all the pegs in the board and the number of pegs dropped during the procedure were noted. Subjects completed the task first with their dominant hand, then the non-dominant hand. The test thus resulted in four scores: time to complete with the dominant hand (dominant time), the number of drops with the dominant hand (dominant drop), time to complete with the non-dominant hand (non-dominant time), and the number of drops with the non-dominant hand (non-dominant drop).
For each of the NP tests administered, a higher score indicates better performance. The one exception is the Grooved Pegboard test, in which a higher score indicates poorer visuomotor coordination.
Imaging was performed with a research-dedicated MRI system at 3.0 T (MR750, GE Healthcare, Waukesha, WI, USA) at Sunnybrook Health Sciences Centre in Toronto. The MRI protocol for each subject (TBI patients and controls) included a three-dimensional (3D) longitudinal relaxation time (T1)-weighted sequence to rule out gross structural abnormalities and for anatomical reference and two-dimensional (2D) FLAIR and GRE sequences to detect microbleeds, along with a 2D DTI sequence. The parameters for these MRI sequences were as follows: 3D axial T1-weighted inversion-recovery prepped fast gradient echo imaging—in-plane resolution = 0.9 mm × 0.9 mm; field of view (FOV) = 220 mm × 165 mm; slice thickness = 1.4 mm; flip angle = 15°; repetition time (TR) = 8.2 ms; echo time (TE) = 3.2 ms, axial FLAIR imaging—resolution = 0.9 mm × 1.1 mm × 3.0 mm; FOV = 220 mm × 200 mm; flip angle = 90°; TR = 9,950 ms; TE = 96 ms; axial GRE imaging—resolution = 0.8 mm × 1.0 mm × 3.0 mm; FOV = 200 mm × 200 mm; flip angle = 20°; TR = 784 ms; TE = 35 ms; and axial DTI imaging—FOV = 240 mm × 240 mm; resolution = 1.9 mm × 1.9 mm; slice thickness = 3 mm; non-collinear diffusion directions = 11 with b = 1,000 s/mm2 and with two volumes applied without diffusion sensitizing gradients (b = 0 s/mm2), averages = 4; TR = 9,400 ms; TE = 89.7 ms. The MRI protocol lasted approximately 40 min.
DTI Data Processing
Diffusion tensor imaging data were processed using the Oxford Centre for Functional MRI of the Brain (FMRIB) Software Library (FSL- version 5.0.5—http://www.fmrib.ox.ac.uk/fsl) (44). The diffusion data were corrected for distortions due to eddy currents and minor head motion using the b = 0 image as a reference. Brain extraction was then performed on these data using the FSL brain extraction tool (45). The diffusion tensor model was fitted to the images of the extracted brain tissue to produce maps of FA, MD, RD, and AD in individual subjects. Voxel-wise analyses of the DTI parameters were performed using the Tract-Based Spatial Statistics (TBSS) package (46). As part of TBSS processing, the FA image of each subject was transformed into the Montreal Neurological Institute (MNI)-152 standard space (47), using non-linear registration (44, 48) to the FMRIB58_FA_1mm image (average high-resolution FA image of 58 subjects in MNI 152 standard space). The mean FA maps of all subjects enrolled in the study were then created. A mean white matter skeleton was subsequently generated using a threshold of FA > 0.2 to include only the centers of white matter tracts common to the entire group. The FA data from each subject were then projected onto this mean FA skeleton to generate a 4D data matrix (x, y, z, FA) characterizing the group. The non-linear transformations used in this procedure were then applied to create analogous 4D data matrices of the skeletonized MD, RD, and AD values, respectively. These matrices were used for the inter-group statistical analyses.
A Student’s t-test was used to compare the mean NP test scores between the two groups, after assessing for normality of the data. Differences were considered statistically significant using a threshold p-value ≤0.05.
A General Linear Model (GLM) (49) was applied using a voxel-based approach to each of the 4D data matrices of FA, MD, AD, and RD within the FSL framework to investigate the relationship of DTI parameters with the group (TBI versus HC). The GLM was applied using non-parametric permutation-based statistics by employing the randomize tool in FSL with threshold-free cluster enhancement and 5,000 permutations (50). Age, sex, and time since injury were included as non-interacting covariates in the GLM. The contrasts were designed to investigate if the two groups had statistically significant differences in DTI parameters (TBI > HC and HC > TBI).
The investigation of whether DTI parameters were correlated with NP test performance was also performed with the randomize tool within the framework of GLM in FSL. For this analysis, a composite score for each category of NP function (memory, executive function, and visuomotor coordination) was calculated using the following procedure. First, the Pearson correlation coefficient (R) was calculated to assess the correlation between scores and sub-scores of various tests within the same category and only those subtests that had significant positive correlation with each other (R > 0.5, p-value ≤0.05) were included in the subsequent analyses. These selected tests/subtests scores for all subjects (TBI and HC) were then converted into ranks. Finally, a composite score for each category was calculated as the mean of ranks of all subtests and tests in that category. The correlation and rank calculations were performed using the Statistical Analysis System, 9.4 (SAS Institute Inc. 2013, Cary, NC, USA). These composites scores for NP tests were entered in the GLM design matrix for investigation of the correlation between DTI parameters and NP test performance. First, an investigation of whether the correlation between NP test performance and DTI parameters varied according to the group was performed. For this, an interaction term for each NP test category was added to the design matrix and contrasts were designed to determine positive and negative interactions between composite test scores in the three NP categories and group predictor variable. Afterward, correlations of the composite scores in the three NP categories were investigated with contrasts designed to calculate positive as well as negative correlations. For these analyses involving NP tests, IQ and education were also included in the GLM as non-interacting covariates.
A p-value ≤0.05, corrected for multiple comparisons within randomize by controlling for the family-wise error rate (51), was considered significant for all statistical comparisons performed within the framework of the GLM in FSL.
Identifying the Regions with Significant Inter-Group Differences/Correlations
The randomize outputs consisting of p-values corrected for multiple comparisons were thresholded to keep voxels with p-value ≤0.05. The FSL atlas query and cluster identification tools (merged as the autoaq algorithm) were used on the thresholded images to identify various white matter tracts using the Johns Hopkins University White Matter Labels atlas (52). To eliminate noise due to averaging and normalization, tracts with <50 significantly different voxels were not considered as tracts exhibiting structural differences between the two groups and were excluded from the reported results.
Post Hoc Analyses
The post hoc analyses, where mentioned, were performed using only those voxels that were found significant at p ≤ 0.05 level. The post hoc analyses were also performed within the framework of the GLM in FSL and the results were corrected for multiple comparisons within randomize by controlling for the family-wise error rate (51).
A total of 52 patients and 18 HC met the inclusion criteria and underwent MRI and standardized NP tests. FLAIR and GRE images of all the subjects were reviewed by two board-certified neuroradiologists to rule out non-hemorrhagic (53) and micro-hemorrhagic (54) lesions. Twelve patients were excluded from the study: five did not complete the MRI exam, four had significant head motion during MRI and three withdrew their consent after MRI. Twenty-one patients and one HC were excluded from the analyses due to gross abnormalities in the central cerebral white matter, visible on FLAIR or GRE images. Most of the TBI subjects had small superficial contusions and superficial hemosiderin staining but were not excluded, as these lesions are not expected to affect the central white matter tracts that this study aimed to analyze. Nineteen patients with a history of trauma from a motor vehicle collision (5 patients) or fall (14 patients) were thus included in the study. Based on the GCS scores at the time of emergency department arrival, 15 patients had mild TBI (reference GCS: 13–15), three patients had moderate TBI (reference GCS: 9–12) and one patient had severe TBI (reference GCS ≤8). Time elapsed between injury and MRI for the present study ranged from 0.2 to 5.1 years (Table 1). Nine subjects were in the sub-acute stage (2 weeks < time since injury < 1 year), whereas the rest of the subjects were in the chronic stage after TBI (time since injury > 1 year) (7).
The group of HC consisted of individuals recruited through online advertisement as well as friends and family members of patients who attended the injury, trauma, or neurosurgery clinics at the hospital. The group of HC was matched by age, sex, and education to the group of TBI patients. Table 1 lists the demographic details of all subjects in the study.
Healthy controls performed better than the TBI group on all NP tests, based on score results. However, only the scores for Digit Span, Phonemic Fluency, the third trial of the HVLT, and Digit Symbol tests were significantly higher in HC at p-value ≤0.05. A significant positive correlation was found between the Digit Span scores (forward and backward), the total correct responses in each of the three trials of HVLT, and the correct delayed recall score in the memory category. In the executive function category, a significant positive correlation was found between the total correct phonemic score, the total correct score in the Digit Symbol test, and the total score in each of the two Token test trials. In the visuomotor coordination category, the dominant and non-dominant times were found to have a significant positive correlation. No significant difference was found in the verbal IQ of the two groups as measured with the NART. The results of NP test performance are summarized in Table 2.
Inter-Group Differences in DTI Parameters
Significantly increased MD and AD values were found in the TBI group compared to HC in several brain regions (controlling for age, sex, and time since injury). Table 3 lists MD and AD values in those tracts where significant inter-group differences in these variables were found. In particular, significant differences in MD and AD were found in the corpus callosum and bilateral internal capsule, external capsule, corona radiata, posterior thalamic radiations, and various other structures. Changes in AD were more widespread compared to those in MD (Figures 1 and 2). No significant differences between the two groups were found in FA and RD. Although trends of increased RD and FA in TBI were visible in several major white matter tracts, notable exceptions included a trend of decreased FA in TBI in the corpus callosum and the bilateral cingulum.
Figure 1. Mean diffusivity (MD) values are significantly higher in the traumatic brain injury (TBI) group with normal-appearing white matter. The MD contrast is overlaid on a standard Montreal Neurological Institute 152 T1 1 mm brain and the mean fractional anisotropy skeleton (green—display threshold 0.2–0.8). The voxels where MD was found to be significantly higher in the TBI group (p-value ≤ 0.05) are shown in red. Every fourth transverse slice only is shown here to approximately cover the entire brain along the superior–inferior axis. Anatomical right side of the brain is shown on the left in this figure.
Figure 2. Axial diffusivity (AD) values are significantly higher in the traumatic brain injury (TBI) group with normal-appearing white matter. The AD contrast is overlaid on a standard Montreal Neurological Institute 152 T1 1 mm brain and the mean fractional anisotropy skeleton (green—display threshold 0.2–0.8). The voxels where mean diffusivity was found to be significantly higher in the TBI group (p-value ≤ 0.05) are shown in red. Every fourth transverse slice only is shown here to approximately cover the entire brain along the superior–inferior axis. Anatomical right side of the brain is shown on the left in this figure.
Post Hoc Analyses of the Inter-Group Differences in DTI Parameters
The post hoc analyses, involving only the significant voxels at p ≤ 0.05 level confirmed the inter-group differences in the DTI variables. The p values of the post hoc analyses are also listed in Table 3.
No significant interactions between the composite test scores in the three NP function categories and group variables were found, indicating that both groups exhibited similar relationships between the NP explanatory variables on the DTI parameters. Therefore, the NP-group interaction terms were excluded, and contrasts were appropriately modified for the subsequent analyses to investigate the effect of group and NP explanatory variables on the DTI parameters.
Correlation of NP Test Scores with DTI Parameters
Continuing with the GLM design matrix with the composite test scores in the three NP function categories added as covariates along with age, sex, time since injury, education and the NART (verbal IQ), the correlation between NP test scores and DTI parameters was investigated with appropriate contrasts. In this analysis, a negative correlation was found between the composite memory score and FA in all subjects, in several tracts widespread in the entire brain. These tracts included bilateral internal and external capsules, corona radiata, cingulum, posterior thalamic radiations and superior longitudinal fasciculi. A positive correlation between composite visuomotor score and RD was also found in both groups, in a few tracts including the corpus callosum and only the right components of corona radiata, cingulum, and superior longitudinal fasciculus. Table 4 lists the tracts where a significant correlation between the DTI parameters and the composite NP test scores was found. No significant correlations between the composite NP function scores and MD or AD were found.
Table 4. Tracts with significant correlations between diffusion tensor imaging parameters and composite neuropsychological scores.
This study was designed to investigate the microstructural abnormalities in NAWM using DTI and to explore their potential impact on cognitive function in post-acute TBI patients compared to HC using a voxel-wise approach. Most often, TBI patients report persisting cognitive deficits and hence are not able to return to their pre-morbid daily activities, sometimes also including the inability to resume work, resulting in financial stress and increased psychological stress. The clinical assessment of this persisting cognitive deficit is challenging, especially when no gross abnormalities indicative of the cognitive or functional loss are detected on routine diagnostic imaging (e.g., structural MRI and CT). The DTI method, however, has been shown to detect micro-abnormalities and hence can be helpful in explaining the functional or cognitive deficit in sub-acute and chronic TBI with NAWM. Voxel-wise analysis of the DTI parameters and of the correlation of cognitive function with DTI parameters performed over the entire brain has the advantage of identifying partial regions or tracts with significant inter-group differences and significant correlations. Such regions may not be successfully detected via analyses which involve ROI selection based on a priori hypotheses. Partial regions or tracts with significant inter-group differences may also be masked by spatial averaging in ROI analyses.
Our neuroradiologists confirmed the absence of gross abnormalities in the central white matter of the brain in only 19 of the 40 TBI patients after reviewing their FLAIR and GRE data. These numbers suggest that some gross abnormalities due to TBI are likely to persist for lengthy time periods and may account for the changes in DTI parameters in chronic TBI, as widely quoted in the literature. It is noteworthy, however, that most of the TBI patients with NAWM that were included in our analyses had mild TBI (15 out of 19 subjects). Mild TBI is more frequent than other forms of TBI world-wide (55), and is less likely to cause gross lesions in the brain.
Our finding of widespread increase in the overall diffusivity of water molecules (as measured with MD) in the post-acute TBI subjects, while controlling for age, sex, and time since injury, robustly reaffirms the findings of previous studies (12, 13, 21, 24, 32, 33, 56, 57) and indicates the persistence of DAI in sub-acute and chronic TBI with NAWM. We also found an accompanying increase in AD in the TBI group. Many DTI studies of TBI have not analyzed AD. However, those studies that explored inter-group differences in AD also found an increase in this parameter in post-acute TBI (21, 32). The changes in AD were more widespread compared to changes in MD, and potentially significantly contributed to increase in MD.
We did not find any significant differences between the two groups in FA. Several DTI studies of TBI reported a decrease of FA in post-acute TBI; however, these studies included more acute patients, up to 1 month from the injury (22, 23, 27, 30). Some other studies with a minimum post-injury time similar to ours (2–3 months) reported variable FA results: some reported a decrease (12, 24, 26, 31, 32); one reported an increase (58); and others reported no change (9, 25, 59). One potential reason for this variability may be the presence of gross abnormalities in the brain, as some of the studies that reported decreased FA either included subjects with gross abnormalities (24, 32) or did not screen subjects for gross abnormalities (26). Furthermore, one study that found reduced FA (12) focused on TBI subjects using prescribed medicines for major depression or post-traumatic stress disorder, which constitute a known confounding factor (60).
We also did not find any significant differences between the two groups in RD; however, trends of increased RD in TBI were evident. Very often, studies of DTI in post- acute TBI have not investigated the RD parameter. A few studies involving RD reported variable findings, with increase (32) or no difference (26) observed in chronic TBI compared to HC.
A survey of the TBI literature suggests that DTI parameters change dynamically with time post-injury. One longitudinal study assessing FA values in TBI patients reported a decrease at the acute and sub-acute stage (range of time since injury = 3–55 days) when compared with controls, an increase at the chronic stage (minimum of 238 days since the first scan) in comparison with patients’ first measurements and no difference at the chronic stage when compared with controls (25). The same study also reported an increase in MD in the TBI group in comparison with controls at both time-points, and an increase in MD with time in the TBI group. Another longitudinal study reported increases in AD, MD, and FA in the acute stage, and no differences in RD; and sustained increase in MD in the chronic stage accompanied by an increase in RD, normal AD levels, and a reduction in FA compared to HC (33). These observations suggest that the absence of statistically significant differences in RD and FA in our study might be due to inclusion of both sub-acute and chronic TBI patients.
Considering these findings collectively with our results, we conclude that DAI, as manifested in elevated diffusivity, potentially persist from the acute stage to well into the chronic stage after TBI, whereas changes in FA may be more sensitive to time since injury and may not remain detectable later on in the chronic stage. The different effects shown by these parameters are interesting, and may reflect their biophysical underpinnings. Whereas MD is a linear combination of diffusivities in the three orthogonal directions and will increase if all three diffusivities are increased, FA is a measure of difference among the three diffusivities and may show no change if the three diffusivities are all modulated proportionally in the same manner. An increase in AD and MD and no difference in FA in our study, therefore, suggest a trend of increased RD values in TBI subjects. Furthermore, AD is thought to be a measure of axonal integrity (61). By extension, our results suggest that the repair of the axonal cytoskeleton, thought to result in rising AD values during recovery (24), happens before the repair of axonal membranes or re-myelination of axons, both of which are thought to restrict perpendicular diffusion and to help normalize increased RD values (62).
Persisting elevated MD values in sub-acute and chronic TBI, encompassing association (e.g., external capsule), projection (e.g., corona radiata, thalamic radiation), and commissural (e.g., corpus callosum) fibers, despite overall NAWM, suggest that microstructural damage or pathological change remains a salient feature of chronic TBI. The significant changes in MD which represent a cumulative effect of AD and RD values are of primary importance, possibly suggesting that MD can be used for assessing the burden of microstructural damage in sub-acute and chronic TBI characterized with NAWM. Our post hoc analyses also confirmed these differences.
Our investigation of the relationship between NP function and DTI parameters primarily revealed effects involving FA and RD, and not AD and MD. The absence of interaction effects between NP test scores and DTI parameters implies that positive and negative correlations hold similarly in all subjects regardless of the group status. A significant negative correlation was found between the composite memory score and FA in several white matter tracts. Notably, the negative correlation was found in the structures of the corpus callosum, bilateral external capsule, bilateral posterior thalamic radiations, and bilateral superior longitudinal fasciculi, which constitute either the commissural, projection, or long-range association fibers. Working and short-term memory is known to be associated with various regions in the frontal, parietal, and anterior-cingulate cortices and basal ganglia (63). The white matter tracts involved in memory are the relatively short-length or U-association fibers connecting these cortical regions. The major white matter tracts that are known to be involved in memory include the fornix (64, 65), which projects to the cingulate cortex, and the internal capsule (66), which connects the cortex with basal ganglia. Furthermore, the DTI literature suggests a positive correlation of FA with function, indicating higher FA (and hence intact structure of axons, with no disruption in myelin or axonal membranes) in subjects with better performance on memory (67) as well as other cognitive tasks (68, 69). Interestingly, our results show a negative correlation of FA in the white tracts which are not directly involved in memory. Keeping in mind the evidence of widespread DAI in the brain, we explain this finding by suggesting that whereas FA increases within the network of interest are correlated with improved performance, the negative correlation of FA in another network or some sub-network might indicate hyper-function resulting in less effective mental processing and thus a decrement in behavioral performance.
We also found a significant positive correlation between RD and visuomotor performance in the corpus callosum and right aspects of corona radiata, cingulum, and superior longitudinal fasciculus. This means that poor visuomotor coordination (higher scores on the Pegboard task) were associated with elevated RD values. Each of the structures where a positive correlation of RD with visuomotor task performance is found is known to be involved in the visual or motor tasks: the corona radiata are associated with the corticospinal tract, the superior longitudinal fasciculi are involved in integrating the motor and decision-making centers with ipsilateral visual centers (70), and the posterior aspect of cingulum is involved in visual skills. Although other studies have reported tracing of superior longitudinal fasciculi associated with investigation of language networks (71), our study is the first to report the involvement of this structure in visuomotor skills using DTI parameters.
Our study has some limitations, which include the collection of DTI data in only 12 diffusion-weighted directions with coarse resolution, the small sample size, and the inclusion of subjects with all TBI severities. Another limitation is that we did not use susceptibility-weighted imaging (SWI) to detect microstructural abnormalities in the NAWM because the protocol was not available on the MRI scanner at the time of data collection. As SWI is known to be superior to GRE in detecting micro-hemorrhages (72), it can be argued that the differences between the two groups might be influenced by undetected micro-hemorrhages in NAWM, which were not visible on GRE. However, the inter-group differences that were found are well supported by earlier studies (12, 24, 32, 73). The similarity of the effect of NP performance on the DTI parameters across both groups and the lack of interaction between NP performance and group variable also indirectly supports that no gross abnormalities were present in the white matter of subjects with TBI.
Conclusion and Future Directions
Overall, this pilot study shows DTI to be helpful in evaluating the microstructural abnormalities in sub-acute and chronic TBI patients even when there are no apparent abnormalities in the central white matter. Elevated diffusivity values appear to persist well into the chronic stage after TBI, and hence MD can be used to assess persisting microstructural damage in chronic TBI. No clear trend for FA values in post-acute TBI is visible. The study also indicates the potential use of DTI in assessing neurocognitive functions in both TBI patients and HC. The authors recommend larger and longitudinal studies for each of mild, moderate, and severe TBI, to establish the use of DTI parameters as biomarkers of DAI and neurocognitive function in patients with NAWM.
Research Ethics Board, St. Michael’s Hospital, Toronto, ON, Canada. M5B 1W8. Telephone: 416-864-6060 extension: 2557. REB # 11-282 (Strategic Teams in Applied Injury Research Protocol). Informed consent was obtained from the participants on a paper form, after a research team member explained the objectives and methods of the study to them.
EH and EC performed the neuroimaging analyses with help from CL and OC. NP provided input for TBSS analyses. RJ helped with the statistical analyses. MC and AB were involved in study design. SZ assisted with patient recruitment and research coordination. EH, EC, NP, and MC wrote the manuscript. All authors participated in discussing the results and editing the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer RF and handling Editor declared their shared affiliation, and the handling Editor states that the process nevertheless met the standards of a fair and objective review.
This research was supported by the Canadian Institutes of Health Research Strategic Team Grant in Applied Injury Research #TIR-103946, and the Ontario Neurotrauma Foundation (www.onf.org). The authors are grateful to Dr. Amy W. Lin (Board-certified radiologist) for screening the MRI scans of all subjects for gross abnormalities in the central white matter of the brain.
2. Arenth PM, Russell KC, Scanlon JM, Kessler LJ, Ricker JH. Corpus callosum integrity and neuropsychological performance after traumatic brain injury: a diffusion tensor imaging study. J Head Trauma Rehabil (2014) 29(2):E1–10. doi:10.1097/HTR.0b013e318289ede5
3. Rugg-Gunn FJ, Symms MR, Barker GJ, Greenwood R, Duncan JS. Diffusion imaging shows abnormalities after blunt head trauma when conventional magnetic resonance imaging is normal. J Neurol Neurosurg Psychiatry (2001) 70(4):530–3. doi:10.1136/jnnp.70.4.530
4. Wang JY, Bakhadirov K, Devous MD Sr, Abdi H, McColl R, Moore C, et al. Diffusion tensor tractography of traumatic diffuse axonal injury. Arch Neurol (2008) 65(5):619–26. doi:10.1001/archneur.65.5.619
5. Marquez de la Plata CD, Yang FG, Wang JY, Krishnan K, Bakhadirov K, Paliotta C, et al. Diffusion tensor imaging biomarkers for traumatic axonal injury: analysis of three analytic methods. J Int Neuropsychol Soc (2011) 17(1):24–35. doi:10.1017/S1355617710001189
6. Perez AM, Adler J, Kulkarni N, Strain JF, Womack KB, Diaz-Arrastia R, et al. Longitudinal white matter changes after traumatic axonal injury. J Neurotrauma (2014) 31(17):1478–85. doi:10.1089/neu.2013.3216
7. Hulkower MB, Poliak DB, Rosenbaum SB, Zimmerman ME, Lipton ML. A decade of DTI in traumatic brain injury: 10 years and 100 articles later. AJNR Am J Neuroradiol (2013) 34(11):2064–74. doi:10.3174/ajnr.A3395
8. Kumar R, Gupta RK, Husain M, Chaudhry C, Srivastava A, Saksena S, et al. Comparative evaluation of corpus callosum DTI metrics in acute mild and moderate traumatic brain injury: its correlation with neuropsychometric tests. Brain Inj (2009) 23(7):675–85. doi:10.1080/02699050903014915
9. Levin HS, Wilde E, Troyanskaya M, Petersen NJ, Scheibel R, Newsome M, et al. Diffusion tensor imaging of mild to moderate blast-related traumatic brain injury and its sequelae. J Neurotrauma (2010) 27(4):683–94. doi:10.1089/neu.2009.1073
10. Rutgers DR, Toulgoat F, Cazejust J, Fillard P, Lasjaunias P, Ducreux D. White matter abnormalities in mild traumatic brain injury: a diffusion tensor imaging study. AJNR Am J Neuroradiol (2008) 29(3):514–9. doi:10.3174/ajnr.A0856
11. Cubon VA, Putukian M, Boyer C, Dettwiler A. A diffusion tensor imaging study on the white matter skeleton in individuals with sports-related concussion. J Neurotrauma (2011) 28(2):189–201. doi:10.1089/neu.2010.1430
12. Yurgelun-Todd DA, Bueler CE, McGlade EC, Churchwell JC, Brenner LA, Lopez-Larson MP. Neuroimaging correlates of traumatic brain injury and suicidal behavior. J Head Trauma Rehabil (2011) 26(4):276–89. doi:10.1097/HTR.0b013e31822251dc
13. Xiong K, Zhu Y, Zhang Y, Yin Z, Zhang J, Qiu M, et al. White matter integrity and cognition in mild traumatic brain injury following motor vehicle accident. Brain Res (2014) 1591:86–92. doi:10.1016/j.brainres.2014.10.030
14. Nakayama N, Okumura A, Shinoda J, Yasokawa YT, Miwa K, Yoshimura SI, et al. Evidence for white matter disruption in traumatic brain injury without macroscopic lesions. J Neurol Neurosurg Psychiatry (2006) 77(7):850–5. doi:10.1136/jnnp.2005.077875
15. Oni MB, Wilde EA, Bigler ED, McCauley SR, Wu TC, Yallampalli R, et al. Diffusion tensor imaging analysis of frontal lobes in pediatric traumatic brain injury. J Child Neurol (2010) 25(8):976–84. doi:10.1177/0883073809356034
16. Goetz P, Blamire A, Rajagopalan B, Cadoux-Hudson T, Young D, Styles P. Increase in apparent diffusion coefficient in normal appearing white matter following human traumatic brain injury correlates with injury severity. J Neurotrauma (2004) 21(6):645–54. doi:10.1089/0897715041269731
17. Rutgers DR, Fillard P, Paradot G, Tadie M, Lasjaunias P, Ducreux D. Diffusion tensor imaging characteristics of the corpus callosum in mild, moderate, and severe traumatic brain injury. AJNR Am J Neuroradiol (2008) 29(9):1730–5. doi:10.3174/ajnr.A1213
18. Bigler ED. Distinguished Neuropsychologist Award Lecture 1999. The lesion(s) in traumatic brain injury: implications for clinical neuropsychology. Arch Clin Neuropsychol (2001) 16(2):95–131. doi:10.1016/S0887-6177(00)00095-0
19. Chang MC, Kim SH, Kim OL, Bai DS, Jang SH. The relation between fornix injury and memory impairment in patients with diffuse axonal injury: a diffusion tensor imaging study. NeuroRehabilitation (2010) 26(4):347–53. doi:10.3233/NRE-2010-0572
20. Newcombe VF, Williams GB, Scoffings D, Cross J, Carpenter TA, Pickard JD, et al. Aetiological differences in neuroanatomy of the vegetative state: insights from diffusion tensor imaging and functional implications. J Neurol Neurosurg Psychiatry (2010) 81(5):552–61. doi:10.1136/jnnp.2009.196246
21. Newcombe V, Chatfield D, Outtrim J, Vowler S, Manktelow A, Cross J, et al. Mapping traumatic axonal injury using diffusion tensor imaging: correlations with functional outcome. PLoS One (2011) 6(5):e19214. doi:10.1371/journal.pone.0019214
22. Niogi SN, Mukherjee P, Ghajar J, Johnson C, Kolster RA, Sarkar R, et al. Extent of microstructural white matter injury in postconcussive syndrome correlates with impaired cognitive reaction time: a 3T diffusion tensor imaging study of mild traumatic brain injury. AJNR Am J Neuroradiol (2008) 29(5):967–73. doi:10.3174/ajnr.A0970
23. Niogi SN, Mukherjee P, Ghajar J, Johnson CE, Kolster R, Lee H, et al. Structural dissociation of attentional control and memory in adults with and without mild traumatic brain injury. Brain (2008) 131(Pt 12):3209–21. doi:10.1093/brain/awn247
24. Bonnelle V, Leech R, Kinnunen KM, Ham TE, Beckmann CF, De Boissezon X, et al. Default mode network connectivity predicts sustained attention deficits after traumatic brain injury. J Neurosci (2011) 31(38):13442–51. doi:10.1523/JNEUROSCI.1163-11.2011
25. Grossman EJ, Jensen JH, Babb JS, Chen Q, Tabesh A, Fieremans E, et al. Cognitive impairment in mild traumatic brain injury: a longitudinal diffusional kurtosis and perfusion imaging study. AJNR Am J Neuroradiol (2013) 34(5):S951–3. doi:10.3174/ajnr.A3358
26. Kraus MF, Susmaras T, Caughlin BP, Walker CJ, Sweeney JA, Little DM. White matter integrity and cognition in chronic traumatic brain injury: a diffusion tensor imaging study. Brain (2007) 130(Pt 10):2508–19. doi:10.1093/brain/awm216
27. Hartikainen KM, Waljas M, Isoviita T, Dastidar P, Liimatainen S, Solbakk AK, et al. Persistent symptoms in mild to moderate traumatic brain injury associated with executive dysfunction. J Clin Exp Neuropsychol (2010) 32(7):767–74. doi:10.1080/13803390903521000
28. Maruta J, Suh M, Niogi SN, Mukherjee P, Ghajar J. Visual tracking synchronization as a metric for concussion screening. J Head Trauma Rehabil (2010) 25(4):293–305. doi:10.1097/HTR.0b013e3181e67936
29. Wozniak JR, Krach L, Ward E, Mueller BA, Muetzel R, Schnoebelen S, et al. Neurocognitive and neuroimaging correlates of pediatric traumatic brain injury: a diffusion tensor imaging (DTI) study. Arch Clin Neuropsychol (2007) 22(5):555–68. doi:10.1016/j.acn.2007.03.004
30. Grossman EJ, Ge Y, Jensen JH, Babb JS, Miles L, Reaume J, et al. Thalamus and cognitive impairment in mild traumatic brain injury: a diffusional kurtosis imaging study. J Neurotrauma (2012) 29(13):2318–27. doi:10.1089/neu.2011.1763
32. Kinnunen KM, Greenwood R, Powell JH, Leech R, Hawkins PC, Bonnelle V, et al. White matter damage and cognitive impairment after traumatic brain injury. Brain (2011) 134(Pt 2):449–63. doi:10.1093/brain/awq347
33. Croall ID, Cowie CJ, He J, Peel A, Wood J, Aribisala BS, et al. White matter correlates of cognitive dysfunction after mild traumatic brain injury. Neurology (2014) 83(6):494–501. doi:10.1212/WNL.0000000000000666
36. Lim KO, Wozniak JR, Mueller BA, Franc DT, Specker SM, Rodriguez CP, et al. Brain macrostructural and microstructural abnormalities in cocaine dependence. Drug Alcohol Depend (2008) 92(1–3):164–72. doi:10.1016/j.drugalcdep.2007.07.019
37. Heinly MT, Greve KW, Bianchini KJ, Love JM, Brennan A. WAIS digit span-based indicators of malingered neurocognitive dysfunction: classification accuracy in traumatic brain injury. Assessment (2005) 12(4):429–44. doi:10.1177/1073191105281099
38. Burton DB, Ryan JJ, Axelrod BN, Schellenberger T. A confirmatory factor analysis of the WAIS-III in a clinical sample with crossvalidation in the standardization sample. Arch Clin Neuropsychol (2002) 17(4):371–87. doi:10.1016/S0887-6177(01)00121-4
39. Shapiro AM, Benedict RH, Schretlen D, Brandt J. Construct and concurrent validity of the Hopkins Verbal Learning Test-revised. Clin Neuropsychol (1999) 13(3):348–58. doi:10.1076/clin.13.3.348.1749
44. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, et al. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage (2004) 23(Suppl 1):S208–19. doi:10.1016/j.neuroimage.2004.07.051
46. Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, Mackay CE, et al. Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. Neuroimage (2006) 31(4):1487–505. doi:10.1016/j.neuroimage.2006.02.024
47. Grabner G, Janke AL, Budge MM, Smith D, Pruessner J, Collins DL. Symmetric atlasing and model based segmentation: an application to the hippocampus in older adults. Med Image Comput Comput Assist Interv (2006) 9(Pt 2):58–66. doi:10.1007/11866763_8
48. Andersson JLR, Jenkinson M, Smith S. Non-Linear Registration, aka Spatial Normalization. Oxford, UK: FMRIB Analysis Group, University of Oxford (2007). 2 p. FMRIB Technical Report TR07JA2. Available from: http://www.fmrib.ox.ac.uk/datasets/techrep/
53. Parizel PM, Van Goethem JW, Ozsarlak O, Maes M, Phillips CD. New developments in the neuroradiological diagnosis of craniocerebral trauma. Eur Radiol (2005) 15(3):569–81. doi:10.1007/s00330-004-2558-z
54. Huisman TA, Sorensen AG, Hergan K, Gonzalez RG, Schaefer PW. Diffusion-weighted imaging for the evaluation of diffuse axonal injury in closed head injury. J Comput Assist Tomogr (2003) 27(1):5–11. doi:10.1097/00004728-200301000-00002
55. Cassidy JD, Carroll LJ, Peloso PM, Borg J, von Holst H, Holm L, et al. Incidence, risk factors and prevention of mild traumatic brain injury: results of the WHO collaborating centre task force on mild traumatic brain injury. J Rehabil Med (2004) 43:28–60. doi:10.1080/16501960410023732
56. Inglese M, Makani S, Johnson G, Cohen BA, Silver JA, Gonen O, et al. Diffuse axonal injury in mild traumatic brain injury: a diffusion tensor imaging study. J Neurosurg (2005) 103(2):298–303. doi:10.3171/jns.2005.103.2.0298
57. Sidaros A, Engberg AW, Sidaros K, Liptrot MG, Herning M, Petersen P, et al. Diffusion tensor imaging during recovery from severe traumatic brain injury and relation to clinical outcome: a longitudinal study. Brain (2008) 131(Pt 2):559–72. doi:10.1093/brain/awm294
58. Lo C, Shifteh K, Gold T, Bello JA, Lipton ML. Diffusion tensor imaging abnormalities in patients with mild traumatic brain injury and neurocognitive impairment. J Comput Assist Tomogr (2009) 33(2):293–7. doi:10.1097/RCT.0b013e31817579d1
59. Lange RT, Iverson GL, Brubacher JR, Madler B, Heran MK. Diffusion tensor imaging findings are not strongly associated with postconcussional disorder 2 months following mild traumatic brain injury. J Head Trauma Rehabil (2012) 27(3):188–98. doi:10.1097/HTR.0b013e318217f0ad
60. Murphy ML, Frodl T. Meta-analysis of diffusion tensor imaging studies shows altered fractional anisotropy occurring in distinct brain areas in association with depression. Biol Mood Anxiety Disord (2011) 1(1):3. doi:10.1186/2045-5380-1-3
61. Di Paola M, Spalletta G, Caltagirone C. In vivo structural neuroanatomy of corpus callosum in Alzheimer’s disease and mild cognitive impairment using different MRI techniques: a review. J Alzheimers Dis (2010) 20(1):67–95. doi:10.3233/JAD-2010-1370
65. Tsivilis D, Vann SD, Denby C, Roberts N, Mayes AR, Montaldi D, et al. A disproportionate role for the fornix and mammillary bodies in recall versus recognition memory. Nat Neurosci (2008) 11(7):834–42. doi:10.1038/nn.2149
66. Tatemichi TK, Desmond DW, Prohovnik I, Cross DT, Gropen TI, Mohr JP, et al. Confusion and memory loss from capsular genu infarction: a thalamocortical disconnection syndrome? Neurology (1992) 42(10):1966–79. doi:10.1212/WNL.42.10.1966
67. Palacios EM, Fernandez-Espejo D, Junque C, Sanchez-Carrion R, Roig T, Tormos JM, et al. Diffusion tensor imaging differences relate to memory deficits in diffuse traumatic brain injury. BMC Neurol (2011) 11:24. doi:10.1186/1471-2377-11-24
68. Grieve SM, Williams LM, Paul RH, Clark CR, Gordon E. Cognitive aging, executive function, and fractional anisotropy: a diffusion tensor MR imaging study. AJNR Am J Neuroradiol (2007) 28(2):226–35.
69. Sugiyama K, Kondo T, Oouchida Y, Suzukamo Y, Higano S, Endo M, et al. Clinical utility of diffusion tensor imaging for evaluating patients with diffuse axonal injury and cognitive disorders in the chronic stage. J Neurotrauma (2009) 26(11):1879–90. doi:10.1089/neu.2008-0839
70. Schmahmann JD, Smith EE, Eichler FS, Filley CM. Cerebral white matter: neuroanatomy, clinical neurology, and neurobehavioral correlates. Ann N Y Acad Sci (2008) 1142:266–309. doi:10.1196/annals.1444.017
71. Kamali A, Flanders AE, Brody J, Hunter JV, Hasan KM. Tracing superior longitudinal fasciculus connectivity in the human brain using high resolution diffusion tensor tractography. Brain Struct Funct (2014) 219(1):269–81. doi:10.1007/s00429-012-0498-y
73. Ljungqvist J, Nilsson D, Ljungberg M, Sorbo A, Esbjornsson E, Eriksson-Ritzen C, et al. Longitudinal study of the diffusion tensor imaging properties of the corpus callosum in acute and chronic diffuse axonal injury. Brain Inj (2011) 25(4):370–8. doi:10.3109/02699052.2011.558038
Keywords: diffusion tensor imaging, chronic traumatic brain injury, white matter microstructure, human, voxel-based analysis, sub-acute traumatic brain injury, normal-appearing white matter
Citation: Hashim E, Caverzasi E, Papinutto N, Lewis CE, Jing R, Charles O, Zhang S, Lin A, Graham SJ, Schweizer TA, Bharatha A and Cusimano MD (2017) Investigating Microstructural Abnormalities and Neurocognition in Sub-Acute and Chronic Traumatic Brain Injury Patients with Normal-Appearing White Matter: A Preliminary Diffusion Tensor Imaging Study. Front. Neurol. 8:97. doi: 10.3389/fneur.2017.00097
Received: 05 August 2016; Accepted: 28 February 2017;
Published: 20 March 2017
Edited by:Renato Anghinah, University of São Paulo, Brazil
Reviewed by:Eric Peter Thelin, University of Cambridge, UK
Renata Areza-Fegyveres, University of São Paulo, Brazil
Copyright: © 2017 Hashim, Caverzasi, Papinutto, Lewis, Jing, Charles, Zhang, Lin, Graham, Schweizer, Bharatha and Cusimano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael D. Cusimano, firstname.lastname@example.org