Abstract
BACKGROUND AND PURPOSE: Diffusion tensor imaging (DTI) can noninvasively detect in vivo white matter (WM) abnormalities on the basis of anisotropic diffusion properties. We analyzed DTI data retrospectively to quantify the abnormalities in different WM regions in children with hydrocephalus during early infancy.
MATERIALS AND METHODS: Seventeen infants diagnosed with hydrocephalus (age range, 0.13–16.14 months) were evaluated with DTI and compared with 17 closely age-matched healthy children (age range, 0.20–16.11 months). Fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity, and radial diffusivity values in 5 regions of interest (ROIs) in the corpus callosum and internal capsule were measured and compared. The correlation between FA and age was also studied and compared by ROI between the 2 study groups.
RESULTS: Infants with hydrocephalus had significantly lower FA, higher MD, and higher radial diffusivity values for all 3 ROIs in the corpus callosum, but not for the 2 ROIs in the internal capsule. In infants with hydrocephalus, the increase of FA with age during normal development was absent in the corpus callosum but was still preserved in the internal capsule. There was also a significant difference in the frequency of occurrence of abnormal FA values in the corpus callosum and internal capsule.
CONCLUSIONS: This retrospective DTI study demonstrated significant WM abnormalities in infants with hydrocephalus in both the corpus callosum and internal capsule. The results also showed evidence that the impact of hydrocephalus on WM was different in the corpus callosum and internal capsule.
Hydrocephalus is a pathologic condition in which excessive CSF accumulates in the ventricular system because of either obstruction along the CSF pathways or an imbalance between CSF production and reabsorption.1 The enlarged ventricles and the associated increased intracranial pressure (ICP) can cause significant damage to various regions of brain, especially to the adjacent white matter (WM) tissue.1–4 The mainstay of treatment for most hydrocephalic patients has been surgical diversion of the excess CSF to a distant location in the body. Although surgical outcomes are generally good, some patients are still at risk for cognitive, motor, and physical developmental delays.5–10 The variability of outcomes in hydrocephalus treated in infancy may reflect the wide spectrum of injury to specific regions of the brain and the variety of recovery mechanisms in the pathophysiology in these locations.
Ventricular size, CSF flow, and ICP are routinely used to guide the treatment of hydrocephalus and have been the basis for diagnostic standards and attempts to predict prognosis. However, neither these nor conventional imaging modalities have been found to be completely accurate. The pathophysiologic mechanism of brain injury has not been clearly elucidated, but measurement of diffusion parameters may provide insight into the mechanisms and/or reversibility of WM injury in childhood hydrocephalus.
Diffusion tensor imaging (DTI) provides quantitative information about anisotropic diffusion properties in WM and has been applied to investigate in vivo WM damage and possible recovery in various neurologic and pathologic disorders.11–13 However, to our knowledge, there are very few published articles on the use of DTI in childhood hydrocephalus. Although Assaf et al14 studied abnormal anisotropic diffusion properties in various WM regions before and after CSF shunt surgery, they did not include patients in early childhood, the common age of hydrocephalus presentation and treatment. In a similar fashion, in a recent DTI study by Hasan et al15 of children with spina bifida and hydrocephalus, the average age of participants was 12.3 ± 2.1 years, well beyond the usual age of diagnosis and treatment of childhood hydrocephalus.
In our study, we used DTI to study WM integrity in infants with hydrocephalus to assess the anisotropic diffusion properties in the corpus callosum and internal capsule preoperatively. We hypothesized that 1) these WM structures would demonstrate abnormal anisotropic diffusion values (FA, MD, axial and radial diffusivity) compared with age-matched healthy children, and 2) the abnormality would be region-specific in the direction and the degree of abnormality.
Materials and Methods
Patient Population
We retrospectively reviewed existing clinical DTI datasets and identified 2 groups of participants for our study: a preshunt hydrocephalus group and an age-matched control group. The Institutional Review Board of Cincinnati Children's Hospital Medical Center approved the study.
The preshunt hydrocephalus group consisted of 17 infants (age range, 0.03–16.14 months; age mean ± SD, 4.65 ± 4.27 months; sex ratio, 7 girls/10 boys) who were diagnosed with hydrocephalus and had MR imaging and DTI performed before shunt surgery as part of their standard clinical care. Demographics and clinical information for these patients are summarized in Table 1. Fifteen of these patients were described as initially presenting with symptoms of accelerated head growth, macrocephaly, enlarged ventricle, or ventriculomegaly. All of the patients demonstrated clinical improvement when evaluated postsurgically: all demonstrated neurologically stable examination or had normal development.
Demographic and clinical information for pediatric patients with HCP
The control group consisted of 17 closely age-matched children (age range, 0.20–16.11 months; age mean ± SD, 4.71 ± 4.17 months; sex ratio, 6 girls/11 boys). These children were selected from a cohort of an ongoing project that aims to establish a frame of reference for the anisotropic diffusion properties throughout development (total n is approximately 250, 0–18 years old; n = 45 for age range 0–16 months). All of the children in the control group met the following criteria: 1) they were scanned for non-central nervous system (CNS)–related problems, 2) they had no previous record of neurologic disorder, 3) they had normal MR imaging results, and 4) they had no record of neurologic disorder for at least 4 months after the MR imaging/DTI scan. Because DTI parameters usually change with age during early childhood, 3 weeks were used as the criteria for the maximal age difference in searching for matching healthy children from the data base. We were only able to find 1 match for each child in the patient group. The age difference in these 17 pairs of children ranges from 0 to 16 days (mean ± SD, 5.35 ± 4.57 days). The 2 groups were not significantly different in age (2-tailed paired t test; P = .32) or sex ratio (Fisher exact test; P = 1).
The fronto-occipital horn ratio (FOHR) was used to assess ventricular size for participants in both groups. This index measures the ratio between the mean of the frontal and occipital horn width to the width of parietal lobe and has been found to correlate well with relative ventricular size.16 Figure 1A shows the methodology for the measurement of FOHR.
A, The measurement and calculation of the FOHR are demonstrated on an axial section of T1-weighted MR imaging study from a child with hydrocephalus in the preshunt group. B, Histogram showing that children in the control group demonstrate a narrow range (0.29∼0.37) in the FOHR, whereas children in the patient group have a distinct range of distribution of the ratio (0.40∼0.65).
Periventricular interstitial edema may be used as a radiologic sign of acute and severe hydrocephalus in children. In our study, only 3 infants with hydrocephalus were found to have abnormal periventricular T2 or fluid-attenuated inversion recovery signals. The number was not sufficient to conduct any statistical comparison and was therefore not examined further. The relative higher water content in young infants may contribute to the low incidence of detectable interstitial edema in periventricular WM in our patient population.
MR Imaging/DTI Scan
All of the MR imaging/DTI scans were performed at Cincinnati Children's Hospital Medical Center between October 2003 and May 2008. Images were acquired clinically either on a 3T Magnetom Trio scanner (Siemens, Erlangen, Germany) or a 1.5T Signa Horizon LX scanner (GE Healthcare, Milwaukee, Wis). All of the DTI images were acquired with diffusion-weighted, spin-echo echo-planar imaging in the axial plane with a b-value of 1000 s/mm2. In the 17 infants with hydrocephalus, 3 infants were scanned on the 3T Siemens scanner with a 12-direction DTI protocol: TR, 6000 ms; TE, 87 ms; resolution, 2 × 2 mm; and section thickness, 2 mm. The other 14 patients were scanned on a 1.5T GE scanner with a 15-direction DTI protocol: TR, 12,000 ms; TE, 81 to 101ms; and resolution, 3 × 3 mm (n = 11) or 1.88 × 1.88 mm (n = 3). In the 17 age-matched control subjects, the parameters were even more inhomogeneous, perhaps because of the variety of protocols used for these children who were referred for MR imaging scanning for more diverse reasons. Of the 17 healthy children, 9 were scanned on the 3T Siemens scanner. Among them, 6 were scanned with a 12-direction DTI protocol: TR, 6000 ms; TE, 87 ms; resolution, 2 × 2 mm; and section thickness = 2 mm. The other 3 healthy children were scanned on the 3T scanner with a 6-direction DTI protocol: TR, 4100 or 5300 ms; TE, 84 ms; resolution, 1.56 × 1.56 mm (n = 1) or 1.72 × 1.72 (n = 2); and section thickness, 4 mm (n = 2) or 3 mm (n = 1). Of the 17 healthy children, 8 were scanned on a 1.5T GE scanner: gradient directions, 15; TR, 12,000 ms; TE, 66 to 97 ms; resolution, between 2.5 × 2.5 and 3 × 3 mm; and section thickness, 3 mm. A matrix of 128 × 128 was used for all patients. The differences in the DTI scan protocols were because of the occasional adjustment in clinical scanning for image quality optimization.
We also conducted a variance ratio test for each pair on the basis of the observed SD and corresponding degree of freedom, with the goal to examine whether there was a systemic image quality bias for the FA value measurement. The number of pairs that demonstrate significant difference at a P level of .05 was found to be small (2/16, 1/7, 1/14, 0/17, and 2/17 for the 5 respective ROIs, respectively). An additional analysis of the effect of field strength and protocol variations on the measured diffusion values found no statistically significant differences for any of the ROIs tested and did not affect the results and conclusion of this study. This issue has been elaborated and discussed elsewhere.17
Data Processing
Image reconstruction, postprocessing, and ROI-based DTI parameter calculations were performed with the software DTIStudio 2.4 software.18 On a pixel-by-pixel basis, the 6 elements (Dxx, Dyy, Dzz, Dxy, Dxz, and Dyz) were calculated and diagonalized to calculate the 3 eigenvalues (λ1, λ2, λ3) corresponding to the 3 eigenvectors in the diffusion tensor matrix. We calculated the FA map using the following formula19,20:
The mean diffusivity (MD) was calculated as the mean of the 3 eigenvalues (MD = (λ1 + λ2 + λ3)/3). Two specific eigenvalue indices, λ‖ (= λ1) and λ⊥(= (λ2+λ3)/2) that represent diffusion properties along axial and radial directions, respectively, are also studied.
On color-coded FA maps, 5 WM regions (Fig 2A and B) were delineated for each participant: 1) genu of the corpus callosum (gCC); 2) body of the corpus callosum (bCC); 3) splenium of the corpus callosum (sCC); 4) anterior limb of the internal capsule (ALIC); and 5) posterior limb of the internal capsule (PLIC). These are all major WM structures that can be easily identified on color-coded FA maps on the basis of their anatomic location and the knowledge of their orientation. Fibers in left-right (eg, gCC, bCC, and sCC), superior-inferior (eg, PLIC), and anteroposterior (eg, ALIC) direction are conventionally coded as red, blue, and green, respectively. We first randomly selected a healthy child and manually drew all of the ROIs on color-coded FA maps as shown in Fig 2. Then the ROIs identified in this subject were used as a guide to manually define ROIs for other subjects as reproducibly as possible. The delineation of these ROIs followed the approach by Hermoye et al21 and has also been described in our previous work.22 In both ALIC and PLIC, because we did not find statistically significant differences between the left and the right side, we used the average DTI value in the analysis. Anatomic distortion caused by hydrocephalus made some ROIs undefinable in some patients. The corpus callosum is expected to be thinner in hydrocephalus. Because the ROIs were all defined on axial images, partial volume effect is likely to occur when CSF is included in the delineation of a very thin bCC. To minimize the potential impact, we determined that the bCC of a patient would be excluded from analysis if the corpus callosum (measured on sagittal T1 images) was thinner than the section thickness (mean ± SD, 3.6 ± 0.65 mm after excluding bCC measurements in 7 patients with hydrocephalus).
Color-coded FA maps for a child in the healthy control group (A and B) and a child with hydrocephalus in the preshunt group (C). The orientation of white matter tracts is coded with different colors with red, green, and blue indicating left-right, anteroposterior, and superior-inferior, respectively. The ROIs used in the study are outlined for the gCC, L_ALIC and R_ALIC, L_PLIC and R_PLIC and sCC in (A), and the bCC in (B). In both ALIC and PLIC, because there were no statistically significant differences between the left and the right side, the average values of DTI measurements were used in the analysis.
Statistical Analysis
We performed statistical analysis using SPSS, version 15 (SPSS, Chicago, Ill). The statistical differences between the 17 children with hydrocephalus and their age-matched control subjects were tested with the paired t test on all the DTI parameters in various ROIs. To control for the expected proportion of incorrectly rejected null hypotheses (type I error rate) in multiple comparison, we made a correction for multiple comparisons using the false discovery rate method.23
As reported in the literature,21,24,25 the developmental trajectory of FA in healthy children can often be modeled by a monoexponential or a biexponential curve. The most drastic increase in values occurs in the first 24 months of life, with the values leveling off and stabilizing before 36 months. The maximal age in our study group was 16 months, and the full age-range over which the normal curve occurs was not represented in this cohort and curve fitting with use of an exponential model was not appropriate. Therefore, the increase of FA with age was fitted with a linear model. The 95% prediction interval was also calculated.
For each ROI, the value of FA was objectively determined to be abnormally high, normal, or abnormally low on the basis of whether the DTI index was above, within, or below the prediction interval at 95% confidence level as derived from the regression analysis in the normal group. The term frequency of occurrence is used to reflect how often the abnormal DTI measurement occurs in a certain cohort. It is equivalent to the percentage of individuals who can be categorized to a certain subgroup. The Freeman-Halton extension26 for the Fisher exact test was used to evaluate the 2 × 3 contingency table and to assess the statistical significance (at a level of P = .05) of the different frequency of occurrence in abnormal FA either across ROIs or across subject groups.
Results
Comparing FOHR between Infants with Hydrocephalus and Age-Matched Control Subjects
The FOHR for the control group followed a normal distribution, with a range from 0.285 to 0.374 (mean ± SD, 0.336 ± 0.023, Fig 1B). The FOHR of the patients with hydrocephalus (range, 0.395–0.648; mean ± SD, 0.530 ± 0.065) demonstrated a wider range of values, with the minimum close to the maximum ratio seen in the control group (paired t test, P < .0001). No statistically significant correlation was found between FOHR and DTI parameters in any of the 5 ROIs examined.
Comparison of DTI Parameters between Infants with Hydrocephalus and Age-Matched Control Subjects
FA, MD, axial diffusivity, and radial diffusivity values for the 2 study groups are presented in Table 2.
Comparison of DTI parameters between children with HCP and age-matched control subjects
FA values in the corpus callosum in children with hydrocephalus were significantly lower than that in age-matched control subjects in all 3 ROIs (2-tailed paired t test controlled for multiple comparison, P < .003, P < .03, and P < .03 for gCC, bCC, and sCC, respectively). MD values in the 3 ROIs in children with hydrocephalus were all higher than those in age-matched control subject. This difference was statistically significant in the gCC (P < .05) and bCC (P < .05); the difference in the splenium showed a trend but did not reach statistical significance (P = .17). No statistically significant differences in axial diffusivity were demonstrated for any ROI in the corpus callosum between the 2 groups. In contrast, all 3 ROIs in the corpus callosum showed significantly higher radial diffusivity in patients compared with age-matched control subjects (P < .03, P < .03, and P = .05 for gCC, bCC, and sCC, respectively). No statistical difference in mean DTI parameters was found between the 2 groups for either the ALIC or PLIC.
Correlation between FA and Age
Figure 3 shows the linear regression between FA and age for sCC (Fig 3A) and PLIC (Fig 3B). In general, FA values of children in the control group increased with age in all 5 ROIs, consistent with the published literature.21,24,25 Within the age range in our study, this increase followed a linear pattern with statistical significance in healthy children for all ROIs (P < .001). The curve-fitting coefficients and the other statistical results are listed in Table 3.
A, sCC. B, PLIC. Scatterplots of FA fitted to a linear model (solid line) in healthy children and in children with hydrocephalus (regression line not shown). The 95% prediction intervals are shown in dashed lines representing age-specific range in which a normal value is supposed to be at 95% certainty.
Results of linear regression between the value of DTI parameters and age
This significant linear correlation between FA and age was not observed in children with hydrocephalus in any of the 3 ROIs in the corpus callosum (Table 3). However, FA and age still correlated highly in children with hydrocephalus in both the ALIC and PLIC. For example, FA values of patients in PLIC (Fig 3B) correlated significantly with age (R2 = 0.549; P = .002). Comparison of residuals from the linear regression found that both groups were unbiased, with a mean value of zero. No increasing or decreasing spread about the regression line was observed as the age increased for any of the ROIs.
The Frequency of Occurrence of Abnormal FA Values in Children with Hydrocephalus
In the corpus callosum, most patients had abnormally low FA values (13/16, 8/10, and 9/15, for gCC, bCC, and sCC, respectively). No child with hydrocephalus was found to have an abnormally high FA value. In the internal capsule, however, only a small portion of the group had abnormally low FA values (2/17 for both ALIC and PLIC), and some had abnormally high FA values (2/13 and 6/17 for ALIC and PLIC, respectively), which was not observed in the corpus callosum. Figure 4 demonstrates the regional difference of the frequency of abnormal FA between the corpus callosum and the internal capsule. The difference is statistically significant (P < .01, with control for multiple comparison) when comparing a ROI in the corpus callosum (gCC, bCC, or sCC) with any ROI in the internal capsule (ALIC or PLIC). On the other hand, no statistical significance was observed between any 2 ROIs within either the corpus callosum or the internal capsule.
Comparisons of frequency of occurrence of FA categorized as abnormally high, normal, and abnormally low in the 5 ROIs. Statistics were conducted on the basis of the Fisher exact test with Freeman-Halton extension.
Discussion
This is the first DTI study of WM abnormalities in children with hydrocephalus during early infancy. Previously published reports have used DTI to examine the effects of hydrocephalus on WM in older children or in young adults.14,15 In our study, the study population ranged from ages 1 day to 16 months at the time of the pretreatment MR imaging/DTI study. Thus, it fills a knowledge gap and provides a radiographic reflection of the understanding for the effects of hydrocephalus on developing WM. The sensitivity of DTI to assess WM diffusion properties in patients with hydrocephalus as demonstrated in our study potentially may help determine whether shunt surgery would be of benefit in borderline ventriculomegaly and benign external hydrocephalus or augment the fetal counseling process for parents. It may also aid in the management of shunted hydrocephalus used at present, allowing a more accurate assessment of structural integrity than the volume measurements used at present.
We included a closely age-matched healthy control group to establish a frame of reference for determining the normality of diffusion properties in patient group. The age difference between the 17 pairs of children ranged from 0 to 16 days. Such a close matching was intended to eliminate the confounding factor of age because DTI parameters have been found to change dramatically during early childhood.24,25,27 Any study of WM in pediatric patients must account for this dramatic change occurring during development, especially the first 36 months of life, to reach any conclusion about pathologic alterations.
Our results demonstrate that the DTI parameters in the corpus callosum in infants with hydrocephalus are abnormal, and this deviation from normal range is region specific. Similar to previous studies in older subjects,14 we found that the corpus callosum in the patients had lower FA and higher MD values. We also found that radial diffusivity in all of the ROIs in the corpus callosum increased significantly, which can explain the above 2 changes. In addition, the correlation of FA with age during normal development seen in our control group and in other studies24,25,27 was absent in the corpus callosum in children with hydrocephalus. Furthermore, no increasing or decreasing spread about the regression line was observed as the age increased for any of the ROIs, suggesting that WM development does not follow a normal trajectory in children with hydrocephalus in regions closest to the ventricles.
In the internal capsule, the mean FA value or other DTI parameters for children with hydrocephalus was not found to be abnormal in either ALIC or PLIC. The trend of linearly increasing of FA with age was also preserved, which may indicate a normal developmental pattern in this brain area. However, further examination demonstrated that, though FA in most patients with hydrocephalus fell within the age-appropriate normal range, there was a greater degree of variability of values in the patient group, with some patients exhibiting abnormally higher FA values and some others having abnormally low FA values. In PLIC, the abnormally high FA values are found mostly in infants older than 3 months, which is in line with the study by Assaf14 of older children with hydrocephalus. Extrapolating from the observations presented in Fig 4, we found that the frequency of occurrence of abnormalities was region specific (ie, FA in patients with hydrocephalus is more often low in the corpus callosum but high in the internal capsule). We can hypothesize that the impact of hydrocephalus is region specific and is likely related to the proximity to the enlarged ventricles. WM structures located farther away from the ventricles, such as the internal capsule, will exhibit less severe damage because the effect is dampened by the intervening compressible deep gray matter and WM structures.
No statistically significant correlation was found between FA and ventricle size. It is believed that there is a wide range of variability in the association between ventricle size and outcomes. Different underlying injury mechanisms may contribute to different FA measurements. The increase of FA as seen in the internal capsules of some patients has been sometimes suggested to be the result of the mechanical compression that leads to increased homogeneity in fiber orientation. The decrease of FA in the corpus callosum, on the other hand, is often regarded as a reflection of permanent WM damage (eg, myelin sheath and axonal cell membrane damage, or demyelination). This latter hypothesis correlates with our observation of increased radial diffusivity in the corpus callosum. However, without knowing the exact underlying mechanism of injury, it is difficult to predict whether ventricle size and FA are inversely or positively correlated.
Discrepancy between the results from our study and published literature exists in various parameters and in various ROIs. For example, increased diffusion coefficient in patients with hydrocephalus has been identified with previous DWI and DTI studies.14,28–30 Our study demonstrates a similar trend of MD change in the corpus callosum but not in the internal capsule. The variability in the FA value in the internal capsule seen in our study is also different from other results.14 These differences may be the result of variations in ROI selection. For example, the study by Uluğ et al28 defined ROIs adjacent to the ventricular horns, which does not correlate with the ROIs in our study. In addition, our cohort included extremely young patients, with rapidly changing WM water content, myelination, and axonal membrane growth in the brain. The difference one may expect during this time of WM development, compared with older children and adolescents, bears significant importance because it may demonstrate the diverse nature of the mechanism of injury seen in hydrocephalus depending on the time during development when this insult has occurred.
This study is based on retrospective analysis of existing patient data in the clinical data base. To accommodate the age range and distribution of the patient group, we were only able to find 1 closely age-matched child from our data base of healthy children. Although our one-to-one matching is no less valid than the one-to-many matching, a larger sample size may increase the efficiency of the analysis. Like most studies of patients with hydrocephalus, we do not have histopathologic correlation of the imaging findings and analysis. It would be ideal to validate the DTI conclusion with the current criterion standard of tissue examination, possibly by studying autopsy-acquired human brain tissue from patients with hydrocephalus or by studying the brain structure in various hydrocephalus animal models. We do not have long-term (≥ 5 years) behavioral and neuropsychological outcome results to relate to the imaging findings in this study. This study was also limited by the heterogeneity in the causes of the patients' hydrocephalus and DTI scanning protocols. The current trend in neurosurgery has been to divert CSF as soon as hydrocephalus is diagnosed to prevent additional injury to the CNS. However, the injury and recovery mechanisms as well as the overall prognosis in patients with hydrocephalus warrant additional investigation. A large-scale prospective longitudinal study may be able to help address the above-noted issues.
Conclusions
This study demonstrated the sensitivity of DTI techniques to investigate WM integrity in pediatric patients with hydrocephalus in infancy. We found significant alterations in diffusion values throughout the corpus callosum in infants with hydrocephalus. In the internal capsule, there was a greater degree of variability in FA values, though FA in most patients fell within the age-appropriate normal range. It is anticipated that DTI may be of value in the management of hydrocephalus, possibly helping to predict long-term outcome on the basis of pretreatment WM diffusion properties.
Footnotes
This study was supported in part by the Robert L. McLaurin, MD, Faculty Development Scholarship in Neurosurgery at Cincinnati Children's Hospital Medical Center.
References
- Received February 23, 2009.
- Accepted after revision April 5, 2009.
- Copyright © American Society of Neuroradiology