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

An Evaluation of 90Y Bremsstrahlung SPECT Image Quality in the Presence of 99mTc: A Technical Perspective on Same-Day Radioembolization

In same-day radioembolization, 99mTc-MAA SPECT/CT, 90Y radioembolization, and post-treatment 90Y SPECT/CT procedures are conducted on the same-day, resulting in a dual-isotope environment of 90Y and 99mTc during post-treatment imaging. This study aimed to quantify the impact of 99mTc on 90Y bremsstrahlung-SPECT/CT image quality and to establish an optimised imaging protocol for both clinical practice, and with advanced reconstruction techniques. Utilising a NEMA IQ phantom, contrast recovery coefficients (CRCs) were measured to evaluate the 90Y image quality degradation caused by 99mTc. SPECT/CT scans of 90Y-only and 90Y with varying amounts of 99mTc were conducted using a standard protocol (90-120 keV energy window, high-energy collimator) and various dual-isotope protocols. The standard protocol resulted in a marked CRC reduction, with the largest sphere's CRC decreasing from 0.21 (90Y-only) to 0.05 when 99mTc activity was 5% of 90Y. For an optimised protocol (160-200 keV energy window, high-energy collimator) CRC values were 0.16 for 90Y-only and 0.15 for 90Y+99mTc. The highest CRC values were achieved with an advanced Monte Carlo-based reconstruction, showing 0.58 for 90Y-only and 0.46 for 90Y+99mTc. Image quality degradation was noted in dual-isotope settings even when using an optimised protocol. Advanced reconstruction techniques markedly improved post-treatment image quality.

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.

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.

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.

EANM guidance document: dosimetry for first-in-human studies and early phase clinical trials

The numbers of diagnostic and therapeutic nuclear medicine agents under investigation are rapidly increasing. Both novel emitters and novel carrier molecules require careful selection of measurement procedures. This document provides guidance relevant to dosimetry for first-in human and early phase clinical trials of such novel agents. The guideline includes a short introduction to different emitters and carrier molecules, followed by recommendations on the methods for activity measurement, pharmacokinetic analyses, as well as absorbed dose calculations and uncertainty analyses. The optimal use of preclinical information and studies involving diagnostic analogues is discussed. Good practice reporting is emphasised, and relevant dosimetry parameters and method descriptions to be included are listed. Three examples of first-in-human dosimetry studies, both for diagnostic tracers and radionuclide therapies, are given.

Energy window optimization in bremsstrahlung imaging after Yttrium-90 microsphere therapy

In imaging of Yttrium-90 patients treated hepatic primary and metastatic cancers, bremsstrahlung photons produced in a wide energy range is used. However, the image quality depends on acquisition energy window. This research aimed energy window optimization for Yttrium-90 bremsstrahlung imaging and 48 patients with various types of cancer received radioembolization therapy were investigated. Patients were imaged using a GE Healthcare Optima NM/CT 640 series gamma camera system with a medium energy general-purpose (MEGP) collimator and planar images were acquired with 8 different energy windows in the 55-400 keV energy range. The data set, formed with the % FOV, contrast, and spatial resolution of image quality parameters calculated from these images, was statistically examined with ANOVA and Tukey tests. According to the visual evaluations and ANOVA/Tukey test results, it was statistically concluded that energy window of 90-110 keV is the optimal energy window while 60-400 keV energy ranges show the lowest image quality for Y-90 bremsstrahlung imaging.

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.

90Y SPECT scatter estimation and voxel dosimetry in radioembolization using a unified deep learning framework

PURPOSE

90Y SPECT-based dosimetry following radioembolization (RE) in liver malignancies is challenging due to the inherent scatter and the poor spatial resolution of bremsstrahlung SPECT. This study explores a deep-learning-based absorbed dose-rate estimation method for 90Y that mitigates the impact of poor SPECT image quality on dosimetry and the accuracy-efficiency trade-off of Monte Carlo (MC)-based scatter estimation and voxel dosimetry methods.

METHODS

Our unified framework consists of three stages: convolutional neural network (CNN)-based bremsstrahlung scatter estimation, SPECT reconstruction with scatter correction (SC) and absorbed dose-rate map generation with a residual learning network (DblurDoseNet). The input to the framework is the measured SPECT projections and CT, and the output is the absorbed dose-rate map. For training and testing under realistic conditions, we generated a series of virtual patient phantom activity/density maps from post-therapy images of patients treated with 90Y-RE at our clinic. To train the scatter estimation network, we use the scatter projections for phantoms generated from MC simulation as the ground truth (GT). To train the dosimetry network, we use MC dose-rate maps generated directly from the activity/density maps of phantoms as the GT (Phantom + MC Dose). We compared performance of our framework (SPECT w/CNN SC + DblurDoseNet) and MC dosimetry (SPECT w/CNN SC + MC Dose) using normalized root mean square error (NRMSE) and normalized mean absolute error (NMAE) relative to GT.

RESULTS

When testing on virtual patient phantoms, our CNN predicted scatter projections had NRMSE of 4.0% ± 0.7% on average. For the SPECT reconstruction with CNN SC, we observed a significant improvement on NRMSE (9.2% ± 1.7%), compared to reconstructions with no SC (149.5% ± 31.2%). In terms of virtual patient dose-rate estimation, SPECT w/CNN SC + DblurDoseNet had a NMAE of 8.6% ± 5.7% and 5.4% ± 4.8% in lesions and healthy livers, respectively; compared to 24.0% ± 6.1% and 17.7% ± 2.1% for SPECT w/CNN SC + MC Dose. In patient dose-rate maps, though no GT was available, we observed sharper lesion boundaries and increased lesion-to-background ratios with our framework. For a typical patient data set, the trained networks took ~ 1 s to generate the scatter estimate and ~ 20 s to generate the dose-rate map (matrix size: 512 × 512 × 194) on a single GPU (NVIDIA V100).

CONCLUSION

Our deep learning framework, trained using true activity/density maps, has the potential to outperform non-learning voxel dosimetry methods such as MC that are dependent on SPECT image quality. Across comprehensive testing and evaluations on multiple targeted lesions and healthy livers in virtual patients, our proposed deep learning framework demonstrated higher (66% on average in terms of NMAE) estimation accuracy than the current "gold-standard" MC method. The enhanced computing speed with our framework without sacrificing accuracy is highly relevant for clinical dosimetry following 90Y-RE.

Predicting the net administered activity in 90 Y-radioembolization patients from post-procedure 90 Y-SPECT/CT

BACKGROUND

The calculation of the net administered activity (Aadmin ) in patients undergoing 90 Y-radioembolization is essential for dosimetry and radiation safety, yet current methods for measuring residual 90 Y activity are often associated with high uncertainty. Therefore, an accurate, robust, and clinically viable method for the determination of Aadmin across approved 90 Y microsphere devices is desirable.

PURPOSE

We report on a novel method to determine Aadmin by leveraging the quantitative capabilities of SPECT/CT to measure 90 Y-emission in vivo from patients following 90 Y-radioembolization with glass or resin microspheres.

METHODS

90 Y-SPECT/CT attenuation-corrected count data from 147 sequential 90 Y-radioembolization patients was used for this analysis. Aadmin was calculated as part of routine clinical practice via the exposure rate differences between the initial 90 Y-vial and the 90 Y-residual jar. This served as our gold standard measure of Aadmin . Patient data for each microsphere device were separated into training and testing cohorts to first develop regression models and then to independently assess model performance. The training cohorts were divided into four groups: first, based on the microsphere device (glass or resin), and second, based on the SPECT volume used to calculate counts (the full SPECT field of view (FOV) or liver only (VOIliver )). Univariate linear regression models were generated for each group to predict Aadmin based on 90 Y-SPECT data from the training cohorts. Leave-one-out cross validation was implemented to estimate variability in model parameters. To assess performance, linear models derived from the training cohort were applied to 90 Y-SPECT data from the testing cohort. A comparison of the models between microspheres devices was also performed.

RESULTS

Linear models derived from the glass and resin training cohorts demonstrated a strong, positive correlation between 90 Y-SPECT image counts and Aadmin for VOIliver and FOV with R2  > 0.98 in all cases. In the glass training cohort, model accuracy (100%-absolute mean prediction error) and precision (95% prediction intervals of mean prediction error) were 99.0% and 15.4% for the VOIliver and 99.7% and 17.5% for the FOV models, respectively. In the resin training cohort, the corresponding values were 98.6% and 16.7% for VOIliver and > 99.9% and 11.4% for the FOV models, respectively. The application of these linear models to 90 Y-SPECT data from the testing cohort showed Aadmin prediction errors to have high accuracy and precision for both microsphere devices. For the glass testing cohort, accuracy (precision) was 96.9% (19.6%) and 98.8% (21.1%) for the VOIliver and FOV models, respectively. The corresponding values for the resin training cohort were 97.3% (26.2%) and 98.5% (25.7%) for the VOIliver and FOV models, respectively. The slope of the linear models between the two microsphere devices was observed to be significantly different with resin microspheres generating 48%-49% more SPECT counts for equivalent 90 Y activity based on each device manufacturer's activity calibration process.

CONCLUSION

90 Y-SPECT image counts can reliably predict (accuracy > 95% and precision < 18%) Aadmin after 90 Y-radioembolization, with performance characteristics essentially equivalent for both glass and resin microspheres. There is a clear indication that activity calibrations are fundamentally different between the two microsphere devices.

For Hepatocellular Carcinoma Treated with Yttrium-90 Microspheres, Dose Volumetrics on Post-Treatment Bremsstrahlung SPECT/CT Predict Clinical Outcomes

In transarterial radioembolization (TARE) of hepatocellular carcinoma (HCC) with Yttrium-90 (Y-90) microspheres, recent studies correlate dosimetry from bremsstrahlung single photon emission tomography (SPECT/CT) with treatment outcomes; however, these studies focus on measures of central tendency rather than volumetric coverage metrics commonly used in radiation oncology. We hypothesized that three-dimensional (3D) isodose coverage of gross tumor volume (GTV) is the driving factor in HCC treatment response to TARE and is best assessed using advanced dosimetry techniques applied to nuclear imaging of actual Y-90 biodistribution. We reviewed 51 lobar TARE Y-90 treatments of 43 HCC patients. Dose prescriptions were 120 Gy for TheraSpheres and 85 Gy for SIR-Spheres. All patients underwent post-TARE Y-90 bremsstrahlung SPECT/CT imaging. Commercial software was used to contour gross tumor volume (GTV) and liver on post-TARE SPECT/CT. Y-90 dose distributions were calculated using the Local Deposition Model based on post-TARE SPECT/CT activity maps. Median gross tumor volume (GTV) dose; GTV receiving less than 100 Gy, 70 Gy and 50 Gy; minimum dose covering the hottest 70%, 95%, and 98% of the GTV (D70, D95, D98); mean dose to nontumorous liver, and disease burden (GTV/liver volume) were obtained. Clinical outcomes were collected for all patients by chart and imaging review. HCC treatment response was assessed according to the modified response criteria in solid tumors (mRECIST) guidelines. Kaplan-Meier (KM) survival estimates and multivariate regression analyses (MVA) were performed using STATA. Median survival was 22.5 months for patients achieving objective response (OR) in targeted lesions (complete response (CR) or partial response (PR) per mRECIST) vs. 7.6 months for non-responders (NR, stable disease or disease progression per mRECIST). On MVA, the volume of underdosed tumor (GTV receiving less than 100 Gy) was the only significant dosimetric predictor for CR (p = 0.0004) and overall survival (OS, p = 0.003). All targets with less than CR (n = 39) had more than 20 cc of underdosed tumor. D70 (p = 0.038) correlated with OR, with mean D70 of 95 Gy for responders and 60 Gy for non-responders (p = 0.042). On MVA, mean dose to nontumorous liver trended toward significant association with grade 3+ toxicity (p = 0.09) and correlated with delivered activity (p < 0.001) and burden of disease (p = 0.05). Dosimetric models supplied area under the curve estimates of > 0.80 predicting CR, OR, and ≥grade 3 acute toxicity. Dosimetric parameters derived from the retrospective analysis of post-TARE Y-90 bremsstrahlung SPECT/CT after lobar treatment of HCC suggest that volumetric coverage of GTV, not a high mean or median dose, is the driving factor in treatment response and that this is best assessed through the analysis of actual Y-90 biodistribution.

Outcomes following transarterial radioembolization with 90Y and nanoparticles loaded resin microspheres

⁹⁰Y bremsstrahlung Single-Photon Emission Tomography (SPECT) imaging is employed to check the possibility of extrahepatic uptake and the quantification of delivered dose in Transarterial Radioembolization (TARE). ⁹⁰Y bremsstrahlung SPECT imaging is challenging due to the nature of bremsstrahlung photons. We reported a Monte Carlo study using the resin microspheres loaded with ⁹⁰Y and Nanoparticles (NPs) in the TARE. By injection of Bismuth (Bi) and Europium (Eu) NPs into the resin microspheres, the sensitivity and the contrast to noise ratio increased for the bremsstrahlung planar images. The highest signal to background ratio was observed in the characteristic X-ray planar images taken with the energy window at the Kα₁ ± 10 keV when Eu NPs were incorporated into the microsphere. The dose enhancement ratio decreased dramatically at NP concentrations >2.4 M. Incorporating NPs into the resin microspheres improves the quality of post-treatment images and establishes a standardized imaging protocol for post-treatment imaging by characteristic X-rays.

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.

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.

Hybrid kernelised expectation maximisation for Bremsstrahlung SPECT reconstruction in SIRT with 90Y micro-spheres

Background

Selective internal radiation therapy with Yttrium-90 microspheres is an effective therapy for liver cancer and liver metastases. Yttrium-90 is mainly a high-energy beta particle emitter. These beta particles emit Bremsstrahlung radiation during their interaction with tissue making post-therapy imaging of the radioactivity distribution feasible. Nevertheless, image quality and quantification is difficult due to the continuous energy spectrum which makes resolution modelling, attenuation and scatter estimation challenging and therefore the dosimetry quantification is inaccurate. As a consequence a reconstruction algorithm able to improve resolution could be beneficial.

Methods

In this study, the hybrid kernelised expectation maximisation (HKEM) is used to improve resolution and contrast and reduce noise, in addition a modified HKEM called frozen HKEM (FHKEM) is investigated to further reduce noise. The iterative part of the FHKEM kernel was frozen at the 72nd sub-iteration. When using ordered subsets algorithms the data is divided in smaller subsets and the smallest algorithm iterative step is called sub-iteration. A NEMA phantom with spherical inserts was used for the optimisation and validation of the algorithm, and data from 5 patients treated with Selective internal radiation therapy were used as proof of clinical relevance of the method.

Results

The results suggest a maximum improvement of 56% for region of interest mean recovery coefficient at fixed coefficient of variation and better identification of the hot volumes in the NEMA phantom. Similar improvements were achieved with patient data, showing 47% mean value improvement over the gold standard used in hospitals.

Conclusions

Such quantitative improvements could facilitate improved dosimetry calculations with SPECT when treating patients with Selective internal radiation therapy, as well as provide a more visible position of the cancerous lesions in the liver.

Role of nanoparticles in transarterial radioembolization with glass microspheres

OBJECTIVE

Transarterial Radioembolization (TARE) with ⁹⁰Y-loaded glass microspheres is a locoregional treatment option for Hepatocellular Carcinoma (HCC). Post-treatment ⁹⁰Y bremsstrahlung imaging using Single-Photon Emission Tomography (SPECT) is currently a gold-standard imaging modality for quantifying the delivered dose. However, the nature of bremsstrahlung photons causes difficulty for dose estimation using SPECT imaging. This work aimed to investigate the possibility of using glass microspheres loaded with ⁹⁰Y and Nanoparticles (NPs) to improve the quantification of delivered doses.

METHODS

The Monte Carlo codes were used to simulate the post-TARE ⁹⁰Y planar imaging. Planar images from bremsstrahlung photons and characteristic X-rays are acquired when 0, 1.2 mol/L, 2.4 mol/L, and 4.8 mol/L of Gold (Au), Hafnium (Hf), and Gadolinium (Gd) NPs are incorporated into the glass microspheres. We evaluated the quality of acquired images by calculating sensitivity and Signal-to-Background Ratio (SBR). Therapeutic effects of NPs were evaluated by calculation of Dose Enhancement Ratio (DER) in tumoral and non-tumoral liver tissues.

RESULTS

The in silico results showed that the sensitivity values of bremsstrahlung and characteristic X-ray planar images increased significantly as the NPs concentration increased in the glass microspheres. The SBR values decreased as the NPs concentration increased for the bremsstrahlung planar images. In contrast, the SBR values increased for the characteristic X-ray planar images when Hf and Gd were incorporated into the glass microspheres. The DER values decreased in the tumoral and non-tumoral liver tissues as the NPs concentration increased. The maximum dose reduction was observed at the NPs concentration of 4.8 mol/L (≈ 7%).

CONCLUSIONS

The incorporation of Au, Hf, and Gd NPs into the glass microspheres improved the quality and quantity of post-TARE planar images. Also, treatment efficiency was decreased significantly at NPs concentration > 4.8 mol/L.

Implications of physics, chemistry and biology for dosimetry calculations using theranostic pairs

Theranostics is an emerging paradigm that combines imaging and therapy in order to personalize patient treatment. In nuclear medicine, this is achieved by using radiopharmaceuticals that target identical molecular targets for both imaging (using emitted gamma rays) and radiopharmaceutical therapy (using emitted beta, alpha or Auger-electron particles) for the treatment of various diseases, such as cancer. If the therapeutic radiopharmaceutical cannot be imaged quantitatively, a “theranostic pair” imaging surrogate can be used to predict the absorbed radiation doses from the therapeutic radiopharmaceutical. However, theranostic dosimetry assumes that the pharmacokinetics and biodistributions of both radiopharmaceuticals in the pair are identical or very similar, an assumption that still requires further validation for many theranostic pairs. In this review, we consider both same-element and different-element theranostic pairs and attempt to determine if factors exist which may cause inaccurate dose extrapolations in theranostic dosimetry, either intrinsic (e.g. chemical differences) or extrinsic (e.g. injecting different amounts of each radiopharmaceutical) to the radiopharmaceuticals. We discuss the basis behind theranostic dosimetry and present common theranostic pairs and their therapeutic applications in oncology. We investigate general factors that could create alterations in the behavior of the radiopharmaceuticals or the quantitative accuracy of imaging them. Finally, we attempt to determine if there is evidence showing some specific pairs as suitable for theranostic dosimetry. We show that there are a variety of intrinsic and extrinsic factors which can significantly alter the behavior among pairs of radiopharmaceuticals, even if they belong to the same chemical element. More research is needed to determine the impact of these factors on theranostic dosimetry estimates and on patient outcomes, and how to correctly account for them.

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.

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.

Post Yttrium-90 Imaging

Transarterial radioembolization with yttrium-90 ( ⁹⁰ Y) is a mainstay for the treatment of liver cancer. Imaging the distribution following delivery is a concept that dates back to the 1960s. As β particles are created during ⁹⁰ Y decay, bremsstrahlung radiation is created as the particles interact with tissues, allowing for imaging with a gamma camera. Inherent qualities of bremsstrahlung radiation make its imaging difficult. SPECT and SPECT/CT can be used but suffer from limitations related to low signal-to-noise bremsstrahlung radiation. However, with optimized imaging protocols, clinically adequate images can still be obtained. A finite but detectable number of positrons are also emitted during ⁹⁰ Y decay, and many studies have demonstrated the ability of commercial PET/CT and PET/MR scanners to image these positrons to understand ⁹⁰ Y distribution and help quantify dose. PET imaging has been proven to be superior to SPECT for quantitative imaging, and therefore will play an important role going forward as we try and better understand dose/response and dose/toxicity relationships to optimize personalized dosimetry. The availability of PET imaging will likely remain the biggest barrier to its use in routine post- ⁹⁰ Y imaging; thus, SPECT/CT imaging with optimized protocols should be sufficient for most posttherapy subjective imaging.

Impacts of acquisition and reconstruction parameters on the absolute technetium quantification of the cadmium-zinc-telluride-based SPECT/CT system: a phantom study

Background

The digital cadmium–zinc–telluride (CZT)-based SPECT system has many advantages, including better spatial and energy resolution. However, the impacts of different acquisition and reconstruction parameters on CZT SPECT quantification might still need to be validated. This study aimed to evaluate the impacts of acquisition parameters (the main energy window and acquisition time per frame) and reconstruction parameters (the number of iterations, subsets in iterative reconstruction, post-filter, and image correction methods) on the technetium quantification of CZT SPECT/CT.

Methods

A phantom (PET NEMA/IEC image quality, USA) was filled with four target-to-background (T/B) ratios (32:1, 16:1, 8:1, and 4:1) of technetium. Mean uptake values (the calculated mean concentrations for spheres) were measured to evaluate the recovery coefficient (RC) changes under different acquisition and reconstruction parameters. The corresponding standard deviations of mean uptake values were also measured to evaluate the quantification error. Image quality was evaluated using the National Electrical Manufacturers Association (NEMA) NU 2–2012 standard.

Results

For all T/B ratios, significant correlations were found between iterations and RCs (r = 0.62–0.96 for 1–35 iterations, r = 0.94–0.99 for 35–90 iterations) as well as between the full width at half maximum (FWHM) of the Gaussian filter and RCs (r = − 0.86 to − 1.00, all P values < 0.05). The regression coefficients of 1–35 iterations were higher than those of 35–90 iterations (0.51–1.60 vs. 0.02–0.19). RCs calculated with AC (attenuation correction) + SC (scatter correction) + RR (resolution recovery correction) combination were more accurate (53.82–106.70%) than those calculated with other combinations (all P values < 0.05). No significant statistical differences (all P values > 0.05) were found between the 15% and 20% energy windows except for the 32:1 T/B ratio (P value = 0.023) or between the 10 s/frame and 120 s/frame acquisition times except for the 4:1 T/B ratio (P value = 0.015) in terms of RCs.

Conclusions

CZT-SPECT/CT of technetium resulted in good quantification accuracy. The favourable acquisition parameters might be a 15% energy window and 40 s/frame of acquisition time. The favourable reconstruction parameters might be 35 iterations, 20 subsets, the AC + SC + RR correction combination, and no filter.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00412-4.

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.

Efficient Monte-Carlo based system modelling for image reconstruction in preclinical pinhole SPECT

The use of multi-pinhole collimation has enabled ultra-high-resolution imaging of SPECT and PET tracers in small animals. Key for obtaining high-quality images is the use of statistical iterative image reconstruction with accurate energy-dependent photon transport modelling through collimator and detector. This can be incorporated in a system matrix that contains the probabilities that a photon emitted from a certain voxel is detected at a specific detector pixel. Here we introduce a fast Monte-Carlo based (FMC-based) matrix generation method for pinhole imaging that is easy to apply to various radionuclides. The method is based on accelerated point source simulations combined with model-based interpolation to straightforwardly change or combine photon energies of the radionuclide of interest. The proposed method was evaluated for a VECTor PET-SPECT system with (i) a HE-UHR-M collimator and (ii) an EXIRAD-3D 3D autoradiography collimator. Both experimental scans with^99mTc,¹¹¹In, and¹²³I, and simulated scans with⁶⁷Ga and⁹⁰Y were performed for evaluation. FMC was compared with two currently used approaches, one based on a set of point source measurements with^99mTc (dubbed traditional method), and the other based on an energy-dependent ray-tracing simulation (ray-tracing method). The reconstruction results show better image quality when using FMC-based matrices than when applying the traditional or ray-tracing matrices in various cases. FMC-based matrices generalise better than the traditional matrices when imaging radionuclides with energies deviating too much from the energy used in the calibration and are computationally more efficient for very-high-resolution imaging than the ray-tracing matrices. In addition, FMC has the advantage of easily combining energies in a single matrix which is relevant when imaging radionuclides with multiple photopeak energies (e.g.⁶⁷Ga and¹¹¹In) or with a continuous energy spectrum (e.g.⁹⁰Y). To conclude, FMC is an efficient, accurate, and versatile tool for creating system matrices for ultra-high-resolution pinhole SPECT.

Stereotactic Inverse Dose Planning After Yttrium-90 Selective Internal Radiation Therapy in Hepatocellular Cancer

PURPOSE

Selective internal radiation therapy (SIRT) is administered to treat tumors of the liver and is generally well tolerated. Although widely adopted for its therapeutic benefits, SIRT is rarely combined with external beam radiation therapy (EBRT) owing to the complexity of the dosimetry resulting from the combination of treatments with distinct radiobiological effects. The purpose of this study was to establish a dosimetric framework for combining SIRT and EBRT using clinical experience derived from representative patients with hepatocellular carcinoma (HCC) who received both therapies.

METHODS AND MATERIALS

Treatments from 10 patients with HCC given EBRT either before or after SIRT were analyzed. The dosimetry framework used here considered differences in the radiobiological effects and fractionation schemes of SIRT versus EBRT, making use of the concepts of biological effective dose (BED) and equivalent dose (EQD). Absorbed dose from SIRT was calculated, converted to BED, and summed with BED from EBRT dose plans. Two of these patients were used in a virtual planning exercise to investigate the feasibility of combining stereotactic body radiation therapy and SIRT.

RESULTS

The combination of EBRT and SIRT in 10 patients with HCC showed no major toxicity. No Child-Pugh scores went above 8 and albumin-bilirubin scores from only 1 patient worsened to grade 3 (> -1.39) from treatment through 3-months follow-up. A framework with radiobiological modeling was developed to manage the combined treatments in terms of their sum BED. The exploratory SIRT plus SABR inverse dose plans for 2 patients, incorporating radiobiologically informed 90Y SIRT dosimetry, achieved dose distributions comparable to SBRT alone.

CONCLUSIONS

Treatment with both EBRT and SIRT can be given safely to patients with HCC. The BED and EQD concepts should be used in combined dosimetry to account for the differing radiobiological effects of EBRT and SIRT. Inverse dose planning of EBRT after SIRT could provide improved dose distributions and flexibility to the clinical workflow. Further research into combination therapy is needed through prospective trials.

Evaluation of qualitative and quantitative data of Y-90 imaging in SPECT/CT and PET/CT phantom studies

Introduction

We aimed to assess the feasibility of SPECT and PET Y-90 imaging, and to compare these modalities by visualizing hot and cold foci in phantoms for varying isotope concentrations.

Materials and methods

The data was acquired from the Jaszczak and NEMA phantoms. In the Jaszczak phantom Y-90 concentrations of 0.1 MBq/ml and 0.2 MBq/ml were used, while higher concentrations, up to 1.0 MBq/ml, were simulated by acquisition time extension with respect to the standard clinical protocol of 30 sec/projection for SPECT and 30 min/bed position for PET imaging. For NEMA phantom, the hot foci had concentrations of about 4 MB/ml and the background 0.1 or 0.0 MBq/ml. All of the acquired data was analysed both qualitatively and quantitatively. Qualitative assessment was conducted by six observers asked to identify the number of visible cold or hot foci. Inter-observer agreement was assessed. Quantitative analysis included calculations of contrast and contrast-to-noise ratio (CNR), and comparisons with the qualitative results.

Results

For SPECT data up to two cold foci were discernible, while for PET four foci were visible. We have shown that CNR (with Rose criterion) is a good measure of foci visibility for both modalities. We also found good concordance of qualitative results for the Jaszczak phantom studies between the observers (corresponding Krippendorf’s alpha coefficients of 0.76 to 0.84).

In the NEMA phantom without background activity all foci were visible in SPECT/CT images. With isotope in the background, 5 of 6 spheres were discernible (CNR of 3.0 for the smallest foci). For PET studies all hot spheres were visible, regardless of the background activity.

Conclusions

PET Y-90 imaging provided better results than Bremsstrahlung based SPECT imaging. This indicates that PET/CT might become the method of choice in Y-90 post radioembolization imaging for visualisation of both necrotic and hot lesions in the liver.

Diagnostic Performance of Theranostic Radionuclides Used in Transarterial Radioembolization for Liver Cancer

Primary liver tumor with hepatocellular carcinoma accounting for 75-80% of all such tumors, is one of the global leading causes of cancer-related death, especially in cirrhotic patients. Liver tumors are highly hypervascularized via the hepatic artery, while normal liver tissues are mainly supplied by the portal vein; consequently, intra-arterially delivered treatment, which includes transarterial chemoembolization (TACE) and transarterial radioembolization (TARE), is deemed as a palliative treatment. With the development of nuclear technology and radiochemistry, TARE has become an alternative for patients with hepatic cancer, especially for patients who failed other therapies, or for patients who need tumor downstaging treatment. In practice, some radionuclides have suitable physicochemical characteristics to act as radioactive embolism agents. Among them, 90Y emits β rays only and is suitable for bremsstrahlung single photon emission computed tomography (BS SPECT) and positron emission tomography (PET); meanwhile, some others, such as 131I, 153Sm, 166Ho, 177Lu, 186Re, and 188Re, emit both β and γ rays, enabling embolism beads to play a role in both therapy and single photon emission computed tomography (SPECT) imaging. During TARE, concomitant imaging provide additive diagnostic information and help to guide the course of liver cancer treatment. Therefore, we review the theranostic radionuclides that have been used or could potentially be used in TARE for liver cancer and focus on the clinical benefits of diagnostic applications, including real-time monitoring of embolism beads, evaluating irradiation dose, predicting therapy effects, and corresponding adjustments to TARE.

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.

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.

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.

Collimator and energy window optimization for YTTRIUM-90 bremsstrahlung SPECT imaging

The optimal collimator and energy window for Yttrium-90 bremsstrahlung SPECT imaging was investigated in the study. Yttrium-90 images were acquired with a dual-head gamma camera, equipped with parallel hole collimators and ⁹⁰Y vial for different energy windows ranging from 56 to 232 keV. Image quality parameters (sensitivity, %FOV, and S/B) were examined for the energy window and collimator combinations. It is concluded that the optimal SPECT imaging was achieved using FBP Method with a HEGP collimator and the energy window of 90-110 keV.

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.

Algorithms and Analyses for Joint Spectral Image Reconstruction in Y-90 Bremsstrahlung SPECT

Quantitative yttrium-90 (Y-90) SPECT imaging is challenging due to the nature of Y-90, an almost pure beta emitter that is associated with a continuous spectrum of bremsstrahlung photons that have a relatively low yield. This paper proposes joint spectral reconstruction (JSR), a novel bremsstrahlung SPECT reconstruction method that uses multiple narrow acquisition windows with accurate multi-band forward modeling to cover a wide range of the energy spectrum. Theoretical analyses using Fisher information and Monte-Carlo (MC) simulation with a digital phantom show that the proposed JSR model with multiple acquisition windows has better performance in terms of covariance (precision) than previous methods using multi-band forward modeling with a single acquisition window, or using a single-band forward modeling with a single acquisition window. We also propose an energy-window subset (ES) algorithm for JSR to achieve fast empirical convergence and maximum-likelihood based initialization for all reconstruction methods to improve quantification accuracy in early iterations. For both MC simulation with a digital phantom and experimental study with a physical multi-sphere phantom, our proposed JSR-ES, a fast algorithm for JSR with ES, yielded higher recovery coefficients (RCs) on hot spheres over all iterations and sphere sizes than all the other evaluated methods, due to fast empirical convergence. In experimental study, for the smallest hot sphere (diameter 1.6cm), at the 20th iteration the increase in RCs with JSR-ES was 66 and 31% compared with single wide and narrow band forward models, respectively. JSR-ES also yielded lower residual count error (RCE) on a cold sphere over all iterations than other methods for MC simulation with known scatter, but led to greater RCE compared with single narrow band forward model at higher iterations for experimental study when using estimated scatter.

Towards standardization of absolute SPECT/CT quantification: a multi-center and multi-vendor phantom study

Abstract

Absolute quantification of radiotracer distribution using SPECT/CT imaging is of great importance for dosimetry aimed at personalized radionuclide precision treatment. However, its accuracy depends on many factors. Using phantom measurements, this multi-vendor and multi-center study evaluates the quantitative accuracy and inter-system variability of various SPECT/CT systems as well as the effect of patient size, processing software and reconstruction algorithms on recovery coefficients (RC).

Methods

Five SPECT/CT systems were included: Discovery™ NM/CT 670 Pro (GE Healthcare), Precedence™ 6 (Philips Healthcare), Symbia Intevo™, and Symbia™ T16 (twice) (Siemens Healthineers). Three phantoms were used based on the NEMA IEC body phantom without lung insert simulating body mass indexes (BMI) of 25, 28, and 47 kg/m². Six spheres (0.5–26.5 mL) and background were filled with 0.1 and 0.01 MBq/mL ^99mTc-pertechnetate, respectively. Volumes of interest (VOI) of spheres were obtained by a region growing technique using a 50% threshold of the maximum voxel value corrected for background activity. RC, defined as imaged activity concentration divided by actual activity concentration, were determined for maximum (RC_max) and mean voxel value (RC_mean) in the VOI for each sphere diameter. Inter-system variability was expressed as median absolute deviation (MAD) of RC. Acquisition settings were standardized. Images were reconstructed using vendor-specific 3D iterative reconstruction algorithms with institute-specific settings used in clinical practice and processed using a standardized, in-house developed processing tool based on the SimpleITK framework. Additionally, all data were reconstructed with a vendor-neutral reconstruction algorithm (Hybrid Recon™; Hermes Medical Solutions).

Results

RC decreased with decreasing sphere diameter for each system. Inter-system variability (MAD) was 16 and 17% for RC_mean and RC_max, respectively. Standardized reconstruction decreased this variability to 4 and 5%. High BMI hampers quantification of small lesions (< 10 ml).

Conclusion

Absolute SPECT quantification in a multi-center and multi-vendor setting is feasible, especially when reconstruction protocols are standardized, paving the way for a standard for absolute quantitative SPECT.

90 Y SPECT/CT quantitative study and comparison of uptake with pretreatment 99 m Tc-MAA SPECT/CT in radiomicrosphere therapy

INTRODUCTION

Yttrium-90 (90 Y) microsphere post-treatment imaging reflects the true distribution characteristics of microspheres in the tumor and liver compartments. However, due to its decay spectra profile lacking a pronounced photopeak, the bremsstrahlung imaging for 90 Y has inherent limitations. The absorbed dose calculations for 90 Y microspheres radiomicrosphere therapy (RMT) sustain a limitation due to the poor quality of 90 Y imaging. The aim of this study was to develop quantitative methods to improve the post-treatment 90 Y bremsstrahlung single photon emission tomography (SPECT)/computed tomography (CT) image analysis for dosimetric purposes and to perform a quantitative comparison with the 99m Tc-MAA SPECT/CT images, which is used for theranostics purposes for liver and tumor dosimetry.

METHODS

Pre and post-treatment SPECT/CT data of patients who underwent RMT for primary or metastatic liver cancer were acquired. A Jasczak phantom with eight spherical inserts of various sizes was used to obtain optimal iteration number for the contrast recovery algorithm for improving 90 Y bremsstrahlung SPECT/CT images. Comparison of uptake on 99m Tc-MAA and 90 Y microsphere SPECT/CT images was assessed using tumor to healthy liver ratios (TLRs). The voxel dosimetry technique was used to estimate absorbed doses. Absorbed doses within the tumor and healthy part of the liver were also investigated for correlation with administered activity.

RESULTS

Improvement in CNR and contrast recovery coefficients on patient and phantom 90 Y bremsstrahlung SPECT/CT images respectively were achieved. The 99m Tc-MAA and 90 Y microspheres SPECT/CT images showed significant uptake correlation (r = 0.9, P = 0.05) with mean TLR of 9.4 ± 9.2 and 5.0 ± 2.2, respectively. The correlation between the administered activity and tumor absorbed dose was weak (r = 0.5, P > 0.05), however, healthy liver absorbed dose increased with administered activity (r = 0.8, P = 0.0).

CONCLUSIONS

This study demonstrated correlation in mean TLR between 99m Tc-MAA and 90 Y microsphere SPECT/CT.

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.

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.

Improved quantitative 90 Y bremsstrahlung SPECT/CT reconstruction with Monte Carlo scatter modeling

PURPOSE

In ⁹⁰ Y microsphere radioembolization (RE), accurate post-therapy imaging-based dosimetry is important for establishing absorbed dose versus outcome relationships for developing future treatment planning strategies. Additionally, accurately assessing microsphere distributions is important because of concerns for unexpected activity deposition outside the liver. Quantitative ⁹⁰ Y imaging by either SPECT or PET is challenging. In ⁹⁰ Y SPECT model based methods are necessary for scatter correction because energy window-based methods are not feasible with the continuous bremsstrahlung energy spectrum. The objective of this work was to implement and evaluate a scatter estimation method for accurate ⁹⁰ Y bremsstrahlung SPECT/CT imaging.

METHODS

Since a fully Monte Carlo (MC) approach to ⁹⁰ Y SPECT reconstruction is computationally very demanding, in the present study the scatter estimate generated by a MC simulator was combined with an analytical projector in the 3D OS-EM reconstruction model. A single window (105 to 195-keV) was used for both the acquisition and the projector modeling. A liver/lung torso phantom with intrahepatic lesions and low-uptake extrahepatic objects was imaged to evaluate SPECT/CT reconstruction without and with scatter correction. Clinical application was demonstrated by applying the reconstruction approach to five patients treated with RE to determine lesion and normal liver activity concentrations using a (liver) relative calibration.

RESULTS

There was convergence of the scatter estimate after just two updates, greatly reducing computational requirements. In the phantom study, compared with reconstruction without scatter correction, with MC scatter modeling there was substantial improvement in activity recovery in intrahepatic lesions (from > 55% to > 86%), normal liver (from 113% to 104%), and lungs (from 227% to 104%) with only a small degradation in noise (13% vs. 17%). Similarly, with scatter modeling contrast improved substantially both visually and in terms of a detectability index, which was especially relevant for the low uptake extrahepatic objects. The trends observed for the phantom were also seen in the patient studies where lesion activity concentrations and lesion-to-liver concentration ratios were lower for SPECT without scatter correction compared with reconstruction with just two MC scatter updates: in eleven lesions the mean uptake was 4.9 vs. 7.1 MBq/mL (P = 0.0547), the mean normal liver uptake was 1.6 vs. 1.5 MBq/mL (P = 0.056) and the mean lesion-to-liver uptake ratio was 2.7 vs. 4.3 (P = 0.0402) for reconstruction without and with scatter correction respectively.

CONCLUSIONS

Quantitative accuracy of ⁹⁰ Y bremsstrahlung imaging can be substantially improved with MC scatter modeling without significant degradation in image noise or intensive computational requirements.

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.