Increased absorbed liver dose in Selective Internal Radiation Therapy (SIRT) correlates with increased sphere-cluster frequency and absorbed dose inhomogeneity

The role of biomarkers and dosimetry parameters in overall and progression free survival prediction for patients treated with personalized 90Y glass microspheres SIRT: a preliminary machine learning study

BACKGROUND

Overall Survival (OS) and Progression-Free Survival (PFS) analyses are crucial metrics for evaluating the efficacy and impact of treatment. This study evaluated the role of clinical biomarkers and dosimetry parameters on survival outcomes of patients undergoing 90Y selective internal radiation therapy (SIRT).

MATERIALS/METHODS

This preliminary and retrospective analysis included 17 patients with hepatocellular carcinoma (HCC) treated with 90Y SIRT. The patients underwent personalized treatment planning and voxel-wise dosimetry. After the procedure, the OS and PFS were evaluated. Three structures were delineated including tumoral liver (TL), normal perfused liver (NPL), and whole normal liver (WNL). 289 dose-volume constraints (DVCs) were extracted from dose-volume histograms of physical and biological effective dose (BED) maps calculated on 99mTc-MAA and 90Y SPECT/CT images. Subsequently, the DVCs and 16 clinical biomarkers were used as features for univariate and multivariate analysis. Cox proportional hazard ratio (HR) was employed for univariate analysis. HR and the concordance index (C-Index) were calculated for each feature. Using eight different strategies, a cross-combination of various models and feature selection (FS) methods was applied for multivariate analysis. The performance of each model was assessed using an averaged C-Index on a three-fold nested cross-validation framework. The Kaplan-Meier (KM) curve was employed for univariate and machine learning (ML) model performance assessment.

RESULTS

The median OS was 11 months [95% CI: 8.5, 13.09], whereas the PFS was seven months [95% CI: 5.6, 10.98]. Univariate analysis demonstrated the presence of Ascites (HR: 9.2[1.8,47]) and the aim of SIRT (segmentectomy, lobectomy, palliative) (HR: 0.066 [0.0057, 0.78]), Aspartate aminotransferase (AST) level (HR:0.1 [0.012-0.86]), and MAA-Dose-V205(%)-TL (HR:8.5[1,72]) as predictors for OS. 90Y-derived parameters were associated with PFS but not with OS. MAA-Dose-V205(%)-WNL, MAA-BED-V400(%)-WNL with (HR:13 [1.5-120]) and 90Y-Dose-mean-TL, 90Y-D50-TL-Gy, 90Y-Dose-V205(%)-TL, 90Y-Dose- D50-TL-Gy, and 90Y-BED-V400(%)-TL (HR:15 [1.8-120]) were highly associated with PFS among dosimetry parameters. The highest C-index observed in multivariate analysis using ML was 0.94 ± 0.13 obtained from Variable Hunting-variable-importance (VH.VIMP) FS and Cox Proportional Hazard model predicting OS, using clinical features. However, the combination of VH. VIMP FS method with a Generalized Linear Model Network model predicting OS using Therapy strategy features outperformed the other models in terms of both C-index and stratification of KM curves (C-Index: 0.93 ± 0.14 and log-rank p-value of 0.023 for KM curve stratification).

CONCLUSION

This preliminary study confirmed the role played by baseline clinical biomarkers and dosimetry parameters in predicting the treatment outcome, paving the way for the establishment of a dose-effect relationship. In addition, the feasibility of using ML along with these features was demonstrated as a helpful tool in the clinical management of patients, both prior to and following 90Y-SIRT.

A microscopic model of the dose distribution in hepatocellular carcinoma after selective internal radiation therapy

The dosimetry evaluation for the selective internal radiation therapy is currently performed assuming a uniform activity distribution, which is in contrast with literature findings. A 2D microscopic model of the perfused liver was developed to evaluate the effect of two different 90Y microspheres distributions: i) homogeneous partitioning with the microspheres equally distributed in the perfused liver, and ii) tumor-clustered partitioning where the microspheres distribution is inferred from the patient specific images.

METHODS

Two subjects diagnosed with liver cancer were included in this study. For each subject, abdominal CT scans acquired prior to the SIRT and post-treatment 90Y positron emission tomography were considered. Two microspheres partitionings were simulated namely homogeneous and tumor-clustered partitioning. The homogeneous and tumor-clustered partitionings were derived starting from CT images. The microspheres radiation is simulated by means of Russell's law.

RESULTS

In homogenous simulations, the dose delivery is uniform in the whole liver while in the tumor-clustered simulations a heterogeneous distribution of the delivered dose is visible with higher values in the tumor regions. In addition, in the tumor-clustered simulation, the delivered dose is higher in the viable tumor than in the necrotic tumor, for all patients. In the tumor-clustered case, the dose delivered in the non-tumoral tissue (NTT) was considerably lower than in the perfused liver.

CONCLUSIONS

The model proposed here represents a proof-of-concept for personalized dosimetry assessment based on preoperative CT images.

MIDOS: a novel stochastic model towards a treatment planning system for microsphere dosimetry in liver tumors

PURPOSE

Transarterial radioembolization (TARE) procedures treat liver tumors by injecting radioactive microspheres into the hepatic artery. Currently, there is a critical need to optimize TARE towards a personalized dosimetry approach. To this aim, we present a novel microsphere dosimetry (MIDOS) stochastic model to estimate the activity delivered to the tumor(s), normal liver, and lung.

METHODS

MIDOS incorporates adult male/female liver computational phantoms with the hepatic arterial, hepatic portal venous, and hepatic venous vascular trees. Tumors can be placed in both models at user discretion. The perfusion of microspheres follows cluster patterns, and a Markov chain approach was applied to microsphere navigation, with the terminal location of microspheres determined to be in either normal hepatic parenchyma, hepatic tumor, or lung. A tumor uptake model was implemented to determine if microspheres get lodged in the tumor, and a probability was included in determining the shunt of microspheres to the lung. A sensitivity analysis of the model parameters was performed, and radiation segmentectomy/lobectomy procedures were simulated over a wide range of activity perfused. Then, the impact of using different microspheres, i.e., SIR-Sphere®, TheraSphere®, and QuiremSphere®, on the tumor-to-normal ratio (TNR), lung shunt fraction (LSF), and mean absorbed dose was analyzed.

RESULTS

Highly vascularized tumors translated into increased TNR. Treatment results (TNR and LSF) were significantly more variable for microspheres with high particle load. In our scenarios with 1.5 GBq perfusion, TNR was maximum for TheraSphere® at calibration time in segmentectomy/lobar technique, for SIR-Sphere® at 1-3 days post-calibration, and regarding QuiremSphere® at 3 days post-calibration.

CONCLUSION

This novel approach is a decisive step towards developing a personalized dosimetry framework for TARE. MIDOS assists in making clinical decisions in TARE treatment planning by assessing various delivery parameters and simulating different tumor uptakes. MIDOS offers evaluation of treatment outcomes, such as TNR and LSF, and quantitative scenario-specific decisions.

EANM procedure guideline for the treatment of liver cancer and liver metastases with intra-arterial radioactive compounds

Primary liver tumours (i.e. hepatocellular carcinoma (HCC) or intrahepatic cholangiocarcinoma (ICC)) are among the most frequent cancers worldwide. However, only 10-20% of patients are amenable to curative treatment, such as resection or transplant. Liver metastases are most frequently caused by colorectal cancer, which accounts for the second most cancer-related deaths in Europe. In both primary and secondary tumours, radioembolization has been shown to be a safe and effective treatment option. The vast potential of personalized dosimetry has also been shown, resulting in markedly increased response rates and overall survival. In a rapidly evolving therapeutic landscape, the role of radioembolization will be subject to changes. Therefore, the decision for radioembolization should be taken by a multidisciplinary tumour board in accordance with the current clinical guidelines. The purpose of this procedure guideline is to assist the nuclear medicine physician in treating and managing patients undergoing radioembolization treatment. PREAMBLE: The European Association of Nuclear Medicine (EANM) is a professional non-profit medical association that facilitates communication worldwide among individuals pursuing clinical and research excellence in nuclear medicine. The EANM was founded in 1985. These guidelines are intended to assist practitioners in providing appropriate nuclear medicine care for patients. They are not inflexible rules or requirements of practice and are not intended, nor should they be used, to establish a legal standard of care. The ultimate judgment regarding the propriety of any specific procedure or course of action must be made by medical professionals taking into account the unique circumstances of each case. Thus, there is no implication that an approach differing from the guidelines, standing alone, is below the standard of care. To the contrary, a conscientious practitioner may responsibly adopt a course of action different from that set out in the guidelines when, in the reasonable judgment of the practitioner, such course of action is indicated by the condition of the patient, limitations of available resources or advances in knowledge or technology subsequent to publication of the guidelines. The practice of medicine involves not only the science but also the art of dealing with the prevention, diagnosis, alleviation and treatment of disease. The variety and complexity of human conditions make it impossible to always reach the most appropriate diagnosis or to predict with certainty a particular response to treatment. Therefore, it should be recognised that adherence to these guidelines will not ensure an accurate diagnosis or a successful outcome. All that should be expected is that the practitioner will follow a reasonable course of action based on current knowledge, available resources and the needs of the patient to deliver effective and safe medical care. The sole purpose of these guidelines is to assist practitioners in achieving this objective.

Radioembolization of Hepatocellular Carcinoma with 90Y Glass Microspheres: No Advantage of Voxel Dosimetry with Respect to Mean Dose in Dose-Response Analysis with Two Radiological Methods

Simple Summary

We confirmed that the non-uniformity of an intra-lesion dose distribution, which was introduced in calculations as voxel dosimetry, did not significantly improve the AUC values of the dose–response relationship with respect to the mean dose. This was probably derived from the strong correlations (all p < 0.0001) among all voxel-based dosimetric variables (minimum Spearman correlation coefficient: 0.67) caused by the limited spatial resolution of nuclear medicine images. Responses were assessed with mRECIST and with an experimental densitometric method with a response threshold optimized at 20% HU variation. Significant dose–response agreement was obtained only with the densitometric method and only with post-therapy 90Y-PET data. More unexpectedly, the injection of Theraspheres™ on day 8 from the reference date rather than on day 4 worsened the dose–response correlation and reduced the efficacy at high doses. This may be explained by the increased non-uniformity following the non-linear mega-clustering effect triggered by the higher number of microspheres/GBq injected on day 8.

Abstract

In this confirmatory study, we tested if a calculation that included the non-uniformity of dose deposition through a voxel-based dosimetric variable Ψ was able to improve the dose–response agreement with respect to the mean absorbed dose D. We performed dosimetry with ^99mTc-MAA SPECT/CT and ⁹⁰Y-PET/CT in 86 patients treated 8 instead of 4 days after the reference date with 2.8 times more ⁹⁰Y glass microspheres/GBq than in our previous study. The lesion-by-lesion response was assessed with the mRECIST method and with an experimental densitometric criterion. A total of 106 lesions were studied. Considering Ψ as a prognostic response marker, having no Ψ provided a significantly higher AUC than D. The correlation, t-test, and AUC values were statistically significant only with the densitometric method and only with post-therapy dosimetry. In comparison with our previous study, the dose–response correlation and AUC values were poorer (maximum r = 0.43, R² = 0.14, maximal AUC = 0.71), and the efficacy at a high dose did not reach 100%. The expected advantages of voxel dosimetry were nullified by the correlation between any Ψ and D due to the limited image spatial resolution. The lower AUC and efficacy may be explained by the mega-clustering effect triggered by the higher number of microspheres/GBq injected on day 8.

Microspheres Used in Liver Radioembolization: From Conception to Clinical Effects

Inert microspheres, labeled with several radionuclides, have been developed during the last two decades for the intra-arterial treatment of liver tumors, generally called Selective Intrahepatic radiotherapy (SIRT). The aim is to embolize microspheres into the hepatic capillaries, accessible through the hepatic artery, to deliver high levels of local radiation to primary (such as hepatocarcinoma, HCC) or secondary (metastases from several primary cancers, e.g., colorectal, melanoma, neuro-endocrine tumors) liver tumors. Several types of microspheres were designed as medical devices, using different vehicles (glass, resin, poly-lactic acid) and labeled with different radionuclides, 90Y and 166Ho. The relationship between the microspheres' properties and the internal dosimetry parameters have been well studied over the last decade. This includes data derived from the clinics, but also computational data with various millimetric dosimetry and radiobiology models. The main purpose of this paper is to define the characteristics of these radiolabeled microspheres and explain their association with the microsphere distribution in the tissues and with the clinical efficacy and toxicity. This review focuses on avenues to follow in the future to optimize such particle therapy and benefit to patients.

Multiscale Computational Fluid Dynamics Modeling for Personalized Liver Cancer Radioembolization Dosimetry

Yttrium-90 (90Y) radioembolization is a minimally invasive procedure increasingly used for advanced liver cancer treatment. In this method, radioactive microspheres are injected into the hepatic arterial bloodstream to target, irradiate, and kill cancer cells. Accurate and precise treatment planning can lead to more efficient and safer treatment by delivering a higher radiation dose to the tumor while minimizing the exposure of the surrounding liver parenchyma. Treatment planning primarily relies on the estimated radiation dose delivered to tissue. However, current methods used to estimate the dose are based on simplified assumptions that make the dosimetry results unreliable. In this work, we present a computational model to predict the radiation dose from the 90Y activity in different liver segments to provide a more realistic and personalized dosimetry. Computational fluid dynamics (CFD) simulations were performed in a 3D hepatic arterial tree model segmented from cone-beam CT angiographic data obtained from a patient with hepatocellular carcinoma (HCC). The microsphere trajectories were predicted from the velocity field. 90Y dose distribution was then calculated from the volumetric distribution of the microspheres. Two injection locations were considered for the microsphere administration, a lobar and a selective injection. Results showed that 22% and 82% of the microspheres were delivered to the tumor, after each injection, respectively, and the combination of both injections ultimately delivered 49% of the total administered 90Y microspheres to the tumor. Results also illustrated the nonhomogeneous distribution of microspheres between liver segments, indicating the importance of developing patient-specific dosimetry methods for effective radioembolization treatment.

Monte Carlo 90Y PET/CT dosimetry of unexpected focal radiation-induced lung damage after hepatic radioembolisation

Transarterial radioembolization (TARE) with ⁹⁰Y-loaded microspheres is an established therapeutic option for inoperable hepatic tumors. Increasing knowledge regarding TARE hepatic dose-response and dose-toxicity correlation is available but few studies have investigated dose-toxicity correlation in extra-hepatic tissues. We investigated absorbed dose levels for the appearance of focal lung damage in a case of off-target deposition of ⁹⁰Y microspheres and compared them with the corresponding thresholds recommended to avoiding radiation induced lung injury following TARE. A 64-year-old male patient received 1.6 GBq of ⁹⁰Y-labelled glass microspheres for an inoperable left lobe hepatocellular carcinoma. A focal off-target accumulation of radiolabeled microspheres was detected in the left lung upper lobe at the post-treatment ⁹⁰Y-PET/CT, corresponding to a radiation-induced inflammatory lung lesion at the 3-months ¹⁸F-FDG PET/CT follow-up. ⁹⁰Y-PET/CT data were used as input for Monte-Carlo based absorbed dose estimations. Dose-volume-histograms were computed to characterize the heterogeneity of absorbed dose distribution. The dose level associated with the appearance of lung tissue damage was estimated as the median absorbed dose measured at the edge of the inflammatory nodule. To account for respiratory movements and possible inaccuracy of image co-registration, three different methods were evaluated to define the irradiated off-target volume. Monte Carlo-derived absorbed dose distribution showed a highly heterogeneous absorbed dose pattern at the site of incidental microsphere deposition (volume = 2.13 ml) with a maximum dose of 630 Gy. Absorbed dose levels ranging from 119 Gy to 133 Gy, were estimated at the edge of the inflammatory nodule, depending on the procedure used to define the target volume. This report describes an original Monte Carlo based patient-specific dosimetry methodology for the study of the radiation-induced damage in a focal lung lesion after TARE. In our patient, radiation-induced focal lung damage occurred at significantly higher absorbed doses than those considered for single administration or cumulative lung dose delivered during TARE.

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

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

Verification of a method to detect glass microspheres via micro-CT

INTRODUCTION

The ability to determine the microscopic distribution of glass microspheres in ⁹⁰ Y radioembolization has important applications in post-treatment microdosimetry and cluster analysis. Current methods are time-intensive and labor-intensive and thus are typically only applied to small samples.

MATERIALS AND METHODS

A high-resolution micro-CT image with a voxel size of 8.74 µm was acquired of phantoms containing ~25 µm-diameter glass microspheres embedded in tissue-equivalent materials that were optically transparent, which allowed true microsphere locations to be determined using transmission light microscopy. A 3-stage algorithm was developed to estimate the number and locations of microspheres in tissue regions. The stages are thresholding the CT image and discarding regions with insufficient voxels, estimating the number of microspheres in each region using the values of the detected and neighboring region voxels and estimating locations for each microsphere using the outputs of the previous two stages. Two different methods for estimating the number of microspheres in each region were derived, as were five methods for localizing microspheres. Metrics for each stage were computed, and the mean absolute error (MAE) between the dose to 72 µm voxels of the true and estimated dose maps created from the microsphere locations was used as the figure of merit for overall algorithm performance. Microsphere locations identified in the optical micrograph were used as the gold standard for the metrics of all stages. The method's utility was then demonstrated using a specimen from a human neuroendocrine tumor (NET) treated with glass ⁹⁰ Y microspheres.

RESULTS

The stage detecting regions containing microspheres found 100% of microspheres inside regions. The number of incorrectly detected regions without microspheres was 1.5% of the total number of regions. In stage 2, with these parameters, nearly 94% of the actual number of spheres in each region was correctly counted, and only 5% of the estimated sphere quantities in each region were false positives. The MAE between the true dose maps and dose maps estimated using the full algorithm with optimal parameter and method choices was 4.2%. A total of 5,713 glass microspheres were identified as being distributed heterogeneously in the NET specimen with a maximum tumor dose of >2500 Gy and 46% of the specimen receiving <20 Gy.

CONCLUSIONS

This work developed and evaluated a method to detect and estimate the three-dimensional locations of glass microspheres in whole tissue samples that require less manual effort than traditional methods. This method could be used to gain important insights into the heterogeneity of microsphere distributions that would be useful for improving radioembolization treatment planning.

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.

Development of a Customizable Hepatic Arterial Tree and Particle Transport Model for Use in Treatment Planning

Optimal treatment planning for radioembolization of hepatic cancers produces sufficient dose to tumors for control and dose to normal liver parenchyma that is below the threshold for toxicity. The non-uniform distribution of particles in liver microanatomy complicates the planning process as different functional regions receive different doses. Having realistic and patient-specific models of the arterial tree and microsphere trapping would be useful for developing more optimal treatment plans. We propose a macrocell-based growth method to generate models of the hepatic arterial tree from the proper hepatic artery to the terminal arterioles supplying the capillaries in the parenchyma. We show how these trees can be adapted to match patient values of pressure, flow, and vessel diameters while still conforming to laws controlling vessel bifurcation, changes in pressure, and blood flow. We also introduce a method to model particle transport within the tree that accounts for vessel and particle diameter distributions and show the non-uniform microsphere deposition pattern that results. Potential applications include investigating dose heterogeneity and microsphere deposition patterns.

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.

Modeling sphere dynamics in blood vessels for SIRT pre-planning - To fathom the potential and limitations

For selective internal radiation therapy (SIRT) the calculation of the 3D distribution of spheres based on individual blood flow properties is still an open and relevant research question. The purpose of this work is to develop and analyze a new treatment planning method for SIRT to calculate the absorbed dose distribution. For this intention, flow dynamics of the SIRT-spheres inside the blood vessels was simulated. The challenge is treatment planning solely using high-resolution imaging data available before treatment. The resolution required to reliably predict the sphere distribution and hence the dose was investigated. For this purpose, arteries of the liver were segmented from a contrast-enhanced angiographic CT. Due to the limited resolution of the given CT, smaller vessels were generated via a vessel model. A combined 1D/3D-flow simulation model was implemented to simulate the final 3D distribution of spheres and dose. Results were evaluated against experimental data from Y90-PET. Analysis showed that the resolution of the vessels within the angiographic CT of about 0.5mm should be improved to a limit of about 150μm to reach a reliable prediction.