Francesca Botta

Quantitative Analysis of 90Y Bremsstrahlung SPECT-CT Images for Application to 3D Patient-Specific Dosimetry

AIM The aim of this study was to evaluate the accuracy of the activity quantification of single-photon emission computed tomography/computed tomography (SPECT-CT) (90)Y-Bremsstrahlung images and to validate the S-voxel method. METHODS An anthropomorphic torso phantom with radioactive inserts ((90)Y) was acquired by SPECT-CT. Constant calibration factors (cps/MBq) for the quantification were evaluated, considering different volume, shape, position inside the phantom, activity concentration and background, and distance from detectors. S-voxel values (EGSnrc) were implemented in MATLAB R0086 USA software. Dose comparisons between S-voxel and the conventional Medical Internal Radiation Dose method were repeated in a group of 11 patients administered with (90)Y-DOTATATE. RESULTS Using the appropriate calibration factors to recover the volume variability, the error about the measurement repeatability and the activity variation was within 4%. The variability of activity quantification, depending on the position in the phantom, detector distance, and background, was <10%, <5%, and <10%, respectively. The absorbed-dose values calculated by OLINDA were in agreement with the mean dose values obtained by the S-voxel method (difference, <10%). CONCLUSIONS The results confirm that, with the hybrid SPECT-CT system, quantitative analysis of SPECT (90)Y-Bremsstrahlung images and the generation of three-dimensional dose distributions are feasible. The improved analysis of Bremsstrahlung images could have a notable clinical impact, allowing to address the dosimetric verification to patients during the course of therapy.

Radioembolization of Hepatic Lesions from a Radiobiology and Dosimetric Perspective

Radioembolization (RE) of liver cancer with 90Y-microspheres has been applied in the last two decades with notable responses and acceptable toxicity. Two types of microspheres are available, glass and resin, the main difference being the activity/sphere. Generally, administered activities are established by empirical methods and differ for the two types. Treatment planning based on dosimetry is a prerogative of few centers, but has notably gained interest, with evidence of predictive power of dosimetry on toxicity, lesion response, and overall survival (OS). Radiobiological correlations between absorbed doses and toxicity to organs at risk, and tumor response, have been obtained in many clinical studies. Dosimetry methods have evolved from the macroscopic approach at the organ level to voxel analysis, providing absorbed dose spatial distributions and dose–volume histograms (DVH). The well-known effects of the external beam radiation therapy (EBRT), such as the volume effect, underlying disease influence, cumulative damage in parallel organs, and different tolerability of re-treatment, have been observed also in RE, identifying in EBRT a foremost reference to compare with. The radiobiological models – normal tissue complication probability and tumor control probability – and/or the style (DVH concepts) used in EBRT are introduced in RE. Moreover, attention has been paid to the intrinsic different activity distribution of resin and glass spheres at the microscopic scale, with dosimetric and radiobiological consequences. Dedicated studies and mathematical models have developed this issue and explain some clinical evidences, e.g., the shift of dose to higher toxicity thresholds using glass as compared to resin spheres. This paper offers a comprehensive review of the literature incident to dosimetry and radiobiological issues in RE, with the aim to summarize the results and to identify the most useful methods and information that should accompany future studies.

A free database of radionuclide voxel S values for the dosimetry of nonuniform activity distributions.

The increasing availability of SPECT/CT devices with advanced technology offers the opportunity for the accurate assessment of the radiation dose to the biological target volume during radionuclide therapy. Voxel dosimetry can be performed employing direct Monte Carlo radiation transport simulations, based on both morphological and functional images of the patient. On the other hand, for voxel dosimetry calculations the voxel S value method can be considered an easier approach than patient-specific Monte Carlo simulations, ensuring a good dosimetric accuracy at least for anatomic regions which are characterized by uniform density tissue. However, this approach has been limited because of the lack of tabulated S values for different voxel dimensions and radionuclides. The aim of this work is to provide a free dataset of values which can be used for voxel dosimetry in targeted radionuclide studies. Seven different radionuclides (89Sr, 90Y, 131I, 153Sm, 177Lu, 186Re, 188Re), and 13 different voxel sizes (2.21, 2.33, 2.4, 3, 3.59, 3.9, 4, 4.42, 4.8, 5, 6, 6.8 and 9.28 mm) are considered. Voxel S values are calculated performing simulations of monochromatic photon and electron sources in two different homogeneous tissues (soft tissue and bone) with DOSXYZnrc code, and weighting the contributions on the basis of the radionuclide emission spectra. The outcomes are validated by comparison with Monte Carlo simulations obtained with other codes (PENELOPE and MCNP4c) performing direct simulation of the radionuclide emission spectra. The differences among the different Monte Carlo codes are of the order of a few per cent when considering the source voxel and the bremsstrahlung tail, whereas the highest differences are observed at a distance close to the maximum continuous slowing down approximation range of electrons. These discrepancies would negligibly affect dosimetric assessments. The dataset of voxel S values can be freely downloaded from the website www.medphys.it.

Calculation of electron and isotopes dose point kernels with fluka Monte Carlo code for dosimetry in nuclear medicine therapy

PURPOSE The calculation of patient-specific dose distribution can be achieved by Monte Carlo simulations or by analytical methods. In this study, FLUKA Monte Carlo code has been considered for use in nuclear medicine dosimetry. Up to now, FLUKA has mainly been dedicated to other fields, namely high energy physics, radiation protection, and hadrontherapy. When first employing a Monte Carlo code for nuclear medicine dosimetry, its results concerning electron transport at energies typical of nuclear medicine applications need to be verified. This is commonly achieved by means of calculation of a representative parameter and comparison with reference data. Dose point kernel (DPK), quantifying the energy deposition all around a point isotropic source, is often the one. METHODS FLUKA DPKS have been calculated in both water and compact bone for monoenergetic electrons (10-3 MeV) and for beta emitting isotopes commonly used for therapy (89Sr, 90Y, 131I 153Sm, 177Lu, 186Re, and 188Re). Point isotropic sources have been simulated at the center of a water (bone) sphere, and deposed energy has been tallied in concentric shells. FLUKA outcomes have been compared to PENELOPE v.2008 results, calculated in this study as well. Moreover, in case of monoenergetic electrons in water, comparison with the data from the literature (ETRAN, GEANT4, MCNPX) has been done. Maximum percentage differences within 0.8.RCSDA and 0.9.RCSDA for monoenergetic electrons (RCSDA being the continuous slowing down approximation range) and within 0.8.X90 and 0.9.X90 for isotopes (X90 being the radius of the sphere in which 90% of the emitted energy is absorbed) have been computed, together with the average percentage difference within 0.9.RCSDA and 0.9.X90 for electrons and isotopes, respectively. RESULTS Concerning monoenergetic electrons, within 0.8.RCSDA (where 90%-97% of the particle energy is deposed), FLUKA and PENELOPE agree mostly within 7%, except for 10 and 20 keV electrons (12% in water, 8.3% in bone). The discrepancies between FLUKA and the other codes are of the same order of magnitude than those observed when comparing the other codes among them, which can be referred to the different simulation algorithms. When considering the beta spectra, discrepancies notably reduce: within 0.9.X90, FLUKA and PENELOPE differ for less than 1% in water and less than 2% in bone with any of the isotopes here considered. Complete data of FLUKA DPKS are given as Supplementary Material as a tool to perform dosimetry by analytical point kernel convolution. CONCLUSIONS FLUKA provides reliable results when transporting electrons in the low energy range, proving to be an adequate tool for nuclear medicine dosimetry.

Correlation between egfr expression and accelerated proliferation during radiotherapy of head and neck squamous cell carcinoma

PurposeTo investigate the correlation between the expression of Epidermal Growth Factor receptor (EGFr) and the reduction of the effective doubling time (TD) during radiotherapy treatment and also to determine the dose per fraction to be taken into account when the overall treatment time (OTT) is reduced in accelerated radiotherapy of head and neck squamous cell carcinoma (HNSCC).MethodsA survey of the published papers comparing 3-years of local regional control rate (LCR) for a total of 2162 patients treated with conventional and accelerated radiotherapy and with a pretreatment assessment of EGFr expression, was made. Different values of TD were obtained by a model incorporating the overall time corrected biologically effective dose (BED) and a 3-year clinical LCR for high and low EGFr groups of patients (HEGFr and LEGFr), respectively. By obtaining the TD from the above analysis and the sub-sites’ potential doubling time (Tpot) from flow cytometry and immunohistochemical methods, we were able to estimate the average TD for each sub-site included in the analysis. Moreover, the dose that would be required to offset the modified proliferation occurring in one day (Dprolif), was estimated.ResultsThe averages of TD were 77 (27-90)95% days in LEGFr and 8.8 (7.3-11.0)95% days in HEGFr, if an onset of accelerated proliferation TK at day 21 was assumed. The correspondent HEGFr sub-sites’ TD were 5.9 (6.6), 5.9 (6.6), 4.6 (6.1), 14.3 (12.9) days, with respect to literature immunohistochemical (flow cytometry) data of Tpot for Oral-Cavity, Oro-pharynx, Hypo-pharynx, and Larynx respectively. The Dprolif for the HEGFr groups were 0.33 (0.29), 0.33 (0.29), 0.42 (0.31), 0.14 (0.15) Gy/day if α = 0.3 Gy-1 and α/β = 10 Gy were assumed.ConclusionsA higher expression of the EGFr leads to enhanced proliferation. This study allowed to quantify the extent of the effect which EGFr expression has in terms of reduced TD and Dprolif for each head and neck sub-site.

Is [18F] fluorodeoxyglucose uptake by the primary tumor a prognostic factor in breast cancer?

BACKGROUND We retrospectively investigated (18)F-FDG uptake by the primary breast tumor as a predictor for relapse and survival. PATIENTS AND METHODS We studied 203 patients with cT1-T3N0 breast cancer. Standardized uptake value (SUVmax), was measured on the primary tumor. After a median follow-up of 68 months (range 22-80), the relation between SUVmax and tumor factors, disease free-survival (DFS) and overall survival (OS) was investigated. RESULTS In the PET-positive patients, the median FDG uptake by the tumor was 4.7. FDG uptake was significantly related to tumor size, number of involved axillary nodes, grade, negative ER, high Ki-67 and HER2 overexpression. No distant metastases or deaths occurred in the PET-negative group. Five-year DFS was 97% and 83%, respectively in the PET-negative and PET-positive groups (P = 0.096). At univariate analysis, DFS was significantly lower in patients with SUVmax >4.7 compared to the patients with negative PET (P = 0.042), but not to the patients with SUVmax ≤4.7 (P = 0.106). At multivariable analysis, among PET-positive patients, SUVmax was not an independent prognostic factor for DFS (HR(>4.7 vs ≤4.7): 1.02 (95% CI 0.45-2.31)). Five-year OS was 100% and 93%, respectively, in the PET-negative and PET-positive groups (P = 0.126). CONCLUSION FDG uptake by the primary lesion was significantly associated with several prognostic variables, but it was not an independent prognostic factor.

Differences in 3D dose distributions due to calculation method of voxel S-values and the influence of image blurring in SPECT.

This study compares 3D dose distributions obtained with voxel S values (VSVs) for soft tissue, calculated by several methods at their current state-of-the-art, varying the degree of image blurring. The methods were: 1) convolution of Dose Point Kernel (DPK) for water, using a scaling factor method; 2) an analytical model (AM), fitting the deposited energy as a function of the source-target distance; 3) a rescaling method (RSM) based on a set of high-resolution VSVs for each isotope; 4) local energy deposition (LED). VSVs calculated by direct Monte Carlo simulations were assumed as reference. Dose distributions were calculated considering spheroidal clusters with various sizes (251, 1237 and 4139 voxels of 3 mm size), uniformly filled with (131)I, (177)Lu, (188)Re or (90)Y. The activity distributions were blurred with Gaussian filters of various widths (6, 8 and 12 mm). Moreover, 3D-dosimetry was performed for 10 treatments with (90)Y derivatives. Cumulative Dose Volume Histograms (cDVHs) were compared, studying the differences in D95%, D50% or Dmax (ΔD95%, ΔD50% and ΔDmax) and dose profiles.For unblurred spheroidal clusters, ΔD95%, ΔD50% and ΔDmax were mostly within some percents, slightly higher for (177)Lu with DPK (8%) and RSM (12%) and considerably higher for LED (ΔD95% up to 59%). Increasing the blurring, differences decreased and also LED yielded very similar results, but D95% and D50% underestimations between 30-60% and 15-50%, respectively (with respect to 3D-dosimetry with unblurred distributions), were evidenced. Also for clinical images (affected by blurring as well), cDVHs differences for most methods were within few percents, except for slightly higher differences with LED, and almost systematic for dose profiles with DPK (-1.2%), AM (-3.0%) and RSM (4.5%), whereas showed an oscillating trend with LED.The major concern for 3D-dosimetry on clinical SPECT images is more strongly represented by image blurring than by differences among the VSVs calculation methods. For volume sizes about 2-fold the spatial resolution, D95% and D50% underestimations up to about 60 and 50% could result, so the usefulness of 3D-dosimetry is highly questionable for small tumors, unless adequate corrections for partial volume effects are adopted.

Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images

Patient-specific absorbed dose calculation for nuclear medicine therapy is a topic of increasing interest. 3D dosimetry at the voxel level is one of the major improvements for the development of more accurate calculation techniques, as compared to the standard dosimetry at the organ level. This study aims to use the FLUKA Monte Carlo code to perform patient-specific 3D dosimetry through direct Monte Carlo simulation on PET-CT and SPECT-CT images. To this aim, dedicated routines were developed in the FLUKA environment. Two sets of simulations were performed on model and phantom images. Firstly, the correct handling of PET and SPECT images was tested under the assumption of homogeneous water medium by comparing FLUKA results with those obtained with the voxel kernel convolution method and with other Monte Carlo-based tools developed to the same purpose (the EGS-based 3D-RD software and the MCNP5-based MCID). Afterwards, the correct integration of the PET/SPECT and CT information was tested, performing direct simulations on PET/CT images for both homogeneous (water) and non-homogeneous (water with air, lung and bone inserts) phantoms. Comparison was performed with the other Monte Carlo tools performing direct simulation as well. The absorbed dose maps were compared at the voxel level. In the case of homogeneous water, by simulating 10(8) primary particles a 2% average difference with respect to the kernel convolution method was achieved; such difference was lower than the statistical uncertainty affecting the FLUKA results. The agreement with the other tools was within 3–4%, partially ascribable to the differences among the simulation algorithms. Including the CT-based density map, the average difference was always within 4% irrespective of the medium (water, air, bone), except for a maximum 6% value when comparing FLUKA and 3D-RD in air. The results confirmed that the routines were properly developed, opening the way for the use of FLUKA for patient-specific, image-based dosimetry in nuclear medicine.

Therapeutic schemes in 177Lu and 90Y-PRRT: radiobiological considerations.

BACKGROUND The purpose of this work is to implement a radiobiological model to compare different treatment schedules for Peptide Receptor Radionuclide Therapy (PRRT) with 177Lu and 90Y. The principal radiobiological quantities were studied as a function of radionuclides, fractionation schemes, activity distribution in kidneys and tumor radiosensitivity. METHODS Clinical data were used to derive representative absorbed doses for several treatment schemes for 177Lu-PRRT and for 90Y-PRRT and considered as input data for the radiobiological model. Both uniform and non-uniform activity distributions were considered for kidneys and cortex; for tumors a possible uptake reduction after each cycle and inter-patient radiosensitivity variability were investigated. Normal-Tissue-Complication-Probability (NTCP) and Tumor-Control-Probability (TCP) were evaluated. RESULTS Hyper-cycling has a limited advantage in terms of BED reduction on kidneys for 177Lu, while for 90Y the effect is sizable and helps in reducing the NTCP. For all 177Lu-schemes the renal toxicity risk is negligible while for some 90Y-schemes the NTCP is not null. In case of tumor uptake reduction with cycles the treatment efficacy is reduced with a BED loss up to 46%. The TCP decreases when assuming normally-distributed tumor radiosensitivity values. CONCLUSIONS This paper discusses how the combination of dosimetry and radiobiological modeling may help in exploring the link between the treatment schedule and the potential clinical outcome. The results highlight the capability of model to reproduce the available clinical data and provide useful qualitative information. Further investigation on dose distribution and dose uptake reduction with accurate clinical data is needed to progress in this field.

Planning Combined Treatments of External Beam Radiation Therapy and Molecular Radiotherapy

Molecular radiotherapy (MRT) with radiolabeled molecules has being constantly evolving, leading to notable results in cancer treatment. In some cases, the absorbed doses delivered to tumors by MRT are sufficient to obtain complete responses; in other cases, instead, to be effective, MRT needs to be combined with other therapeutic approaches. Recently, several studies proposed the combination of MRT with external beam radiation therapy (EBRT). Some describe the theoretical basis within radiobiological models, others report the results of clinical phase I-II studies aimed to assess the feasibility and tolerability. The latter includes the treatment of various tumors, such as meningiomas, paragangliomas, non-Hodgkins lymphomas, bone, brain, hepatic, and breast lesions. The underlying principle of combined MRT and EBRT is the possibility of exploiting the full potential of each modality, given the different organs at risk. Target tissues can indeed receive a higher irradiation, while respecting the threshold limits of more than one critical tissue. Nevertheless, clinical trials are empirical and optimization is still a theoretical issue. This article describes the state of the art of combined MRT and EBRT regarding the rationale and the results of clinical studies, with special focus on the possibility of treatment improvement.