Prospective longitudinal MRI study of brain volumes and diffusion changes during the first year after moderate to severe traumatic brain injury.

The objectives of this prospective study in 62 moderate-severe TBI patients were to investigate volume change in cortical gray matter (GM), hippocampus, lenticular nucleus, lobar white matter (WM), brainstem and ventricles using a within subject design and repeated MRI in the early phase (1-26 days) and 3 and 12 months postinjury and to assess changes in GM apparent diffusion coefficient (ADC) in normal appearing tissue in the cortex, hippocampus and brainstem. The impact of Glasgow Coma Scale (GCS) score at admission, duration of post-traumatic amnesia (PTA), and diffusion axonal injury (DAI) grade on brain volumes and ADC values over time was assessed. Lastly, we determined if MRI-derived brain volumes from the 3-month scans provided additional, significant predictive value to 12-month outcome classified with the Glasgow Outcome Scale-Extended after adjusting for GCS, PTA and age. Cortical GM loss was rapid, largely finished by 3 months, but the volume reduction was unrelated to GCS score, PTA, or presence of DAI. However, cortical GM volume at 3 months was a significant independent predictor of 12-month outcome. Volume loss in the hippocampus and lenticular nucleus was protracted and statistically significant first at 12 months. Slopes of volume reduction over time for the cortical and subcortical GGM were significantly different. Hippocampal volume loss was most pronounced and rapid in individuals with PTA > 2 weeks. The 3-month volumes of the hippocampus and lentiform nucleus were the best independent predictors of 12-month outcome after adjusting for GCS, PTA and age. In the brainstem, volume loss was significant at both 3 and 12 months. Brainstem volume reduction was associated with lower GCS score and the presence of DAI. Lobar WM volume was significantly decreased first after 12 months. Surprisingly DAI grade had no impact on lobar WM volume. Ventricular dilation developed predominantly during the first 3 months, and was strongly associated with volume changes in the brainstem and cortical GM, but not lobar WM volume. Higher ADC values were detected in the cortex in individuals with severe TBI, DAI and PTA > 2 weeks, from 3 months. There were no associations between ADC values and brain volumes, and ADC values did not predict outcome.

Introduction

General cerebral atrophy is a common consequence of moderate to severe traumatic brain injury (TBI) (Trivedi et al., 2007; Ding et al., 2008; Sidaros et al., 2009; Warner et al., 2010a). This general atrophy is not predicted by focal lesion volume (Marcoux et al., 2008), and develops over at least 3 years with the bulk loss occurring in the first year (Sidaros et al., 2009; Blatter et al., 1997). However, there is limited knowledge on the trajectories of volume changes for different brain structures, especially in humans. There is evidence from rodent models of closed head injury that brain region and injury type determine the histopathological response, and hence degree of tissue loss and outcome (Fox et al., 1998; Immonen et al., 2009; Meaney and Smith, 2011). Longitudinal studies of structural changes following TBI in humans have so far only contained two, often poorly defined, time points of MRI (Ross, 2011), which is insufficient to describe trajectories of volume loss in different brain structures. Furthermore, the impact of different clinical characteristics related to TBI subtype classification, e.g. injury severity and presence of diffuse axonal injury (DAI)/traumatic axonal injury, on trajectory of volume changes in different brain structures has not been explored systematically.

Based on the animal literature we predicted that cortical volume declines first, mainly between the early and 3-month scans. To the best of our knowledge there are no human studies documenting longitudinal cortical volume changes postinjury. Rodent TBI studies consistently show that the histopathological response including neuronal necrosis and apoptosis is largely completed within 3 months after injury in the cortex (Conti et al., 1998; Rodriguez-Paez et al., 2005; Sato et al., 2001), although volume loss may continue for 12 months (Smith et al., 1997).

In the hippocampus, on the other hand, protracted neuronal loss lasting for at least 12 months is described in animals (Conti et al., 1998; Smith et al., 1997; Williams et al., 2001). In addition, prominent changes in synapse morphology and reduced hippocampal neurogenesis are reported (Gao et al., 2011; Scheff et al., 2005; Gao et al., 2008). Hippocampal volume loss is reported to be greatest during the first weeks, but continues for at least 12 months in rodents (Immonen et al., 2009; Smith et al., 1997). Thus, we predicted that hippocampal volume loss in humans following moderate–severe TBI would be protracted, and continue during the entire observation period.

For the lenticular nucleus less information is available, but quite disparate histopathological results are reported for other subcortical GM nuclei. In animals and humans, both very early and slowly evolving histopathological changes have been shown (Hicks et al., 1996; Greer et al., 2011; Tong et al., 2002). Patient studies have documented both volume loss (Warner et al., 2010b; Bendlin et al., 2008; Wilde et al., 2007) and no volume change in the lenticular nucleus (Salmond et al., 2005). Taken together the data suggest that lenticular nucleus volume most likely declines slowly, and we expected to see a significant reduction at 12 months.

Lobar WM volume changes are described as the slowest and take place over a very protracted time period due to Wallerian/Wallerian-like degradation being lengthy and incomplete in the brain compared to the peripheral nervous system (Vargas and Barres, 2007; Coleman, 2005; Coleman and Freeman, 2010). Significant lobar WM loss was therefore expected to develop between 3 and 12 months after injury.

Human studies have shown that the brainstem volume is reduced in the chronic phase after TBI (Sidaros et al., 2009; Kim et al., 2008), but animal studies have not focused on this region. We anticipated that brainstem volumes decreased slowly, similar to the lobar WM.

TBI leads to ventricular dilation in both humans and animals. One patient study reports that the period of ventricular dilation lasts until 2–7 months postinjury with no significant increase thereafter (Blatter et al., 1997). In rodents ventricular enlargement is described to take place until 6 months postinjury (Bouilleret et al., 2009; Liu et al., 2010). We therefore expected the ventricular volume increase to take place mainly between the early and 3-month scans. Ventricular dilation is considered to result primarily from WM volume loss (Anderson and Bigler, 1995), but if ventricular dilation develops early after TBI, GM volume loss may also contribute significantly. We therefore investigated the relationship between ventricular volume change over time and the volume changes in cortical GM, lobar WM and brainstem to establish the relative importance of volume loss in these structures to ventricular dilation.

Furthermore, the impact of clinical measures of TBI mechanism and severity, i.e. Glasgow Coma Scale (GCS), duration of post-traumatic amnesia (PTA), and DAI, on volume loss in different brain structures was assessed. Lower GCS scores were expected to lead to more severe loss of lobar WM and brainstem volumes, as unconsciousness postinjury is associated with deeper lesions (Jenkins et al., 1986) and increasing mechanical forces converging in the midbrain (Meaney and Smith, 2011). Longer duration of PTA (>2 weeks) was predicted to be associated with lower hippocampal volumes since PTA is a manifestation of hippocampal dysfunction and hence likely degree of injury. DAI was predicted to result in increased loss of lobar WM volumes, and increasing DAI severity was predicted to lead to greater brainstem volume loss. Based on postmortem evidence, DAI also gives rise to increased cortical thinning and could thus reduce cortical volume (Maxwell et al., 2010).

GM volume changes could be accompanied by changes in GM microstructure. The apparent diffusion coefficient (ADC) obtained from diffusion weighted imaging (DWI) is a measure of the degree of random isotropic diffusion. In GM diffusion is not restricted to a certain direction as in WM, and ADC values are considered to reflect cell number (neurons and for instance inflammatory cells), cell size (swelling vs. shrinkage), as well as other cytoarchitectonic features such as the composition and fraction of neuropil and extracellular matrix. It has previously been demonstrated that TBI patients with lesions visible on conventional MRI have increased ADC values in the brainstem the first month after TBI (Brandstack et al., 2011), while ADC values in cortex and subcortical GM increased between ~2 and ~3 months in moderate–severe TBI (Bendlin et al., 2008). Thus we expected the severe TBI patients to have increased ADC values in the normal appearing midbrain at 3 months. On the other hand, the ADC values in normal appearing cortical GM should increase more between 3 and 12 months in both moderate and severe TBI. Since DAI leads to loss of neurons, neuronal soma shrinkage, axonal changes and myelin loss in cortex (Greer et al., 2011; Coleman and Freeman, 2010; Maxwell et al., 2010), we expected that cortical ADC values would be higher in DAI than in non-DAI patients after 3 and 12 months. Moreover, we predicted that longer PTA would be associated with increased hippocampal ADC values due to increased cell loss.

The primary objectives of this prospective study in moderate and severe TBI patients were to investigate the trajectories of volume change in total cortical GM, hippocampus, lenticular nucleus, lobar WM, brainstem and ventricles using a within subject design and repeated MRI in the early phase and 3 and 12 months after injury and to relate the volume changes to clinical measures of TBI mechanism and severity; GCS score, PTA, and DAI. Changes in local diffusion properties in GM were studied to identify any associations between changes in tissue microstructure as described with ADC values and the clinical measures of TBI  mechanism and severity; GCS score, PTA, and DAI. In addition the relationship between ADC values and brain volumes was examined. Our final aim was to determine if MRI-derived brain volumes and ADC values obtained after 3 months provided additional, significant predictive value to 12-month outcome after adjusting for GCS score, PTA and age.

Materials and methods

The study was approved by the Regional Committee for Medical Research Ethics and the Norwegian Social Science Data Services. Written consent was obtained from the patient or, for individuals underaged or incapacitated, their next of kin. Permission was obtained from the Norwegian Directorate of Health to use data from the deceased without consent from their next of kin.

Subjects

124 moderate to severe (GCS score ≤ 13) TBI patients aged 14–65 years were admitted to the Neurosurgical Department, St. Olav's Hospital, Trondheim University Hospital, Norway, in the period October 10, 2004–October 10, 2007. St. Olav's Hospital is the only level 1 trauma center in a region of 680,000 inhabitants. 17 of these patients did not participate in MRI study (did not consent, could not be examined with MRI or did not participate for other reasons). An additional 20 patients died before acquisition of follow-up MRI scans. Longitudinal MRI data was thus available for 87 TBI patients. Only patients with at least 2 successful segmentations of the T1-weighted 3D brain volumes were included in the study, therefore another 25 patients were excluded.

For each T1-weighted 3D scan the segmentation quality was assessed visually and scans were excluded if hematomas were included as brain tissue, or GM, WM or ventricles were incorrectly delineated.

See Fig. 1 for details about patient inclusion and exclusion and Fig. 2 for an example of successful NeuroQuant segmentation. Six scans were excluded due to insufficient segmentation quality. Five of these were acquired in the early phase, and segmentation failed due to hematomas segmented as brain tissue (n = 3), incorrect delineation of WM and GM of injured frontal lobe (n = 1), and excessive motion (n = 1). One 3-month scan failed due to an MRI technical mistake.

Fig. 1

TBI patients with GCS ≤ 13 that could not be explained by other factors than traumatic head injury, admitted to the Neurosurgical Department, St. Olav's Hospital, Trondheim University Hospital, Norway in the period 2004–2007. Flow chart of the inclusion and exclusion criteria, and the final number of patients participating in the current study.

Fig. 2

Successful NeuroQuant segmentation of repeated MRIs obtained in a 21-year-old patient with moderate traumatic brain injury at a) 9 days postinjury and b) 3 months and c) 12 months postinjury.

Unsuccessful segmentations resulted in exclusion of two TBI patients. Furthermore, in four TBI patients the number of scans was reduced from three to two.

In total 161 segmented T1-weighted 3D scans from 62 moderate to severe TBI patients were included in the final sample. 37 patients had three and 25 patients had two scans successfully segmented.

MRI and image processing

Scanning was performed at three different time points: early phase (0–26 days postinjury) and 3 months and 12 months postinjury. The early phase scan was acquired as soon as the clinical condition of the patients allowed it.

Scan protocol

MRI acquisition was performed on a 1.5 T Siemens Symphony Sonata scanner (Siemens, Erlangen, Germany) with an eight-channel head coil. Imaging was performed following a standardized protocol, where slice positioning for transversal scans was parallel with the anterior–posterior commissural line. A T1-weighted 3D magnetization prepared rapid acquisition gradient echo (MPRAGE) sequence was obtained in the sagittal plane with TR = 7.1 ms; TE = 3.45 ms; flip angle = 7°; TI = 1000 ms; FOV = 256 × 256; acquisition matrix of 256 × 192 × 128, reconstructed to 256 × 256 × 128, giving a reconstructed voxel resolution of 1.00 × 1.00 × 1.33 mm. For diffusion-weighted imaging a single-shot spin echo planar imaging sequence with 19 slices, slice thickness 5 mm, TR = 3300 ms; TE = 110 ms; NEX = 4; acquisition matrix size 256 × 256; FoV = 230; in-plane resolution 0.9 × 0.9 mm and baseline images (b = 0 s/mm2) plus varying diffusion gradient strength along each of three orthogonal directions with b = 500 and 1000 s/mm2 was used. Diffusion trace maps were computed from the isotropic diffusion images and used to estimate ADC values. FLAIR, T2, and T2*-weighted hemosequences were also acquired.

Segmentation of the T1-weighted 3D brain scans

Fully automated segmentation was performed using NeuroQuant (CorTechs Labs, La Jolla, CA, US) which is a FDA-cleared tool for clinical evaluation of hippocampal atrophy in mild cognitive impairment and Alzheimer's disease (Brewer et al., 2012). The output from NeuroQuant segmentation contains volumes of the total lobar WM, total cortical GM, lateral, third, fourth, and inferior lateral ventricles, cerebellum, hippocampus, amygdala, caudate, putamen, pallidum, thalamus and brainstem given as both absolute volume in cm3 and relative to the intracranial volume (ICV). In this study only volumes relative to ICV (%ICV) are used. The ICV-corrected volumes of the segmented structures from the left and right hemispheres were combined. The ICVs in the early phase and 12-month scan were significantly different (linear mixed model; mean ± standard error: early phase 1740.06 ± 19.98 cm3; at 12 months 1706.10 ± 19.74 cm3; p < 0.001). ICV is estimated from the total brain volume in NeuroQuant, and may therefore be influenced by brain edema in the early phase and brain atrophy in the chronic phase. Differences between ICV measurements based on total brain volume in the early and chronic phases are not surprising in TBI patients (Ding et al., 2008; Sidaros et al., 2009). Therefore, volumes of segmented structures were normalized and expressed as proportions of the ICV estimated from the 3-month scan.

Because of the large number of segmented brain structures, we chose to analyze the brain structures with low frequency of segmentation mistakes. The ICV-corrected volumes of cortical GM, hippocampus, brainstem, lenticular nucleus, lobar WM, and ventricles were compared between the early, 3-month and 12-month scans. Putamen and pallidum were evaluated together as they form the lenticular nucleus. The lenticular nucleus was used since it was better delineated than the thalamus in the TBI patients. Total ventricular volume was calculated as the sum of third, fourth and lateral ventricular volumes. We evaluated the ventricular volumes together since all ventricles have been shown to dilate in TBI (Warner et al., 2010b; Kim et al., 2008).

In addition the volume change over time was calculated for all structures. For the solid structures the following was used: early phase volume minus volume at 3 months, and early phase volume minus volume at 12 months. For the ventricles: volume at 3 months minus early phase volume, and volume at 12 months minus early phase volume. This was done in order for all mean changes to be positive.

Total brain atrophy from the early phase to 3 and 12 months was estimated using linear mixed model based on the sum of the volume changes of the solid structures (cortical GM, hippocampus, lenticular nucleus, lobar WM, and brainstem) subtracting the ventricular volume change (as ventricular volume change is defined as reversed, please see above).

ADC measurements in GM

ADC measurements were performed in PACS using Spectra Workstation IDS5 11.4.1 in all 62 patients. Ten circle-shaped ROIs with radius of 2.7 mm were positioned in apparently healthy GM bilaterally in the cortex, hippocampus and the brainstem (Fig. 3). These ROIs were chosen because they are located to regions where ADC value changes have been reported earlier (see Introduction), the regions are readily identifiable, and the ROIs are located to parts of GM with minimal WM contamination. All ROIs were positioned in apparently healthy tissue within the structure of interest. In patients with visible focal pathology on MRI the ROI was moved to the closest healthy tissue if necessary. The cortical ROIs in the superior frontal sulcus and postcentral sulcus were positioned at the sulcal depth in order to avoid measuring ADC in GM close to the skull with possible minor contusions. The cortical ROIs were placed on the same slice, just above the most dorsal slice displaying the lateral ventricles. In the posterior insular gyrus ROIs were set on the slice where the head of the caudate nucleus had its greatest extent. ROIs in the hippocampus and brainstem were positioned on the most caudal slice where mesencephalon was still visible. In the hippocampus the ROIs were located as laterally as possible in the anterior part of the hippocampal body in order to avoid placing the ROIs within the hippocampal head which frequently suffers from contusions in TBI. For all ROIs, the ADC values from the left and right hemispheres were averaged. The image analyst (V.B.) was blinded to clinical information. Intra-rater reliability with regard to ROI positioning and subsequent ADC measurements was assessed in 14 scans which were reanalyzed after >1 month from the first assessment. Intraclass correlation coefficients (ICC) using two-way random single measures were used to calculate reliability.

Fig. 3

ROI placement for ADC measurements in a) superior frontal sulcus ROI, postcentral sulcus ROI; b) insular gyrus ROI; c) hippocampal ROI, brainstem ROI. Circle-shaped ROIs with radius of 2.7 mm were positioned in apparently healthy tissue. In TBI patients with visible focal pathology, the ROI was moved to the closest apparently healthy GM tissue within the structure of interest.

Injury-related variables

GCS score was recorded at hospital admittance, or before intubation in case of a pre-hospital intubation. All patients included had GCS ≤ 13 that could not be explained by other factors than head injury, as explained in previous publications from this cohort (Skandsen et al., 2010; Moen et al., 2012). The head injury severity scale (HISS) is based on GCS, where GCS ≤ 8 is considered severe while GCS 9–13 as moderate head injury (Stein and Spettell, 1995).

GCS scores were available for all included patients.

PTA was recorded and dichotomized into PTA > 2 weeks and PTA ≤ 2 weeks which corresponds to moderate and moderate–severe Mississippi PTA intervals, respectively. Mississippi PTA intervals have been shown to be more accurate predictors of TBI outcome than the traditional Russell intervals (Nakase-Richardson et al., 2011). Data on PTA were missing in 4 patients. These were excluded from analyses including PTA.

DAI classification was performed by experienced senior neuroradiologists based on the early FLAIR, diffusion-, T2- and T2*-weighted scans, for details see Skandsen et al. (2010). DAI was classified into grade 1; traumatic lesions confined to lobar WM, grade 2; lesions also detected in the corpus callosum, and in grade 3; presence of brainstem lesions (Gentry, 1994). DAI classification was available for all included patients.

Outcome assessment

Global outcome was assessed by telephone interview 12 months after injury using the structured interview for Glasgow Outcome Scale—Extended (GOSE) (Jennett and Bond, 1975; Jennett et al., 1981). To reduce the potential error associated with the telephone setting (LeGrand et al., 2007), relatives or caregivers also provided information, in two-thirds of the cases, and the best source of information was used (based on the judgment of the authors) for details see Skandsen et al. (2010). Outcome data were available for all included individuals, but missing in 5 non-included patients. Patients were dichotomized into disability (GOSE < 7) and good recovery (GOSE ≥ 7) after 12 months.

Statistical analyses

Characteristics of included and excluded TBI patients.

The Mann–Whitney U-test, χ2-test, independent samples T-test and the Wilcoxon U-test were used to compare demographic and clinical data.

Longitudinal brain volume changes

Significant changes in the volume of each of the brain structures at each time point was calculated with a linear mixed effects model.

The relationships between ventricular volume change and volume changes of lobar WM, cortical GM and brainstem were explored using linear regression separately for the periods early to 3 months and early to 12 months. There were separate regression models for lobar WM, cortical GM, and brainstem volume changes as the volume changes were expected to be correlated.

Impact of TBI subgroup on brain volumes was investigated in the following dichotomized groups; moderate vs. severe TBI, PTA ≤ 2 weeks vs. > 2 weeks, and DAI vs. non-DAI. The volumes were compared over time within each group and between the dichotomized groups. For DAI a subdivision into four DAI groups, i.e. non-DAI = 0, and DAI grades 1, 2, and 3 were performed and the effect of degree of DAI was assessed with one-way ANOVA.

A significant F-value was followed by a paired Student's T-tests within the group, and with independent samples T-test between the dichotomized groups. The resulting p-values were Bonferroni corrected for number of time points evaluated for each structure separately.

The following analyses of volume differences between the different injury mechanism groups were performed:

Moderate vs. severe TBI

Differences between cortical GM, hippocampal, lenticular nucleus, lobar WM, brainstem and ventricular volumes in the early phase and at 3 and 12 months were investigated. Since DAI was more frequent in severe TBI, the presence of DAI was included in the analysis of lobar WM volume differences between moderate and severe TBI.

Differences between cortical GM, hippocampal, and ventricular volumes in the early phase and at 3 and 12 months were investigated.

DAI vs. non-DAI

Differences in cortical GM, lobar WM, brainstem, and ventricular volumes in the early phase and at 3 and 12 months were investigated.

DAI grades 0–3

The effect of degree of DAI (no DAI, DAI grade 1, 2 or 3) as an ordinal variable on brainstem and total ventricular volumes was assessed in the early phase and at 3 and 12 months.

Differences in volume changes over time between brain structures

The presence of statistically significant differences in rate of volume loss between different brain structures was evaluated in a linear mixed model with random intercept since this model fitted the data best, i.e. had the lowest Akaike's information criterion. There were not enough observations to build a model exploring non-linear changes over time. The following comparisons were done: cortical GM vs. hippocampus, cortical GM vs. lobar WM, hippocampus vs. lentiform nucleus, and hippocampus vs. lobar WM matter. For these analyses significant differences in the interaction term between volume loss over time and brain structure were calculated, and the estimated difference in slope with 95% confidence intervals was reported.

Outcome prediction

Hierarchical regression analysis was performed to determine if 3-month MRI-derived volumes contributed significantly to 12-month GOSE outcome after taking into account the established outcome predictors GCS score, PTA (PTA ≤ 2 weeks = 0, PTA > 2 weeks = 1) and age which were entered in the first step as the base model with which all subsequent models were compared. The MRI measures were cortical GM, hippocampal, lenticular nucleus and ventricular volumes at 3 months. In addition we used the Rotterdam CT scores (Maas et al., 2005) obtained from the worst CT scan, to see if this neuroimaging modality provided additional predictive value (Yuh et al., 2012). Overall variance explained (R2) and statistical significance were calculated for every model. R2-change and significance of F-value change were calculated as differences between the base model and all subsequent models.

Longitudinal analysis of ADC changes in GM

Intraclass correlation coefficients (ICC) using two way random single measures were used to calculate reliability of the ROI ADC values (Gwet, 2008).

A linear mixed effects model was used in the statistical comparisons of the  ADC values in the different ROIs over time in the entire TBI group, and between the TBI subgroups; moderate vs. severe TBI, PTA ≤ 2 weeks vs. >2 weeks, and DAI vs. non-DAI. The effect of degree of DAI (non-DAI = 0, DAI grades 1, 2, and 3) was assessed with ANOVA. Significant F-values were followed by paired Student's T-test or independent sample T-test. The resulting p-values were Bonferroni corrected as described earlier. Associations between ADC values and corresponding volumes at 3 and 12 months were examined using Pearson correlation.

A similar hierarchical regression analysis as described for the volumes above was performed to establish if cortical ADC values from the superior frontal and post-central sulci ROIs at 3 months were independent significant predictors of 12 month outcome.

Clinical and demographic data (age, GCS and GOSE scores) are given as mean ± standard deviation, and volumes and ADC values as mean ± standard error. The threshold for statistical significance was p < 0.05 for all analyses, including analyses applying Bonferroni correction for multiple comparisons.

Results

TBI patient characteristics (Table 1)

There were no significant differences with regard to the clinical variables between included and excluded TBI patients who survived for 12 months. Of the included patients 58% suffered moderate and 42% severe TBI. See Table 1 for details.

Table 1

TBI patient characteristics.

		All (n = 62)	Moderatea (n = 36)	Severeb (n = 26)
Age (years)		32.5 ± 15.2	34.3 ± 15.8	30.1 ± 14.4
GCS		9.1 ± 3.4	11.5 ± 1.6*	5.7 ± 1.8*
Mechanism of injury	Road accident	34 (55%)	19 (53%)	15 (58%)
	Fall injury	23 (37%)	13 (36%)	10 (38%)
	Other	4 (6%)	3 (8%)	1 (4%)
	Unknown	1 (2%)	1 (3%)	0
Post-traumatic amnesia	>2 weeks	19 (31%)	5 (14%)*	14 (54%)*
	≤2 weeks	39 (63%)	30 (83%)*	9 (35%)*
Diffuse axonal injury	No	19 (31%)	17 (47%)*	2 (8%)*
	Yes	41 (66%)	19 (53%)*	22 (85%)*
	Grade 1	13 (21%)	5 (14%)*	8 (31%)*
	Grade 2	17 (27%)	11 (31%)*	6 (23%)*
Outcome at 12 months	Grade 3	11 (18%)	3 (8%)*	8 (31%)*
With disability (GOSE < 7)		27 (44%)	12 (33%)*	15 (58%)*
Good recovery (GOSE 7–8)		35 (56%)	24 (67%)*	11 (42%)*
Days between injury and MRIc	Scan 1	10.2 ± 6.7	9.0 ± 6.2	12.1 ± 7.3
	Scan 2	98.4 ± 15.2	98.2 ± 14.3	98.6 ± 16.7
	Scan 3	373.5 ± 26.3	377.1 ± 30.1	369.0 ± 20.3

GCS 9–13.

GCS 3–8.

In days.

Significant difference between moderate and severe TBI.

Numbers of patients (percentages) or means ± standard deviations are shown. Percentages were calculated considering the total number of patients in each column.

GCS, Glasgow Coma Scale; GOSE, Glasgow Outcome Scale—Extended; MRI, magnetic resonance imaging.

Comparison between the included males (n = 45) and females (n = 17) showed that GCS score, HISS, and PTA and GOSE score at 12 months were similar between the sexes.

Injury mechanism, age and number of days between the head injury and the MRI scans were similar in the moderate and severe TBI patients. There were significantly more patients with PTA > 2 weeks, any grade of DAI and poor 12-month outcome in the severe than in the moderate TBI group.

Longitudinal brain structure volume changes in the entire TBI group (Figs. 4 and 5)

During the first 12 months after TBI significant volume differences were found for all analyzed structures (Fig. 4). Cortical GM experienced a rapid volume reduction and was significantly reduced at 3 months compared to the early scan, but did not change significantly after that. The hippocampal and lenticular nucleus volumes decreased more slowly and the volume was not significantly reduced before 12 months postinjury. Lobar WM volume was significantly reduced first at 12 months compared to both the early and the 3-month scan. Brainstem volumes decreased significantly between the early and 3-month scans and also between the 3-month and 12-month scans. The ventricular volume expanded from the early phase, but a significant increase was present first at 12-months. Statistical comparison of the interaction between slope of volume loss over time for the different brain structures demonstrated that cortical GM loss followed a significantly different trajectory compared to the hippocampus (estimated difference in slope was 1.2 (CI: 0.7–1.8), p < 0.0001), but not compared to lobar WM (estimated difference in slope was −0.04 (CI: −0.9–0.8); p = 0.9). The hippocampus and lenticular nucleus had similar trajectories of volume loss (estimated difference in slope was −0.01 (CI: −0.3–0.0); p = 0.1). For lobar WM and hippocampus significantly different trajectories of volume loss were present (estimated difference in slope was 1.2 (CI: 0.7–1.7); p < 0.0001).

Fig. 4

Longitudinal changes in mean ICV-corrected brain volumes with standard errors in the early phase (0–26 days post-injury, scan 1) and 3 months (scan 2), and 12 months (scan 3) postinjury. Volumes of the left and right sides of the structures were summed and evaluated together. ***p < 0.001, **p < 0.01, *p < 0.05.

Ventricular dilation was strongly associated with change in brainstem volumes at 3 months (R2 = 0.246, p = 0.001) and 12 months (R2 = 0.488, p < 0.001), and total cortical GM volume at 3 months (R2 = 0.143, p = 0.012) and 12 months (R2 = 0.331, p < 0.001). Ventricular dilation was not significantly associated with lobar WM volume changes at any time point.

The total brain volume decline was 3.8 ± 0.8% of ICV between the early and 3-month scans, and 5.7 ± 0.8% of ICV between the early and 12-month scans.

The impact of clinical TBI subtypes on brain volumes (Figs. 6 and 7)

TBI severity

Lobar WM and brainstem volume at 12 months were smaller in the severe compared to the moderate TBI group (lobar WM: 26.34 ± 0.58% ICV vs. 28.60 ± 0.51% ICV, p = 0.012; brainstem: 1.36 ± 0.03% ICV vs. 1.47 ± 0.02% ICV, p = 0.005) while cortical GM, ventricular, hippocampal and lenticular nucleus volumes were similar. The difference in lobar WM volume at 12 months was significant also after correcting for the presence of DAI (p = 0.012).

Moreover, hippocampal volume change between the early and 3-month scans was larger in the severe group compared to the moderate group (0.03 ± 0.01 %ICV vs. 0.00 ± 0.01%ICV, p = 0.046). For the other structures, there were no differences in volume changes over time between moderate and severe TBI.

PTA duration (Fig. 6a, b)

Patients with PTA > 2 weeks had smaller hippocampi at 3 and 12 months, and larger ventricles at 3 and 12 months compared to patients with PTA ≤ 2 weeks. Cortical GM volumes did not differ significantly between the PTA groups.

Hippocampal volume change over time was larger in the PTA > 2 weeks than in the PTA ≤ 2 weeks group between the early and 3-month scans (0.05 ± 0.01 %ICV vs. 0.00 ± 0.01 %ICV, p = 0.003). Cortical GM and ventricular volume changes over time did not differ between the PTA groups.

DAI groups (Figs. 6c, d and 7)

Smaller brainstem volumes were present at all time points in the DAI group compared to the non-DAI group. No differences in cortical GM, lobar WM or ventricular volumes were present between the DAI and non-DAI groups.

The non-DAI group had larger brainstem volumes than the DAI grade 2 group in the early phase; than the DAI grade 2 and 3 groups at 3 months, and the DAI grade 1, 2 and 3 groups at 12 months. In DAI 3 brainstem volumes decreased significantly between the early and 3 month scans, and further till 12 months. In DAI 2 brainstem volume decrease was present first after 12 months. In DAI 1 brainstem volumes decreased between the early scan and 12 months and between the 3 and 12 month scans. In the non-DAI group, brainstem volumes did not change with time.

No significant differences between the DAI groups were observed in ventricular volumes.

For the non-DAI, DAI 2 and DAI 3 groups the ventricular volume decreased between the early and 3 month scans, and between the early and 12 month scans. DAI 1 ventricular volume decreased between the early and 12 month scans.

Outcome prediction based on 3-month MRI volumes (Table 2)

The hierarchical regression analyses with PTA, GCS scores and age in the base model demonstrated that the MRI-derived brain volumes of cortical GM, hippocampus, and lenticular nucleus at 3 months explained significantly more of the final outcome than the base model (Table 2). The ventricular 3-month volume did not add to the 12-month outcome prediction, and neither did the Rotterdam CT score.

Table 2

Hierarchical regression analysis of 12-month outcome prediction using the 3-month brain structure volumes with GCS and PTA as base model.

Models	R2	R2-change	Significance of F change
Base model	0.385
CortGM 3 months	0.515	0.053	0.023
Ventricles 3 months	0.467	0.005	0.494
Hippocampus 3 months	0.535	0.073	0.007
Lenticular ncl. 3 months	0.529	0.067	0.010
Rotterdam CT score	0.465	0.003	0.596

Base predictors: GCS, PTA duration (0 = PTA ≤ 2 weeks; 1 = PTA > 2 weeks), and age

CortGM predictors: GCS, PTA duration, age, and cortical GM volume at 3 months

Ventricular predictors: GCS, PTA duration, age, and ventricular volume at 3 months

Hippocampal predictors: GCS, PTA duration, age, and hippocampal volume at 3 months

Lenticular nucleus predictors: GCS, PTA duration, age, and lenticular nucleus volume at 3 months

Rotterdam CT score predictors: GCS, PTA duration, age, and Rotterdam CT score

CortGM, cortical gray matter; GCS, Glasgow Coma Scale; PTA, posttraumatic amnesia duration.

Analysis of ADC value changes in GM

Intra-rater reliability of ADC measurements in the ROIs

ICC of repeated measurements were 0.74 for superior frontal sulcus, 0.86 for postcentral sulcus, 0.69 for insular gyrus, 0.88 for hippocampus and 0.89 for brainstem.

Overall changes in entire TBI group

No ADC value changes over time were found in any ROI in the entire TBI group.

The impact of clinical TBI subtypes on ADC value changes

TBI severity (Table 3a)

The severe TBI group had higher ADC values in the brainstem at 3 months, the superior frontal sulcus at 3 and 12 months, and the postcentral sulcus at 12 months compared to the moderate TBI group. No significant differences in ADC values between moderate and severe TBI groups were found at any time point for the insula or hippocampus.

In the severe TBI group the ADC values increased over time in the superior frontal sulcus from the early phase to 3 months (p = 0.043), and the postcentral sulcus from the early phase to both 3 and 12 months (early phase < 3 months, p = 0.028; early phase < 12 months, p = 0.027). There were no differences in brainstem, insula or hippocampal ADC values between any time points. There were no differences in ADC values over time in any ROI in the moderate TBI group.

PTA duration (Table 3b)

At 12 months the ADC values in the superior frontal sulcus were higher in patients with PTA > 2 weeks than in those with PTA ≤ 2 weeks, while no differences were found in the other ROIs. In the PTA > 2 weeks group ADC values in insula increased between 3 and 12 months (p = 0.046). In the group with PTA ≤ 2 weeks no differences in ADC values overtime were found in any ROI.

DAI groups (Figs. 8 and 9)

In the superior frontal sulcus and postcentral sulcus ADC values were higher in TBI patients with DAI than in non-DAI patients at 3 and 12 months. The ADC values increased significantly from the early phase to 3 months in the superior frontal (p = 0.009) and the postcentral sulcus (p = 0.005) in the DAI group. No differences were found in ADC values or ADC value change over time in the insula, hippocampus or brainstem ROIs between the non-DAI and the DAI groups, or within the DAI grades.

The different DAI grades were associated with significant differences in cortical GM ADC values (Fig. 9). In DAI 3, the superior frontal sulcus ADC values were higher than those in the non-DAI group at 3 months, and compared to both non-DAI and DAI 1 at 12 months. In DAI 2 the ADC values were higher in the superior frontal sulcus at 3 and 12 months than those in the non-DAI group. The DAI 3 group had higher ADC values in the postcentral sulcus at both 3 and 12 months than non-DAI and DAI 1 groups. In the superior frontal sulcus ADC values decreased between 3 and 12 months (p = 0.018) in DAI 1, and increased between the early phase and both 3 and 12 months in DAI 2 (early phase < 3 months, p = 0.036; early phase < 12 months, p = 0.018). There were no significant changes over time within the non-DAI and the DAI 3 groups.

Fig. 9

The associations between longitudinal changes in mean ADC values for the non-DAI group (non-DAI = DAI 0) and the three DAI groups (DAI grades 1, 2, 3). Longitudinal analysis within each DAI group was performed separately for both ROIs using a linear mixed model. Significant F-values were followed by the Student's T-test. The resulting p-values were Bonferroni corrected separately for every ROI. At each time point, differences between DAI subgroups were evaluated using one-way ANOVA, and significant F-values were followed by independent samples T-tests, the resulting p-values were Bonferroni corrected for the three time points separately. Significant differences between different time points for each DAI group were confirmed and marked as: @—significantly different from scan 2; $—significantly different from scan 3 (***p < 0.001, **p < 0.01, *p < 0.05).

The ADC values increased between the early phase and 3 months in the DAI 2 group (p = 0.020) in the postcentral sulcus. In DAI 3, the ADC values were significantly increased both at 3 and 12 months compared to the early phase (early phase < 3 months, p = 0.003; early phase < 12 months, p = 0.003). ADC values in postcentral sulcus did not change with time in the non-DAI and DAI 1 groups.

ROI ADC values vs. brain structure volumes.

There were no correlations between the ADC values in the cortex, hippocampus or brainstem and the volumes of these structures at 3 and 12 months.

ADC and outcome prediction

ADC values of superior frontal and postcentral sulci from 3 months did not add to the 12-month outcome prediction (Table 4).

Table 4

Hierarchical regression analysis of 12-month outcome prediction using the 3-month ADC values with GCS, PTA and age as base model.

Models	R2	R2-change	Significance of F change
Base model	0.462
Superior frontal sulcus model	0.479	0.017	0.208
Postcentral sulcus model	0.480	0.018	0.193

Base predictors: GCS, PTA duration, and age.

Superior frontal sulcus model: GCS, PTA duration, age, and superior frontal sulcus ADC value at 3 months

Postcentral sulcus model: GCS, PTA duration, age, and postcentral sulcus ADC value at 3 months

ADC, apparent diffusion coefficient; GCS, Glasgow Coma Scale; PTA, posttraumatic amnesia duration.

Discussion

This is to the best of our knowledge the largest number of moderate and severe TBI patients included in a prospective, longitudinal MRI study of brain volumes and diffusion changes where data have been systematically collected at 3 different time points, i.e. early phase and 3 and 12 months postinjury. Moreover, the patient group included in this study was representative for the entire group of TBI patients surviving for 1 year. The results not only confirmed most of our predictions, but also revealed new, unexpected findings.

Longitudinal volume changes

Significant cortical GM volume loss took place between the early and 3-month scans without any further significant decline at 12 months. This is the first direct evidence for cortical GM volume loss being an early occurrence after TBI in humans, similar to results in animal models (Rodriguez-Paez et al., 2005; Smith et al., 1997; Liu et al., 2010). The rapid cortical volume loss probably stems from a combination of neuronal necrosis and apoptosis initiated at time of injury and peaking early postinjury (Conti et al., 1998). Supporting these observations in animals, human PET studies show early and widespread neurodegeneration in intact cortex following TBI (Wu et al., 2004; Bergsneider et al., 2001). We did not find any evidence for GCS scores, duration of PTA or presence of DAI being associated with cortical volume loss. The effect of DAI and PTA on cortical volumes may, however, be localized to specific cortical regions for instance connected to affected WM tracts or the medial temporal lobe, and hence not be reflected in total cortical volumes.

The volume loss in the hippocampus was significantly slower than that in the cortex and not significantly different from the early and 3 month scans before 12 months postinjury. This slow, continuous loss of hippocampal volume after TBI in humans concurs with animal studies (Immonen et al., 2009; Smith et al., 1997). However, both animal data (Immonen et al., 2009; Smith et al., 1997) and autopsy data from humans (Maxwell et al., 2003) demonstrate significant early neuronal loss in the hippocampus. The lack of a significant decline in hippocampal volume between the early and 3 month-scans in the current study may be due to the early scans being obtained at ~10.2  days after injury. At that time point some volume loss may already have taken place. Alternatively, MRI-based hippocampal volume loss may take longer to develop than histological identified cell loss. Still, our result concurs with findings in cross-sectional studies where decreased hippocampal volume is detected first after 3 months in moderate–severe TBI patients compared to matched controls (Wilde et al., 2007; Bigler et al., 1997). The current data also agree with results from longitudinal studies with two time points, which show that hippocampal volumes decline over a protracted period from 1 week to 2.5 years (Warner et al., 2010b; Ng et al., 2008). The slowly evolving hippocampal volume loss is consistent with the protracted apoptotic neuronal death described in this region (Williams et al., 2001), coupled with glial proliferation (Gao and Chen, 2013) which may counteract some of the volume change, and impaired neurogenesis (Gao et al., 2008) which may affect volume when accumulated over time. The hippocampus is critical for declarative memory (Squire et al., 2004; Eichenbaum et al., 2007), and hippocampal damage is known to produce symptoms of anterograde and retrograde amnesia (Cipolotti and Bird, 2006). It was thus not surprising to find that patients with PTA > 2 weeks had significantly lower hippocampal volumes at both 3 and 12 months than patients with PTA ≤ 2 weeks. Moreover, hippocampal volume loss was markedly more pronounced from the early to 3 month scans in the PTA > 2 weeks group. These results demonstrated that PTA reflects hippocampal neuronal injury and that more severe PTA is associated with more rapid and extensive neuronal death in this region.

Significant volume loss developed over 12 months in the lenticular nucleus. Decreased lenticular nucleus volume has been reported in cross-sectional and longitudinal studies of TBI in children and adults ≥1 year after injury (Wilde et al., 2007; Anderson and Bigler, 1995; Primus et al., 1997). The histopathological changes in the lenticular nucleus include apoptosis, neuronal loss and/or shrinkage due to transneuronal degeneration and loss of myelinated axons as described in animal studies (Tong et al., 2002; Ng et al., 1994). Based on the current results these processes take place over a protracted period. The present data do not support the notion that basal ganglia respond more rapidly and briefly to TBI than the hippocampus (Hicks et al., 1996). Rather, the current study showed that the trajectory of volume loss was similar in the subcortical GM structures, i.e. the lenticular nucleus and hippocampus. Moreover, the subcortical GM structures were shown to differ significantly from the cortex with regard to the slopes of volume loss suggesting important differences in the histopathological responses between these GM structures to TBI.

As predicted lobar WM volume loss decreased mainly in the late phase, being significant first at 12 months. This concurs with the slowly evolving pathological response in WM (Vargas and Barres, 2007; Coleman, 2005; Coleman and Freeman, 2010). The extended period of volume loss found in lobar WM in humans in the current study is most likely the macroscopic consequence of the protracted histological changes in WM. Previous longitudinal studies have only been able to show a significant decline in WM volumes between an early (1–2 months) and a later time point (6–12 months) after mild to severe TBI (Ding et al., 2008; Sidaros et al., 2009; Warner et al., 2010b). Surprisingly the trajectory of WM loss and cortical GM loss did not differ significantly although the slope of WM loss was significantly different from the hippocampus, lenticular nucleus, and brainstem. This finding could reflect a significant interaction between WM volume loss and cortical volume loss over time at the individual level. Another explanation could be that volume losses in cortical GM and lobar WM volumes were the least linear. The lack of a significant interaction between volume losses over time between these two regions could hence be due to the model not being able to take nonlinearities into account.

Lobar WM volume loss was associated with injury severity as patients with severe TBI had significantly reduced lobar WM volumes compared to the moderate TBI group. This difference was present also after correcting for the presence of DAI. Unexpectedly, DAI was not associated with decreased lobar WM volume. This may be due to DAI classifications being based on visible lesions on conventional MRI which may not fully describe WM involvement in TBI.

The brainstem was the only structure in which volume loss was significant between every scan session, thus demonstrating particularly pervasive and enduring consequences of TBI to this region. Previous MRI studies have reported brainstem volume loss in cross-sectional and longitudinal studies with two time points in both children and adults with moderate and/or severe TBI (Sidaros et al., 2009; Warner et al., 2010b; Fearing et al., 2008). WM injury probably contributed to the early brainstem volume decline since patients with DAI grade 3 had the fastest and largest reduction of brainstem volume, present from 3 months. Also the group with severe TBI had smaller brainstem volumes, but this was detectable first after 12 months, appearing at the same time as reduced lobar WM volumes in this group. The WM tracts in the brain stem thus appear to be particularly sensitive to TBI.

Ventricular volumes were significantly increased at 3 months, followed by a non-significant increase at 12 months. This result is a refinement of the previously described time frame of significant ventricular dilation which was shown to last until ˜2–7 months postinjury (Blatter et al., 1997). Ventricular dilation at 3 and 12 months was associated with brainstem and cortical GM volumes. Neither lobar WM volumes nor the presence of DAI was associated with larger ventricular volumes. These results demonstrate the importance of deep WM loss in particular, followed by cortical GM to ventricular dilation in TBI, rather than hemispheric WM loss.

For the entire TBI group mean brain atrophy was ~4% of ICV after 3 months and ~6% of ICV after 12 months, demonstrating that the bulk volume loss occurred within the first 3 months. Still, a notable protracted component to the atrophy was present. No directly comparable human studies exist, but endpoint brain volume atrophy in the chronic phase after TBI is reported to be from 1.4% to 8% in cross-sectional and longitudinal studies in mild to severe TBI (Trivedi et al., 2007; Ding et al., 2008; Sidaros et al., 2009).

Longitudinal ADC changes

ADC values increased over time in the individuals with the most serious brain injuries, i.e. severe TBI, PTA > 2 weeks, and DAI, and there were no associations between ADC values and volumes. Participants with severe TBI had higher ADC values in deep, radiologically unaffected sulcal cortical GM, which increased significantly throughout the 12 month period, even when total cortical volume remained similar between 3 and 12 months. Thus it appears that small-scale remodeling of cortex after TBI persisted for a longer period than volumetric changes. Moreover, a significant contribution from DAI to the increase in cortical ADC values was present. More severe DAI grades gave rise to increased cortical ADC values developing at a more rapid timescale. In DAI the cortical ADC changes may be connected to delayed neuronal death, soma shrinkage, loss of cortical axons and myelination developing over an extended time period in this condition (Goetz et al., 2004; Hergan et al., 2002).

As expected there were increased ADC values at 3 months in the brainstem in severe TBI, but not in moderate TBI. This finding lends further support to the notion that deep brain injury is significantly more pronounced in severe than moderate TBI. PTA > 2 weeks was not associated with increased hippocampal ADC values. Both increased and decreased hippocampal ADC values are reported after TBI in animals and humans (Bendlin et al., 2008; Obenaus et al., 2007; Van putten et al., 2005). We predicted that PTA > 2 weeks would lead to increased ADC values due to increased cell loss, but did not find support for this assumption in the data. It is possible that the pathophysiological mechanisms leading to PTA are more heterogeneous than we expected and/or that volume loss combined with glial proliferation, inflammation (Woodcock and Morganti-Kossmann, 2013; Lin and Wen, 2013; Giunta et al., 2012) and/or neurogenesis (Gao et al., 2008; Greer et al., 2011; Gao and Chen, 2013) counteracted any effect of neuronal loss on ADC values.

MRI-derived 3-month volumes for 12-month outcome prediction

The hierarchical regression analysis demonstrated that MRI-derived volumes at 3 months contributed independently to 12-month outcome prediction even adjusting for duration of PTA, GCS scores and age. This finding demonstrates a significant clinical potential for MRI volumes obtained at 3 months for outcome prediction. Rotterdam CT score, on the other hand, did not have additional predictive value. Furthermore, cortical ADC values at 3 months did not improve outcome prediction.

Methodological considerations

This is the first study using a fully automated method, NeuroQuant, for deriving brain volumes from MRI scans in patients with TBI. Our results show that NeuroQuant performed very well in moderate–severe TBI patients even in the presence of different types of injuries and postoperative changes. Its limitations were restricted to segmentation problems in the early phase due to hematomas and contusions, and excessive movement during scanning. For the 3-month scans no segmentations failed. Since several volumes derived from the 3-month scan were independent predictors of outcome at 12 months, MRI at 3 months appears to be of clinical value in the assessment of the brain injury and prediction of later function. However, for the individual patient the presence of focal injury (contusions, hemorrhages) also plays a significant role for outcome, but it was beyond the scope of the present paper to evaluate the interaction between focal injuries and the atrophy in different brain structures.

It should be noted that the volume change measures obtained using independent NeuroQuant segmentations would not be expected to have the spatial specificity and power to detect subtle change that many across time point registration methods might provide (Holland et al., 2009; Fox et al., 2001; Thompson et al., 2004). On the other hand, NeuroQuant has previously been shown to have results comparable to that of hand segmentation of subregional volumes by anatomical experts (Brewer et al., 2009), and thus provides an identifiable and translatable measure of structure volume. As a FDA-cleared tool used at more than 100 clinical sites and applied to over 30,000 individuals in routine clinical practice, the tool yields standardized data that are directly comparable to results obtained in the clinic. Further, by providing independent segmentations of each time point, this approach is inverse-consistent and not subject to registration bias that may inflate measures of change provided by registration-based techniques (Holland et al., 2012). Thus, the approach yields a conservative measure of the regional brain volume changes associated with TBI.

Reliability of the ADC measurements was high to very high, but only performed by a single rater and interrater variability could not be assessed as Ozturk et al. (2008) recommend. There are several ways of analyzing ADC values, and we choose manually positioned ROIs to try to optimize ROI localization to avoid partial volume effects as slice thickness was 5 mm. However, more isotropic voxels and/or smaller voxels would be preferable to avoid partial volume effects, which potentially could have made it easier to detect statistical significant differences between groups and over time. For all the GM ROIs, the effect sizes of the differences between ADC values for the different TBI groups were very large. This was also seen in the regions where there were no statistical significant differences. Taken together, our results show that that variability in ADC measures is substantial in the material, and that ROI positioning probably plays some part in this, but that the largest contribution comes from individual differences within and between subjects. ADC values in GM may change slightly during normal aging. However, according to earlier studies, brain ADC values of healthy subjects remain constant during the age period 20–60 years (Klimas et al., 2013; Watanabe et al., 2013). As most of the subjects included in the present study are in this age range, we assumed that the diffusion changes revealed are TBI-associated. Including a matched control group may have helped in the interpretation of the ADC values, which relies on inferences from histopathology. Newer diffusion based MRI methods such as restriction spectrum imaging (White et al., 2013) may be better suited for characterizing tissue microstructural changes in TBI in the future.

The large number of analyses conducted across multiple brain regions increased the susceptibility to type I errors in the current study.

5. Conclusion

In summary, the effect of TBI on brain volumes differed between brain structures, and injury subgroups (GCS score, PTA, DAI) were associated with specific patterns of volume loss. Changes in cortical diffusion properties were detected in the more severe TBI subgroups in particular in patients with DAI grade 3, where diffusion properties increased in value and encompassed more cortical regions with time. Volume changes and diffusion changes were not associated. Cortical GM, hippocampus and lenticular nucleus volumes at 3 months were all very significant independent predictors of 12-month outcome even after adjusting for PTA duration, GCS scores and age. Thus TBI induced early and widespread structural brain volume loss which has significant impact on outcome in the chronic phase in moderate and severe TBI patients.

Table 3

Longitudinal ADC values of TBI subgroups and comparison between the clinical TBI subtypes.

a) Moderate vs. severe TBI.
		Moderate TBI	Severe TBI	p-Value
Superior frontal sulcus
	Early phase	68.2 ± 1.1	68.1 ± 1.3	NS
	3 months	66.9 ± 1.0	71.5 ± 1.2	0.011
	12 months	65.0 ± 1.1	70.7 ± 1.2	0.005
Postcentral sulcus
	Early phase	70.2 ± 0.9	69.6 ± 1.1	NS
	3 months	69.9 ± 0.9	73.2 ± 1.0	NS
	12 months	68.9 ± 1.0	73.2 ± 1.0	0.009
Insular gyrus
	Early phase	81.5 ± 1.1	81.9 ± 1.3	NS
	3 months	80.7 ± 1.0	81.9 ± 1.2	NS
	12 months	81.4 ± 1.1	83.9 ± 1.2	NS
Hippocampus
	Early phase	86.1 ± 1.4	88.0 ± 1.7	NS
	3 months	85.8 ± 1.3	87.8 ± 1.5	NS
	12 months	85.1 ± 1.4	88.7 ± 1.5	NS
Brainstem
	Early phase	80.6 ± 1.7	81.5 ± 2.1	NS
	3 months	80.5 ± 1.6	87.5 ± 1.9	0.016
	12 months	80.3 ± 1.7	83.5 ± 1.9	NS
b) Patients with PTA ≤ 2 weeks vs. >2 weeks
		PTA ≤ 2 weeks	PTA > 2 weeks	p-Value
Superior frontal sulcus
	Early phase	68.2 ± 1.0	67.9 ± 1.5	NS
	3 months	68.0 ± 1.0	70.4 ± 1.5	NS
	12 months	66.3 ± 1.1	70.8 ± 1.5	0.037
Postcentral sulcus
	Early phase	70.0 ± 0.9	69.8 ± 1.4	NS
	3 months	70.4 ± 0.9	72.9 ± 1.3	NS
	12 months	69.9 ± 1.0	72.6 ± 1.3	NS
Insular gyrus
	Early phase	81.7 ± 0.9	81.5 ± 1.4	NS
	3 months	81.9 ± 0.9	79.5 ± 1.3	NS
	12 months	81.8 ± 1.0	84.0 ± 1.3	NS
Hippocampus
	Early phase	87.0 ± 1.1	86.2 ± 1.6	NS
	3 months	86.0 ± 1.1	85.0 ± 1.5	NS
	12 months	86.1 ± 1.1	87.8 ± 1.5	NS
Brainstem
	Early phase	80.8 ± 1.5	81.3 ± 2.2	NS
	3 months	81.8 ± 1.4	83.8 ± 2.1	NS
	12 months	81.2 ± 1.5	82.6 ± 2.1	NS

ADC in 10−5 m2/s. Values given as means ± standard deviations.

Early phase scans were acquired 0–26 days postinjury.

Statistical differences between moderate vs. severe or PTA ≤ 2 vs. >2 weeks were Bonferroni-corrected for every structure separately.

ADC, apparent diffusion coefficient; NS, non-significant; TBI, traumatic brain injury.