Treatment planning for prostate implants using magnetic-resonance spectroscopy imaging

MR Imaging and MR Spectroscopy in Prostate Cancer Management

Currently, endorectal coil MR imaging has the ability to improve accuracy in staging of localized prostate cancer. The addition of MR spectroscopic imaging has further improved the sensitivity of MR imaging for intraprostatic tumor localization. Additional refinements and techniques are expected to further improve the performance of MR imaging for prostate cancer imaging and to aid in patient management. Further studies are required to identify the ideal role for MR imaging in the diagnosis and management of prostate cancer.

IMRT boost dose planning on dominant intraprostatic lesions: Gold marker-based three-dimensional fusion of CT with dynamic contrast-enhanced and 1H-spectroscopic MRI

PURPOSE

To demonstrate the theoretical feasibility of integrating two functional prostate magnetic resonance imaging (MRI) techniques (dynamic contrast-enhanced MRI [DCE-MRI] and 1H-spectroscopic MRI [MRSI]) into inverse treatment planning for definition and potential irradiation of a dominant intraprostatic lesion (DIL) as a biologic target volume for high-dose intraprostatic boosting with intensity-modulated radiotherapy (IMRT).

METHODS AND MATERIALS

In 5 patients, four gold markers were implanted. An endorectal balloon was inserted for both CT and MRI. A DIL volume was defined by DCE-MRI and MRSI using different prostate cancer-specific physiologic (DCE-MRI) and metabolic (MRSI) parameters. CT-MRI registration was performed automatically by matching three-dimensional gold marker surface models with the iterative closest point method. DIL-IMRT plans, consisting of whole prostate irradiation to 70 Gy and a DIL boost to 90 Gy, and standard IMRT plans, in which the whole prostate was irradiated to 78 Gy were generated. The tumor control probability and rectal wall normal tissue complication probability were calculated and compared between the two IMRT approaches.

RESULTS

Combined DCE-MRI and MRSI yielded a clearly defined single DIL volume (range, 1.1-6.5 cm3) in all patients. In this small, selected patient population, no differences in tumor control probability were found. A decrease in the rectal wall normal tissue complication probability was observed in favor of the DIL-IMRT plan versus the plan with IMRT to 78 Gy.

CONCLUSION

Combined DCE-MRI and MRSI functional image-guided high-dose intraprostatic DIL-IMRT planned as a boost to 90 Gy is theoretically feasible. The preliminary results have indicated that DIL-IMRT may improve the therapeutic ratio by decreasing the normal tissue complication probability with an unchanged tumor control probability. A larger patient population, with more variations in the number, size, and localization of the DIL, and a feasible mechanism for treatment implementation has to be studied to extend these preliminary tumor control and toxicity estimates.

Treatment of neuroblastoma using the fused imaging guided radiotherapy (FIGURA) system

PURPOSE

The purpose of this study was to describe our department's experience with the fused imaging-guided radiotherapy (FIGURA) system for planning radiation treatment of high-risk neuroblastoma.

PATIENTS AND METHODS

Between 1999 and 2002, 11 patients received radiation therapy as consolidation after chemotherapy in 9 and for palliation in 2. Diagnostic metaiodobenzylguanidine (MIBG) imaging was used, which is specific for neuroblastoma, to identify the residual tumor, followed by computed tomography scanning in the radiation treatment position. The FIGURA software fused the images obtained by the 2 modalities and transferred the result to a 3-dimensional radiation treatment planning system. Radiation was delivered at a total dose of 25.2 Gy according to the FIGURA.

RESULTS

Five patients achieved complete remission and 2 partial remission; 3 were stabilized. One child with a highly rapid progressive course died of the disease.

CONCLUSION

FIGURA is a new, feasible technique for defining target volumes. By using standard hospital equipment, it is possible to treat residual disease identified by sensitive metaiodobenzylguanidine imaging and localized with the anatomic computed tomography scan. Treating a more accurate target volume spares normal tissue and organs and minimizes side effects.

[MR spectroscopic imaging for evaluation of prostate cancer]

MR spectroscopic imaging (MRSI) provides a noninvasive method of evaluating metabolic markers of prostate cancer or healthy prostatic tissue (the metabolites choline and citrate), and is performed in conjunction with high-resolution MR anatomic imaging. Multiple studies have showed the incremental role of MRSI combined with the anatomical information provided by MRI for assessment of cancer location and extent within the prostate, staging, and cancer aggressiveness. In addition, MRSI has a potential role for pre- and post-treatment evaluation in non surgical patients. Ongoing technical developments show the potential role of MRSI for guidance of biopsies or focal treatment. Further developments - including new 3T technology - will likely provide improved spectral resolution for better prostate cancer detection and characterization.

Quantification of prostate MRSI data by model-based time domain fitting and frequency domain analysis

This paper compares two spectral processing methods for obtaining quantitative measures from in vivo prostate spectra, evaluates their effectiveness, and discusses the necessary modifications for accurate results. A frequency domain analysis (FDA) method based on peak integration was compared with a time domain fitting (TDF) method, a model-based nonlinear least squares fitting algorithm. The accuracy of both methods at estimating the choline + creatine + polyamines to citrate ratio (CCP:C) was tested using Monte Carlo simulations, empirical phantom MRSI data and in vivo MRSI data. The paper discusses the different approaches employed to achieve the quantification of the overlapping choline, creatine and polyamine resonances. Monte Carlo simulations showed induced biases on the estimated CCP:C ratios. Both methods were successful in identifying tumor tissue, provided that the CCP:C ratio was greater than a given (normal) threshold. Both methods predicted the same voxel condition in 94% of the in vivo voxels (68 out of 72). Both TDF and FDA methods had the ability to identify malignant voxels in an artifact-free case study using the estimated CCP:C ratio. Comparing the ratios estimated by the TDF and the FDA, the methods predicted the same spectrum type in 17 out of 18 voxels of the in vivo case study (94.4%).

Image-fusion of MR spectroscopic images for treatment planning of gliomas

1H magnetic resonance spectroscopic imaging (MRSI) can improve the accuracy of target delineation for gliomas, but it lacks the anatomic resolution needed for image fusion. This paper presents a simple protocol for fusing simulation computer tomography (CT) and MRSI images for glioma intensity-modulated radiotherapy (IMRT), including a retrospective study of 12 patients. Each patient first underwent whole-brain axial fluid-attenuated-inversion-recovery (FLAIR) MRI (3 mm slice thickness, no spacing), followed by three-dimensional (3D) MRSI measurements (TE/TR: 144/1000 ms) of a user-specified volume encompassing the extent of the tumor. The nominal voxel size of MRSI ranged from 8 x 8 x 10 mm3 to 12 x 12 x 10 mm3. A system was developed to grade the tumor using the choline-to-creatine (Cho/Cr) ratios from each MRSI voxel. The merged MRSI images were then generated by replacing the Cho/Cr value of each MRSI voxel with intensities according to the Cho/Cr grades, and resampling the poorer-resolution Cho/Cr map into the higher-resolution FLAIR image space. The FUNCTOOL processing software was also used to create the screen-dumped MRSI images in which these data were overlaid with each FLAIR MRI image. The screen-dumped MRSI images were manually translated and fused with the FLAIR MRI images. Since the merged MRSI images were intrinsically fused with the FLAIR MRI images, they were also registered with the screen-dumped MRSI images. The position of the MRSI volume on the merged MRSI images was compared with that of the screen-dumped MRSI images and was shifted until agreement was within a predetermined tolerance. Three clinical target volumes (CTVs) were then contoured on the FLAIR MRI images corresponding to the Cho/Cr grades. Finally, the FLAIR MRI images were fused with the simulation CT images using a mutual-information algorithm, yielding an IMRT plan that simultaneously delivers three different dose levels to the three CTVs. The image-fusion protocol was tested on 12 (six high-grade and six low-grade) glioma patients. The average agreement of the MRSI volume position on the screen-dumped MRSI images and the merged MRSI images was 0.29 mm with a standard deviation of 0.07 mm. Of all the voxels with Cho/Cr grade one or above, the distribution of Cho/Cr grade was found to correlate with the glioma grade from pathologic finding and is consistent with literature results indicating Cho/Cr elevation as a marker for malignancy. In conclusion, an image-fusion protocol was developed that successfully incorporates MRSI information into the IMRT treatment plan for glioma.

Three-dimensional conformal external beam radiotherapy compared with permanent prostate implantation in low-risk prostate cancer based on endorectal magnetic resonance spectroscopy imaging and prostate-specific antigen level

PURPOSE

To evaluate the metabolic response by comparing the time to resolution of spectroscopic abnormalities (TRSA) and the time to prostate-specific antigen level in low-risk prostate cancer patients after treatment with three-dimensional conformal external beam radiotherapy (3D-CRT) compared with permanent prostate implantation (PPI). Recent studies have suggested that the treatment of low-risk prostate cancer yields similar results for patients treated with 3D-CRT or PPI.

METHODS AND MATERIALS

A total of 50 patients, 25 in each group, who had been treated with 3D-CRT or PPI, had undergone endorectal magnetic resonance spectroscopy imaging before and/or at varying times after therapy. The 3D-CRT patients had received radiation doses of > or =72 Gy compared with 144 Gy for the PPI patients. The spectra from all usable voxels were examined for detectable levels of metabolic signal, and the percentages of atrophic and cancerous voxels were tabulated.

RESULTS

The median time to resolution of the spectroscopic abnormalities was 32.2 and 24.8 months and the time to the nadir prostate-specific antigen level was 52.4 and 38.0 months for the 3D-CRT and PPI patients, respectively. Of the 3D-CRT patients, 92% achieved negative endorectal magnetic resonance spectroscopy imaging findings, with 40% having complete metabolic atrophy. All 25 PPI patients had negative endorectal magnetic resonance spectroscopy imaging findings, with 60% achieving complete metabolic atrophy.

CONCLUSION

The results of this study suggest that metabolic and biochemical responses of the prostate are more pronounced after PPI. Our results have not proved PPI is more effective at curing prostate cancer, but they have demonstrated that it may be more effective at destroying prostate metabolism.

Registration of MR prostate images with biomechanical modeling and nonlinear parameter estimation

Magnetic resonance imaging (MRI) and magnetic resonance spectroscopic imaging (MRSI) have been shown to be very useful for identifying prostate cancers. For high sensitivity, the MRI/MRSI examination is often acquired with an endorectal probe that may cause a substantial deformation of the prostate and surrounding soft tissues. Such a probe is removed prior to radiation therapy treatment. To register diagnostic probe-in magnetic resonance (MR) images to therapeutic probe-out MR images for treatment planning, a new deformable image registration method is developed based on biomechanical modeling of soft tissues and estimation of uncertain tissue parameters using nonlinear optimization. Given two-dimensional (2-D) segmented probe-in and probe-out images, a finite element method (FEM) is used to estimate the deformation of the prostate and surrounding tissues due to displacements and forces resulting from the endorectal probe. Since FEM requires tissue stiffness properties and external force values as input, the method estimates uncertain parameters using nonlinear local optimization. The registration method is evaluated using images from five balloon and five rigid endorectal probe patient cases. It requires on average 37 s of computation time on a 1.6 GHz Pentium-M PC. Comparing the prostate outline in deformed probe-out images to corresponding probe-in images, the method obtains a mean Dice Similarity Coefficient (DSC) of 97.5% for the balloon probe cases and 98.1% for the rigid probe cases. The method improves significantly over previous methods (P < 0.05) with greater improvement for balloon probe cases with larger tissue deformations.

Simultaneous beam geometry and intensity map optimization in intensity-modulated radiation therapy

PURPOSE

In current intensity-modulated radiation therapy (IMRT) plan optimization, the focus is on either finding optimal beam angles (or other beam delivery parameters such as field segments, couch angles, gantry angles) or optimal beam intensities. In this article we offer a mixed integer programming (MIP) approach for simultaneously determining an optimal intensity map and optimal beam angles for IMRT delivery. Using this approach, we pursue an experimental study designed to (a) gauge differences in plan quality metrics with respect to different tumor sites and different MIP treatment planning models, and (b) test the concept of critical-normal-tissue-ring--a tissue ring of 5 mm thickness drawn around the planning target volume (PTV)--and its use for designing conformal plans.

METHODS AND MATERIALS

Our treatment planning models use two classes of decision variables to capture the beam configuration and intensities simultaneously. Binary (0/1) variables are used to capture "on" or "off" or "yes" or "no" decisions for each field, and nonnegative continuous variables are used to represent intensities of beamlets. Binary and continuous variables are also used for each voxel to capture dose level and dose deviation from target bounds. Treatment planning models were designed to explicitly incorporate the following planning constraints: (a) upper/lower/mean dose-based constraints, (b) dose-volume and equivalent-uniform-dose (EUD) constraints for critical structures, (c) homogeneity constraints (underdose/overdose) for PTV, (d) coverage constraints for PTV, and (e) maximum number of beams allowed. Within this constrained solution space, five optimization strategies involving clinical objectives were analyzed: optimize total intensity to PTV, optimize total intensity and then optimize conformity, optimize total intensity and then optimize homogeneity, minimize total dose to critical structures, minimize total dose to critical structures and optimize conformity simultaneously. We emphasize that the objectives that include optimizing conformity make use of the critical-normal-tissue-ring. Three tumor sites: head-and-neck, pediatric brain, and prostate are used for comparison.

RESULTS

The critical-normal-tissue-ring acts as a good device for enforcing conformity. Trends in the characteristics and quality of plans resulting from each model were observed. Attempts to reduce dose to critical structures tend to worsen PTV conformity (1.542 to 3.092) and homogeneity (1.223 to 1.984), depending on the relative size and spatial distance of the critical structures to the PTV. When the critical structures are relatively small compared with the PTV (as in the case for the pediatric brain tumor, where each is less than 15% in volume), dose reduction to critical structures is accompanied by much worse scores in conformity (2.482) and homogeneity (1.984). When the critical structures are larger, as in the case of head-and-neck (approximately 50%), the conformity and homogeneity deterioration is less significant (1.542 and 1.239, respectively). There is a clear tradeoff between homogeneity, conformity, and minimum dose to organs at risk (OARs). For head-and-neck and pediatric brain tumor, the model that minimizes total dose to critical structures and optimizes conformity simultaneously offers a compromise among these factors, resulting in reduced critical structure dose with conformal and homogeneous plans. In the prostate case, the tumor is smaller than the two large nearby critical structures, and all models provide very homogeneous PTV dose distribution. However, minimizing dose to critical structures worsens conformity, as it spreads the radiation to the area surrounding the PTV. The maximum dose to the critical structures also increases slightly. Compared with plans used in the clinic which generally have uniformly spaced beam angles, the optimal clinically acceptable plans obtained via the methods herein do not have equispaced beams. The optimal beam angles returned appear to be nonintuitive, and depend on PTV size and geometry and the spatial relationship between the tumor and critical structures.

CONCLUSIONS

The MIP model described allows simultaneous optimization over the space of beamlet fluence weights and beam and couch angles. Based on experiments with tumor data, this approach can return good plans that are clinically acceptable and practical. This work distinguishes itself from recent IMRT research in several ways. First, in previous methods beam angles are selected before intensity map optimization. Herein, we employ 0/1 variables to model the set of candidate beams, and thereby allow the optimization process itself to select optimal beams. Second, instead of incorporating dose-volume criteria within the objective function as in previous work, herein, a combination of discrete and continuous variables associated with each voxel provides a mechanism to strictly enforce dose-volume criteria within the constraints. Third, using the construct of critical-normal-tissue-ring within the objective function can enhance the achievement of conformal plans. Based on the three tumor sites considered, it appears that volume and spatial geometry with respect to the PTV are important factors to consider when selecting objectives to optimize, and in estimating how well suited a particular model is for achieving a specified goal.

Biological imaging in radiation oncology

The goal of this study was to discuss the value of integrating biological imaging (PET, SPECT, MRS etc.) in radiation treatment planning and monitoring. Studies in patients with brain tumors have shown that, compared to CT and MRI alone, the image fusion of CT/MRI and amino acid SPECT or PET allows a more correct delineation of gross tumor volume (GTV) and planning target volume (PTV). For FDG-PET comparable results with different techniques are reported in the literature also for bronchial carcinoma, ear-nose-and-throat tumors, and cervical carcinoma, or, in the case of MRS, for prostate cancer. Imaging of hypoxia, cell proliferation, apoptosis, tumor angiogenesis, and gene expression leads to the identification of differently aggressive areas of a biologically inhomogeneous tumor mass that can be individually and more appropriately targeted using innovative IMRT. Thus, a biological, inhomogeneous dose distribution can be generated, the so-called dose painting. In addition, the biological imaging can play a significant role in the evaluation of the therapy response after radiochemotherapy. Clinical studies in ear-nose-and-throat tumors, bronchial carcinoma, esophagus carcinoma, and cervical carcinoma suggest that the sensitivity and specificity of FDG-PET for the therapy response are higher compared to anatomical imaging (CT and MRI). Clinical and experimental studies are required to define the real impact of these investigations in radiation treatment planning, and especially in the evaluation of therapy response.

Narrow band deformable registration of prostate magnetic resonance imaging, magnetic resonance spectroscopic imaging, and computed tomography studies

PURPOSE

Endorectal (ER) coil-based magnetic resonance imaging (MRI) and magnetic resonance spectroscopic imaging (MRSI) is often used to obtain anatomic and metabolic images of the prostate and to accurately identify and assess the intraprostatic lesions. Recent advancements in high-field (3 Tesla or above) MR techniques affords significantly enhanced signal-to-noise ratio and makes it possible to obtain high-quality MRI data. In reality, the use of rigid or inflatable endorectal probes deforms the shape of the prostate gland, and the images so obtained are not directly usable in radiation therapy planning. The purpose of this work is to apply a narrow band deformable registration model to faithfully map the acquired information from the ER-based MRI/MRSI onto treatment planning computed tomography (CT) images.

METHODS AND MATERIALS

A narrow band registration, which is a hybrid method combining the advantages of pixel-based and distance-based registration techniques, was used to directly register ER-based MRI/MRSI with CT. The normalized correlation between the two input images for registration was used as the metric, and the calculation was restricted to those points contained in the narrow bands around the user-delineated structures. The narrow band method is inherently efficient because of the use of a priori information of the meaningful contour data. The registration was performed in two steps. First, the two input images were grossly aligned using a rigid registration. The detailed mapping was then modeled by free form deformations based on B-spline. The limited memory Broyden-Fletcher-Goldfarb-Shanno algorithm (L-BFGS), which is known for its superior performance in dealing with high-dimensionality problems, was implemented to optimize the metric function. The convergence behavior of the algorithm was studied by self-registering an MR image with 100 randomly initiated relative positions. To evaluate the performance of the algorithm, an MR image was intentionally distorted, and an attempt was then made to register the distorted image with the original one. The ability of the algorithm to recover the original image was assessed using a checkerboard graph. The mapping of ER-based MRI onto treatment planning CT images was carried out for two clinical cases, and the performance of the registration was evaluated.

RESULTS

A narrow band deformable image registration algorithm has been implemented for direct registration of ER-based prostate MRI/MRSI and CT studies. The convergence of the algorithm was confirmed by starting the registration experiment from more than 100 different initial conditions. It was shown that the technique can restore an MR image from intentionally introduced deformations with an accuracy of approximately 2 mm. Application of the technique to two clinical prostate MRI/CT registrations indicated that it is capable of producing clinically sensible mapping. The whole registration procedure for a complete three-dimensional study (containing 256 x 256 x 64 voxels) took less than 15 min on a standard personal computer, and the convergence was usually achieved in fewer than 100 iterations.

CONCLUSIONS

A deformable image registration procedure suitable for mapping ER-based MRI data onto planning CT images was presented. Both hypothetical tests and patient studies have indicated that the registration is reliable and provides a valuable tool to integrate the ER-based MRI/MRSI information to guide prostate radiation therapy treatment.

Accuracy of finite element model-based multi-organ deformable image registration

As more pretreatment imaging becomes integrated into the treatment planning process and full three-dimensional image-guidance becomes part of the treatment delivery the need for a deformable image registration technique becomes more apparent. A novel finite element model-based multiorgan deformable image registration method, MORFEUS, has been developed. The basis of this method is twofold: first, individual organ deformation can be accurately modeled by deforming the surface of the organ at one instance into the surface of the organ at another instance and assigning the material properties that allow the internal structures to be accurately deformed into the secondary position and second, multi-organ deformable alignment can be achieved by explicitly defining the deformation of a subset of organs and assigning surface interfaces between organs. The feasibility and accuracy of the method was tested on MR thoracic and abdominal images of healthy volunteers at inhale and exhale. For the thoracic cases, the lungs and external surface were explicitly deformed and the breasts were implicitly deformed based on its relation to the lung and external surface. For the abdominal cases, the liver, spleen, and external surface were explicitly deformed and the stomach and kidneys were implicitly deformed. The average accuracy (average absolute error) of the lung and liver deformation, determined by tracking visible bifurcations, was 0.19 (s.d.: 0.09), 0.28 (s.d.: 0.12) and 0.17 (s.d.: 0.07) cm, in the LR, AP, and IS directions, respectively. The average accuracy of implicitly deformed organs was 0.11 (s.d.: 0.11), 0.13 (s.d.: 0.12), and 0.08 (s.d.: 0.09) cm, in the LR, AP, and IS directions, respectively. The average vector magnitude of the accuracy was 0.44 (s.d.: 0.20) cm for the lung and liver deformation and 0.24 (s.d.: 0.18) cm for the implicitly deformed organs. The two main processes, explicit deformation of the selected organs and finite element analysis calculations, require less than 120 and 495 s, respectively. This platform can facilitate the integration of deformable image registration into online image guidance procedures, dose calculations, and tissue response monitoring as well as performing multi-modality image registration for purposes of treatment planning.

Current imaging paradigms in radiation oncology

The technological revolution in imaging during recent decades has transformed the way image-guided radiation therapy is performed. Anatomical imaging (plain radiography, computed tomography, magnetic resonance imaging) greatly improved the accuracy of delineating target structures and has formed the foundation of 3D-based radiation treatment. However, the treatment planning paradigm in radiation oncology is beginning to shift toward a more biological and molecular approach as advances in biochemistry, molecular biology, and technology have made functional imaging (positron emission tomography, nuclear magnetic resonance spectroscopy, optical imaging) of physiological processes in tumors more feasible and practical. This review provides an overview of the role of current imaging strategies in radiation oncology, with a focus on functional imaging modalities, as it relates to staging and molecular profiling (cellular proliferation, apoptosis, angiogenesis, hypoxia, receptor status) of tumors, defining radiation target volumes, and assessing therapeutic response. In addition, obstacles such as imaging-pathological validation, optimal timing of post-therapy scans, spatial and temporal evolution of tumors, and lack of clinical outcome studies are discussed that must be overcome before a new era of functional imaging-guided therapy becomes a clinical reality.

Prostate Depiction at Endorectal MR Spectroscopic Imaging: Investigation of a Standardized Evaluation System

PURPOSE

To investigate the accuracy and interobserver variability of a standardized evaluation system for endorectal three-dimensional (3D) magnetic resonance (MR) spectroscopic imaging of the prostate.

MATERIALS AND METHODS

The human research committee approved the study, and all patients provided written informed consent. Endorectal MR imaging and MR spectroscopic imaging were performed in 37 patients before they underwent radical prostatectomy. For the 22 patients with good or excellent MR spectroscopic imaging data, step-section histopathologic tumor maps were used to identify spectroscopic voxels of unequivocally benign (n = 306) or malignant (n = 81) peripheral zone tissue. Two independent spectroscopists, unaware of all other findings, scored the spectra of the selected voxels by using a scale of 1 (benign) to 5 (malignant) that was based on standardized metabolic criteria. Descriptive statistical, receiver operating characteristics (ROC), and kappa statistical analyses of the data obtained by both readers were performed by using two definitions of cancer: one based on a voxel score of 3-5 and the other based on a score of 4 or 5.

RESULTS

The scoring system had good accuracy (74.2%-85.0%) in the differentiation between benign and malignant tissue voxels, with areas under the ROC curve of 0.89 for reader 1 and 0.87 for reader 2. Specificities of 84.6% and 89.3% were achieved when a voxel score of 4 or 5 was used to identify cancer, and sensitivities of 90% and 93% were achieved when a score of 3-5 was used to identify cancer. Readers demonstrated excellent interobserver agreement (kappa values, 0.79 and 0.80).

CONCLUSION

The good accuracy and excellent interobserver agreement achieved by using the standardized five-point scale to interpret peripheral zone metabolism demonstrate the potential effectiveness of using metabolic information to identify prostate cancer, and the clinical usefulness of this system warrants testing in prospective clinical trials of MR imaging combined with MR spectroscopic imaging.

Mapping of the prostate in endorectal coil-based MRI/MRSI and CT: A deformable registration and validation study

The endorectal coil is being increasingly used in magnetic resonance imaging (MRI) and MR spectroscopic imaging (MRSI) to obtain anatomic and metabolic images of the prostate with high signal-to-noise ratio (SNR). In practice, however, the use of endorectal probe inevitably distorts the prostate and other soft tissue organs, making the analysis and the use of the acquired image data in treatment planning difficult. The purpose of this work is to develop a deformable image registration algorithm to map the MRI/MRSI information obtained using an endorectal probe onto CT images and to verify the accuracy of the registration by phantom and patient studies. A mapping procedure involved using a thin plate spline (TPS) transformation was implemented to establish voxel-to-voxel correspondence between a reference image and a floating image with deformation. An elastic phantom with a number of implanted fiducial markers was designed for the validation of the quality of the registration. Radiographic images of the phantom were obtained before and after a series of intentionally introduced distortions. After mapping the distorted phantom to the original one, the displacements of the implanted markers were measured with respect to their ideal positions and the mean error was calculated. In patient studies, CT images of three prostate patients were acquired, followed by 3 Tesla (3 T) MR images with a rigid endorectal coil. Registration quality was estimated by the centroid position displacement and image coincidence index (CI). Phantom and patient studies show that TPS-based registration has achieved significantly higher accuracy than the previously reported method based on a rigid-body transformation and scaling. The technique should be useful to map the MR spectroscopic dataset acquired with ER probe onto the treatment planning CT dataset to guide radiotherapy planning.

Role of ProstaScint for brachytherapy in localized prostate adenocarcinoma

ProstaScint (CYT-356 or capromab pendetide, Cytogen) is an 111In-labeled monoclonal mouse antibody specific for prostate-specific membrane antigen, a prostate transmembrane glycoprotein that is upregulated in prostate adenocarcinoma. ProstaScint scans are US Food and Drug Administration approved for pretreatment evaluation of metastatic disease in high-risk patients. They are also approved for post-prostatectomy assessment of recurrent disease in patients with a rising prostate-specific antigen level. This review explores the literature on ProstaScint and its use in guiding the treatment of prostate cancer. A novel technique for identifying areas of cancer within the prostate using ProstaScint images fused with pelvic computed tomography scans is also described. The identification of areas of high antibody signal provides targets for radiotherapeutic dose escalation, with the overall goals of improving treatment outcome while preserving adjacent tissue structures and decreasing treatment morbidity.

MRI-guided HDR prostate brachytherapy in standard 1.5T scanner

PURPOSE

Magnetic resonance imaging (MRI) provides superior visualization of the prostate and surrounding anatomy, making it the modality of choice for imaging the prostate gland. This pilot study was performed to determine the feasibility and dosimetric quality achieved when placing high-dose-rate prostate brachytherapy catheters under MRI guidance in a standard "closed-bore" 1.5T scanner.

METHODS AND MATERIALS

Patients with intermediate-risk and high-risk localized prostate cancer received MRI-guided high-dose-rate brachytherapy boosts before and after a course of external beam radiotherapy. Using a custom visualization and targeting program, the brachytherapy catheters were placed and adjusted under MRI guidance until satisfactory implant geometry was achieved. Inverse treatment planning was performed using high-resolution T(2)-weighted MRI.

RESULTS

Ten brachytherapy procedures were performed on 5 patients. The median percentage of volume receiving 100% of prescribed minimal peripheral dose (V(100)) achieved was 94% (mean, 92%; 95% confidence interval, 89-95%). The urethral V(125) ranged from 0% to 18% (median, 5%), and the rectal V(75) ranged from 0% to 3.1% (median, 0.3%). In all cases, lesions highly suspicious for malignancy could be visualized on the procedural MRI, and extracapsular disease was identified in 2 patients.

CONCLUSION

High-dose-rate prostate brachytherapy in a standard 1.5T MRI scanner is feasible and achieves favorable dosimetry within a reasonable period with high-quality image guidance. Although the procedure was well tolerated in the acute setting, additional follow-up is required to determine the long-term safety and efficacy of this approach.

Inverse planning for HDR prostate brachytherapy used to boost dominant intraprostatic lesions defined by magnetic resonance spectroscopy imaging

PURPOSE

To dose escalate selected regions inside the prostate without compromising the dose coverage of the prostate and the protection to the urethra, rectum, and bladder for prostate cancer patients treated with high-dose-rate brachytherapy.

METHODS AND MATERIALS

Magnetic resonance imaging combined with magnetic resonance spectroscopy imaging was used to differentiate between normal and malignant prostate and define cancer-validated dominant intraprostatic lesions (DIL) on 10 patients. The DILs were then contoured on the planning scans (CT or MRI based, 5 patients each), and our inverse planning dose optimization algorithm (called IPSA) was used to generate dose distributions for 3 different boost levels. Dose-volume histograms of the target and each organ at risk were compared with optimized plans without DIL boost.

RESULTS

Combined MRI/magnetic resonance spectroscopic imaging identified 2 DILs in 8/10 of the 10 patients studied and a single DIL in the remaining 2 patients. The average prostate dose coverage V100 was 97% (sigma = 1.0%). When the minimum DIL dose requested was 120% of the prescribed dose, the average DIL V120 was 97.1% (sigma = 1.8%). For a boost value of 150%, the average V150 ranged from 77.8% to 86.1%, depending on the upper limit of the dose constraints. The bladder V50 increased by 1%, independently of the boost levels. The absolute increases in V50 for the rectum varied from 1% to 3%, depending on the boost level. The urethra V120 were increased by 13.4% and 32.5% for the lowest and highest boost levels, respectively.

CONCLUSION

The DIL dose can be escalated to a minimum of 120% while the entire prostate is treated simultaneously, without increasing the dose to surrounding normal tissues. Higher boost levels between 150% and 170% are feasible, but with slightly larger doses delivered to the rectum and urethra.

Optimization of tumour control probability for heterogeneous tumours in fractionated radiotherapy treatment protocols

We find the dose distribution that maximizes the tumour control probability (TCP) for a fixed mean tumour dose per fraction. We consider a heterogeneous tumour volume having a radiation response characterized by the linear quadratic model with heterogeneous radiosensitivity and repopulation rate that may vary in time. Using variational calculus methods a general solution is obtained. We demonstrate the spatial dependence of the optimal dose distribution by explicitly evaluating the solution for different functional forms of the tumour properties. For homogeneous radiosensitivity and growth rate, we find that the dose distribution that maximizes TCP is homogeneous when the clonogen cell density is homogeneous, while for a heterogeneous initial tumour density we find that the first dose fraction is inhomogeneous, which homogenizes the clonogen cell density, and subsequent dose fractions are homogeneous. When the tumour properties have explicit spatial dependence, we show that the spatial variation of the optimized dose distribution is insensitive to the functional form. However, the dose distribution and tumour clonogen density are sensitive to the value of the repopulation rate. The optimized dose distribution yields a higher TCP than a typical clinical dose distribution or a homogeneous dose distribution.

Deformable image registration for the use of magnetic resonance spectroscopy in prostate treatment planning

PURPOSE

There is now convincing evidence that prostate cancer cells lack the ability to produce and accumulate citrate. Using magnetic resonance spectroscopy imaging (MRSI), regions of absent or low citrate concentration in the prostate can be visualized at a resolution of a few mm. This new advancement provides not only a tool for early diagnosis and screening but also the opportunity for preferential targeting of radiation to regions of high tumor burden in the prostate. The differences in the shape and location of the prostate between MRSI imaging and treatment have been the major obstacle in integrating MRSI in radiation therapy treatment planning. The purpose of this study is to develop a reliable method for deforming the prostate and surrounding regions from the geometry of MRSI imaging to the geometry of treatment planning, so that the regions of high tumor burden identified by the MRSI study can be faithfully transferred to the images used for treatment planning.

METHODS AND MATERIALS

Magnetic resonance spectroscopy imaging studies have been performed on 2 prostate cancer patients using a commercial MRSI system with an endorectal coil and coupling balloon. At the end of each study, we also acquired the MRI of the pelvic region at both the deformed state where the prostate is distorted by the endorectal balloon and the resting state with the endorectal balloon deflated and removed. The task is to find a three-dimensional matrix of transformation vectors for all volume elements that links the two image sets. We have implemented an optimization method to iteratively optimize the transformation vectors using a Newton-Ralphson algorithm. The objective function is based on the mutual information. The distorted images using the transformation vectors are compared with the images acquired at the resting conditions.

RESULTS AND DISCUSSION

The algorithm is capable of performing the registration automatically without the need for intervention. It does not require manual contouring of the organs. By applying the algorithm to multiple image sets of different patients, we found a good agreement between the images transformed from those acquired at the deformed state and those acquired at resting conditions. The computation time required for achieving the registration is in the range of a half-hour (for image size: 256 pixels x 256 pixels x 25 slices). However, the space of registration can be restricted to speed up the process.

CONCLUSION

In this article, we described a three-dimensional deformable image registration method to automatically transform images from the deformed imaging state to resting state. Our examples show that this method is feasible and useful to the treatment planning system.

The robustness of dose distributions to displacement and migration of 125I permanent seed implants over a wide range of seed number, activity, and designs

PURPOSE

To investigate the robustness of permanent prostate implant dosimetry for various (125)I seed activities and various seed models. The dosimetric impact of seed misplacement and seed migration (seed loss) is also taken into account using various standard dose indices.

METHODS AND MATERIALS

A dose-based inverse planning algorithm is used for automated dosimetric plan creation (45-60 s per plan) and provides an unbiased way to compare the robustness of various optimal dosimetric plans. Seed misplacement and seed migration are simulated by way of Monte Carlo, based on the measured displacement distributions from clinical postimplant cases. Plans were generated for seed activities between 0.2 and 1.4 mCi (0.25 to 1.78 U) and for 11 different seed models.

RESULTS

The numbers of seeds and needles are shown to decrease rapidly for a seed activity between 0.3 mCi and 0.6 mCi (0.38 and 0.76 U). The loss in V100, from 100%, because of seed misplacement is below 10% for an apparent activity ranging from 0.2 to 0.9 mCi (0.25 to 1.14 U). A minimum degradation in V100 is observed around 0.6-0.7 mCi (0.76-0.89 U). D90 increases from 150 to 170 Gy between 0.3 and 0.7 mCi (0.38 and 0.89 U) and decreases afterward to fall below 140 Gy at higher activity. V200 and D10 to the target volume both show an increase in hot spots up to 0.7 mCi, and then decrease linearly at higher activities for all seed models. V200 and D10 to the urethra remain about constant for all seed activities up to 0.8 mCi (1.02 U), at which point they start to decrease. All seed models follow this general trend.

CONCLUSIONS

Plans were shown to be robust to misplacement and migration of seeds over a wide range of seed activity and for various seed models. With a properly tuned inverse planning algorithm able to ensure the dose coverage and protection for the organs at risk in the presence of placement errors (displacement and migration), the choice of a preferred seed activity, in a range up to about 0.7 mCi (0.89 U), is open. The upper part of this range offers the opportunity to significantly reduce the number of seeds and needles, thus reducing surgical trauma to the patient, saving time in an operating room planning setting, and reducing the cost of a permanent prostate implant procedure.

Proton high-resolution magic angle spinning NMR analysis of fresh and previously frozen tissue of human prostate

The previously observed improvement in spectral resolution of tissue proton NMR with high-resolution magic angle spinning (HRMAS) was speculated to be due largely to freeze-thawing artifacts resulting from tissue storage. In this study, 12 human prostate samples were analyzed on a 14.1T spectrometer at 3 degrees C, with HRMAS rates of 600 and 700 Hz. These samples were measured fresh and after they were frozen for 12-16 hr prior to thawing. The spectral linewidths measured from fresh and previously frozen samples were identical for all metabolites except citrate and acetate. The metabolite intensities of fresh and freeze-thawed samples depend on the quantification procedures used; however, in this experiment the differences of means were <30%. As expected, it was found that tissue storage impacts tissue quality for pathological analysis, and HRMAS conditions alone are not sufficiently destructive to impair pathological evaluation. Furthermore, although storage conditions affect absolute metabolite concentrations in NMR analysis, relative metabolite concentrations are less affected.

Combined quantitative dynamic contrast-enhanced MR imaging and1H MR spectroscopic imaging of human prostate cancer

PURPOSE

To differentiate prostate carcinoma from healthy peripheral zone and central gland using quantitative dynamic contrast-enhanced (DCE) magnetic resonance (MR) imaging and two-dimensional (1)H MR spectroscopic imaging (MRSI) combined into one clinical protocol.

MATERIALS AND METHODS

Twenty-three prostate cancer patients were studied with a combined DCE-MRI and MRSI protocol. Cancer regions were localized by histopathology of whole mount sections after radical prostatectomy. Pharmacokinetic modeling parameters, K(trans) and k(ep), as well as the relative levels of the prostate metabolites citrate, choline, and creatine, were determined in cancer, healthy peripheral zone (PZ), and in central gland (CG).

RESULTS

K(trans) and k(ep) were higher (P < 0.05) in cancer and in CG than in normal PZ. The (choline + creatine)/citrate ratio was elevated in cancer compared to the PZ and CG (P < 0.05). While a (choline + creatine)/citrate ratio above 0.68 was found to be a reliable indicator of cancer, elevated K(trans) was only a reliable cancer indicator in the diagnosis of individual patients. K(trans) and (choline + creatine)/citrate ratios in cancer were poorly correlated (Pearson r(2) = 0.07), and thus microvascular and metabolic abnormalities may have complementary value in cancer diagnosis.

CONCLUSION

The combination of high-resolution spatio-vascular information from dynamic MRI and metabolic information from MRSI has excellent potential for improved localization and characterization of prostate cancer in a clinical setting. J. Magn. Reson. Imaging 2004;20:279-287.

Intraoperative dynamic dose optimization in permanent prostate implants

PURPOSE

With the advent of intraoperative optimized planning, the treatment of prostate cancer with permanent implants has reached an unprecedented level of dose conformity. However, because of well-documented (and unavoidable) inaccuracies in seed placement into the gland, carrying out a plan results in a large degree of variability relative to the intended dose distribution. This brings forth the need to periodically readjust the plan to allow for the real positions of seeds already implanted. In this paper, an algorithm for performing this task, hereby described as intraoperative dynamic dose optimization (IDDO), is presented and assessed.

METHODS AND MATERIALS

The general scheme for performing IDDO consists of three steps: (1) at some point during the implant, coordinates of implanted seeds are identified; (2) seed images are projected onto the reference frame of the ultrasound images for planning; and (3) the plan is reoptimized. Work on the first two steps is reported elsewhere. Here, we focus on the strategy for implementing the reoptimization step. An optimal treatment plan is first obtained based on initial operating room-acquired ultrasound images. We analyze the sensitivity and effect of the IDDO procedure with respect to the total number of reoptimizations performed. Specifically, we consider reoptimizing 2, 3, and 4 times. When two reoptimizations are used, half of the seeds from the initial optimal plan are implanted. The first reoptimization is performed on the remaining possible seed positions, and all the seeds designated in this reoptimized plan are implanted. The second (final) reoptimization is done on the remaining unused seed positions to ensure 100% coverage of the gland and to eliminate possible cold spots in the gland. Similarly, when three reoptimization steps are used, one-third of the seeds from the initial optimized plan, one-half of the seeds from the first reoptimization, and all seeds from the second reoptimization are implanted. The third (final) reoptimization is performed to assist in eliminating possible cold spots. Reoptimizing four times proceeds in a like manner. Fifteen patient cases are used for comparison. Strict dose bounds of 100% and 120% of the prescription dose are imposed on the urethra, and 100% coverage is imposed on the prostate volume. To assist in achieving good conformity, prostate contour points are assigned a target upper dose bound of 150% of the prescription dose.

RESULTS

A two-way comparison is performed: (a) initial optimized plan, (b) IDDO plan. Postimplant dose analysis, coverage and conformity measures, as well as actual dose received by urethra and rectum are used to gauge the results. The initial optimized plan consistently provides 93% prescription dose coverage to the gland with average conformity index of 1.32. The urethra dose ranges within 100% to 150%, and the maximum dose delivered to the rectum reaches 91% of the prescription dose. On average, about 50% of the urethra receives more than 120% of the prescription dose, and 19% of the rectum volume receives more than the 78% upper dose limit. For the IDDO plan, 100% postimplant coverage with 1.16 conformity is achieved. Urethra and rectum dose is maintained within the prescribed 100% to 120% range and 78% upper bound, respectively.

CONCLUSIONS

With real-time treatment planning, it is possible to dynamically reoptimize treatment plans to account for actual seed positions (as opposed to planned positions) and needle-induced swelling to the gland during implantation. Postimplant analysis shows that the final seed configuration resulting from the IDDO method yields improved dosimetry. The algorithmic design ensures that one can achieve complete coverage while maintaining good conformity, thus sparing excess radiation to external tissue. The study also provides evidence of the possibility of morbidity reduction to urethra and rectum (because of reduced dose delivered to these structures) via the use of IDDO planning. Clinical studies are needed to validate the importance of our approach.

Intraoperative dynamic dosimetry for prostate implants

This paper describes analytic tools in support of a paradigm shift in brachytherapy treatment planning for prostate cancer--a shift from standard pre-planning to intraoperative planning using dosimetric feedback based on the actual deposited seed positions within the prostate. The method proposed is guided by several desiderata: (a) bringing both planning and evaluation in the operating room (i.e. make post-implant evaluation superfluous) therefore making rectifications--if necessary--still achievable; (b) making planning and implant evaluation consistent by using the same imaging system (ultrasound); and (c) using only equipment commonly found in a hospital operating room. The intraoperative dosimetric evaluation is based on the fusion between ultrasound images and 3D seed coordinates reconstructed from fluoroscopic projections. Automatic seed detection and registration of the fluoroscopic and ultrasound information, two of the three key ingredients needed for the intraoperative dynamic dosimetry optimization (IDDO), are explained in detail. The third one, the reconstruction of 3D coordinates from projections, was reported in a previous article. The algorithms were validated using a custom-designed phantom with non-radioactive (dummy) seeds. Also, fluoroscopic images were taken at the conclusion of an actual permanent prostate implant and compared with data on the same patient obtained from radiographic-based post-implant evaluation. To offset the effect of organ motion the comparison was performed in terms of the proximity function of the two seed distributions. The agreement between the intra- and post-operative seed distributions was excellent.

Image guidance for precise conformal radiotherapy

PURPOSE

To review the state of the art in image-guided precision conformal radiotherapy and to describe how helical tomotherapy compares with the image-guided practices being developed for conventional radiotherapy.

MATERIALS AND METHODS

Image guidance is beginning to be the fundamental basis for radiotherapy planning, delivery, and verification. Radiotherapy planning requires more precision in the extension and localization of disease. When greater precision is not possible, conformal avoidance methodology may be indicated whereby the margin of disease extension is generous, except where sensitive normal tissues exist. Radiotherapy delivery requires better precision in the definition of treatment volume, on a daily basis if necessary. Helical tomotherapy has been designed to use CT imaging technology to plan, deliver, and verify that the delivery has been carried out as planned. The image-guided processes of helical tomotherapy that enable this goal are described.

RESULTS

Examples of the results of helical tomotherapy processes for image-guided intensity-modulated radiotherapy are presented. These processes include megavoltage CT acquisition, automated segmentation of CT images, dose reconstruction using the CT image set, deformable registration of CT images, and reoptimization.

CONCLUSIONS

Image-guided precision conformal radiotherapy can be used as a tool to treat the tumor yet spare critical structures. Helical tomotherapy has been designed from the ground up as an integrated image-guided intensity-modulated radiotherapy system and allows new verification processes based on megavoltage CT images to be implemented.

Towards integrating functional imaging in the treatment of prostate cancer with radiation: the registration of the MR spectroscopy imaging to ultrasound/CT images and its implementation in treatment planning

PURPOSE

Dose-escalation to intraprostatic tumor deposits detected by magnetic resonance spectroscopy (MRS) is an example of tumor-targeted radiation therapy. Because treatment planning for prostate brachytherapy is performed based on ultrasound (US)/computed tomography (CT) images, a sine qua non of this technique is the ability to map MRS-positive volumes (obtained in a gland deformed by the endorectal balloon coil) to the US/CT images. An empirical algorithm designed to perform this function, and its validation, are described.

METHODS AND MATERIALS

Mathematically, the problem of mapping points between the MR and US/CT domains comes to: (a) ascertaining that the position of any point in the interior of the prostate is uniquely determined by the shape of the gland, and (b) finding an algorithm that describes this relationship. The image registration algorithm described here is based on the assumption that points within the gland maintain the same relative position with respect to both the axial contours of the prostate and the center of the prostate along the superior-inferior direction. Relative positions of MRS-positive voxels are calculated with this method in both MR and US/CT space. For a particular voxel in the MR space, one obtains first the z coordinate in the US/CT space, that is, along the superior-inferior direction. This determines the axial slice in the US/CT frame of reference where the other two coordinates (x, y) will be calculated. The validity of this algorithm was examined with the aid of a pelvic phantom built to simulate realistically the prostate and its surrounding bony and tissue structures and with CT scans of implanted patients obtained, at several weeks' intervals, as part of an edema-resolution study. Seventy-five "dummy" seeds were placed in the phantom, within the simulated prostate gland, in a quasi-regular pattern. The coordinates of these seeds were determined and thus served as markers of prostate deformation when an inflated rectal probe was introduced in the phantom. CT images of this phantom were taken for different volumes of the MR rectal probe and in each case the prostate outlines were contoured and seed coordinates calculated. Using these data, the predictions of the mapping algorithm could be directly verified.

RESULTS

Absolute values of the 3D-positional errors in this algorithm were 2.2 mm +/- 1.2 mm (average +/- SD). Only 6 of 75 seeds had positional displacement of 4 mm or more. Similar results were obtained in the patient analysis.

CONCLUSIONS

In comparison to the MRS voxel size (6.25 x 6.25 x 3.0 mm3), the present algorithm achieves the desired clinical accuracy. As well, with this 3D algorithm seed positions are reconstructed with an uncertainty that, along the z direction, is less than half the thickness of the typical US slice (0.5 cm).

Comparison of MRI pulse sequences in defining prostate volume after permanent implantation

PURPOSE

To determine the relative value of three MRI pulse sequences in defining the prostate volume after permanent implantation.

METHODS AND MATERIALS

A total of 45 patients who received a permanent 125I implant were studied. Two weeks after implantation, an axial CT scan (2 mm thickness) and T1-weighted, T1-weighted fat saturation, and T2-weighted axial MRI (3-mm) studies were obtained. The prostate volumes were compared with the initial ultrasound planning volumes, and subsequently the CT, T1-weighted, and T1-weighted fat saturation MRI volumes were compared with the T2-weighted volumes. Discrepancies in volume were evaluated by visual inspection of the registered axial images and the registration of axial volumes on the sagittal T2-weighted volumes. In a limited set of patients, pre- and postimplant CT and T2-weighted MRI studies were available for comparison to determine whether prostate volume changes after implant were dependent on the imaging modality.

RESULTS

T1-weighted and T1-weighted fat saturation MRI and CT prostate volumes were consistently larger than the T2-weighted MRI prostate volumes, with a volume on average 1.33 (SD 0.24) times the T2-weighted volume. This discrepancy was due to the superiority of T2-weighted MRI for prostate definition at the following critical interfaces: membranous urethra, apex, and anterior base-bladder and posterior base-seminal vesicle interfaces. The differences in prostate definition in the anterior base region suggest that the commonly reported underdose may be due to overestimation of the prostate in this region by CT. The consistent difference in volumes suggests that the degree of swelling observed after implantation is in part a function of the imaging modality. In patients with pre- and postimplant CT and T2-weighted MRI images, swelling on the T2-weighted images was 1.1 times baseline and on CT was 1.3 times baseline, confirming the imaging modality dependence of prostate swelling.

CONCLUSION

Postimplant T2-weighted MRI images provided superior prostate definition in all critical regions of the prostate compared with CT and the other MRI sequences tested. In addition to defining an optimal technique, these findings call two prior observations into question. Under dosing at the anterior base region may be overestimated because of poor definition of the prostate-bladder muscle interface. The swelling observed after implantation was lower on T2-weighted images as well, suggesting that a fraction of postimplant swelling is a function of the imaging modality. These findings have implications for preimplant planning and postimplant evaluation. As implant planning techniques become more conformal, and registration methods become more efficient, T2-weighted MRI after implantation will improve the accuracy of postimplant dosimetry.

Inverse planning for functional image-guided intensity-modulated radiation therapy

Radiation therapy is an image-guided process whose success critically depends on the imaging modality used for treatment planning and the level of integration of the available imaging information. In this work, we establish a dose optimization framework for incorporating metabolic information from functional imaging modalities into the intensity-modulated radiation therapy (IMRT) inverse planning process and to demonstrate the technical feasibility of planning deliberately non-uniform dose distributions in accordance with functional imaging data. For this purpose, a metabolic map from functional images is discretized into a number of abnormality levels (ALs) and then fused with CT images. To escalate dose to the metabolically abnormal regions, we assume, for a given spatial point, a linear relation between the AL and the prescribed dose. But the formalism developed here is independent of the assumption and any other relation between AL and prescription is applicable. For a given AL and prescription relation, it is only necessary to prescribe the dose to the lowest AL in the target and the desired doses to other regions with higher AL values are scaled accordingly. To accomplish differential sparing of a sensitive structure when its functional importance (FI) distribution is known, we individualize the tolerance doses of the voxels within the structure according to their Fl levels. An iterative inverse planning algorithm in voxel domain is used to optimize the system with in homogeneous dose prescription. To model intra-structural trade-off, a mechanism is introduced through the use of voxel-dependent weighting factors, in addition to the conventional structure specific weighting factors which model the inter-structural trade-off. The system is used to plan a phantom case with a few hypothetical functional distributions and a brain tumour treatment with incorporation of magnetic resonance spectroscopic imaging data. The results indicated that it is technically feasible to produce deliberately non-uniform dose distributions according to the functional imaging requirements. Integration of functional imaging information into radiation therapy dose optimization allows for consideration of patient-specific biologic information and provides a significant opportunity to truly individualize radiation treatment. This should enhance our capability to safely and intelligently escalate dose and lays the technical foundation for future clinical studies of the efficacy of functional imaging-guided IMRT.

Combined magnetic resonance imaging and spectroscopic imaging approach to molecular imaging of prostate cancer

Magnetic resonance spectroscopic imaging (MRSI) provides a noninvasive method of detecting small molecular markers (historically the metabolites choline and citrate) within the cytosol and extracellular spaces of the prostate, and is performed in conjunction with high-resolution anatomic imaging. Recent studies in pre-prostatectomy patients have indicated that the metabolic information provided by MRSI combined with the anatomical information provided by MRI can significantly improve the assessment of cancer location and extent within the prostate, extracapsular spread, and cancer aggressiveness. Additionally, pre- and post-therapy studies have demonstrated the potential of MRI/MRSI to provide a direct measure of the presence and spatial extent of prostate cancer after therapy, a measure of the time course of response, and information concerning the mechanism of therapeutic response. In addition to detecting metabolic biomarkers of disease behavior and therapeutic response, MRI/MRSI guidance can improve tissue selection for ex vivo analysis. High-resolution magic angle spinning ((1)H HR-MAS) spectroscopy provides a full chemical analysis of MRI/MRSI-targeted tissues prior to pathologic and immunohistochemical analyses of the same tissue. Preliminary (1)H HR-MAS spectroscopy studies have already identified unique spectral patterns for healthy glandular and stromal tissues and prostate cancer, determined the composition of the composite in vivo choline peak, and identified the polyamine spermine as a new metabolic marker of prostate cancer. The addition of imaging sequences that provide other functional information within the same exam (dynamic contrast uptake imaging and diffusion-weighted imaging) have also demonstrated the potential to further increase the accuracy of prostate cancer detection and characterization.

Intraoperative planning and evaluation of permanent prostate brachytherapy: report of the American Brachytherapy Society

PURPOSE

The preplanned technique used for permanent prostate brachytherapy has limitations that may be overcome by intraoperative planning. The goal of the American Brachytherapy Society (ABS) project was to assess the current intraoperative planning process and explore the potential for improvement in intraoperative treatment planning (ITP).

METHODS AND MATERIALS

Members of the ABS with expertise in ITP performed a literature review, reviewed their clinical experience with ITP, and explored the potential for improving the technique.

RESULTS

The ABS proposes the following terminology in regard to prostate planning process: *Preplanning--Creation of a plan a few days or weeks before the implant procedure. *Intraoperative planning--Treatment planning in the operating room (OR): the patient and transrectal ultrasound probe are not moved between the volume study and the seed insertion procedure. * Intraoperative preplanning--Creation of a plan in the OR just before the implant procedure, with immediate execution of the plan. *Interactive planning--Stepwise refinement of the treatment plan using computerized dose calculations derived from image-based needle position feedback. *Dynamic dose calculation--Constant updating of dose distribution calculations using continuous deposited seed position feedback. Both intraoperative preplanning and interactive planning are currently feasible and commercially available and may help to overcome many of the limitations of the preplanning technique. Dosimetric feedback based on imaged needle positions can be used to modify the ITP. However, the dynamic changes in prostate size and shape and in seed position that occur during the implant are not yet quantifiable with current technology, and ITP does not obviate the need for postimplant dosimetric analysis. The major current limitation of ITP is the inability to localize the seeds in relation to the prostate. Dynamic dose calculation can become a reality once these issues are solved. Future advances can be expected in methods of enhancing seed identification, in imaging techniques, and in the development of better source delivery systems. Additionally, ITP should be correlated with outcome studies, using dosimetric, toxicity, and efficacy endpoints.

CONCLUSION

ITP addresses many of the limitations of current permanent prostate brachytherapy and has some advantages over the preplanned technique. Further technologic advancement will be needed to achieve dynamic real-time calculation of dose distribution from implanted sources, with constant updating to allow modification of subsequent seed placement and consistent, ideal dose distribution within the target volume.

In vivo proton (H1) magnetic resonance spectroscopy for cervical carcinoma

Proton magnetic resonance spectroscopy (MRS) may be a useful tool in both the initial diagnosis of cervical carcinoma and the subsequent surveillance after radiation therapy, particularly when other standard diagnostic methods are inconclusive. Single voxel magnetic resonance (MR) spectral data were acquired from 8 normal volunteers, 16 patients with cervical cancer before radiation therapy, and 18 patients with cervical cancer after radiation therapy using an external pelvic coil at a 1.5-T on a Signa system. The presence or absence of various resonances within each spectrum was evaluated for similarities within each patient group and for spectral differences between groups. Resonances corresponding to lipid and creatine dominated the spectrum for the eight normal volunteers without detection of a choline resonance. Spectra from 16 pretreatment patients with biopsy-proven cervical cancer revealed strong resonances at a chemical shift of 3.25 ppm corresponding to choline. Data acquired from the 18 posttreatment setting studies was variable, but often correlated well with the clinical findings. Biopsy confirmation was obtained in seven patients. H1 MRS of the cervix using a noninvasive pelvic coil consistently demonstrates reproducible spectral differences between normal and neoplastic cervical tissue in vivo. However, signal is still poor for minimal disease recurrence. Further study is needed at intervals before, during, and after definitive irradiation with biopsy confirmation to validate the accuracy of MRS in distinguishing persistence or recurrence of disease from necrosis and fibrosis.

Magnetic resonance spectroscopy of the malignant prostate gland after radiotherapy: a histopathologic study of diagnostic validity

PURPOSE

Accurate spatial representation of tumor clearance after conformal radiotherapy is an endpoint of clinical importance. Magnetic resonance spectroscopy (MRS) can diagnose malignancy in the untreated prostate gland through measurements of cellular metabolites. In this study we sought to describe spectral metabolic changes in prostatic tissue after radiotherapy and validate a multivariate analytic strategy (based on MRS) that could identify viable tumor.

METHODS AND MATERIALS

Transrectal ultrasound-guided prostate biopsies from 35 patients were obtained 18-36 months after external beam radiotherapy. One hundred sixteen tissue specimens were subjected to 1H MRS, submitted to histopathology, and analyzed for correlation with a multivariate strategy specifically developed for biomedical spectra.

RESULTS

The sensitivity and specificity of MRS in identifying a malignant biopsy were 88.9% and 92% respectively, with an overall classification accuracy of 91.4%. The diagnostic spectral regions identified by our algorithm included those due to choline, creatine, glutamine, and lipid. Citrate, an important discriminating resonance in the untreated prostate gland, was invisible in all spectra, regardless of histology.

CONCLUSIONS

Although the spectral features of prostate tissue markedly change after radiotherapy, MRS combined with multivariate methods of analysis can accurately identify histologically malignant biopsies. MRS shows promise as a modality that could integrate three-dimensional measures of tumor response.

Radioimmunoguided imaging of prostate cancer foci with histopathological correlation

PURPOSE

We have previously presented a technique that fuses ProstaScint and pelvic CT images for the purpose of designing brachytherapy that targets areas at high risk for treatment failure. We now correlate areas of increased intensity seen on ProstaScint-CT fusion images to biopsy results in a series of 7 patients to evaluate the accuracy of this technique in localizing intraprostatic disease.

METHODS AND MATERIALS

The 7 patients included in this study were evaluated between June 1998 and March 29, 1999 at Metrohealth Medical Center and University Hospitals of Cleveland in Cleveland, Ohio. ProstaScint and CT scans of each patient were obtained before transperineal biopsy and seed implantation. Each patient's prostate gland was biopsied at 12 separate sites determined independently of Prostascint-CT scan results.

RESULTS

When correlated with biopsy results, our method yielded an overall accuracy of 80%: with a sensitivity of 79%, a specificity of 80%, a positive predictive value of 68%, and a negative predictive value of 88%.

CONCLUSION

The image fusion of the pelvic CT scan and ProstaScint scan helped identify foci of adenocarcinoma within the prostate that correlated well with biopsy results. These data may be useful to escalate doses in regions containing tumor by either high-dose rate or low-dose rate brachytherapy, as well as by external beam techniques such as intensity modulated radiotherapy (IMRT).