|
|
||||||||
Neuroradiology |
1 From the Department of Neurosurgery, Neurocenter (A.S., O.G., M.B., C.N.), Departments of Neuropathology (R.B.) and Neurology (T.H.), University of Erlangen-Nuremberg, Erlangen, Germany; Department of Radiology, Landesklinikum St. Poelten, Propst Fuehrer Strasse 4, A-3100 St. Poelten, Austria (A.S., E.S.); and MR Centre of Excellence, Medical University of Vienna, Vienna, Austria (S.G., E.M.). Received June 4, 2005; revision requested August 1; revision received September 2; accepted September 22; final version accepted November 17. A.S. supported by the German Research Society (SFB 603, C9). Address correspondence to E.S. (e-mail: erich.salomonowitz{at}stpoelten.lknoe.at).
| ABSTRACT |
|---|
|
|
|---|
Materials and Methods: The institutional review board required neither ethics committee approval nor patient informed consent for this study. Twenty patients (eight women, 12 men; age range, 1853 years) with glioma (seven World Health Organization grade II and 13 grade III tumors) underwent diffusion-tensor MR imaging at 1.5 T. Diffusion-tensor data were obtained with an echo-planar imaging sequence with six diffusion directions (b = 1000 sec/mm2), isotropic 1.9-mm voxels, and five averages. FA and mean diffusivity values were calculated from diffusion-tensor data. Coregistration with a three-dimensional MR imaging data set (used for stereotactic brain biopsies) enabled correlation of FA and mean diffusivity values with the histopathologic findings total cell number (CN), tumor CN, and percentage tumor infiltration (TI) by using linear, exponential, and logarithmic models. Student t and Mann-Whitney U tests were performed.
Results: Histopathologic findings of 77 MR imageguided stereotactic biopsies in all 20 patients were correlated with FA and mean diffusivity values at the biopsy locus. For FA and mean diffusivity, a logarithmic model showed strongest correlation with tumor CN and total CN; a linear model showed strongest correlation with percentage TI. For FA there were negative logarithmic (R = 0.802, P < .001) and linear (R = 0.796, P < .001) correlations with tumor CN and percentage TI, respectively. For mean diffusivity there were positive logarithmic (R = 0.557, P < .001) and linear correlations (R = 0.521, P < .001) with tumor CN and percentage TI, respectively. Differences between correlations for FA and mean diffusivity versus tumor CN (P < .001) and percentage TI (P < .001) were significant.
Conclusion: FA is better than mean diffusivity for assessment and delineation of different degrees of pathologic changes (ie, TI) in glioma.
© RSNA, 2006
| INTRODUCTION |
|---|
|
|
|---|
Pathologic processes (eg, tumors, stroke, multiple sclerosis) result in changes in diffusion because of several reasons, such as loss of tissue organization or changes in extracellular space. Cerebral gliomas cause not only disruption or displacement of white matter structures but also widening of fiber bundles due to tumor infiltration (TI) or edema (310). Information on white matter involvement and glioma microstructure is essential for neurosurgical planning with the aim of maximum tumor resection and minimum damage to healthy brain. Routine MR imaging methods do not enable these questions to be answered adequately for the care of patients with glioma.
Diffusion-tensor MR imaging enables the macroscopic detection of anisotropic water diffusion due to elongated structures such as white matter bundles (11,12). Diffusion-tensor MR imaging measures differences in the diffusion displacement of water in three dimensionsdepending on image positionby employing at least six diffusion gradient directions. The diffusion tensor D for each image voxel is estimated from this series of experiments, which include an MR imaging acquisition performed without diffusion weighting (b = 0 sec/mm2). The magnitude and directionality of water diffusionreferred to as mean diffusivity and fractional anisotropy (FA), respectivelycan be calculated by using D. Mean diffusivity is mathematically equivalent to the well-known apparent diffusion coefficient obtained with standard diffusion-weighted imaging methods (2,13).
Several researchers (5,7,10,13,14) have recommended histologic or clinical evaluation of the results of diffusion-tensor MR imaging methods used for the examination of patients with tumors. The purpose of our study was to retrospectively correlate the changes in FA and mean diffusivity in gliomas at diffusion-tensor MR imaging with the degree of tumor cell infiltration determined histologically.
| MATERIALS AND METHODS |
|---|
|
|
|---|
MR Imaging Methods
MR imaging was performed with a 1.5-T clinical whole-body unit (Magnetom Sonata; Siemens, Erlangen, Germany) equipped with the standard head coil. For diffusion-tensor imaging, a diffusion-weighted echo-planar imaging sequence was used with the following parameters: repetition time msec/echo time msec, 9200/86; matrix size, 128 x 128; field of view, 240 x 240 mm; section thickness, 1.9 mm; and bandwidth, 1502 Hz/pixel. Sixty sections were measured with no intersection gap and an isotropic voxel size of 1.9 x 1.9 x 1.9 mm3. Diffusion-gradient encoding in six directions with a b value of 1000 sec/mm2 and an additional measurement without diffusion-gradient encoding (b = 0 sec/mm2) were performed.
The sequence design was based on the use of balanced diffusion gradients to minimize eddy current artifacts. For sufficient signal-to-noise ratios, five averages were applied, requiring a total diffusion-tensor imaging data acquisition time of 5 minutes 31 seconds. Immobilization of the patient's head to minimize artifacts due to patient motion was achieved by fixation of the head in a headrest.
Five adhesive skin fiducial marks were placed in a scattered pattern on the head surface for registration to the frameless stereotactic system. To obtain an MR data set for neuronavigation, a 3D anatomic magnetization-prepared rapid-acquisition gradient-echo sequence was performed in a single session 1 day before surgery with the following parameters: 2020/4.38; field of view, 25 x 25 cm; 1-mm isotropic voxels; and 160 sections (15).
Data Processing
Diffusion-tensor imaging data were transferred off-line to a workstation (Inspiron 8200; Dell, Round Rock, Tex) for analysis. Data processing was performed with software (DTI-Studio, version 2.3; H. Jiang, S. Mori; Johns Hopkins University, Baltimore, Md). The diffusion tensor (D) was calculated by one author (A.S., with 4 years of experience with brain MR imaging) for each voxel by using a multivariate linear fitting method (16,17). Diagonalization of the tensor yielded three eigenvalues (
1,
2, and
3) and three eigenvectors (17,18). Parametric maps of mean diffusivity (MD) and FA (19) were calculated by using the following equations:
|
|
|
Mean diffusivity has the unit square millimeters per second, while FA values are unitless and vary between 0 (isotropic diffusion) and 1 (complete anisotropy). To compare means of mean diffusivity and FA in tumor and contralateral normal-appearing white matter (NAWM), regions of interest were drawn to cover approximately 90% of the lesion and an area equal in size on the contralateral side. Region-of-interest drawing was performed in consensus by a neurosurgeon (O.G., with 12 years of experience in brain tumor imaging) and a radiologist (E.S., with 31 years of neuroradiology experience). The average mean diffusivity and FA values for the lesion and the contralateral NAWM were calculated by averaging them over the regions of interest.
MR Imageguided Stereotactic Biopsy
Before tumor resection, biopsies were performed by one of two neurosurgeons (O.G. or C.N., with 13 and 15 years of neurosurgery experience, respectively) in several regions according to the information from the 3D magnetization-prepared rapid-acquisition gradient-echo examination by using a stereotactic needle tracked by a navigation system (VectorVision Sky; Brain LAB, Heimstetten, Germany). This procedure was performed before the actual tumor resection and thus ensured a minimal interference of brain shift, which would render neuronavigation inaccurate. The coordinates of each biopsy locus were labeled and documented within the 3D MR imaging data set by one of two neurosurgeons (O.G. or C.N.). No perioperative adverse effects were observed during or after biopsy sampling.
Diffusion-tensor imaging data sets were transferred to the planning workstation of the navigation system for retrospective registration to the 3D MR imaging data sets by using image fusion software (VV2 Planning, version 1.3; BrainLAB). This automated registration process uses a rigid registration algorithm that applies an intensity-based pyramidal approach by using mutual information (20). Histopathologic findings for each specimen were traced back to the exact voxel positions in the diffusion-tensor imaging data set (by A.S.). To increase statistical accuracy and to account for possible minor matching inaccuracies, the averages of the mean diffusivity and FA values were calculated from a 3 x 3 x 3-voxel cube around the determined voxel position. The volume of the voxel cube thus evaluated was 5.7 x 5.7 x 5.7 mm3.
Histopathologic Evaluation
All glioma specimens were histologically examined by a neuropathologist (R.B.). The tumors were classified according to the WHO guidelines for tumors of the nervous system by using hematoxylin-eosin staining, as well as a standardized panel of immunohistochemical markers (MAP2, GFAP, p53, and Ki67).
The following antibodies were used: glial fibrillary acidic protein (GFAP clone 6F2; DAKO, Glostrup, Denmark), Ki67 epitope (Mib1; Dianova, Hamburg, Germany), and p53 (DAKO). For detection of the microtubule-associated protein MAP2, we used the monoclonal antibody clone C (supplied by Beat Riederer, MD, Lausanne, Switzerland). Commercially available antibodies reacting with either high- or low-molecular-weight isoforms of MAP2 were used (clone HM-210, Sigma-Aldrich, Munich, Germany; and MAB364, Chemicon, Hampshire, England). The specificity of MAP2 in human brain tumor samples has been extensively studied with in situ hybridization, Western blotting experiments, and immunohistochemical analysis (21,22).
Quantitative assessment of glial tumor cells and preexisting brain parenchyma was microscopically performed by using analysis software (Soft Imaging System, Leinenfeld-Echterdingen, Germany) at x200 magnification averaged over five 348 x 261-µm subfields. Only MAP2 immunoreactive cells with a distinct nucleus were taken into account. The averaged number of glial tumor cells was defined as mean tumor cell number (CN). The sum of the averaged number of glial tumor cells and preexisting brain parenchymal cells was defined as mean total CN. The ratio of mean tumor CN to mean total CN was defined as percentage TI.
Statistical Analysis
Data were analyzed by one author (A.S.) with statistical software (SPSS, version 12; SPSS, Chicago, Ill). For comparisons of mean FA and mean diffusivity values averaged over tumor regions of interest and contralateral control regions, a two-sided paired Student t test, with significant differences indicated by P < .05, was used for statistical analysis. The differences in mean FA and average mean diffusivity between tumor and contralateral NAWM for each patient were calculated. The means and the 95% confidence intervals (CIs) for these differences were determined. Standard error estimates for sample means were denoted by the standard error of the mean (SEM). Comparisons of diffusion-tensor imaging measurements between groups (patients with WHO grade II glioma and those with WHO grade III glioma) were performed by using the Mann-Whitney U test.
Linear, exponential, and logarithmic regression analyses were performed for correlations of diffusion-tensor imaging metrics (FA and mean diffusivity) with the three histologic parameters (mean tumor CN, mean total CN, and percentage TI). The significance of differences between correlations for FA and mean diffusivity were tested by using a t test (23). Linear regressions were performed to determine intercept points, ranges of intersection determined as the CI, and SEMs at a percentage TI of 0 for FA and mean diffusivity.
Interpretations of the correlation coefficients were made by A.S., O.G., and E.S., with consideration of the specifications for interpretation of correlation coefficients given by Zou et al (24).
| RESULTS |
|---|
|
|
|---|
The average mean diffusivity for WHO grade II and III gliomas was 1.557 x 103 mm2/sec (standard deviation, 0.353; SEM, 0.079) compared with an average mean diffusivity for the equal-sized region of interest in contralateral NAWM of 0.826 x 103 mm2/sec (standard deviation, 0.089; SEM, 0.020). The mean difference between mean diffusivity values in gliomas and contralateral NAWM was 0.731 x 103 mm2/sec (95% CI: 0.574, 0.888). These differences were significant (P < .001, t test).
Grade II versus Grade III Glioma
The mean FA for the subgroup of patients with a WHO grade II glioma was 0.196 (standard deviation, 0.032; SEM, 0.012), while that for the patients with a WHO grade III glioma was 0.170 (standard deviation, 0.051; SEM, 0.014). The average mean diffusivity for WHO grade II glioma was 1.504 x 103 mm2/sec (standard deviation, 0.242; SEM, 0.091), and that for WHO grade III glioma was 1.585 x 103 mm2/sec (standard deviation, 0.407; SEM, 0.113). Comparisons of FA and mean diffusivity values between patients with WHO grade II glioma and patients with WHO grade III glioma with the Mann-Whitney U test revealed significant differences only for FA in tumor (P = .046). All other comparisons revealed no significant differences.
Histopathologic Correlation with Stereotactic Biopsy Specimens
In this study, 77 stereotactic biopsy specimens were evaluated histopathologically and correlated with diffusion-tensor imaging metrics (FA and mean diffusivity) by means of coregistration of patient anatomy and diffusion-tensor imaging data by using image fusion by semiautomatic rigid registration (25) (Fig 1).
|
|
|
|
|
|
|
|
|
|
|
|
Regression analyses of FA and mean diffusivity revealed strongly negative and moderately positive linear correlations versus percentage TI (Fig 4), respectively. Intercept values at a percentage TI of 0% were obtained at an FA value of 0.278 (95% CI: 0.258, 0.294; SEM, 0.034) and a mean diffusivity value of 1.172 x 103 mm2/sec (95% CI: 1.021, 1.322; SEM, 0.287), respectively, which might be interpreted as ranges of threshold values for brain tissue that was not infiltrated by tumor and/or peritumoral edema.
|
|
| DISCUSSION |
|---|
|
|
|---|
Mean diffusivity values obtained in our study are in agreement with those in earlier published studies (7,13,14,26) for both tumor area and contralateral NAWM. However, regarding FA values, the results are not unequivocal. Our FA values are in agreement with data published by Provenzale et al (13) and Lu et al (14). Differing values were communicated by Tropine et al (26) (0.123 for tumor; 0.504 for contralateral NAWM) and Sinha et al (7) (0.11 for tumor core; 0.48 for contralateral NAWM). This difference may be attributed to the higher number of glioblastomas in their patient groups.
Our results regarding subgroups of WHO grade II and III gliomas showed only a weakly significant (P = .046) difference for FA in tumor. Owing to the small number of patients with WHO grade II glioma in our study (ie, eight), this result is not representative. Hence, we conclude that diffusion-tensor imaging may not be useful for the preoperative differentiation of WHO grade II and III gliomas. This supports the findings of Provenzale et al regarding low-grade glioma and glioblastoma.
Although numerous publications exist in which the role of diffusion-tensor imaging was investigated for characterization of peritumoral infiltration area and solid tumor, an often suggested investigation of a correlation between diffusion-tensor imaging metrics and histopathologic findings (5,7,10,13,14) has not yet been sufficiently addressed. In two publications by Beppu and colleagues (3,27), an attempt was made to predict cell density and proliferation in glioblastoma through correlation with computed tomographyguided stereotactic samples. Their finding of a linear increase in FA with increasing tumor CN is in contrast to the commonly obtained finding of a decrease in FA from normal to pathologic tissue. The authors did not describe how the correlation of the stereotactic specimens to the diffusion-tensor imaging metrics was established. A reason for the different results may be that the regions of interest chosen for the measurement of the FA values that corresponded to the biopsy site were rather large in the studies by Beppu et al. Also, they investigated a rather heterogeneous patient population in one study and only glioblastoma in the other study while performing only one biopsy per case, which may not yield sufficient information on the pathologic situation.
In our study, the correlation of FA and mean diffusivity versus total CN and tumor CN, respectively, showed a logarithmic relationship between histologic parameters and diffusion properties. For small tumor CNs, the logarithmic regression lines demonstrate the typical coursenamely, a strong decrease in FA and an increase in mean diffusivity for low tumor CNs followed by an asymptotic plateauing trend for high tumor CNs.
We hypothesize that these findings reflect the infiltrative invasion mechanism of gliomas: enlargement of the extracellular space but widely preserved normal microstructure of the fiber bundles in the border zone. Increased TI causes a decrease in the extracellular space (eg, edema) and a decrease in directionality of diffusion because of derangement of the microstructures.
The key finding of our study is that the comparisons of mean diffusivity and FA in correlation with histopathologic findings demonstrate that the FA values render a better correlation with all histopathologic parameters (mean tumor CN, mean total CN, and percentage TI). On the other hand, the ranges of interception for FA and mean diffusivity values at a percentage TI of 0% might lead to the interpretation that these values could serve as a threshold for the diagnosis of nontumor-infiltrated brain tissue and/or peritumoral edema.
Limitations of this study were related to the spatial distortions inherent in the echo-planar imaging sequence, so registration with the standard anatomic image data used for navigation, which is necessary to locate the biopsy position in the diffusion-tensor imaging data, may have caused some inaccuracy. Furthermore, the patient registration error of the navigation system, which is typically on the order of 12 mm, also influences this correct spatial correlation. We tried to consider these errors by averaging the FA and mean diffusivity values of a cube measuring 3 x 3 x 3 voxels around the biopsy position.
In conclusion, FA is better than mean diffusivity for assessment and delineation of different degrees of pathologic changes (ie, TI) in glioma.
| ADVANCES IN KNOWLEDGE |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Abbreviations: CI = confidence interval CN = cell number FA = fractional anisotropy NAWM = normal-appearing white matter SEM = standard error of the mean 3D = three-dimensional TI = tumor infiltration WHO = World Health Organization
Authors stated no financial relationship to disclose.
Author contributions: Guarantors of integrity of entire study, A.S., O.G., E.S., C.N.; study concepts/study design or data acquisition or data analysis/interpretation, all authors; manuscript drafting or manuscript revision for important intellectual content, all authors; approval of final version of submitted manuscript, all authors; literature research, O.G., R.B., S.G., E.M., E.S., C.N.; clinical studies, O.G., C.N.; experimental studies, A.S., O.G., R.B., C.N.; statistical analysis, A.S., T.H., E.S.; and manuscript editing, all authors
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L.S. Hu, L.C. Baxter, K.A. Smith, B.G. Feuerstein, J.P. Karis, J.M. Eschbacher, S.W. Coons, P. Nakaji, R.F. Yeh, J. Debbins, et al. Relative Cerebral Blood Volume Values to Differentiate High-Grade Glioma Recurrence from Posttreatment Radiation Effect: Direct Correlation between Image-Guided Tissue Histopathology and Localized Dynamic Susceptibility-Weighted Contrast-Enhanced Perfusion MR Imaging Measurements AJNR Am. J. Neuroradiol., March 1, 2009; 30(3): 552 - 558. [Abstract] [Full Text] [PDF] |
||||
![]() |
C.-H. Toh, M. Castillo, A.M.-C. Wong, K.-C. Wei, H.-F. Wong, S.-H. Ng, and Y.-L. Wan Primary Cerebral Lymphoma and Glioblastoma Multiforme: Differences in Diffusion Characteristics Evaluated with Diffusion Tensor Imaging AJNR Am. J. Neuroradiol., March 1, 2008; 29(3): 471 - 475. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Sadeghi, N. D'Haene, C. Decaestecker, M. Levivier, T. Metens, C. Maris, D. Wikler, D. Baleriaux, I. Salmon, and S. Goldman Apparent Diffusion Coefficient and Cerebral Blood Volume in Brain Gliomas: Relation to Tumor Cell Density and Tumor Microvessel Density Based on Stereotactic Biopsies AJNR Am. J. Neuroradiol., March 1, 2008; 29(3): 476 - 482. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Stadlbauer, C. Nimsky, S. Gruber, E. Moser, T. Hammen, T. Engelhorn, M. Buchfelder, and O. Ganslandt Changes in Fiber Integrity, Diffusivity, and Metabolism of the Pyramidal Tract Adjacent to Gliomas: A Quantitative Diffusion Tensor Fiber Tracking and MR Spectroscopic Imaging Study AJNR Am. J. Neuroradiol., March 1, 2007; 28(3): 462 - 469. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |