Regional brain morphometry predicts memory rehabilitation outcome after traumatic brain injury.

Cognitive deficits following traumatic brain injury (TBI) commonly include difficulties with memory, attention, and executive dysfunction. These deficits are amenable to cognitive rehabilitation, but optimally selecting rehabilitation programs for individual patients remains a challenge. Recent methods for quantifying regional brain morphometry allow for automated quantification of tissue volumes in numerous distinct brain structures. We hypothesized that such quantitative structural information could help identify individuals more or less likely to benefit from memory rehabilitation. Fifty individuals with TBI of all severities who reported having memory difficulties first underwent structural MRI scanning. They then participated in a 12 session memory rehabilitation program emphasizing internal memory strategies (I-MEMS). Primary outcome measures (HVLT, RBMT) were collected at the time of the MRI scan, immediately following therapy, and again at 1-month post-therapy. Regional brain volumes were used to predict outcome, adjusting for standard predictors (e.g., injury severity, age, education, pretest scores). We identified several brain regions that provided significant predictions of rehabilitation outcome, including the volume of the hippocampus, the lateral prefrontal cortex, the thalamus, and several subregions of the cingulate cortex. The prediction range of regional brain volumes were in some cases nearly equal in magnitude to prediction ranges provided by pretest scores on the outcome variable. We conclude that specific cerebral networks including these regions may contribute to learning during I-MEMS rehabilitation, and suggest that morphometric measures may provide substantial predictive value for rehabilitation outcome in other cognitive interventions as well.

Introduction

Traumatic brain injury (TBI) is a heterogeneous phenomenon, with a variety of causes, severities, and anatomical injury distributions. Despite this heterogeneity, TBI has been associated with a constellation of typical anatomical injuries as well as a prototypical set of functional deficits. Anatomical injuries typically include diffuse axonal injury (DAI), traumatic axonal injury (TAI), focal contusions in frontotemporal areas, and regions particularly vulnerable to excitotoxicity (Povlishock and Katz, 2005). Damage may be concentrated in acetylcholine dense pathways (Salmond et al., 2005), but has the potential for a wide spatial distribution (Merkley et al., 2008). In terms of long-term functional consequences of TBI, deficits are commonly observed in memory, attention, and in tasks relying heavily on executive functioning (Cicerone et al., 2000, 2006; McAllister et al., 2006). It is presumed that the consistency, albeit modest, of anatomical alterations partially accounts for the commonly observed functional deficits.

One approach to investigate the contribution of different anatomical regions to functional deficits is to utilize quantitative brain morphological measures (Bigler, 2001b). Such measures are made by dividing, or parcellating, brain tissue into separate regions using either manual processes (Kennedy et al., 1989) or computational ones (Fischl et al., 2002, 2004; Chen et al., 2004; Mega et al., 2005). Computational processes have the advantages of reducing subjectivity, analysis time, and cost, making them feasible to apply even on large cohorts of subjects. Importantly, computational approaches have exhibited accuracy similar to manual parcellation for many brain structures (Fischl et al., 2002, 2004; Desikan et al., 2009). One approach for parcellating brain tissue on magnetic resonance imaging (MRI) scans is to fit a mesh of points to the pial surface of the brain and a second mesh to the gray–white matter border. These two computed surfaces are then parcellated into regions (precentral gyrus, central sulcus, and so forth) based on regional curvature plus structural priors (Fischl et al., 2002; Han and Fischl, 2007). Regional measurements can then be made on the resulting parcellations, including regional cortical thickness, regional gray matter volume, and local surface curvature (Fjell et al., 2009; Jovicich et al., 2009).

Such morphometric measures have been used in studies of aging (Salat et al., 2004), finding regionally specific cortical thinning with age, and in pediatric TBI, finding relatively widespread cortical thinning post-injury (Merkley et al., 2008; Bigler et al., 2010). Semi-automatic parcellation techniques have demonstrated widespread gray matter volume loss in TBI (Levine et al., 2008). In adults, such volume loss has been further associated with injury severity, using both parcellation techniques (Levine et al., 2008) and voxel-based morphometry (Gale et al., 2005). Recent clinical studies have demonstrated an association between prospective memory and regional cortical thickness in children (McCauley et al., 2010), and an association between regional cortical thickness and behavioral performance in Alzheimer's disease (Dickerson et al., 2009). These results naturally lead to the hypothesis of a link between regional cortical gray matter volumes and TBI outcome.

Indeed, a number of studies have demonstrated a positive association between the extent of anatomical injury and general TBI outcome (Azouvi, 2000; Wilde et al., 2005; Ariza et al., 2006; Jorge et al., 2007; Kraus et al., 2007). In these, more extensive or severe injury to specific brain regions is associated with poorer performance on neuropsychological and cognitive tests. Less well understood, however, is the relationship between anatomical injury and outcome from any specific rehabilitation program, particularly in chronic TBI. While TBI rehabilitation programs can help improve memory and executive function (Cicerone et al., 2000), attention (Kim et al., 2009), community integration (Cicerone et al., 2004), and functional outcomes (Giles, 2001), little work has been done to predict outcomes from specific rehabilitation programs prior to their initiation. In a previous study of chronic TBI, we showed that functional brain activity may help predict memory rehabilitation outcome in this way (Strangman et al., 2008b). However, no studies have yet examined the ability of brain morphological measures to predict outcome from a TBI rehabilitation program, despite the widespread availability of structural MRI in clinical settings as compared to fMRI. We therefore sought to test whether regional volumetric measures derived from structural MRI scans could help predict outcomes from a TBI memory rehabilitation program.

The rehabilitation program we employed was developed as part of a recent study and was designed to improve participants’ memory by training individuals to use internally based semantic association and related strategies (O'Neil-Pirozzi et al., 2010). Prior neuroimaging research has demonstrated that such internal semantic and strategic cognitive manipulations typically involve dorsolateral prefrontal cortex (DLPFC), ventrolateral prefrontal cortex (VLPFC), posterior parietal cortex (PPC), and the hippocampal region (McAllister et al., 2001; Logan et al., 2002; Strangman et al., 2008a). We therefore hypothesized that decreased gray matter volumes in these regions would predict less favorable outcomes from our rehabilitation program. Given the complex nature of TBI, however, we also planned to examine the remaining brain regions from the parcellation process as an exploratory analysis.

Materials and Methods

Participants

Participants with TBI were recruited via mailings to clients of The Commonwealth of Massachusetts Brain Injury and Statewide Specialized Community Services Department, members of the Brain Injury Association of Massachusetts, local support groups, and to patients of study-affiliated physicians. Fifty-eight individuals were initially enrolled, following a procedure approved by the Human Research Committee at the Massachusetts General Hospital. All participants (1) were at least 18 years of age at the time of injury, (2) had sustained a TBI of any severity at least 12 months prior to the study, (3) were right-handed (Oldfield, 1971) and fluent in English, and (4) self-reported as having difficulty with memory following their injury. Exclusion criteria included a score of less than 4 on either the expression or comprehension items of the FIM, non-traumatic etiology of cerebral dysfunction in addition to TBI, active major illnesses, pre-injury history of psychiatric disease, inability to read single words at an eighth-grade reading level, and current drug or alcohol dependence (criteria from the Diagnostic and Statistical Manual of Mental Disorders, 4th Edition).

From this group, data was excluded from five individuals who failed to complete at least eight of 12 memory intervention sessions and/or return for posttesting. In addition, technical difficulties for two participants’ structural scans led to insufficient image quality to attempt parcellation. And in one participant, the cortical parcellation process could not run to completion, likely due to two particularly large focal regions of cortical tissue loss (approximately 8.95 and 3.98 cm3, respectively). We proceeded with analysis on the remaining 50 participants with TBI (14 female, 36 male). Demographic characteristics for these participants appear in Table 1.

Table 1.

Demographic characterization of study participants.

Characteristics	TBI Participants (n = 50)
	Mean (SD)	Min.	Max.
Age (years)	47.2 (11.4)	25.2	65.4
Education (years)	14.6 (2.2)	10	20
Time since injury (years)	11.5 (9.3)	1.3	37.6
Loss of consciousness (days)	12.0 (23.5)	0	105

Experimental procedures

Session 1 consisted of a series of standard neuropsychological tests of memory, executive function, and language. Outcome variables included the Hopkins Verbal Learning Test-Revised (HVLT-R) and the Rivermead Behavioral Memory Test II (RBMT-II). In the HVLT, a person is verbally presented a list of 12 words (four words from each of three semantically related groups), in a semantically unrelated order. The delayed correct recall score measures the ability to remember three presentations of these words after a 20-min delay. HVLT was a primary outcome because it is widely used, it emphasizes semantic association memory (a focus of the I-MEMS intervention), and it has multiple parallel versions to help minimize practice effects. RBMT test items include asking a patient to remember a hidden personal belonging, to remember to ask an appointment-related question, and to deliver a message. The RBMT is an ecologically valid, broad measure of impairment in everyday memory functioning. It has four parallel versions with good alternate form reliability. Other tests included the Boston Naming Test (BNT) 2nd Edition Short Form, the Boston Diagnostic Aphasia Examination (BDAE) 3rd Edition Animal Naming, Trail Making Test Part A and B, and Wechsler Memory Scale-Revised Digit Span (Forward and Backward). In Session 2, approximately 3 days later, each study participant underwent structural MRI scanning as described below.

Sessions 3 through 14 were memory intervention sessions led by one of the investigators that emphasized internal memory strategy training and were designed around evidence-based approaches; details are described in (O'Neil-Pirozzi et al., 2010). In brief, we utilized a group intervention with 3–6 members per group cycle and 2–3 group facilitators. Each group cycle ran twice weekly for 6 weeks, 90 min per session, for a total of 12 sessions. Sessions emphasized semantic organization and other internally based strategies (e.g., elaboration and imagery) from encoding, storage, and retrieval perspectives.

Session 15 (posttest 1) was conducted within 3 days after the final memory intervention session and repeated all of the tests from session 1, except for the BNT and BDAE. Session 16 (posttest 2) was conducted 1 month following session 15 and again repeated the same tests. No MRI scanning was performed at either posttest session. Forms of the various tests were counterbalanced across participants and sessions. Pre- and posttest evaluations were conducted by different individuals, with regular cross-validation of test administrators.

MRI data acquisition

MRI was performed with a Siemens Avanto 1.5 Tesla scanner with eight-channel Tim head-coil. Two MPRAGE sequences were collected for high-resolution anatomy: relaxation time (TR) = 1.91 s, echo time (TE) = 4.13 ms, inversion time (TI) = 1.1 s, flip angle = 15°, 120 slices, matrix 128 × 128, 1 × 1 × 1 mm resolution), which was designed to provide optimal differentiation of gray and white matter tissue for Freesurfer analysis. The following scans were also performed for neuroradiological evaluation: T2-weighted spin-echo (TR = 4.54 s, TE = 110 ms, 21 slices, 0.8 × 0.4 × 5.0 mm), FLAIR (TR = 9.0 s, TE = 89 ms, TI = 2500 ms, 21 slices, 1.0 × 0.8 × 5.0 mm), T2-weighted hemosiderin scan (TR = 800 ms, TE = 26 ms, flip angle = 30°, 21 slices, 1.1 × 0.9 × 5.0 mm) and one diffusion-weighted echo planar scan (TR = 5.0 s, TE = 84 ms, 23 slices, 1.7 × 1.7 × 5.0 mm, 30 directions, six B0 images).

Data analysis

Neuropsychological tests were scored per respective standards. The two rehabilitation outcome measures we considered were (1) the HVLT delayed correct recall score (range: 0–12), and (2) the RBMT total standardized profile score (range: 0–24). Injury severity was categorized as mild, moderate or severe based on a method used in our previous studies (Strangman et al., 2008b). In brief, 85% of cases were based on duration of loss of consciousness (LOC), (Ruff et al., 1993; Katz and Alexander, 1994; Dikmen et al., 2003) where mild was 0–30 min LOC (unless post traumatic amnesia exceeded 24 h); moderate was 30 min < LOC ≤ 24 h, and >24 h LOC was designated severe. In the absence of LOC data, the determination was made based on Glasgow coma scale scores (<9 = severe, 9–12 = moderate, 13–15 = mild) (Katz and Alexander, 1994; Temkin et al., 1995).

Structural MRI data was analyzed using Freesurfer v4.5 (Fischl et al., 2004), with the associated recon-all processing stream applied to each participant's pair of MPRAGE scans. In brief, the preprocessing stream consists of co-registering the two scans, non-uniform intensity normalization, Talairach transformation, skull stripping, volumetric labeling, tissue type segmentation, fitting of the two cortical surfaces (pial and white matter), and cortical parcellation. From this, data summaries for each cortical region are computed, including regional volume, surface area, curvature, and thickness; for details, including information on the 2009 version of the Destrieux atlas used for parcellation (see Fischl et al., 1999; Fischl et al., 2004; Destrieux et al., 2010; Schmansky, 2010). The Destrieux atlas was selected, as opposed to the Desikan atlas (Desikan et al., 2006), because it differentiates gyral and sulcal tissue. We deemed this difference potentially important in TBI where many gyri (as opposed to sulci) may contact the skull upon impact. In all, we collected – for each subject – gray matter volumes for 75 labeled cortical regions per hemisphere (see Figure 1) plus gray matter volume in 23 distinct subcortical regions. Of these latter 23 regions, we retained nine – specifically excluding total hemisphere measures, ventricles, vessels, image hypo-intensities, and the optic chiasm (Fischl et al., 2002). Upon completion of the automated processing, results were manually examined for clear failures to fit cortical surfaces (e.g., crossover of pial and white matter surfaces, errors related to inadequate skull stripping), major topological defects (holes or handles) in these surfaces, and failures in subcortical segmentation (e.g., distinct “islands” for contiguous regions). We sought to identify only gross errors for two reasons: (1) manual editing can cause bias, and (2) such editing is impractical in many clinical settings. Beyond the one failure to reach process completion (indicated earlier), no defects or failures were identified as requiring correction.

Figure 1

Brain morphometric parcellation results for an example subject, identifying key regions from this study. (A,B) Lateral and medial representations of the inflated cortical surface showing the 75 parcellated cortical regions per hemisphere. Only the left hemisphere is shown, laterally (A) and medially (B). (C) Coronal slice highlighting several subcortical structures automatically identified by the Freesurfer algorithms (thalamus, ventral diencephalon, hippocampus, basal ganglia, etc).

Regions of interest

As discussed, we hypothesized that the volume of the dorsolateral and ventrolateral prefrontal cortex (DLPFC, VLPFC), the PPC and the hippocampus were likely to provide significant predictive value for rehabilitation outcome. The Freesurfer region most coincident with the DLPFC region previously identified as involved in strategic memory processing (Strangman et al., 2008b) was the middle frontal gyrus (mFG). Our VLPFC region, also derived from (Strangman et al., 2008b), spanned two cortical labels, so we examined both the triangular and opercular parts of the inferior frontal gyrus. The predicted region of the PPC also straddled two cortical labels, so we examined both the angular and supramarginal portions of the inferior parietal gyrus. Each subject contributed two volume measurements per region, one from each hemisphere, as we hypothesized that the integrity of any particular region could be important bilaterally and in fact one hemisphere may “compensate” for damage in the opposite hemisphere (Cramer, 2004).

In addition to analysis of our a priori regions of interest, the remaining 70 cortical regions as well as eight of the subcortical brain regions (excluding ventricles, and other non-gray matter tissue types) were retained for further analysis, for a total of 84 regions analyzed.

Constructing the prediction model

Selecting a prediction model to use for all brain regions was achieved via a two-pass procedure. In the first pass, we began with a linear mixed effects regression model to predict posttest outcome from the following predictors: pretest outcome, age, gender, time since injury, education, moderateTBI, severeTBI, regional volume, TIV, and regional volume*TIV, where posttest outcome was either the HVLT delayed correct recall score or RBMT total standardized profile score at posttest 1 or 2 (i.e., immediately or 1 month after completion of therapy), pretest outcome was the same outcome measure at pretest, age was the participant's age in decades at the time of the MRI, gender was a binary variable, education and time since injury were in years, regional volume was either the left or right hemisphere volume of the brain region being tested, total intracranial volume (TIV) was the reciprocal of the volume of the total brain mask (gray matter plus white matter plus ventricles and cerebrospinal fluid, as generated by Freesurfer), and finally the interaction of regional volume with TIV. These last three terms were needed to compensate for the relationship between regional volumes and TIV, given that regional volume tends to be positively associated with TIV (Kronmal, 1993). ModerateTBI and severeTBI were binary variables indicating whether a participant had sustained a moderate or severe injury, as determined from Glasgow Coma Scale and length of loss of consciousness data (Strangman et al., 2008b). This model was fit for each region and outcome, and we then counted the number of models where each term survived a simple p < 0.05 significance threshold. Using this process, education, time since injury, and gender terms never reached significance in any of the models and were discarded from further consideration.

In the second pass, we computed the following multiple linear regression model for each brain region, which retained the remaining terms from our initial fits:

posttest_outcome = pretest_outcome + age + moderateTBI + severeTBI + regional_volume + TIV + regional_volume*TIV

The separate regional volume measurements are largely independent measures, but not entirely so due to part–whole relationships, so the optimal adjustment for multiple comparisons is not known. Due to the exploratory component of the study, we report all models where the final volume coefficient exhibited p < 0.05 uncorrected, discussing in detail those that that meet a conservative Bonferroni correction: for our primary hypotheses this is 0.5/6 = 0.0083, and for the whole-brain search this is p < 0.05/84 regions = 0.0006.

Supplementary analyses

Two additional analyses were also conducted. First, we sought to directly test for laterality effects to identify any potential hemispheric asymmetries. To do so, we utilized the same final regression model as described in the previous section, but added a binary variable to indicate whether the regional volume data came from the participant's left or right hemisphere.

Second, all candidate final models were checked for non-additivity, non-normality of residuals, quadratic components of regional volume and age, and effects of leverage. Leverage was identified by computing Cook's distance for every data point in the regression and determining whether the largest and fourth largest distance values were separated by a factor greater than 10 (Draper and Smith, 1998). All regression analyses were performed using R v2.8.1.

Results

Participant characterization and behavioral results

Demographics for study participants appear in Table 1. Injury evaluations appear in Table 2, revealing that many of the brain injuries were remote (>5 years), and nearly half were severe. The outcome measures pre- and post-therapy for TBI participants appear in Table 3, along with the results of simple paired t-tests of posttest versus pretest scores. As a whole, participants with TBI demonstrated significant improvements for both semantic memory (HVLT; up to 2.6 points) and everyday memory function (RBMT; 2.1 points), and at both posttest 1 and 2 (i.e., immediate and delayed time points; see also O'Neil-Pirozzi et al. (2010).

Table 2.

Injury severity, patients taking anticonvulsant medications, and cause of injury information for participants with TBI.

Condition	Participants
INJURY SEVERITY
Mild	12 (24%)
Moderate	12 (24%)
Severe	24 (48%)
Unverifiable	2 (4%)
Anticonvulsant medications	13 (26%)
CAUSE OF INJURY
Vehicle accident	38 (76%)
Blunt force trauma	8 (16%)
Fall	6 (12%)
Mixed	5 (10%)

Table 3.

Outcome scores for TBI participants (.

Time point	HVLT-R delayed recall		RBMT-II total score
	Mean (SD)	t (p)	Mean (SD)	t (p)
Pretest	6.2 (3.6)	–	16.3 (4.8)	–
Posttest 1 (immediate)	8.1 (3.2)	4.4 (<0.0001)	18.4 (4.8)	3.4 (0.001)
Posttest 2 (1-month follow-up)	8.8 (2.8)	7.4 (<0.0001)	18.4 (4.3)	3.2 (0.002)

Morphometric parcellation

The fully automated Freesurfer cortical reconstruction and parcellation process completed successfully for 50 participants. Even in cases of significant tissue degeneration (n = 7), the reconstruction process was able to identify all sulci and gyri, and provide reasonable parcellations and regional volumes. The results from a participant with a particularly large region of right temporal lobe encephalomalacia appear in Figure 2.

Figure 2

Example surface reconstruction in a case of substantial right temporal lobe encephalomalacia. (A) T1-weighted axial MRI image through the region of degeneration (oval). (B) Close up of T1 with overlay of the surfaces generated at the white-/gray-matter interface (green) and the pial surface (red). (C) Inflated right hemisphere showing curvature (reds = depth of sulci, greens = peak of gyri), with minimal curvature in the region of degeneration. (D) Re-folded surface with individual regional cortical labels (colors).

Prediction results

From our a priori regions of interest, we found a significant positive relationship between hippocampal volume (pooling across left and right hemispheres) and HVLT outcome at posttest 2 (p = 0.007; cf. Table 4 and Figure 3A). The regression coefficient for hippocampal volume was 1.056 per 1 cm3 (1,000 mm3) of tissue, whereas that for the HVLT pretest score was 0.53; thus a 1-unit change in pretest was roughly equivalent to a 0.5 cm3 change in hippocampal volume in terms of predicting HVLT at posttest 2. Over the full range of hippocampal volumes (1.335–5.182 cm3), the hippocampal volume term predicted a range of 4.1 points on the HVLT, or 34% of the full-scale (12 point) range. For comparison, the HVLT pretest score predicted a range of 6.38 points, or 53% of the full-scale HVLT range. The hippocampal prediction model also included a significant effect of age, whereby older individuals had better rehabilitation outcomes (1.6 points on HVLT across the span of four decades), and marginal effects of injury severity (severities greater than mild tended towards an average 0.85 point decrease in HVLT at outcome). Other results from our six regions of interest with volume terms exhibiting p < 0.05 appear in Table 5. This included the mFG which exhibited modest predictive value for HVLT at posttest 2 (coefficient = 0.27 on a 7.633 cm3 range) as well as at posttest 1, and the supramarginal gyrus, which predicted both HVLT and RBMT outcomes, but only at posttest 1.

Table 4.

All model terms for the primary prediction analysis, pooling across left and right hippocampal gray matter volume.

		HVLT outcome, Posttest 2
Region (volume range)	Variable	Coef.	CI-low	CI-high	t	p
Hippocampus (1.335–5.182)	Volume	1.06	0.29	1.82	2.74	0.007
	Prehvlt	0.54	0.42	0.65	9.26	0.000
	Age	0.39	0.03	0.76	2.13	0.035
	ModerateTBI	−0.91	−1.95	0.12	−1.76	0.082
	SevereTBI	−0.81	−1.75	0.13	−1.71	0.091
	TIV	4537	−2861	11936	1.22	0.226
	Volume*TIV	−12240	−25409	930	−1.85	0.068
	Intercept	9.27	8.56	9.97	26.14	0.000

Columns include regression coefficient, 95% confidence interval for the coefficient, T score and associated p-value.

Figure 3

Depiction of the relationship between HVLT delayed correct recall score at posttest 2 versus the volume of (A) the hippocampus, and (B) the posterior dorsal cingulate gyrus. Simple least-squares lines are plotted; detailed regression results appear in Tables 4 and 5.

Table 5.

Positive prediction regression results for the regional gray matter volume term in .

Region (volume range)		HVLT outcome					RBMT outcome
		Coef.	CI-low	CI-high	t	p	Coef.	CI-low	CI-high	t	p
Frontal: middle	post1	0.41	0.07	0.75	2.39	0.0188	0.49	−0.04	1.01	1.86	(0.07)
(4.137–11.770)	post2	0.27	0	0.53	2.01	0.0472	0.36	−0.14	0.85	1.43	n.s.
Parietal: supramarginal	post1	0.92	0.28	1.57	2.86	0.0053	1.25	0.25	2.25	2.48	0.015
(3.764–9.162)	post2	0.34	−0.16	0.84	1.35	n.s.	0.53	−0.43	1.48	1.1	n.s.
Hippocampus	post1	0.6	−0.44	1.65	1.15	n.s.	0.77	−0.78	2.31	0.98	n.s.
(1.335–5.182)	post2	1.06	0.29	1.82	2.74	0.0073	−0.18	−1.64	1.27	−0.25	n.s.

We next performed our exploratory analysis investigating the predictive value of all remaining brain regions generated by the Freesurfer parcellation (results in Table 6). The dorsal posterior cingulate cortex passed a strict Bonferroni correction (p < 0.0006) for predicting HVLT at posttest 2 (Figure 3B). The remaining entries in Table 6 are for those regions passing a p < 0.05 cutoff for at least one of the outcome measures. These identified regions, rather than being randomly distributed around the brain, fell into one of four geographic clusters. The first cluster was on the medial wall surrounding the corpus callosum, including four of five cingulate subregions and the subcallosal gyrus (immediately inferior to the anterior cingulate). The second cluster included the intraparietal sulcus and postcentral sulcus, which are immediately adjacent to our supramarginal and angular gyrus regions of interest. The third cluster included the occipital–temporal and medial lingual gyrus and sulcus, plus the inferior temporal sulcus, all on the inferior surface of the temporal lobe. The fourth cluster was in the diencephalon (thalamus and ventral diencephalon). A final region, the inferior frontal sulcus, was immediately adjacent to our two a priori VLPFC regions of interest. In essentially all regions, the regression coefficients were similar in magnitude to those discussed above, resulting in ranges of predicted outcome scores similar to those for outcome pretest scores. Only one region, the subcallosal gyrus, exhibited negative regression coefficients for RBMT – that is, increasing volume tended to predict poorer RBMT outcome scores.

Table 6.

Prediction results for the regional volume term in the brain-wide exploratory analysis.

Region (volume range)		HVLT outcome					RBMT outcome
		Coef.	CI-low	CI-high	t	p	Coef.	CI-low	CI-high	t	p
GYRUS PLUS SULCUS
Cingulate: anterior	post1	1.08	0.52	1.64	3.81	0.0003	0.67	−0.23	1.58	1.47	n.s.
(2.892–6.931)	post2	0.62	0.19	1.05	2.87	0.0051	0.27	−0.6	1.13	0.61	n.s.
Cingulate: mid-anterior	post1	1.53	0.54	2.53	3.06	0.0029	1.3	−0.31	2.91	1.6	n.s.
(1.523–3.887)	post2	0.73	−0.03	1.49	1.91	(0.06)	0.59	−0.93	2.11	0.77	n.s.
Cingulate: mid-posterior	post1	2.1	0.97	3.23	3.69	0.0004	1.86	0.08	3.63	2.07	0.041
(1.246–3.814)	post2	1.35	0.51	2.19	3.19	0.002	2.38	0.75	4.01	2.9	0.0047
GYRI
Cingulate: posterior dorsal	post1	3.02	1.27	4.77	3.43	0.0009	2.95	0.25	5.65	2.17	0.0329
(0.783–2.230)	post2	2.61	1.33	3.88	4.06	0.0001	1.86	−0.72	4.45	1.43	n.s.
Occ. temporal/medial lingual	post1	0.62	0.03	1.21	2.09	0.0399	0.32	−0.59	1.23	0.7	n.s.
(1.664–7.118)	post2	0.09	−0.37	0.54	0.38	n.s.	0.29	−0.56	1.14	0.68	n.s.
Subcallosal	post1	1.62	−0.96	4.21	1.25	n.s.	−4.75	−8.4	−1.09	−2.58	0.0116
(0.056–1.052)	post2	0.39	−1.57	2.36	0.4	n.s.	−3.24	−6.72	0.25	−1.84	(0.07)
SULCI
Frontal: inferior	post1	0.88	−0.17	1.93	1.67	(0.10)	1.61	0.05	3.17	2.05	0.0432
(1.508–4.133)	post2	0.51	−0.29	1.32	1.26	n.s.	−0.14	−1.64	1.36	−0.19	n.s.
Intraparietal/parietal transv.	post1	0.97	0.06	1.88	2.12	0.0371	1.21	−0.17	2.59	1.74	(0.08)
(2.188–5.254)	post2	1	0.34	1.65	3.01	0.0034	1.11	−0.18	2.41	1.71	(0.09)
Occ. temporal/medial lingual	post1	0.61	−0.24	1.46	1.42	n.s.	1.3	0.04	2.56	2.06	0.0427
(1.245–4.827)	post2	0.62	−0.02	1.26	1.94	(0.06)	0.72	−0.47	1.91	1.21	n.s.
Postcentral	post1	0.42	−0.19	1.03	1.38	n.s.	0.93	0.02	1.84	2.03	0.0456
(1.711–5.543)	post2	0.53	0.06	0.99	2.26	0.0263	0.96	0.11	1.81	2.24	0.0272
Temporal: inferior	post1	0.81	−0.15	1.76	1.68	(0.10)	1.63	0.23	3.02	2.31	0.0229
(0.490–3.535)	post2	0.7	−0.03	1.42	1.91	(0.06)	1.21	−0.11	2.53	1.82	(0.07)
SUBCORTICAL
Thalamus	post1	0.53	−0.06	1.13	1.78	(0.08)	1.03	0.18	1.89	2.4	0.0186
(4.267–9.933)	post2	0.68	0.26	1.11	3.18	0.002	0.13	−0.7	0.96	0.32	n.s.
Ventral diencephalon	post1	0.49	−0.67	1.65	0.84	n.s.	2.58	1.02	4.14	3.3	0.0014
(2.266–5.236)	post2	1.23	0.4	2.06	2.94	0.0041	1.02	−0.51	2.55	1.33	n.s.

Notes: p-values in brackets indicate trends (p < 0.1). Designations of post1 and post2 refer to posttest 1 (immediately following therapy) and posttest 2 (1-month follow-up).

The strongest findings in Table 6 appeared in the cingulate gyrus, including the anterior cingulate, anterior mid-cingulate, and dorsal posterior cingulate for HVLT, and the posterior mid-cingulate cortex (pMCC) for both HVLT and RBMT (see Table 6). The pMCC produced the strongest effect of all regions at posttest 2 for RBMT, the most functionally relevant outcome (Figure 4). With a regression coefficient estimate of 2.378 and a volume range of 2.568 cm3, this region had a prediction range of 6.1 points on the RBMT 24 point scale, compared to a prediction range of 6.7 points over the full observed range for RBMT pretest scores. As with all other models, the positive regression coefficient implies that larger volumes were associated with better outcome scores.

Figure 4

Relationship between our functionally relevant outcome (RBMT) 1-month post-therapy and the volume of the pMCC. A simple least-squares regression line is plotted. Results from the multiple regression model appears in Table 6.

Our test for laterality differences only revealed one significant effect, in the anterior cingulate region. In particular, for the HVLT outcome at posttest 1, the left anterior cingulate predicted an outcome score 1.3 points higher than the right (t = 2.73, p = 0.02).

Finally, our Cook's distance check for outliers proved negative in all cases where significant relationships with volume were identified, indicating that no individual data points had substantial influence on the regression results. Each regression with a significant regional volume term was also checked for quadratic effects of age and regional volume. Two volume-squared terms were found significant: there were negative quadratic effects of volume in the thalamus and ventral diencephalon when predicting HVLT at posttest 2.

Discussion

Using a fully automated morphometric brain parcellation technique, we measured regional gray matter volumes in individuals with chronic TBI prior to their participation in a standardized memory rehabilitation program. We then identified brain regions for which the regional brain volume predicted post-therapy outcome. Substantial predictive value was found both for an outcome closely related to the rehabilitation training (i.e., HVLT), as well as for a more functionally relevant outcome (i.e., RBMT). Predictive regions included the hippocampus, the DLPFC and the PPC, which had been previously demonstrated as important to task performance and exhibited decreased activity in individuals with TBI versus controls (Strangman et al., 2009). In exploratory analysis, we also observed substantial predictive value from various subregions of the cingulate gyrus, areas adjacent to the predicted VLPFC and PPC regions, and the thalamus. Importantly, the magnitude of the regression coefficients for regional volumes rivaled those for pretest scores on the outcome variable – that is, knowledge of regional brain volumes tended to account for similar amounts of variance on the posttest score as did the pretest score. Finally, we found very little evidence for hemispheric specialization. Only in the anterior cingulate did we find modest support for the left hemisphere predicting significantly higher rehabilitation outcome scores than the right hemisphere.

Predicting outcome from memory rehabilitation

Given the hippocampus’ well known role in learning and memory (Squire, 1992; Burgess et al., 2002), it is not surprising that the volume of the hippocampus exhibited a positive relationship with memory improvement following memory rehabilitation. This is also consistent with previous findings of degeneration in the medial temporal lobe following TBI (Ng et al., 2008), as well as the association between hippocampal damage and TBI outcome (Umile et al., 2002).

The findings in the mFG and supramarginal gyrus, as well as the regions adjacent to the supramarginal and angular gyri and in the VLPFC, are all regions previously identified as functionally important in HVLT-like tasks (Strangman et al., 2008b, 2009). This adds support to the notion that the volume of regions functionally relevant to task performance is important to rehabilitation outcomes. The same argument also extends to the cingulate regions identified in this study. While in previous work the cingulate did not exhibit a main effect difference in activation between controls and individuals with TBI, functional activation in this region held a key role in various interactions between population and task performance (Strangman et al., 2009).

While the location of identified regions was consistent with previous work, the large magnitude of the regression coefficients for regional volume was less expected. Prediction ranges afforded by regional volume often rivaled those of pretest scores, even after accounting for variables such as pretest, age, and injury severity. This strongly supports the notion that regional volumes can provide substantial added value for predicting rehabilitation outcome, at least in individuals one or more years post-injury. The generalizability of these findings to individuals less than a year post-injury remains to be investigated. However, we do note that structural MRI scans are commonly obtained on patients with TBI, such scans are the only data required for the Freesurfer parcellation programs, and our results involved no manual editing of the parcellation process, all of which support the clinical feasibility of investigating individuals closer to their date of injury.

The reason for the extent of the added value remains to be understood. The regional specificity and lack of significant predictive value for global measures such as all cerebral gray matter, argues strongly against a mass effect. It is likely, however, that the anatomical specificity provided by regional morphometric measurements, and/or the sensitivity of such measurements as compared to radiological reports (Bigler, 2001b), provides key information about the disruptions to functional brain networks. Such information is presumably lost when simply considering, for example, overall injury severity or total lesion load. Other imaging modalities that provide detailed anatomical information about individual injuries, such as diffusion tensor imaging (Wilde et al., 2006), may prove equally or more useful in predicting outcomes from cognitive rehabilitation programs, or may complement morphological measures such as those used here. Regardless of imaging modality, we speculate that optimizations can still be made in terms of the level of anatomical detail considered. Our morphological measures lie in a middle ground between gross pathology and voxel-wise assessments. It is not yet clear which level of analysis may be optimally effective at predicting TBI rehabilitation outcomes.

One particularly interesting finding from the exploratory analysis was in the pMCC, which exhibited substantial predictive capability for the more functionally relevant RBMT outcome at 1-month post-therapy. The RBMT is comprised of a number of quite disparate memory tasks, each of which has considerable ecological validity (e.g., remembering where a personal item is, remembering to ask an appointment-related question, and remembering to deliver a message). While the prefrontal cortex in general is known to be involved in working memory (Hillary et al., 2006), long-term episodic memory (Ranganath et al., 2003), and executive control (Messinger et al., 2009; Rossi et al., 2009), it is still unexpected that any one brain area might be particularly important for successful performance across the broad range of RBMT tasks. While predicting RBMT at 1-month post-therapy is of high clinical importance, we nevertheless anticipated only modest success in predicting functionally relevant outcomes such as RBMT, particularly a month post-intervention. Identification the pMCC as having substantial predictive range for this outcome was thus quite promising. It provides evidence in favor of predicting outcome not just for laboratory tasks like the HVLT, but also for more ecologically valid outcomes. It further provides guidance on what brain regions and networks to focus on.

Unlike more anterior portions of the cingulate cortex, relatively little is known about the pMCC and its involvement in cognitive processing. However, the pMCC is known to interact with the PPC (Vogt, 2005), a key area for multimodal integration. It was also recently found to be functionally linked with various medial temporal lobe regions as well as the insular cortex (Taylor et al., 2009). One hypothesized role is orienting the individual to sensory events (Vogt, 2005). Collectively, the evidence suggests that pMCC interacts with areas involved in multimodal integration and declarative and autobiographical memory processing. It is therefore plausible that the pMCC is a key component of a functional network involving stimulus orientation, integration, and transfer to memory. If sufficiently general purpose, such a network would be expected to have predictive value for specific outcomes (e.g., HVLT) as well as more general outcomes (e.g., RBMT), as we observed. It remains to be determined whether this posited integrative role of the pMCC is specific to memory rehabilitation or extends in some way to the rehabilitation of attentional or executive functioning following TBI.

Finally, previous TBI studies have demonstrated negative relationships between age and rehabilitation outcome. In addition, a number of researchers have suggested that TBI may have a particularly negative effect on older brains (Rapoport and Feinstein, 2000; Bigler, 2001a; Marquez de la Plata et al., 2008), implying that age at time of injury may be an important variable in rehabilitation outcome. In this study, we found age to be positively related to rehabilitation outcome when predicting HVLT from the volume of the hippocampus. This positive relationship with age held when controlling for age at time of injury, and continued to hold in supplementary analyses where we excluded all volume terms, and even when predicting HVLT outcome solely from age. Thus, we had no evidence of a negative relationship between age and outcome. While these findings contradict previous results, we note that (1) our therapy was atypical relative to other studies in that it focused exclusively on internal strategies for memory improvement, and (2) the coefficient on age was quite modest at approximately 0.4 points on HVLT per decade, meaning our oldest participant was predicted to score 1.6 points higher than our youngest, across our 40 year participant age range. Thus, while significant and positive, the effect of age was substantially less than half the effect contributed by volume, and in a range below what would typically be considered a clinically meaningful change.

Limitations

Our study included individuals with a wide range of injury severities. Thus, there may have been differential responses from mild or severe cases that did not correspond to the “average” prediction coefficients we observed, or the group as a whole may have been skewed by one particular subpopulation. However, we consider our severity range to be a strength of the study: among individuals post-TBI who present with memory difficulties, the ability to predict memory rehabilitation outcome remained strong in spite of such heterogeneity, and hence the findings may apply to a relatively broad clinical population. Similarly, although the MRI scans were all chronic and were collected no sooner than 1-year post-injury, the time since injury and rehabilitation histories differed widely across participants. This additional variability almost certainly contributed to reduced (rather than inflated) prediction sensitivity. In addition, our approach may be less sensitive if applied exclusively in mild TBI where neurostructural changes may be less evident, or exclusively in severe TBI where MRI may be contraindicated. It remains to be determined whether our approach would be suitable for individuals who sustained a TBI less than 1 year prior, particularly given the variability of injury evolution in the first few months immediately post-injury.

We recognize that automated approaches may contain segmentation or parcellation errors, may be fooled by hemorrhagic lesions, and may miss subtle variations or complexities that can be more appropriately handled by an expert neuroanatomist. For a relatively large study such as this, we believe the benefits of automation actually make not only the study, but potential clinical translation of the approach, realistically feasible. Related, it is difficult to judge volumetric accuracy in face of encephalomalacia or other gross structural injury. Scar tissue may lead to over-estimation of gray matter volume, white matter may degenerate sufficiently to be classified as gray matter, or a major loss of tissue may result in unrecoverable segmentation failures. For this study, we assumed that such effects contribute to the noise in our data and prediction models. Finally, there are various alternatives to compensate for TIV, a common approach being to divide all regional volumes in an individual by their TIV. While the best approach is debated, the divisor approach is has certain statistical shortcomings (Kronmal, 1993), and hence we utilized an interaction approach to help control for differing brain volumes.

Conclusions

We identified a number of specific regions that provided substantial predictive value for memory rehabilitation outcomes. While this represents only an early stage in predicting outcomes from specific cognitive rehabilitation programs, and is specific to just one such program, our work suggests that the approach is viable, and provides a benchmark with which to compare other neuroimaging prediction approaches. Since our study was conducted at the chronic stage (all participants being more than 1-year post-injury), further work needs to be done to determine if more acute imaging provides equivalent, better, or worse predictive value. If the prediction of cognitive rehabilitation outcomes can be similarly achieved with other interventions, structural neuroimaging could eventually become a tool to help individuals with TBI and their caregivers more objectively select among candidate rehabilitation programs based on post injury structural scans. Based on our prior work, we further anticipate that combining results from multiple imaging modalities – for example, T1-, T2-, hemosiderin-, and diffusion-weighted MRI, spectroscopy, resting functional MRI and task-related functional MRI – may help achieve even greater sensitivity, specificity and predictive value. Such an approach must always balance the costs and benefits of scanning and analysis relative to the therapy. However, for expensive and human-intensive efforts such as TBI rehabilitation, neuroimaging may provide a reliable method for better targeting and customizing TBI rehabilitation therapies.

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.