| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Full paper |
1 Division of Imaging Science and Biomedical Engineering, Faculty of Medical and Human Sciences, The University of Manchester, Stopford Building, Oxford Road, Manchester M13 9PT, UK, 2 Department of Radiology, School of Medicine, University of California at San Francisco, San Francisco, California 94143, USA
Correspondence: Professor Alan Jackson, Division of Imaging Science and Biomedical Engineering, Faculty of Medical and Human Sciences, The University of Manchester, Stopford Building, Oxford Road, Manchester M13 9PT, UK. E-mail: alan.jackson{at}man.ac.uk
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
The most commonly applied method for mapping CBV is based on dynamic contrast-enhanced imaging using susceptibility (T2*) weighted images (DSCE-MRI) [18, 23, 24]. This approach was first used to study the perfusion of normal cerebral capillary beds in grey and white matter where the bloodbrain barrier (BBB) is intact and gadolinium based contrast agents can be assumed to act as purely intravascular markers. In these applications susceptibility based contrast mechanisms have the advantage over other MR techniques of increased signal-to-noise ratio in areas of low blood volume such as capillary beds [2426]. However, susceptibility based techniques have significant disadvantages related to spatial distortion occurring in areas of paramagnetic variation and residual relaxivity effects in areas of extravascular contrast leakage [27]. This has led several workers to try to derive CBV maps from T1 weighted data [2832]. Although these methods showed some promise they are not routinely used, largely because of their tendency to give rise to erroneously high estimates of CBV in areas of extravascular contrast leakage.
In 2000, Li et al [33] described a novel method for the analysis of relaxivity (T1) based dynamic contrast enhanced MRI (DRCE-MRI) which uses a shape-based analysis to decompose the dynamic contrast concentration time course data into intravascular and extravascular components. One consequence of this analytical approach is the generation of an estimate of the intravascular contrast concentration time course curve that is free of leakage effects and therefore ideal for the calculation of leakage-free estimates of CBV. Examination of the method using Monte Carlo simulation techniques has shown that the representation of CBV values can be expected to be accurate across a wide range of CBV values [34].
The purpose of the study presented here is to compare this novel DRCE-MRI technique with conventional DSCE-MRI approaches. The CBV maps generated from both methods have been compared in terms of image quality (spatial information, distortions etc.), pixel value distributions within normal brain and measured values within intraaxial cerebral tumours. The aim is to determine the suitability of these CBV maps for diagnosis, grading and the planning and guidance of procedures such as surgery and radiotherapy.
| Methods and materials |
|---|
|
|
|---|
|
Routine clinical T1 and T2 weighted anatomical imaging was performed in all patients prior to dynamic studies. Post-contrast T1 weighted volume images (TR 24 ms/TE 11 ms) were acquired at the end of the dynamic studies and used for qualitative comparison with parametric maps of CBV.
For acquisition of T1 weighted dynamic contrast-enhanced data, a 3D T1 FFE (T1 fast field echo) scanning sequence was applied, with an image matrix of 128 x 128 pixels and 25 slices in the axial plane. The field of view used was 230 mm (square), with a slice thickness of 6 mm with 3 mm overlap (Fourier interpolation), resulting in an effective slice thickness of 3 mm. TR was set at 4.2 ms, and TE at 1.2 ms. Three pre-contrast data sets were acquired for the baseline T1 calculation using flip angles of 2°, 10° and 35°. This was followed by a dynamic contrast-enhanced acquisition series at a flip angle of 35°, consisting of 120 scans with a temporal spacing of approximately 5 s. Gadolinium-based contrast agent (Gd-DTPA-BMA; OmniscanTM, Amersham Health AS, Oslo, Norway) was injected as a bolus over 4 s at a dose of 0.1 mmol kg1 of body weight.
Acquisition of T2* weighted data was carried out immediately after completion of T1 weighted data acquisition, with the time delay set at a fixed duration, so that pre-enhancement would minimize residual relaxivity effects ("T1 shine through") in the dynamic data [35]. T2* weighted (T2*W) dynamic contrast-enhanced data were acquired using a multislice 2D T2*W-FEEPI (field echo (
gradient echo) EPI) multi-shot sequence. The image matrix was 128 x 128 pixels, with nine slices in the axial plane. The field of view was again 230 mm (square), with slice thickness of 5 mm with 1 mm interslice gap, giving an effective slice thickness of 6 mm. TR was set at 440 ms, and TE at 30 ms. The dynamic contrast-enhanced series consisted of 52 scans at a flip angle of 35° and a temporal spacing of 1.8 s. A second dose of contrast agent (Gd-DTPA) was injected after the fifth dynamic scan using the same protocol described above.
For all scans the tumour was centred in the imaging volume, the imaging volume included the superior sagittal sinus to provide a vascular input function for analysis and the geometry of the T2*W acquisitions was selected from the preceding T1W acquisition so that images were spatially coincident with the central images in the T1W dynamic data-set. Due to its volumetric coverage the temporal resolution of the T1W dynamic sequence is almost three times slower than the multi-slice 2D T2*W dynamic sequence. Parallel imaging was not available at the time of these scans. The pre-enhancement required for the T2*W sequence meant that the T1W DCE acquisition was always carried out using the first contrast injection before the T2*W DCE acquisition, and it was therefore not possible to alternate the sequences.
Image analysis
Dynamic contrast-enhanced T1W data were analysed using the First Pass Leakage Profile (FPLP) method described by Li et al [33, 36]. This technique uses a shape analysis to decompose contrast concentration time course data into two separate components representing intravascular and extravascular contrast agent, allowing calculation of CBV free from the effects of contrast leakage (T1 CBV), and of the transfer coefficient (Kfp) for the passage of contrast between the plasma and the extracellular extravascular space (EES).
T2*W dynamic contrast-enhanced data were analysed using the technique described by Zhu et al [14], based on the techniques of Kassner et al [35]. Multi-slice maps of rate of change of T2* (
R2*) were calculated from the T2*W-FEEPI dynamic data signals for each dynamic phase and a gamma variate model [37] was used to fit the first pass
R2*(t) data. Relative cerebral blood volume (rCBV) maps were calculated by pixel-by-pixel integration of the resulting gamma variate curves.
Comparison of T1 CBV and T2* CBV
Visual comparison was performed between maps of T1 CBV, T2* CBV and standard post-contrast T1 weighted anatomical images. Qualitative analysis included calculation of lesion conspicuity, a calculation of the correlation coefficient between median CBV values from each technique and an attempt to produce pixel-by-pixel correlation values from one representative case.
Due to spatial distortions resulting from susceptibility effects in the T2*W acquisition it was not possible to reliably apply automated coregistration of T1W and T2*W images. Separate corresponding regions of interest (ROIs) were therefore manually defined on matched T1W and T2*W images. ROIs were drawn by an experienced neuroradiologist (TAP) and were defined on contrast-enhanced images from the late phase dynamic acquisitions. ROI definition included all enhancing components of the tumour but excluded areas of non-enhancement. Calculated values of CBV from T2*W images were normalized to a value of 1 for a group of voxels (to increase signal-to-noise) with an intensity value of greater than 99% of the maximum (assumed to be 100% CBV) to obtain an estimate of absolute CBV. Median tumour CBV values from T1 CBV and T2* CBV were calculated from each ROI, the values were plotted on a scattergram and the correlation between the values tested using Pearson's correlation coefficient. In addition direct subjective visual comparison of T1-CBV and T2*-CBV maps was performed using a standardized colour map to aid comparison.
ROIs were also drawn on normal appearing white matter on both T1 CBV and T2* CBV maps, and the mean (µ) and standard deviation (
) from these were compared with values from the enhancing tumour ROIs in the same maps, to ascertain an index of tumour conspicuity using the following expression:
|
|
The square-root of the CBV (sqrt (CBV)) maps was taken in order to produce estimates with uniform error and to expand the dynamic range of the data to make it easier to see any trends in the correlation between the two data sets. The volume of T1 sqrt(CBV) maps from a representative case was manually registered, using a linear affine transform, to the corresponding volume of T2* sqrt(CBV) maps. A single T2*W slice was chosen and the T1 sqrt(CBV) volume was resliced using a sinc 5 kernel for interpolation. The T1 sqrt(CBV) map was then normalized to the (already normalized) T2* sqrt(CBV) by finding the peak in the distribution of the logged CBV ratio over all voxels. A scattergram of the sqrt(CBV) of the two data sets in normal vascularized brain only was produced in the range 00.6 sqrt(CBV). This range was chosen for clarity, as the majority of sqrt(CBV) values fell within this range. A second scattergram was obtained using data from the same region of interest as the first, by taking the most similar value from the T1 sqrt(CBV) image to the pixel of interest in the T2* sqrt(CBV) image, corresponding to a half pixel linear interpolation in one of four directions (anterior, posterior, left, right) a pixel shuffle technique, as has previously been used to compare pre- and post-contrast MR images [38]. The technique allows the selection of the most likely match between the two data sets in order to minimize the errors introduced by the non-rigidity between the two volumetric acquisitions and broadening of large blood vessels in T2* data due to susceptibility effects.
| Results |
|---|
|
|
|---|
Figure 1
shows a typical representative axial slice through the centre of a high grade glioma (patient 8). A standard high resolution T1W volume post-contrast anatomical image showing the enhancing portion of the tumour and enhanced vasculature (Figure 1a
) is provided for comparison with parametric maps. Figure 1b
is a T2* CBV map at the same location as that of Figure 1a
and displays the distribution of CBV values within the brain and tumour tissue (note that the skull and scalp have been stripped from this image during processing). Normal vasculature, such as the branches of the middle cerebral artery and the cerebral veins, show high values of CBV as expected, and the enhancing portion of the tumour shows a heterogeneous distribution of CBV values. There is good differentiation of normal grey and white matter. Note that blood vessels appear much broader and smoother-edged on the T2* CBV map than on the high resolution anatomical image. There is also considerable signal drop-out and distortion in the basal prefrontal cortex due to susceptibility artefact resulting from air in the paranasal sinuses. Figure 1d
is a T1 CBV map at the same location as Figure 1a
. In contrast to Figure 1b
the map shows improved demonstration of high spatial frequency features and anatomical details, particularly vascular structure, corresponding closely to the anatomical image (Figure 1a
). The signal-to-noise ratio in normal grey and white matter in Figure 1d
is poorer than in the T2* CBV image but still allows clear subjective visual discrimination. The spatial distribution of CBV values in Figure 1b and 1d
is similar in normal brain, tumour and vessels. The FPLP method also generates maps of Kfp as shown in Figure 1c
indicating an intact BBB where Kfp is zero or consistent with noise, and areas of higher Kfp representing contrast leakage only within the tumour and choroid plexus.
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
Instead of using T2* weighted gradient echo sequences data may also be collected with T2 weighted (T2W) spin echo sequences. T2W spin echo techniques are relatively more sensitive to the small vessels than T2*W gradient-echo imaging techniques [42, 43], while being less prone to the distortion associated with susceptibility effects. On the other hand, T2*W imaging, which represents the effects of total blood volume in capillaries up to large vessels and weights all vessels approximately equally [44], is expected to be more suitable for evaluating brain tumours [45], and more useful for differentiating low-grade from high-grade gliomas than the T2W spin echo techniques [8]. The current study thus focused on the comparison between T2* CBV and T1 CBV in intra-axial cerebral tumours.
In view of the problems associated with DSCE-MRI techniques several groups have attempted to use relaxivity based dynamic images (DRCE-MRI) to derive estimates of CBV. The potential advantages of this approach include the lack of spatial distortion associated with T1W images and the potential for fast 3D imaging approaches without the need for EPI collections. The major problem with DRCE-MRI in enhancing tissues is the inability to separate the effects of intravascular and extravasated contrast on the observed signal change. This means that signal changes in areas of high capillary permeability will result to a large degree from contrast leakage and estimates of CBV will therefore be erroneously high.
Despite this disadvantage, Hacklander's group in the late nineties described the use of DRCE-MRI in cerebral tumours in a series of publications [2932]. This group dealt with the leakage problem by assuming that contrast leakage is so slow that it can be ignored over the time course of a single passage of a contrast bolus. They stated that: "gadopentate dimeglumine can almost be regarded to be an intravascular contrast agent, even in cases of a disturbed blood brain barrier". They concluded that T1 CBV measurements could be used in enhancing tissues. However, it is interesting to note that they were unable to demonstrate significant differences between grade III and grade IV glioma even though this has been shown by several groups using susceptibility based techniques [10, 14, 28]. In addition, careful comparison of the results from the T2* and T1 based techniques shows systematic overestimation of T1 CBV in tumours with high values (Figure 7 in [32]). Although the authors do not comment on this a similar observation was highlighted by Bruening et al [28], who described areas of apparently very high CBV on T1 CBV maps which were not seen on T2* CBV or T2 CBV. They commented: "In the high grade group, different values between T1 and T2 CBV maps were apparent", concluding: "Theoretically one would anticipate that in settings of blood brain barrier disruption there would be a tendency for T1 rCBV maps to cause elevated rCBV measurements. In contrast T2 rCBV maps tend to underestimate the apparent rCBV values in the presence of a blood brain barrier breakdown and may show false negative findings in the event of an active tumour recurrence."
In practice, as described above, the tendency for T2* CBV maps to underestimate CBV in areas of contrast leakage can be minimized by a variety of acquisition strategies which are discussed in detail by Kassner et al [35], so that the major residual disadvantages of T2* based methods are the distortion associated with susceptibility effects and the restrictions on imaging time. The sequences used here represent a compromise between spatial resolution and susceptibility sensitivity using a segmented EPI acquisition protocol to reduce scan time whilst minimizing the degree of spatial distortion. The use of optimized T2* based sequences such as PRESTO offers the opportunity to improve on this performance by combining echo shifted, segmented EPI collection and volume acquisition techniques in order to maximize temporal and spatial resolution whilst limiting sensitivity to susceptibility artefacts and distortions. Although these sequences are still not widely available and we have not used them in this study our experience with them indicates that the improvements in susceptibility distortion are relatively limited and spatial distortion remains a significant problem [35, 46].
The use of DRCE-MRI avoids problems associated with susceptibility-based spatial distortion and, in addition, it is possible to use reduced contrast doses compared with DSCE-MRI, although we have not explored that aspect in this study [28, 32]. The major problem with DRCE-MRI is that the signal changes resulting from intravascular contrast and from extravasated contrast occur in the same direction and, since they result from the same physical contrast mechanism, they cannot be separated by modifications of the acquisition technique. We have used the first pass leakage profile model described by Li et al [34, 36] to separate these effects and to generate leakage free T1 CBV maps. This analysis technique decomposes the contrast concentration time course data from DRCE-MRI into two components: the first due to intravascular contrast agent, and the second due to contrast agent leakage into the extravascular extracellular space. The technique works by using a constrained shape model of each component and decomposing the data to produce the optimized fit to these constraints. Although the technique was designed to produce estimates of transendothelial contrast transfer coefficient (Ktrans [47], or more specifically Kfp [48]) which are free of intravascular contrast effects, it also allows calculation of leakage-free CBV maps. The theoretical advantages in terms of temporal resolution, tissue coverage and freedom from image distortion have been the subject of the present study. However, it is equally important to know that CBV estimates from DRCE-MRI provide comparable biological information to those derived from conventional DSCE-MRI techniques. This study has confirmed that this is the case with close correlation between median values and also excellent pixel-by-pixel agreement, particularly when spatial distortion effects are reduced. The advantages of T1 techniques allow confident use of parametric maps for image guided surgical procedures and radiotherapy planning without any risk of error due to spatial distortion.
In conclusion we would recommend further investigation of T1-weighted DCE-MRI for the routine measurement of CBV in cerebral tumours. The technique avoids the risk of significant spatial distortion, provides biologically equivalent data to conventional T2* weighted DCE-MRI methods and has the added advantage of providing high-quality maps of the transfer coefficient [48, 49].
Received for publication June 14, 2006. Accepted for publication July 13, 2006.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L.-G. Vezina Imaging of Central Nervous System Tumors in Children: Advances and Limitations J Child Neurol, October 1, 2008; 23(10): 1128 - 1135. [Abstract] [PDF] |
||||
![]() |
BJR review of the year -- 2007 Br. J. Radiol., April 1, 2008; 81(964): 265 - 269. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| BJR | DMFR | IMAGING | ALL BIR JOURNALS |