Quantified MRI and cognition in TBI with diffuse and focal damage☆

In patients with chronic-phase traumatic brain injury (TBI), structural MRI is readily attainable and provides rich anatomical information, yet the relationship between whole-brain structural MRI measures and neurocognitive outcome is relatively unexplored and can be complicated by the presence of combined focal and diffuse injury. In this study, sixty-three patients spanning the full range of TBI severity received high-resolution structural MRI concurrent with neuropsychological testing. Multivariate statistical analysis assessed covariance patterns between volumes of grey matter, white matter, and sulcal/subdural and ventricular CSF across 38 brain regions and neuropsychological test performance. Patients with diffuse and diffuse + focal injury were analyzed both separately and together. Tests of speeded attention, working memory, and verbal learning and memory robustly covaried with a distributed pattern of volume loss over temporal, ventromedial prefrontal, right parietal regions, and cingulate regions. This pattern was modulated by the presence of large focal lesions, but held even when analyses were restricted to those with diffuse injury. Effects were most consistently observed within grey matter. Relative to regional brain volumetric data, clinically defined injury severity (depth of coma at time of injury) showed only weak relation to neuropsychological outcome. The results showed that neuropsychological test performance in patients with TBI is related to a distributed pattern of volume loss in regions mediating mnemonic and attentional processing. This relationship holds for patients with and without focal lesions, indicating that diffuse injury alone is sufficient to cause significant neuropsychological disability in relation to regional volume loss. Quantified structural brain imaging data provides a highly sensitive index of brain integrity that is related to cognitive functioning in chronic phase TBI.

Introduction

Disability due to traumatic brain injury (TBI) is estimated to affect 2% of U.S. population (Thurman et al., 1999). While the neurocognitive effects of TBI on speeded attention, memory, and executive functioning are well known (Levin et al., 1982), the nature of the relationship between TBI neuropathology and neuropsychological outcome remains unclear. Certain brain–behavior relationships have been described for focal lesions in ventral frontal and anterior temporal regions (Mattson and Levin, 1990; Stuss and Gow, 1992), yet the primary neuropathology of TBI is diffuse axonal injury (DAI) (Gennarelli et al., 1982; Povlishock and Katz, 2005). Relative to focal lesion effects, less is known about the effects of DAI on behavior, and it is not known how these effects interact with focal lesions.

Quantified structural MRI can reveal volume loss not appreciated in clinical interpretation. Yet the volume of specific structures visualized on MRI (e.g., the corpus callosum, fornix, and hippocampus) is only inconsistently correlated with neuropsychological outcome (Gale et al., 1993; Serra-Grabulosa et al., 2005; Tomaiuolo et al., 2004; Yount et al., 2002), possibly because additional damage outside these regions of interest is unaccounted for. Whole brain analysis with voxel-based morphometry has been applied in TBI, with variable relationships to behavior (Bendlin et al., 2008; Gale et al., 2005; Tomaiuolo et al., 2005; Warner et al., 2010a; Wilde et al., 2011) or with comparison to global outcome rather than to cognitive test performance (Sidaros et al., 2009; Warner et al., 2010b).

Many of the foregoing studies are complicated by small sample sizes, variance in imaging methods, and confounding the effects of focal and diffuse injury (Levine et al., 2006). Automated voxel-based methods pose significant problems with registration (Bookstein, 2001), especially given the distorted brain morphology in TBI (Kim et al., 2008). Using non-linear deformation field matching to preserve patients' anatomy, robust tissue compartment segmentation, regional parcellation of the cerebrum, and multivariate analysis, we previously found systematic relationships between quantified volume loss and TBI severity (Levine et al., 2008). Here we investigate the relationship of this volume loss to behavior as measured by concurrent neuropsychological testing.

Methods

Participants

Sixty-three TBI patients were recruited from consecutive admissions to a level I trauma center at approximately one year post-injury as part of the Toronto TBI Study (Levine et al., 2008). Injury severity was determined by the GCS at the time of discharge from the Trauma Unit (corresponding to the 6-hour GCS score; Teasdale and Jennett, 1974). Severity classification (mild/moderate/severe) was upgraded in eight cases where extended loss of consciousness (> 2 h), post-traumatic amnesia (> 48 h) or focal lesions suggested more severe injury than indicated by the GCS.

To minimize confounds from psychosocial effects specific to the TBI cohort, 27 comparison participants were recruited from friends and family members of the TBI patients. Exclusion criteria for all subjects included significant psychiatric history, substance abuse, learning disability, multiple traumatic brain injuries, severe injury to organs other than the brain, prior history of neurological disease, major medical problems or medications impacting cognition, and lack of English proficiency. As the patients were recruited from consecutive hospital admissions, compensation-related symptom exaggeration was minimized; all but one patient with large bifrontal lesions performed normally on a test of symptom validity (Green et al., 1996).

As seen in Table 1, participants in different groups were well matched for sex, age, education, vocabulary knowledge (Zachary, 1986), and time since injury. All participants were informed of the experimental aim of the study and gave their consent to participate.

Table 1

Characteristics of TBI patients and controls.

Group	Na	Age		Education (yrs)		Vocabularyb		GCSc		LOC (h)d		PTA (d)e		TSI (yr)f
		Mean	SD	Mean	SD	Mean	SD	Mean	SD	Median	IQR	Median	IQR	Mean	SD
Mild	13 (8)	33.7	13.1	13.2	2.0	28.8	3.6	14.6	0.7	0.00	0–0.13	0.0	0–0.30	1.19	0.42
Moderate	29 (15)	32.2	11.2	15.0	2.3	29.6	6.2	11.0	2.1	26.0	2–144	10.0	4.5–21	1.12	0.38
Severe	21 (14)	28.5	8.4	14.6	2.7	27.8	5.2	6.3	2.5	122	63–336	36	18–60	1.03	0.2
Controls	27 (11)	27.7	7.6	15.1	1.8	30.4	4.0	n/a	n/a	n/a	n/a	n/a	n/a	n/a	n/a

Number of males in parentheses.

Raw score on the vocabulary subtest of the Shipley Institute of Living Scale (Zachary, 1986).

Glasgow Coma Scale score (Teasdale and Jennett, 1974).

Duration of loss of consciousness.

Duration of post-traumatic amnesia.

Time since injury.

Neuropsychological tests

All participants received a battery of standardized neuropsychological tests sensitive to information processing deficits in attention, memory, and executive functioning caused by TBI. These included the Symbol-Digit Modalities Test (oral and written versions, SDMTo and SDMTw) (Smith, 1978), assessing clerical speed and accuracy; the Trail-Making Test, parts A and B (TRa, TRb) (Army Individual Test Battery, 1944), a paper-and-pencil test of speeded information processing; the Self-Ordered Pointing Test (six, eight, ten, and twelve item forms, SOP6, SOP8, SOP10, and SOP12) (Petrides and Milner, 1982), a measure of working memory; the Hopkins Verbal Learning Test—Revised (Benedict et al., 1998), involving verbal learning acquisition (HVLTln) and delayed recall (HVLTdl), Phonemic Word List Generation (FAS) (Spreen and Strauss, 1998) requiring speeded word retrieval to phonemic cues; and the Wisconsin Card Sorting Test (categories, criterion perseverations, response perseverations, and set loss errors, WCSTc, WCSTpc, WCSTpr, and WCSTs) (Stuss et al., 2000), a test of executive and frontal lobe functioning. Neuropsychological test data were scaled so that high scores indicate intact performance.

MRI acquisition and image processing

TBI patients were scanned with a 1.5 T-MR system (Signa, CV/i hardware, LX software, General Electric Healthcare, Waukesha, WI) at the time of testing. A sagittal T1-weighted 3D volume technique produced 124 1.3 mm slices (TR/TE of 35/5 ms, flip angle of 35°, 1.0 NEX, and FOV of 22 cm, acquisition matrix: 256 × 256 X 124). Proton density and T2-weighted images with a slice thickness of 3 mm were obtained using an interleaved sequence (TR/TE of 3000/30, 80 ms, 0.5 NEX, and FOV of 22 cm, acquisition matrix: 256 × 256). For TBI patients, gradient echo T2 sequences with a slice thickness of 6 mm were obtained to emphasize hemosiderin deposits (TR/TE of 750/35 ms, flip angle of 20°, 2.0 NEX, and FOV of 22 cm, acquisition matrix: 256 × 256). For technical reasons, gradient echo images were unavailable for 21 patients.

The image-processing pipeline is described in our earlier publication (Levine et al., 2008) (see also Fujiwara et al., 2008; McKinnon et al., 2008; Rosenbaum et al., 2008). In brief, each patient's T1-weighted image was registered to a template brain (Woods et al., 1998a, 1998b), preserving the original size of the brain while standardizing the position and orientation. Template matching was accomplished via non-linear registration of T1-weighted images to the template image (Collins and Evans, 1997). Removal of non-brain tissue from the image incorporated thresholding information derived from the PD- and T2-weighted images, facilitating the distinction between dura mater and gray matter (Kovacevic et al., 2002), optimizing estimation of subdural CSF.

Focal cortical contusions, appearing on at least two slices with a minimal diameter of 3 mm were manually defined in the axial plane in 20 of the 63 patients following a ventral frontal/anterior temporal distribution typical for TBI (Courville, 1937; Gentry et al., 1988) (see Fig. 1). Scans were also clinically interpreted by a board-certified neuroradiologist specializing in TBI. In addition to the large focal lesions, evidence of DAI (e.g., hemosiderin deposits) was present for 91% of the moderate and severe TBI patients and 15% of the mild TBI patients. These clinically evident lesions represent only the regions where the confluence of damage is sufficient to be detected by the naked eye (Gentry, 1990) and are provided for reference as volumetric measures are considered to be more sensitive markers of diffuse injury (Levine et al., 2008).

Fig. 1

Distribution of focal lesions (contusions) in 20 TBI patients. Lesion tracings are projected on selected axial slices of a template brain derived from 12 healthy control subjects. The color scale indicates degree of lesion overlap across patients (max = 5).

The voxels on the T1-image were then classified as representing gray matter, white matter, or CSF using a robust automated tissue classification method (Kovacevic et al., 2002). A trained operator reclassified CSF-segmented voxels inside the ventricles, allowing for the separate assessment of ventricular and sulcal/subdural CSF (vCSF and sCSF, respectively). A modified Semi-Automated Brain Region Extraction (SABRE) (Dade et al., 2004) method was then used to create 38 ROIs on the template brain (19 per hemisphere; see Fig. 2). Non-linear deformation field matching of the template to individual images was used to customize these regions to fit each patient's brain anatomy, again preserving inter-individual topographical variability. Regional volumes were adjusted for total intracranial capacity using a regression-based method (Arndt et al., 1991).

Fig. 2

SABRE regional cortical divisions in axial and sagittal views. Abbreviations: LSF: lateral superior frontal, MSF: medial superior frontal, LMF: lateral middle frontal, MMF: medial middle frontal, LVF: lateral ventral frontal, MVF: medial ventral frontal, GCG: genual cingulate gyrus, ACG: anterior cingulate gyrus, MCG: middle cingulate gyrus, PCG: posterior cingulate gyrus, AT: anterior temporal, MT: medial temporal, PT: posterior temporal, O: occipital, ABGT: anterior basal ganglia/thalamus, PBGT: posterior basal ganglia/thalamus, EC: external capsule/corona radiata, IP: inferior parietal, SP: superior parietal.

Image analysis

Partial Least Squares (PLS), a flexible, widely used multivariate neuroimage analysis technique (McIntosh and Lobaugh, 2004; McIntosh et al., 1996) was used to identify patterns of regional volume loss related to performance on neuropsychological tests. PLS uses singular value decomposition (SVD) applied to the brain–behavior correlation matrix to identify latent variables [LVs] that indicate optimal relations between patterns of regional brain volume loss and the profile of test scores (McIntosh and Lobaugh, 2004; McIntosh et al., 1996). The statistical significance of each LV is assessed by 1500 permutation tests (Edgington, 1980). The stability of each brain region's contribution to the LV is determined through bootstrap resampling (subjects were resampled 500 times) (Wasserman and Bockenholt, 1989). Brain regions were considered reliable if they had a ratio of salience to tandard error (hereafter referred to as the bootstrap ratio, interpreted similar to a Z-score McIntosh and Lobaugh, 2004) greater than 3, corresponding to 99% confidence limits.

Two PLS analyses were conducted. The first set incorporated the full sample of 63 TBI patients, providing a global picture of brain–behavior relationships in this sample. The second set was restricted to 43 patients without traceable focal lesions, hereafter referred to as the “diffuse injury” group. Each analysis included CSF, grey matter, and white matter volumes. Additional analyses were conducted on separate vCSF and sCSF volumes. We also used PLS to assess the relationship of TBI presence and severity classification and test performance (i.e., without considering brain imaging data). This provided a comparison of injury severity measures to brain imaging data in their relationship to cognitive functioning.

Results

Relationship of regional brain volumes to test performance

The PLS analysis for the full sample indicated one significant latent variable for each of the full and diffuse injury samples, ps = .0007 and .0133, respectively, by permutation test. These latent variables accounted for 70% and 58% of the covariance between test scores and regional volumes, respectively.

Contribution of test scores to patterns of brain volume changes

Figs. 3 and 4 display the latent variables for the full sample and diffuse injury PLS analyses. The contribution of each test to the patterns of volume changes described below can be determined by examining the correlations and their 99% confidence intervals as plotted in the upper panels of Figs. 3 and 4. For the full sample, correlations for all test scores except those associated with the WCST were significantly different from 0, with the highest correlations noted for tests of speeded information processing (SDMTo, SDMTw, TRa, TRb), although tests of mnemonic function were also significant. This pattern held for the diffuse injury sample, except that SOP10 was no longer significant, and WCSTpc and WCSTpr were marginally significant.

Patterns of volumetric changes

The color-coded plots of bootstrap ratios in Figs. 3 and 4 indicate the patterns of CSF, white, and grey matter volume changes associated with the above-described pattern of test scores. For the full sample, significant bootstrap ratios for CSF volumes were noted in a right-lateralized pattern, including the right inferior parietal and posterior temporal regions. The left inferior parietal lobe was also significant. The middle cingulate region was significant, right greater than left, with additional involvement of the right anterior cingulate region. All basal ganglia/thalamic regions were significant. In the frontal lobes, the medial sector of the middle frontal region was significant. Within grey matter, the highest bootstrap ratio was observed in the right posterior cingulate region. The medial temporal lobes were significant bilaterally, as was the right posterior temporal lobe. The ventromedial frontal regions were significant bilaterally.

The regional pattern for the diffuse injury sample was limited to grey matter, where the highest bootstrap ratio was observed for the right posterior temporal lobe and the right anterior cingulate gyrus, with additional involvement of the right posterior cingulate gyrus. Other posterior regions included the medial temporal lobes bilaterally and the right inferior parietal lobe. Within frontal lobes, significant bootstrap ratios were observed for the lateral sectors of the middle frontal regions and the ventromedial sectors bilaterally.

Ventricular vs. subdural/sulcal CSF

As above, PLS analyses for the full sample and diffuse injury group yielded a single latent variable for each compartment. For the full sample, p values for permutation tests for sCSF and vCSF were 0.009 and 0, respectively, accounting for 73 and 96% of the covariance between imaging and neuropsychological data, respectively. For the diffuse injury group, these values were .03 and .003, respectively, with 61 and 94% of the covariance accounted for, respectively. Overall, vCSF was more closely related to behavior than was sCSF. For the full sample, all tests but the WCST measures, certain SOP measures, and FAS (vCSF only) contributed to the pattern of volume loss. All regions of vCSF contributed to the pattern, with peaks in the right medial and left posterior temporal lobe, but for sCSF only the right superior parietal and left inferior parietal regions were significant. For the diffuse injury group, all tests but the HVLTd, SOP, and WCST contributed to a generalized pattern of subtle sCSF volume changes, with no specific region of sCSF emerging as significant in this analysis. All tests but the HVLTln, certain SOP measures, and the WCST covaried with vCSF loss in a single region: the left posterior temporal lobe.

Relationship of TBI severity classification to test performance

PLS analysis assessing the relationship between TBI presence/severity and test performance for the full sample and non-injured comparison subjects indicated one significant latent variable. The permutation test yielded a p value of 0, indicating that none of the 1500 permuted orderings of observations produced a latent variable exceeding the value of the observed one. A single latent variable also emerged from the analysis of the diffuse injury sample, p = .0007 by permutation test. The latent variables accounted for 82% and 78% of the covariance between test scores and severity for the full sample and diffuse injury samples, respectively. For both samples, the non-injured and severe TBI groups contributed to the identified pattern of test scores; the mild and moderate TBI groups did not significantly contribute to the test score pattern. In other words, with respect to their test performance, patients with mild or moderate TBI could not be statistically distinguished from non-injured or severe TBI patients. By comparison, in this same sample, mild and moderate TBI were associated with significant volume loss relative to non-injured comparison subjects (Levine et al., 2008).

Relative to the brain imaging data, injury severity classification was related to a comparatively constrained profile of test performance. Bootstrap ratios for SDMTo, SDMTw, TMTb, HVLTln and HVLTdl were statistically significant for the full sample. The diffuse injury group showed a similar pattern, with SDMTo, HVLTln, HVLTdl, and SOP10 identified as significant.

Discussion

Understanding the nature of brain–behavior relationships in TBI has implications for both management of patients with TBI and theory. Clinical interpretation of brain MRI findings has been inconsistently related to behavior as measured by test performance (Levin et al., 1992; Scheid et al., 2006). While quantified MRI has yielded correlations between structural brain volumes and neuropsychological test performance (Bendlin et al., 2008; Gale et al., 1993; Tomaiuolo et al., 2004; Warner et al., 2010a; Wilde et al., 2011), these correlations can be elusive or lacking in neuroanatomical specificity (Serra-Grabulosa et al., 2005; Yount et al., 2002), possibly due to heterogeneity in patient ascertainment, sample size, test selection, image acquisition and analysis methods, and confounding focal and diffuse injury effects.

In this study, we investigated the relationship between structural brain integrity and cognition in a large sample of chronic-phase TBI patients drawn from consecutive admissions to a major trauma center. We found systematic relationships between quantified measures of regional brain volume and neuropsychological test performance, even when patients with focal cortical contusions were excluded. The relationships to brain volume were spatially distributed, indicating that interpretation of brain–behavior relationships in TBI requires whole-brain assessment. Regional brain volume measures were more closely related to behavior than indicators of acute consciousness disruption routinely used to classify TBI severity.

Whereas our previous study in this cohort related regional brain volumes to TBI severity (Levine et al., 2008), in this study we considered regional brain volumes expressing focal or diffuse lesion load as more precise measures of TBI severity effects at the chronic stage. Accordingly, these regional volumetric measures demonstrated a finer-grained relationship to concurrent neuropsychological test performance than did acute TBI severity measures, where patients with mild and moderate TBI did not significantly contribute to the pattern of test scores even though these groups were differentiated in terms of patterns of regional brain volume changes (Levine et al., 2008).

The measures most robustly associated with TBI severity and to regional brain volumes involved speeded attention, working memory, and verbal learning and memory, which covaried with frontal, temporal, and parietal volumes in patients with and without focal lesions. Grey matter volumes in the medial and posterior temporal lobe, the ventromedial prefrontal cortex, and the right inferior parietal lobe covaried with test performance. These regions are classically associated with learning, memory, attention, and self-regulation as well as complaints of working memory deficits. Cingulate regions also figured prominently in the pattern, particularly in the posterior region that is vulnerable in TBI (Gentry et al., 1988; Kim et al., 2008; Levine et al., 2008). This region, closely interconnected to the medial temporal lobes, frontal lobes, and anterior cingulate gyrus (Morris et al., 1999) is considered a connector hub region crucial to integration of information across the cortex (Hagmann et al., 2008; Raja Beharelle et al., 2012).

The oral version of the SDMT, which requires rapid shifting of visual attention and maintenance and manipulation of information on-line, is an excellent measure of such integrated brain function that is highly sensitive to brain disease in other patient populations (Christodoulou et al., 2003; Parmenter et al., 2007; Turken et al., 2008). This test showed the strongest relation to volume loss in both samples. On the other hand, a canonical measure of executive/frontal function, the WCST, did not significantly covary with regional volumetric measures in the full sample. This may reflect the fact that lesions in this group, concentrated in ventral frontal regions, are not in the critical dorsal regions sensitive to WCST (Stuss et al., 2000). In the diffuse injury group, WCST performance subtly contributed to a pattern of test scores that covaried with lateral frontal grey matter (among other regions). Tasks sensitive to ventral frontal damage may show greater specificity in TBI. Indeed, in this same sample, we found a relationship between the tests of ventral frontal function, particularly the Smell Identification Test (Doty et al., 1984) and ventral frontal integrity (Fujiwara et al., 2008). Our data suggest that executive control deficits classically associated with TBI may be secondary to slowed information processing and impaired working memory, which were more robustly related to structural changes following TBI.

In relation to TBI severity, volume loss was demonstrated in grey matter, white matter, and CSF compartments (Levine et al., 2008). In relation to behavior, however, only CSF and grey matter contributed. Thus white matter volume loss is sensitive to TBI, yet the functional consequences as measured by the tests in this study are related to grey matter loss, which is assumed to reflect cell death due to delayed axotomy (Povlishock and Christman, 1995). Indeed, the contribution of white matter volume loss to increased CSF may have attenuated relationships between CSF and behavior. The lack of covariance between test behavior and white matter volumes may relate to limitations of functional localization within white matter as assessed by the T1-weighted image. This should not imply that white matter damage does not contribute to cognitive disability in TBI. On the contrary, axonal damage due to DAI leads to neuronal death and deafferentation of connected fields, affecting neuronal processing in select cortical regions (Maxwell et al., 2010; Povlishock and Christman, 1995). Accordingly, specific effects have been noted using diffusion tensor imaging (DTI) (Bendlin et al., 2008; Kinnunen et al., 2011; Wilde et al., 2006), which can show axonal damage effects in normally appearing white matter on the T1-weighted image. Nonetheless, DTI and the associated processing pipelines have yet to achieve the accessibility of standard T1-weighted imaging.

The relationship between vCSF and behavior in this study replicates prior work (Blatter et al., 1997). This finding, however, may appear difficult to reconcile with our finding that sCSF is more strongly related to injury severity and to whole brain grey matter volume (Levine et al., 2008) (see also Kim et al., 2008). This discrepancy may relate to the regional pattern of brain–behavior relationships that included midline regions with greater vCSF volume, particularly in the temporal lobes, which has been related to mnemonic deficits in patients with TBI (Gale et al., 1994).

Conclusions

TBI causes widely distributed brain damage with associated impairments on measures of integrated brain function mediating attentional, executive, and mnemonic capacities. In a large sample of patients with chronic-phase TBI drawn from consecutive admissions, variance in temporal, frontal, parietal, and cingulate regional volume covaried with performance on concurrently administered neuropsychological tests of speeded attention, working memory, and verbal learning and memory. Although both white and grey matter are affected in TBI, grey matter volumes more strongly covaried with behavior. The brain–behavior covariance patterns were modified by presence of focal lesions, yet regionally-specific effects held even when patients with focal lesions were excluded, illustrating that diffuse injury alone can cause significant neuropsychological disability. Quantified structural MRI provides a highly sensitive index of brain damage in TBI that is more closely related to behavior than are indicators of acute consciousness disruption routinely used to classify TBI patients into severity groups.