Comparing voxel-based absorbed dosimetry methods in tumors, liver, lung, and at the liver-lung interface for (90)Y microsphere selective internal radiation therapy

Averaging of absorbed doses: How matter matters

BACKGROUND

Dosimetry in radionuclide therapy often requires the calculation of average absorbed doses within and between spatial regions, for example, for voxel-based dosimetry methods, for paired organs, or across multiple tumors. Formation of such averages can be made in different ways, starting from different definitions.

PURPOSE

The aim of this study is to formally specify different averaging strategies for absorbed doses, and to compare their results when applied to absorbed dose distributions that are non-uniform within and between regions.

METHODS

For averaging within regions, two definitions of the average absorbed dose are considered: the simple average over the region (the region average) and the average when weighting by the mass density (density-weighted region average). The latter is shown to follow from the definition of mean absorbed dose according to the ICRU, and to be consistent with the MIRD formalism. For averaging between different spatial regions, three definitions follow: the volume-weighted, the mass-weighted, and the unweighted average. With respect to characterizing non-uniformity, the different average definitions lead to the use of dose-volume histograms (DVHs) (region average), dose-mass histograms (DMHs) (density-weighted region average), and unweighted histograms (unweighted average). Average absorbed doses are calculated for three worked examples, starting from the different definitions. The first, schematic, example concerns the calculation of the average absorbed dose between two regions with different volumes or mass densities. The second, stylized, example concerns voxel-based dosimetry, for which the average absorbed-dose rate within a region is calculated. The geometries studied include three 177 Lu-filled voxelized spheres, where the sphere masses are held constant while the material compositions, densities, and volumes are varied. For comparison, the mean absorbed-dose rates obtained using unit-density sphere S-values are also included. The third example concerns SPECT/CT-based tumor dosimetry for five patients undergoing therapy with 177 Lu-PSMA and six patients undergoing therapy with 177 Lu-DOTA-TATE, for which the average absorbed-dose rates across multiple tumors are calculated. For the second and third examples, analyses also include representations by histograms.

RESULTS

Example 1 shows that the average absorbed doses, calculated using different definitions, can differ considerably if the masses and absorbed doses for two regions are markedly different. From example 2 it is seen that the density-weighted region average is stable under different activity and density distributions and is also in line with results using S-values. In contrast, the region average varies as function of the activity distribution. In example 3, the absorbed dose rates for individual tumors differ by (1.1 ± 4.3)% and (-0.1 ± 0.4)% with maximum deviations of +34.4% and -1.4% for 177 Lu-PSMA and 177 Lu-DOTA-TATE, respectively, when calculated as region averages or density-weighted region averages, with largest deviations obtained when the density is non-uniform. The average absorbed doses calculated across all tumors are similar when comparing mass-weighted and volume-weighted averages but these differ substantially from unweighted averages.

CONCLUSION

Different strategies for averaging of absorbed doses within and between regions can lead to substantially different absorbed-dose estimates. At reporting of radionuclide therapy dosimetry, it is important to specify the averaging strategy applied.

Influence on voxel-based dosimetry: noise effect on absorbed dose dosimetry at single time-point versus sequential single-photon emission computed tomography

OBJECTIVE

The purpose of this study was to evaluate how statistical fluctuation in single-photon emission computed tomography (SPECT) images propagate to absorbed dose maps.

METHODS

SPECT/computed tomography (CT) images of iodine-131 filled phantoms, using different acquisition and processing protocols, were evaluated using STRATOS software to assess the absorbed dose distribution at the voxel level. Absorbed dose values and coefficient of variation (COV) were analyzed for dosimetry based on single time-point SPECT images and time-integrated activities of SPECT sequences with low and high counts.

RESULTS

Considering dosimetry based on a single time-point, the mean absorbed dose was not significantly affected by total counts or reconstruction parameters, but the uniformity of the absorbed dose maps had an almost linear correlation with SPECT noise. When high- and low-count SPECT sequences were used to generate an absorbed dose map, the absorbed dose COV for each of the temporal sequences was slightly lower than the absorbed dose COV based on the single SPECT image with the highest count included in the sequence.

CONCLUSION

The impact of changes in SPECT counts and reconstruction parameters is almost linear when dosimetry is based on isolated SPECT images, but less pronounced when dosimetry is based on sequential SPECTs.

Single-Compartment Dose Prescriptions for Ablative 90Y-Radioembolization Segmentectomy

BACKGROUND

Yttrium-90 (⁹⁰Y) radioembolization is increasingly being utilized with curative intent. While single-compartment doses with respect to the perfused volume for the complete pathologic necrosis (CPN) of tumors have been reported, the actual doses delivered to the tumor and at-risk margins that leads to CPN have hitherto not been estimated. We present an ablative dosimetry model that calculates the dose distribution for tumors and at-risk margins based on numerical mm-scale dose modeling and the available clinical CPN evidence and report on the necessary dose metrics needed to achieve CPN following ⁹⁰Y-radioembolization.

METHODS

Three-dimensional (3D) activity distributions (MBq/voxel) simulating spherical tumors were modeled with a 121 × 121 × 121 mm³ soft tissue volume (1 mm³ voxels). Then, 3D dose distributions (Gy/voxel) were estimated by convolving 3D activity distributions with a ⁹⁰Y 3D dose kernel (Gy/MBq) sized 61 × 61 × 61 mm³ (1 mm³ voxels). Based on the published data on single-compartment segmental doses for the resected liver samples of HCC tumors showing CPN after radiation segmentectomy, the nominal voxel-based mean tumor dose (DmeanCPN), point dose at tumor rim (DrimCPN), and point dose 2 mm beyond the tumor boundary (D2mmCPN), which are necessary to achieve CPN, were calculated. The single-compartment dose prescriptions to required achieve CPN were then analytically modeled for more general cases of tumors with diameters dt = 2, 3, 4, 5, 6, and 7 cm and with tumor-to-normal-liver uptake ratios T:N = 1:1, 2:1, 3:1, 4:1, and 5:1.

RESULTS

The nominal case defined to estimate the doses needed for CPN, based on the previously published clinical data, was a single hyperperfused tumor with a diameter of 2.5 cm and T:N = 3:1, treated with a single-compartment segmental dose of 400 Gy. The voxel-level doses necessary to achieve CPN were 1053 Gy for the mean tumor dose, 860 Gy for the point dose at the tumor boundary, and 561 Gy for the point dose at 2 mm beyond the tumor edge. The single-compartment segmental doses necessary to satisfy the criteria for CPN in terms of the mean tumor dose, point dose at the tumor boundary, and the point dose at 2 mm beyond the tumor edge were tabulated for a range of tumor diameters and tumor-to-normal-liver uptake ratios.

CONCLUSIONS

The analytical functions that describe the relevant dose metrics for CPN and, more importantly, the single-compartment dose prescriptions for the perfused volume needed to achieve CPN are reported for a large range of conditions in terms of tumor diameters (1-7 cm) and T:N uptake ratios (2:1-5:1).

Voxel-S-Value based 3D treatment planning methods for Y-90 microspheres radioembolization based on Tc-99m-macroaggregated albumin SPECT/CT

Partition model (PM) for Y-90 microsphere radioembolization is limited in providing 3D dosimetrics. Voxel-S-Values (VSV) method has good agreement with Monte Carlo (MC) simulations for 3D absorbed dose conversion. We propose a new VSV method and compare its performance along with PM, MC and other VSV methods for Y-90 RE treatment planning based on Tc-99m MAA SPECT/CT. Twenty Tc-99m-MAA SPECT/CT patient data are retrospectively analyzed. Seven VSV methods are implemented: (1) local energy deposition; (2) liver kernel; (3) liver kernel and lung kernel; (4) liver kernel with density correction (LiKD); (5) liver kernel with center voxel scaling (LiCK); (6) liver kernel and lung kernel with density correction (LiLuKD); (7) proposed liver kernel with center voxel scaling and lung kernel with density correction (LiCKLuKD). Mean absorbed dose and maximum injected activity (MIA) obtained by PM and VSV are evaluated against MC results, and 3D dosimetrics generated by VSV are compared with MC. LiKD, LiCK, LiLuKD and LiCKLuKD have the smallest deviation in normal liver and tumors. LiLuKD and LiCKLuKD have the best performance in lungs. MIAs are similar by all methods. LiCKLuKD could provide MIA consistent with PM, and precise 3D dosimetrics for Y-90 RE treatment planning.

Voxel-S-value methods adapted to heterogeneous media for quantitative Y-90 microsphere radioembolization dosimetry

PURPOSE

The absorbed dose estimation from Voxel-S-Value (VSV) method in heterogeneous media is suboptimal as VSVs are calculated in homogeneous media. The aim of this study is to develop and evaluate new VSV methods in order to enhance the accuracy of Y-90 microspheres absorbed dose estimation in liver, lungs, tumors and lung-liver interface regions.

METHODS

Ten patients with Y-90 microspheres SPECT/CT and PET/CT data, six of whom had additional Tc-99m-macroaggregated albumin SPECT/CT data, were analyzed from the Deep Blue Data Repository. Seven existing VSV methods along with three newly proposed VSV methods were evaluated: liver and lung kernel with center voxel scaling (LiLuCK), liver kernel with density correction and lung kernel with center voxel scaling (LiKDLuCK), liver kernel with center voxel scaling and lung kernel with density correction (LiCKLuKD). Monte Carlo (MC) results were regarded as the gold standard. Absolute absorbed dose errors (%AADE) of these methods for the liver, lungs, tumors, upper liver, and lower lungs were assessed.

RESULTS

Liver and tumor's median %AADE of all methods were <3% for three types of imaging data. In the lungs, however, three recently proposed VSV methods provided median %AADEs of less than 7%, whereas the differences exceeded 20% for existing methods that did not use a lung kernel. LiCKLuKD could achieve median %AADE <2% in the liver, upper liver and tumors, and median %AADE <7% in the lungs and lower lungs in three types of data.

CONCLUSION

All methods are consistent with MC in the liver and tumors. Methods with tissue-specific kernel and effective correction achieve smaller errors in lungs. LiCKLuKD has comparable results with MC in absorbed dose estimation of Y-90 radioembolization for all target regions.

The American Brachytherapy Society consensus statement for permanent implant brachytherapy using Yttrium-90 microsphere radioembolization for liver tumors

PURPOSE

To develop a multidisciplinary consensus for high quality multidisciplinary implementation of brachytherapy using Yttrium-90 (90Y) microspheres transarterial radioembolization (90Y TARE) for primary and metastatic cancers in the liver.

METHODS AND MATERIALS

Members of the American Brachytherapy Society (ABS) and colleagues with multidisciplinary expertise in liver tumor therapy formulated guidelines for 90Y TARE for unresectable primary liver malignancies and unresectable metastatic cancer to the liver. The consensus is provided on the most recent literature and clinical experience.

RESULTS

The ABS strongly recommends the use of 90Y microsphere brachytherapy for the definitive/palliative treatment of unresectable liver cancer when recommended by the multidisciplinary team. A quality management program must be implemented at the start of 90Y TARE program development and follow-up data should be tracked for efficacy and toxicity. Patient-specific dosimetry optimized for treatment intent is recommended when conducting 90Y TARE. Implementation in patients on systemic therapy should account for factors that may enhance treatment related toxicity without delaying treatment inappropriately. Further management and salvage therapy options including retreatment with 90Y TARE should be carefully considered.

CONCLUSIONS

ABS consensus for implementing a safe 90Y TARE program for liver cancer in the multidisciplinary setting is presented. It builds on previous guidelines to include recommendations for appropriate implementation based on current literature and practices in experienced centers. Practitioners and cooperative groups are encouraged to use this document as a guide to formulate their clinical practices and to adopt the most recent dose reporting policies that are critical for a unified outcome analysis of future effectiveness studies.

Relevance of artefacts in 99mTc-MAA SPECT scans on pre-therapy patient-specific 90Y TARE internal dosimetry: a GATE Monte Carlo study

Objective. The direct Monte Carlo (MC) simulation of radiation transport exploiting morphological and functional tomographic imaging as input data is considered the gold standard for internal dosimetry in nuclear medicine, and it is increasingly used in studies regarding trans-arterial radio-embolization (TARE). However, artefacts affecting the functional scans, such as reconstruction artefacts and motion blurring, decrease the accuracy in defining the radionuclide distribution in the simulations and consequently lead to errors in absorbed dose estimations. In this study, the relevance of such artefacts in patient-specific three-dimensional MC dosimetry was investigated in three cases of 90Y TARE. Approach. The pre-therapy 99mTc MacroAggregate Albumin (Tc-MAA) SPECTs and CTs of patients were used as input for simulations performed with the GEANT4-based toolkit GATE. Several pre-simulation SPECT-masking techniques were implemented, with the aim of zeroing the decay probability in air, in lungs, or in the whole volume outside the liver. Main results. Increments in absorbed dose up to about +40% with respect to the native-SPECT simulations were found in liver-related volumes of interest (VOIs), depending on the masking procedure adopted. Regarding lungs-related VOIs, decrements in absorbed doses in right lung as high as −90% were retrieved. Significance. These results highlight the relevant influence of SPECT artefacts, if not properly treated, on dosimetric outcomes for 90Y TARE cases. Well-designed SPECT-masking techniques appear to be a promising way to correct for such misestimations.

Organ-level internal dosimetry for intra-hepatic-arterial administration of 99m Tc-macroaggregated albumin

PURPOSE

There are no published data on organ doses following intra-hepatic-arterial administration of ^99m Tc-macroaggregated-albumin (IHA ^99m Tc-MAA) routinely used in ⁹⁰ Y-radioembolization-treatment planning to assess intra- and extra-hepatic depositions and calculate lung-shunt-fraction (LSF). We propose a method to model the organ doses following IHA ^99m Tc-MAA that incorporates three in vivo constituent biodistributions, the ^99m Tc-MAA that escape the liver due to LSF, and the ^99m Tc-MAA disassociation fraction (DF).

METHODS

The potential in vivo biodistributions for IHA ^99m Tc-MAA are: Liver-Only MAA with all activity sequestered in the liver (LSF = 0&DF = 0), Intravenous MAA with all activity transferred intravenously as ^99m Tc-MAA (LSF = 1&DF = 0), and Intravenous Pertechnetate with all activity is transferred intravenously as ^99m Tc-pertechnetate (LSF = 0&DF = 1). Organ doses for Liver-Only MAA were determined using OLINDA/EXM 2.2, where liver was modeled as the source organ containing ^99m Tc-MAA, while those for Intravenous MAA and Intravenous Pertechnetate were from ICRP 128. Organ doses for the general case can be determined as a weighted-linear-combination of the three constituent biodistributions depending on the LSF and DF. The maximum-dose scenario was modeled by selecting the highest dose rate for each organ amongst the three constituent cases.

RESULTS

For Liver-Only MAA, the liver as source organ received the highest dose at 98.6 and 126 mGy/GBq for the Adult Male and Adult Female phantoms, respectively; all remaining organs received <27 and <32 mGy/GBq. For Intravenous MAA, the lung as source organ received the highest dose at 66 and 97 mGy/GBq; all remaining organs received <16 and <21 mGy/GBq. The organ with the highest dose for Intravenous Pertechnetate was the upper-large-intestinal wall at 56 and 73 mGy/GBq; all remaining organs received <26 and <34 mGy/GBq. The liver and lung doses for the maximum-dose scenario with 5 mCi (185 MBq) ^99m Tc-MAA were estimated at 18.2 and 12.2 mGy, and 23.3 and 17.9 mGy, for the Adult Male and Adult Female phantoms, respectively.

CONCLUSION

Organ dose estimates following IHA ^99m Tc-MAA based on constituent biodistribution models and patient-specific LSF and DF values have been derived. Liver and lung were the organs with highest dose, receiving at most 15 - 25 mGy in the maximum-dose scenario, following 5 mCi IHA ^99m Tc-MAA. This article is protected by copyright. All rights reserved.

Pre- and post-treatment image-based dosimetry in90Y-microsphere radioembolization using the TOPAS Monte Carlo toolkit

Objective. To evaluate the pre-treatment and post-treatment imaging-based dosimetry of patients treated with 90Y-microspheres, including accurate estimations of dose to tumor, healthy liver and lung. To do so, the Monte Carlo (MC) TOPAS platform is in this work extended towards its utilization in radionuclide therapy.Approach. Five patients treated at the Massachusetts General Hospital were selected for this study. All patients had data for both pre-treatment SPECT-CT imaging using 99mTc-MAA as a surrogate of the 90Y-microspheres treatment and SPECT-CT imaging immediately after the 90Y activity administration. Pre- and post-treatment doses were computed with TOPAS using the SPECT images to localize the source positions and the CT images to account for tissue inhomoegeneities. We compared our results with analytical calculations following the voxel-based MIRD scheme.Main results. TOPAS results largely agreed with the MIRD-based calculations in soft tissue regions: the average difference in mean dose to the liver was 0.14 Gy GBq^-1(2.6%). However, dose distributions in the lung differed considerably: absolute differences in mean doses to the lung ranged from 1.2 to 6.3 Gy GBq^-1and relative differences from 153% to 231%. We also found large differences in the intra-hepatic dose distributions between pre- and post-treatment imaging, but only limited differences in the pulmonary dose.Significance. Doses to lung were found to be higher using TOPAS with respect to analytical calculations which may significantly underestimate dose to the lung, suggesting the use of MC methods for 90Y dosimetry. According to our results, pre-treatment imaging may still be representative of dose to lung in these treatments.

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.

Precision Interventional Brachytherapy: A Promising Strategy Toward Treatment of Malignant Tumors

Precision interventional brachytherapy is a radiotherapy technique that combines radiation therapy medicine with computer network technology, physics, etc. It can solve the limitations of conventional brachytherapy. Radioactive drugs and their carriers change with each passing day, and major research institutions and enterprises worldwide have conducted extensive research on them. In addition, the capabilities of interventional robotic systems are also rapidly developing to meet clinical needs for the precise delivery of radiopharmaceuticals in interventional radiotherapy. This study reviews the main radiopharmaceuticals, drug carriers, dispensing and fixation technologies, and interventional robotic precision delivery systems used in precision brachytherapy of malignant tumors. We then discuss the current needs in the field and future development prospects in high-precision interventional brachytherapy.

A Prospective Phase II Study of Safety and Efficacy of Sorafenib Followed by 90Y Glass Microspheres for Patients with Advanced or Metastatic Hepatocellular Carcinoma

Purpose

The most common cause of death in advanced/metastatic hepatocellular carcinoma (HCC) is liver failure due to tumor progression. While retrospective studies and meta-analyses of systemic therapy combined with liver-directed therapy have been performed, prospective studies of safety/efficacy of antiangiogenesis followed by intra-arterial therapies are lacking. We tested our hypothesis that sorafenib followed by yttrium 90 glass microspheres (⁹⁰Y GMs) is safe and that survival outcomes may improve by controlling hepatic tumors.

Methods

We enrolled 38 Child–Pugh A patients with advanced/metastatic HCC. In sum, 34 received sorafenib, followed after 4 weeks by ⁹⁰Y GMs. Analysis of safety and survival outcomes was performed to assess adverse events, median progression-free survival, and overall survival.

Results

A total of 34 patients were evaluable: 14 (41.2%) with chronic hepatitis, nine (26.5%) with vascular invasion, and eleven (32.4%) with extrahepatic diseases. Safety analysis revealed that the combination therapy was well tolerated. Grade III–IV adverse events comprised fatigue (n=3), diarrhea (n=2), nausea (n=1), vomiting (n=2), hypertension (n=4), thrombocytopenia (n=1), hyperbilirubinemia (n=1), proteinuria (n=1), hyponatremia (n=1), and elevated alanine or aspartate aminotransferase (n=5). Median progression-free and overall survival were 10.4 months (95% CI 5.8–14.4) and 13.2 months (95% CI 7.9–18.9), respectively. Twelve patients (35.3%) achieved partial responses and 16 (47.0%) stable disease. Median duration of sorafenib was 20 (3–90) weeks, and average dose was 622 (466–800) mg daily. Dosimetry showed similar mean doses between planned and delivered calculations to normal liver and tumor:normal liver uptake ratio, with no significant correlation with adverse events at 3 and 6 months post-⁹⁰Y treatment.

Conclusion

This is the first prospective study to evaluate sorafenib followed by ⁹⁰Y in patients with advanced HCC. The study validated our hypothesis of safety with encouraging efficacy signals of the sequencing treatment, and provides proof of concept for future combination modalities for patients with advanced or metastatic HCC.

Clinical Trial Registration Number

NCT01900002.

Post-administration dosimetry in yttrium-90 radioembolization through micro-CT imaging of radiopaque microspheres in a porcine renal model

The purpose of this study is to perform post-administration dosimetry in yttrium-90 radioembolization through micro-CT imaging of radiopaque microsphere distributions in a porcine renal model and explore the impact of spatial resolution of an imaging system on the extraction of specific dose metrics. Following the administration of radiopaque microspheres to the kidney of a hybrid farm pig, the kidney was explanted and imaged with micro-CT. To produce an activity distribution, 400 MBq of yttrium-90 activity was distributed throughout segmented voxels of the embolized vasculature based on an established linear relationship between microsphere concentration and CT voxel value. This distribution was down-sampled to coarser isotropic grids ranging in voxel size from 2.5 to 15 mm to emulate nominal resolutions comparable to those found in yttrium-90 PET and Bremsstrahlung SPECT imaging. Dose distributions were calculated through the convolution of activity distributions with dose-voxel kernels generated using the GATE Monte Carlo toolkit. Contours were computed to represent normal tissue and target volumes. Dose-volume histograms, dose metrics, and dose profiles were compared to a ground truth dose distribution computed with GATE. The mean dose to the target for all studied voxel sizes was found to be within 5.7% of the ground truth mean dose.D70was shown to be strongly correlated with image voxel size of the dose distribution (r² = 0.90).D70is cited in the literature as an important dose metric and its dependence on voxel size suggests higher resolution dose distributions may provide new perspectives on dose-response relationships in yttrium-90 radioembolization. This study demonstrates that dose distributions with large voxels incorrectly homogenize the dose by attributing escalated doses to normal tissues and reduced doses in high-dose target regions. High-resolution micro-CT imaging of radiopaque microsphere distributions can provide increased confidence in characterizing the absorbed dose heterogeneity in yttrium-90 radioembolization.

Clinical Evaluation of a Three-Dimensional Internal Dosimetry Technique for Liver Radioembolization with 90Y Microspheres Using Dose Voxel Kernels

Background: The purpose of this study was to develop a rapid, reliable, and efficient tool for three-dimensional (3D) dosimetry treatment planning and post-treatment evaluation of liver radioembolization with ⁹⁰Y microspheres, using tissue-specific dose voxel kernels (DVKs) that can be used in everyday clinical practice. Materials and Methods: Two tissue-specific DVKs for ⁹⁰Y were calculated through Monte Carlo (MC) simulations. DVKs for the liver and lungs were generated, and the dose distribution was compared with direct MC simulations. A method was developed to produce a 3D dose map by convolving the calculated DVKs with the activity biodistribution derived from clinical single-photon emission computed tomography (SPECT) or positron emission tomography (PET) images. Image registration for the SPECT or PET images with the corresponding computed tomography scans was performed before dosimetry calculation. The authors first compared the DVK convolution dosimetry with a direct full MC simulation on an XCAT anthropomorphic phantom. They then tested it in 25 individual clinical cases of patients who underwent ⁹⁰Y therapy. All MC simulations were carried out using the GATE MC toolkit. Results: Comparison of the measured absorbed dose using tissue-specific DVKs and direct MC simulation on 25 patients revealed a mean difference of 1.07% ± 1.43% for the liver and 1.03% ± 1.21% for the tumor tissue, respectively. The largest difference between DVK convolution and full MC dosimetry was observed for the lung tissue (10.16% ± 1.20%). The DVK statistical uncertainty was <0.75% for both media. Conclusions: This semiautomatic algorithm is capable of performing rapid, accurate, and efficient 3D dosimetry. The proposed method considers tissue and activity heterogeneity using tissue-specific DVKs. Furthermore, this method provides results in <1 min, making it suitable for everyday clinical practice.

High quality imaging and dosimetry for yttrium-90 (90Y) liver radioembolization using a SiPM-based PET/CT scanner

PURPOSE

Transarterial radioembolization (TARE) with yttrium-90 (⁹⁰Y) microspheres is a liver-directed treatment for primary and secondary hepatic malignancies. Personalized dosimetry aims for maximum treatment effect and reduced toxicity. We aimed to compare pre-treatment voxel-based dosimetry from ^99mTc macroaggregated albumin (MAA) SPECT/CT with post-treatment ⁹⁰Y PET/CT for absorbed dose values, and to evaluate image quality of ⁹⁰Y SiPM-based PET/CT.

METHODS

Forty-two patients (28 men, 14 women, mean age: 67 ± 11 years) with advanced hepatic malignancies were prospectively enrolled. Twenty patients were treated with glass and 22 with resin microspheres. Radiation absorbed doses from planning ^99mTc-MAA SPECT/CT and post-therapy ⁹⁰Y PET/CT were assessed. ⁹⁰Y PET/CT images were acquired for 20 min and reconstructed to produce 5-, 10-, 15-, and 20-min datasets, then evaluated using the 5-point Likert scale.

RESULTS

The mean administered activity was 3.44 ± 1.5 GBq for glass and 1.62 ± 0.7 GBq for resin microspheres. The mean tumor absorbed doses calculated from ^99mTc-MAA SPECT/CT and ⁹⁰Y PET/CT were 175.69 ± 113.76 Gy and 193.58 ± 111.09 Gy (P = 0.61), respectively for glass microspheres; they were 60.18 ± 42.20 Gy and 70.98 ± 49.65 Gy (P = 0.37), respectively for resin microspheres. The mean normal liver absorbed doses from ^99mTc-MAA SPECT/CT and ⁹⁰Y PET/CT were 32.70 ± 22.25 Gy and 30.62 ± 20.09 Gy (P = 0.77), respectively for glass microspheres; they were 18.33 ± 11.08 Gy and 24.32 ± 15.58 Gy (P = 0.17), respectively for resin microspheres. Image quality of ⁹⁰Y PET/CT at 5-, 10-, 15-, and 20-min scan time showed a Likert score of 3.6 ± 0.54, 4.57 ± 0.58, 4.84 ± 0.37, and 4.9 ± 0.3, respectively.

CONCLUSIONS

^99mTc-MAA SPECT/CT demonstrated great accuracy for treatment planning dosimetry. SiPM-based PET/CT scanner showed good image quality at 10-min scan time, acquired in one bed position. A PET/CT scan time of 5 min showed acceptable image quality and suffices for dosimetry and treatment verification. This allows for inclusion of ⁹⁰Y PET/CT in busy routine clinical workflows. Studies with larger patient cohorts are needed to confirm these findings.

Lung shunt and lung dose calculation methods for radioembolization treatment planning

Radioembolization, also known as selective internal radiation therapy (SIRT), is firmly established in the management of patients with unresectable liver cancers. Advances in normal and tumor liver dosimetry and new knowledge about tumor dose response relationships have helped promote the safe use of higher prescribed doses, consequently transitioning radioembolization from palliative to curative therapy. The lungs are considered a critical organ of risk for radioembolization treatment planning. Unfortunately, lung dosimetry has not achieved similar advances in dose calculation methodology as liver dosimetry. Current estimations of lung dose are dependent on a number of parameters associated with data acquisition and processing algorithms, leading to poor accuracy and precision. Therefore, the efficacy of curative radioembolization may be compromised in patients for whom the lung dose derived using currently available methods unnecessarily limits the desired administered activity to the liver. We present a systematic review of the various methods of determining the lung shunt fraction (LSF) and lung mean dose (LD). This review encompasses pretherapy estimations and post-therapy assessments of the LSF and LD using both 2D planar and 3D SPECT/CT based calculations. The advantages and limitations of each of these methods are deliberated with a focus on accuracy and practical considerations. We conclude the review by presenting a lexicon to precisely describe the methodology used for the estimation of LSF and LD; specifically, category, agent, modality, contour and algorithm, in order to aid in their interpretation and standardization in routine clinical practice.

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.

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.

Transarterial Radioembolization for Hepatocellular Carcinoma and Hepatic Metastases: Clinical Aspects and Dosimetry Models

Transarterial radioembolization (TARE) with Yttrium-90 (⁹⁰Y) microspheres is a liver-directed therapy for primary and metastatic disease. This manuscript provides a review of the clinical literature on TARE indications and efficacy with overviews of patient-selection and toxicity. Current dosimetry models used in practice are safe, relatively simple, and easy for clinicians to use. Planning currently relies on the imperfect surrogate, ^99mTc macroaggregated albumin. Post-therapy quantitative imaging (⁹⁰Y SPECT/CT or ⁹⁰Y PET/CT) of microspheres can be used to calculate the macroscopic in vivo absorbed dose distribution. Similar to the evolution of other brachytherapy dose calculations, TARE is moving toward more patient-specific dosimetry that includes calculating and reporting nonuniform dose distributions throughout tumors and normal uninvolved liver.

Respiratory motion compensation in interventional liver SPECT using simultaneous fluoroscopic and nuclear imaging

PURPOSE

Quantitative accuracy of the single photon emission computed tomography (SPECT) reconstruction of the pretreatment procedure of liver radioembolization is crucial for dosimetry; visual quality is important for detecting doses deposited outside the planned treatment volume. Quantitative accuracy is limited by respiratory motion. Conventional gating eliminates motion by count rejection but increases noise, which degrades the visual reconstruction quality. Motion compensation using all counts can be performed if the motion signal and motion vector field over time are known. The measurement of the motion signal of a patient currently requires a device (such as a respiratory belt) attached to the patient, which complicates the acquisition. The motion vector field is generally extracted from a previously acquired four-dimensional scan and can differ from the motion in the scan performed during the intervention. The simultaneous acquisition of fluoroscopic and nuclear projections can be used to obtain both the motion vector field and the projections of the corresponding (moving) activity distribution. This eliminates the need for devices attached to the patient and provides an accurate motion vector field for SPECT reconstruction. Our approach to motion compensation would primarily be beneficial for interventional SPECT because the time-critical setting requires fast scans and no inconvenience of an external apparatus. The purpose of this work is to evaluate the performance of the motion compensation approach for interventional liver SPECT by means of simulations.

METHODS

Nuclear and fluoroscopic projections of a realistic digital human phantom with respiratory motion were generated using fast Monte Carlo simulators. Fluoroscopic projections were sampled at 1-5 Hz. Nuclear data were acquired continuously in list mode. The motion signal was extracted from the fluoroscopic projections by calculating the center-of-mass, which was then used to assign each photon to a corresponding motion bin. The fluoroscopic projections were reconstructed per bin and coregistered, resulting in a motion vector field that was used in the SPECT reconstruction. The influence of breathing patterns, fluoroscopic imaging dose, sampling rate, number of bins, and scanning time was studied. In addition, the motion compensation method was compared with conventional gating to evaluate the detectability of spheres with varying uptake ratios.

RESULTS

The liver motion signal was accurately extracted from the fluoroscopic projections, provided the motion was stable in amplitude and the sampling rate was greater than 2 Hz. The minimum total fluoroscopic dose for the proposed method to function in a 5-min scan was 10 µGy. Although conventional gating improved the quantitative reconstruction accuracy, substantial background noise was observed in the short scans because of the limited counts available. The proposed method similarly improved the quantitative accuracy, but generated reconstructions with higher visual quality. The proposed method provided better visualization of low-contrast features than when using gating.

CONCLUSION

The proposed motion compensation method has the potential to improve SPECT reconstruction quality. The method eliminates the need for external devices to measure the motion signal and generates an accurate motion vector field for reconstruction. A minimal increase in the fluoroscopic dose is required to substantially improve the results, paving the way for clinical use.

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.

Assessing Spatial Concordance Between Theranostic Pairs Using Phantom and Patient-Specific Acceptance Criteria: Application to 99mTc-MAA SPECT/90Y-Microsphere PET

PURPOSE

Predictive 3-dimensional dosimetry requires spatial concordance between diagnostic and therapeutic activity distributions. We assess similarity between theranostic pairs (^99mTc-macroaggregated albumin [MAA] single photon emission computed tomography [SPECT] and ⁹⁰Y microsphere positron emission tomography [PET]) in patients using criteria that account for spatial resolution differences and misregistration.

METHODS AND MATERIALS

Phantom-based acceptance criteria were determined using a liver phantom filled with ^99mTc and ⁹⁰YCl₃ and scanned with SPECT/computed tomography [CT] and PET/CT, respectively. Gaussian blurring was applied to PET to match ^99mTc phantom scan image quality. After rigid registration between SPECT/CT and PET/CT, perturbations up to ±3 voxels were applied to determine the similarity metric (SM) sensitivity. ^99mTc-MAA SPECT/CT and ⁹⁰Y microsphere PET/CT image pairs/patients (n = 23) were processed analogously. SMs calculated included the Pearson correlation coefficient (ρ_r), Lin's concordance correlation coefficient (ρ_c), Spearman's rank correlation coefficient (ρ_s), the mean squared difference, and the Dice similarity coefficient (DSC). Patient-specific acceptance criteria were determined by evaluating the SMs of the blurred PET compared with itself misregistered.

RESULTS

After transforming PET to SPECT resolution, high similarity was found in phantom, with ρ_c, ρ_r, ρ_s > 0.98 ± 0.01, a mean squared difference of (4.1 ± 0.3) × 10^-4 and DSC > 0.85 ± 0.01 for investigated thresholds (5%, 30%, and 50%). SMs for patients varied from poor to good. A small percentage (13%-30%) of patient scans were acceptable using phantom-based acceptance criteria. The percentage increased slightly (17%-35%) using patient-specific acceptance criteria. DSC for most patients were substantially lower (average 0.95 vs 0.61 for 5% threshold) than phantom values.

CONCLUSIONS

At best, 35% of patients had an SM within the acceptance criteria established to account for imaging-related effects impacting spatial concordance between ^99mTc-MAA SPECT and ⁹⁰Y PET. Additional clinical factors should be evaluated in the future. The procedure of accounting for image-related effects when assessing spatial concordance can be applied to other theranostic pairs.

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.

Multi-modal image analysis for semi-automatic segmentation of the total liver and liver arterial perfusion territories for radioembolization

Purpose

We have developed a multi-modal imaging approach for SIRT, combining ^99mTc-MAA SPECT/CT and/or ⁹⁰Y PET, ¹⁸F-FDG PET/CT, and contrast-enhanced CBCT for voxel-based dosimetry, as a tool for treatment planning and verification. For radiation dose prediction calculations, a segmentation of the total liver volume and of the liver perfusion territories is required.

Method

In this paper, we proposed a procedure for multi-modal image analysis to assist SIRT treatment planning. The pre-treatment ¹⁸F-FDG PET/CT, ^99mTc-MAA SPECT/CT, and contrast-enhanced CBCT images were registered to a common space using an initial rigid, followed by a deformable registration. The registration was scored by an expert using Likert scores. The total liver was segmented semi-automatically based on the PET/CT and SPECT/CT images, and the liver perfusion territories were determined based on the CBCT images. The segmentations of the liver and liver lobes were compared to the manual segmentations by an expert on a CT image.

Result

Our methodology showed that multi-modal image analysis can be used for determination of the liver and perfusion territories using CBCT in SIRT using all pre-treatment studies. The results for image registration showed acceptable alignment with limited impact on dosimetry.

The image registration performs well according to the expert reviewer (scored as perfect or with little misalignment in 94% of the cases). The semi-automatic liver segmentation agreed well with manual liver segmentation (dice coefficient of 0.92 and an average Hausdorff distance of 3.04 mm). The results showed disagreement between lobe segmentation using CBCT images compared to lobe segmentation based on CT images (average Hausdorff distance of 14.18 mm), with a high impact on the dosimetry (differences up to 9 Gy for right and 21 Gy for the left liver lobe).

Conclusion

This methodology can be used for pre-treatment dosimetry and for SIRT planning including the determination of the activity that should be administered to achieve the therapeutical goal. The inclusion of perfusion CBCT enables perfusion-based definition of the liver lobes, which was shown to be markedly different from the anatomical definition in some of the patients.

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.

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.

Characterization of 90 Y-SPECT/CT self-calibration approaches on the quantification of voxel-level absorbed doses following 90 Y-microsphere selective internal radiation therapy

PURPOSE

⁹⁰ Y-microsphere selective internal radiation therapy (⁹⁰ Y-SIRT or ⁹⁰ Y-radioembolization) is used in the management of unresectable liver tumors. ⁹⁰ Y-SIRT presents a unique situation where the total ⁹⁰ Y activity inside the liver can be determined with high accuracy (> 95%). ⁹⁰ Y bremsstrahlung single-photon emission computed tomography (SPECT)/computed tomography (CT) can be self-calibrated to provide quantitative images that facilitate voxel-level absorbed dose calculations. We investigated the effects of different approaches for ⁹⁰ Y-SPECT self-calibration on the quantification of absorbed doses following ⁹⁰ Y-SIRT.

METHODS

⁹⁰ Y bremsstrahlung SPECT/CT images of 31 patients with hepatocellular carcinoma, collected following ⁹⁰ Y-SIRT, were analyzed, yielding 48 tumor and 31 normal liver contours. We validated the accuracy of absorbed doses calculated by a commercial software against those calculated using Monte Carlo-based radiation transport. The software package was used to analyze the following definitions of SPECT volume of interest used for ⁹⁰ Y-SPECT self-calibration: (a) SPECT field-of-view (FOV), (b) chest-abdomen contour, (c) total liver contour, (d) total liver contour expanded by 5 mm, and (e) total liver contour contracted by 5 mm. Linear correlation and Bland-Altman analysis were performed for tumor and normal liver tissue absorbed dose volume histogram metrics between the five different approaches for ⁹⁰ Y-SPECT self-calibration.

RESULTS

The mean dose calculated using the commercial software was within 3% of Monte Carlo for tumors and normal liver tissues. The tumor mean dose calculated using the chest-abdomen calibration was within 2% of that calculated using the SPECT FOV, whereas the doses calculated using the total liver contour, expanded total liver contour, and contracted total liver contour were within 68%, 47%, and 107%, respectively, of doses calculated using the SPECT FOV. The normal liver tissue mean dose calculated using the chest-abdomen contour was within 1.3% of that calculated using the SPECT FOV, whereas the doses calculated using the total liver contour, expanded total liver contour, and contracted total liver contour were within 73%, 50%, and 114%, respectively, of doses calculated using the SPECT FOV.

CONCLUSIONS

The mean error of < 3% for commercial software can be considered clinically acceptable for ⁹⁰ Y-SIRT dosimetry. Absorbed dose quantification using ⁹⁰ Y-SPECT self-calibration with the chest-abdomen contour was equivalent to that calculated using the SPECT FOV, but self-calibration with the total liver contour yielded substantially higher (~70%) dose values. The large biases revealed by our study suggest that consistent absorbed dose calculation approaches are essential when comparing ⁹⁰ Y-SIRT dosimetry between different clinical studies.

Phantom validation of quantitative Y-90 PET/CT-based dosimetry in liver radioembolization

Background

PET/CT has recently been shown to be a viable alternative to traditional post-infusion imaging methods providing good quality images of ⁹⁰Y-laden microspheres after selective internal radiation therapy (SIRT). In the present paper, first we assessed the quantitative accuracy of ⁹⁰Y-PET using an anthropomorphic phantom provided with lungs, liver, spine, and a cylindrical homemade lesion located into the hepatic compartment. Then, we explored the accuracy of different computational approaches on dose calculation, including (I) direct Monte Carlo radiation transport using Raydose, (II) Kernel convolution using Philips Stratos, (III) local deposition algorithm, (IV) Monte Carlo technique (MCNP) considering a uniform activity distribution, and (V) MIRD (Medical Internal Radiation Dose) analytical approach. Finally, calculated absorbed doses were compared with those obtained performing measurements with LiF:Mg,Cu,P TLD chips in a liquid environment.

Results

Our results indicate that despite ⁹⁰Y-PET being likely to provide high-resolution images, the ⁹⁰Y low branch ratio, along with other image-degrading factors, may produce non-uniform activity maps, even in the presence of uniform activity. A systematic underestimation of the recovered activity, both for the tumor insert and for the liver background, was found. This is particularly true if no partial volume correction is applied through recovery coefficients. All dose algorithms performed well, the worst case scenario providing an agreement between absorbed dose evaluations within 20%. Average absorbed doses determined with the local deposition method are in excellent agreement with those obtained using the MIRD and the kernel-convolution dose calculation approach.

Finally, absorbed dose assessed with MC codes are in good agreement with those obtained using TLD in liquid solution, thus confirming the soundness of both calculation approaches. This is especially true for Raydose, which provided an absorbed dose value within 3% of the measured dose, well within the stated uncertainties.

Conclusions

Patient-specific dosimetry is possible even in a scenario with low true coincidences and high random fraction, as in ⁹⁰Y–PET imaging, granted that accurate absolute PET calibration is performed and acquisition times are sufficiently long. Despite Monte Carlo calculations seeming to outperform all dose estimation algorithms, our data provide a strong argument for encouraging the use of the local deposition algorithm for routine ⁹⁰Y dosimetry based on PET/CT imaging, due to its simplicity of implementation.

Impact of respiratory motion and acquisition settings on SPECT liver dosimetry for radioembolization

PURPOSE

Respiratory motion may impose significant inaccuracies on emission activity estimation in quantitative SPECT. This effect may be a major issue in dosimetry as used in the management of liver radioembolization. The purpose of this study was to investigate the impact of respiratory motion on radioembolization liver dosimetry for different SPECT acquisition settings.

METHODS

In a series of SPECT/CT Monte Carlo simulations using several digital XCAT phantoms, the following parameters were varied: breathing/nonbreathing, liver tumor size (0.3-35 ml) and location, patient properties (body mass index ranging from underweight to obese; male and female), acquisition time (10-30 s/view), collimator setup (High Sensitivity, High Resolution, Ultra High Resolution) and tumor VOI. The effect of applying a respiratory gating scheme was examined as well.

RESULTS

Breathing decreased activity recovery and tumor/non-tumor (T/N) ratios on average from 90% to 66%. VOIs based on SPECT images instead of breath-hold CT improved T/N values significantly. The most accurate results were obtained using a gating scheme combined with SPECT-based VOIs. Scan duration, body mass index, sex, and location all had a minor effect. Lung shunt fraction estimations were relatively unaffected by any of the varied parameters.

CONCLUSIONS

Respiratory motion has a large effect on SPECT activity quantitation of liver tumors as used in radioembolization treatment planning and assessment. As compared with the other parameters that were varied in this study, respiration is the predominant degrading effect on image quantitation. Gating alleviates much of this detrimental effect.

Practical reconstruction protocol for quantitative (90)Y bremsstrahlung SPECT/CT

PURPOSE

To develop a practical background compensation (BC) technique to improve quantitative (90)Y-bremsstrahlung single-photon emission computed tomography (SPECT)/computed tomography (CT) using a commercially available imaging system.

METHODS

All images were acquired using medium-energy collimation in six energy windows (EWs), ranging from 70 to 410 keV. The EWs were determined based on the signal-to-background ratio in planar images of an acrylic phantom of different thicknesses (2-16 cm) positioned below a (90)Y source and set at different distances (15-35 cm) from a gamma camera. The authors adapted the widely used EW-based scatter-correction technique by modeling the BC as scaled images. The BC EW was determined empirically in SPECT/CT studies using an IEC phantom based on the sphere activity recovery and residual activity in the cold lung insert. The scaling factor was calculated from 20 clinical planar (90)Y images. Reconstruction parameters were optimized in the same SPECT images for improved image quantification and contrast. A count-to-activity calibration factor was calculated from 30 clinical (90)Y images.

RESULTS

The authors found that the most appropriate imaging EW range was 90-125 keV. BC was modeled as 0.53× images in the EW of 310-410 keV. The background-compensated clinical images had higher image contrast than uncompensated images. The maximum deviation of their SPECT calibration in clinical studies was lowest (<10%) for SPECT with attenuation correction (AC) and SPECT with AC + BC. Using the proposed SPECT-with-AC + BC reconstruction protocol, the authors found that the recovery coefficient of a 37-mm sphere (in a 10-mm volume of interest) increased from 39% to 90% and that the residual activity in the lung insert decreased from 44% to 14% over that of SPECT images with AC alone.

CONCLUSIONS

The proposed EW-based BC model was developed for (90)Y bremsstrahlung imaging. SPECT with AC + BC gave improved lesion detectability and activity quantification compared to SPECT with AC only. The proposed methodology can readily be used to tailor (90)Y SPECT/CT acquisition and reconstruction protocols with different SPECT/CT systems for quantification and improved image quality in clinical settings.

Evaluating the Application of Tissue-Specific Dose Kernels Instead of Water Dose Kernels in Internal Dosimetry: A Monte Carlo Study

PURPOSE

The aim of this work is to evaluate the application of tissue-specific dose kernels instead of water dose kernels to improve the accuracy of patient-specific dosimetry by taking tissue heterogeneities into consideration.

MATERIALS AND METHODS

Tissue-specific dose point kernels (DPKs) and dose voxel kernels (DVKs) for yttrium-90 (⁹⁰Y), lutetium-177 (¹⁷⁷Lu), and phosphorus-32 (³²P) are calculated using the Monte Carlo (MC) simulation code GATE (version 7). The calculated DPKs for bone, lung, adipose, breast, heart, intestine, kidney, liver, and spleen are compared with those of water. The dose distribution in normal and tumorous tissues in lung, liver, and bone of a Zubal phantom is calculated using tissue-specific DVKs instead of those of water in conventional methods. For a tumor defined in a heterogeneous region in the Zubal phantom, the absorbed dose is calculated using a proposed algorithm, taking tissue heterogeneity into account. The algorithm is validated against full MC simulations.

RESULTS

The simulation results indicate that the highest differences between water and other tissue DPKs occur in bone for ⁹⁰Y (12.2% ± 0.6%), ³²P (18.8% ± 1.3%), and ¹⁷⁷Lu (16.9% ± 1.3%). The second highest discrepancy corresponds to the lung for ⁹⁰Y (6.3% ± 0.2%), ³²P (8.9% ± 0.4%), and ¹⁷⁷Lu (7.7% ± 0.3%). For ⁹⁰Y, the mean absorbed dose in tumorous and normal tissues is calculated using tissue-specific DVKs in lung, liver, and bone. The results are compared with doses calculated considering the Zubal phantom water equivalent and the relative differences are 4.50%, 0.73%, and 12.23%, respectively. For the tumor in the heterogeneous region of the Zubal phantom that includes lung, liver, and bone, the relative difference between mean calculated dose in tumorous and normal tissues based on the proposed algorithm and the values obtained from full MC dosimetry is 5.18%.

CONCLUSIONS

A novel technique is proposed considering tissue-specific dose kernels in the dose calculation algorithm. This algorithm potentially enables patient-specific dosimetry and improves estimation of the average absorbed dose of ⁹⁰Y in a tumor located in lung, bone, and soft tissue interface by 6.98% compared with the conventional methods.

Selective Internal Radiation Therapy With Yttrium-90 Glass Microspheres: Biases and Uncertainties in Absorbed Dose Calculations Between Clinical Dosimetry Models

PURPOSE

To quantify differences that exist between dosimetry models used for ⁹⁰Y selective internal radiation therapy (SIRT).

METHODS AND MATERIALS

Retrospectively, 37 tumors were delineated on 19 post-therapy quantitative ⁹⁰Y single photon emission computed tomography/computed tomography scans. Using matched volumes of interest (VOIs), absorbed doses were reported using 3 dosimetry models: glass microsphere package insert standard model (SM), partition model (PM), and Monte Carlo (MC). Univariate linear regressions were performed to predict mean MC from SM and PM. Analysis was performed for 2 subsets: cases with a single tumor delineated (best case for PM), and cases with multiple tumors delineated (typical clinical scenario). Variability in PM from the ad hoc placement of a single spherical VOI to estimate the entire normal liver activity concentration for tumor (T) to nontumoral liver (NL) ratios (TNR) was investigated. We interpreted the slope of the resulting regression as bias and the 95% prediction interval (95%PI) as uncertainty. MC_NL^single represents MC absorbed doses to the NL for the single tumor patient subset; other combinations of calculations follow a similar naming convention.

RESULTS

SM was unable to predict MC_T^single or MC_T^multiple (p>.12, 95%PI >±177 Gy). However, SM^single was able to predict (p<.012) MC_NL^single, albeit with large uncertainties; SM^single and SM^multiple yielded biases of 0.62 and 0.71, and 95%PI of ±40 and ± 32 Gy, respectively. PM_T^single and PM_T^multiple predicted (p<2E-6) MC_T^single and MC_T^multiple with biases of 0.52 and 0.54, and 95%PI of ±38 and ± 111 Gy, respectively. The TNR variability in PM_T^single increased the 95%PI for predicting MC_T^single (bias = 0.46 and 95%PI = ±103 Gy). The TNR variability in PM_T^multiple modified the bias when predicting MC_T^multiple (bias = 0.32 and 95%PI = ±110 Gy).

CONCLUSIONS

The SM is unable to predict mean MC tumor absorbed dose. The PM is statistically correlated with mean MC, but the resulting uncertainties in predicted MC are large. Large differences observed between dosimetry models for ⁹⁰Y SIRT warrant caution when interpreting published SIRT absorbed doses. To reduce uncertainty, we suggest the entire NL VOI be used for TNR estimates when using PM.