A Comparison of Techniques for (90)Y PET/CT Image-Based Dosimetry Following Radioembolization with Resin Microspheres

Translating contrast enhanced computed tomography images to liver radioembolization dose distribution for more comprehensively indicating patients

Objective.Contrast-enhanced computed tomography (CECT) is commonly used in the pre-treatment evaluation of liver Y-90 radioembolization feasibility. CECT provides detailed imaging of the liver and surrounding structures, allowing healthcare providers to assess the size, location, and characteristics of liver tumors prior to the treatment. Here we propose a method for translating CECT images to an expected dose distribution for tumor(s) and normal liver tissue.Approach.A pre-procedure CECT is used to obtain an iodine arterial-phase distribution by subtracting the non-contrast CT from the late arterial phase. The liver segments surrounding the targeted tumor are selected using Couinaud's method. The resolution of the resulting images is then degraded to match the resolution of the positron emission tomography (PET) images, which can image the Y-90 activity distribution post-treatment. The resulting images are then used in the same way as PET images to compute doses using the local deposition method. CECT images from three patients were used to test this method retrospectively and were compared with Y-90 PET-based dose distributions through dose volume histograms.Main results.Results show a concordance between predicted and delivered Y-90 dose distributions with less than 10% difference in terms of mean dose, for doses greater than 10% of the 98th percentile (D2%).Significance.CECT-derived predictions of Y-90 radioembolization dose distributions seem promising as a supplementary tool for physicians when assessing treatment feasibility. This dosimetry prediction method could provide a more comprehensive pre-treatment evaluation-offering greater insights than a basic assessment of tumor opacification on CT images.

Establishing Updated Safety Standards for Independent 99mTc-MAA SPECT/CT Treatment Planning in Radioembolization

PURPOSE

Significant improvements within radioembolization imaging and dosimetry permit the development of an accurate and personalized pretreatment plan using technetium 99m-labeled macroaggregated albumin (99mTc-MAA) and single-photon emission computed tomography (SPECT) combined with anatomical CT (SPECT/CT). Despite these potential advantages, the clinical transition to pretreatment protocols with SPECT/CT is hindered by their unknown safety constraints. This study aimed to address this issue by establishing novel dose limits for 99mTc-MAA SPECT/CT to enable quantitative pretreatment planning.

METHODS AND MATERIALS

Stratification criteria to determine images most viable for dosimetry analysis were created from a cohort of 85 patients. SPECT/CT, cone beam CT, and activity calculations derived from the local deposition method were used to create an accurate pretreatment protocol. Planar and SPECT/CT images were compared using linear regression and modified Bland-Altman analyses to convert accepted planar dose limits to SPECT/CT. To validate these new dose limits, activity calculations based on SPECT/CT were compared with those calculated with the body surface area and planar methods for three treatment plans.

RESULTS

A total of 38 of 85 patients were deemed viable for dosimetry analysis. SPECT yielded greater lung shunt fractions (LSFs) than planar imaging when LSFs were <4.89%, whereas SPECT yielded lower LSFs than planar imaging when LSFs were >4.89%. Planar to SPECT/CT dose conversions were 0.76×, 0.70×, and 0.55× for the whole liver, normal liver, and lungs, respectively. Patients with SPECT LSFs ≤4.89% were safely treated with the direct application of planar lung dose limits. Activity calculations with the newly established SPECT/CT dose limits were greater than those of the body surface area method by a median range of 33.1% to 61.9% and were lower than planar-based activity calculations by a median range of 12.5% to 13.7% for the whole liver and by 29.4% to 32.2% for the normal liver.

CONCLUSIONS

This study demonstrated a safe method for translating dose limits from 99mTc-MAA planar imaging to SPECT/CT. A robust pretreatment protocol was further developed guided by the current knowledge in the field. Established SPECT/CT dose limits safely treated 97.5% of patients and permitted the application of independent pretreatment planning with 99mTc-MAA SPECT/CT.

Characterizing the voxel-based approaches in radioembolization dosimetry with reDoseMC

BACKGROUND

Yttrium-90 (90 Y $^{90}{\rm {Y}}$ ) represents the primary radioisotope used in radioembolization procedures, while holmium-166 (166 Ho $^{166}{\rm {Ho}}$ ) is hypothesized to serve as a viable substitute for90 Y $^{90}{\rm {Y}}$ due to its comparable therapeutic potential and improved quantitative imaging. Voxel-based dosimetry for these radioisotopes relies on activity images obtained through PET or SPECT and dosimetry methods, including the voxel S-value (VSV) and the local deposition method (LDM). However, the evaluation of the accuracy of absorbed dose calculations has been limited by the use of non-ideal reference standards and investigations restricted to the liver. The objective of this study was to expand upon these dosimetry characterizations by investigating the impact of image resolutions, voxel sizes, target volumes, and tissue materials on the accuracy of90 Y $^{90}{\rm {Y}}$ and166 Ho $^{166}{\rm {Ho}}$ dosimetry techniques.

METHODS

A specialized radiopharmaceutical dosimetry software called reDoseMC was developed using the Geant4 Monte Carlo toolkit and validated by benchmarking the generated90 Y $^{90}{\rm {Y}}$ kernels with published data. The decay spectra of both90 Y $^{90}{\rm {Y}}$ and166 Ho $^{166}{\rm {Ho}}$ were also compared. Multiple VSV kernels were generated for the liver, lungs, soft tissue, and bone for isotropic voxel sizes of 1 mm, 2 mm, and 4 mm. Three theoretical phantom setups were created with 20 or 40 mm activity and mass density inserts for the same three voxel sizes. To replicate the limited spatial resolutions present in PET and SPECT images, image resolutions were modeled using a 3D Gaussian kernel with a Full Width at Half Maximum (FWHM) ranging from 0 to 16 mm and with no added noise. The VSV and LDM dosimetry methods were evaluated by characterizing their respective kernels and analyzing their absorbed dose estimates calculated on theoretical phantoms. The ground truth for these estimations was calculated using reDoseMC.

RESULTS

The decay spectra obtained through reDoseMC showed less than a 1% difference when compared to previously published experimental data for energies below 1.9 MeV in the case of90 Y $^{90}{\rm {Y}}$ and less than 1% for energies below 1.5 MeV for166 Ho $^{166}{\rm {Ho}}$ . Additionally, the validation kernels for90 Y $^{90}{\rm {Y}}$ VSV exhibited results similar to those found in published Monte Carlo codes, with source dose depositions having less than a 3% error margin. Resolution thresholds (FWHM thresh s ${\rm {FWHM}}_\mathrm{thresh}{\rm {s}}$ ), defined as resolutions that resulted in similar dose estimates between the LDM and VSV methods, were observed for90 Y $^{90}{\rm {Y}}$ . They were 1.5 mm for bone, 2.5 mm for soft tissue and liver, and 8.5 mm for lungs. For166 Ho $^{166}{\rm {Ho}}$ , the accuracy of absorbed dose deposition was found to be dependent on the contributions of absorbed dose from photons. Volume errors due to variations in voxel size impacted the final dose estimates. Larger target volumes yielded more accurate mean doses than smaller volumes. For both radioisotopes, the radial dose profiles for the VSV and LDM approximated but never matched the reference standard.

CONCLUSIONS

reDoseMC was developed and validated for radiopharmaceutical dosimetry. The accuracy of voxel-based dosimetry was found to vary widely with changes in image resolutions, voxel sizes, chosen target volumes, and tissue material; hence, the standardization of dosimetry protocols was found to be of great importance for comparable dosimetry analysis.

Quantitative evaluation of 90Y-PET/CT and 90Y-SPECT/CT-based dosimetry following Yttrium-90 radioembolization

BACKGROUND

Following yttrium-90 radioembolization (90Y-RE), 90Y-PET/CT and 90Y-SPECT/CT imaging provide the means to calculate the voxelized absorbed dose distribution. Given the widespread use of the two imaging modalities and lack of well-established standardized dosimetry protocols for 90Y-RE, there is a clinical need to systematically investigate and evaluate differences in the performance of voxel-based dosimetry between 90Y-PET/CT and 90Y-SPECT/CT.

PURPOSE

To quantitatively analyze and compare 90Y-PET/CT and 90Y-SPECT/CT-based dosimetry following 90Y-RE.

METHODS

90Y-PET/CT and 90Y-SPECT/CT imaging was acquired for 35 patients following 90Y-RE with TheraSphere for the treatment of unresectable hepatocellular carcinoma. Dosimetry was performed using the local deposition method with known activity and the mean dose (Dmean) was calculated for perfused liver volumes (PV), tumors (T), and perfused normal livers (NL). Additionally, the absorbed dose to x% of the volume (Dx, x ∈ $ \in $ [5%, 10%, …, 90%, 95%]) and the volume receiving y Gy (Vy, y ∈ $ \in $ [10 Gy, 20 Gy, …, 190 Gy, 200 Gy]) were calculated for T and NL, respectively. Dose metrics were compared using linear regression, Bland-Altman analysis, and statistical testing.

RESULTS

Both 90Y-SPECT/CT and 90Y-PET/CT-based tumor Dmean were strongly correlated (R2 ≥ 0.90) with Dx, excluding metrics on the extrema. Intra-modality comparisons of various Dx and Vy metrics yielded statistically significant differences (ANOVA, p < 0.001) for both90Y-PET/CT and 90Y-SPECT/CT. Based on statistical testing, only Dx metrics separated by greater than 20%-30% coverage, and only Vy metrics separated by greater than 40-70 Gy, reported significant differences. For PV, there was a strong correlation (R2 ≥ 0.99) between Dmean derived separately from 90Y-PET/CT and 90Y-SPECT/CT imaging. The strength of the correlation was slightly reduced for T and NL with R2 = 0.91 and R2 = 0.95, respectively. For PV, the mean bias ± standard error (SE) and 95% limits of agreement (LOA) between Dmean from the two modalities was effectively zero with -0.8 ± 0.4% (± 2.5%). For T and NL, the mean bias ± SE (± LOA) was -14.5 ± 3.7% (± 24%) and 9.4 ± 4.7% (± 27%), respectively.

CONCLUSION

The strong correlation between Dmean and Dx suggests information from multiple dose metrics (e.g., D70 and Dmean) is largely redundant when establishing dose-response relationships in 90Y-RE. Dmean is highly correlated between 90Y-PET/CT and 90Y-SPECT/CT-based dosimetry, for all liver VOIs. Relative to 90Y-SPECT/CT, 90Y-PET/CT, on average, yielded higher Dmean to tumors (14%) and lower Dmean to perfused normal livers (9%). Absorbed dose differences for perfused liver volumes between 90Y-SPECT/CT and 90Y-PET/CT were negligible.

A comparison of methods for in vivo activity and absorbed dose quantification with PET/CT following yttrium-90 radioembolization

BACKGROUND

Yttrium-90 (90Y) positron emission tomography (PET)/computed tomography (CT) imaging is increasingly being used to perform tumor (T) and normal liver (NL) voxel dosimetry after 90Y-radioembolization (90Y-RE). Yet, the accuracy of in vivo 90Y-PET/CT imaging, subject to motion blur and co-registration inaccuracies, and 90Y-PET/CT dose quantification, subject to availability of different voxel dosimetry algorithms, are not well understood.

PURPOSE

The purpose of this study was to investigate the accuracy of 90Y-PET/CT-based activity estimates following 90Y-RE and characterize differences between 90Y-PET/CT-based voxel dosimetry algorithms.

METHODS

Thirty-five patients underwent 90Y-PET/CT imaging after 90Y-RE with TheraSphere. The net administered 90Y activity (Aadmin) was determined using a dose calibrator and pre- and post-procedure exposure rate measurements. The summation of image-based activity (Aimage) was extracted from perfused volume (PV) and 3D-isotropically 2-cm expanded PV contour (PV+2 cm). Absorbed doses were calculated using voxel S-value (VSV), local deposition method (LDM), and LDM with known activity (LDMKA) dosimetry algorithms. Linear regression and Bland-Altman analysis quantified the relationship between Aimage and Aadmin and between mean dose estimates (DLDM, DVSV, DLDM-KA) for PV, T, and perfused NL volumes.

RESULTS

While Aadmin and Aimage in PV were highly correlated (R2 > 0.95), the mean bias ± standard error (SE) and (95% limits of agreement, LOA) was significantly non-zero with -22.7 ± 4.7% (± 28.4%). In PV+2 cm, the mean bias ± SE (± LOA) decreased to 1.3 ± 3.4% (± 18.0%) consistent with zero mean error. DLDM and DVSV were highly correlated (R2 > 0.99) for all volumes of interest (VOIs) and the mean bias ± SE (± LOA) was 2.2 ± 0.2% (± 1.0%), 0.7 ± 0.4% (± 2.8%), and 3.2 ± 0.5% (± 2.8%) for PV, T, and NL, respectively. DLDM-KA and DVSV were correlated with R2 = 0.86, 0.80, and 0.86 for PV, T, and NL, respectively. The mean bias ± SE (± LOA) between DLDM-KA and DVSV was significantly non-zero with -19.6 ± 5.1% (± 31.0%), -20.8 ± 4.4% (± 29.0%), and -18.1 ± 5.3% (± 31.1%) for PV, T, and NL, respectively.

CONCLUSIONS

The summation of Aimage in PV was underestimated relative to Aadmin. Only by accounting for respiratory motion, limited spatial resolution, and PET/CT co-registration errors through VOI expansion was Aimage, on average, equal to Aadmin. The differences between DLDM and DVSV were not clinically relevant, though DLDM-KA was approximately 20% greater than DVSV. Given the high quantitative accuracy of dose calibrators and challenges associated with accurate 90Y-PET/CT quantification, LDMKA is the preferred algorithm for accurate 90Y-PET/CT-based dosimetry following 90Y-RE.

Intraprocedural C-arm dual-phase cone-beam enhancement patterns correlate with tumor absorbed dose after radioembolization

BACKGROUND

Recent studies have shown a clear relationship between absorbed dose and tumor response to treatment after hepatic radioembolization. These findings help to create more personalized treatment planning and dosimetry. However, crucial to this goal is the ability to predict the dose distribution prior to treatment. The microsphere distribution is ultimately determined by (i) the hepatic vasculature and the resulting blood flow dynamics and (ii) the catheter position.

PURPOSE

To show that pretreatment, intra-procedural imaging of blood flow patterns, as quantified by catheter-directed intra-arterial contrast enhancement, correlate with posttreatment microsphere accumulation and, consequently, absorbed dose.

MATERIALS AND METHODS

Patients who participated in a clinical trial (NCT01177007) and for whom both a pretreatment dual-phase contrast-enhanced cone-beam CT (CBCT) and a posttreatment 90Y PET/CT scan were available were included in this retrospective study. Tumors and perfused volumes were manually delineated on the CBCT by an experienced radiologist. The mean, sum, and standard deviation of the voxels in each volume were recorded. The delineations were transferred to the PET-based absorbed dose maps by coregistration of the corresponding CTs. Linear multiple regression was used to correlate pretreatment CBCT enhancement to posttreatment 90Y PET/CT-based absorbed dose in each region. Leave-one-out cross-validation and Bland-Altman analyses were performed on the predicted versus measured absorbed doses.

RESULTS

Nine patients, with a total of 23 tumors were included. All presented with hepatocellular carcinoma (HCC). Visually, all patients had a clear correspondence between CBCT enhancement and absorbed dose. The correlation between CBCT enhancement and posttherapy absorbed tumor dose based was strong (R2 = 0.91), and moderate for the non-tumor liver tissue (R2 = 0.61). Limits of agreement were approximately ±55 Gray for tumor tissue.

CONCLUSION

There is a linear relationship between pretreatment blood dynamics in HCC tumors and posttreatment absorbed dose, which, if shown to be generalizable, allows for pretreatment tumor absorbed dose prediction.

AAPM Medical Physics Practice Guideline 14.a: Yttrium-90 microsphere radioembolization

Radioembolization using Yttrium-90 (90 Y) microspheres is widely used to treat primary and metastatic liver tumors. The present work provides minimum practice guidelines for establishing and supporting such a program. Medical physicists play a key role in patient and staff safety during these procedures. Products currently available are identified and their properties and suppliers summarized. Appropriateness for use is the domain of the treating physician. Patient work up starts with pre-treatment imaging. First, a mapping study using Technetium-99m (Tc-99m ) is carried out to quantify the lung shunt fraction (LSF) and to characterize the vascular supply of the liver. An MRI, CT, or a PET-CT scan is used to obtain information on the tumor burden. The tumor volume, LSF, tumor histology, and other pertinent patient characteristics are used to decide the type and quantity of 90 Y to be ordered. On the day of treatment, the appropriate dose is assayed using a dose calibrator with a calibration traceable to a national standard. In the treatment suite, the care team led by an interventional radiologist delivers the dose using real-time image guidance. The treatment suite is posted as a radioactive area during the procedure and staff wear radiation dosimeters. The treatment room, patient, and staff are surveyed post-procedure. The dose delivered to the patient is determined from the ratio of pre-treatment and residual waste exposure rate measurements. Establishing such a treatment modality is a major undertaking requiring an institutional radioactive materials license amendment complying with appropriate federal and state radiation regulations and appropriate staff training commensurate with their respective role and function in the planning and delivery of the procedure. Training, documentation, and areas for potential failure modes are identified and guidance is provided to ameliorate them.

A single centre intercomparison between commercial treatment planning systems for 90Y radioembolization using virtual and experimental phantoms

INTRODUCTION

Dedicated Treatment Planning Systems (TPSs) were developed to personalize 90Y-transarterial radioembolization. This study evaluated the agreement among four commercial TPSs assessing volumes of interest (VOIs) volumes and dose metrics.

METHODS

A homogeneous (EH) and an anthropomorphic phantom with hot and cold inserts (EA) filled with 99mTc-pertechnetate were acquired with a SPECT/CT scanner. Their virtual versions (VH and VA, respectively) and a phantom with activity inside a single voxel (VK) were generated by an in-house MATLAB script. Images and delineated VOIs were imported into the TPSs to compute voxel-based absorbed dose distributions with various dose deposition approaches: local deposition method (LDM) and dose kernel convolution (DKC) with/without local density correction (LDC). VOI volumes and mean absorbed doses were assessed against their median value across TPSs. Dose-volume histograms (DVHs) and VK-derived dose profiles were evaluated.

RESULTS

Small (<2.1 %) and large (up to 42.4 %) relative volume differences were observed on large (>500 ml) and small VOIs, respectively. Mean absorbed doses relative differences were < 3 % except for small VOIs with steep dose gradients (up to 89.1 % in the VA Cold Sphere VOI). Within the same TPS, LDC negligibly affected the mean absorbed dose, while DKC and LDM showed differences up to 63 %. DHVs were mostly overlapped in experimental phantoms, with some differences in the virtual versions. Dose profiles agreed within 1 %.

CONCLUSION

TPSs showed an overall good agreement except for small VOI volumes and mean absorbed doses of VOIs with steep dose gradients. These discrepancies should be considered in the dosimetry uncertainty assessment, thus requiring an appropriate harmonization.

Pre-treatment dosimetry in90Y-SIRT: Is it possible to optimise SPECT reconstruction parameters and calculation methods for accurate dosimetry?

PURPOSE

The aim of this study was (a) to optimise the99mTc-SPECT reconstruction parameters for the pre-treatment dosimetry of90Y-selective internal radiation therapy (SIRT) and (b) to compare the accuracy of clinical dosimetry methods with full Monte-Carlo dosimetry (fMCD) performed with Gate.

METHODS

To optimise the reconstruction parameters, two hundred reconstructions with different parameters were performed on a NEMA phantom, varying the number of iterations, subsets, and post-filtering. The accuracy of the dosimetric methods was then investigated using an anthropomorphic phantom. Absorbed dose maps were generated using (1) the Partition Model (PM), (2) the Dose Voxel Kernel (DVK) convolution, and (3) the Local Deposition Method (LDM) with known activity restricted to the whole phantom (WP) or to the liver and lungs (LL). The dose to the lungs was calculated using the "multiple DVK" and "multiple LDM" methods.

RESULTS

Optimal OSEM reconstruction parameters were found to depend on object size and dosimetric criterion chosen (Dmean or DVH-derived metric). The Dmean of all three dosimetric methods was close (≤ 10%) to the Dmean of fMCD simulations when considering large segmented volumes (whole liver, normal liver). In contrast, the Dmean to the small volume (∅=31) was systemically underestimated (12%-25%). For lungs, the "multiple DVK" and "multiple LDM" methods yielded a Dmean within 20% for the WP method and within 10% for the LL method.

CONCLUSIONS

All three methods showed a substantial degradation of the dose-volume histograms (DVHs) compared to fMCD simulations. The DVK and LDM methods performed almost equally well, with the "multiple DVK" method being more accurate in the lungs.

99mTc-macroaggregated albumin SPECT/CT predictive dosimetry and dose-response relationship in uveal melanoma liver metastases treated with first-line selective internal radiation therapy

First-line selective internal radiation therapy (SIRT) showed promising outcomes in patients with uveal melanoma liver metastases (UMLM). Patient survival depends on liver's disease control. SIRT planning is essential and little is known about dosimetry. We investigated whether 99mTc-MAA-SPECT/CT dosimetry could predict absorbed doses (AD) evaluated on 90Y-PET/CT and assess the dose-response relationship in UMLM patients treated with first-line SIRT. This IRB-approved, single-center, retrospective analysis (prospectively collected cohort) included 12 patients (median age 63y, range 43-82). Patients underwent MRI/CT, 18F-FDG-PET/CT before and 3-6 months post-SIRT, and 90Y-PET/CT immediately post-SIRT. Thirty-two target lesions were included. AD estimates in tumor and non-tumor liver were obtained from 99mTc-MAA-SPECT/CT and post-SIRT 90Y-PET/CT, and assessed with Lin's concordance correlation coefficients (ρc and Cb), Pearson's coefficient correlation (ρ), and Bland-Altman analyses (mean difference ± standard deviation; 95% limits-of-agreement (LOA)). Influence of tumor characteristics and microsphere type on AD was analyzed. Tumor response was assessed according to size-based, enhancement-based and metabolic response criteria. Mean target lesion AD was 349 Gy (range 46-1586 Gy). Concordance between 99mTc-MAA-SPECT/CT and 90Y-PET/CT tumor dosimetry improved upon dose correction for the recovery coefficient (RC) (ρ = 0.725, ρc = 0.703, Cb = 0.969) with good agreement (mean difference: - 4.93 ± 218.3 Gy, 95%LOA: - 432.8-422.9). Without RC correction, concordance was better for resin microspheres (ρ = 0.85, ρc = 0.998, Cb = 0.849) and agreement was very good between predictive 99mTc-MAA-SPECT/CT and 90Y-PET/CT dosimetry (mean difference: - 4.05 ± 55.9 Gy; 95%LOA: - 113.7-105.6). After RC correction, 99mTc-MAA-SPECT/CT dosimetry overestimated AD (- 70.9 ± 158.9 Gy; 95%LOA: - 382.3-240.6). For glass microspheres, concordance markedly improved with RC correction (ρ = 0.790, ρc = 0.713, Cb = 0.903 vs without correction: ρ = 0.395, ρc = 0.244, Cb = 0.617) and 99mTc-MAA-SPECT/CT dosimetry underestimated AD (148.9 ± 267.5 Gy; 95%LOA: - 375.4-673.2). For non-tumor liver, concordance was good between 99mTc-MAA-SPECT/CT and 90Y-PET/CT dosimetry (ρ = 0.942, ρc = 0.852, Cb = 0.904). 99mTc-MAA-SPECT/CT slightly overestimated liver AD for resin (3.4 ± 3.4 Gy) and glass (11.5 ± 13.9 Gy) microspheres. Tumor AD was not correlated with baseline or post-SIRT lesion characteristics and no dose-response threshold could be identified. 99mTc-MAA-SPECT/CT dosimetry provides good estimates of AD to tumor and non-tumor liver in UMLM patients treated with first-line SIRT.

Lung Dose Measured on Postradioembolization 90Y PET/CT and Incidence of Radiation Pneumonitis

Radiation pneumonitis is a rare but possibly fatal side effect of ⁹⁰Y radioembolization. It may occur 1-6 mo after therapy, if a significant part of the ⁹⁰Y microspheres shunts to the lungs. In current clinical practice, a predicted lung dose greater than 30 Gy is considered a criterion to exclude patients from treatment. However, contrasting findings regarding the occurrence of radiation pneumonitis and lung dose were previously reported in the literature. In this study, the relationship between the lung dose and the eventual occurrence of radiation pneumonitis after ⁹⁰Y radioembolization was investigated. Methods: We retrospectively analyzed 317 ⁹⁰Y liver radioembolization procedures performed during an 8-y period (February 2012 to September 2020). We calculated the predicted lung mean dose (LMD) using ^99mTc-MAA planar scintigraphy (LMD_MAA) acquired during the planning phase and left LMD (LMD_Y-90) using the ⁹⁰Y PET/CT acquired after the treatment. For the lung dose computation, we used the left lung as the representative lung volume, to compensate for scatter from the liver moving in the craniocaudal direction because of breathing and mainly affecting the right lung. Results: In total, 272 patients underwent ⁹⁰Y procedures, of which 63% were performed with glass microspheres and 37% with resin microspheres. The median injected activity was 1,974 MBq (range, 242-9,538 MBq). The median LMD_MAA was 3.5 Gy (range, 0.2-89.0 Gy). For 14 procedures, LMD_MAA was more than 30 Gy. Median LMD_Y-90 was 1 Gy (range, 0.0-22.1 Gy). No patients had an LMD_Y-90 of more than 30 Gy. Of the 3 patients with an LMD_Y-90 of more than 12 Gy, 2 patients (one with an LMD_Y-90 of 22.1 Gy and an LMD_MAA of 89 Gy; the other with an LMD_Y-90 of 17.7 Gy and an LMD_MAA of 34.1 Gy) developed radiation pneumonitis and consequently died. The third patient, with an LMD_Y-90 of 18.4 Gy (LMD_MAA, 29.1 Gy), died 2 mo after treatment, before the imaging evaluation, because of progressive disease. Conclusion: The occurrence of radiation pneumonitis as a consequence of a lung shunt after ⁹⁰Y radioembolization is rare (<1%). No radiation pneumonitis developed in patients with a measured LMD_Y-90 lower than 12 Gy.

Accuracy of two dosimetry software programs for 177Lu radiopharmaceutical therapy using voxel-based patient-specific phantoms

Purpose

Virtual dosimetry using voxel-based patient-specific phantoms and Monte Carlo (MC) simulations offer the advantage of having a gold standard against which absorbed doses may be benchmarked to establish the dosimetry accuracy. Furthermore, these reference values assist in investigating the accuracy of the absorbed dose methodologies from different software programs. Therefore, this study aimed to compare the accuracy of the absorbed doses computed using LundADose and OLINDA/EXM 1.0.

Methods

The accuracy was based on ¹⁷⁷Lu-DOTATATE distributions of three voxel-based phantoms. SPECT projection images were simulated for 1, 24, 96, and 168 h post-administration and reconstructed with LundADose using 3D OS-EM reconstruction. Mono-exponential curves were fitted to the bio-kinetic data for the kidneys, liver, spleen, and tumours resulting in SPECT time-integrated activity (SPECT-TIA). The SPECT-TIA were used to compute mean absorbed doses using LundADose (LND-D_SPECT) and OLINDA (OLINDA-D_SPECT) for the organs. Pre-defined true activity images, were used to obtain TRUE-TIA and, together with full MC simulations, computed the true doses (MC-D_True). The dosimetry accuracy was assessed by comparing LND-D_SPECT and OLINDA-D_SPECT to MC-D_True.

Results

Overall, the results presented an overestimation of the mean absorbed dose by LND-D_SPECT compared to the MC-D_True with a dosimetry accuracy ≤6.6%. This was attributed to spill-out activity from the reconstructed LND-D_SPECT, resulting in a higher dose contribution than the MC-D_True. There was a general underestimation (<8.1%) of OLINDA-D_SPECT compared to MC-D_True attributed to the geometrical difference in shape between the voxel-based phantoms and the OLINDA models. Furthermore, OLINDA-D_SPECT considers self-doses while MC-D_True reflects self-doses plus cross-doses.

Conclusion

The better than 10% accuracy suggests that the mean dose values obtained with LND-D_SPECT and OLINDA-D_SPECT approximate the true values. The mean absorbed doses of the two software programs, and the gold standard were comparable. This work shall be of use for optimising ¹⁷⁷Lu dosimetry for clinical applications.

Precision dosimetry in yttrium-90 radioembolization through CT imaging of radiopaque microspheres in a rabbit liver model

PURPOSE

To perform precision dosimetry in yttrium-90 radioembolization through CT imaging of radiopaque microspheres in a rabbit liver model and to compare extracted dose metrics to those produced from conventional PET-based dosimetry.

MATERIALS AND METHODS

A CT calibration phantom was designed containing posts with nominal microsphere concentrations of 0.5 mg/mL, 5.0 mg/mL, and 25.0 mg/mL. The mean Hounsfield unit was extracted from the post volumes to generate a calibration curve to relate Hounsfield units to microsphere concentration. A nominal bolus of 40 mg of microspheres was administered to the livers of eight rabbits, followed by PET/CT imaging. A CT-based activity distribution was calculated through the application of the calibration curve to the CT liver volume. Post-treatment dosimetry was performed through the convolution of yttrium-90 dose-voxel kernels and the PET- and CT-based cumulated activity distributions. The mean dose to the liver in PET- and CT-based dose distributions was compared through linear regression, ANOVA, and Bland-Altman analysis.

RESULTS

A linear least-squares fit to the average Hounsfield unit and microsphere concentration data from the calibration phantom confirmed a strong correlation (r2 > 0.999) with a slope of 14.13 HU/mg/mL. A poor correlation was found between the mean dose derived from CT and PET (r2 = 0.374), while the ANOVA analysis revealed statistically significant differences (p < 10-12) between the MIRD-derived mean dose and the PET- and CT-derived mean dose. Bland-Altman analysis predicted an offset of 15.0 Gy between the mean dose in CT and PET. The dose within the liver was shown to be more heterogeneous in CT than in PET with an average coefficient of variation equal to 1.99 and 1.02, respectively.

CONCLUSION

The benefits of a CT-based approach to post-treatment dosimetry in yttrium-90 radioembolization include improved visualization of the dose distribution, reduced partial volume effects, a better representation of dose heterogeneity, and the mitigation of respiratory motion effects. Post-treatment CT imaging of radiopaque microspheres in yttrium-90 radioembolization provides the means to perform precision dosimetry and extract accurate dose metrics used to refine the understanding of the dose-response relationship, which could ultimately improve future patient outcomes.

Personalized Dosimetry in Targeted Radiation Therapy: A Look to Methods, Tools and Critical Aspects

Targeted radiation therapy (TRT) is a strategy increasingly adopted for the treatment of different types of cancer. The urge for optimization, as stated by the European Council Directive (2013/59/EURATOM), requires the implementation of a personalized dosimetric approach, similar to what already happens in external beam radiation therapy (EBRT). The purpose of this paper is to provide a thorough introduction to the field of personalized dosimetry in TRT, explaining its rationale in the context of optimization and describing the currently available methodologies. After listing the main therapies currently employed, the clinical workflow for the absorbed dose calculation is described, based on works of the most experienced authors in the literature and recent guidelines. Moreover, the widespread software packages for internal dosimetry are presented and critical aspects discussed. Overall, a selection of the most important and recent articles about this topic is provided.

Quantitative SPECT/CT for dosimetry of peptide receptor radionuclide therapy

Neuroendocrine tumors (NETs) are uncommon malignancies of increasing incidence and prevalence. As these slow growing tumors usually overexpress somatostatin receptors (SSTRs), the use of ⁶⁸Ga-DOTA-peptides (gallium-68 chelated with dodecane tetra-acetic acid to somatostatin), which bind to the SSTRs, allows for PET based imaging and selection of patients for peptide receptor radionuclide therapy (PRRT). PRRT with radiolabeled somatostatin analogues such as ¹⁷⁷Lu-DOTATATE (lutetium-177-[DOTA,Tyr³]-octreotate), is mainly used for the treatment of metastatic or inoperable NETs. However, PRRT is generally administered at a fixed injected activity in order not to exceed dose limits in critical organs, which is suboptimal given the variability in radiopharmaceutical uptake among patients. Advances in SPECT (single photon emission computed tomography) imaging enable the absolute quantitative measure of the true radiopharmaceutical distribution providing for PRRT dosimetry in each patient. Personalized PRRT based on patient-specific dosimetry could improve therapeutic efficacy by optimizing effective tumor absorbed dose while limiting treatment related radiotoxicity.

Overview of commercial treatment planning systems for targeted radionuclide therapy

INTRODUCTION

Targeted Radionuclide Therapy (TRT) is a branch of cancer medicine dealing with the therapeutic use of radioisotopes associated with biological vectors accumulating in the tumors/targets, indicated as Molecular Radiotherapy (MRT), or directly injected into the arteries that supply blood to liver tumour vasculature, indicated as Selective RT (SRT). The aim of this work is to offer a panoramic view on the increasing number of commercially-available TRT treatment planning systems (TPSs).

MATERIALS AND METHODS

A questionnaire was sent to manufacturers' representatives. Academic software were not considered. Questions were grouped as follows: general information, clinical workflow, calibration procedure, image processing/reconstruction, image registration and segmentation tools, time-activity curve (TAC) fitting and absorbed dose calculation.

RESULTS

All software reported have CE-marking. TPSs were divided between SRT-dedicated software [4] and MRT [5] dosimetry software. In SRT, since no kinetic process is involved, absorbed dose calculation does not require TAC fitting, and image registration is not fully developed in all TPS. All software requires a radionuclide-specific calibration. In SRT, a relative image calibration can be obtained by scaling the counts to a known activity. Automated VOI contouring and rigid/deformable propagation between different acquisitions time-points is implemented in most TPSs, although DICOM export is rare. Different TAC fits are available depending on the number of time-points. Voxel S-value and Local deposition methods are the most frequent dosimetric approaches; dose-voxel kernel convolution and semi-Monte Carlo method are also available.

CONCLUSIONS

Available TPSs allows performing personalized dosimetry in clinical practice. Individual variations in methodology/algorithms must be considered in the standardisation/harmonization processes.

Efficacy of voxel-based dosimetry map for predicting response to trans-arterial radioembolization therapy for hepatocellular carcinoma: a pilot study

OBJECTIVE

Typical clinical dosimetry models for trans-arterial radioembolization (TARE) assume uniform dose distribution in each tissue compartment. We performed simple voxel-based dosimetry using post-treatment 90Y PET following TARE with 90Y-resin microspheres and investigated its prognostic value in a pilot cohort.

METHOD

Ten patients with 14 hepatocellular carcinoma lesions who underwent TARE with 90Y-resin microspheres were retrospectively included. The partition model-based expected target tumor dose (TDp) was calculated using a pretreatment 99mTc-macroaggregated albumin scan. From post-treatment 90Y-microsphere PET and voxel-wise S-value kernels, voxel-based dose maps were produced and the absorbed dose of each lesion (TDv) was calculated. Heterogeneity of intratumoral absorbed doses was assessed using the SD and coefficient of variation of voxel doses. The response of each lesion was determined based on contrast-enhanced MRI or CT, or both. Lesion responses were classified as local control success or failure. Prognostic values of dosimetry parameters and clinicopathological factors were evaluated in terms of progression-free survival (PFS) of each lesion.

RESULTS

TDv was significantly different between local control success and failure groups, whereas tumor size, TDp and intratumoral dose heterogeneity were not. Univariate survival analysis identified serum aspartate transaminase level ≥40 IU/L, tumor size ≥66 mm and TDv <81 Gy as significant prognostic factors for PFS. However, only TDv was an independent predictive factor in the multivariate analysis (P = 0.022). There was a significant correlation between TDv and PFS (P = 0.009; r = 0.669).

CONCLUSIONS

In TARE, voxel-based dose index TDv can be estimated on post-treatment 90Y PET using a simple method. TDv was a more effective prognostic factor for TARE than TDp and clinicopathologic factors in this pilot study. Further studies are warranted on the role of voxel-based dose and dose distribution in TARE.

EANM dosimetry committee series on standard operational procedures: a unified methodology for 99mTc-MAA pre- and 90Y peri-therapy dosimetry in liver radioembolization with 90Y microspheres

The aim of this standard operational procedure is to standardize the methodology employed for the evaluation of pre- and post-treatment absorbed dose calculations in 90Y microsphere liver radioembolization. Basic assumptions include the permanent trapping of microspheres, the local energy deposition method for voxel dosimetry, and the patient-relative calibration method for activity quantification.The identity of 99mTc albumin macro-aggregates (MAA) and 90Y microsphere biodistribution is also assumed. The large observed discrepancies in some patients between 99mTc-MAA predictions and actual 90Y microsphere distributions for lesions is discussed. Absorbed dose predictions to whole non-tumoural liver are considered more reliable and the basic predictors of toxicity. Treatment planning based on mean absorbed dose delivered to the whole non-tumoural liver is advised, except in super-selective treatments.Given the potential mismatch between MAA simulation and actual therapy, absorbed doses should be calculated both pre- and post-therapy. Distinct evaluation between target tumours and non-tumoural tissue, including lungs in cases of lung shunt, are vital for proper optimization of therapy. Dosimetry should be performed first according to a mean absorbed dose approach, with an optional, but important, voxel level evaluation. Fully corrected 99mTc-MAA Single Photon Emission Computed Tomography (SPECT)/computed tomography (CT) and 90Y TOF PET/CT are regarded as optimal acquisition methodologies, but, for institutes where SPECT/CT is not available, non-attenuation corrected 99mTc-MAA SPECT may be used. This offers better planning quality than non dosimetric methods such as Body Surface Area (BSA) or mono-compartmental dosimetry. Quantitative 90Y bremsstrahlung SPECT can be used if dedicated correction methods are available.The proposed methodology is feasible with standard camera software and a spreadsheet. Available commercial or free software can help facilitate the process and improve calculation time.

Standardizing SPECT/CT dosimetry following radioembolization with yttrium-90 microspheres

Background

Multiple post-treatment dosimetry methods are currently under investigation for Yttrium-90 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{90}\hbox {Y}$$\end{document}90Y) radioembolization. Within each methodology, a variety of dosimetric inputs exists that affect the final dose estimates. Understanding their effects is essential to facilitating proper dose analysis and crucial in the eventual standardization of radioembolization dosimetry. The purpose of this study is to investigate the dose differences due to different self-calibrations and mass density assignments in the non-compartmental and local deposition methods. A practical mean correction method was introduced that permits dosimetry in images where the quality is compromised by patient motion and partial volume effects.

Methods

Twenty-one patients underwent \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{90}\hbox {Y}$$\end{document}90Y radioembolization and were imaged with SPECT/CT. Five different self-calibrations (FOV, Body, OAR, Liverlung, and Liver) were implemented and dosimetrically compared. The non-compartmental and local deposition method were used to perform dosimetry based on either nominal- or CT calibration-based mass densities. A mean correction method was derived assuming homogeneous densities. Cumulative dose volume histograms, linear regressions, boxplots, and Bland Altman plots were utilized for analysis.

Results

Up to 270% weighted dose difference was found between self-calibrations with mean dose differences up to 50 Gy in the liver and 23 Gy in the lungs. Between the local deposition and non-compartmental methods, the liver and lung had dose differences within 0.71 Gy and 20 Gy, respectively. The local deposition method’s nominal and CT calibration-based mass density implementations dosimetric metrics were within 1.4% in the liver and 24% in the lungs. The mean lung doses calculated with the CT method were shown to be inflated. The mean correction method demonstrated that the corrected mean doses were greater by up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 5$$\end{document}∼5 Gy in the liver and lower by up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 12$$\end{document}∼12 Gy in the lungs.

Conclusions

The OAR calibration may be utilized as a potentially more accurate and precise self-calibration. The non-compartmental method was found more comparable to the local deposition method in organs that were more homogeneous in mass densities. Due to the potential for inflated lung mean doses, the non-compartmental and local deposition method implemented with nominal mass densities is recommended for more consistent dosimetric results. If patient motion and partial volume effects are present in the liver, our practical correction method will calculate more representative doses in images suboptimal for dosimetry.

Yttrium-90 TOF-PET-Based EUD Predicts Response Post Liver Radioembolizations Using Recommended Manufacturer FDG Reconstruction Parameters

Purpose

Explaining why ⁹⁰Y TOF-PET based equivalent uniform dose (EUD) using recommended manufacturer FDG reconstruction parameters has been shown to predict response.

Methods

The hot rods insert of a Jaszczak deluxe phantom was partially filled with a 2.65 GBq ⁹⁰Y - 300ml DTPA water solution resulting in a 100 Gy mean absorbed dose in the 6 sectors. A two bed 20min/position acquisition was performed on a 550ps- and on a 320ps- TOF-PET/CT and reconstructed with recommended manufacturer FDG reconstruction parameters, without and with additional filtering. The whole procedure was repeated on both PET after adding 300ml of water (50Gy setup). The phantom was acquired again after decay by a factor of 10 (5Gy setup), but with 200min per bed position. For comparison, the phantom was also acquired with ¹⁸F activity corresponding to a clinical FDG whole body acquisition.

Results

The 100Gy-setup provided a hot rod sectors image almost as good as the ¹⁸F phantom. However, despite acquisition time compensation, the 5Gy-setup provides much lower quality imaging. TOF-PET based sectors EUDs for the three large rod sectors agreed with the actual EUDs computed with a radiosensitivity of 0.021Gy^-1 well in the range observed in external beam radiotherapy (EBRT), i.e. 0.01-0.04Gy^-1. This agreement explains the reunification of the dose-response relationships of the glass and resin spheres in HCC using the TOF-PET based EUD. Additional filtering reduced the EUDs agreement quality.

Conclusions

Recommended manufacturer FDG reconstruction parameters are suitable in TOF-PET post ⁹⁰Y liver radioembolization for accurate tumour EUD computation. The present results rule out the use of low specific activity phantom studies to optimize reconstruction parameters.

Lung shunt fraction calculation using 99mTc-MAA SPECT/CT imaging for 90Y microsphere selective internal radiation therapy of liver tumors

BACKGROUND

^99mTc-macroaggregated albumin (^99mTc-MAA) scintigraphy is utilized in treatment planning for Yttrium-90 (⁹⁰Y) Selective Internal Radiation Therapy (SIRT) of liver tumors to evaluate hepatopulmonary shunting by calculating the lung shunt fraction (LSF). The purpose of this study was to evaluate if LSF calculation using SPECT/CT instead of planar gamma camera imaging is more accurate and if this can potentially lead to more effective treatment planning of hepatic lesions while avoiding excessive pulmonary irradiation.

RESULTS

LSF calculation was obtained using two different methodologies in 85 cases from consecutive patients intended to receive ⁹⁰Y SIRT. The first method was based on planar gamma camera imaging in the anterior and posterior views with geometric mean calculation of the LSF from regions of interest of the liver and lungs. The second method was based on segmentation of the liver and lungs from SPECT/CT images of the thorax and abdomen. The differences in planar imaging versus SPECT/CT derived LSF values along with the estimated absorbed lung mean dose (LMD) were evaluated. The LSF values were higher in planar imaging versus SPECT/CT in 81/85 cases, with a mean value of 8.5% vs. 4.6% respectively; the difference was statistically significant using a paired t-test (alpha = 0.05). In those patients who received SIRT, the estimated absorbed LMD calculated with planar imaging was significantly higher than with SPECT/CT (t-test, P < 0.005). Repeated phantom experiments using an anthropomorphic torso phantom with variable ^99mTc activity concentrations for the liver and lungs were performed with the standard patient protocol, demonstrated improved accuracy of the LSF calculation based on SPECT/CT than planar imaging (mean overestimated value of 6% vs. 26%).

CONCLUSIONS

This study demonstrates that LSF calculation using planar imaging can be significantly overestimated while calculation using SPECT/CT imaging and appropriate segmentation tools can be more accurate. Minimizing the errors in obtaining the LSF can lead to more effective ⁹⁰Y SIRT treatment planning for hepatic tumors while ensuring the lung dose will not exceed the standard acceptable safety thresholds.

Yttrium-90 quantitative phantom study using digital photon counting PET

BACKGROUND

PET imaging of 90Y-microsphere distribution following radioembolisation is challenging due to the count-starved statistics from the low branching ratio of e+/e- pair production during 90Y decay. PET systems using silicon photo-multipliers have shown better 90Y image quality compared to conventional photo-multiplier tubes. The main goal of the present study was to evaluate reconstruction parameters for different phantom configurations and varying listmode acquisition lengths to improve quantitative accuracy in 90Y dosimetry, using digital photon counting PET/CT.

METHODS

Quantitative PET and dosimetry accuracy were evaluated using two uniform cylindrical phantoms specific for PET calibration validation. A third body phantom with a 9:1 hot sphere-to-background ratio was scanned at different activity concentrations of 90Y. Reconstructions were performed using OSEM algorithm with varying parameters. Time-of-flight and point-spread function modellings were included in all reconstructions. Absorbed dose calculations were carried out using voxel S-values convolution and were compared to reference Monte Carlo simulations. Dose-volume histograms and root-mean-square deviations were used to evaluate reconstruction parameter sets. Using listmode data, phantom and patient datasets were rebinned into various lengths of time to assess the influence of count statistics on the calculation of absorbed dose. Comparisons between the local energy deposition method and the absorbed dose calculations were performed.

RESULTS

Using a 2-mm full width at half maximum post-reconstruction Gaussian filter, the dosimetric accuracy was found to be similar to that found with no filter applied but also reduced noise. Larger filter sizes should not be used. An acquisition length of more than 10 min/bed reduces image noise but has no significant impact in the quantification of phantom or patient data for the digital photon counting PET. 3 iterations with 10 subsets were found suitable for large spheres whereas 1 iteration with 30 subsets could improve dosimetry for smaller spheres.

CONCLUSION

The best choice of the combination of iterations and subsets depends on the size of the spheres. However, one should be careful on this choice, depending on the imaging conditions and setup. This study can be useful in this choice for future studies for more accurate 90Y post-dosimetry using a digital photon counting PET/CT.

Impact of contouring methods on pre-treatment and post-treatment dosimetry for the prediction of tumor control and survival in HCC patients treated with selective internal radiation therapy

INTRODUCTION

The aim of this study was to evaluate the impact of the contouring methods on dose metrics and their predictive value on tumor control and survival, in both situations of pre-treatment and post-treatment dosimetry, for patients with advanced HCC treated with SIRT.

METHODS

Forty-eight patients who underwent SIRT between 2012 and 2020 were retrospectively included in this study. Target volumes were delineated using two methods: MRI-based contours manually drawn by a radiologist and then registered on SPECT/CT and PET/CT via deformable registration (Pre-C_MRI and Post-C_MRI), ^99mTc-MAA-SPECT and ⁹⁰Y-microspheres-PET 10% threshold contouring (Pre-C_SPECT and Post-C_PET). The mean absorbed dose (Dm) and the minimal absorbed dose delivered to 70% of the tumor volume (D70) were evaluated with both contouring methods; the tumor-to-normal liver uptake ratio (TNR) was evaluated with MRI-based contours only. Tumor response was assessed using the mRECIST criteria on the follow-up MRIs.

RESULTS

No significant differences were found for Dm and TNR between pre- and post-treatment. TNR evaluated with radiologic contours (Pre-C_MRI and Post-C_MRI) were predictive of tumor control at 6 months on pre- and post-treatment dosimetry (OR 5.9 and 7.1, respectively; p = 0.02 and 0.01). All dose metrics determined with both methods were predictive of overall survival (OS) on pre-treatment dosimetry, but only Dm with MRI-based contours was predictive of OS on post-treatment images with a median of 23 months for patients with a supramedian Dm versus 14 months for the others (p = 0.04).

CONCLUSION

In advanced HCC treated with SIRT, Dm and TNR determined with radiologic contours were predictive of tumor control and OS. This study shows that a rigorous clinical workflow (radiologic contours + registration on scintigraphic images) is feasible and should be prospectively considered for improving therapeutic strategy.

Relationship of radiation dose to efficacy of radioembolization of liver metastasis from breast cancer

PURPOSE

To determine the relationship of tumoral and nontumoral radiation dose to response and toxicity after transarterial radioembolization (TARE) of breast cancer liver metastasis.

METHODS

This retrospective study evaluated all patients with breast cancer liver metastases treated with TARE (2/2011-6/2019). Extent of disease was measured as unilobar or bilobar on baseline PET/CT prior to TARE. Response was assessed for targeted regions with modified PERCIST criteria on first follow-up PET/CT. Tumoral and nontumoral liver dosimetry was evaluated by performing volumetric segmentation on post-TARE Bremsstrahlung SPECT/CT. ≥Grade 3 hepatotoxicity was defined as ≥grade 3 bilirubin/AST/ALT elevation or ascites requiring intervention. Fisher's exact tests, Wilcoxon rank sum tests, and Kaplan-Meier survival analysis were performed.

RESULTS

Among 64 women, 60 patients had pre- and post-TARE PET/CT, of whom 46/60 (77 %) achieved objective response (OR). Responders received higher tumoral dose with a median (interquartile range) of 167 (96-217) vs. 54 (45-62) Gy (p < 0.001). ≥Grade 3 hepatotoxicity occurred in 8/64 (12.5 %) and was associated with higher pre-treatment bilirubin levels of 0.9 (0.9-1.1) vs. 0.5 (0.4-0.7) mg/dL (p = 0.013). Median overall survival (OS) was 11 (95 % CI 10-19) months. Bilobar disease (Hazard Ratio [HR]: 2.77, 95 % CI 1.11-6.89, p = 0.028) and elevated pre-TARE AST (HR 1.02, 95 % CI 1.01-1.03, p < 0.001) were independently associated with shorter survival. ≥Grade 3 hepatotoxicity was associated with reduced survival (p < 0.001). OR was associated with longer OS of 17 months, compared with 10 months (p = 0.027).

CONCLUSION

In TARE for breast cancer liver metastasis, higher tumoral radiation dose (>79.5 Gy) was associated with OR, which was associated with longer survival. Pre-existing liver dysfunction was associated with hepatotoxicity, which was associated with decreased survival.

Clinical implementation of PLANET® Dose for dosimetric assessment after [177Lu]Lu-DOTA-TATE: comparison with Dosimetry Toolkit® and OLINDA/EXM® V1.0

Background

The aim of this study was to compare a commercial dosimetry workstation (PLANET® Dose) and the dosimetry approach (GE Dosimetry Toolkit® and OLINDA/EXM® V1.0) currently used in our department for quantification of the absorbed dose (AD) to organs at risk after peptide receptor radionuclide therapy with [¹⁷⁷Lu]Lu-DOTA-TATE.

Methods

An evaluation on phantom was performed to determine the SPECT calibration factor variations over time and to compare the Time Integrated Activity Coefficients (TIACs) obtained with the two approaches. Then, dosimetry was carried out with the two tools in 21 patients with neuroendocrine tumours after the first and second injection of 7.2 ± 0.2 GBq of [¹⁷⁷Lu]Lu-DOTA-TATE (40 dosimetry analyses with each software). SPECT/CT images were acquired at 4 h, 24 h, 72 h and 192 h post-injection and were reconstructed using the Xeleris software (General Electric). The liver, spleen and kidneys masses and TIACs were determined using Dosimetry Toolkit® (DTK) and PLANET® Dose. The ADs were calculated using OLINDA/EXM® V1.0 and the Local Deposition Method (LDM) or Dose voxel-Kernel convolution (DK) on PLANET® Dose.

Results

With the phantom, the 3D calibration factors showed a slight variation (0.8% and 3.3%) over time, and TIACs of 225.19 h and 217.52 h were obtained with DTK and PLANET® Dose, respectively. In patients, the root mean square deviation value was 8.9% for the organ masses, 8.1% for the TIACs, and 9.1% and 7.8% for the ADs calculated with LDM and DK, respectively. The Lin’s concordance correlation coefficient was 0.99 and the Bland–Altman plot analysis estimated that the AD value difference between methods ranged from − 0.75 to 0.49 Gy, from − 0.20 to 0.64 Gy, and from − 0.43 to 1.03 Gy for 95% of the 40 liver, kidneys and spleen dosimetry analyses. The dosimetry method had a minor influence on AD differences compared with the image registration and organ segmentation steps.

Conclusions

The ADs to organs at risk obtained with the new workstation PLANET® Dose are concordant with those calculated with the currently used software and in agreement with the literature. These results validate the use of PLANET® Dose in clinical routine for patient dosimetry after targeted radiotherapy with [¹⁷⁷Lu]Lu-DOTA-TATE.

Dosimetry for Optimized, Personalized Radiopharmaceutical Therapy

Radiopharmaceutical therapy (RPT) has grown rapidly over the last decade for treatment of numerous cancer types. Dosimetric guidance, as with other radiotherapy modalities, has benefitted patients by reducing the incidence of side effects and improving overall survival in populations treated under this paradigm. Development of tools and techniques for dosimetry-guided therapy is ongoing, with numerous the Food and Drug Administration-cleared products reaching the U.S. market in 2019. Safe use of commercial dosimetry platforms requires a deep understanding of the underlying physical principles and thoroughly vetted input data. Likewise, interpretation of dosimetry results relies on an understanding of radiobiological principles, and the principles of uncertainty propagation. In this article, we review strategies commonly employed for dosimetry-guided RPT - including quantitative imaging, dose calculation methods, and modeling of dose across time-points. Additionally, we review recent literature evidence (2013-2020) demonstrating the efficacy of personalized RPT.

IPEM topical report: current molecular radiotherapy service provision and guidance on the implications of setting up a dosimetry service

Despite a growth in molecular radiotherapy treatment (MRT) and an increase in interest, centres still rarely perform MRT dosimetry. The aims of this report were to assess the main reasons why centres are not performing MRT dosimetry and provide advice on the resources required to set-up such a service. A survey based in the United Kingdom was developed to establish how many centres provide an MRT dosimetry service and the main reasons why it is not commonly performed. Twenty-eight per cent of the centres who responded to the survey performed some form of dosimetry, with 88% of those centres performing internal dosimetry. The survey showed that a 'lack of clinical evidence', a 'lack of guidelines' and 'not current UK practice' were the largest obstacles to setting up an MRT dosimetry service. More practical considerations, such as 'lack of software' and 'lack of staff training/expertise', were considered to be of lower significance by the respondents. Following on from the survey, this report gives an overview of the current guidelines, and the evidence available demonstrating the benefits of performing MRT dosimetry. The resources required to perform such techniques are detailed with reference to guidelines, training resources and currently available software. It is hoped that the information presented in this report will allow MRT dosimetry to be performed more frequently and in more centres, both in routine clinical practice and in multicentre trials. Such trials are required to harmonise dosimetry techniques between centres, build on the current evidence base, and provide the data necessary to establish the dose-response relationship for MRT.

Impact of the dosimetry approach on the resulting 90Y radioembolization planned absorbed doses based on 99mTc-MAA SPECT-CT: is there agreement between dosimetry methods?

BACKGROUND

Prior radioembolization, a simulation using ^99mTc-macroaggregated albumin as ⁹⁰Y-microspheres surrogate is performed. Gamma scintigraphy images (planar, SPECT, or SPECT-CT) are acquired to evaluate intrahepatic ⁹⁰Y-microspheres distribution and detect possible extrahepatic and lung shunting. These images may be used for pre-treatment dosimetry evaluation to calculate the ⁹⁰Y activity that would get an optimal tumor response while sparing healthy tissues. Several dosimetry methods are available, but there is still no consensus on the best methodology to calculate absorbed doses. The goal of this study was to retrospectively evaluate the impact of using different dosimetry approaches on the resulting ⁹⁰Y-radioembolization pre-treatment absorbed dose evaluation based on ^99mTc-MAA images.

METHODS

Absorbed doses within volumes of interest resulting from partition model (PM) and 3D voxel dosimetry methods (3D-VDM) (dose-point kernel convolution and local deposition method) were evaluated. Additionally, a new "Multi-tumor Partition Model" (MTPM) was developed. The differences among dosimetry approaches were evaluated in terms of mean absorbed dose and dose volume histograms within the volumes of interest.

RESULTS

Differences in mean absorbed dose among dosimetry methods are higher in tumor volumes than in non-tumoral ones. The differences between MTPM and both 3D-VDM were substantially lower than those observed between PM and any 3D-VDM. A poor correlation and concordance were found between PM and the other studied dosimetry approaches. DVH obtained from either 3D-VDM are pretty similar in both healthy liver and individual tumors. Although no relevant global differences, in terms of absorbed dose in Gy, between both 3D-VDM were found, important voxel-by-voxel differences have been observed.

CONCLUSIONS

Significant differences among the studied dosimetry approaches for ⁹⁰Y-radioembolization treatments exist. Differences do not yield a substantial impact in treatment planning for healthy tissue but they do for tumoral liver. An individual segmentation and evaluation of the tumors is essential. In patients with multiple tumors, the application of PM is not optimal and the 3D-VDM or the new MTPM are suggested instead. If a 3D-VDM method is not available, MTPM is the best option. Furthermore, both 3D-VDM approaches may be indistinctly used.

3D image-based dosimetry for Yttrium-90 radioembolization of hepatocellular carcinoma: Impact of imaging method on absorbed dose estimates

BACKGROUND

To improve therapy outcome of Yttrium-90 selective internal radiation therapy (⁹⁰Y SIRT), patient-specific post-therapeutic dosimetry is required. For this purpose, various dosimetric approaches based on different available imaging data have been reported. The aim of this work was to compare post-therapeutic 3D absorbed dose images using Technetium-99m (^99mTc) MAA SPECT/CT, Yttrium-90 (⁹⁰Y) bremsstrahlung (BRS) SPECT/CT, and ⁹⁰Y PET/CT.

METHODS

Ten SIRTs of nine patients with unresectable hepatocellular carcinoma (HCC) were investigated. The ^99mTc SPECT/CT data, obtained from ^99mTc-MAA-based treatment simulation prior to ⁹⁰Y SIRT, were scaled with the administered ⁹⁰Y therapy activity. 3D absorbed dose images were generated by dose kernel convolution with scaled ^99mTc/⁹⁰Y SPECT/CT, ⁹⁰Y BRS SPECT/CT, and ⁹⁰Y PET/CT data of each patient. Absorbed dose estimates in tumor and healthy liver tissue obtained using the two SPECT/CT methods were compared against ⁹⁰Y PET/CT.

RESULTS

The percentage deviation of tumor absorbed dose estimates from ⁹⁰Y PET/CT values was on average -2 ± 18% for scaled ^99mTc/⁹⁰Y SPECT/CT, whereas estimates from ⁹⁰Y BRS SPECT/CT differed on average by -50 ± 13%. For healthy liver absorbed dose estimates, all three imaging methods revealed comparable values.

CONCLUSION

The quantification capabilities of the imaging data influence ⁹⁰Y SIRT tumor dosimetry, while healthy liver absorbed dose values were comparable for all investigated imaging data. When no ⁹⁰Y PET/CT image data are available, the proposed scaled ^99mTc/⁹⁰Y SPECT/CT dosimetry method was found to be more appropriate for HCC tumor dosimetry than ⁹⁰Y BRS SPECT/CT based dosimetry.

Correlation and Agreement of Yttrium-90 Positron Emission Tomography/Computed Tomography with Ex Vivo Radioembolization Microsphere Deposition in the Rabbit VX2 Liver Tumor Model

PURPOSE

To demonstrate a stronger correlation and agreement of yttrium-90 (⁹⁰Y) positron emission tomography (PET)/computed tomography (CT) measurements with explant liver tumor dosing compared with the standard model (SM) for radioembolization.

MATERIALS AND METHODS

Hepatic VX2 tumors were implanted into New Zealand white rabbits, with growth confirmed by 7 T magnetic resonance imaging. Seventeen VX2 rabbits provided 33 analyzed tumors. Treatment volumes were calculated from manually drawn volumes of interest (VOI) with three-dimensional surface renderings. Radioembolization was performed with glass ⁹⁰Y microspheres. PET/CT imaging was completed with scatter and attenuation correction. Three-dimensional ellipsoid VOI were drawn to encompass tumors on fused images. Tumors and livers were then explanted for inductively coupled plasma (ICP)-optical emission spectroscopy (OES) analysis of microsphere content. ⁹⁰Y PET/CT and SM measurements were compared with reference standard ICP-OES measurements of tumor dosing with Pearson correlation and Bland-Altman analyses for agreement testing with and without adjustment for tumor necrosis.

RESULTS

The median infused activity was 33.3 MBq (range, 5.9-152.9). Tumor dose was significantly correlated with ⁹⁰Y PET/CT measurements (r = 0.903, P < .001) and SM estimates (r = 0.607, P < .001). Bland-Altman analyses showed that the SM tended to underestimate the tumor dosing by a mean of -8.5 Gy (CI, -26.3-9.3), and the degree of underestimation increased to a mean of -18.3 Gy (CI, -38.5-1.9) after the adjustment for tumor necrosis.

CONCLUSIONS

⁹⁰Y PET/CT estimates were strongly correlated and had better agreement with reference measurements of tumor dosing than SM estimates.

The Impact of Radiobiologically Informed Dose Prescription on the Clinical Benefit of 90Y SIRT in Colorectal Cancer Patients

The purpose of this study was to establish the dose-response relationship of selective internal radiation therapy (SIRT) in patients with metastatic colorectal cancer (mCRC), when informed by radiobiological sensitivity parameters derived from mCRC cell lines exposed to 90Y. Methods: Twenty-three mCRC patients with liver metastases refractory to chemotherapy were included. 90Y bremsstrahlung SPECT images were transformed into dose maps assuming the local dose deposition method. Baseline and follow-up CT scans were segmented to derive liver and tumor volumes. Mean, median, and D70 (minimum dose to 70% of tumor volume) values determined from dose maps were correlated with change in tumor volume and volumetric RECIST response using linear and logistic regression, respectively. Radiosensitivity parameters determined by clonogenic assays of mCRC cell lines HT-29 and DLD-1 after exposure to 90Y or external beam radiotherapy (EBRT; 6 MV photons) were used in biologically effective dose (BED) calculations. Results: Mean administered radioactivity was 1,469 ± 428 MBq (range, 847-2,185 MBq), achieving a mean absorbed radiation dose to tumor of 35.5 ± 9.4 Gy and mean normal liver dose of 26.4 ± 6.8 Gy. A 1.0 Gy increase in mean, median, and D70 absorbed dose was associated with a reduction in tumor volume of 1.8%, 1.8%, and 1.5%, respectively, and an increased probability of a volumetric RECIST response (odds ratio, 1.09, 1.09, and 1.10, respectively). Threshold mean, median and D70 doses for response were 48.3, 48.8, and 41.8 Gy, respectively. EBRT-equivalent BEDs for 90Y are up to 50% smaller than those calculated by applying protraction-corrected radiobiological parameters derived from EBRT alone. Conclusion: Dosimetric studies have assumed equivalence between 90Y SIRT and EBRT, leading to inflation of BED for SIRT and possible undertreatment. Radiobiological parameters for 90Y were applied to a BED model, providing a calculation method that has the potential to improve assessment of tumor control.

Planning dosimetry for 90 Y radioembolization with glass microspheres: Evaluating the fidelity of 99m Tc-MAA and partition model predictions

PURPOSE

^99m Tc-MAA-SPECT/CT may be used in ⁹⁰ Y-glass microsphere radioembolization treatment planning to assess perfused liver volumes and absorbed dose distributions. The partition model (PM) offers a more detailed planning dosimetry option beyond the single-compartment model more traditionally used in ⁹⁰ Y radioembolization. As ⁹⁰ Y radioembolization treatments shift toward activities and doses that aim to achieve tumor control, accurate and reliable treatment planning dosimetry for both tumors and normal liver (NL) becomes more critical. In this work, we explore the accuracy and precision of ⁹⁰ Y dosimetry predictions from pretherapy ^99m Tc-MAA and PM.

METHODS

Both PM and voxel dosimetry models were used to calculate tumor and NL mean doses using both planning ^99m Tc-MAA and verification ⁹⁰ Y-SPECT/CT in this retrospective analysis of hepatocellular carcinoma cases treated with glass microspheres (NCT01900002, n = 32). Linear regression models were developed at first access, and then later correct, the estimates by (a) ^99m Tc-MAA for ⁹⁰ Y voxel dosimetry and (b) ^99m Tc-MAA PM for voxel dosimetry, separately for both tumors and NL. Bland-Altman analysis was then used to evaluate the accuracy and precision of the regression model predictions with the mean bias and 95% prediction intervals (PI, ±1.96σ). Two categories of cases were stratified (catheter matched vs catheter unmatched) by establishing the level of ^99m Tc-MAA and ⁹⁰ Y catheter position alignment. Only catheter-matched cases were included in the ^99m Tc-MAA vs ⁹⁰ Y voxel dosimetry comparison, while all cases were used to compare dosimetry models (PM vs voxel).

RESULTS

Half (16/32) of cases were deemed catheter matched. ^99m Tc-MAA could reliably predict NL doses in catheter-matched cases after application of the linear model, with mean bias (PI) of -1% (±31%). PM was equivalent to voxel dosimetry for NL doses with mean bias (PI) of 0% (±1%). Even among catheter-matched cases, ^99m Tc-MAA planning for ⁹⁰ Y tumor voxel doses was poor, overestimating dose by an average of nearly 40%. Upon application of the linear model, ^99m Tc-MAA predictions for ⁹⁰ Y tumor voxel dose were only minimally biased (-4%) but possessed very large PI (±104%). PM predictions for tumor voxel dose using the linear model also showed small bias (-6%) but maintained similarly high PI of ±90%. Cases with tumors representing a large majority (>80%) of the total tumor volume demonstrated the best scenarios for ^99m Tc-MAA and PM tumor dose predictions, with mean biases (PI) of -3% (±53%) and -4% (±21%), respectively.

CONCLUSION

The unconditional use of ^99m Tc-MAA to predict ⁹⁰ Y dosimetry across all cases is not recommended due to: (a) demonstrated the risk of unmatched catheter positions between procedures, and (b) large bias and uncertainty in ^99m Tc-MAA predictions in cases with matched catheter locations. However, NL voxel dose predictions with ^99m Tc-MAA are clinically viable and either PM or voxel dosimetry can be used to produce equivalent predictions. Both ^99m Tc-MAA and PM can provide tumor dose predictions with potential clinical utility, but only in catheter-matched cases and with tumors comprising a clear majority (>80%) of the total tumor volume. These findings stratify the predictive fidelity of ^99m Tc-MAA- and PM-based treatment planning for ⁹⁰ Y dosimetry in improving treatment outcomes.

Yttrium-90 radioembolization as a possible new treatment for brain cancer: proof of concept and safety analysis in a canine model

Purpose

To evaluate the safety, feasibility, and preliminary efficacy of yttrium-90 (⁹⁰Y) radioembolization (RE) as a minimally invasive treatment in a canine model with presumed spontaneous brain cancers.

Materials

Three healthy research dogs (R1–R3) and five patient dogs with spontaneous intra-axial brain masses (P1–P5) underwent cerebral artery RE with ⁹⁰Y glass microspheres (TheraSphere). ⁹⁰Y-RE was performed on research dogs from the unilateral internal carotid artery (ICA), middle cerebral artery (MCA), and posterior cerebral artery (PCA) while animals with brain masses were treated from the ICA. Post-treatment ⁹⁰Y PET/CT was performed along with serial neurological exams by a veterinary neurologist. One month after treatment, research dogs were euthanized and the brains were extracted and sent for microdosimetric and histopathologic analyses. Patient dogs received post-treatment MRI at 1-, 3-, and 6-month intervals with long-term veterinary follow-up.

Results

The average absorbed dose to treated tissue in R1–R3 was 14.0, 30.9, and 73.2 Gy, respectively, with maximum doses exceeding 1000 Gy. One month after treatment, research dog pathologic analysis revealed no evidence of cortical atrophy and rare foci consistent with chronic infarcts, e.g., < 2-mm diameter. Absorbed doses to masses in P1–P5 were 45.5, 57.6, 58.1, 45.4, and 64.1 Gy while the dose to uninvolved brain tissue was 15.4, 27.6, 19.2, 16.7, and 33.3 G, respectively. Among both research and patient animals, 6 developed acute neurologic deficits following treatment. However, in all surviving dogs, the deficits were transient resolving between 7 and 33 days post-therapy. At 1 month post-therapy, patient animals showed a 24–94% reduction in mass volume with partial response in P1, P3, and P4 at 6 months post-treatment. While P2 initially showed a response, by 5 months, the mass had advanced beyond pre-treatment size, and the dog was euthanized.

Conclusion

This proof of concept demonstrates the technical feasibility and safety of ⁹⁰Y-RE in dogs, while preliminary, initial data on the efficacy of ⁹⁰Y-RE as a potential treatment for brain cancer is encouraging.

Quantitative comparison of pre-treatment predictive and post-treatment measured dosimetry for selective internal radiation therapy using cone-beam CT for tumor and liver perfusion territory definition

Background

Selective internal radiation therapy (SIRT) is a promising treatment for unresectable hepatic malignancies. Predictive dose calculation based on a simulation using ^99mTc-labeled macro-aggregated albumin (^99mTc-MAA) before the treatment is considered as a potential tool for patient-specific treatment planning. Post-treatment dose measurement is mainly performed to confirm the planned absorbed dose to the tumor and non-tumor liver volumes. This study compared the predicted and measured absorbed dose distributions.

Methods

Thirty-one patients (67 tumors) treated by SIRT with resin microspheres were analyzed. Predicted and delivered absorbed dose was calculated using ^99mTc-MAA-SPECT and ⁹⁰Y-TOF-PET imaging. The voxel-level dose distribution was derived using the local deposition model. Liver perfusion territories and tumors have been delineated on contrast-enhanced CBCT images, which have been acquired during the ^99mTc-MAA work-up. Several dose-volume histogram (DVH) parameters together with the mean dose for liver perfusion territories and non-tumoral and tumoral compartments were evaluated.

Results

A strong correlation between the predicted and measured mean dose for non-tumoral volume was observed (r = 0.937). The ratio of measured and predicted mean dose to this volume has a first, second, and third interquartile range of 0.83, 1.05, and 1.25. The difference between the measured and predicted mean dose did not exceed 11 Gy. The correlation between predicted and measured mean dose to the tumor was moderate (r = 0.623) with a mean difference of − 9.3 Gy. The ratio of measured and predicted tumor mean dose had a median of 1.01 with the first and third interquartile ranges of 0.58 and 1.59, respectively. Our results suggest that ^99mTc-MAA-based dosimetry could predict under or over dosing of the non-tumoral liver parenchyma for almost all cases. For more than two thirds of the tumors, a predictive absorbed dose correctly indicated either good tumor dose coverage or under-dosing of the tumor.

Conclusion

Our results highlight the predictive value of ^99mTc-MAA-based dose estimation to predict non-tumor liver irradiation, which can be applied to prescribe an optimized activity aiming at avoiding liver toxicity. Compared to non-tumoral tissue, a poorer agreement between predicted and measured absorbed dose is observed for tumors.

Comparison of commercial dosimetric software platforms in patients treated with 177 Lu-DOTATATE for peptide receptor radionuclide therapy

Purpose

The aim of this study was to quantitatively compare five commercial dosimetric software platforms based on the analysis of clinical datasets of patients who benefited from peptide receptor radionuclide therapy (PRRT) with ¹⁷⁷Lu‐DOTATATE (LUTATHERA^®).

Methods

The dosimetric analysis was performed on two patients during two cycles of PRRT with ¹⁷⁷Lu. Single photon emission computed tomography/computed tomography images were acquired at 4, 24, 72, and 192 h post injection. Reconstructed images were generated using Dosimetry Toolkit^® (DTK) from Xeleris™ and HybridRecon‐Oncology version_1.3_Dicom (HROD) from HERMES. Reconstructed images using DTK were analyzed using the same software to calculate time‐integrated activity coefficients (TIAC), and mean absorbed doses were estimated using OLINDA/EXM V1.0 with mass correction. Reconstructed images from HROD were uploaded into PLANET® OncoDose from DOSIsoft, STRATOS from Phillips, Hybrid Dosimetry Module™ from HERMES, and SurePlan™ MRT from MIM. Organ masses, TIACs, and mean absorbed doses were calculated from each application using their recommendations.

Results

The majority of organ mass estimates varied by <9.5% between all platforms. The highest variability for TIAC results between platforms was seen for the kidneys (28.2%) for the two patients and the two treatment cycles. Relative standard deviations in mean absorbed doses were slightly higher compared with those observed for TIAC, but remained of the same order of magnitude between all platforms.

Conclusions

When applying a similar processing approach, results obtained were of the same order of magnitude regardless of the platforms used. However, the comparison of the performances of currently available platforms is still difficult as they do not all address the same parts of the dosimetric analysis workflow. In addition, the way in which data are handled in each part of the chain from data acquisition to absorbed doses may be different, which complicates the comparison exercise. Therefore, the dissemination of commercial solutions for absorbed dose calculation calls for the development of tools and standards allowing for the comparison of the performances between dosimetric software platforms.

Quantitative and Qualitative Improvement of Low-Count [68Ga]Citrate and [90Y]Microspheres PET Image Reconstructions Using Block Sequential Regularized Expectation Maximization Algorithm

PURPOSE

There are several important positron emission tomography (PET) imaging scenarios that require imaging with very low photon statistics, for which both quantitative accuracy and visual quality should not be neglected. For example, PET imaging with the low photon statistics is closely related to active efforts to significantly reduce radiation exposure from radiopharmaceuticals. We investigated two examples of low-count PET imaging: (a) imaging [90Y]microsphere radioembolization that suffers the very small positron emission fraction of Y-90's decay processes, and (b) cancer imaging with [68Ga]citrate with uptake time of 3-4 half-lives, necessary for visualizing tumors. In particular, we investigated a type of penalized likelihood reconstruction algorithm, block sequential regularized expectation maximization (BSREM), for improving both image quality and quantitative accuracy of these low-count PET imaging cases.

PROCEDURES

The NEMA/IEC Body phantom filled with aqueous solution of Y-90 or Ga-68 was scanned to mimic the low-count scenarios of corresponding patient data acquisitions on a time-of-flight (TOF) PET/magnetic resonance imaging system. Contrast recovery, background variation, and signal-to-noise ratio were evaluated in different sets of count densities using both conventional TOF ordered subset expectation (TOF-OSEM) and TOF-BSREM algorithms. The regularization parameter, beta, in BSREM that controls the tradeoff between image noise and resolution was evaluated to find a value for improved confidence in image interpretation. Visual quality assessment of the images obtained from patients administered with [68Ga]citrate (n = 6) was performed. We also made preliminary visual image quality assessment for one patient with [90Y]microspheres. In Y-90 imaging, the effect of 511-keV energy window selection for minimizing the number of random events was also evaluated.

RESULTS

Quantitatively, phantom images reconstructed with TOF-BSREM showed improved contrast recovery, background variation, and signal-to-noise ratio values over images reconstructed with TOF-OSEM. Both phantom and patient studies of delayed imaging of [68Ga]citrate show that TOF-BSREM with beta = 500 gives the best tradeoff between image noise and image resolution based on visual assessment by the readers. The NEMA-IQ phantom study with [90Y]microspheres shows that the narrow energy window (460-562 keV) recovers activity concentrations in small spheres better than the regular energy window (425-650 keV) with the beta value of 2000 using the TOF-BSREM algorithm. For the images obtained from patients with [68Ga]citrate using TOF-BSREM with beta = 500, the visual analogue scale (VAS) was improved by 17 % and the Likert score was increased by 1 point on average, both in comparison to corresponding scores for images reconstructed using TOF-OSEM.

CONCLUSION

Our investigation shows that the TOF-BSREM algorithm improves the image quality and quantitative accuracy in low-count PET imaging scenarios. However, the beta value in this algorithm needed to be adjusted for each radiopharmaceutical and counting statistics at the time of scans.

The number of microspheres in Y90 radioembolization directly affects normal tissue radiation exposure

PURPOSE

In Y90 radioembolization, the number of microspheres infused varies by more than a factor of 20 over the shelf-life of the glass radioembolization device. We investigated the effect of the number of Y90 microspheres on normal liver tissue.

METHOD

Healthy pigs received lobar radioembolization with glass Y90 microspheres at 4, 8, 12, and 16 days post-calibration, representing a > 20× range in the number of microspheres deposited per milliliter in tissue. Animals were survived for 1-month post-treatment and the livers were explanted and scanned on a micro CT system to fully characterize the microscopic distribution of individual microspheres. A complete 3D microdosimetric evaluation of each liver was performed with a spatially correlated analysis of histopathologic effect.

RESULTS

Through whole-lobe microscopic identification of each microsphere, a consistent number of microspheres per sphere cluster was found at 4, 8, and 12 days postcalibration, despite an 8-fold increase in total microspheres infused from days 4 to 12. The additional microspheres instead resulted in more clusters formed and, therefore, a more homogeneous microscopic absorbed dose. The increased absorbed-dose homogeneity resulted in a greater volume fraction of the liver receiving a potentially toxic absorbed dose based on radiobiologic models. Histopathologic findings in the animals support a possible increase in normal liver toxicity in later treatments with more spheres (i.e., ≥ day 12) compared to early treatments with less spheres (i.e., ≤ day 8).

CONCLUSION

The microdosimetric evidence presented supports a recommendation of caution when treating large volumes (e.g., right lobe) using glass ⁹⁰Y microspheres at more than 8 days post-calibration, i.e., after "2nd week" Monday. The favorable normal tissue microscopic distribution and associated low toxicity of first week therapies may encourage opportunities for dose escalation with glass microspheres and could also be considered for patients with decreased hepatic reserve.

Comparative dosimetry between 99mTc-MAA SPECT/CT and 90Y PET/CT in primary and metastatic liver tumors

INTRODUCTION

The aim of this study is to determine whether ^99mTc-MAA SPECT/CT-based dosimetry could predict the actual absorbed dose in hepatocellular carcinoma (HCC) or liver metastases, treated by glass or resin microspheres.

MATERIAL AND METHODS

Fifty-seven patients who underwent selective internal radiation therapy (SIRT) were retrospectively included in the study, for a total of 59 treatments. Nineteen HCC were treated by resin microspheres (HCC-SIR), 20 HCC with glass microspheres (HCC-Thera), and 20 liver metastases with resin microspheres (Metastases-SIR). The mean absorbed doses in tumoral liver (Dm) and non-tumoral liver (DmNTL) were determined on the ^99mTc-MAA SPECT/CT and the ⁹⁰Y PET/CT, and compared with each other.

RESULTS

DmNTL was < 50 Gy in the 3 groups, with a strong correlation in all population, albeit slightly lower in Metastases-SIR than HCC-SIR and HCC-Thera (CCC 0.8, 0.94 and 0.96, respectively). In tumoral liver, Dm was higher in HCC than metastases (159 ± 117 Gy versus 63 ± 31 Gy). ^99mTc-MAA SPECT/CT proved to be a better indicator of Dm in HCC compared with metastases, with similar ^99mTc-MAA-⁹⁰Y concordance in resin and glass microspheres (CCC HCC-SIR 0.82, CCC HCC-Thera 0.82, and CCC Metastases-SIR 0.52).

CONCLUSION

^99mTc-MAA SPECT/CT is a reasonably reliable tool for predicting the dose to the non-tumoral liver in both HCC and metastases, regardless of the type of microspheres. It is also fairly reliable for predicting the tumor dose in HCC, again regardless of the type of spheres, although individual variations are observed.

Prediction of Tumor Control in 90Y Radioembolization by Logit Models with PET/CT-Based Dose Metrics

The aim of this work was to develop models for tumor control probability (TCP) in radioembolization with ⁹⁰Y PET/CT-derived radiobiologic dose metrics. Methods: Patients with primary liver cancer or liver metastases who underwent radioembolization with glass microspheres were imaged with ⁹⁰Y PET/CT for voxel-level dosimetry to determine lesion absorbed dose (AD) metrics, biological effective dose (BED) metrics, equivalent uniform dose, and equivalent uniform BED for 28 treatments (89 lesions). The lesion dose-shrinkage correlation was assessed on the basis of RECIST and, when available, modified RECIST (mRECIST) at first follow-up. For a subset with mRECIST, logit regression TCP models were fit via maximum likelihood to relate lesion-level binary response to the dose metrics. As an exploratory analysis, the nontumoral liver dose-toxicity relationship was also evaluated. Results: Lesion dose-shrinkage analysis showed that there were no significant differences between model parameters for primary and metastatic subgroups and that correlation coefficients were superior with mRECIST. Therefore, subsequent TCP analysis was performed for the combined group using mRECIST only. The overall lesion-level mRECIST response rate was 57%. The AD and BED metrics yielding 50% TCP were 292 and 441 Gy, respectively. All dose metrics considered for TCP modeling, including mean AD, were significantly associated with the probability of response, with high areas under the curve (0.87-0.90, P < 0.0001) and high sensitivity (>0.75) and specificity (>0.83) calculated using a threshold corresponding to 50% TCP. Because nonuniform AD deposition by microspheres cannot be determined by PET at a microscopic scale, radiosensitivity values extracted here by fitting models to clinical response data were substantially lower than reported for in vitro cell cultures or for external-beam radiotherapy clinical studies. There was no correlation between nontumoral liver AD and toxicity measures. Conclusion: Despite the heterogeneous patient cohort, logistic regression TCP models showed a strong association between various dose metrics and the probability of response. The performance of mean AD was comparable to that of radiobiologic dose metrics that involve more complex calculations. These results demonstrate the importance of considering TCP in treatment planning for radioembolization.

Validation of post-treatment PET-based dosimetry software for hepatic radioembolization of Yttrium-90 microspheres

PURPOSE

Yttrium-90 (⁹⁰ Y) microsphere radioembolization enables selective internal radiotherapy for hepatic malignancies. Currently, there is no standard postdelivery imaging and dosimetry of the microsphere distribution to verify treatment. Recent studies have reported utilizing the small positron yield of ⁹⁰ Y (32 ppm) with positron emission tomography (PET) to perform treatment verification and dosimetry analysis. In this study, we validated a commercial dosimetry software, MIM SurePlan™ LiverY90 (MIM Software Inc., Cleveland, OH), for clinical use.

METHODS

A MATLAB-based algorithm for ⁹⁰ Y PET-based dosimetry was developed in-house and validated for the purpose of commissioning the commercial software. The algorithm is based on voxel S values and dosimetry formalism reported in MIRD Pamphlet 17. We validated the in-house algorithm to establish it as the ground truth by comparing results from a digital point phantom and a digital uniform cylinder to manual calculations. Once we validated our in-house MATLAB-based algorithm, we used it to perform acceptance testing and commissioning of the commercial dosimetry software, MIM SurePlan, which uses the same dosimetry formalism. A 0.4 cm/5% gamma test was performed on PET-derived dose maps from each algorithm of uniform digital and nonuniform physical phantoms filled with ⁹⁰ Y chloride solution. Average dose (D_avg ) and minimum dose to 70% (D₇₀ ) of a given volume of interest (VOI) were compared for the digital phantom, the physical phantom, and five patient cases (27 tumor VOIs), representing different clinical scenarios.

RESULTS

The gamma-pass rates were 97.26% and 97.66% for the digital and physical phantoms, respectively. The differences between D_avg and D₇₀ were 0.076% and 0.10% for the digital phantom, respectively, and <5.2% for various VOIs in the physical phantom. In the clinical cases, 96.3% of the VOIs had a difference <5% for D_avg , and 88.9% of the VOIs had a difference <5% for D₇₀ .

CONCLUSIONS

Dose calculation results from MIM SurePlan were found to be in good agreement with our in-house algorithm. This indicates that MIM SurePlan performs as it should and, hence, can be deemed accepted and commissioned for clinical use for post-implant PET-based dosimetry of ⁹⁰ Y radioembolization.

Dose volume histogram-based optimization of image reconstruction parameters for quantitative 90 Y-PET imaging

PURPOSE

⁹⁰ Y-microsphere radioembolization or selective internal radiation therapy is increasingly being used as a treatment option for tumors that are not candidates for surgery and external beam radiation therapy. Recently, volumetric ⁹⁰ Y-dosimetry techniques have been implemented to explore tumor dose-response on the basis of 3D ⁹⁰ Y-activity distribution from PET imaging. Despite being a theranostic study, the optimization of quantitative ⁹⁰ Y-PET image reconstruction still uses the mean activity concentration recovery coefficient (RC) as the objective function, which is more relevant to diagnostic and detection tasks than is to dosimetry. The aim of this study was to optimize ⁹⁰ Y-PET image reconstruction by minimizing errors in volumetric dosimetry via the dose volume histogram (DVH). We propose a joint optimization of the number of equivalent iterations (the product of the iterations and subsets) and the postreconstruction filtration (FWHM) to improve the accuracy of voxel-level ⁹⁰ Y dosimetry.

METHODS

A modified NEMA IEC phantom was used to emulate clinically relevant ⁹⁰ Y-PET imaging conditions through various combinations of acquisition durations, activity concentrations, sphere-to-background ratios, and sphere diameters. PET data were acquired in list mode for 300 min in a single-bed position; we then rebinned the list mode PET data to 60, 45, 30, 15, and 5 min per bed, with 10 different realizations. Errors in the DVH were calculated as root mean square errors (RMSE) of the differences in the image-based DVH and the expected DVH. The new optimization approach was tested in a phantom study, and the results were compared with the more commonly used objective function of the mean activity concentration RC.

RESULTS

In a wide range of clinically relevant imaging conditions, using 36 equivalent iterations with a 5.2-mm filtration resulted in decreased systematic errors in volumetric ⁹⁰ Y dosimetry, quantified as image-based DVH, in ⁹⁰ Y-PET images reconstructed using the ordered subset expectation maximization (OSEM) iterative reconstruction algorithm with time of flight (TOF) and point spread function (PSF) modeling. Our proposed objective function of minimizing errors in DVH, which allows for joint optimization of ⁹⁰ Y-PET iterations and filtration for volumetric quantification of the ⁹⁰ Y dose, was shown to be superior to conventional RC-based optimization approaches for image-based absorbed dose quantification.

CONCLUSION

Our proposed objective function of minimizing errors in DVH, which allows for joint optimization of iterations and filtration to reduce errors in the PET-based volumetric quantification ⁹⁰ Y dose, is relevant to dosimetry in therapy procedures. The proposed optimization method using DVH as the objective function could be applied to any imaging modality used to assess voxel-level quantitative information.

The physics of radioembolization

Radioembolization is an established treatment for chemoresistant and unresectable liver cancers. Currently, treatment planning is often based on semi-empirical methods, which yield acceptable toxicity profiles and have enabled the large-scale application in a palliative setting. However, recently, five large randomized controlled trials using resin microspheres failed to demonstrate a significant improvement in either progression-free survival or overall survival in both hepatocellular carcinoma and metastatic colorectal cancer. One reason for this might be that the activity prescription methods used in these studies are suboptimal for many patients.In this review, the current dosimetric methods and their caveats are evaluated. Furthermore, the current state-of-the-art of image-guided dosimetry and advanced radiobiological modeling is reviewed from a physics' perspective. The current literature is explored for the observation of robust dose-response relationships followed by an overview of recent advancements in quantitative image reconstruction in relation to image-guided dosimetry.This review is concluded with a discussion on areas where further research is necessary in order to arrive at a personalized treatment method that provides optimal tumor control and is clinically feasible.

Comparison of PET/CT and PET/MR imaging and dosimetry of yttrium-90 (90Y) in patients with unresectable hepatic tumors who have received intra-arterial radioembolization therapy with 90Y microspheres

BACKGROUND

The aim of our study was to compare 90Y dosimetry obtained from PET/MRI versus PET/CT post-therapy imaging among patients with primary or metastatic hepatic tumors. First, a water-filled Jaszczak phantom containing fillable sphere with 90Y-chloride was acquired on both the PET/CT and PET/MRI systems, in order to check the cross-calibration of the modalities. Following selective internal radiation therapy (SIRT) with 90Y microspheres, 32 patients were imaged on a PET/CT system, immediately followed by a PET/MRI study. Reconstructed images were transferred to a common platform and used to calculate 90Y dosimetry. A Passing-Bablok regression scatter diagram and the Bland and Altman method were used to analyze the difference between the dosimetry values.

RESULTS

The phantom study showed that both modalities were calibrated with less than 1% error. The mean liver doses for the 32 subjects calculated from PET/CT and PET/MRI were 51.6 ± 24.7 Gy and 46.5 ± 22.7 Gy, respectively, with a mean difference of 5.1 ± 5.0 Gy. The repeatability coefficient was 9.0 (18.5% of the mean). The Spearman rank correlation coefficient was very high, ρ = 0.97. Although the maximum dose to the liver can be significantly different (up to 40%), mean liver doses from each modalities were relatively close, with a difference of 18.5% or less.

CONCLUSIONS

The two main contributors to the difference in 90Y dosimetry calculations using PET/CT versus PET/MRI can be attributed to the differences in regions of interest (ROIs) and differences attributed to attenuation correction. Due to the superior soft-tissue contrast of MRI, liver contours are usually better seen than in CT images. However, PET/CT provides better quantification of PET images, due to better attenuation correction. In spite of these differences, our results demonstrate that the dosimetry values obtained from PET/MRI and PET/CT in post-therapy 90Y studies were similar.

Feasibility assessment of yttrium-90 liver radioembolization imaging using amplitude-based gated PET/CT

Purpose

The usage of PET/computed tomography (CT) to monitor hepatocellular carcinoma patients following yttrium-90 (⁹⁰Y) radioembolization has increased. Respiratory motion causes liver movement, which can be corrected using gating techniques at the expense of added noise. This work examines the use of amplitude-based gating on ⁹⁰Y-PET/CT and its potential impact on diagnostic integrity.

Patients and methods

Patients were imaged using PET/CT following ⁹⁰Y radioembolization. A respiratory band was used to collect respiratory cycle data. Patient data were processed as both standard and motion-corrected images. Regions of interest were drawn and compared using three methods. Activity concentrations were calculated and converted into dose estimates using previously determined and published scaling factors. Diagnostic assessments were performed using a binary scale created from published ⁹⁰Y-PET/CT image interpretation guidelines.

Results

Estimates of radiation dose were increased (P<0.05) when using amplitude-gating methods with ⁹⁰Y PET/CT imaging. Motion-corrected images show increased noise, but the diagnostic determination of success, using the Kao criteria, did not change between static and motion-corrected data.

Conclusion

Amplitude-gated PET/CT following ⁹⁰Y radioembolization is feasible and may improve ⁹⁰Y dose estimates while maintaining diagnostic assessment integrity.

SPECT/CT image-based dosimetry for Yttrium-90 radionuclide therapy: Application to treatment response

This work demonstrates the efficacy of voxel‐based ⁹⁰Y microsphere dosimetry utilizing post‐therapy SPECT/CT imaging and applies it to the prediction of treatment response for the management of patients with hepatocellular carcinoma (HCC). A ⁹⁰Y microsphere dosimetry navigator (RapidSphere) within a commercial platform (Velocity, Varian Medical Systems) was demonstrated for three microsphere cases that were imaged using optimized bremsstrahlung SPECT/CT. For each case, the ⁹⁰Y SPECT/CT was registered to follow‐up diagnostic MR/CT using deformable image registration. The voxel‐based dose distribution was computed using the local deposition method with known injected activity. The system allowed the visualization of the isodose distributions on any of the registered image datasets and the calculation of dose‐volume histograms (DVHs). The dosimetric analysis illustrated high local doses that are characteristic of blood‐flow directed brachytherapy. In the first case, the HCC mass demonstrated a complete response to treatment indicated by a necrotic region in follow‐up MR imaging. This result was dosimetrically predicted since the gross tumor volume (GTV) was well covered by the prescription isodose volume (V150 Gy = 85%). The second case illustrated a partial response to treatment which was characterized by incomplete necrosis of an HCC mass and a remaining area of solid enhancement in follow‐up MR imaging. This result was predicted by dosimetric analysis because the GTV demonstrated incomplete coverage by the prescription isodose volume (V470 Gy = 18%). The third case demonstrated extrahepatic activity. The dosimetry indicated that the prescription (125 Gy) isodose region extended outside of the liver into the duodenum (178 Gy maximum dose). This was predictive of toxicity as the patient later developed a duodenal ulcer. The ability to predict outcomes and complications using deformable image registration, calculated isodose distributions, and DVHs, points to the clinical utility of patient‐specific dose calculations for ⁹⁰Y radioembolization treatment planning.

Hepatic selective internal radiation therapy: how well does pretreatment 99mTc-macroaggregated albumin activity distribution and quantification agree with post-therapy bremsstrahlung imaging?

AIM

The aim of this study was to examine the agreement of pretreatment Tc-macroaggregated albumin imaging performed for selective internal radiation therapy (SIRT) workup with Y percentage lung shunt (PLS) and regional hepatic distribution in subsequent post-therapy bremsstrahlung imaging.

PATIENTS AND METHODS

Planar images were used to calculate PLS. The significant Y bremsstrahlung scatter required background correction. Results using both Y lung background regions of interest (ROI) reported in previous studies and extended ROIs (reflecting lung background variation) were compared with Tc-MAA PLS. Lesion and healthy liver volumes were outlined on diagnostic computed tomography scans and registered to Tc-MAA and Y single-photon emission computed tomography/computed tomographydata. Single-photon emission computed tomography voxel values were normalized to injected Y activity. Volume mean activities were calculated, and converted into the mean absorbed dose. Agreement was quantified using Bland-Altman analysis.

RESULTS

PLS: The bias using previous studies' lung background ROIs was -10.71%, with a 95% confidence interval (CI) of -18.79 to -2.64%. The extended ROI yielded a bias of 0.77% (95% CI: -2.23 to 3.70%). Liver: The healthy liver bias was 0.01 MBq/ml (0.17 Gy), with a -0.05 to 0.06 MBq/ml (95% CI:0.80 -1.93 Gy). The lesion mean activity/ml bias was -0.02 MBq/ml (3.71 Gy), with a -0.81 to 0.76 MBq/ml (95% CI: -35.49 to 28.07 Gy).

CONCLUSIONS

The PLS agreement was sensitive to the Y lung background correction ROI, potentially explaining a previously published controversy. The mean activity and absorbed dose agreement for the metastatic lesions was poorer than the healthy liver volumes studied here.

A Woman with Black Beads in Her Stomach: Severe Gastric Ulceration Caused by Yttrium-90 Radioembolization

Radioembolization (RE) is a selective internal radiation therapy (SIRT) delivering targeted, high-dose, intra-arterial radiation directly to the vascular supply of liver tumors. Complications can occur due to aberrant deposition or migration of radiation microspheres into nontarget locations, including normal hepatic parenchyma, lungs, pancreas, and upper gastrointestinal (UGI) tract. We report a case of gastric ulcers due to yttrium-90 (⁹⁰Y) seed migration to the stomach to alert clinicians to this rare cause of gastric injury. A 57-year-old woman with stage IV breast cancer with liver and lung metastases presented to the hospital with 2 months of worsening nausea and vomiting. Two months prior, she had received SIRT with ⁹⁰Y microspheres without complications. Upper GI endoscopy showed diffuse gastritis and extensive antral ulceration. Biopsies revealed black, spherical foreign bodies, consistent with ⁹⁰Y microspheres, documenting radiation injury. Radiation-induced UGI ulceration is caused by direct radiation injury from beta-radiation. Delay in diagnosis may be due to the nonspecificity of symptoms and temporal delay of symptom onset from SIRT, which was 2 months in our patient. Also, complaints may be attributed erroneously to adjuvant chemotherapy or widespread metastatic disease. Clinicians must consider radiation-associated toxicity in any SIRT-treated patient developing abdominal symptoms.