Radioembolization Dosimetry: The Road Ahead

MIRD Pamphlet No. 29: MIRDy90-A 90Y Research Microsphere Dosimetry Tool

90Y-microsphere radioembolization has become a well-established treatment option for liver malignancies and is one of the first U.S. Food and Drug Administration-approved unsealed radionuclide brachytherapy devices to incorporate dosimetry-based treatment planning. Several different mathematical models are used to calculate the patient-specific prescribed activity of 90Y, namely, body surface area (SIR-Spheres only), MIRD single compartment, and MIRD dual compartment (partition). Under the auspices of the MIRDsoft initiative to develop community dosimetry software and tools, the body surface area, MIRD single-compartment, MIRD dual-compartment, and MIRD multicompartment models have been integrated into a MIRDy90 software worksheet. The worksheet was built in MS Excel to estimate and compare prescribed activities calculated via these respective models. The MIRDy90 software was validated against available tools for calculating 90Y prescribed activity. The results of MIRDy90 calculations were compared with those obtained from vendor and community-developed tools, and the calculations agreed well. The MIRDy90 worksheet was developed to provide a vetted tool to better evaluate patient-specific prescribed activities calculated via different models, as well as model influences with respect to varying input parameters. MIRDy90 allows users to interact and visualize the results of various parameter combinations. Variables, equations, and calculations are described in the MIRDy90 documentation and articulated in the MIRDy90 worksheet. The worksheet is distributed as a free tool to build expertise within the medical physics community and create a vetted standard for model and variable management.

Navigating the nuances: comparative analysis and hyperparameter optimisation of neural architectures on contrast-enhanced MRI for liver and liver tumour segmentation

In medical imaging, accurate segmentation is crucial to improving diagnosis, treatment, or both. However, navigating the multitude of available architectures for automatic segmentation can be overwhelming, making it challenging to determine the appropriate type of architecture and tune the most crucial parameters during dataset optimisation. To address this problem, we examined and refined seven distinct architectures for segmenting the liver, as well as liver tumours, with a restricted training collection of 60 3D contrast-enhanced magnetic resonance images (CE-MRI) from the ATLAS dataset. Included in these architectures are convolutional neural networks (CNNs), transformers, and hybrid CNN/transformer architectures. Bayesian search techniques were used for hyperparameter tuning to hasten convergence to the optimal parameter mixes while also minimising the number of trained models. It was unexpected that hybrid models, which typically exhibit superior performance on larger datasets, would exhibit comparable performance to CNNs. The optimisation of parameters contributed to better segmentations, resulting in an average increase of 1.7% and 5.0% in liver and tumour segmentation Dice coefficients, respectively. In conclusion, the findings of this study indicate that hybrid CNN/transformer architectures may serve as a practical substitute for CNNs even in small datasets. This underscores the significance of hyperparameter optimisation.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Personalized 90 Y-resin microspheres dose determination: a retrospective study on the impact of dosimetry software on the treatment of patients with selective internal radiotherapy

INTRODUCTION

The calculation of resin yttrium-90 ( 90 Y-) microspheres activity for selective internal radiotherapy (SIRT) needs to be investigated.

METHODS AND MATERIALS

Analyses using Simplicit 90 Y (Boston Scientific, Natick, Massachusetts, USA) dosimetry software were performed to determine the concordance between the absorbed doses to the tumor (D T1 and D T2 ) and healthy liver (D N1 and D N2 ) during the pre-treatment and the post-treatment phases. An optimized calculation of the activity of 90 Y-microspheres performed using dosimetry software was applied retrospectively to assess the impact of this calculation method on the treatment.

RESULTS

D T1 ranged from 38.8 to 372 Gy, with a mean value of 128.9 ± 73.6 Gy and median of 121.2 Gy [interquartile range (IQR): 81.7-158.8 Gy]. The median D N1 and D N2 was 10.5 Gy (IQR: 5.8-17.6). A significant correlation was between D T1 and D T2 ( r  = 0.88, P  < 0.001) and D N1 and D N2 ( r  = 0.96, P  < 0.001). The optimized activities were calculated; the target dose to the tumor compartment was 120 Gy. No activity reduction was applied in accordance with the tolerance of the healthy liver. Optimization of the microspheres dosages would have resulted in a significant increase in activity for nine treatments (0.21-2.54 GBq) and a reduction for seven others (0.25-0.76 GBq).

CONCLUSIONS

The development of customized dosimetry software adapted to clinical practice makes it possible to use dosimetry to optimize the dosage for each patient.

Dose-response Analysis in Hepatic Tumors Treated with 90Y-TARE According to a Personalized Dosimetric Workflow: Preliminary Results

BACKGROUND

Transarterial Radioembolization (TARE) is a widespread radiation therapy for unresectable hepatic lesions, but a clear understanding of the dose-response link is still missing. The aim of this preliminary study is to investigate the role of both dosimetric and clinical parameters as classifiers or predictors of response and survival for TARE in hepatic tumors and to present possible response cut-off.

METHODS

20 patients treated with glass or resin microspheres according to a personalized workflow were included. Dosimetric parameters were extracted from personalized absorbed dose maps obtained from the convolution of 90Y PET images with 90Y voxel S-values.

RESULTS

D95 ≥ 104 Gy and tumor mean absorbed dose MADt ≥ 229 Gy were found to be optimal cut-off values for complete response, while D30 ≥ 180 Gy and MADt ≥ 117 Gy were selected as cut-off values for at least partial response and predicted better survival. Clinical parameters Alanine Transaminase (ALT) and Model for End-Stage Liver Disease (MELD) didn't show sufficient classification capability for response or survival.

CONCUSION

These preliminary results highlight the importance of an accurate dosimetric evaluation and suggest a cautious approach when considering clinical indicators. Dosimetric cut-off values could be a support tool in both planning and post-treatment phases. Larger multi-centric randomized trials, with standardized methods regarding patient selection, response criteria, Regions of Interest definition, dosimetric approach and activity planning are needed to confirm these promising results.

The Impact of Injection Distance to Bifurcations on Yttrium-90 Distribution in Liver Cancer Radioembolization

PURPOSE

To model the effect of the injection location on the distribution of yttrium-90 (⁹⁰Y) microspheres in the liver during radioembolization using computational simulation and to determine the potential effects of radial movements of the catheter tip.

MATERIALS AND METHODS

Numerical studies were conducted using images from a representative patient with hepatocellular carcinoma. The right hepatic artery (RHA) was segmented from contrast-enhanced cone-beam computed tomography scans. The blood flow was investigated in the trunk of the RHA using numerical simulations for 6 injection position scenarios at 2 sites located at a distance of approximately 5 and 20 mm upstream of the first bifurcation (RHA diameters of approximately 4.6 mm). The ⁹⁰Y delivery to downstream vessels was calculated from the simulated hepatic artery hemodynamics.

RESULTS

Varying the injection location along the RHA and across the vessel cross-section resulted in different simulated microsphere distributions in the downstream vascular bed. When the catheter tip was 5 mm upstream of the bifurcation, ⁹⁰Y distribution in the downstream branches varied by as much as 53% with a 1.5-mm radial movement of the tip. However, the catheter radial movement had a weaker effect on the microsphere distribution when the injection plane was farther from the first bifurcation (20 mm), with a maximum delivery variation of 9% to a downstream branch.

CONCLUSIONS

An injection location far from bifurcations is recommended to minimize the effect of radial movements of the catheter tip on the microsphere distribution.

Inflammatory markers and long term hematotoxicity of holmium-166-radioembolization in liver-dominant metastatic neuroendocrine tumors after initial peptide receptor radionuclide therapy

Purpose

In patients with neuroendocrine tumor liver metastases, additional tumor reduction can be achieved by sequential treatment with [¹⁶⁶Ho]-radioembolization after peptide receptor radionuclide therapy (PRRT). The aim of this study was to analyze hematotoxicity profiles, (i.e. lymphocyte and neutrophile toxicity) and the prognostic value of neutrophil-to-lymphocyte ratio (NLR) and thrombocyte-to-lymphocyte ratio (TLR).

Methods

All patients included in the prospective HEPAR PLuS study were included in this study. Blood testing was performed at baseline (before radioembolization) and at regular intervals during 1-year follow-up. Radiological response was assessed at 3, 6, 9, and 12 months according to RECIST 1.1. Logistic regression was used to analyze the prognostic value of NLR and TLR on response.

Results

Thirty-one patients were included in the toxicity analysis; thirty were included in the response analysis. Three weeks after radioembolization, a significant decrease in lymphocyte count (mean change − 0.26 × 10⁹/L) was observed. Ten patients (32.2%) experienced grade 3–4 lymphocyte toxicity. This normalized at 6 weeks and 3 months after treatment, while after 6 months a significant increase in lymphocyte count was observed. An increase in NLR and TLR at 3 weeks, compared to baseline, significantly predicted response at 3 months (AUC = 0.841 and AUC = 0.839, respectively) and at 6 months (AUC = 0.779 and AUC = 0.765). No significant relation with survival was found.

Conclusions

Toxicity after sequential treatment with PRRT and [¹⁶⁶Ho]-radioembolization is limited and temporary, while significant additional benefit can be expected. Change in NLR and TLR at 3-weeks follow-up may be valuable early predictors of response.

Trial registration ClinicalTrials.gov, NCT02067988. Registered 20 February 2014, https://clinicaltrials.gov/ct2/show/record/NCT02067988.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13550-022-00880-4.

Diagnostic and prognostic value of 99mTc-MAA SPECT/CT for treatment planning of 90Y-resin microsphere radioembolization for hepatocellular carcinoma: comparison with planar image

^99mTc-macroaggregated albumin (MAA) imaging is performed before transarterial radioembolization (TARE), in which SPECT/CT is presumed more precise than planar image. However, additive role of SPECT/CT has not been well established. Thirty-four consecutive hepatocellular carcinoma patients of intermediate and advanced stages who underwent ⁹⁰Y-microsphere TARE were recruited. On pre-treatment planning scan using ^99mTc-MAA, image characteristics and absorbed dose for target tumors calculated by partition model methods were estimated on planar image and SPECT/CT, respectively. The measurements were repeated on post-treatment ⁹⁰Y PET/CT, as the reference standard. Treatment response was assessed and predictive values of image parameters were analyzed. The image characteristics including heterogeneity, necrosis and thrombosis uptake were better delineated on SPECT/CT than planar scan. The agreement and correlation of TNr between SPECT/CT and PET/CT were stronger than those between planar scan and PET/CT. Tumor dose estimated on ^99mTc-MAA SPECT/CT was more effective than planar image for prediction of treatment response, with cutoff value 125 Gy (sensitivity of 86% and specificity of 75%). In conclusion, ^99mTc-MAA SPECT/CT is more closely correlated with post-treatment ⁹⁰Y PET/CT, and is more effective for predicting treatment response than planar scan. SPECT/CT is superior to planar image in simulation before ⁹⁰Y TARE.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Basics to advances in nanotherapy of colorectal cancer

Colorectal cancer (CRC) is the third most common cancer existing across the globe. It begins with the formation of polyps leading to the development of metastasis, especially in advanced stage patients, who necessitate intensive chemotherapy that usually results in a poor response and high morbidity owing to multidrug resistance and severe untoward effects to the non-cancerous cells. Advancements in the targeted drug delivery permit the targeting of tumor cells without affecting the non-tumor cells. Various nanocarriers such as liposomes, polymeric nanoparticles, carbon nanotubes, micelles, and nanogels, etc. are being developed and explored for effective delivery of cytotoxic drugs to the target site thereby enhancing the drug distribution and bioavailability, simultaneously subduing the side effects. Moreover, immunotherapy for CRC is being explored for last few decades. Few clinical trials have even potentially benefited patients suffering from CRC, still immunotherapy persists merely an experimental alternative. Assessment of the ongoing and completed trials is to be warranted for effective treatment of CRC. Scientists are paying efforts to develop novel carrier systems that may enhance the targeting potential of low therapeutic index chemo- and immune-therapeutics. Several preclinical studies have revealed the superior efficacy of nanotherapy in CRC as compared to conventional approaches. Clinical trials are being recruited to ascertain the safety and efficacy of CRC therapies. The present review discourses in a nutshell the molecular interventions including the genetics, signaling pathways involved in CRC, and advances in various strategies explored for the treatment of CRC with a special emphasis on nanocarriers based drug targeting.

Analysis of differences between 99mTc-MAA SPECT- and 90Y-microsphere PET-based dosimetry for hepatocellular carcinoma selective internal radiation therapy

Background

The aim of this study was to compare predictive and post-treatment dosimetry and analyze the differences, investigating factors related to activity preparation and delivery, imaging modality used, and interventional radiology.

Methods

Twenty-three HCC patients treated by selective internal radiation therapy with ⁹⁰Y glass microspheres were included in this study. Predictive and post-treatment dosimetry were calculated at the voxel level based on ^99mTc-MAA SPECT/CT and ⁹⁰Y-microsphere PET/CT respectively. Dose distribution was analyzed through mean dose, metrics extracted from dose-volume histograms, and Dice similarity coefficients applied on isodoses. Reproducibility of the radiological gesture and its influence on dose deviation was evaluated.

Results

⁹⁰Y delivered activity was lower than expected in 67% (16/24) of the cases mainly due to the residual activity. A mean deviation of − 6 ± 11% was observed between the delivered activity and the ⁹⁰Y PET’s FOV activity. In addition, a substantial difference of − 20 ± 8% was measured on ⁹⁰Y PET images between the activity in the liver and in the whole FOV. After normalization, ^99mTc-MAA SPECT dosimetry was highly correlated and concordant with ⁹⁰Y-microsphere PET dosimetry for all dose metrics evaluated (ρ = 0.87, ρ_c = 0.86, P = 3.10⁻⁸ and ρ = 0.91, ρ_c = 0.90, P = 7.10⁻¹⁰ for tumor and normal liver mean dose respectively for example). Besides, mean tumor dose deviation was lower when the catheter position was identical than when it differed (16 Gy vs. 37 Gy, P = 0.007). Concordance between predictive and post-treatment dosimetry, evaluated with Dice similarity coefficients applied on isodoses, significantly correlated with the distance of the catheter position from artery bifurcation (P = 0.04, 0.0004, and 0.05, for 50 Gy, 100 Gy, and 150 Gy isodoses respectively).

Conclusions

Discrepancies between planned activity and activity measured on ⁹⁰Y PET images were observed and seemed to be mainly related to clinical hazards and equipment issues. Predictive vs. post-treatment comparison of relative dose distributions between tumor and normal liver showed a good correlation and no significant difference highlighting the predictive value of ^99mTc MAA SPECT/CT-based dosimetry. Besides, the reproducibility of catheter tip position appears critical in the agreement between predictive and actual dose distribution.

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.

Fast technetium-99m liver SPECT for evaluation of the pretreatment procedure for radioembolization dosimetry

PURPOSE

The efficiency of radioembolization procedures could be greatly enhanced if results of the 99m Tc-MAA pretreatment procedure were immediately available in the interventional suite, enabling 1-day procedures as a result of direct estimation of the hepatic radiation dose and lung shunt fraction. This would, however, require a relatively fast, but still quantitative, SPECT procedure, which might be achieved with acquisition protocols using nonuniform durations of the projection images.

METHODS

SPECT liver images of the 150-MBq 99m Tc-MAA pretreatment procedure were simulated for eight different lesion locations and two different lesion sizes using the digital XCAT phantom for both single- and dual-head scanning geometries with respective total acquisition times of 1, 2, 5, 10, and 30 min. Three nonuniform projection-time acquisition protocols ("half-circle SPECT (HCS)," "nonuniform SPECT (NUS) I," and "NUS II") for fast quantitative SPECT of the liver were designed and compared with the standard uniform projection-time protocol. Images were evaluated in terms of contrast-to-noise ratio (CNR), activity recovery coefficient (ARC), tumor/non-tumor (T/N) activity concentration ratio, and lung shunt fraction (LSF) estimation. In addition, image quality was verified with a physical phantom experiment, reconstructed with both clinical and Monte Carlo-based reconstruction software.

RESULTS

Simulations showed no substantial change in image quality and dosimetry by usage of a nonuniform projection-time acquisition protocol. Upon shortening acquisition times, CNR dropped, but ARC, T/N ratio, and LSF estimates were stable across all simulated acquisition times. Results of the physical phantom were in agreement with those of the simulations.

CONCLUSION

Both uniform and nonuniform projection-time acquisition liver SPECT protocols yield accurate dosimetric metrics for radioembolization treatment planning in the interventional suite within 10 min, without compromising image quality. Consequently, fast quantitative SPECT of the liver in the interventional suite is feasible.

Radionuclide therapy: current status and prospects for internal dosimetry in individualized therapeutic planning

The efficacy and toxicity of radionuclide therapy are believed to be directly related to the radiation doses received by target tissues; however, nuclear medicine therapy continues to be based primarily on the administration of empirical activities to patients and less frequently on the use of internal dosimetry for individual therapeutic planning. This review aimed to critically describe the techniques and clinical evidence of dosimetry as a tool for therapeutic planning and the main limitations to its implementation in clinical practice. The present article is a nonsystematic review of voxel-based dosimetry. Clinical evidence pointing to a correlation between the radiation dose and therapeutic response in various diseases, such as thyroid carcinoma, neuroendocrine tumors and prostate cancer, is reviewed. Its limitations include technical aspects related to image acquisition and processing and the lack of randomized clinical trials demonstrating the impact of dosimetry on patient therapy. A more widespread use of dosimetry in therapeutic planning involves the development of user-friendly dosimetric protocols and confirmation that dose estimation implies good efficacy and low treatment-related toxicity.

Impact of 90Y PET gradient-based tumor segmentation on voxel-level dosimetry in liver radioembolization

BACKGROUND

The purpose was to validate ⁹⁰Y PET gradient-based tumor segmentation in phantoms and to evaluate the impact of the segmentation method on reported tumor absorbed dose (AD) and biological effective dose (BED) in ⁹⁰Y microsphere radioembolization (RE) patients. A semi-automated gradient-based method was applied to phantoms and patient tumors on the ⁹⁰Y PET with the initial bounding volume for gradient detection determined from a registered diagnostic CT or MR; this PET-based segmentation (PS) was compared with radiologist-defined morphologic segmentation (MS) on CT or MRI. AD and BED volume histogram metrics (D90, D70, mean) were calculated using both segmentations and concordance/correlations were investigated. Spatial concordance was assessed using Dice similarity coefficient (DSC) and mean distance to agreement (MDA). PS was repeated to assess intra-observer variability.

RESULTS

In phantoms, PS demonstrated high accuracy in lesion volumes (within 15%), AD metrics (within 11%), high spatial concordance relative to morphologic segmentation (DSC > 0.86 and MDA < 1.5 mm), and low intra-observer variability (DSC > 0.99, MDA < 0.2 mm, AD/BED metrics within 2%). For patients (58 lesions), spatial concordance between PS and MS was degraded compared to in-phantom (average DSC = 0.54, average MDA = 4.8 mm); the average mean tumor AD was 226 ± 153 and 197 ± 138 Gy, respectively for PS and MS. For patient AD metrics, the best Pearson correlation (r) and concordance correlation coefficient (ccc) between segmentation methods was found for mean AD (r = 0.94, ccc = 0.92), but worsened as the metric approached the minimum dose (for D90, r = 0.77, ccc = 0.69); BED metrics exhibited a similar trend. Patient PS showed low intra-observer variability (average DSC = 0.81, average MDA = 2.2 mm, average AD/BED metrics within 3.0%).

CONCLUSIONS

⁹⁰Y PET gradient-based segmentation led to accurate/robust results in phantoms, and showed high concordance with MS for reporting mean tumor AD/BED in patients. However, tumor coverage metrics such as D90 exhibited worse concordance between segmentation methods, highlighting the need to standardize segmentation methods when reporting AD/BED metrics from post-therapy ⁹⁰Y PET. Estimated differences in reported AD/BED metrics due to segmentation method will be useful for interpreting RE dosimetry results in the literature including tumor response data.

The physics of radioembolization

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

Additional hepatic 166Ho-radioembolization in patients with neuroendocrine tumours treated with 177Lu-DOTATATE; a single center, interventional, non-randomized, non-comparative, open label, phase II study (HEPAR PLUS trial)

Background

Neuroendocrine tumours (NET) consist of a heterogeneous group of neoplasms with various organs of origin. At diagnosis 21% of the patients with a Grade 1 NET and 30% with a Grade 2 NET have distant metastases. Treatment with peptide receptor radionuclide therapy (PRRT) shows a high objective response rate and long median survival after treatment. However, complete remission is almost never achieved. The liver is the most commonly affected organ in metastatic disease and is the most incriminating factor for patient survival. Additional treatment of liver disease after PRRT may improve outcome in NET patients. Radioembolization is an established therapy for liver metastasis. To investigate this hypothesis, a phase 2 study was initiated to assess effectiveness and toxicity of holmium-166 radioembolization (¹⁶⁶Ho-RE) after PRRT with lutetium-177 (¹⁷⁷Lu)-DOTATATE.

Methods

The HEPAR PLUS trial (“HolmiumEmbolizationParticles forArterialRadiotherapyPlus¹⁷⁷Lu-DOTATATE inSalvage NET patients”) is a single centre, interventional, non-randomized, non-comparative, open label study. In this phase 2 study 30–48 patients with > 3 measurable liver metastases according to RECIST 1.1 will receive additional ¹⁶⁶Ho-RE within 20 weeks after the 4th and last cycle of PRRT with 7.4 GBq ¹⁷⁷Lu-DOTATATE. Primary objectives are to assess tumour response, complete and partial response according to RECIST 1.1, and toxicity, based on CTCAE v4.03, 3 months after ¹⁶⁶Ho-RE. Secondary endpoints include biochemical response, quality of life, biodistribution and dosimetry.

Discussion

This is the first prospective study to combine PRRT with ¹⁷⁷Lu-DOTATATE and additional ¹⁶⁶Ho-RE in metastatic NET. A radiation boost on intrahepatic disease using ¹⁶⁶Ho-RE may lead to an improved response rate without significant additional side-effects.

Trial registration

Clinicaltrials.gov NCT02067988, 13 February 2014. Protocol version: 6, 30 november 2016.

Y-90 SPECT ML image reconstruction with a new model for tissue-dependent bremsstrahlung production using CT information: a proof-of-concept study

While the yield of positrons used in Y-90 PET is independent of tissue media, Y-90 SPECT imaging is complicated by the tissue dependence of bremsstrahlung photon generation. The probability of bremsstrahlung production is proportional to the square of the atomic number of the medium. Hence, the same amount of activity in different tissue regions of the body will produce different numbers of bremsstrahlung photons. Existing reconstruction methods disregard this tissue-dependency, potentially impacting both qualitative and quantitative imaging of heterogeneous regions of the body such as bone with marrow cavities. In this proof-of-concept study, we propose a new maximum-likelihood method that incorporates bremsstrahlung generation probabilities into the system matrix, enabling images of the desired Y-90 distribution to be reconstructed instead of the 'bremsstrahlung distribution' that is obtained with existing methods. The tissue-dependent probabilities are generated by Monte Carlo simulation while bone volume fractions for each SPECT voxel are obtained from co-registered CT. First, we demonstrate the tissue dependency in a SPECT/CT imaging experiment with Y-90 in bone equivalent solution and water. Visually, the proposed reconstruction approach better matched the true image and the Y-90 PET image than the standard bremsstrahlung reconstruction approach. An XCAT phantom simulation including bone and marrow regions also demonstrated better agreement with the true image using the proposed reconstruction method. Quantitatively, compared with the standard reconstruction, the new method improved estimation of the liquid bone:water activity concentration ratio by 40% in the SPECT measurement and the cortical bone:marrow activity concentration ratio by 58% in the XCAT simulation.

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.

Estimation of lung shunt fraction from simultaneous fluoroscopic and nuclear images

Radioembolisation with yttrium-90 (⁹⁰Y) is increasingly used as a treatment of unresectable liver malignancies. For safety, a scout dose of technetium-99m macroaggregated albumin (^99mTc-MAA) is used prior to the delivery of the therapeutic activity to mimic the deposition of ⁹⁰Y. One-day procedures are currently limited by the lack of nuclear images in the intervention room. To cope with this limitation, an interventional simultaneous fluoroscopic and nuclear imaging device is currently being developed. The purpose of this simulation study was to evaluate the accuracy of estimating the lung shunt fraction (LSF) of the scout dose in the intervention room with this device and compare it against current clinical methods.

METHODS

A male and female XCAT phantom, both with two respiratory profiles, were used to simulate various LSFs resulting from a scout dose of 150 MBq ^99mTc-MAA. Hybrid images were Monte Carlo simulated for breath-hold (5 s) and dynamic breathing (10 frames of 0.5 s) acquisitions. Nuclear images were corrected for attenuation with the fluoroscopic image and for organ overlap effects using a pre-treatment CT-scan. For comparison purposes, planar scintigraphy and mobile gamma camera images (both 300 s acquisition time) were simulated. Estimated LSFs were evaluated for all methods and compared to the phantom ground truth.

RESULTS

In the clinically relevant range of 10-20% LSF, hybrid imaging overestimated LSF with approximately 2 percentage points (pp) and 3 pp for the normal and irregular breathing phantoms, respectively. After organ overlap correction, LSF was estimated with a more constant error. Errors in planar scintigraphy and mobile gamma camera imaging were more dependent on LSF, body shape and breathing profile.

CONCLUSION

LSF can be estimated with a constant minor error with a hybrid imaging device. Estimated LSF is highly dependent on true LSF, body shape and breathing pattern when estimated with current clinical methods. The hybrid imaging device is capable of accurately estimating LSF within a few seconds in an interventional setting.

Efficacy of Radioembolization with 166Ho-Microspheres in Salvage Patients with Liver Metastases: A Phase 2 Study

Radioembolization of liver malignancies with ¹⁶⁶Ho-microspheres has been shown to be safe in a phase 1 dose-escalation study. The purpose of this study was to investigate the efficacy of ¹⁶⁶Ho radioembolization. Methods: In this prospective single-arm study, 56 patients were enrolled, all with liver metastases refractory to systemic therapy and ineligible for surgical resection. The primary outcome was a response by 2 target lesions on triphasic liver CT scans 3 mo after therapy, as assessed using RECIST, version 1.1. Secondary outcomes included overall tumor response, time to imaging progression, overall survival, toxicity, quality of life, and quantification of the microspheres on SPECT and MRI. Results: Between May 2012 and March 2015, 38 eligible patients were treated, one of whom was not evaluable. In 27 (73%) of 37 patients, the target lesions showed complete response, partial response, or stable disease (disease control) at 3 mo (95% confidence interval [CI], 57%-85%). The median overall survival was 14.5 mo (95% CI, 8.6-22.8 mo). For colorectal cancer patients (n = 23), the median overall survival was 13.4 mo (95% CI, 8.2-15.7 mo). Grade 3 or 4 toxic events after treatment (according to the Common Terminology Criteria for Adverse Events, version 4.03) included abdominal pain (in 18% of patients), nausea (8%), ascites (3%), fatigue (3%), gastric stenosis (3%), hepatic failure (3%), liver abscesses (3%), paroxysmal atrial tachycardia (3%), thoracic pain (3%), upper gastrointestinal hemorrhage (3%), and vomiting (3%). On SPECT, ¹⁶⁶Ho could be quantified with high accuracy and precision, with a mean overestimation of 9.3% ± 7.1% in the liver. Conclusion: Radioembolization with ¹⁶⁶Ho-microspheres induced a tumor response with an acceptable toxicity profile in salvage patients with liver metastases.

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.

Yttrium-90 hepatic radioembolization: clinical review and current techniques in interventional radiology and personalized dosimetry

In recent years, yttrium-90 ((90)Y) microsphere radioembolization has been establishing itself as a safe and efficacious treatment for both primary and metastatic liver cancers. This extends to both first-line therapies as well as in the salvage setting. In addition, radioembolization appears efficacious for patients with portal vein thrombosis, which is currently a contraindication for surgery, transplantation and transarterial chemoembolization. This article reviews the efficacy and expanding use of (90)Y microsphere radioembolization with an added emphasis on recent advances in personalized dosimetry and interventional radiology techniques. Directions for future research into combination therapies with radioembolization and expansion into sites other than the liver are also explored.

A descriptive analysis of remnant activity during (90)Y resin microspheres radioembolization of hepatic tumors: technical factors and dosimetric implications

OBJECTIVE

Activity planning for (90)Y radioembolization aims to maximize the effect of the treatment while keeping toxicity acceptably low. Our aim was to describe the amount of residual activity in post-treatment v-vials and tubing and analyze the possible factors affecting it (total activity administered, number of split activity injection(s), previous treatments, administration artery and microcatheter size), as these may influence dosimetric planning and treatment.

METHODS

This was a retrospective review using case records of patients who received (90)Y-radioembolization for hepatic tumors at a single tertiary center. From August 2013 to September 2015, seventy-seven out of one hundred and fifty patients who received radioembolization with (90)Y resin microspheres due to inoperable Hepatocellular Carcinoma (HCC) or liver metastases were included. The rest were mainly excluded due to incomplete data sets. The number of split activities (injections) for the radioembolization could be: one single injection, two or three. The remnant activity in post-treatment v-vials and tubing were measured for every patient. The administration arteries evaluated were: proper hepatic artery (PHA), right hepatic artery (RHA), middle hepatic artery (MHA), left hepatic artery (LHA) and small caliber branch arteries. The sizes of the microcatheters (2.2 or 2.7 Fr) used to administer the dose were also evaluated.

RESULTS

In total, 77 out of 150 patients were included in the final analysis. There were 59 men of median age 64.0 years old. The total median dose loss was 0.10 GBq. The total dose loss increased 0.244 GBq [95 % CI = (0.169, 0.318)] when three split activities were given compared to single activity injection. Activity loss for each injection increased 0.0297 GBq [95 % CI = (0.0151, 0.0443)] for every 1.0 GBq increase of split activity injection. There were no significant statistical differences in the rest of patient characteristics.

CONCLUSIONS

There is significant loss of activity observed during radioembolization, which can have a major dosimetric impact. The total administered activity and the number of split injections during radioembolization are the main influencing factors. Further prospective studies as well as measures of clinical outcome are warranted.

⁹⁰Y Hepatic Radioembolization: An Update on Current Practice and Recent Developments

Radioembolization is an established treatment modality that has been subjected to many improvements over the last decade. Developments are occurring at a high pace, affecting patient selection and treatment. The aim of this review is therefore to provide an overview of current practice, with a focus on recent developments in the field of radioembolization. Several practical issues and recommendations in the application of radioembolization will be discussed, ranging from patient selection to treatment response and future applications.