Abstract
We assessed the performance of parametric fusion PET/MRI based on 11C-choline PET/CT and apparent diffusion coefficient (ADC) maps derived from diffusion-weighted MRI for the identification of primary prostate cancer. Methods: 11C-choline PET/CT and MRI were performed in 17 patients with untreated primary prostate cancer, followed by prostatectomy. Registration of in vivo imaging with histology was achieved using a mutual-information objective function and by performing ex vivo MRI of the prostatectomy specimen (obtained at 3 T) and whole-mount sectioning with block-face photography as intermediate steps. Data analysis included volumetrically registered whole-mount histology with Gleason scoring, 11C-choline, and ADC data (obtained at 1.5 T). Volumes of interest were defined on the basis of histologically proven tumor tissue to calculate tumor–to–benign prostate background ratios (TBRs) for 11C-choline, ADC, and a derived fusion PET/MRI parameter calculating the quotient of 11C-choline over ADC (PCHOL/ADC). Results: Fifty-one tumor nodules were identified at pathology. The TBRs for 11C-choline (P < 0.05) and PCHOL/ADC (P < 0.005) were significantly higher in prostate cancers with a Gleason score of ≥3 + 4 than with a Gleason score of ≤3 + 3 disease and controls. For Gleason ≥ 3 + 4, the ADC TBRs were significantly lower than controls and Gleason ≤ 3 + 3 disease (P < 0.05). The absolute value of TBRs obtained from Gleason ≥ 3 + 4 cancers increased from ADC to 11C-choline PET/CT and from 11C-choline PET/CT to PCHOL/ADC, with each step being statistically significant. Conclusion: Our data indicate that parametric PET/MRI using PCHOL/ADC improves lesion-to-background contrast (TBRs) of Gleason ≥ 3 + 4 disease, compared with 11C-choline PET/CT or diffusion-weighted MRI, and thus hold promise that parametric imaging performed on hybrid PET/MRI may further improve identification and localization of significant primary prostate cancer.
In the past 20 y, imaging of prostate cancer with CT, hybrid PET and CT (PET/CT), and MRI has rapidly developed for disease identification, staging, and management of treatment. Each modality has, however, well-known limitations inherent to anatomic and molecular imaging approaches that limit their diagnostic accuracy for the assessment of tumor stage and therapy response. For primary prostate cancer, several key clinical imaging scenarios exist: the identification of suspected prostate cancer (elevated prostate-specific antigen with previous negative biopsy), localization of clinically significant (more aggressive) primary disease, and detection of local and distant metastatic disease.
PET/CT with 11C- and 18F-choline has shown promise in the detection of prostate cancer, especially in the setting of recurrent and metastatic disease (1–4). Studies evaluating 11C-choline for the detection of primary prostate cancer have provided mixed results (5–8), potentially because of overlap with benign prostate hyperplasia and prostatitis (9) and a considerable uncertainty about the true location of disease in relation to the imaging findings they were felt to represent. To overcome these technical limitations, we recently developed a registration method using a standard mutual-information objective function and thin-plate spline warping to register PET/CT and histology in prostate cancer (10). This registration technique—using high-resolution ex vivo MRI of the prostatectomy specimen and whole-mount histology—allowed the assessment of imaging findings associated with various histologic features, most importantly the local Gleason score. Using this technique in a patient population specifically selected to minimize potential inflammation from recent prostate biopsies or clinical signs of prostatitis, we identified that 11C-choline preferentially localizes to aggressive primary prostate cancer when using tumor-to-normal background ratios, however, with considerable overlap to low-grade disease and benign tissues (11).
Anatomic MRI of primary prostate cancer has rapidly developed in the last decade. More recently, diffusion-weighted MRI (DWI) using apparent diffusion coefficient (ADC) maps has been applied. DWI techniques are based on the molecular mobility of water, which is sensitive to interaction with intracellular elements, macromolecules, cell membranes, cell density, and microstructural organization (12). DWI has identified restricted water diffusion in several malignancies, including prostate cancer (13,14). Recent studies performed in untreated primary prostate cancer suggested a subtle but significant negative correlation of ADCs with Gleason scores; thus, DWI ADC may be helpful in identifying peripheral zone cancers of intermediate and high risk (14–18).
To overcome specific limitations of morphologic and molecular imaging, hybrid PET/MRI scanners have most recently been introduced into clinical practice (19,20). Such systems promise optimal spatial and temporal coregistration of PET and MRI data, with the added benefit of performing functional MRI. It is still unknown whether such devices, compared with their established individual components performed separately, will indeed hold the promise of improving clinical imaging (21). Here, we investigated for the first time whether the combination of 11C-choline PET/CT and DWI may improve lesion detection of primary prostate cancer. In addition to using the mutual-information objective function, PET/MRI fusion also included registration of whole mount pathology with established registration errors (10). Although we are aware that PET/MRI fusion imaging as performed in this study is limited because of unavoidable misregistration, this study may still prove valuable in assessing the clinical potential of future hybrid PET/MRI in prostate cancer.
MATERIALS AND METHODS
Patient Population
The study was conducted with 17 men (mean age ± SD, 60.5 ± 7.2 y) with biopsy-proven, untreated localized adenocarcinoma of the prostate scheduled to undergo prostatectomy; 14 patients had been included in our prior report on 11C-choline (11). All patients received 12-sample prostate biopsies within 16 wk before enrollment, resulting in 3 or more cores positive for a summed Gleason score of 6 or greater from at least 1 lobe of the prostate. Patients with prior androgen ablation treatment, prostate biopsies performed less than 6 wk before PET/CT and MRI, previous external radiation treatment of the pelvic region, any prior malignancies, active inflammatory bowel disease, or evidence of prostatitis were not eligible.
In Vivo Imaging
In vivo imaging was performed within 14 d before radical prostatectomy (mean, 7.2 d; range, 3–13 d). In 14 cases, MRI and 11C-choline PET/CT were performed on the same day, whereas in the remaining 3 cases imaging was completed 1 d apart. Anatomic in vivo MRI was performed either on a 3-T system (Achieva; Philips Healthcare) using a torso–cardiac surface coil or on a 1.5-T system (Signa; GE Healthcare) using a pelvic surface coil combined with an endorectal coil. Parameters for the axial T2-weighted scans on the 3-T system were fast spin echo (echo train length, 15); repetition time/echo time, 2,800/80 ms; field of view (FOV), 280 × 280 mm; matrix, 388 × 307; and twenty-eight 3-mm slices. Parameters on the 1.5-T system for the axial T2-weighted scans were fast spin echo (echo train length, 16); repetition time/echo time, 3,050/120 ms; FOV, 140 × 140 mm; matrix, 256 × 194; and twenty-eight 3-mm slices.
Axial diffusion scans obtained on the 1.5-T system only were single-shot echo planar imaging; repetition time/echo time, 6,000/65 ms; FOV, 240 × 240 mm; 128 × 128 acquired matrix reconstructed to a 256 × 256 matrix; 4-mm slices through the prostate; and isotropic diffusion weighting at b-values 0 and 800 s/mm2. ADC maps were calculated by the standard formula (22):
PET/CT was performed on a Siemens Biograph classic scanner, which incorporates an ECAT HR+ PET camera (Siemens) with a 2-channel helical CT component (n = 12). In 5 cases, PET/CT was performed on a Biograph TrueV HiRez scanner with an extended FOV and a 6-channel helical CT component (Siemens Medical Solutions). The intrinsic axial resolutions were approximately 5.2 and 4.5 mm in full width at half maximum, respectively, for the Biograph classic scanner (23) and Biograph TrueV HiRez scanner (11). Five minutes after injection of 690 ± 26 MBq of 11C-choline, a limited body scan of the lower abdomen and pelvis was obtained (7 min per bed position), with the first bed position at the level of the prostate. A total of 2 bed positions were included on the Biograph TrueV and 3 on the Biograph classic PET scanner. Images were reconstructed in a 128 × 128 × 210 matrix, resulting in a 5.2 × 5.2 × 2.4 mm voxel dimension. Image reconstruction methods were identical on both scanners. Applying iterative (ordered-subset expectation maximization) algorithms (Fourier rebinning, 4 iterations, 8 subsets) with a 7-mm gaussian filter using a low-dose CT scan without intravenous or oral contrast for attenuation correction resulted in a similar effective axial reconstructed resolution on both scanners (11).
Image Coregistration
As described earlier, a standard registration method based on a mutual-information cost function and thin-plate spline deformation was applied (10,24). Although registration of in vivo modalities (T2-weighted MRI, DWI, PET/CT) is generally achieved without difficulties, our approach separated the difficult direct registration of histology and in vivo imaging (PET/CT, MRI) into achievable subregistration tasks involving intermediate ex vivo modalities (block-face photography and specimen MRI). Histology sections were corrected for 2-dimensional microtome-induced deformations by registering each section with its respective block-face photograph acquired immediately before cutting. Volumetric stacks of histology sections were then registered onto high-resolution (and high-information-content) ex vivo specimen MRI volumes. The final registration of the ex vivo specimen MRI volume onto the in vivo MRI reference integrated multiple discrete 2-dimensional histology slice volumes into the 3-dimensional in vivo reference space (T2-weighted MRI), with errors in the range of 2.26–3.74 mm (10).
Histologic Assessment
Fresh prostate specimens were fixed overnight in 10% neutral formalin. After specimen MRI at 3 T (10), the prostate was sectioned using a whole-mount technique and stained (hematoxylin and eosin [HE]) (25,26). Tumor maps were generated for each whole-mount slice to determine individual tumor foci in multifocal prostate cancers as previously described (27). Each tumor focus was assigned a primary and secondary Gleason grade and staged according to the 1998 guidelines of the American Joint Committee on Cancer (28). Tumors were graded according to Gleason scores, with ≤3 + 3 as low-, 3 + 4 as intermediate-, and ≥ 4 + 3 as high-grade tumors. Tumor borders were highlighted by the pathologist on HE slides. Benign hyperplasia was inhomogeneously present in prostatic tissues, favoring the central gland.
Image Data Analysis
Volumes of interest (VOI) were defined on the basis of registered HE histology, contouring the borders of each tumor lesion individually according to the pathologist's indications. Once HE slices were registered with T2-weighted MRI volumes, the tumor volume (cm3) was determined. Depending on the tumor size per lobe, generally 2 (in 3 cases, just 1) separate standardized VOIs (1.5 cm3) were defined in benign tissue in the peripheral zone and central gland, avoiding proximity to cancerous lesions. Those VOIs were then overlaid with the registered image volume (11C-choline or ADC); thus, the definition of tumor lesions and benign reference background was independent from particular imaging features identified visually on each modality. Standardized uptake value (SUV) and ADC values were obtained on a voxel-by-voxel basis from every VOI by calculating the mean ± SD of SUVs and ADCs.
For parametric fusion PET/MRI, coregistered 11C-choline SUV and ADC voxels were used to create parametric voxel ratio data with a simple quotient of SUV over ADC values.
A small value of 0.001 was added to the ADC value in the denominator of Equation 2, thus avoiding division by zero. Because the quotient of 11C-choline over ADC (PCHOL/ADC) was computed in the reference MRI space, voxel dimensions and image matrix were identical to the T2-weighted MRI.
Further analyses were performed after calculating the mean tumor–to–benign prostate background ratio (TBR) using contralateral histologically benign tissue (which included areas of benign hyperplasia) as a reference. To account for variations of signal intensities on ADC maps and 11C-choline uptake pattern within the prostate gland, smaller peripheral zone tumors were normalized by nonmalignant contralateral peripheral zone tissue, whereas larger tumors involving significant portions of the prostate lobe were normalized by the mean of the benign peripheral and central gland (11). For each modality (PET [SUV], MRI [ADC], and parametric fusion PET/MRI [PCHOL/ADC], the same coregistered VOI ratios (TBRs) were used.
Statistics
Analyses were conducted using the R statistical package (29). For all analyses, we used a linear mixed-effects (LME) model. A Box-Cox transformation (a log transformation was suggested) was performed on all outcomes (11C-choline, ADC, and PCHOL/ADC). To minimize the impact of partial-volume effects on the data analysis, we used the tumor volume as a weighting factor. Tumors with a measured tumor volume below 0.7 cm3 were down-weighted, whereas tumors ≥ 0.7 cm3 were assessed without weighting. The level of the threshold was set at 0.7 as a compromise between image resolution of PET and MRI data. The unit of correlation for the LME models was the patient because several VOIs were measured within each subject. Contrasts (differences in model parameter estimates associated with covariates) were used within the LME models for assessment of statistical significance using log-transformed TBRs to first compare each imaging method at 3 histologic levels (Gleason ≥ 3 + 4 vs. Gleason ≤ 3 + 3 vs. benign) and then the 3 imaging methods (PET vs. ADC vs. PCHOL/ADC) at a Gleason ≥ 3 + 4 level. Multiple comparisons were adjusted by Bonferroni method.
RESULTS
All patients were staged T2b to T3b at pathology. The mean weight of the prostate specimen was 47.6 ± 12.1 g. The mean prostate-specific antigen level before biopsy was 7.6 ± 3.8 ng/mL.
A total of 51 tumor nodules were identified at HE histology. In all but 1 patient, multifocal disease was present. Gleason scoring of the largest (primary) tumor lesion per patient varied among 3 + 4 (n = 7), 4 + 3 (n = 8), and 4 + 4 (n = 2). Of 34 multifocal (secondary) tumor lesions, 28 were classified as Gleason ≤ 3 + 3, 4 as Gleason 3 + 4, and 2 as 4 + 4 disease. Tumor volumes, as identified on registered histology on T2-weighted MRI, varied with Gleason score categories. Low-grade tumors were considerably smaller (mean volume, 0.44 cm3; range, 0.03–1.53 cm3; n = 29) than intermediate- and high-grade disease. The average volume of intermediate-grade lesions was 4.00 cm3 (range, 0.22–25.74 cm3; n = 10), whereas high-grade lesions had similar average volume (4.09 cm3) and range (0.34–12.63 cm3; n = 12). Thus, many tumor volumes, mainly low-grade tumors, were below the reconstructed resolution of 11C-choline PET/CT and DWI. In the statistical analyses, we accounted for partial-volume effects using the tumor volume (≤0.7 cm3) as weighting factor.
For 11C-choline, the mean TBR was significantly increased in Gleason ≥ 3 + 4 disease, compared with controls (P < 0.05) and significantly higher than Gleason ≤ 3 + 3 disease (P < 0.01) (Fig. 1A). Figure 1B shows that the mean ADC TBR decreased in Gleason ≥ 3 + 4 disease, compared with controls (P < 0.05), and was significantly lower than Gleason ≤ 3 + 3 disease (P < 0.05). For parametric PET/MRI using PCHOL/ADC, the mean TBR was significantly elevated in Gleason ≥ 3 + 4 disease, compared with controls (P < 0.005), and also significantly higher than Gleason ≤ 3 + 3 disease (P < 0.001) (Fig. 1C).
Figure 2 displays imaging results of all modalities in high-grade disease using an overlay of PET, DWI, histology, and PET/MRI fusion imaging (PCHOL/ADC) with T2-weighted MRI. A larger Gleason 4 + 3 lesion (volume, 6.78 cm3) was identified visually on the left on ADC and 11C-choline images. 11C-choline PET images represent SUVs, whereas DWI displays ADC values. Accordingly, the PET/MRI fusion images represent PCHOL/ADC. In addition, 1 separate low-volume Gleason 3 + 3 satellite lesion (volume, 0.36 cm3) found in the peripheral zone on the right side was missed on all imaging.
Lesion identification on ADC and 11C-choline imaging is strongly influenced by the image contrast between lesion and surrounding (benign) tissues. We therefore assessed whether the lesion-to-background contrast of Gleason ≥ 3 + 4 cancers differed for the 3 imaging modalities (PET, ADC, and fusion PET/MRI) using the imaging method as contrast in the LME model. We identified that the absolute TBRs of Gleason ≥ 3 + 4 cancers increased from ADC to 11C-choline PET/CT and from 11C-choline PET/CT to PCHOL/ADC, with each step being statistically significant (ADC vs. 11C-choline PET/CT, P < 0.02; ADC vs. PCHOL/ADC, P < 0.05; and 11C-choline PET/CT vs. PCHOL/ADC, P < 0.01). The positive effect on lesion identification is demonstrated in Figure 2.
DISCUSSION
The clinical behavior of primary prostate cancer ranges from microscopic, well-differentiated tumors to aggressive cancer with a high likelihood of invasion and metastasis. Clinical management can be most diverse and may range from definitive treatment options over localized treatments to watchful waiting. In the developed world, prostate-specific antigen testing has resulted in a dramatic stage-migration toward detection of lower-grade disease. In this situation, a reliable noninvasive detection of significant disease (30) could reduce unnecessary morbidity and treatment costs in patients with clinically indolent disease, thereby addressing the pressing issue of overtreatment of prostate cancer (31).
Imaging competes with the diagnostic transrectal prostate biopsy procedure, often performed repeatedly. Although the value of Gleason scoring based on final pathology obtained from prostatectomy specimens is undisputed as a reliable predictor of outcome, the accuracy of biopsy-based diagnosis is severely hampered by considerable sampling errors, resulting in significant under- or overestimation of the true Gleason score with limited concordance rates between 43% and 67% in contemporary studies (32–34).
We have recently shown in a partially overlapping patient population (n = 14) that—when TBRs are applied—11C-choline preferentially localizes to aggressive prostate cancer (11). Here, we demonstrate in 17 patients that similar results can be obtained from DWI using ADC maps. Our data and several other recent studies indicate that ADCs are decreasing with increasing prostate cancer aggressiveness (14–17,35,36). However, similar to our findings regarding 11C-choline, a considerable overlap between low-, intermediate-, and high-grade disease limits classification of individual tumor nodules.
Fusion of separately acquired PET and MRI has rarely been attempted in humans (37,38). Such approaches are severely hampered by unavoidable registration errors and temporal changes of functional image signals. Here, intermodality fusion between DWI and PET was assisted by high-resolution ex vivo MRI of the prostate specimen, which also offered the integration of registered histologic information. We have previously shown that registration between MRI and histology is possible with small errors (10), but with significant uncertainty in regard to lesion boundaries when using automated thresholding techniques (39). Fusion of in vivo imaging (MRI vs. PET/CT) is generally possible with sufficient precision; however, definitive measurements for the registration accuracy between in vivo PET/CT and MRI are lacking. In this study, in vivo PET/CT and MRI was performed mostly on the same day within a narrow time window, without endorectal coil and with voided urinary bladder. These precautions limited the potential for significant deformation of the prostate gland between scans.
Because Gleason ≥ 3 + 4 tumors, compared with low-grade cancers and benign prostate tissues, demonstrated significantly increased 11C-choline uptake and decreased ADCs, we evaluated the potential value of fusion PET/MRI for lesion detection and characterization. Our results suggest that parametric PET/MRI using PCHOL/ADC significantly increased the lesion-to-benign background contrast (expressed by TBRs) for Gleason ≥ 3 + 4 tumors, compared with those obtained from 11C-choline or ADC as a single modality. Despite improvements in characterization of Gleason ≥ 3 + 4 tumors, TBRs obtained from PCHOL/ADC in Gleason ≥ 3 + 4 still showed considerable overlap with low-grade disease and benign tissues, although to a lesser degree than TBRs from 11C-choline and ADC (Fig. 1). We attribute such overlap—at least in part—to an increased variability of fusion data, compared with single-modality data, which is related to unavoidable registration errors. In the case of small lesions, even relatively minute registration errors may influence parameter estimates from PET/MRI fusion imaging. We accounted for this uncertainty by down-weighting low-volume disease for the statistical analyses.
Performing registration of histology with in vivo imaging has profound effects on the visual appearance of image data. It limits the FOV to the smallest component of the registration process (histology); thus, in vivo imaging registered to the reference space (T2-weighted MRI) can be visualized only for a confined anatomic space around the prostate gland as seen in Figure 2. Reformatting (warping) lower-resolution DWI and PET volumes to match the higher resolution of the reference space reduces voxel sizes for DWI and PET while maintaining their respective numeric data content (10), but reformatted DWI and PET images display an inferior visual clarity. Reformatting histology forces many high-resolution voxels to be replaced with 1 low-resolution voxel bearing the average intensity. As a result, the reformatted histology, compared with original HE slices, is blurred but sufficient to identify the exact lesion borders for VOI definition based on the pathologist's indications.
CONCLUSION
These pilot data are the first, to our knowledge, to suggest that parametric PET/MRI is feasible in humans and may improve identification of significant primary prostate cancer, a relevant question in the clinical management of primary disease. Obviously, parametric PET/MRI is not limited to the particular parameter (PCHOL/ADC) presented here. Various other metabolic markers could be combined using diffusion- or perfusion-weighted MRI or MR spectroscopy to improve lesion detection or to enhance characterization for PET markers generally resulting in low contrast (such as hypoxia markers). Therefore, we anticipate parametric hybrid PET/MRI to be a versatile and potentially powerful imaging tool for prostate cancer and beyond.
DISCLOSURE STATEMENT
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
We thank the staff of the PET and MRI suites for their excellent technical support. This research was supported by grants 1P01CA87634 and P50CA069568 (Prostate Cancer SPORE) and 1P01CA87634 (National Institutes of Health). No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online Mar. 14, 2012.
- © 2012 by the Society of Nuclear Medicine, Inc.
REFERENCES
- Received for publication April 11, 2011.
- Accepted for publication December 8, 2011.