This single-center retrospective study was performed in accordance with the national ethical guidelines and regulations. The national Ethics Committee has approved the methods of this study. MRI and CT images of patients were used in this study without any intervention in the diagnostic or treatment procedures. In addition, gathering the informed consent was waived because of the retrospective nature of the study.
Imaging data (CT and MRI scans) of 20 male patients having local prostate cancer who had no previous surgery, hormone therapy (AST or ADT), and prostate radiation therapy with at least one non-high risk intraprostatic lesion (IL) in stages of T1 to T3a, were used in this study. Patients’ ages ranged from 54 to 85 years, with a mean age of 69.4.
The CT (Matrix size: 512*512; Slice thickness: 3 to 5mm), T2w-MRI (fast spin echo pulse sequence with TE: 80 ms and TR: 7800ms), diffusion-weighted MRI (echo planar imaging with TE: 88ms, and TR: 4600ms) and apparent diffusion coefficient (fast spin echo pulse sequence with TE: 100ms and TR: 3000ms) images were taken using a Siemens 16-slice Emotion CT and a 1.5 Tesla Avanto MRI machine (Siemens Healthcare GmbH, Germany). The patients were placed in supine positions for both imaging procedures. Diffusion-weighted images (DWI) were gathered with three signals per image with a scattering-sensitive gradient in three orthogonal planes and b-values of 0, 250, 500, and 1000 seconds per square. DW-MRI images have a resolution of 1.64×1.64×3 mm and a FOV of 210×210 mm, a matrix size of 128×128 pixels, and a NEX (number of excitation) parameter equal to four. Apparent diffusion coefficient (ADC) maps (images) were automatically calculated from DW-MRI images.
CT and MRI images were combined using a rigid registration algorithm in the treatment planning software (TPS) based on bone landmarks, gold markers implanted in the prostate and skin surfaces, and then verified by a specialized physician. The MRI and CT registration were used to contour the lesion volumes inside the patients’ prostate and radiation-sensitive organs. It also allows for more precise target volume delineation in prostate cancer patients (31).
An in-house MATLAB program was developed to automatically identify lesion regions on ADC images based on ADC values. The MATLAB code is available in the “supplementary materials” section. Two types of predominant lesions were considered in the prostate, one related to lesions with ADC values lower than 750 mm2/s (PTV-1), and the other was related to lesions with ADC values higher than 750 and less than 1500 mm2 /s (PTV-2). The upper limit of ADC values in tumor tissues varies in different studies but usually was considered more than 1300 mm2/s (32). In this study, the apparent diffusion coefficient threshold for distinguishing tumor tissue from normal prostate tissue was 1500 mm2/s. The highest measurable value of ADC in MATLAB software was 5000 mm2/s. The ADC images were imported to the MATLAB program, and the voxels in a specific range of ADC values were determined. The related voxels for determining areas must at least have a minimum number (400 voxels) located next to each other so that the software can identify those areas separately. Considering the relationship between Gleason score (GS) and ADC cut-off value based on a study by Pepe et al. (33), different target volumes were identified within the prostate. Each patient's output DICOM RS file was then transferred to the treatment planning system to contour these new structures (prostatic lesions) on the CT images. Contouring of other organs at risk (OARs) was performed by an experienced radiation oncologist in the treatment planning system.
Prescription dose levels, except in dominant intraprostatic lesions (DILs), were taken from the Jereczek-Fossa et al. study (34). The upper limit for the prescribed dose for DIL was considered 70 Gy in 27 fractions for high risk DILs. This hypofractionated dose escalation was used and evaluated in many studies (35–42). According to a study by Onjukka et al. (43), this hypofractionated dose was equivalent to 86 Gy in 37 sessions (used in the study of Uzan et al. (44). A prescribed dose for DILs with lower risks was considered 66 Gy. The clinical target volume of the base of the seminal vesicles was considered CTV53Gy. The planning target volume for this target (PTV53Gy) was formed by adding eight millimeter isotropic margins to CTV53Gy in order to account for patient and equipment placement errors. The whole prostate volume (except the DILs) was considered CTV60; similarly, PTV60Gy was formed by adding 5 mm margins to the CTV60Gy. The margin was reduced to zero in the posterior region, where the target volume overlaps with the rectum. Two millimeters margins were added to CTV66Gy and CTV70Gy to create PTV66Gy and PTV70Gy without extending beyond CTV60Gy or overlapping with the rectum, bladder, and urethra due to the uncertainty in defining the DILs (44). Furthermore, planning at risk volumes (PRV) were created for high-risk organs, including the rectum, urethra, and bladder, with margins of two millimeters. The dose escalated DIL regions with the whole prostate were presented in figure 1. Furthermore, the procedure used to automatically contour the intraprostatic DILs is illustrated in figures 1-a and 1-b.
The IMRT plans were designed with Eclipse software (version 11, Varian Corporation, USA) for each patient. IMRT plan with nine coplanar fields in gantry angles of 0, 30, 60, 105, 140, 220, 260, 300, and 330 was designed to irradiate PTVs with prescribed doses. All the plans were interactively optimized based on our institutional planning protocol derived from a previous study by Pollak et al. (45). The planning optimization objectives are presented in table 1. An experienced physicist evaluated all the treatment plans to ensure compliance with reported dose constraints (46).
Table 1- The planning optimization objectives used for prostate dose painting IMRT
Structures
|
Objectives
|
Bladder
|
V40.8Gy<50%
|
V48.6Gy<25%
|
V60Gy<5%
|
Max dose < 65Gy
|
Rectum
|
V40.8Gy<50%
|
V48.6Gy<35%
|
Max dose < 65Gy
|
V60Gy<3%
|
Femoral heads
|
Max dose < 40Gy
|
Bowel
|
V50Gy<17cc
|
Max dose < 60Gy
|
PTVs
|
V98%>98%
|
V105%<2%
|
*VxGy represents the percentage of the structure volume received at least x Gy.
*Vx% represents the percentage of the structure volume received at least x% of the prescribed dose.
After treatment planning optimization and final dose calculation, dose volume histograms (DVHs) of the CTVs for each patient were entered in BioSuite software (47). The tumor control probability (TCP) values were calculated using the Poisson model (48), based on radiobiological model parameters proposed by Deb and Fielding (49).
We introduced two indices of effectiveness (IOE) for evaluating IMRT dose painting plan dose distribution. One IOE can evaluate the conformity of CTVs, IOE(C), and another IOE can assess the overall dose distribution homogeneity of target volumes, IOE(H). The previous equation proposed by Park et al. (28) has relative volume coefficients, and these coefficients were included in our IOE equations accounting for the effect of each target (DIL) on the overall value of IOE. Furthermore, cell density values obtained from ADC maps were used in the IOE equations. The cell density is a measure of the clonogenicity level for each of the tumor volumes. The equations of IOE(H) and IOE(C) are as follows:
