Comparison of quantitative Y-90 SPECT and non-time-of-flight PET imaging in post-therapy radioembolization of liver cancer

Using Voxel-Based Dosimetry to Evaluate Sphere Concentration and Tumor Dose in Hepatocellular Carcinoma Treated with Yttrium-90 Radiation Segmentectomy with Glass Microspheres

PURPOSE

To utilize voxel-based dosimetry following radiation segmentectomy (RS) to understand microsphere distribution and validate current literature regarding radiologic and pathologic outcomes.

MATERIALS AND METHODS

A retrospective, single-center analysis of patients with solitary hepatocellular carcinoma (N = 56) treated with yttrium-90 (90Y) RS with glass microspheres (Therasphere; Boston Scientific, Marlborough, Massachusetts) from 2020 to 2022 was performed. Posttreatment voxel-based dosimetry was evaluated using Mirada DBx Build 1.2.0 Simplicit90Y software (Boston Scientific) and utilized to calculate sphere concentration to tumor as well as D70 (minimum dose to 70% total tumor volume), D90, and D99. Time to progression (TTP), treatment response, and adverse events were studied.

RESULTS

Fifty-six solitary tumors were analyzed with a median tumor diameter of 3.4 cm (range, 1.2-6.8 cm) and median tumor absorbed dose of 732 Gy (range, 252-1,776 Gy). Median sphere activity (SA) at the time of delivery was 1,446 Bq (range, 417-2,621 Bq). Median tumor sphere concentration was 12,868 spheres/mL (range, 2,655-37,183 spheres/mL). Sphere concentration into tumor and normal tissue was inversely correlated with perfused treatment volume (R2 = 0.21 and R2 = 0.39, respectively). Of the 51 tumors with posttreatment imaging, objective response was noted in 49 patients (96%) and complete response in 42 patients (82%). The median TTP was not reached with a 2-year progression rate of 11%. Fifteen patients underwent liver transplant. Median tumor necrosis was 99% (range, 80%-100%). Lower tumor volumes and higher D99 were associated with complete pathologic necrosis (P < .001 and P = .022, respectively).

CONCLUSIONS

Voxel-based dosimetry following 90Y radioembolization can be utilized to account for sphere deposition and distribution into tumor. Ablative RS with high SA yields durable radiologic and pathologic outcomes.

Enhancing voxel-based dosimetry accuracy with an unsupervised deep learning approach for hybrid medical image registration

BACKGROUND

Deformable registration is required to generate a time-integrated activity (TIA) map which is essential for voxel-based dosimetry. The conventional iterative registration algorithm using anatomical images (e.g., computed tomography (CT)) could result in registration errors in functional images (e.g., single photon emission computed tomography (SPECT) or positron emission tomography (PET)). Various deep learning-based registration tools have been proposed, but studies specifically focused on the registration of serial hybrid images were not found.

PURPOSE

In this study, we introduce CoRX-NET, a novel unsupervised deep learning network designed for deformable registration of hybrid medical images. The CoRX-NET structure is based on the Swin-transformer (ST), allowing for the representation of complex spatial connections in images. Its self-attention mechanism aids in the effective exchange and integration of information across diverse image regions. To augment the amalgamation of SPECT and CT features, cross-stitch layers have been integrated into the network.

METHODS

Two different 177 Lu DOTATATE SPECT/CT datasets were acquired at different medical centers. 22 sets from Seoul National University and 14 sets from Sunway Medical Centre are used for training/internal validation and external validation respectively. The CoRX-NET architecture builds upon the ST, enabling the modeling of intricate spatial relationships within images. To further enhance the fusion of SPECT and CT features, cross-stitch layers have been incorporated within the network. The network takes a pair of SPECT/CT images (e.g., fixed and moving images) and generates a deformed SPECT/CT image. The performance of the network was compared with Elastix and TransMorph using L1 loss and structural similarity index measure (SSIM) of CT, SSIM of normalized SPECT, and local normalized cross correlation (LNCC) of SPECT as metrics. The voxel-wise root mean square errors (RMSE) of TIA were compared among the different methods.

RESULTS

The ablation study revealed that cross-stitch layers improved SPECT/CT registration performance. The cross-stitch layers notably enhance SSIM (internal validation: 0.9614 vs. 0.9653, external validation: 0.9159 vs. 0.9189) and LNCC of normalized SPECT images (internal validation: 0.7512 vs. 0.7670, external validation: 0.8027 vs. 0.8027). CoRX-NET with the cross-stitch layer achieved superior performance metrics compared to Elastix and TransMorph, except for CT SSIM in the external dataset. When qualitatively analyzed for both internal and external validation cases, CoRX-NET consistently demonstrated superior SPECT registration results. In addition, CoRX-NET accomplished SPECT/CT image registration in less than 6 s, whereas Elastix required approximately 50 s using the same PC's CPU. When employing CoRX-NET, it was observed that the voxel-wise RMSE values for TIA were approximately 27% lower for the kidney and 33% lower for the tumor, compared to when Elastix was used.

CONCLUSION

This study represents a major advancement in achieving precise SPECT/CT registration using an unsupervised deep learning network. It outperforms conventional methods like Elastix and TransMorph, reducing uncertainties in TIA maps for more accurate dose assessments.

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.

Y-90 PET/MR imaging optimization with a Bayesian penalized likelihood reconstruction algorithm

Positron Emission Tomography (PET) imaging after90 Y liver radioembolization is used for both lesion identification and dosimetry. Bayesian penalized likelihood (BPL) reconstruction algorithms are an alternative to ordered subset expectation maximization (OSEM) with improved image quality and lesion detectability. The investigation of optimal parameters for90 Y image reconstruction of Q.Clear, a commercial BPL algorithm developed by General Electric (GE), in PET/MR is a field of interest and the subject of this study. The NEMA phantom was filled at an 8:1 sphere-to-background ratio. Acquisitions were performed on a PET/MR scanner for clinically relevant activities between 0.7 and 3.3 MBq/ml. Reconstructions with Q.Clear were performed varying the β penalty parameter between 20 and 6000, the acquisition time between 5 and 20 min and pixel size between 1.56 and 4.69 mm. OSEM reconstructions of 28 subsets with 2 and 4 iterations with and without Time-of-Flight (TOF) were compared to Q.Clear with β = 4000. Recovery coefficients (RC), their coefficient of variation (COV), background variability (BV), contrast-to-noise ratio (CNR) and residual activity in the cold insert were evaluated. Increasing β parameter lowered RC, COV and BV, while CNR was maximized at β = 4000; further increase resulted in oversmoothing. For quantification purposes, β = 1000-2000 could be more appropriate. Longer acquisition times resulted in larger CNR due to reduced image noise. Q.Clear reconstructions led to higher CNR than OSEM. A β of 4000 was obtained for optimal image quality, although lower values could be considered for quantification purposes. An optimal acquisition time of 15 min was proposed considering its clinical use.

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.

Optimization of Y-90 Radioembolization Imaging for Post-Treatment Dosimetry on a Long Axial Field-of-View PET/CT Scanner

BACKGROUND

PET imaging after yttrium-90 (Y-90) radioembolization is challenging because of the low positron fraction of Y-90 (32 × 10-6). The resulting low number of events can be compensated by the high sensitivity of long axial field-of-view (LAFOV) PET/CT scanners. Nevertheless, the reduced event statistics require optimization of the imaging protocol to achieve high image quality (IQ) and quantification accuracy sufficient for post-treatment dosimetry.

METHODS

Two phantoms (NEMA IEC and AbdoMan phantoms, mimicking human liver) filled with Y-90 and a 4:1 sphere (tumor)-to-background ratio were scanned for 24 h with the Biograph Vision Quadra (Siemens Healthineers). Eight patients were scanned after Y-90 radioembolization (1.3-4.7 GBq) using the optimized protocol (obtained by phantom studies). The IQ, contrast recovery coefficients (CRCs) and noise were evaluated for their limited and full acceptance angles, different rebinned scan durations, numbers of iterations and post-reconstruction filters. The s-value-based absorbed doses were calculated to assess their suitability for dosimetry.

RESULTS

The phantom studies demonstrate that two iterations, five subsets and a 4 mm Gaussian filter provide a reasonable compromise between a high CRC and low noise. For a 20 min scan duration, an adequate CRC of 56% (vs. 24 h: 62%, 20 mm sphere) was obtained, and the noise was reduced by a factor of 1.4, from 40% to 29%, using the full acceptance angle. The patient scan results were consistent with those from the phantom studies, and the impacts on the absorbed doses were negligible for all of the studied parameter sets, as the maximum percentage difference was -3.89%.

CONCLUSIONS

With 2i5s, a 4 mm filter and a scan duration of 20 min, IQ and quantification accuracy that are suitable for post-treatment dosimetry of Y-90 radioembolization can be achieved.

Predictive value of 99mTc-MAA-based dosimetry in personalized 90Y-SIRT planning for liver malignancies

BACKGROUND

Selective internal radiation therapy with ⁹⁰Y radioembolization aims to selectively irradiate liver tumours by administering radioactive microspheres under the theragnostic assumption that the pre-therapy injection of ^99mTc labelled macroaggregated albumin (^99mTc-MAA) provides an estimation of the ⁹⁰Y microspheres biodistribution, which is not always the case. Due to the growing interest in theragnostic dosimetry for personalized radionuclide therapy, a robust relationship between the delivered and pre-treatment radiation absorbed doses is required. In this work, we aim to investigate the predictive value of absorbed dose metrics calculated from ^99mTc-MAA (simulation) compared to those obtained from ⁹⁰Y post-therapy SPECT/CT.

RESULTS

A total of 79 patients were analysed. Pre- and post-therapy 3D-voxel dosimetry was calculated on ^99mTc-MAA and ⁹⁰Y SPECT/CT, respectively, based on Local Deposition Method. Mean absorbed dose, tumour-to-normal ratio, and absorbed dose distribution in terms of dose-volume histogram (DVH) metrics were obtained and compared for each volume of interest (VOI). Mann-Whitney U-test and Pearson's correlation coefficient were used to assess the correlation between both methods. The effect of the tumoral liver volume on the absorbed dose metrics was also investigated. Strong correlation was found between simulation and therapy mean absorbed doses for all VOIs, although simulation tended to overestimate tumour absorbed doses by 26%. DVH metrics showed good correlation too, but significant differences were found for several metrics, mostly on non-tumoral liver. It was observed that the tumoral liver volume does not significantly affect the differences between simulation and therapy absorbed dose metrics.

CONCLUSION

This study supports the strong correlation between absorbed dose metrics from simulation and therapy dosimetry based on ⁹⁰Y SPECT/CT, highlighting the predictive ability of ^99mTc-MAA, not only in terms of mean absorbed dose but also of the dose distribution.

Triple modality image reconstruction of PET data using SPECT, PET, CT information increases lesion uptake in images of patients treated with radioembolization with [Formula: see text] micro-spheres

PURPOSE

Nuclear medicine imaging modalities like computed tomography (CT), single photon emission CT (SPECT) and positron emission tomography (PET) are employed in the field of theranostics to estimate and plan the dose delivered to tumors and the surrounding tissues and to monitor the effect of the therapy. However, therapeutic radionuclides often provide poor images, which translate to inaccurate treatment planning and inadequate monitoring images. Multimodality information can be exploited in the reconstruction to enhance image quality. Triple modality PET/SPECT/CT scanners are particularly useful in this context due to the easier registration process between images. In this study, we propose to include PET, SPECT and CT information in the reconstruction of PET data. The method is applied to Yttrium-90 ([Formula: see text]Y) data.

METHODS

Data from a NEMA phantom filled with [Formula: see text]Y were used for validation. PET, SPECT and CT data from 10 patients treated with Selective Internal Radiation Therapy (SIRT) were used. Different combinations of prior images using the Hybrid kernelized expectation maximization were investigated in terms of VOI activity and noise suppression.

RESULTS

Our results show that triple modality PET reconstruction provides significantly higher uptake when compared to the method used as standard in the hospital and OSEM. In particular, using CT-guided SPECT images, as guiding information in the PET reconstruction significantly increases uptake quantification on tumoral lesions.

CONCLUSION

This work proposes the first triple modality reconstruction method and demonstrates up to 69% lesion uptake increase over standard methods with SIRT [Formula: see text]Y patient data. Promising results are expected for other radionuclide combination used in theranostic applications using PET and SPECT.

Comparison of PMT-based TF64 and SiPM-based Vereos PET/CT systems for 90 Y imaging and dosimetry optimization: A quantitative study

BACKGROUND

Selective internal radiotherapy based on transarterial radioembolization (TARE) with yttrium-90 (⁹⁰ Y) microspheres is an established treatment for primary or metastatic liver disease.

PURPOSE

The objective of this work is to optimize the dosimetry of patients treated with ⁹⁰ Y TARE, using positron emission tomography (PET) images.

METHODS

The NEMA 2012 PET phantom was filled with nearly 3.9 GBq of ⁹⁰ Y activity and acquired at days 0, 3, 5, 7, and 9 on a classic time-of-flight PET/computed tomography (CT) scanner (Philips TF64) and on a silicon photomultiplier (SiPM)-based PET/CT scanner (Philips Vereos). Acquisitions were carried on following the guidelines proposed in a previously published multicentric trial and images were reconstructed by varying and combining the available parameters. Comparisons were performed to identify the best set(s) of parameters leading to the most accurate ⁹⁰ Y-PET image(s), in terms of activity distribution. Then, for both scanners, the best images were analyzed with Simplicit⁹⁰ Y, a personalized dosimetry software using multicompartmental Medical Internal Radiation Dose model. The comparison between measured and true doses allowed to identify the image granting the most consistent dose estimations and, therefore, to designate the set of parameters to be applied on patients' data for the reconstruction of optimized clinical images. Posttreatment dosimetry of four patients was then realized with Simplicit⁹⁰ Y using optimized imaging datasets.

RESULTS

Based on activity distribution comparisons and dose estimations over phantom and patients data, the SiPM-based PET/CT system appeared more suitable than the photomultiplier tube-based TF64 for ⁹⁰ Y-PET imaging. With the SiPM-based PET/CT system, reconstructed images with a 2-mm voxel size combined with the application of the point spread function correction led to the most accurate results for quantitative ⁹⁰ Y measures.

CONCLUSIONS

For the SiPM-based PET/CT scanner, an optimized set of reconstruction parameters has been identified and applied on patients' data in order to generate the most accurate image to be used for an improved personalized ⁹⁰ Y-PET dosimetry, ensuring a reliable evaluation of the delivered doses.

Beyond the MAA-Y90 Paradigm: The Evolution of Radioembolization Dosimetry Approaches and Scout Particles

Radioembolization is a well-established treatment for primary and metastatic liver cancer. There is increasing interest in personalized treatment planning supported by dosimetry, as it provides an opportunity to optimize dose delivery to tumor and minimize nontarget deposition, which demonstrably increases the efficacy and safety of this therapy. However, the optimal dosimetry procedure in the radioembolization setting is still evolving; existing data are limited as few trials have prospectively tailored dose based on personalized planning and predominantly semi-empirical methods are used for dose calculation. Since the pretreatment or "scout" procedure forms the basis of dosimetry calculations, an accurate and reliable technique is essential. ^99m Tc-MAA SPECT constitutes the current accepted standard for pretreatment imaging; however, inconsistent patterns in published data raise the question whether this is the optimal agent. Alternative particles are now being introduced to the market, and early indications suggest use of an identical scout and treatment particle may be superior to the current standard. This review will undertake an evaluation of the increasingly refined dosimetric methods driving radioembolization practices, and a horizon scanning exercise identifying alternative scout particle solutions. Together these constitute a compelling vision for future treatment planning methods that prioritize individualized care.

Evaluation of different CT maps for attenuation correction and segmentation in static 99m Tc-MAA SPECT/CT for 90 Y radioembolization treatment planning: A simulation study

PURPOSE

Conventional ^99m Tc-macroaggregated albumin (^99m Tc-MAA) planar scintigraphy overestimates lung shunt fraction (LSF) compared to SPECT/CT. However, the respiratory motion artifact due to the temporal mismatch between static SPECT and helical CT (HCT) may compromise the SPECT quantitation accuracy by incorrect attenuation correction (AC) and volume-of-interest (VOI) segmentation. This study aims to evaluate AC and VOI segmentation effects systematically and to propose a CT map for LSF and tumor-to-normal liver ratio (TNR) estimation in static ^99m Tc-MAA SPECT/CT.

METHODS

The 4D XCAT phantom was used to simulate a phantom population of 120 phantoms, modeling 10 different anatomical variations, nine TNRs (2-13.2), nine tumor sizes (2-6.7 cm diameter), eight tumor locations, three axial motion amplitudes of 1, 1.5, and 2 (cm), and four LSFs of 5%, 10%, 15%, and 20%. An analytical projector for low-energy high-resolution parallel-hole collimator was used to simulate 60 noisy projections over 360°, modeling attenuation and geometric collimator-detector response (GCDR). AC and VOI mismatch effects were investigated independently and together, using cine average CT (CACT), HCT at end-inspiration (HCT-IN), mid-respiration (HCT-MID), and end-expiration (HCT-EX) respectively as attenuation and segmentation maps. SPECT images without motion, AC, and VOI errors were also generated as reference. LSF and TNR errors were measured as compared to the ground truth.

RESULTS

HCT-MID has slightly better performance for AC effect compared with other CT maps in LSF and TNR estimation, while HCT-EX and HCT-MID perform better for VOI effect. For a respiratory motion amplitude of 1.5 cm and a LSF of 5%, the LSF errors are 19.56 ± 4.58%, -6.79 ± 1.74%, 77.29 ± 14.74%, and 111.25 ± 18.29% corresponding to HCT-MID, HCT-EX, HCT-IN, and CACT in static SPECT. The TNR errors are -12.38 ± 6.42%, -20.55 ± 11.25%, -20.89 ± 9.98%, and -22.89 ± 14.38% respectively. HCT-MID has the best performance for LSF estimation for LSF > 10% and TNR estimation, followed by HCT-EX, HCT-IN, and CACT.

CONCLUSIONS

The HCT-MID is recommended for AC and segmentation to alleviate respiratory artifacts and improve quantitation accuracy in ⁹⁰ Y radioembolization treatment planning. HCT-EX would also be a recommended choice if HCT-MID is not available.

International recommendations for personalised selective internal radiation therapy of primary and metastatic liver diseases with yttrium-90 resin microspheres

PURPOSE

A multidisciplinary expert panel convened to formulate state-of-the-art recommendations for optimisation of selective internal radiation therapy (SIRT) with yttrium-90 (90Y)-resin microspheres.

METHODS

A steering committee of 23 international experts representing all participating specialties formulated recommendations for SIRT with 90Y-resin microspheres activity prescription and post-treatment dosimetry, based on literature searches and the responses to a 61-question survey that was completed by 43 leading experts (including the steering committee members). The survey was validated by the steering committee and completed anonymously. In a face-to-face meeting, the results of the survey were presented and discussed. Recommendations were derived and level of agreement defined (strong agreement ≥ 80%, moderate agreement 50%-79%, no agreement ≤ 49%).

RESULTS

Forty-seven recommendations were established, including guidance such as a multidisciplinary team should define treatment strategy and therapeutic intent (strong agreement); 3D imaging with CT and an angiography with cone-beam-CT, if available, and 99mTc-MAA SPECT/CT are recommended for extrahepatic/intrahepatic deposition assessment, treatment field definition and calculation of the 90Y-resin microspheres activity needed (moderate/strong agreement). A personalised approach, using dosimetry (partition model and/or voxel-based) is recommended for activity prescription, when either whole liver or selective, non-ablative or ablative SIRT is planned (strong agreement). A mean absorbed dose to non-tumoural liver of 40 Gy or less is considered safe (strong agreement). A minimum mean target-absorbed dose to tumour of 100-120 Gy is recommended for hepatocellular carcinoma, liver metastatic colorectal cancer and cholangiocarcinoma (moderate/strong agreement). Post-SIRT imaging for treatment verification with 90Y-PET/CT is recommended (strong agreement). Post-SIRT dosimetry is also recommended (strong agreement).

CONCLUSION

Practitioners are encouraged to work towards adoption of these recommendations.

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.

Tumor-to-Normal Ratio Relationship between Planning Technetium-99 Macroaggregated Albumin and Posttherapy Yttrium-90 Bremsstrahlung SPECT/CT

PURPOSE

To quantify the relationship of the tumor-to-normal ratio (TNR) attained from the technetium-99m macroaggregated albumin (MAA) and posttreatment yttrium-90 bremsstrahlung (Y90-Brem) single-photon emission computerized tomography (SPECT)/computer tomography (CT) studies in patients with hepatocellular carcinoma (HCC) treated with glass microspheres.

MATERIALS AND METHODS

Retrospectively, a total of 190 consecutive patients with HCC who underwent 204 MAA and Y90-Brem SPECT/CT for glass microsphere Y90 radiation segmentectomy (Y90-RS) or lobar treatment (Y90-RLT) between 2013 and 2018 were included. Semi-automated regions-of-interests were drawn around the targeted tumor and nontumoral liver tissue on the SPECT/CT studies. TNR values from MAA and Y90-Brem SPECT/CT were compared using paired t-tests, Pearson correlation, and median with interquartile ranges (IQR).

RESULTS

The mean TNR for MAA and Y90-Brem SPECT/CT was 2.96 ± 1.86 (median, 2.64; IQR, 2.50) and 2.29 ± 1.10 (median, 2.06; IQR, 1.05), respectively (P < .0001). The mean Y90-RLT TNR was 2.88 ± 1.67 (median, 2.59; IQR, 0.83) and 2.17 ± 0.89 (median, 1.98; IQR, 0.81) for MAA and Y90-Brem SPECT/CT, respectively (P < .0001). The mean Y90-RS TNR was 3.02 ± 2.01 (median, 2.87; IQR, 3.01) and 2.39 ± 1.25 (median, 2.11; IQR, 1.28) for MAA and Y90-Brem SPECT/CT, respectively (P = .0003). TNR attained from MAA and Y90 SPECT/CT studies showed a moderate correlation in a positive linear fashion for the overall (r = 0.54; P < .001), Y90-RLT (r = 0.66, P < .001), and Y90-RS cohorts (r = 0.48, P < .001).

CONCLUSIONS

The TNR attained from Y90-Brem SPECT/CT is often underestimated, positively correlated, and less variable than that attained from MAA SPECT/CT.

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.

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.

Radiopharmaceutical therapy in cancer: clinical advances and challenges

Radiopharmaceutical therapy (RPT) is emerging as a safe and effective targeted approach to treating many types of cancer. In RPT, radiation is systemically or locally delivered using pharmaceuticals that either bind preferentially to cancer cells or accumulate by physiological mechanisms. Almost all radionuclides used in RPT emit photons that can be imaged, enabling non-invasive visualization of the biodistribution of the therapeutic agent. Compared with almost all other systemic cancer treatment options, RPT has shown efficacy with minimal toxicity. With the recent FDA approval of several RPT agents, the remarkable potential of this treatment is now being recognized. This Review covers the fundamental properties, clinical development and associated challenges of RPT.

Current Status of Radiopharmaceutical Therapy

In radiopharmaceutical therapy (RPT), a radionuclide is systemically or locally delivered with the goal of targeting and delivering radiation to cancer cells while minimizing radiation exposure to untargeted cells. Examples of current RPTs include thyroid ablation with the administration of ¹³¹I, treatment of liver cancer with ⁹⁰Y microspheres, the treatment of bony metastases with ²²³Ra, and the treatment of neuroendocrine tumors with ¹⁷⁷Lu-DOTATATE. New RPTs are being developed where radionuclides are incorporated into systemic targeted therapies. To assure that RPT is appropriately implemented, advances in targeting need to be matched with advances in quantitative imaging and dosimetry methods. Currently, radiopharmaceutical therapy is administered by intravenous or locoregional injection, and the treatment planning has typically been implemented like chemotherapy, where the activity administered is either fixed or based on a patient's body weight or body surface area. RPT pharmacokinetics are measurable by quantitative imaging and are known to vary across patients, both in tumors and normal tissues. Therefore, fixed or weight-based activity prescriptions are not currently optimized to deliver a cytotoxic dose to targets while remaining within the tolerance dose of organs at risk. Methods that provide dose estimates to individual patients rather than to reference geometries are needed to assess and adjust the injected RPT dose. Accurate doses to targets and organs at risk will benefit the individual patients and decrease uncertainties in clinical trials. Imaging can be used to measure activity distribution in vivo, and this information can be used to determine patient-specific treatment plans where the dose to the targets and organs at risk can be calculated. The development and adoption of imaging-based dosimetry methods is particularly beneficial in early clinical trials. In this work we discuss dosimetric accuracy needs in modern radiation oncology, uncertainties in the dosimetry in RPT, and best approaches for imaging and dosimetry of internal radionuclide therapy.

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.

Comparison of the Biograph Vision and Biograph mCT for quantitative 90Y PET/CT imaging for radioembolisation

BACKGROUND

New digital PET scanners with improved time of flight timing and extended axial field of view such as the Siemens Biograph Vision have come on the market and are expected to replace current generation photomultiplier tube (PMT)-based systems such as the Siemens Biograph mCT. These replacements warrant a direct comparison between the systems, so that a smooth transition in clinical practice and research is guaranteed, especially when quantitative values are used for dosimetry-based treatment guidance. The new generation digital PET scanners offer increased sensitivity. This could particularly benefit 90Y imaging, which tends to be very noisy owing to the small positron branching ratio and high random fraction of 90Y. This study aims to determine the ideal reconstruction settings for the digital Vision for quantitative 90Y imaging and to evaluate the image quality and quantification of the digital Vision in comparison with its predecessor, the PMT-based mCT, for 90Y imaging in radioembolisation procedures.

METHODS

The NEMA image quality phantom was scanned to determine the ideal reconstruction settings for the Vision. In addition, an anthropomorphic phantom was scanned with both the Vision and the mCT, mimicking a radioembolisation patient with lung, liver, tumour, and extrahepatic deposition inserts. Image quantification of the anthropomorphic phantom was assessed by the lung shunt fraction, the tumour to non-tumour ratio, the parenchymal dose, and the contrast to noise ratio of extrahepatic depositions.

RESULTS

For the Vision, a reconstruction with 3 iterations, 5 subsets, and no post-reconstruction filter is recommended for quantitative 90Y imaging, based on the convergence of the recovery coefficient. Comparing both systems showed that the noise level of the Vision is significantly lower than that of the mCT (background variability of 14% for the Vision and 25% for the mCT at 2.5·103 MBq for the 37 mm sphere size). For quantitative 90Y measures, such as needed in radioembolisation, both systems perform similarly.

CONCLUSIONS

We recommend to reconstruct 90Y images acquired on the Vision with 3 iterations, 5 subsets, and no post-reconstruction filter for quantitative imaging. The Vision provides a reduced noise level, but similar quantitative accuracy as compared with its predecessor the mCT.

Melanoma-Targeting Property of Y-90-Labeled Lactam-Cyclized α-Melanocyte-Stimulating Hormone Peptide

Purpose: The purpose of this study was to evaluate melanoma-targeting property of ⁹⁰Y-DOTA-GGNle-CycMSH_hex to facilitate its potential therapeutic application. Materials and Methods: DOTA-GGNle-CycMSH_hex was synthesized and readily labeled with ⁹⁰Y in 0.25 M NH₄Ac-buffered solution to generate ⁹⁰Y-DOTA-GGNle-CycMSH_hex. The specific receptor binding, internalization, and efflux of ⁹⁰Y-DOTA-GGNle-CycMSH_hex were determined on B16/F10 murine melanoma cells. The biodistribution property of ⁹⁰Y-DOTA-GGNle-CycMSH_hex was examined on B16/F10 melanoma-bearing C57 mice. Results: ⁹⁰Y-DOTA-GGNle-CycMSH_hex displayed receptor-specific binding, rapid internalization, and prolonged efflux on B16/F10 melanoma cells. ⁹⁰Y-DOTA-GGNle-CycMSH_hex exhibited high uptake and prolonged retention in melanoma, and fast urinary clearance on B16/F10 melanoma-bearing C57 mice. The B16/F10 tumor uptake was 20.73% ± 7.99%, 19.93% ± 5.73%, 14.8% ± 4.61%, and 6.69% ± 1.85% ID/g at 0.5, 2, 4, and 24 h postinjection, respectively. Conclusions: ⁹⁰Y-DOTA-GGNle-CycMSH_hex displayed melanocortin-1 receptor (MC1R) targeting and specificity on B16/F10 melanoma cells and tumors. The favorable melanoma-targeting property and fast urinary clearance of ⁹⁰Y-DOTA-GGNle-CycMSH_hex warranted its evaluation for melanoma therapy in future studies.

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.

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.

Dosimetry Estimate and Initial Clinical Experience with 90Y-PSMA-617

Because of different physical properties, the β-emitters ¹⁷⁷Lu and ⁹⁰Y offer specific radiologic-biologic advantages in dedicated clinical situations. Our objective was to introduce ⁹⁰Y-labeled prostate-specific membrane antigen (PSMA)-617 to clinical application, providing additional avenues for personalized medicine. Here, we present our dosimetry estimate for ⁹⁰Y-PSMA-617, report first clinical experiences, and discuss the advantages and drawbacks of varying the β-emitter in PSMA-targeting radioligand therapy. Methods: To approximate radiation dosimetry, 4 patients with metastatic castration-resistant prostate cancer underwent serially performed imaging up to 1 wk after ¹⁷⁷Lu-PSMA-617 therapy. Time-activity curves were extrapolated to the half-life of ⁹⁰Y, and OLINDA was used to calculate the dosimetry estimate. In clinical practice, 11 patients with PSMA-positive lymph-nodal bulk disease were stratified to receive ⁹⁰Y-PSMA-617 radioligand therapy (mean, 3.2 GBq; range, 2.8-3.7 GBq); afterward, safety lab tests, prostate-specific antigen (PSA) response, and clinical findings were thoroughly followed. Results: The projected dosimetry for ⁹⁰Y-PSMA-617 estimated a mean kidney dose of 3.47 ± 1.40 Gy/GBq, red marrow dose of 0.11 ± 0.04 Gy/GBq, and salivary gland dose of 5.57 ± 1.34 Gy/GBq; randomly chosen metastases were approximated with 22.8 ± 16.10 Gy/GBq. The observed acute hematologic toxicity (5 cases of leukopenia and 2 of thrombocytopenia, all grade 1 or 2) and clinical side effects (2 cases of transient xerostomia and 1 of nausea, all grade 1 or 2), as well as PSA response (any PSA response, 7/11 patients; >50% PSA decline, 5/11 patients), were comparable to ¹⁷⁷Lu-PSMA-617 literature data. Conclusion: A factor 3-4 lower treatment activity for ⁹⁰Y-PSMA-617 translates into a comparable dosimetry estimate and clinical findings similar to those of ¹⁷⁷Lu-PSMA-617. However, safety was demonstrated only for patients with oligometastatic disease. Further studies are needed to evaluate its potential in patients with more disseminated bone involvement or visceral metastasis.

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.