|
|
||||||||
Neuroradiology |
1 From the Department of Radiology, Landesklinikum St Poelten, Propst Fuehrer Strasse 4, A-3100 St Poelten, Austria. Received April 20, 2007; revision requested June 11; revision received July 27; final version accepted September 12. Address correspondence to E.S. (e-mail: erich.salomonowitz{at}stpoelten.lknoe.at).
| ABSTRACT |
|---|
|
|
|---|
Materials and Methods: This study was approved by the institutional review board, and informed consent was obtained. Diffusion-tensor imaging data with an isotropic voxel size of 1.9 mm3 were acquired at 3 T in 38 healthy volunteers (age range, 18–88 years; 18 women). Quantitative fiber tracking was used to calculate fractional anisotropy (FA) and mean diffusivity values, eigenvalues (
1,
2, and
3), the number of fiber projections, and the number of fiber projections per voxel for three-dimensional reconstructed association, callosal, projection, and total brain fibers. Bivariate linear regression models were used to analyze correlations. Significant differences between correlations were assessed with the Hotelling-Williams test.
Results: For FA, the strongest degradation in association fibers and no significant changes in projection fibers were observed. The difference in correlation was significant (P = .002). The number of fiber projections and the number of fiber projections per voxel showed strong to moderate negative correlations that were dependent on age (P < .001) in the three fiber structures and total brain fibers, with the exception of the number of fiber projections per voxel in projection fibers, which showed no significant correlation. The decrease in the number of fiber projections was significantly greater (P = .043) in projection fibers than in total brain fibers, whereas the decrease in the number of fiber projections per voxel was significantly weaker (P = .005). Association fibers showed the largest changes per decade of age for FA (–1.13%) and for the number of fiber projections per voxel (–4.7%), whereas callosal fibers showed the largest changes per decade of age for the number of fiber projections (–10.4%).
Conclusion: Quantitative fiber tracking enables identification of differences in diffusivity and fiber characteristics due to normal aging.
© RSNA, 2008
| INTRODUCTION |
|---|
|
|
|---|
The three major fiber structures of functional cortical regions are association, callosal, and projection fibers. Association fibers are cortex-to-cortex connections in the same hemisphere. Callosal fibers consist of an array of bundles that connect cortical areas of the hemispheres via the corpus callosum. Projection fibers connect the cortex with the thalamus, brainstem, and spinal cord (8,9). These highly organized formations of white matter affect the molecular movement of water, which is fastest when it is parallel to myelinated axonal fibers (10). Thus, the anisotropy of water diffusion reflects the microstructure of white matter tissue. Processes related to normal aging lead to changes in diffusion that result from loss of tissue organization and alterations in extracellular spaces (11,12).
Diffusion-tensor (DT) imaging enables in vivo detection of the anisotropy of water diffusion (13–16) with use of diffusion gradients in at least six directions. In general, DT data are evaluated by calculating parametric maps of fractional anisotropy (FA) and mean diffusivity (MD) that describe the directionality and magnitude of water diffusion, respectively. To determine potential regional differences, two-dimensional regions of interest (ROIs) can be defined manually within these maps. Fiber tracking enables reconstruction of white matter pathways in three dimensions with use of a tracking algorithm for comparison of orientations of water diffusion anisotropy on a voxel-by-voxel basis. Reconstructed fiber bundles may be evaluated visually or with qualitative assessment. Quantitative fiber tracking (17–20) is a technique used to determine diffusivity parameters (FA, MD, and eigenvalues) and fiber characteristics (number of fiber projections and mean number of fiber projections per voxel). The individual fiber bundles are defined by using a tracking algorithm for three-dimensional segmentation.
The purpose of our study was to prospectively quantify differences in age-related changes in diffusivity parameters and fiber characteristics between association, callosal, and projection fibers.
| MATERIALS AND METHODS |
|---|
|
|
|---|
0.5) between the diffusivity parameters (both diffusivity parameters and fiber characteristics) and age, and a power analysis (1 – β = 0.8,
= .1) revealed 36 subjects needed to be included in this study for it to have statistical significance. We examined 38 healthy volunteers (18 women, 20 men; mean age, 49.6 years ± 20.1 [standard deviation]; age range, 18–88 years). There were no significant differences in demographic variables between women and men. The mean ages of women and men were 51.6 years ± 22.1 and 47.8 years ± 18.6, respectively. The age difference between women and men was not significant (t = –0.569, df = 36, P = .573). Also, the age distribution of the participants, as assessed with the Kolmogorov-Smirnov z test, did not differ significantly (z = 0.855, P = .458). An experienced radiologist (E.S., 32 years of experience in neuroradiology) interpreted conventional MR images to exclude morphologic abnormalities. Age-related white matter changes on MR images were not an exclusion criterion; however, neurologic or psychiatric disorders, head trauma, or any history of loss of consciousness was an exclusion criterion. Participants aged 50 years or older were screened for dementia by a neurologist (T.H., 13 years of experience in neurology) using the Mini-Mental State Examination (MMSE), as is usual for this age group (21–23). A minimum MMSE score of 27 points (maximum MMSE score, 30 points) was an inclusion criterion (24) (Table 1).
|
Conventional MR imaging performed to define anatomic details consisted of (a) a coronal T2-weighted fast spin-echo sequence (7720/85, 5-mm section thickness), (b) a transverse fluid-attenuated inversion-recovery sequence (9000/133, 5-mm section thickness), (c) a sagittal T2-weighted fast spin-echo sequence (2800/83, 2-mm section thickness), and (d) a three-dimensional T1-weighted magnetization-prepared rapid acquisition gradient-echo sequence (1550/2.6, 0.47 x 0.47 x 0.94-mm voxel size, 176 sections).
We used T1-weighted magnetization-prepared rapid acquisition gradient-echo data sets and a region-growing process of the BrainVoyager QX software (version 1.7; Brain Innovation, Maastricht, the Netherlands) for volumetry of white matter (25). Talairach normalization was performed to transform subcortical brain structures to standardized spatial positions. This approach provided a priori positional knowledge, which was needed to segment white and gray brain matter from the ventricles and the surrounding head tissue, including the cerebellum.
To determine a threshold to automatically separate white matter from gray matter, a series of 10 signal intensity histograms across several transverse sections were analyzed. The curves showed two major peaks: The first peak corresponded to the intensity of gray matter, whereas the second peak corresponded to the intensity of white matter. The aforementioned software program was used to compute the threshold value that separated white matter from gray matter. This value was used for the segmentation procedure, which was used to perform a region-growing process to label white matter voxels.
Additionally, the region-growing process needed a seed voxel to start the growing process. The seed voxel was determined with an imaging process that started in the center of the brain and continued until a voxel that had an intensity higher than that of the determined white matter–gray matter threshold was found. From the determined seed voxel, neighboring voxels (ie, the six voxels in the direct neighborhood) are considered for inclusion in white matter segmentation. A neighboring voxel is added to the white matter segmentation only if its intensity is higher than the intensity threshold. The neighboring voxels of all included voxels were checked for inclusion with use of the same threshold. This process continued until no more neighbor voxels with intensities that exceeded the intensity threshold were found. After completion of the region-growing process, the gray matter voxels and remaining head tissue were set to zero; this resulted in a representation containing only the labeled white matter voxels. The labeled voxels were counted, and the white matter volume was calculated.
DT Data Processing and Quantitative Fiber Tracking
DT data were transferred offline to a workstation (Inspiron 8200; Dell, Round Rock, Tex) for analysis. One author (A.S., 6 years of experience in brain MR imaging) using DtiStudio software (version 2.4; H. Jiang, S. Mori, Department of Radiology, Johns Hopkins University, Baltimore, Md) (26) performed data processing. The DT was calculated for each voxel with a multivariate linear fitting method (27). Diagonalization of the DT was used to calculate three eigenvectors and three eigenvalues (
1,
2,
3) that corresponded to the main diffusion directions and associated diffusivities. The concept of diffusion ellipsoids enables three-dimensional representation of the diffusion distance in a given diffusion time. The eigenvalues are the axes of these ellipsoids and can be calculated for each image voxel (10,28,29). Parametric maps of FA and MD were calculated (30).
The three major white matter tracts that were investigated in this study were defined according to the classification of functional categories described by Wakana et al (27), as follows: (a) association fibers, including the superior longitudinal, inferior longitudinal, and inferior fronto-occipital fasciculi; (b) callosal fibers; and (c) projection fibers—including the corticobulbar and corticospinal tracts—and thalamic fibers—including anterior, superior, and posterior thalamic radiations.
The Fiber Assignment by Continuous Tracking, or FACT, algorithm was used for three-dimensional reconstruction of fiber tracts. With this algorithm, one starts tracking by tracking every voxel, and all voxels of the image volume are encompassed (16,31,32) so that each voxel has multiple fibers passing through it. The fibers of interest were selected by designating manually defined multiple ROIs and using the three logical operators AND, OR, and NOT (26). ROIs were determined by a radiologist (E.S.) and neurosurgeron (O.G., 13 years of experience in brain MR imaging) working in consensus and by using T2-weighted MR images or DT imaging–based color-coded maps (for the superior longitudinal fasciculus) according to the guidelines given by Wakana et al (27) and Mori et al (8). The tracking procedure was stopped when a track-turning angle greater than 60° was encountered. An FA threshold of 0.2 was used as the indicator for termination of tract elongation.
For visualization of the superior longitudinal fasciculus, one ROI was defined on a coronal section of the DT imaging–based color-coded maps at the posterior tip of the putamen. One ROI at the parieto-occipital sulcus, which was identified at the middle of the coronal section along the superior-inferior axis, was used for both the inferior longitudinal fasciculus and the inferior fronto-occipital fasciculus. For the inferior longitudinal fasciculus, an additional ROI was defined on a coronal section at the midtemporal lobe at the section level of the posterior tip of the putamen. For the inferior fronto-occipital fasciculus, an additional ROI was defined on a coronal section at the frontal lobe at the section level where the frontal and temporal lobes were separated. The ROIs for the inferior longitudinal fasciculus and the inferior fronto-occipital fasciculus were combined by using the operator AND (8).
The callosal fibers were reconstructed in three sections. The first ROI was placed in the corpus callosum at the midsagittal level. Two additional ROIs were combined with the operator AND and were placed on coronal sections in the corona radiata. One was placed anterior to the genu, and one was placed posterior to the splenium of the corpus callosum. One ROI was placed on transverse sections superior to the body of the corpus callosum. The three reconstructed fiber bundles were combined with the operator OR. The operator NOT was used to remove thalamic fibers, which were reconstructed occasionally (27).
The corticobulbar and corticospinal tracts were reconstructed by placing one ROI at the cerebral peduncle and one ROI at the internal capsule. Both ROIs were combined by using the operator AND. For reconstruction of the thalamic fibers, the entire thalamus was defined as the first ROI. The second ROI for the anterior and posterior thalamic radiations was placed on coronal sections to define the frontal lobe at the section level where the frontal and temporal lobes were separated and to define the occipital lobe at the section level of the posterior tip of the putamen, respectively. For superior thalamic radiation, the second ROI was defined on a transverse section above the corpus callosum, and it occupied the entire hemisphere. Again, the first and the second ROIs were combined by using the operator AND (8,27).
Quantitative fiber tracking of the reconstructed fibers was performed by using the tract statistics function of the DtiStudio software, which enables statistical evaluation of pixels occupied by reconstructed fibers. Diffusivity parameters, as well as fiber characteristics, were calculated for the selected fiber bundles. For total brain fibers, the diffusivity parameters (FA, MD, and eigenvalues) cannot be calculated with the tract statistics function of the DtiStudio software.
The number of fiber projections is the number of streamlines that were reconstructed and that penetrate the ROI. Hence, this parameter includes not only the fibers started from the points in the ROI but also the streamlines started from other points that pass through the ROI. Both the number of streamlines within each voxel of the reconstructed fiber structure and the number of voxels occupied by the fiber structure (measured in cubic millimeters) were counted and used to calculate the mean number of fiber projections per voxel. Therefore, the number of fiber projections per voxel was used to measure the fiber density of the selected fiber bundle. The number of fiber projections and the number of fiber projections per voxel are limited by the spatial resolution achieved and the technique applied in this study. Both parameters have to be interpreted as relative parameters and are not measures of the true number of axons passing through a voxel.
Statistical Methods
Data were analyzed by two authors (A.S., G.S.; each with 15 years of experience in statistics) by using statistical software (SPSS 12; SPSS, Chicago, Ill). A two-sided paired Student t test was used to compare differences between hemispheres. A Kolmogorov-Smirnov z test was used to compare age distributions between women and men. Bivariate linear regression models were used to analyze the influence of age on diffusivity parameters and fiber characteristics. The significance of differences between correlations was tested with the Hotelling-Williams test (33). Three authors (A.S., G.S., and O.G.) interpreted the correlation coefficients with consideration of the specifications for interpretation of correlation coefficients given by Zou et al (34). For all tests, the level of significance was set at P < .05.
To assess the influence of age and sex on diffusivity parameters and fiber characteristics, multiple linear regression models were used. Multiple linear correlations with sex as the control variable revealed that sex did not have a significant influence, and results were qualitatively similar to those obtained with the bivariate linear regression model. Thus, the authors thought it was justifiable to present the results for only the bivariate model.
| RESULTS |
|---|
|
|
|---|
1, P = .128;
2, P = .446;
3, P = .058; number of fiber projections, P = .349; and number of fiber projections per voxel, P = .073) and projection fibers (FA, P = .073; MD, P = .074;
1, P = .128;
2, P = .446;
3, P = .058; number of fiber projections, P = .349; and number of fiber projections per voxel, P = .073). For further statistical evaluation, we used (a) mean values from the left and right hemispheres for FA, MD, eigenvalues, and number of fiber projections per voxel and (b) the sum from the left and right hemispheres for the number of fiber projections.
|
|
1 of the three fiber structures showed moderately positive correlations that were dependent on age (association fibers: r = 0.349, P = .032; callosal fibers: r = 0.377, P = .02; and projection fibers: r = 0.39, P = .016). For
2, we found weakly positive correlations that were significant for association fibers (r = 0.315, P = .044) but not for callosal (r = 0.239, P = .148) or projection (r = 0.251, P = .128) fibers. The
3 showed moderately significant positive correlations that depended on age for all three fiber structures (association fibers: r = 0.317, P = .042; callosal fibers: r = 0.363, P = .025; and projection fibers: r = 0.378, P = .019). There were no significant differences in correlations of eigenvalues between the fiber structures. For
1, P values were as follows: .846 for association versus callosal fibers, .8 for association versus projection fibers, and .927 for callosal versus projection fibers. For
2, P values were as follows: .88 for association versus callosal fibers, .811 for association versus projection fibers, and .946 for callosal versus projection fibers. For
3, P values were as follows: .887 for association versus callosal fibers, .792 for association versus projection fibers, and .937 for callosal versus projection fibers.
Age-related Changes
Age-related changes in the number of fiber projections for association, callosal, projection, and total brain fibers are shown in Figure 3. For total brain fibers and all three fiber structures, a bivariate linear regression model showed there was a strong negative correlation between the number of fiber projections and age (P < .001). Only the difference between the correlations of projection fibers (r = –0.834) and total brain fibers (r = –0.688) was significant (P = .043) at Hotelling-Williams testing (ie, the number of fiber projections of the projection fibers decreased significantly more with age than did the number of fiber projections of the total brain fibers).
|
|
|
Hotelling-Williams testing revealed significant differences between correlations for the number of fiber projections for fiber structures and the white matter volume (P = .01 for association fibers, P = .039 for callosal fibers, P < .001 for projection fibers, and P < .001 for total brain fibers). However, the differences between correlations for the number of fiber projections per voxel for the fiber structures and white matter volume were not significant (P = .227 for association fibers, P = .792 for callosal fibers, P = .27 for projection fibers, and P = .192 for total brain fibers).
For diffusivity parameters, only FA in projection fibers had a significantly weaker correlation compared with the correlation of white matter volume (P = .027). All other correlations for diffusivity parameters were not significantly different from the correlation of white matter volume with age. For association fibers, P values were as follows: .303 for FA, .525 for MD, .4 for
1, .155 for
2, and .171 for
3. For callosal fibers, P values were as follows: .874 for FA, .975 for MD, .505 for
1, .196 for
2, and .439 for
3. For projection fibers, P values were as follows: .853 for MD, .545 for
1, .224 for
2, and .38 for
3.
| DISCUSSION |
|---|
|
|
|---|
Our results for the age-related changes in fiber characteristics in projection fibers—a large significant decrease in the number of fiber projections and a small insignificant decrease in the number of fiber projections per voxel—may be explained by the fact that an age-related decrease in the number of fiber projections may be associated with a decrease in the volume of the fiber structure that resulted in a smaller reduction in fiber density (the number of fiber projections per voxel) and FA. On the other hand, it appears that the volume of association fiber structures decreases the least compared with the volume of other white matter structures investigated in this study. The largest decrease was registered for the number of fiber projections per voxel. This finding is supported by the fact that the decrease in FA was largest in association fibers. This decrease was related to the largest and only significant increase in the perpendicular component of diffusion
2, and it may have been due to increased extracellular space (11,12).
Age-related changes in callosal fibers occupied an intermediate position: The most pronounced decrease in the number of fiber projections was associated with a large decrease in both the number of fiber projections per voxel and the FA, indicating a similar decrease in callosal volume. However, in our estimation, the most important parameter was the number of fiber projections, which was most affected by age-related degradation of fiber structures, including frontal white matter (callosal and projection fibers). Furthermore, there was a significantly larger decrease in the number of fiber projections in the three fiber structures and the total brain fibers compared with the decrease in total white matter volume, indicating additional age-related changes beyond the well-known brain atrophies.
In a number of studies, researchers used DT imaging and manually defined ROIs to determine age-related changes in the diffusivity of water in different regions of white matter. Study results have shown a significant age-related decrease in tissue volume in frontal white matter (35), the genu of the corpus callosum (36–38), and the posterior limb of the internal capsule (39,40). These results are in agreement with our findings, as we registered a significant correlation between FA and age in callosal fibers. The use of manually defined two-dimensional ROIs in the evaluation of DT data limits the detection of a three-dimensional process, such as that which occurs with aging in white matter structures of the brain. Advanced and interesting approaches for three-dimensional evaluation of age-related changes in DT data have been presented in two relatively recent studies (12,21).
Nusbaum et al (12) co-registered DT data to T1-weighted MR images and standardized maps of relative anisotropy to normal coordinates for statistical probability mapping. Pixels with which there was a significant correlation between relative anisotropy and age were calculated with a linear regression model. Significantly decreased anisotropy with increasing age was found in the frontal white matter, the genu and splenium of the corpus callosum, and the periventricular white matter. Significant increases in relative anisotropy were found in the internal capsules on both sides. Nusbaum et al (12) concluded that their technique may permit global assessment of changes in the organization of white matter pathways that occur with normal aging.
Pfefferbaum et al (21) presented an anteroposterior profile analysis of the integrity of the white matter microstructure across supratentorial regions by using DT data collected at 3 T. The DT imaging profile analysis revealed that the frontal distribution of low white matter FA was more robust in healthy older adults than in healthy younger adults. Sullivan et al (23) described the use of a quantitative fiber-tracking approach to examine age-related degradation in diffusivity parameters (FA and apparent diffusion coefficient) and the number of fiber projections in six regions of the corpus callosum. In the frontal callosal fibers particularly, older participants had lower FA values, higher apparent diffusion coefficient values (apparent diffusion coefficient is comparable to MD), and fewer fibers than did younger participants. Sullivan et al (23) concluded that the quantitative fiber-tracking approach enables one to confirm the selective vulnerability of frontal callosal fibers to normal aging as the mechanism underlying age-related decreases in cognition. The results of a study by Ota et al (41), in which the authors used quantitative fiber tracking, confirmed these results for age-related FA and MD changes in frontal and callosal fibers. However, unlike Sullivan et al (23) and us, Ota et al (41) did not evaluate fiber characteristics. The evaluation of fiber characteristics yields detailed information about the mechanisms of age-related changes in the fiber structures of white matter.
Our study was limited by the rather small number of participants examined. However, the number of participants in each age group was more or less equal. A further limitation might have been the fact that we did not include other fiber structures and did not investigate differences between the substructures of association, callosal, and projection fibers. These issues should be addressed in future studies (eg, for association fibers, projection fibers, or short-range U fibers in different lobes), with a view toward identifying potential differences between the substructures.
In conclusion, we found that quantitative evaluation of fiber-tracking results enabled us to identify differences in age-related changes in diffusivity parameters and fiber characteristics between different fiber structures. The degradation in the number of fiber projections was stronger compared with the decrease in total white matter volume, suggesting that a process more complex than the loss of white matter volume occurs in the normal aging brain.
| ADVANCE IN KNOWLEDGE |
|---|
|
|
|---|
| IMPLICATION FOR PATIENT CARE |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Abbreviations: DT = diffusion tensor FA = fractional anisotropy MD = mean diffusivity ROI = region of interest
Author contributions: Guarantors of integrity of entire study, A.S., E.S., T.H., O.G.; study concepts/study design or data acquisition or data analysis/interpretation, all authors; manuscript drafting or manuscript revision for important intellectual content, all authors; manuscript final version approval, all authors; literature research, A.S., E.S., G.S., O.G.; clinical studies, A.S., E.S., T.H.; statistical analysis, A.S., E.S., G.S., O.G.; and manuscript editing, all authors
Authors stated no financial relationship to disclose.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. A. Chaddock, G. J. Barker, N. Marshall, K. Schulze, M. H. Hall, A. Fern, M. Walshe, E. Bramon, X. A. Chitnis, R. Murray, et al. White matter microstructural impairments and genetic liability to familial bipolar I disorder The British Journal of Psychiatry, June 1, 2009; 194(6): 527 - 534. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |