OUP user menu

Biomarkers of increased diffusion anisotropy in semi-acute mild traumatic brain injury: a longitudinal perspective

Josef M. Ling, Amanda Peña, Ronald A. Yeo, Flannery L. Merideth, Stefan Klimaj, Charles Gasparovic, Andrew R. Mayer
DOI: http://dx.doi.org/10.1093/brain/aws073 1281-1292 First published online: 14 April 2012

Summary

Mild traumatic brain injury is the most prevalent neurological insult and frequently results in neurobehavioural sequelae. However, little is known about the pathophysiology underlying the injury and how these injuries change as a function of time. Although diffusion tensor imaging holds promise for in vivo characterization of white matter pathology, both the direction and magnitude of anisotropic water diffusion abnormalities in axonal tracts are actively debated. The current study therefore represents both an independent replication effort (n = 28) of our previous findings (n = 22) of increased fractional anisotropy during semi-acute injury, as well as a prospective study (n = 26) on the putative recovery of diffusion abnormalities. Moreover, new analytical strategies were applied to capture spatially heterogeneous white matter injuries, which minimize implicit assumptions of uniform injury across diverse clinical presentations. Results indicate that whereas a general pattern of high anisotropic diffusion/low radial diffusivity was present in various white matter tracts in both the replication and original cohorts, this pattern was only consistently observed in the genu of the corpus callosum across both samples. Evidence for a greater number of localized clusters with increased anisotropic diffusion was identified across both cohorts at trend levels, confirming heterogeneity in white matter injury. Pooled analyses (50 patients; 50 controls) suggested that measures of diffusion within the genu were predictive of patient classification, albeit at very modest levels (71% accuracy). Finally, we observed evidence of recovery in lesion load in returning patients across a 4-month interval, which was correlated with a reduction in self-reported post-concussive symptomatology. In summary, the corpus callosum may serve as a common point of injury in mild traumatic brain injury secondary to anatomical (high frequency of long unmyelinated fibres) and biomechanics factors. A spatially heterogeneous pattern of increased anisotropic diffusion exists in various other white matter tracts, and these white matter anomalies appear to diminish with recovery. This macroscopic pattern of diffusion abnormalities may be associated with cytotoxic oedema following mechanical forces, resulting in changes in ionic homeostasis, and alterations in the ratio of intracellular and extracellular water. Animal models more specific to the types of mild traumatic brain injury typically incurred by humans are needed to confirm the histological correlates of these macroscopic markers of white matter pathology.

  • mild traumatic brain injury
  • DTI
  • longitudinal

Introduction

In spite of the large number of individuals who sustain a mild traumatic brain injury (TBI) each year (Faul et al., 2010), the pathophysiology of mild TBI and its role in producing neurobehavioural symptomatology is still actively debated. Since complex cognitive processes such as attention, executive functions and memory are dependent on intact white matter tracts between frontal, parietal and medial temporal lobes (Grieve et al., 2007), putative white matter pathology provides an intriguing link to the neurobehavioural symptomatology that is frequently observed during the semi-acute stage (i.e. first weeks) of mild TBI (Belanger et al., 2007; Bigler, 2008). Indeed, accumulating evidence from both animal (Mac Donald et al., 2007a; Spain et al., 2010) and human (Blumbergs et al., 1994; Bigler, 2010; Niogi and Mukherjee, 2010) studies suggest subtle white matter abnormalities that are better captured by diffusion tensor imaging (DTI) than by conventional imaging. However, the direction of tissue water diffusivity abnormalities (as measured by DTI) has proven difficult to characterize (Niogi and Mukherjee, 2010). Specifically, different studies report increased (Bazarian et al., 2007; Wilde et al., 2008; Mayer et al., 2010; Henry et al., 2011), reduced (Arfanakis et al., 2002; Inglese et al., 2005; Miles et al., 2008; Smits et al., 2011) or unchanged (Zhang et al., 2010; Messe et al., 2011) white matter diffusivity during semi-acute mild TBI.

A commonly applied DTI metric is fractional anisotropy, which indirectly measures preferential water diffusion along white matter tracts and thus serves as a surrogate marker of tissue integrity. Animal studies of healthy, non-myelinated fibres and computer simulations suggest a more prominent role for axonal membranes in determining anisotropic diffusion rather than myelin (Beaulieu, 2002), with the ratio of intracellular to extracellular fluid also a contributing factor (Peled, 2007). Moreover, axonal pathology is more pronounced in the acute phase of injury (Mac Donald et al., 2007a; Dikranian et al., 2008; Spain et al., 2010), with evidence of injury progression over a 4–6-week period that correlates with cognitive dysfunction (Spain et al., 2010). These traumatic axonal injuries can occur either as a direct result of rotational forces or from complex secondary cellular processes involving mechanoreceptor dysregulation, ionic flux and alterations to the cytoskeleton (Albensi et al., 2000; Povlishock and Katz, 2005; Van Putten et al., 2005).

In contrast to animal work, differences in injury severity, injury type, measurement time post-injury, use of representative versus persistently symptomatic samples, and sample type (e.g. athlete versus military versus typical emergency room patient) all potentially contribute to the heterogeneity of findings in human diffusivity studies. For example, studies examining patients with chronic (months to years post-injury) mild TBI generally report decreased fractional anisotropy (Niogi and Mukherjee, 2010). However, patients with chronically symptomatic mild TBI reflect only a minority (3–20%) of the overall mild TBI population and thus represent a specific sample with potentially different clinical characteristics (Dikmen et al., 2001; Iverson, 2005). Similarly, sample-specific attributes (athletes or combat-exposed military personnel) may increase the likelihood of temporally proximal repetitive concussive or sub-concussive injuries (Zhang et al., 2006; Levin et al., 2010; Mac Donald et al., 2011), as well as introduce differences in general health and other complicating psychiatric factors (Hoge et al., 2008). As a result, directly recruiting sequential patients from the emergency room and scanning during the semi-acute injury phase may be the best strategy for collecting a sample that best generalizes to the overall population of mild brain injuries (Mayer et al., 2010).

As recently demonstrated with other neuroimaging modalities (Vagnozzi et al., 2010; Mayer et al., 2011; Yeo et al., 2011), imaging studies with longitudinal components are also critical for clarifying the recovery course of diffusivity abnormalities and the resolution of clinical symptoms. Animal models suggest that anisotropic measures are a good predictor of time post-injury (Mac Donald et al., 2007b), but these findings have not been easily replicated in human ‘cross-sectional’ studies (Inglese et al., 2005; Rutgers et al., 2008). Preliminary longitudinal work by our group (Mayer et al., 2010) and others (Arfanakis et al., 2002) suggests partial normalization of diffusion indices, although in each study sample sizes were small (n < 10) and the direction of ‘normalization’ varied. A recent longitudinal study of a military sample that included both semi-acute and post-semi-acute patients reported that while lower relative anisotropy in blast-related TBI persisted into follow-up, initial mean and axial diffusivity scores normalized across a 6–12-month interval (Mac Donald et al., 2011).

In addition to variability introduced by clinical factors, the analytical approach adopted for quantifying white matter abnormalities may also contribute to discrepancies in semi-acute mild TBI DTI studies. Specifically, both region of interest and voxel-based analyses have previously been used to characterize white matter pathology during semi-acute mild TBI. Voxel-based analyses emphasizing the alignment of white matter tracts [tract-based spatial statistics (TBSS)] have also been introduced. In general, both region of interest and voxel-wise analyses are based on the implicit assumption that clinically heterogeneous patients have a homogenous (i.e. high degree of spatial overlap) pattern of white matter abnormalities, although this assumption is much stronger for voxel-wise approaches. However, the validity of the spatial homogeneity assumption is questionable, which motivated us in the current study to adopt a new approach for quantifying diffusion abnormalities through a metric similar to lesion-load (White et al., 2009). Specifically, clusters of abnormally high or low anisotropic diffusion were determined on a voxel-wise basis and then summed to represent total burden of distributed pathology.

The current study thus had three primary goals. First, we attempted to independently replicate our previous finding of increased fractional anisotropy in a new cohort of patients with mild TBI (n = 28) using identical experimental protocols. Replication efforts are critical for clarifying the distinction between sample characteristics discussed earlier and the reproducibility of the methodologies used to examine their effects, especially given that previously published studies on semi-acute mild TBI have included relatively modest sample sizes (5 ≤ n ≤ 22). Secondly, data were combined across our original and replication samples (hereafter referred to as the pooled sample) to examine questions of recovery of function (through both cross-sectional and longitudinal analyses) and whether differences in fractional anisotropy were associated with indices of injury severity or type. Finally, a new metric for quantifying diffusion abnormalities that is less reliant on common injury sites across patients with mild TBI was introduced. Based on prior results, we predicted that patients with mild TBI in the replication sample would again show increased anisotropy in region of interest but not voxel-based analyses, and that there would be evidence of recovery of function at both short-term (i.e. a negative correlation with days post-injury in cross-sectional analyses) and longitudinal analyses.

Materials and Methods

Participants and clinical data: replication sample

A total of 29 patients with mild TBI and 29 gender-, age- and education-matched healthy controls were recruited through a case-control design. Specifically, the majority of controls were directly recruited to match patients in terms of gender, age (± 2 years) and education (± 2 years). One patient and one control were eliminated during quality assurance protocols (see following sections) leaving a total of 28 patients with mild TBI (16 males, 12 females; 28.21 ± 10.58 years old; 13.11 ± 2.04 years of education) and 28 healthy controls (27.36 ± 10.22 years old; 13.86 ± 1.63 years of education). All patients were recruited from local Albuquerque Emergency Rooms and were evaluated clinically (clinical exam = 15.71 ± 4.45 days post-injury) and with neuroimaging (imaging exam = 15.57 ± 4.29 days post-injury) within 21 days of injury. Clinical and imaging protocols occurred within 3 days of each other for the majority (82%) of patients (see Supplementary Table 1).

Inclusion criteria for the mild TBI sample were based on the American Congress of Rehabilitation Medicine [Glasgow Coma Score of 13–15, loss of consciousness (if present) < 30 min, post-traumatic amnesia (if present) < 24 h]. At minimum, all participants experienced an alteration in mental status, although the majority of the sample (18/28) also experienced a loss of consciousness. Mild TBI participants and controls were excluded if there was a history of neurological disease, psychiatric disturbance, additional closed-head injuries with >5 min loss of consciousness, a head injury within the last year, learning disorder, ADHD, or a history of substance or alcohol abuse. Informed consent was obtained from all participants according to institutional guidelines at the University of New Mexico.

Our clinical protocol was identical to our previous publication (please see Mayer et al., 2010, for details), and included a full neuropsychological and symptom assessment. This included measures of effort (Test of Memory Malingering), premorbid intelligence (Wechsler Test of Adult Reading), attention, working memory, processing speed, executive function and memory (composite indices). Self-reported symptomatology was also assessed with the Neurobehavioural Symptom Checklist, the Beck Depression Inventory—Second Edition and the State-Trait Anxiety Index. A modified Rivermead scale was also administered to capture current and previous head injuries and associated symptomatology. The neuropsychological data from one healthy control were excluded secondary to potential confounds (practice effects from an unrelated experiment). Additionally, one healthy control had abnormally low scores on the Test of Memory Malingering (T-score < 0) in conjunction with normal performance on other more demanding neuropsychological measures. Therefore, this subject’s Test of Memory Malingering score was excluded from further analyses.

Quality assurance

Data from the replication sample were combined with the data from our original sample (22 patients with mild TBI; 21 healthy controls) along with an additional unmatched healthy control to create a total sample of 51 mild TBI and 51 healthy controls. As head motion also introduces bias among DTI scalar measurements (Storey et al., 2007; Ling et al., 2012), quality assurance measures were also performed on all subjects (original and replication samples). One patient with mild TBI and one healthy control from the replication sample were identified as motion outliers [3 standard deviations (SD) greater than their cohorts] and were therefore eliminated from all (replication and the pooled) DTI analyses. There were no differences in head motion between mild TBI and healthy control groups in the replication (all P > 0.10) or original cohort (Ling et al., 2012). For the complete quality assurance protocol and results, see Supplementary material.

The current protocol did not include phantom data collection. Therefore, the number of days between each subject’s first and second scans in the study was used as a covariate in all DTI analyses, comparing original and replication cohorts to control for any potential changes in scanner drift.

Image acquisition and statistical analyses

Image acquisition was also identical to the previous protocol (Mayer et al., 2010). High-resolution T1-weighted anatomic images were acquired with a 5-echo multi-echo MPRAGE sequence (voxel size = 1.00 × 1.00 × 1.00 mm). T2-weighted images were collected with a fast spin echo sequence (voxel size = 1.15 × 1.15 × 1.5 mm). Two DTI scans (b = 800 s/mm2; voxel size = 2 × 2 × 2 mm) were acquired using a twice-refocused spin echo sequence with 30 diffusion gradients and the b = 0 experiment repeated five times.

DTI image analyses were performed using a mixture of AFNI and FSL software, and analyses were identical to our previous publication with the following best-practice updates: (i) the two DTI image data sets were concatenated prior to estimating the tensor; (ii) images were co-registered using an alignment of 12 rather than six degrees of freedom and (iii) the rotational effects of this alignment were applied to the gradient table. Both the original and replication samples were analysed with this new pipeline for the pooled data. Region of interest analyses were performed for the genu, splenium and body of the corpus callosum, as well as the superior longitudinal fasciculus, the anterior and posterior corona radiata, the superior corona radiata, the uncinate fasciculus and the anterior and posterior limbs of the internal capsule for both the right and left hemispheres based on the intersection of the Johns Hopkins University (JHU) white matter labels atlas from FSL (Mori and van Zijl, 2007) and each participant’s white matter tracts (i.e. a T1 image-segmented white matter mask). Multivariate tests were selected for region of interest analyses as a result of known correlations both within and across hemispheric white matter tracts.

In addition to region of interest analyses, two different voxel-wise analyses were performed. In the first analysis, each subject’s fractional anisotropy data were non-linearly aligned (using FNIRT from the FSL package) to a whole-brain fractional anisotropy template (FMRIB58_FA) and then blurred with a 6-mm full-width at half-maximum kernel prior to group analyses. To reduce the number of comparisons and restrict the analysis to white matter, all data were masked by the group fractional anisotropy mask (values ≥ 0.25). For the second voxel-wise analyses (TBSS), each subject’s fractional anisotropy data were aligned to a skeleton and all group-wise comparisons were also restricted to the skeleton (Smith et al., 2006). For both voxel-wise analyses, a permutation-based method (FSL Randomize) was used for group comparisons in conjunction with Threshold-Free Cluster Enhancement to correct for family-wise error (P < 0.05).

Finally, we also conducted a cluster analysis to capture the overall burden of white matter diffusion abnormalities, separately summing clusters of abnormally high and low anisotropy, regardless of specific location (White et al., 2009). Specifically, the mean and SD of fractional anisotropy was first calculated for each voxel from the spatially normalized (whole-brain fractional anisotropy template) pooled sample of matched healthy controls (n = 50). Next, individual subject data (both mild TBI and matched healthy controls) were transformed to signed z-scores on a voxel-wise basis using the calculated means and standard deviations from the overall healthy control sample. Individual subject z-score maps were then thresholded based on 2 SD above (z ≥ 2) or below (z ≤ −2) the healthy control mean, and masked by both the JHU atlas and the subject’s fractional anisotropy (fractional anisotropy ≥ 0.25). To further reduce false positives, contiguous voxels were required to form a cluster with a minimum of 16 native voxels (0.128 ml).

To derive an index score of diffusion abnormalities similar to lesion load, the total number of abnormal clusters and the associated total voxel count (total volume across all clusters) was then calculated by summing across all white matter tracts separately for abnormally high or low diffusion criteria. The square root of both values was then derived to more closely approximate a normal distribution.

Results

Original versus replication cohorts

There were no significant differences between healthy controls from the original and replication samples (Supplementary Material). The patient replication sample was evaluated at more days post-injury for both imaging [3.62 days difference; t(1,48) = 2.55, P = 0.014] and clinical [3.81 days difference; t(1,47) = 2.81, P = 0.007] exams, but was otherwise similar.

Replication sample results

Neuropsychological and clinical measures

Summary results (major demographics and domain scores) are presented in Table 1. There were no differences (P > 0.10) between patients with mild TBI and helathy controls on major demographic variables, indicating that the matching procedure was effective. A total of 19/28 healthy controls and 21/28 patients with mild TBI completed the modified Rivermead questionnaire, with a higher proportion (χ2 = 4.61, P = 0.032) of patients with mild TBI (10/21) reporting previous mild injury/injuries relative to the healthy controls (3/19). Healthy controls performed worse than patients with mild TBI on a measure of effort [Test of memory Malingering; t(1,32.36) = −2.09, P = 0.045], although performance across both groups was in the average range or above (T-score = 47.5 for healthy controls). Independent samples t-test indicated that patients had significantly lower estimates of premorbid intelligence Wechsler Test of Adult Reading than controls [t(1,54) = 3.02, P = 0.004]. Given the known relationship between intelligence and other neuropsychological measures as well as diffusion metrics (Grieve et al., 2007), Wechsler Test of Adult Reading was used as a covariate in all subsequent analyses (replication and pooled).

View this table:
Table 1

Neuropsychological and clinical summary measures

ReplicationPooled
Mild TBIControls
Mean (SD)Mean (SD)Cohen’s d
Demographic
    Age28.2 (10.6)27.4 (10.2)0.051
    Education13.1 (2)13.9 (1.6)−0.362
    HQ81.2 (39.7)79.2 (44.7)−0.029
Neuropsychology
    Attentiona52.8 (4.38)53.61 (6.13)−0.234
    Memorya52.2 (6.96)51.96 (7.34)−0.189
    Working memorya52.13 (6.22)53.4 (5.47)−0.048
    Processing speeda43.37 (5.47)46.74 (7.07)−0.322
    Executive functiona50.41 (5.55)49.4 (4.37)−0.055
    WTAR50.2 (8.6)56.4 (6.5)−0.690
    TOMM54.9 (4.3)47.5 (17.1)0.429
Self-report
    Emotional47.8 (8.5)42.8 (6.2)0.787
    NBSI-Som8.6 (7.9)1.6 (2.5)1.226
    NBSI-Cog4.5 (3.9)1.5 (2.5)0.975
Days post-injury
    Imaging15.6 (4.3)N/AN/A
    Neuropsychology15.6 (4.2)N/AN/A
  • a Means, standard deviations and effect sizes for neuropsychological indices reported following correction for Wechsler Test of Adult Reading as covariate.

  • HQ = handedness quotient; N/A = not applicable; NBSI-Cog = Neurobehavioural Symptom Inventory cognitive complaints (Som = somatic complaints); TOMM = Test of Memory Malingering; WTAR = Wechsler Test of Adult Reading.

The multivariate effect of group was not significant for the multivariate analysis of covariance (MANCOVA) analyses of cognitive composites (P > 0.10), although the univariate effect for processing speed showed a trend for slowed processing speed in the mild TBI group [F(1,53) = 3.33, P = 0.074]. In contrast, a significant multivariate effect of group [F(3,52) = 6.50, P = 0.001] was observed during a MANCOVA examining self-report measures (Table 1), with univariate effects suggesting significantly greater complaints for patients with mild TBI in the domains of emotional [F(1,54) = 6.36, P = 0.015], cognitive [F(1,54) = 12.02, P = 0.001] and somatic [F(1,54) = 20.15, P = 0.000] symptoms. Supplementary Table 2 provides a complete listing of neuropsychological data for the replication sample.

Structural imaging data

Eighteen of 28 patients with mild TBI had a CT scan at the time of their emergency room visit, but none of the CT scans were deemed to contain trauma-related pathology by a non-blinded neuroradiologist. In addition, T1- and T2-weighted MRI images were reviewed by a neuroradiologist blinded to patient diagnosis. Three patients (negative CT scan) had a finding that may have been secondary to trauma. Therefore, there were a total of 3/28 patients with visible lesions using conventional neuroimaging techniques.

Diffusion tensor imaging scalar analyses

Three MANCOVAs were performed to examine group differences in fractional anisotropy within five left and right hemisphere tracts and three divisions of the corpus callosum (Fig. 1). Significant univariate results are also reported for comparison purposes with our original sample. In the replication sample, the multivariate effect of group was not significant for the corpus callosum or the right or left hemisphere white matter tracts (P > 0.10). Univariate results indicated a significant effect of group within the genu [F(1,53) = 5.28, P = 0.026, Cohen’s d = 0.67] of the corpus callosum, whereas non-significant trends were observed for the right corona radiata (F(1,53) = 3.88, P = 0.054; d = 0.57) and right internal capsule [F(1,53) = 3.87, P = 0.054; d = 0.57]. For all structures, examination of the means indicated that fractional anisotropy was elevated for the mild TBI group, with follow-up tests indicating significantly decreased radial diffusivity (RD) for subjects with mild TBI within the genu [F(1,53) = 4.60, P = 0.037, d = −0.38] with a non-significant trend in the right internal capsule [F(1,53) = 3.52, P = 0.066, d = −0.32]. There were no differences in these tracts for axial diffusivity (P > 0.10).

Figure 1

Fractional anisotropy (FA) values from all regions of interest. This figure presents the mean FA values from the replication cohort during their first visit. Fractional anisotropy values are corrected for differences in premorbid intelligence, patients with mild TBI (mTBI) represented by black bars and healthy controls (HC) by grey bars. Error bars represent the standard deviation of the sample. (A) Regions of interest include (A) the genu (GNU), body (BDY) and splenium (SPL) of the corpus callosum (CC), as well as the superior corona radiata (SCR), the superior longitudinal fasciculus (SLF), the uncinate fasciculus (UF), the corona radiata (CR), and the internal capsule (IC) from the left (B) and right (C) hemispheres. Significant effects are denoted with double asterisks, statistical trends with a single asterisk.

Similar to our previous supplementary analyses (Mayer et al., 2010), there were no significant group differences in fractional anisotropy based on either of the two voxel-wise analyses (TBSS or voxel-wise white matter comparison) following recommended correction for family wise error.

The total number of clusters with high (z ≥ 2) or low (z ≤ −2) anisotropic diffusion and associated total cluster volume were also directly compared across patients with mild TBI and controls with ANCOVAs. Results indicated non-significant trends for both an increased number of clusters with high anisotropy [F(1,53) = 3.87, P = 0.055, d = 0.58] as well as total cluster volume [F(1,53) = 3.64, P = 0.062, d = 0.56] for patients with mild TBI in the replication sample relative to controls (Fig. 2B). In contrast, there were no significant differences between the groups on either cluster measure corresponding to low anisotropy (P > 0.10). An identical series of analyses were also conducted with the original cohort data (Fig. 2A), with results showing similar effects in terms of both increased number of clusters with high anisotropy [F(1,40) = 4.06, P = 0.051, d = 0.64] as well as total cluster volume [F(1,40) = 4.68, P = 0.037, d = 0.69] and non-significant (P > 0.10) results for clusters with low anisotropy.

Figure 2

Cluster analysis. This figure shows the mean square root of the number of clusters and their associated volumes for which individual subject z-scores exhibited low (z ≤ −2 SD below the pooled healthy controls mean) or high (z ≥ 2 SD above the healthy controls mean) anisotropy. All values are corrected for differences in premorbid intelligence, with mild traumatic brain injury patients (mild TBI) represented by black bars and healthy controls (HC) by grey bars (error bar = SD). (A) Results for the original cohort of patients and controls, whereas (B) represents the replication cohort. Significant effects are denoted with double asterisks, statistical trends with a single asterisk.

To ensure that these findings were not solely related to our minimal cluster volume (0.128 ml), the data were reanalysed using spatial volumes ranging from 0.032 to 0.320 ml in incrementing steps of 0.016 ml (Fig. 3). However, results indicated direction and magnitude of increased anisotropic diffusion were generally consistent regardless of the minimal cluster volume threshold. A full range of spatial volumes is presented in Supplementary Fig. 1.

Figure 3

Cluster results as a function of voxel threshold. The y-axis represents the square root of the mean number (top) and mean volume (bottom panel) for clusters exhibiting high anisotropy (z ≥ 2.0) in patients with mild traumatic brain injury (mild TBI; black bars) and healthy controls (HC; grey bars). The x-axis represents different minimal cluster volume thresholds, ranging between 32 and 320 µl. Error bars represent the standard deviation. A full range of voxel thresholds (16–1008 µl) are shown in Supplementary Fig. 1.

Pooled sample results

Region of interest analyses

Pooled analyses (ANCOVAs) were restricted to summary results from the cluster analyses (cluster number/total cluster volume with high anisotropic diffusion) or region of interest that showed significant effects or trends in either the original (splenium and left superior corona radiata, uncinate fasciculus, internal capsule and corona radiata), replication (right internal capsule and corona radiata) or both original and replication cohorts (genu of corpus callosum). Significantly higher fractional anisotropy in the genu [F(1,97) = 12.02, P = 0.001, Cohen’s d = 0.74] and left corona radiata [F(1,97) = 4.80, P = 0.031, d = 0.47] was observed for patients with mild TBI relative to the healthy control group, with non-significant trends in the left superior corona radiata [F(1,97) = 3.16, P = 0.079, d = 0.38] and right corona radiata [F(1,97) = 2.92, P = 0.091, d = 0.36]. Additionally, the number of clusters exhibiting high anisotropy [F(1,97) = 6.41, P = 0.013, d = 0.54] as well as total cluster volume [F(1,97) = 6.63, P = 0.012, d = 0.55] was also significant for pooled sample results.

Effects of clinical characteristics

Several cross-sectional analyses were conducted to examine the effects of recovery (days post-injury), injury severity (based on American Academy of Neurology criteria) or type (motor vehicle accident versus falls/assaults/struck by object), and self-reported post-concussive symptoms for the cluster statistics (cluster number and total cluster volume) and four regions of interest (genu, left and right corona radiata, and left superior corona radiata) that showed significant or trend differences for fractional anisotropy in the pooled sample ANCOVAs.

The first series of analyses examined the correlation between days post-injury and fractional anisotropy. Although all six variables showed the predicted negative relationship between anisotropy and days post-injury (−0.200 ≤ r ≤ −0.022), none of the findings were statistically significant (P > 0.10). Moreover, voxel-wise analyses also indicated that there was no relationship between time post-injury and fractional anisotropy (Supplementary material).

To examine effects of injury severity on fractional anisotropy, we first reclassified the subjects with mild TBI into two groups based on American Academy of Neurology classification (37 mild TBI Grade 3; 13 Grade 1 or 2). The 13 Grade 1 and 2 subjects were then matched with 13 subjects from the Grade 3 group on age, gender, and education and compared. However, there were no significant differences (P > 0.10) between injury severity and measures of increased anisotropic diffusion in region of interest or cluster indices. A similar procedure was utilized to match the 13 patients with mild TBI from a motor vehicle accident to 13 mild TBI from other injury mechanisms (n = 37); however, these results were also negative (all P > 0.10).

Linear regressions were conducted to investigate potential relationships between self-reported concussive symptoms and intelligence on measures of increased anisotropic diffusion. The overall model was not significant for any of the regressions, although cognitive (t = 1.77, P = 0.084) and somatic (t = −1.88, P = 0.066) complaints were weakly associated with fractional anisotropy in the genu.

Finally, we attempted to objectively classify patients and controls using metrics of anisotropic diffusion through binary logistic regression. In the first model, estimates of premorbid intelligence were entered first followed by fractional anisotropy from the genu, bilateral corona radiata and left superior corona radiata. Results indicated that the estimate of premorbid intelligence (Wald = 9.86, P = 0.002) was significantly able to discriminate between healthy controls (66% accuracy) and patients with mild TBI (64%). The inclusion of fractional anisotropy from regions of interest significantly improved classification accuracy across both groups, although the increase in accuracy was marginal (healthy controls = 70%; mild TBI = 72%) and only significant for the genu (Wald = 7.35, P = 0.007). The second model utilized a similar hierarchical structure (premorbid intelligence entered first) with the cluster statistics (positive cluster number and total cluster volume) entered as a second step. However, neither metric significantly improved prediction of patient/control status (Wald P > 0.10).

Longitudinal comparisons

Twenty-six patients with mild TBI from the pooled sample returned for their longitudinal visit. Returning patients were more educated [t(1,48) = −2.40, P = 0.020] and performed better on the memory composite index [t(1,47) = −2.40, P = 0.020] relative to the non-returning cohort. There were no other significant differences in key demographic variables, clinical variables, cluster statistics or fractional anisotropy of the genu for returning and non-returning patients (P > 0.10). Twenty-six returning healthy controls (total returning controls = 42) were selected on gender, age and education to best match with returning patients, with no significant differences in key demographic variables across the groups.

Longitudinal analyses were first conducted to examine changes in self-reported symptomatology as a function of recovery. The main effects of group and time were significant across all composite indices (all P < 0.05). The Group × Time interaction was also significant for cognitive [F(1,49) = 7.84, P = 0.007], somatic [F(1,49) = 6.64, P = 0.013] and emotional [F(1,49) = 5.19, P = 0.027] complaints, with simple effects testing suggested that the level of post-concussive complaints decreased as a function of time (Visit 1 to Visit 2) within the mild TBI subjects (all P < 0.05) but not healthy controls (all P > 0.10) group.

Three 2 × 2 (Group × Time) mixed measures ANCOVAs were then performed to examine how the number of clusters with high anisotropy, total cluster volume and fractional anisotropy within the genu changed as a function of recovery (Fig. 4). Longitudinal analyses were restricted to these metrics as they exhibited reliable differences across both the original and replication cohorts. Results indicated that the main effects of group and time were not significant for any of the analyses (P > 0.10). A trend for the Group × Time interaction [F(1,49) = 3.68, P = 0.061] was observed for the total number of positive clusters. Simple effect testing indicated that whereas healthy control subjects exhibited no change in the total number of clusters with high anisotropy (P > 0.10), clusters with high anisotropy were significantly reduced at Visit 2 for the patients with mild TBI [t(1,25) = 2.40, P = 0.024]. Even though the omnibus interaction term was not significant, the total cluster volume also decreased in Visit 2 for patients with mild TBI [t(1,25) = 3.07, P = 0.005] but not controls (P > 0.10).

Figure 4

Longitudinal measures of high anisotropy. This figure presents longitudinal changes in anisotropy measures for patients with mild traumatic brain injury (mild TBI; black bars) and healthy controls (HC; grey bars). Solid bars indicate data from Visit 1 whereas bars with diagonals show Visit 2 data, with all error bars equivalent to the standard deviation. (A) Mean fractional anisotropy (FA) of the genu of the corpus callosum, whereas the number of clusters with high anisotropic diffusion (z ≥ 2.0) and their respective volumes are presented in (B) and (C), respectively. Significant effects are denoted with double asterisks.

Finally, a multiple regression was conducted to determine whether changes (Visit 1 to Visit 2) in the number of clusters with increased anisotropy was associated with self-reported cognitive, somatic or emotional symptomatology for the subjects with mild TBI. Although the overall model was not significant (P > 0.10), a reduction in severity of somatic complaints significantly predicted a reduction in the number of clusters with increased anisotropy (t = 2.33, P = 0.029).

Exploratory analyses

A primary aim of the current study was to replicate our previous findings of increased fractional anisotropy in a limited selection of regions of interest. However, several other white matter tracts have also been implicated in studies of both semi-acute and chronic mild TBI. Therefore, exploratory analyses were also conducted to compare fractional anisotropy from all 48 regions of interest in the JHU atlas across the patients and controls in the pooled sample. Results indicated significantly higher fractional anisotropy in the mild TBI group for 7/48 white matter tracts (uncorrected P < 0.05), with an additional 4/48 tracts showing a non-significant trend (uncorrected P < 0.10). Fractional anisotropy was increased for mild TBI subjects in varying degrees of magnitude for 32/37 remaining tracts, resulting in a significant signed test (P = 0.001). Sample means, standard deviations, effect sizes, uncorrected P-values and number of patients/controls with significant clusters of high anisotropy are presented for all 48 regions of interest (Supplementary material).

Discussion

A clear picture of the diffusion abnormalities associated with mild TBI has not yet emerged from human neuroimaging studies (Niogi and Mukherjee, 2010), perhaps because of different sample characteristics and the implicit assumption that white matter injuries will occur in the same location across patients whose injuries reflect diverse mechanical deformations. As noted above, previous studies utilizing region of interest techniques have reported both reduced fractional anisotropy in the internal capsule and splenium (Huisman et al., 2004), centrum semiovale, corpus callosum and internal capsule (Miles et al., 2008), as well as increased fractional anisotropy in the corpus callosum (Bazarian et al., 2007; Mayer et al., 2010; Henry et al., 2011), left superior corona radiata, left corona radiata and left uncinate fasciculus (Mayer et al., 2010) during the semi-acute stage of mild TBI. Similarly, findings have also been mixed for studies utilizing voxel-based techniques (Mayer et al., 2010; Cubon et al., 2011; Messe et al., 2011; Smits et al., 2011). Therefore, a principal goal of the current study was to replicate our previous finding of increased anisotropic diffusion (fractional anisotropy) and reduced radial diffusion (Mayer et al., 2010) in an independent cohort of patients using nearly identical clinical and experimental procedures.

Results from our replication cohort confirmed a general pattern of increased fractional anisotropy and reduced radial diffusivity during the semi-acute phase of mild TBI; however, the spatial location of these diffusion abnormalities varied across cohorts. Specifically, several left hemisphere tracts exhibited increased diffusion anisotropy in the original but not replication cohort, whereas two right hemisphere tracts exhibited trends for increased diffusion in the replication but not original cohort. Only the genu of the corpus callosum exhibited a significant effect across both cohorts, suggesting that this structure may serve as a common point of injury regardless of initial biomechanical factors. Approximately 80% of fibres in the corpus callosum are unmyelinated, rendering them more vulnerable to traumatic axonal injury (Reeves et al., 2005), which likely contributes to the frequent occurrence of histopathological abnormalities in this structure following TBI (Gentry et al., 1988). In addition, modelling studies indicate that the inter-hemispheric fibres of the corpus callosum may be more susceptible to mechanical strain as a result of brain deformation compared to long intra-hemispheric fibre bundles (Viano et al., 2005). Finally, individual subject-based finite element modelling has recently been used to show that mean and maximum strain rates directly correlate with changes in anisotropic diffusion in the corpus callosum of 10 recently concussed athletes (McAllister et al., 2012).

In contrast to our region of interest analyses, two separate voxel-wise comparisons (TBSS and a voxel-based white matter analysis) did not reveal any group differences in anisotropic diffusion following recommended corrections for false positives. Null findings in voxel-wise comparisons have been previously reported (Levin et al., 2010; Mayer et al., 2010; Zhang et al., 2010; Cubon et al., 2011; Messe et al., 2011), although others have reported significant group differences by using alternative thresholding strategies (Henry et al., 2011; Smits et al., 2011). These findings suggest that the putative combination of heterogeneity in injury location and strict corrections for reducing false positives may limit the utility of voxel-based approaches in mild TBI research. In addition, known differences in individual white matter morphology must be corrected through spatial smoothing and complex registration algorithms, both of which can introduce variability in the data (Jones et al., 2005; Van Hecke et al., 2009; Moraschi et al., 2010). TBSS (Smith et al., 2006) attempts to mitigate these two issues through registration to a tract-based skeleton, but in doing so limits the number of voxels examined to the core of white tracts. This may be particularly detrimental for detecting abnormalities in grey-white matter boundaries, where pathology is frequently secondary to differences in tissue rigidity (Ducreux et al., 2005; Le et al., 2005).

From a methodological perspective, region of interest analyses somewhat relax the assumption of spatial homogeneity of injury by smoothing (i.e. averaging) over a priori defined white matter tracts. However, the manual selection of regions of interest can affect the reliability and reproducibility of results as well as be quite time-consuming. Our region of interest methodology mitigates the first two concerns by utilizing an automatic atlas-driven approach and reduces partial-voluming effects (intersection with individual subject’s white matter mask). A larger concern is that region of interest analyses may reduce the sensitivity for detecting the relatively sporadic lesions that exist even within tracts exposed to similar deformation forces (Browne et al., 2011; Greer et al., 2011).

To overcome the respective limitations of both region of interest and voxel-based methods, an approach for identifying diffuse traumatic axonal injuries on a voxel-wise basis was developed (White et al., 2009). Specifically, we defined individual clusters of both increased and decreased diffusion anisotropy based on our healthy control data to more accurately model a prediction of spatially heterogeneous diffuse lesions. Results indicated trends for a larger number of clusters with increased diffusion and greater total volume of these clusters for patients with mild TBI across both the original and replication cohorts. These findings are similar to a co-occurring publication demonstrating spatially heterogeneous white matter injuries following blast injury (Davenport et al., 2012). However, the Davenport and colleagues report regions of reduced anisotropic diffusion, potentially reflecting a relatively more chronic injury. Exploratory analyses confirmed that 42 out of the 48 regions of interest showed increased anisotropy for patients with mild TBI relative to controls to some degree, with six tracts exhibiting a moderate or greater effect size (Cohen’s d ≥ 0.40). Thus, current results provide evidence of diffusion abnormalities in several other tracts that may be targeted for future studies.

As previously discussed by our group and others (Bazarian et al., 2007; Wilde et al., 2008; Mayer et al., 2010), increased fractional anisotropy in white matter may be reflective of changes in the concentrations of intra and extracellular water (cytotoxic oedema). Specifically, the mechanical forces of mild TBI result in the stretching of axons and related supporting structures such as oligodendrocytes (Povlishock and Katz, 2005; Browne et al., 2011), altering the function of gated ion channels, and resulting in increased intracellular and decreased extracellular water (Rosenblum, 2007). The decrease in extracellular water may lead to a decrease in diffusivity perpendicular to the axon (radial diffusivity), secondary to more tightly compacted axons and the resultant changes in tortuosity of water (Sotak, 2002; Rosenblum, 2007). Other potential physiological bases for increased fractional anisotropy include changes in myelin structure. However, myelin plays a relatively minor role in anisotropic diffusion in comparison with axonal membranes (Beaulieu, 2002) and animal models of concussion more consistently indicate axonal rather than myelin pathology (Spain et al., 2010).

Human diffusion studies (i.e. DTI) are also influenced by a variety of other non-trauma-related factors including head motion (Ling et al., 2012), individual differences such as intelligence level (Grieve et al., 2007), smoking history (Gons et al., 2011) and emotional status (Korgaonkar et al., 2011), all factors that likely account for the variability in anisotropic diffusion that is observed within both patient and control groups in many studies (Ling et al., 2012). The contributions of some of these factors can partially be mitigated through careful selection of appropriate control populations through careful case-control. For example, in the current study all controls were directly matched on gender, and 47/51 controls were matched to within 2 years (age and education) of their respective patients. As discussed in the ‘Introduction’ section, sample characteristics also likely contribute to the differences in findings that are observed across different studies of mild TBI (Niogi and Mukherjee, 2010). Patients with positive CT scans (complicated mild TBI) tend to experience more neurobehavioural symptoms and poorer prognoses relative to non-complicated patients (Iverson, 2006; Lange et al., 2009). Future studies are needed to determine whether poor outcome is related to the presence of a CT detectable lesion or may result from more diffuse injuries such as those reported in the current study.

DTI also measures underlying pathology on a relatively macroscopic level (i.e. standard voxel volume of 0.008 ml) and there are several morphological changes, metabolic processes and inflammatory responses that occur following TBI on a cellular level (Albensi et al., 2000; Van Putten et al., 2005; Barkhoudarian et al., 2011). Therefore, it is a considerable challenge to directly relate macroscopic measures of diffusion (DTI) with cell histology. The current findings of increased anisotropic diffusivity appear to be inconsistent with animal models of TBI, which have consistently indicated reduced anisotropic white matter water diffusion in the acute and semi-acute injury stages (Mac Donald et al., 2007a, b; Budde et al., 2011; van de Looij et al., 2012). Although other animal studies of blast (Rubovitch et al., 2011) and cortical impact injuries (Budde et al., 2011; Xu et al., 2011) indicate ‘increased’ fractional anisotropy, these changes are more prevalent in grey rather than white matter regions. It is important to note that animal injury models frequently induce cortical contusions or other alterations (e.g. enlarged ventricles) of sufficient severity that they are visible with magnetic resonance, and these lesions are rarely observed in human mild TBI (Hughes et al., 2004). This suggests the need for animal models that more closely parallel the types of injuries typically observed in human mild TBI (Rubovitch et al., 2011).

For example, a recent study carefully attempted to replicate the mechanical forces necessary to result in axonal deformations that are typically present in human mild TBI in a swine model (Browne et al., 2011). Although DTI metrics were not obtained, histopathological results indicated widespread axonal swelling secondary to a disruption in axonal transport following mild TBI. It has been demonstrated that fast axonal transport and associated microtubule structure have little influence on anisotropic diffusion after pharmacological lesions (Beaulieu and Allen, 1994); however, the pathology of mechanical trauma is different from pharmacological injury. Seminal studies of Weiss and Hiscoe (1948) suggest that axonal constriction results in disrupted slow axonal transport, with accumulating neurofilament proteins and axoplasm ultimately resulting in cell distension upstream of the injury site (cited in Koehnle and Brown, 1999; Miller and Heidemann, 2008). Of note, cell distension should theoretically result in increased radial diffusivity (and thus reduced fractional anisotropy) ‘within a single axon’. It is not clear how axonal swelling affects macroscopic diffusion measurements that may be influenced by partial compression of adjacent axons and the reduction of extracellular fluid. Clearly, more research is needed to understand disruptions in slow axonal transport and how it relates to the macroscopic measures of diffusion used in human studies. This may be especially relevant given the emerging literature on the role of disrupted slow axonal transport in neurodegenerative disorders and tauopathies (Roy et al., 2005), with equivalent findings of tauopathies in human histopathological studies of repetitive mild TBI (McKee et al., 2009).

A critical goal of all biomarker research is to provide an objective, sensitive and specific measure that is associated with symptomatology following injury. Similar to our previous finding (Mayer et al., 2010), we were able to objectively classify mild TBI patients from controls based on measures of anisotropic diffusion from the genu of the corpus callosum in our pooled sample. However, the sensitivity and specificity of classification were relatively modest (healthy controls = 70%; mild TBI = 72%), likely a result of the considerable variability that was present across both patient and control groups. Finally, the relationship between our metrics of increased anisotropic diffusion (both cluster metrics and region of interest) and clinical measures of injury severity (based on American Academy of Neurology), injury type and post-concussive symptomatology were relatively weak. Collectively, current and previous (Niogi and Mukherjee, 2010) results suggest that although diffusion abnormalities may provide unique objective information about mild TBI that is not available through other standard assessment techniques, by itself it may only provide a weak biomarker of injury.

The second major goal of the current study was to examine how diffusion abnormalities changed as patients transitioned from the semi-acute to more chronic injury stages. Previous cross-sectional studies examining the evolution of anisotropic diffusion across multi-year periods have not reported significant differences between semi-acute and patients with chronic mild TBI (Inglese et al., 2005). Similarly, our moderately sized cohort of 50 patients did not exhibit any evidence of changes in increased anisotropic diffusion during the semi-acute injury phase (2–21-day post-injury). The lack of correlation with days post-injury suggests that white matter abnormalities may not evolve as quickly as alterations in neurometabolic functioning (Vagnozzi et al., 2010; Yeo et al., 2011), and hence, may represent part of the neuroanatomical substrate for the repetitive injuries.

The current study also incorporated a more powerful longitudinal design to examine the recovery of diffusion abnormalities. In our initial study we showed evidence for a partial normalization of fractional anisotropy within the splenium and corona radiata for the mild TBI group across a 4–5-month interval. Other studies have reported a continued reduction in anisotropic values during longitudinal studies in more severely injured (Bendlin et al., 2008; Kumar et al., 2009) and mild TBI (Arfanakis et al., 2002; Mac Donald et al., 2011) populations. The current study provided mixed evidence of recovery in a cohort of 26 patients, with the number of clusters and total volume of voxels with increased diffusion decreasing as patients recovered from their injury. Importantly, these reductions were also associated with a reduction in self-reported post-concussive symptomatology (somatic complaints). However, there were no significant changes in fractional anisotropy values of the genu, indicating that measuring the longitudinal recovery of anisotropic diffusion may be challenging in the presence of considerable individual heterogeneity.

The current study had several limitations. First, despite the implementation of several protocols to maximize participant retention (i.e. regular telephone and mail contact; alternative contact information), a significant number of patients did not return for their follow-up visit. The primary reasons for attrition were inability to contact patients (71%), patient no-shows (21%) and inability of patient to complete second visit secondary to personal circumstances (8%). Our attrition rate is similar to previously published TBI studies (Jorge et al., 2004; Rapoport et al., 2006; Jacobs et al., 2010) and there were no differences between returning and non-returning patients on key demographic, clinical and imaging variables. However, patient attrition may have biased the sample to a more or less symptomatic population. Secondly, our anatomical imaging protocol was rather limited for classifying lesions. Future studies should include a more detailed anatomical protocol (e.g. susceptibility weighted imaging) to better characterize additional lesion types such as petechial haemorrhages that are also seen in mild TBI. Thirdly, we did not assess patients for litigation, although the rate of litigation is likely to be very low among patients with mild TBI prospectively recruited from the emergency room.

In conclusion, we independently replicated our previous finding of increased fractional anisotropy and reduced radial diffusivity during the semi-acute stage of mild TBI in the genu of the corpus callosum. The corpus callosum has been previously established as a common point of trauma, and may be more susceptible to pathology regardless of actual injury focus secondary to anatomical (high frequency of unmyelinated fibres) and biomechanical factors. In contrast, the lack of replication between our previous (Mayer et al., 2010) and our current results across several other white matter tracts argues forcefully for the heterogeneous nature of diffusion abnormalities, potentially resulting from the unique circumstances that characterize each injury. This heterogeneity can be partially overcome by using methods that minimize the assumption of spatial homogeneity of injury, resulting in a finding of increased anisotropic diffusion across both our original and replication samples.

Funding

National Institutes of Health (R24-HD050836, R21-NS064464-01A1 and 3 R21 NS064464-01S1 to A.M.).

Supplementary material

Supplementary material is available at Brain online.

Acknowledgements

Special thanks to Diana South and Cathy Smith for assistance with data collection.

Abbreviations
DTI
diffusion tensor imaging
TBI
traumatic brain injury
TBSS
tract-based spatial statistics

References

View Abstract