Personalized Dosimetry in the Context of Radioiodine Therapy for Differentiated Thyroid Cancer.

The most frequent thyroid cancer is Differentiated Thyroid Cancer (DTC) representing more than 95% of cases. A suitable choice for the treatment of DTC is the systemic administration of 131-sodium or potassium iodide. It is an effective tool used for the irradiation of thyroid remnants, microscopic DTC, other nonresectable or incompletely resectable DTC, or all the cited purposes. Dosimetry represents a valid tool that permits a tailored therapy to be obtained, sparing healthy tissue and so minimizing potential damages to at-risk organs. Absorbed dose represents a reliable indicator of biological response due to its correlation to tissue irradiation effects. The present paper aims to focus attention on iodine therapy for DTC treatment and has developed due to the urgent need for standardization in procedures, since no unique approaches are available. This review aims to summarize new proposals for a dosimetry-based therapy and so explore new alternatives that could provide the possibility to achieve more tailored therapies, minimizing the possible side effects of radioiodine therapy for Differentiated Thyroid Cancer.

Yttrium-90 Radioembolization of Colorectal Cancer Liver Metastases: A Systematic Review of Clinical Evidence.

The present work aims to review the clinical evidence of survival outcome after treatment of colorectal cancer liver metastases using yttrium-90 radioembolization, existence, and relevance of clinical, imaging, and genomic predictors of treatment efficacy and the amount of administered activity. Publications listed in PubMed between July 2016 and May 2021 were screened. Predictors of overall survival were reported and distinguished in clinical, imaging, and genomic variables. Administered activity is reported as median and mean value; overall survival is presented as a median value from the treatment. Fourteen papers resulted to be eligible for this systematic review, 11 retrospective, and 3 prospective studies. Ten studies reported administered activity data, with an average mean administered activity of 1.63 GBq and an average median administered activity of 1.53 GBq. Many clinical, imaging and genomic variables have been identified or defined as predictors of efficacy, leading to the possibility of improvements in patient selection criteria. The overall survival resulted to be about 9 months. The clinical evidence on the application of radioembolization with yittrium-90 resumed in this work underlines the importance to analyze several clinical, imaging, and genomic variables to predict the outcome of the therapy. The overall survival has not improved significantly with respect to older studies. Further developments on treatment planning and patient selection could lead to better clinical outcomes.

Technical note: A wearable radiation measurement system for collection of patient-specific time-activity data in radiopharmaceutical therapy: system design and Monte Carlo simulation results.

PURPOSE: A high level of personalization in Molecular Radiotherapy (MRT) could bring advantages in terms of treatment effectiveness and toxicity reduction. Individual organ-level dosimetry is crucial to describe the radiopharmaceutical biodistribution expressed by the patient, to estimate absorbed doses to normal organs and target tissue(s). This paper presents a proof-of-concept Monte Carlo simulation study of "WIDMApp" (Wearable Individual Dose Monitoring Apparatus), a multi-channel radiation detector and data processing system for in vivo patient measurement and collection of radiopharmaceutical biokinetic data (i.e., time-activity data). Potentially, such a system can increase the amount of such data that can be collected while reducing the need to derive it via nuclear medicine imaging. METHODS: a male anthropomorphic MIRD phantom was used to simulate photons (i.e., gamma-rays) propagation in a patient undergoing a 131 I thyroid treatment. The administered activity was set to the amount usually administered for the treatment of differentiated carcinoma while its initial distribution in different organs was assigned following the ICRP indications for the 131 I biokinetics. Using this information, the simulation computes the Time-dependent Counts Curves (TCCs) that would have been measured by seven WIDMApp-like sensors placed and oriented to face each one of five emitting organs plus two thyroid lobes. A deconvolution algorithm was then applied on this simulated data set to reconstruct the Time-Activity Curve (TAC) of each organ. Deviations of the reconstructed TACs parameters from values used to generate them were studied as a function of the deconvolution algorithm initialization parameters and assuming non-Poisson fluctuation of the TCCs data points. RESULTS: This study demonstrates that it is possible, at least in the simple simulated scenario, to reconstruct the organ cumulated activity by measuring the time dependence of counts recorded by several detectors placed at selected positions on the patient's body. The ability to perform in vivo sampling more frequently than conventional biokinetic studies increases the number of time points and therefore the accuracy in TAC estimates. In this study, an accuracy on cumulated activity of 5% is obtained even with a 20% error on the TCC data points and a 50% error on the initial guess on the parameters of the deconvolution algorithm. CONCLUSIONS: the WIDMApp approach could provide an effective tool to characterize more accurately the radiopharmaceutical biokinetics in MRT patients, reducing the need of resources of nuclear medicine departments, such as technologist and scanner time, to perform individualized biokinetics studies. The relatively simple hardware for the approach proposed would allow its application to large numbers of patients. The results obtained justify development of an actual prototype system to characterize this technique under realistic conditions.

Dosimetric optimization of nuclear medicine therapy based on the Council Directive 2013/59/EURATOM and the Italian law N. 101/2020. Position paper and recommendations by the Italian National Associations of Medical Physics (AIFM) and Nuclear Medicine (AIMN).

This recommendation by the Italian Associations of Nuclear Medicine (AIMN) and Medical Physics (AIFM) focuses on the dosimetric optimization of Nuclear Medicine Therapy (NMT) as clearly requested by the article 56 of the EURATOM Directive 2013/59 and its consequent implementation in article 158 in the Italian Law n. 101/2020. However, this statement must deal with scientific and methodological limits that still exist and, above all, with the currently available limited resources. This paper addresses these specific issues. It distinguishes among many possible kinds of NMT. For each type, dosimetric optimization is recommended or considered optional, according to the general criteria adopted in any human choice, i.e. a check of technical feasibility first, followed by a cost/benefit argument. The classification of therapies as standardized or non-standardized is presented. This is based on the complexity of the type of pathology, on the variability of the treatment outcome, and on the risks involved. According to the present document, which was officially delivered to Italian Health Ministry as necessary interpretation of the law, a therapeutic team can, in science and consciousness, overcome the indications of posology, to optimize and tailoring a treatment with dosimetry, on the basis of published national or international data or guidelines, without need of an Ethics Committee approval. Data collected in this way will provide additional evidence about optimal dosimetric reference values. As conclusion, a formal appeal is made to the European and National regulatory agencies for pharmaceuticals to obtain the official acknowledgment of this principle.

Diagnostic Value of the Early Heart-to-Mediastinum Count Ratio in Cardiac 123I-mIBG Imaging for Parkinson's Disease.

BACKGROUND: Early diagnosis of Parkinson's disease (PD) is of primary importance. The delayed (3-4 h after injection) Iodine-123-Metaiodobenzylguanidine (123I-mIBG) scintigraphy has been proven to be effective in early differential diagnosis for Lewy body disease. But early imaging (15-30 min after injection) has only been marginally studied for its possible diagnostic role. In this prospective study, a threshold for the early Heart-to-Mediastinum (H/M) count ratio has been investigated, obtaining a diagnostic accuracy analogous to conventional, delayed imaging. METHODS: One hundred and eight patients with suspected Parkinson's disease (PD) were acquired after 15 and 240 minutes from the injection of 150-185 MBq of 123I-mIBG. The early and late H/M (He/Me and Hl/Ml) were evaluated by drawing Region-of-Interests on the heart and the upper half of the mediastinum. Optimal threshold (Youden index) and overall predictive performance were determined by receiver operating characteristic curve, classifying tentatively patients having an Hl/Ml lower than 1.6 as suffering from PD. RESULTS: He/Me was not significantly different from Hl/Ml (p-value=0.835). The Area-under-curve was 0.935 (95%CI: 0.845-1.000). The He/Me optimal threshold was 1.66, with sensitivity, specificity, and diagnostic accuracy of 95.5%, 85.7 and 90.7% respectively. CONCLUSION: The He/Me Ratio is almost as accurate as the widely used delayed 123I-mIBG imaging, reducing the burden of delayed imaging but preserving the diagnostic accuracy of the method. Moreover the differential diagnosis in Parkinson's disease can be made in just 25 minutes against the 4 hours currently needed, lowering costs of the healthcare system and improving patients compliance.

Deep Seated Tumour Treatments With Electrons of High Energy Delivered at FLASH Rates: The Example of Prostate Cancer.

Different therapies are adopted for the treatment of deep seated tumours in combination or as an alternative to surgical removal or chemotherapy: radiotherapy with photons (RT), particle therapy (PT) with protons or even heavier ions like 12C, are now available in clinical centres. In addition to these irradiation modalities, the use of Very High Energy Electron (VHEE) beams (100-200 MeV) has been suggested in the past, but the diffusion of that technique was delayed due to the needed space and budget, with respect to standard photon devices. These disadvantages were not paired by an increased therapeutic efficacy, at least when comparing to proton or carbon ion beams. In this contribution we investigate how recent developments in electron beam therapy could reshape the treatments of deep seated tumours. In this respect we carefully explored the application of VHEE beams to the prostate cancer, a well-known and studied example of deep seated tumour currently treated with high efficacy both using RT and PT. The VHEE Treatment Planning System was obtained by means of an accurate Monte Carlo (MC) simulation of the electrons interactions with the patient body. A simple model of the FLASH effect (healthy tissues sparing at ultra-high dose rates), has been introduced and the results have been compared with conventional RT. The study demonstrates that VHEE beams, even in absence of a significant FLASH effect and with a reduced energy range (70-130 MeV) with respect to implementations already explored in literature, could be a good alternative to standard RT, even in the framework of technological developments that are nowadays affordable.

PET-CT post therapy dosimetry in radioembolization with resin 90Y microspheres: Comparison with pre-treatment SPECT-CT 99mTc-MAA results.

Resin microspheres radioembolization is an effective treatment for liver tumors when the surgical option is not feasible. Doses delivered to tumor and normal liver can be assess in the pre-therapy phase by means of a 99mTc-MAA SPECT-CT simulation and after the treatment with 90Y PET-CT acquisition. The optimal therapeutic 90Y activity is determined on 99mTc-MAA SPECT-CT dose results in order to avoid healthy parenchyma toxicity and to effectively irradiate the tumor. The assumption of identical radiopharmaceutical distribution between simulation and verification is still under debate and literature data showed controversial results. In this study 10 HCC patient's dosimetry performed on 99mTc SPECT-CT and 90Y PET-CT were compared. Patients were selected when a good agreement between the pre and post-therapy distribution was observed in order to investigate the intrinsic dosimetric variations between the two imaging modalities. Mean doses (MIRD and Voxel approaches) showed a good correlation (Pearson's coefficient r > 0.90) both for tumor and normal liver. Dose Volume Histogram curves were compared with a good agreement particularly for normal liver (D50). Goal doses were achieved for 90% of patients. Bland-Altman analysis indicates lower variations for healthy parenchyma than for tumor (1.96 SD equal to 9.1 Gy and 68 Gy respectively) confirming the robustness of the dose-toxicity approach. PET-CT dosimetry well correlates with SPECT-CT doses (under assumption of same catheter position and 90Y activity). Better agreement was showed for 7/10 and 8/10 patients for T and NL respectively, confirming dosimetry as effective tool to optimize and individualize the treatment.

T1 and extracellular volume fraction mapping in cardiac magnetic resonance: estimation of accuracy and precision of a novel algorithm.

Magnetic resonance imaging is effective for non-invasive detection of myocardial diseases by extracellular volume fraction (ECV) estimation. A new methodology for T1 and ECV mapping is tested in this work, comparing results with other well-consolidated methods. The associated level of uncertainty for data was also estimated, to assess the reliability of the technique. A phantom with known T1 values was used for reference, and 60 subjects (40 controls, 20 diseased patients) were examined, using the modified look-locker inversion-recovery (MOLLI) method. Obtained T1 data were studied in terms of accuracy (difference with reference T1), paired comparisons with other methods, and Gamma-tool analysis with tolerances criteria of 4.21 mm for distance-to-agreement, and between 2%-5% for T1 difference. Accuracy and precision of the T1 mapping was analysed by phantom measurements, and the uncertainty of the ECV was estimated by T1 error propagation. Differences (in paired comparisons) of T1 datasets were not significant neither for statistical tests, nor for Bland-Altman analysis. T1 accuracy was between -12% and -1% across methods, slightly better for the tested method (mean accuracy in the T1 range of interest better than 2%). The Gamma analysis confirm slightly better results for the tested method than other methodologies. The new method improves the computational efficiency by a factor of 25 (at least), revealing to be also more suitable for Big Data-related applications.

Correlation of dose with toxicity and tumour response to 90Y- and 177Lu-PRRT provides the basis for optimization through individualized treatment planning.

PURPOSE: Peptide receptor radionuclide therapy (PRRT) with 90Y-labelled and 177Lu-labelled peptides is an effective strategy for the treatment of metastatic/nonresectable neuroendocrine tumours (NETs). Dosimetry provides important information useful for optimizing PRRT with individualized regimens to reduce toxicity and increase tumour responses. However, this strategy is not applied in routine clinical practice, despite the fact that several dosimetric studies have demonstrated significant dose-effect correlations for normal organ toxicity and tumour response that can better guide therapy planning. The present study reviews the key relationships and the radiobiological models available in the literature with the aim of providing evidence that optimization of PRRT is feasible through the implementation of dosimetry. METHODS: The MEDLINE database was searched combining specific keywords. Original studies published in the English language reporting dose-effect outcomes in patients treated with PRRT were chosen. RESULTS: Nine of 126 studies were selected from PubMed, and a further five were added manually, reporting on 590 patients. The studies were analysed and are discussed in terms of weak and strong elements of correlations. CONCLUSION: Several studies provided evidence of clinical benefit from the implementation of dosimetry in PRRT, indicating the potential contribution of this approach to reducing severe toxicity and/or reducing undertreatment that commonly occurs. Prospective trials, possibly multicentre, with larger numbers of patients undergoing quantitative dosimetry and with standardized methodologies should be carried out to definitively provide robust predictive paradigms to establish effective tailored PRRT.

Dosimetry-Based Consideration on Remission and Relapse after Therapy with 223Ra-Dichloride in Castration-Resistant Prostate Cancer (CRPC) with Bone Metastases. A Case Report.

Here, we present the case of a 64-year-old male patient diagnosed with castration-resistant prostate cancer (CRPC) with bone metastasis, treated with abiraterone prednisone/prednisolone in combination with 223Ra-dichloride therapy, who had remission and a subsequent relapse of bone metastasis on repeated bone scans after therapy. We also discuss the possibility of continuing the 223Ra-dichloride therapy over the six planned administrations by administering other cycles at the same dose or at higher doses, if shown to be devoid of a significant increase in side effects, based on dosimetry considerations.

Impact of missing attenuation and scatter corrections on 99m Tc-MAA SPECT 3D dosimetry for liver radioembolization using the patient relative calibration methodology: A retrospective investigation on clinical images.

PURPOSE: To investigate the clinical implication of performing pre-treatment dosimetry for 90 Y-microspheres liver radioembolization on 99m Tc-MAA SPECT images reconstructed without attenuation or scatter correction and quantified with the patient relative calibration methodology. METHODS: Twenty-five patients treated with SIR-Spheres® at Istituto Europeo di Oncologia and 31 patients treated with TheraSphere® at Istituto Nazionale Tumori were considered. For each acquired 99m Tc-MAA SPECT, four reconstructions were performed: with attenuation and scatter correction (AC_SC), only attenuation (AC_NoSC), only scatter (NoAC_SC) and without corrections (NoAC_NoSC). Absorbed dose maps were calculated from the activity maps, quantified applying the patient relative calibration to the SPECT images. Whole Liver (WL) and Tumor (T) regions were drawn on CT images. Injected Liver (IL) region was defined including the voxels receiving absorbed dose >3.8 Gy/GBq. Whole Healthy Liver (WHL) and Healthy Injected Liver (HIL) regions were obtained as WHL = WL - T and HIL = IL - T. Average absorbed dose to WHL and HIL were calculated, and the injection activity was derived following each Institute's procedure. The values obtained from AC_NoSC, NoAC_SC and NoAC_NoSC images were compared to the reference value suggested by AC_SC images using Bland-Altman analysis and Wilcoxon paired test (5% significance threshold). Absorbed-dose maps were compared to the reference map (AC_SC) in global terms using the Voxel Normalized Mean Square Error (%VNMSE), and at voxel level by calculating for each voxel the normalized difference with the reference value. The uncertainty affecting absorbed dose at voxel level was accounted for in the comparison; to this purpose, the voxel counts fluctuation due to Poisson and reconstruction noise was estimated from SPECT images of a water phantom acquired and reconstructed as patient images. RESULTS: NoAC_SC images lead to activity prescriptions not significantly different from the reference AC_SC images; the individual differences (<0.1 GBq for all IEO patients, <0.6 GBq for all but one INT patients) were comparable to the uncertainty affecting activity measurement. AC_NoSC and NoAC_NoSC images, instead, yielded significantly different activity prescriptions and wider 95% confidence intervals in the Bland-Altman analysis. Concerning the absorbed dose map, AC_NoSC images had the smallest %VNMSE value and the highest fraction of voxels differing less than 2 standard deviations from AC_SC. CONCLUSIONS: The patient relative calibration methodology can compensate for the missing attenuation correction when performing healthy liver pre-treatment dosimetry: safe treatments can be planned even on NoAC_SC images, suggesting activities comparable to AC_SC images. Scatter correction is recommended due to its heavy impact on healthy liver dosimetry.

Myocardial blood flow estimates from dynamic contrast-enhanced magnetic resonance imaging: three quantitative methods.

Dynamic contrast-enhanced cardiovascular magnetic resonance imaging can be used to quantitatively assess the myocardial blood flow (MBF), recovering the tissue impulse response function for the transit of a gadolinium bolus through the myocardium. Several deconvolution techniques are available, using various models for the impulse response. The method of choice may influence the results, producing differences that have not been deeply investigated yet. Three methods for quantifying myocardial perfusion have been compared: Fermi function modelling (FFM), the Tofts model (TM) and the gamma function model (GF), with the latter traditionally used in brain perfusion MRI. Thirty human subjects were studied at rest as well as under cold pressor test stress (submerging hands in ice-cold water), and a single bolus of gadolinium weighing 0.1 ± 0.05 mmol kg-1 was injected. Perfusion estimate differences between the methods were analysed by paired comparisons with Student's t-test, linear regression analysis, and Bland-Altman plots, as well as also using the two-way ANOVA, considering the MBF values of all patients grouped according to two categories: calculation method and rest/stress conditions. Perfusion estimates obtained by various methods in both rest and stress conditions were not significantly different, and were in good agreement with the literature. The results obtained during the first-pass transit time (20 s) yielded p-values in the range 0.20-0.28 for Student's t-test, linear regression analysis slopes between 0.98-1.03, and R values between 0.92-1.01. From the Bland-Altman plots, the paired comparisons yielded a bias (and a 95% CI)-expressed as ml/min/g-for FFM versus TM, -0.01 (-0.20, 0.17) or 0.02 (-0.49, 0.52) at rest or under stress respectively, for FFM versus GF, -0.05 (-0.29, 0.20) or -0.07 (-0.55, 0.41) at rest or under stress, and for TM versus GF, -0.03 (-0.30, 0.24) or -0.09 (-0.43, 0.26) at rest or under stress. With the two-way ANOVA, the results were p = 0.20 for the method effect (not significant), p < 0.0001 for the rest/stress condition effect (highly significant, as expected), whereas no interaction resulted between the rest/stress condition and method (p = 0.70, not significant). Considering a wider time-frame (60 s), the estimates for both rest and stress conditions were 25%-30% higher (p in the range 0.016-0.025) than those obtained in the 20 s time-frame. MBF estimates obtained by various methods under rest/stress conditions were not significantly different in the first-pass transit time, encouraging quantitative perfusion estimates in DCE-CMRI with the used methods.

A 3-variable prognostic score (3-PS) for overall survival prediction in metastatic castration-resistant prostate cancer treated with 223Radium-dichloride.

OBJECTIVE: In mCRPC patients treated with 223Ra, a major issue is the validation of reliable prognostic and predictive biomarkers to maximize clinical benefit and minimize toxicities and costs. Bearing in mind how changes in tALP did not meet statistical requirements as surrogate marker for survival, aim of this single-center retrospective study was to characterize the prognostic and predictive role of baseline clinical variables associated with overall survival in patients receiving 223Ra treatment. METHODS: 92 consecutive CRPC patients with symptomatic bone metastases receiving 223Ra treatment were included. Available baseline clinical data relevant to the survival analysis were retrospectively collected. The primary end-point of the study was overall survival, which was established from the first 223Ra administration until date of death from any cause. RESULTS: Median follow-up time from the first 223Ra administration was 6 months (range 1-31 months). The univariate analysis evaluating the prognostic value of all baseline clinical variables showed that patients' weight, BMI, ECOG PS, Hb and tALP values were independently associated with OS. On multivariable analysis only baseline Hb value and ECOG PS remained significantly correlated with OS. To determine reliable baseline predictive factors for survival in patients receiving 223Ra treatment, we produced a predictive score. We tried all possible variable combinations, and found that the best score was obtained by combining baseline ECOG PS with Hb < 12 g/dl and PSA ≥ 20 ng/ml. This resulted in a score ranging from 0 to 4, with AUC 78.4% (p < 0.001). CONCLUSIONS: We propose a multidimensional clinical evaluation to select those mCRPC subjects suitable to receive the maximum benefit from 223Ra treatment.

Gamma camera calibrations for the Italian multicentre study for lesion dosimetry in 223Ra therapy of bone metastases.

PURPOSE: The aim was to calibrate gamma cameras in the framework of the Italian multicentre study for lesion dosimetry in 223Ra therapy of bone metastases. Equipments of several manufacturers and different models were used. METHODS: Eleven gamma cameras (3/8- and 5/8-inch crystal) were used, acquiring planar static images with double-peak (82 and 154keV, 20% wide) and MEGP collimator. The sensitivity was measured in air, varying source-detector distance and source size. Transmission curves were measured, calculating the parameters used for attenuation/scatter correction with the pseudo-extrapolation number method, and assessing their variations with the source size. RESULTS: Values of the calibration factor (geometric mean of both detector sensitivities) ranged from 41.1 to 113.9cps/MBq. For the smallest source (diameter of 3.5cm), the calibration factor decrease ranged from -30% to -4%, highlighting the importance of partial volume effects according to the equipment involved. The sensitivity variation with the source-detector distance, with respect to the 15cm-value, reached 10% (in absolute value) in the range 5-30cm, but fixing the distance between the two heads, the calibration factor variation with the distance from the midline was within 3.6%. Appreciable variation of the transmission curves with the source size were observed, examining the results obtained with six gamma cameras. CONCLUSION: Assessments of sensitivity and transmission curve variations with source size should be regularly implemented in calibration procedures. The results of this study represent a useful compendium to check the obtained calibrations for dosimetric purposes.

The Italian multicentre dosimetric study for lesion dosimetry in 223Ra therapy of bone metastases: Calibration protocol of gamma cameras and patient eligibility criteria.

PURPOSE: The aims of this work were to explore patient eligibility criteria for dosimetric studies in 223Ra therapy and evaluate the effects of differences in gamma camera calibration procedures into activity quantification. METHODS: Calibrations with 223Ra were performed with four gamma cameras (3/8-inch crystal) acquiring planar static images with double-peak (82 and 154keV, 20% wide) and MEGP collimator. The sensitivity was measured in air by varying activity, source-detector distance, and source diameter. Transmission curves were measured for attenuation/scatter correction with the pseudo-extrapolation number method, varying the experimental setup. 223Ra images of twenty-five patients (69 lesions) were acquired to study the lesions visibility. Univariate ROC analysis was performed considering visible/non visible lesions on 223Ra images as true positive/true negative group, and using as score value the lesion/soft tissue contrast ratio (CR) derived from 99mTc-MDP WB scan. RESULTS: Sensitivity was nearly constant varying activity and distance (maximum s.d.=2%). Partial volume effects were negligible for object area ⩾960mm2. Transmission curve measurements are affected by experimental setup and source size, leading to activity quantification errors up to 20%. The ROC analysis yielded an AUC of 0.972 and an optimal threshold of CR of 10, corresponding to an accuracy of 92%. CONCLUSION: The minimum calibration protocol requires sensitivity and transmission curve measurements varying the object size, performing a careful procedure standardisation. Lesions with 99mTc-MDP CR higher than 10, not overlapping the GI tract, are generally visible on 223Ra images acquired at 24h after the administration, and possibly eligible for dosimetric studies.

Impact of SPECT corrections on 3D-dosimetry for liver transarterial radioembolization using the patient relative calibration methodology.

PURPOSE: Many centers aim to plan liver transarterial radioembolization (TARE) with dosimetry, even without CT-based attenuation correction (AC), or with unoptimized scatter correction (SC) methods. This work investigates the impact of presence vs absence of such corrections, and limited spatial resolution, on 3D dosimetry for TARE. METHODS: Three voxelized phantoms were derived from CT images of real patients with different body sizes. Simulations of (99m)Tc-SPECT projections were performed with the SIMIND code, assuming three activity distributions in the liver: uniform, inside a "liver's segment," or distributing multiple uptaking nodules ("nonuniform liver"), with a tumoral liver/healthy parenchyma ratio of 5:1. Projection data were reconstructed by a commercial workstation, with OSEM protocol not specifically optimized for dosimetry (spatial resolution of 12.6 mm), with/without SC (optimized, or with parameters predefined by the manufacturer; dual energy window), and with/without AC. Activity in voxels was calculated by a relative calibration, assuming identical microspheres and (99m)Tc-SPECT counts spatial distribution. 3D dose distributions were calculated by convolution with (90)Y voxel S-values, assuming permanent trapping of microspheres. Cumulative dose-volume histograms in lesions and healthy parenchyma from different reconstructions were compared with those obtained from the reference biodistribution (the "gold standard," GS), assessing differences for D95%, D70%, and D50% (i.e., minimum value of the absorbed dose to a percentage of the irradiated volume). γ tool analysis with tolerance of 3%/13 mm was used to evaluate the agreement between GS and simulated cases. The influence of deep-breathing was studied, blurring the reference biodistributions with a 3D anisotropic gaussian kernel, and performing the simulations once again. RESULTS: Differences of the dosimetric indicators were noticeable in some cases, always negative for lesions and distributed around zero for parenchyma. Application of AC and SC reduced systematically the differences for lesions by 5%-14% for a liver segment, and by 7%-12% for a nonuniform liver. For parenchyma, the data trend was less clear, but the overall range of variability passed from -10%/40% for a liver segment, and -10%/20% for a nonuniform liver, to -13%/6% in both cases. Applying AC, SC with preset parameters gave similar results to optimized SC, as confirmed by γ tool analysis. Moreover, γ analysis confirmed that solely AC and SC are not sufficient to obtain accurate 3D dose distribution. With breathing, the accuracy worsened severely for all dosimetric indicators, above all for lesions: with AC and optimized SC, -38%/-13% in liver's segment, -61%/-40% in the nonuniform liver. For parenchyma, D50% resulted always less sensitive to breathing and sub-optimal correction methods (difference overall range: -7%/13%). CONCLUSIONS: Reconstruction protocol optimization, AC, SC, PVE and respiratory motion corrections should be implemented to obtain the best possible dosimetric accuracy. On the other side, thanks to the relative calibration, D50% inaccuracy for the healthy parenchyma from absence of AC was less than expected, while the optimization of SC was scarcely influent. The relative calibration therefore allows to perform TARE planning, basing on D50% for the healthy parenchyma, even without AC or with suboptimal corrections, rather than rely on nondosimetric methods.

A case report of image-based dosimetry of bone metastases with Alpharadin ((223)Ra-dichloride) therapy: inter-fraction variability of absorbed dose and follow-up.

A 70-year-old man affected by bone metastases from castration resistant prostate cancer underwent Alpharadin ((223)Ra-dichloride) therapy (6 administrations of 50 kBq per kg i.v., once every 4 weeks). The inter-fraction variability of the absorbed dose to lesions was evaluated for four injections. Dosimetric assessments were performed following the MIRD approach and a recently published methodology. The mean absorbed dose and standard deviation for 4 lesions [mean (σ %)] were: 434 mGy (15%) and 516 mGy (21%) for the right and left humeral head, 1205 mGy (14%) and 781 mGy (8%) for the right and left glenoid. The estimated total absorbed dose after the whole treatment, considering also the relative-biological effectiveness of alpha particles (RBE = 5), yielded a D RBE range of 13-36 Gy. A good correlation between (99m)Tc and (223)Ra uptake was obtained (R (2) = 0.7613). The tumour-non-tumour (TNT) ratio of 8 lesions (those above, plus 4 additional), monitored by six (99m)Tc-MDP bone scans over a period of about 10 months, evidenced a TNT reduction in two lesions (-42 and -48 %), but in most lesions the TNT remained fairly constant, evidencing that (223)Ra-dichloride therapy tends to prevent further progression of osseous disease, leading to chronicity of the metastatic status.

Dosimetry of bone metastases in targeted radionuclide therapy with alpha-emitting (223)Ra-dichloride.

PURPOSE: Ra-dichloride is an alpha-emitting radiopharmaceutical used in the treatment of bone metastases from castration-resistant prostate cancer. Image-based dosimetric studies remain challenging because the emitted photons are few. The aim of this study was to implement a methodology for in-vivo quantitative planar imaging, and to assess the absorbed dose to lesions using the MIRD approach. METHODS: The study included nine Caucasian patients with 24 lesions (6 humeral head lesions, 4 iliac wing lesions, 2 scapular lesions, 5 trochanter lesions, 3 vertebral lesions, 3 glenoid lesions, 1 coxofemoral lesion). The treatment consisted of six injections (one every 4 weeks) of 50 kBq per kg body weight. Gamma-camera calibrations for (223)Ra included measurements of sensitivity and transmission curves. Patients were statically imaged for 30 min, using an MEGP collimator, double-peak acquisition, and filtering to improve the image quality. Lesions were delineated on (99m)Tc-MDP whole-body images, and the ROIs superimposed on the (223)Ra images after image coregistration. The activity was quantified with background, attenuation, and scatter correction. Absorbed doses were assessed deriving the S values from the S factors for soft-tissue spheres of OLINDA/EXM, evaluating the lesion volumes by delineation on the CT images. RESULTS: In 12 lesions with a wash-in phase the biokinetics were assumed to be biexponential, and to be monoexponential in the remainder. The optimal timing for serial acquisitions was between 1 and 5 h, between 18 and 24 h, between 48 and 60 h, and between 7 and 15 days. The error in cumulated activity neglecting the wash-in phase was between 2 % and 12 %. The mean effective half-life (T 1/2eff) of (223)Ra was 8.2 days (range 5.5-11.4 days). The absorbed dose (D) after the first injection was 0.7 Gy (range 0.2-1.9 Gy. Considering the relative biological effectiveness (RBE) of alpha particles (RBE = 5), D RBE = 899 mGy/MBq (range 340-2,450 mGy/MBq). The percent uptake of (99m)Tc and (223)Ra (activity extrapolated to t = 0) were significantly correlated. CONCLUSION: The feasibility of in vivo quantitative imaging in (223)Ra therapy was confirmed. The lesion uptake of (223)Ra-dichloride was significantly correlated with that of (99m)Tc-MDP. The D RBE to lesions per unit administered activity was much higher than that of other bone-seeking radiopharmaceuticals, but considering a standard administration of 21 MBq (six injections of 50 kBq/kg to a 70-kg patient), the mean cumulative value of D RBE was about 19 Gy, and was therefore in the range of those of other radiopharmaceuticals. The macrodosimetry of bone metastases in treatments with (223)Ra-dichloride is feasible, but more work is needed to demonstrate its helpfulness in predicting clinical outcomes.

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.

Improving the dose-myelotoxicity correlation in radiometabolic therapy of bone metastases with 153Sm-EDTMP.

PURPOSE: (153)Sm-ethylene diamine tetramethylene phosphonic acid ((153)Sm-EDTMP) is widely used to palliate pain from bone metastases, and is being studied for combination therapy beyond palliation. Conceptually, red marrow (RM) dosimetry allows myelotoxicity to be predicted, but the correlation is poor due to dosimetric uncertainty, individual sensitivity and biological effects from previous treatments. According to EANM guidelines, basic dosimetric procedures have been studied to improve the correlation between dosimetry and myelotoxicity in (153)Sm-EDTMP therapy. METHODS: RM dosimetry for 33 treatments of bone metastases from breast, prostate and lung tumours was performed prospectively (with (99m)Tc-MDP) and retrospectively, acquiring whole-body scans early and late after injection. The (153)Sm-EDTMP activity was calculated by prospective dosimetry based on measured skeletal uptake and full physical retention, with the RM absorbed dose not exceeding 3.8 Gy. Patient-specific RM mass was evaluated by scaling in terms of body weight (BW), lean body mass (LBM) and trabecular volume (TV) estimated from CT scans of the L2­L4 vertebrae. Correlations with toxicity were determined in a selected subgroup of 27 patients, in which a better correlation between dosimetry and myelotoxicity was expected. RESULTS: Skeletal uptakes of (99m)Tc and (153)Sm (Tc% and Sm%) were well correlated. The median Sm% was higher in prostate cancer (75.3 %) than in lung (60.5%, p = 0.005) or breast (60.8%, p = 0.008). PLT and WBC nadirs were not correlated with administered activity, but were weakly correlated with uncorrected RM absorbed doses, and the correlation improved after rescaling in terms of BW, LBM and TV. Most patients showed transient toxicity (grade 1­3), which completely and spontaneously recovered over a few days. Using TV, RM absorbed dose was in the range 2­5 Gy, with a median of 312 cGy for PLT in patients with toxicity and 247 cGy in those with no toxicity (p = 0.019), and 312 cGy for WBC in those with toxicity and 232 cGy in those with no toxicity (p = 0.019). ROC curves confirmed the correlations, yielding toxicity absorbed dose thresholds of 265 cGy for PLT and 232 cGy for WBC. CONCLUSION: The best predictor of myelotoxicity and blood cells nadir was obtained scaling the RM absorbed dose in terms of the estimated TV. It seems clear that the increase in skeletal uptake due to the presence of bone metastases and the assumption of full physical retention cause an overestimation of the RM absorbed dose. Nevertheless, an improvement of the dose­toxicity correlation is easily achievable by simple methods, also leading to possible improvement in multifactorial analyses of myelotoxicity.