|
|
||||||||
Genitourinary Imaging |
1 From Imaging Science and Biomedical Engineering, University of Manchester, Stopford Building, Oxford Road, Manchester M13 9PT, England (D.L.B., C.R., G.J.M.P., C.E.H.); and Department of Clinical Oncology, Christie Hospital, Manchester, England (J.P.L.). Received December 23, 2003; revision requested February 24, 2004; revision received March 18; accepted April 8. Supported by a grant from the U.S. Department of Defense Prostate Cancer Research Program (PC991154). Address correspondence to D.L.B. (e-mail: david.buckley@man.ac.uk).
| ABSTRACT |
|---|
|
|
|---|
MATERIALS AND METHODS: Twenty-two patients between 57 and 76 years of age (mean age, 67 years) with histologically proved adenocarcinoma of the prostate were examined by using three-dimensional dynamic contrast-enhanced T1-weighted MR imaging at 1.5 T. The local research ethics committee approved this study, and written consent was obtained from all patients. Data from regions of interest drawn in tumor, normal-appearing peripheral zone tissue, and muscle were analyzed to provide estimates of perfusion, blood volume, interstitial volume, and microvascular permeabilitysurface area product. These estimates were compared by using the nonparametric Wilcoxon signed rank test.
RESULTS: Mean blood flow was significantly (P < .001) higher in 22 prostate tumors than in 20 contralateral peripheral zones (66 vs 32 mL/100 mL/min). Similarly, the interstitial distribution volume in tumors was enlarged compared with the interstitial distribution volume in normal peripheral zones (42 vs 27 mL/100 mL). Blood volume and microvascular permeabilitysurface area product values in tumors (1.0 mL/100 mL and 22 mL/100 mL/min, respectively) were similar to estimated values in peripheral zone tissue (1.5 mL/100 mL and 21 mL/100 mL/min, respectively).
CONCLUSION: These findings show considerable promise for isolating vascular characteristics of prostate cancer.
© RSNA, 2004
Index terms: Magnetic resonance (MR), tissue characterization, 844.12117, 844.12146 Neoplasms, blood supply, 844.316, 844.32 Prostate neoplasms, MR, 844.12146
| INTRODUCTION |
|---|
|
|
|---|
It is likely that the early enhancement seen in cancer is due in part to angiogenesis associated with early tumor growth. The growth of prostate cancer is associated with the development of a rich blood supply fed by a large network of immature, leaky blood vessels (8). The density of this network is related to tumor grade and metastatic potential (9). Because conventional methods of assessing angiogenesis require the processing of pathologic specimens, there is considerable interest in developing noninvasive tools to assess the microvascular characteristics of prostate cancer in vivo (10).
A link between contrast material uptake in solid tumors and their microvascular characteristics certainly exists, but the precise relationship between tracer kinetics, which are often measured in terms of heuristic measurements, and microvessel density, which is measured at pathologic examination, remains unclear (11,12). The assessment of tracer kinetics in absolute physiologic terms may help to address this question (13). Increasing microvessel density will lead to an increase in blood flow, blood volume, and the surface area of vessel walls. An upturn in vascular endothelial growth factor production is likely to increase the permeability of these vessel walls.
Blood flow, blood volume, and microvascular permeabilitysurface area product are all, in principle, quantifiable through analysis of contrast-enhanced MR imaging data by using distributed-parameter tracer kinetics models (10,14). To our knowledge, these have yet to be applied in vivo with MR imaging in a substantial patient group, although results of a recent preliminary computed tomography (CT)-based study (15) are promising.
The use of tracer kinetics modeling in the analysis of tomographic images of the prostate is widespread, from the earliest study involving positron emission tomography (PET) (16) through a series of studies involving MR imaging (1719) and CT (15,20). Although the results of the studies that involved MR imaging indicated that MR imaging has sufficient spatial and contrast resolution to depict cancer, normal peripheral zone tissue, and benign prostatic hyperplasia within individual glands, there is limited information about the microvascular characteristics of these tissues.
Results of studies involving compartmental models have shown that peripheral zone cancers have higher extraction-flow products (transfer constants [Ktrans]) (21) than normal peripheral zone tissue (17,18). This finding is supported by the generally higher blood flow measured in tumor-containing prostates than in prostates containing benign prostatic hyperplasia alone (16). Differences between tumors and normal tissues that contain benign prostatic hyperplasia are less clear cut in the central portion of the gland, although a recent MR imagingbased study revealed potential differences in flow-related parameters (22). Thus, the purpose of our study was to use contrast-enhanced MR imaging and a distributed-parameter tracer kinetics model for prospectively evaluating the vascular characteristics of prostate cancer.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Johnson and Wilson (25) introduced a simple distributed-parameter model in which the internal capillary space contains a tracer, the concentration of which varies in both time and space (along the capillary length). Tracer in the capillary space is in exchange with an external interstitial space that behaves as a well-mixed compartment (ie, there is no spatial variation in tracer concentration). In 1998, St Lawrence and Lee (24) introduced an adiabatic approximation to this model, providing the following time-domain solution suitable for the analysis of imaging data:
|
is the mean microvascular transit time, which is equal to the blood volume, Vb, divided by F. Ve is the interstitial tracer distribution volume, and u is a dummy integration variable. The extraction fraction E is calculated with the following equation: 1 exp[PS/F(1 hct)], where PS is the microvascular permeabilitysurface area product. The product EF(1 hct) is also referred to as Ktrans (21).
This model contains four unknown parameters (F, E, Ve, and
) that may be estimated with curve fitting. From these, the parameters PS and Vb may be evaluated. This model has since been used successfully for the analysis of both MR imaging and CT data (15,26,27).
Patients
Twenty-two consecutive men with histologically proved adenocarcinoma of the prostate who were to be treated with neoadjuvant hormonal treatment and conformal radiation therapy were recruited into this study between June 2001 and February 2003. The local research ethics committee approved the study, and written consent was obtained from all men. Histologic diagnosis was rendered with tissue obtained at transrectal ultrasonographically guided biopsy in 21 men and with tissue obtained at the time of a transurethral resection of the prostate in one man. The clinical stages at presentation were T1c through T3b, with the modal T stage being T3 in 14 patients. All men had had negative results at isotope bone scanning. The mean Gleason score was 6 (range, 58), the mean prostate-specific antigen level was 27.5 ng/mL (range, 674 ng/mL), and the mean age of the patients was 67 years (range, 5776 years). After MR imaging and in accordance with local protocol, all men were treated with neoadjuvant hormonal therapy involving a luteinizing hormonereleasing hormone analogue for 34 months and three-dimensional conformal radiation therapy.
MR Imaging
MR imaging was performed at 1.5 T (Gyroscan NT/Intera; Philips Medical Systems, Best, the Netherlands) by using a four-element phased-array receiver coil wrapped around the pelvis. These men were expected to undergo follow-up MR imaging after radiation therapy; thus, the use of an endorectal receiver was precluded (28). After the acquisition of scout images, data encompassing the abdomen and pelvis in both the transverse and coronal planes were acquired by using T1-weighted spin-echo sequences. These images were used to assess the possibility of extraprostatic spread. Subsequently, T2-weighted fast spin-echo images of the prostate were obtained in the sagittal and transverse planes with the following parameters: field of view, 200 x 200 mm; section thickness, 3 mm; intersection gap, 0.3 mm; matrix, 256 x 190; repetition time msec/effective echo time msec, 4111/135; echo train length, 19; and number of signals acquired, six. After the acquisition of quantitative data during contrast agent administration (detailed below), postcontrast T1-weighted spin-echo images of the abdomen and pelvis were again acquired in the transverse plane.
Quantitative data encompassing the pelvic region from the apex of the prostate gland to the bladder were acquired by using a transverse three-dimensional spoiled gradient-echo sequence (ie, a T1-weighted fast-field-echo sequence with an elliptical k-space shutter, a field of view of 300 x 300 x 100 mm, a matrix [after interpolation] of 128 x 128 x 20, and a repetition time msec/echo time msec of 2.5/0.86). Before contrast agent administration, the native T1 relaxation time was mapped by using a series of volume acquisitions (each obtained with four signals acquired) with flip angles of 2°, 10°, 20°, and 30°. After acquisition of the T1 relaxation map data, a dynamic series was acquired by using a 30° flip angle and no averaging over a period of approximately 4 minutes with a total of 100106 volumes acquired; one volume was acquired every 2.3 seconds. Contrast agent was administered through a cannula that had been placed in the antecubital vein before imaging. After the acquisition of the first five volumes of the dynamic series, 0.1 mmol of gadodiamide (Omniscan; Amersham-Health, Little Chalfont, Buckinghamshire, England) per kilogram of body weight was injected at 3 mL/sec by using an MR imagingcompatible power injector (Spectris; Medrad, Indianola, Pa). The injection of contrast material was immediately followed by injection of a similar volume of normal saline.
Data Collection and Statistical Analysis
Data were transferred off-line for quantitative analysis. For each patient, 4-pixel (22-mm2) regions of interest were placed (by C.R.) in the external iliac or femoral arteries (to provide a vascular input function). With the aid of the biopsy results and T2-weighted MR images, a radiologist (C.E.H.) also drew regions over muscle (the internal obturator muscle, which was selected as a control tissue) that ranged in size from 148 to 956 mm2. The radiologist then drew regions over hypointense areas within the peripheral zone that were suspected of containing prostate tumor, and, where possible, regions over normal contralateral peripheral zone tissue. These regions were drawn to encompass as much tumor and normal contralateral peripheral zone as was visible in the selected section. Signal intensity variations in these regions were converted to temporal changes in gadodiamide concentration by using estimates of native T1 calculated by using the precontrast data (29,30).
Finally, the distributed-parameter model described by St Lawrence and Lee (24) was fitted to the data. Multiple fits were performed to each tissue region by using a range of starting points (systematic variation in the initial estimates of
). The final solution selected was that which minimized the sum of squared differences between the data and the model fit (26,27). To determine whether differences existed between the vascular characteristics of tumor, normal prostate, and muscle in the patients examined, parameter estimates were compared by using the nonparametric Wilcoxon signed rank test (SPSS for Windows, version 11.5; SPSS, Chicago, Ill). P < .05 was considered to indicate a statistically significant difference.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Validation
Our estimates of native tissue T1 calculated by using variable-flip-angle three-dimensional gradient-echo acquisitions (29) compare well with previous measurements obtained in both muscle (32) and the prostate (33). Single-capillary distributed-parameter models have previously been used in the analysis of PETdata (23), and the St Lawrence and Lee model specifically has been used in both MR imagingbased studies (26) and CT-based studies (15). Use of a variant of the model for measuring blood flow at CT in comparison with microsphere measurements has been validated in two separate studies (34,35). Although the St Lawrence and Lee model has not yet been formally validated for use in evaluating the prostate, results of simulation studies indicate that this model performs well in comparison to both independent multiple path distributed-parameter models and commonly used compartmental models (27).
In accordance with results of previous work (17,18), the combination of E and F (Ktrans) measured in the current study was significantly higher in tumor than in normal peripheral zone tissue. The estimates of blood flow in the present study compare well with estimates obtained by using PET (16). Even the small decrease in blood volume that was associated with cancer in our study follows the trend observed in that previous study (16).
Finally, the model-based estimates of blood flow and microvessel permeabilitysurface area product correlate well with the predictions made by pathologists in previous studies (36). Although tumors contain a greater number of vessels, specifically in angiogenic "hot spots" (8), than do normal tissues, the vessel surface area density in tumors is actually reduced (36). This reflects the size and shape of tumor microvessels, which tend to be broad and less branched than those in normal tissue or low-grade tumors. This reduced vessel surface area in high-grade cancers is offset by increased blood flow and vessel wall permeability (36)physiologic changes that are obscured in the estimate of permeabilitysurface area product and may explain the similarities between estimates made in tumors and in normal tissue.
Study Limitations
Although concordant imaging and biopsy results provide a highly specific marker of cancer in the peripheral zone (37), this study lacked a standard of reference such as whole-mount prostatectomy specimens (38). Essential to the analysis of any experiment is data quality. The data acquired in the external iliac and femoral arteries and the prostate gland were of good quality, but the signal obtained from the internal obturator muscle was characterized by a poor signal-to-noise ratio and sporadic motion artifacts. Motion interfered with a number of the signal-time courses analyzed, and future studies would probably benefit if image registration was performed before data analysis. The signal-to-noise ratio is conventionally addressed at the expense of temporal resolution.
A reduction in temporal resolution was not an option in the present study because it was necessary to acquire data very rapidly so that a distributed-parameter model could be used. Indeed, the acquisition time that we used was already longer than that recommended in similar studies (39). The limited signal-to-noise ratio is a common problem in prostate MR imaging and may be partly addressed by the use of an endorectal coil (38). Although the use of an endorectal coil is certain to improve the signal from the prostate, it is unlikely to increase the signal from the artery or muscle. Moreover, use of the endorectal coil may be contraindicated in the assessment of patients, such as the patients in our study, who have undergone radiation therapy (28). Alternatively, the signal-to-noise ratio could be improved by using a high-field-strength MR imaging system (40). These limitations may have influenced the accuracy and precision of the parameter estimates that we obtained.
A previous study revealed that estimates of blood volume obtained by using this implementation of the St Lawrence and Lee model were low (27), and the results in the present study fall below previous estimates obtained by using PET (16). This may reflect noise or an inadequate temporal resolution in the experimental data (39) that led to limitations in definition and integration of the arterial input (20). It may also indicate the existence of a limited transendothelial water exchange rate (41,42). The influence of such effects may be reduced by further decreasing the sensitivity of the imaging sequence to water exchange (by increasing the flip angle with so-called exchange minimization [41]). Such an approach further compromises the signal-to-noise ratio of the sequence, and a practical trade-off must be achieved.
Clinical Implications
Our results provide new quantitative physiologic data concerning the prostate gland. These methods, which include measurement of an arterial input function, three-dimensional coverage of the entire prostate gland, and absolute estimates of both flow and permeabilitysurface area product, represent an advance on those used in previous MR imagingbased studies (7,18). Moreover, similar CT examinations lack tissue coverage (only one or two sections may be examined) and the intrinsic tissue contrast at MR imaging that enables discrimination of tumor from normal peripheral zone tissue (15,20).
These techniques will allow further investigators to characterize tissue types, improve diagnosis, and categorize disease prognosis. The promising results obtained by using qualitative or heuristic methods could be assessed to identify their underlying physiologic basis (7,22). Tissue characterization in the present study was limited by the method of region selection, which was dictated by the patient group examined. Nevertheless, concordant imaging and biopsy results providing a highly specific marker of cancer in the peripheral zone (37) guided selection of prostate tumor regions. Future work would benefit if MR imaging findings were compared with findings in whole-mount prostatectomy specimens, which would constitute an improved standard of reference (38).
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Author contributions: Guarantor of integrity of entire study, D.L.B.; study concepts and design, D.L.B., J.P.L., C.E.H.; literature research, D.L.B.; clinical studies, C.E.H., J.P.L.; data acquisition, D.L.B., C.E.H., J.P.L.; data analysis/interpretation, D.L.B., C.R., G.J.M.P., J.P.L.; statistical analysis, D.L.B.; manuscript preparation, D.L.B.; manuscript definition of intellectual content, D.L.B., G.J.M.P., C.E.H.; manuscript editing, D.L.B., G.J.M.P., J.P.L., C.E.H.; manuscript revision/review and final version approval, all authors
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N VAN AS, E CHARLES-EDWARDS, A JACKSON, S JHAVAR, S REINSBERG, N DESOUZA, D DEARNALEY, M BAILEY, A THOMPSON, T CHRISTMAS, et al. Correlation of diffusion-weighted MRI with whole mount radical prostatectomy specimens Br. J. Radiol., June 1, 2008; 81(966): 456 - 462. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. K. Kim, B. K. Park, H. M. Lee, S. S. Kim, and E. Kim MRI Techniques for Prediction of Local Tumor Progression After High-Intensity Focused Ultrasonic Ablation of Prostate Cancer Am. J. Roentgenol., May 1, 2008; 190(5): 1180 - 1186. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Casciani, E. Polettini, E. Carmenini, I. Floriani, G. Masselli, L. Bertini, and G. F. Gualdi Endorectal and Dynamic Contrast-Enhanced MRI for Detection of Local Recurrence After Radical Prostatectomy Am. J. Roentgenol., May 1, 2008; 190(5): 1187 - 1192. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. A. Mullani, R. S. Herbst, R. G. O'Neil, K. L. Gould, B. J. Barron, and J. L. Abbruzzese Tumor Blood Flow Measured by PET Dynamic Imaging of First-Pass 18F-FDG Uptake: A Comparison with 15O-Labeled Water-Measured Blood Flow J. Nucl. Med., April 1, 2008; 49(4): 517 - 523. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. N. Bloch, E. Furman-Haran, T. H. Helbich, R. E. Lenkinski, H. Degani, C. Kratzik, M. Susani, A. Haitel, S. Jaromi, L. Ngo, et al. Prostate Cancer: Accurate Determination of Extracapsular Extension with High-Spatial-Resolution Dynamic Contrast-enhanced and T2-weighted MR Imaging Initial Results Radiology, October 1, 2007; 245(1): 176 - 185. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Jackson, J. P.B. O'Connor, G. J.M. Parker, and G. C. Jayson Imaging Tumor Vascular Heterogeneity and Angiogenesis using Dynamic Contrast-Enhanced Magnetic Resonance Imaging Clin. Cancer Res., June 15, 2007; 13(12): 3449 - 3459. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Hricak, P. L. Choyke, S. C. Eberhardt, S. A. Leibel, and P. T. Scardino Imaging Prostate Cancer: A Multidisciplinary Perspective Radiology, April 1, 2007; 243(1): 28 - 53. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. J. Choi, J. K. Kim, N. Kim, K. W. Kim, E. K. Choi, and K.-S. Cho Functional MR Imaging of Prostate Cancer RadioGraphics, January 1, 2007; 27(1): 63 - 75. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Futterer, S. W. T. P. J. Heijmink, T. W. J. Scheenen, J. Veltman, H. J. Huisman, P. Vos, C. A. H.-V. de Kaa, J. A. Witjes, P. F. M. Krabbe, A. Heerschap, et al. Prostate Cancer Localization with Dynamic Contrast-enhanced MR Imaging and Proton MR Spectroscopic Imaging Radiology, November 1, 2006; 241(2): 449 - 458. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |