Heterogeneity of microsphere distribution in resected liver and tumour tissue following selective intrahepatic radiotherapy

Spatial distribution of fractionally administered holmium microspheres in non-tumorous human liver tissue: how livers survive transarterial radioembolisation

BACKGROUND

Relatively high mean absorbed doses to the non-tumorous liver tissue (NTLT) are generally well tolerated in transarterial radioembolisation (TARE), potentially due to a heterogeneous dose distribution. This study investigates the macroscopic and microscopic distribution of fractionally administered TARE holmium microspheres in NTLT using an experimental setup of ex vivo perfused human donor livers under magnetic resonance imaging (MRI), and validates these findings through a comparison with MRI data from TARE-treated patients.

RESULTS

MRI-based dose maps of the TARE-treated ex vivo livers and patients revealed a heterogeneous dose distribution pattern throughout the NTLT (heterogeneity index (HI) range 2.96-10.11). Microscopic analysis confirmed this, as a wide variation in the percentage of tissue within 2.1 mm of microspheres (5.4%-84.3%) was observed. Microspheres administered in consecutive fractions decreased the heterogeneity, which was observed macroscopically by a decreased HI, and microscopically by the formation of new microsphere clusters. However, this HI decrease appeared finite, and new clusters formed near existing clusters, maintaining the overall distribution pattern.

CONCLUSIONS

TARE induces a heterogeneous dose distribution pattern in human NTLT. This heterogeneous dose distribution pattern persists across additional microsphere fractions, leaving parts of the NTLT unexposed to lethal doses of ionising radiation. Combined with the regenerative capacity of the liver, this may explain why relatively high mean absorbed doses to the NTLT are generally well tolerated in TARE.

REGISTRATION

For validation purposes, clinical data from patients who participated in a previous study (ClinicalTrials.gov, identifier NCT04269499, registered on February 13, 2020) was analysed in the current study.

Sequencing microsphere selective internal radiotherapy after external beam radiotherapy for hepatocellular carcinoma: proof of concept of a synergistic combination

OBJECTIVES

This study aims to explore the synergistic effects of combining stereotactic body radiation therapy (SBRT) and selective internal radiation therapy (SIRT) in that specific sequence for treating hepatocellular carcinoma (HCC), particularly in patients at high risk of radiation-induced liver disease (RILD).

METHODS

We analysed a case of a patient with HCC who was treated with SBRT at our institution. A virtual 90Y dose distribution was added using our in-house MIDOS model to keep a minimum dose to the healthy liver tissue. BED and EUD metrics were calculated to harmonize the dose distributions of SBRT and SIRT.

RESULTS

Our radiation biology-based models suggest that the combination of SBRT and SIRT could maintain effective tumour control while reducing the dose to normal liver tissue. Specifically, an SBRT plan of 10 Gy×3 fractions combined with SIRT yielded comparable tumour control probability to an SBRT-only plan of 10 Gy×5 fractions.

CONCLUSIONS

The combination of SBRT and SIRT is a promising treatment strategy for HCC patients at high risk of RILD. While the LQ model and associated formalisms provide a useful starting point, further studies are needed to account for factors beyond these models. Nonetheless, the potential for significant dose reduction to normal tissue suggests that this combination therapy could offer substantial clinical benefits.

ADVANCES IN KNOWLEDGE

This article presents a proposal to combine SBRT and SIRT, in this specific order, for HCC, discussing its advantages. A framework for future research into optimizing combination therapy for HCC is provided, utilizing a novel HCC vascular model to simulate 90Y doses.

Imageable Radioembolization Microspheres for Treatment of Unresectable Hepatocellular Carcinoma: Interim Results from a First-in-Human Trial

PURPOSE

To determine 6-month interim safety, effectiveness, and multimodal imageability of imageable glass microsphere yttrium-90 (90Y) radioembolization for unresectable hepatocellular carcinoma (HCC) in a first-in-human trial.

MATERIALS AND METHODS

Imageable microspheres (Eye90 Microspheres; ABK Biomedical, Halifax, Nova Scotia, Canada), a U.S. Food and Drug Administration (FDA) Breakthrough-Designated Device consisting of glass radiopaque 90Y microspheres visible on computed tomography (CT) and single photon emission CT (SPECT), were used to treat 6 subjects with unresectable HCC. Patients underwent selective (≤2 segments) treatment in a prospective open-label pilot trial. Key inclusion criteria included liver-only HCC, performance status ≤1, total lesion diameter ≤9 cm, and Child-Pugh A status. Prospective partition dosimetry was utilized. Safety (measured by Common Terminology Criteria for Adverse Events [CTCAE] v5), multimodal imageability on CT and SPECT, and 3- and 6-month imaging response by modified Response Evaluation Criteria in Solid Tumors on magnetic resonance (MR) imaging were evaluated.

RESULTS

Seven tumors in 6 subjects were treated and followed to 180 days. Administration success was 100%. Microsphere distribution measured by radiopacity on CT correlated with SPECT. Ninety-day target lesion complete response (CR) was observed in 3 of 6 subjects (50%) and partial response (PR) in 2 (33.3%). At 180 days, target lesion CR was maintained in 3 subjects (50%) and PR in 1 (16.7%). Two subjects could not be reassessed, having undergone intervening chemoembolization. All subjects reported adverse events (AEs), and 5 reported AEs related to treatment. There were no treatment-related Grade ≥3 AEs.

CONCLUSIONS

Radioembolization using imageable microspheres was safe and effective in 6 subjects with unresectable HCC at 6-month interim analysis. Microsphere distribution by radiopacity on CT correlated with radioactivity distribution by SPECT, providing previously unavailable CT-based tumor targeting information.

Transcatheter pseudo-vascular isolation for localization and concentration of a large molecule theranostic probe into a transgenic OncoPIG kidney tumor

INTRODUCTION

Great strides have been made identifying molecular and genetic changes expressed by various tumor types. These molecular and genetic changes are used as pharmacologic targets for precision treatment using large molecule (LM) proteins with high specificity. Theranostics exploits these LM biomolecules via radiochemistry, creating sensitive diagnostic and therapeutic agents. Intravenous (i.v.) LM drugs have an extended biopharmaceutical half-life thus resulting in an insufficient therapeutic index, permitting only palliative brachytherapy due to unacceptably high rates of systemic nontarget radiation doses to normal tissue. We employ tumor arteriole embolization isolating a tumor from the systemic circulation, and local intra-arterial (i.a.) infusion to improve uptake of a LM drug within a porcine renal tumor (RT).

METHODS

In an oncopig RT we assess the in vivo biodistribution of 99mTc-labeled macroaggregated albumin (MAA) a surrogate for a LM theranostics agent in the RT, kidney, liver, spleen, muscle, blood, and urine. Control animals underwent i.v. infusion and experimental group undergoing arteriography with pseudovascular isolation (PVI) followed by direct i.a. injection.

RESULTS

Injected dose per gram (%ID/g) of the LM at 1 min was 86.75 ± 3.76 and remained elevated up to 120 min (89.35 ± 5.77) with i.a. PVI, this increase was statistically significant (SS) compared to i.v. (13.38 ± 1.56 and 12.02 ± 1.05; p = 0.0003 p = 0.0006 at 1 and 120 min respectively). The circulating distribution of LM in the blood was less with i.a. vs i.v. infusion (2.28 ± 0.31 vs 25.17 ± 1.84 for i.v. p = 0.033 at 1 min). Other organs displayed a trend towards less exposure to radiation for i.a. with PVI compared to i.v. which was not SS.

CONCLUSION

PVI followed by i.a. infusion of a LM drug has the potential to significantly increase the first pass uptake within a tumor. This minimally invasive technique can be translated into clinical practice, potentially rendering monoclonal antibody based radioimmunotherapy a viable treatment for renal tumors.

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.

Concepts and methods for the dosimetry of radioembolisation of the liver with Y-90-loaded microspheres

This article aims at presenting in a didactic way, dosimetry concepts and methods that are relevant for radio-embolization of the liver with 90Y-microspheres. The application of the medical internal radiation dose formalism to radio-embolization is introduced. This formalism enables a simplified dosimetry, where the absorbed dose in a given tissue depends on only its mass and initial activity. This is applied in the single-compartment method, partition model, for the liver, tumour and lung dosimetry, and multi-compartment method, allowing identification of multiple tumours. Voxel-based dosimetry approaches are also discussed. This allows taking into account the non-uniform uptake within a compartment, which translates into a non-uniform dose distribution, represented as a dose-volume histogram. For this purpose, dose-kernel convolution allows propagating the energy deposition around voxel-sources in a computationally efficient manner. Alternatively, local-energy deposition is preferable when the spatial resolution is comparable or larger than the beta-particle path. Statistical tools may be relevant in establishing dose-effect relationships in a given population. These include tools such as the logistic regression or receiver operator characteristic analysis. Examples are given for illustration purpose. Moreover, tumour control probability modelling can be assessed through the linear-quadratic model of Lea and Catcheside and its counterpart, the normal-tissue complication probability model of Lyman, which is suitable to the parallel structure of the liver. The selectivity of microsphere administration allows tissue sparing, which can be considered with the concept of equivalent uniform dose, for which examples are also given. The implication of microscopic deposition of microspheres is also illustrated through a liver toxicity model, even though it is not clinically validated. Finally, we propose a reflection around the concept of therapeutic index (TI), which could help tailor treatment planning by determining the treatment safety through the evaluation of TI based on treatment-specific parameters.

Precision dosimetry in yttrium-90 radioembolization through CT imaging of radiopaque microspheres in a rabbit liver model

PURPOSE

To perform precision dosimetry in yttrium-90 radioembolization through CT imaging of radiopaque microspheres in a rabbit liver model and to compare extracted dose metrics to those produced from conventional PET-based dosimetry.

MATERIALS AND METHODS

A CT calibration phantom was designed containing posts with nominal microsphere concentrations of 0.5 mg/mL, 5.0 mg/mL, and 25.0 mg/mL. The mean Hounsfield unit was extracted from the post volumes to generate a calibration curve to relate Hounsfield units to microsphere concentration. A nominal bolus of 40 mg of microspheres was administered to the livers of eight rabbits, followed by PET/CT imaging. A CT-based activity distribution was calculated through the application of the calibration curve to the CT liver volume. Post-treatment dosimetry was performed through the convolution of yttrium-90 dose-voxel kernels and the PET- and CT-based cumulated activity distributions. The mean dose to the liver in PET- and CT-based dose distributions was compared through linear regression, ANOVA, and Bland-Altman analysis.

RESULTS

A linear least-squares fit to the average Hounsfield unit and microsphere concentration data from the calibration phantom confirmed a strong correlation (r2 > 0.999) with a slope of 14.13 HU/mg/mL. A poor correlation was found between the mean dose derived from CT and PET (r2 = 0.374), while the ANOVA analysis revealed statistically significant differences (p < 10-12) between the MIRD-derived mean dose and the PET- and CT-derived mean dose. Bland-Altman analysis predicted an offset of 15.0 Gy between the mean dose in CT and PET. The dose within the liver was shown to be more heterogeneous in CT than in PET with an average coefficient of variation equal to 1.99 and 1.02, respectively.

CONCLUSION

The benefits of a CT-based approach to post-treatment dosimetry in yttrium-90 radioembolization include improved visualization of the dose distribution, reduced partial volume effects, a better representation of dose heterogeneity, and the mitigation of respiratory motion effects. Post-treatment CT imaging of radiopaque microspheres in yttrium-90 radioembolization provides the means to perform precision dosimetry and extract accurate dose metrics used to refine the understanding of the dose-response relationship, which could ultimately improve future patient outcomes.

Yttrium-90 TOF-PET-Based EUD Predicts Response Post Liver Radioembolizations Using Recommended Manufacturer FDG Reconstruction Parameters

Purpose

Explaining why ⁹⁰Y TOF-PET based equivalent uniform dose (EUD) using recommended manufacturer FDG reconstruction parameters has been shown to predict response.

Methods

The hot rods insert of a Jaszczak deluxe phantom was partially filled with a 2.65 GBq ⁹⁰Y - 300ml DTPA water solution resulting in a 100 Gy mean absorbed dose in the 6 sectors. A two bed 20min/position acquisition was performed on a 550ps- and on a 320ps- TOF-PET/CT and reconstructed with recommended manufacturer FDG reconstruction parameters, without and with additional filtering. The whole procedure was repeated on both PET after adding 300ml of water (50Gy setup). The phantom was acquired again after decay by a factor of 10 (5Gy setup), but with 200min per bed position. For comparison, the phantom was also acquired with ¹⁸F activity corresponding to a clinical FDG whole body acquisition.

Results

The 100Gy-setup provided a hot rod sectors image almost as good as the ¹⁸F phantom. However, despite acquisition time compensation, the 5Gy-setup provides much lower quality imaging. TOF-PET based sectors EUDs for the three large rod sectors agreed with the actual EUDs computed with a radiosensitivity of 0.021Gy^-1 well in the range observed in external beam radiotherapy (EBRT), i.e. 0.01-0.04Gy^-1. This agreement explains the reunification of the dose-response relationships of the glass and resin spheres in HCC using the TOF-PET based EUD. Additional filtering reduced the EUDs agreement quality.

Conclusions

Recommended manufacturer FDG reconstruction parameters are suitable in TOF-PET post ⁹⁰Y liver radioembolization for accurate tumour EUD computation. The present results rule out the use of low specific activity phantom studies to optimize reconstruction parameters.

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.

Correlation and Agreement of Yttrium-90 Positron Emission Tomography/Computed Tomography with Ex Vivo Radioembolization Microsphere Deposition in the Rabbit VX2 Liver Tumor Model

PURPOSE

To demonstrate a stronger correlation and agreement of yttrium-90 (⁹⁰Y) positron emission tomography (PET)/computed tomography (CT) measurements with explant liver tumor dosing compared with the standard model (SM) for radioembolization.

MATERIALS AND METHODS

Hepatic VX2 tumors were implanted into New Zealand white rabbits, with growth confirmed by 7 T magnetic resonance imaging. Seventeen VX2 rabbits provided 33 analyzed tumors. Treatment volumes were calculated from manually drawn volumes of interest (VOI) with three-dimensional surface renderings. Radioembolization was performed with glass ⁹⁰Y microspheres. PET/CT imaging was completed with scatter and attenuation correction. Three-dimensional ellipsoid VOI were drawn to encompass tumors on fused images. Tumors and livers were then explanted for inductively coupled plasma (ICP)-optical emission spectroscopy (OES) analysis of microsphere content. ⁹⁰Y PET/CT and SM measurements were compared with reference standard ICP-OES measurements of tumor dosing with Pearson correlation and Bland-Altman analyses for agreement testing with and without adjustment for tumor necrosis.

RESULTS

The median infused activity was 33.3 MBq (range, 5.9-152.9). Tumor dose was significantly correlated with ⁹⁰Y PET/CT measurements (r = 0.903, P < .001) and SM estimates (r = 0.607, P < .001). Bland-Altman analyses showed that the SM tended to underestimate the tumor dosing by a mean of -8.5 Gy (CI, -26.3-9.3), and the degree of underestimation increased to a mean of -18.3 Gy (CI, -38.5-1.9) after the adjustment for tumor necrosis.

CONCLUSIONS

⁹⁰Y PET/CT estimates were strongly correlated and had better agreement with reference measurements of tumor dosing than SM estimates.

Personalized Dosimetry for Liver Cancer Y-90 Radioembolization Using Computational Fluid Dynamics and Monte Carlo Simulation

Yttrium-90 (Y-90) transarterial radioembolization uses radioactive microspheres injected into the hepatic artery to irradiate liver tumors internally. One of the major challenges is the lack of reliable dosimetry methods for dose prediction and dose verification. We present a patient-specific dosimetry approach for personalized treatment planning based on computational fluid dynamics (CFD) simulations of the microsphere transport combined with Y-90 physics modeling called CFDose. The ultimate goal is the development of a software to optimize the amount of activity and injection point for optimal tumor targeting. We present the proof-of-concept of a CFD dosimetry tool based on a patient's angiogram performed in standard-of-care planning. The hepatic arterial tree of the patient was segmented from the cone-beam CT (CBCT) to predict the microsphere transport using multiscale CFD modeling. To calculate the dose distribution, the predicted microsphere distribution was convolved with a Y-90 dose point kernel. Vessels as small as 0.45 mm were segmented, the microsphere distribution between the liver segments using flow analysis was predicted, the volumetric microsphere and resulting dose distribution in the liver volume were computed. The patient was imaged with positron emission tomography (PET) 2 h after radioembolization to evaluate the Y-90 distribution. The dose distribution was found to be consistent with the Y-90 PET images. These results demonstrate the feasibility of developing a complete framework for personalized Y-90 microsphere simulation and dosimetry using patient-specific input parameters.

Significant artefactual noise in 90Y TOF-PET imaging of low specific activity phantoms arises despite increased acquisition time

Volumes of usual PET phantoms are about four to sixfold that of a human liver. In order to avoid count rate saturation and handling of very high ⁹⁰Y activity, reported TOF-PET phantom studies are performed using specific activities lower than those observed in liver radioembolization.

However, due to the constant random coincidence rate induced by the natural crystal radioactivity, reduction of ⁹⁰Y specific activity in TOF-PET imaging cannot be counterbalanced by increasing the acquisition time. As a result, most ⁹⁰Y phantom studies reported images noisier than those obtained in whole-body ¹⁸F-FDG, and thus advised to use dedicated noise control in TOF-PET imaging post ⁹⁰Y liver radioembolization.

We performed acquisitions of the Jaszczak Deluxe phantom in which the hot rod insert was only partially filled with 2.6 GBq of ⁹⁰Y. Standard reconstruction parameters recommended by the manufacturer for whole-body ¹⁸F-FDG PET were used.

Low specific activity setups, although exactly compensated by increasing the acquisition time in order to get the same number of detected true coincidences per millilitre, were impacted by significant noise. On the other hand, specific activity and acquisition time setup similar to that used in post ⁹⁰Y liver radioembolization provided image quality very close to that of whole-body ¹⁸F-FDG.

This result clearly discards the use of low specific activity phantoms intended to TOF-PET reconstruction parameter optimization. Volume reduction of large phantoms can be achieved by vertically setting the phantoms or by adding Styrofoam inserts.

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.

90Y TOF-PET based EUD reunifies patient survival prediction in resin and glass microspheres radioembolization of HCC tumours

Clinical studies reported a twofold ratio between the efficacies per Gy of resin versus glass spheres. Our aim is to investigate whether this difference could result from the different degrees of heterogeneity in sphere distribution between the two medical devices. The ⁹⁰Y TOF-PET based equivalent uniform doses (EUD) was used for this purpose. 58 consecutive HCC radioembolizations were retrospectively analyzed. Absorbed doses D and Jones-Hoban EUD in lesions were computed. Radioembolization efficacy was assessed using Kaplan-Meier survival curves. In order to match together the glass and resin spheres survival curves using a 40 Gy-threshold, an efficacy factor of 0.73 and 0.36 has to be applied on their absorbed dose, respectively. Using EUD, a nice matching between glass and resin survival curves was obtained with a better separation of the responding and not responding survival curves. The results clearly support the fact that the activity heterogeneity observed in ⁹⁰Y TOF-PET post radioembolization does not only result from statistical noise, but also reflects the actual heterogeneity of the spheres distribution. Use of EUD reunifies the efficacy of the two medical devices.

Correlation of radiation dose and activity with clinical outcomes in metastatic colorectal cancer after selective internal radiation therapy using yttrium-90 resin microspheres

PURPOSE

Yttrium-90 (Y)-resin microspheres are prescribed using activity. We evaluated overall survival (OS) and radiographic tumor response after selective internal radiation therapy (SIRT) with resin microspheres in patients with liver metastases from colorectal cancer.

PATIENTS AND METHODS

We retrospectively reviewed 60 metastatic colorectal cancer patients treated at our institution with SIRT using Y-resin microspheres. Each patient underwent pre-SIRT MRI or computed tomography imaging of the liver with intravenous contrast. Patients underwent post-treatment imaging at 2-3-month intervals with response assessed according to unidimensional Response Evaluation Criteria in Solid Tumors (RECIST) criteria as well as published three-dimensional volumetric criteria. We then related the prescribed activity established by the body surface area method and the corresponding prescribed dose to radiographic treatment response and OS.

RESULTS

The median follow-up after the first SIRT treatment was 8.9 months. The mean prescribed activity and the prescribed dose were 26.6 mCi and 52.8 Gy, respectively. OS was not significantly associated with either prescribed activity or prescribed dose. Prescribed dose was also not related to response. However, a significant relationship was found between a higher prescribed activity and an improved radiographic response by RECIST (P=0.04) at the second follow-up.

CONCLUSION

The prescribed activity of Y-resin microspheres may be correlated with radiographic response by RECIST criteria at 4-6 months post-treatment. For a more accurate prediction of response, a valid dose calculation model based on post-Y PET dosimetry is likely needed given the heterogeneous dose delivery seen in SIRT.

Quantitative biomarkers for liver metastases: comparison of MRI diffusion-weighted imaging heterogeneity index and fluorine-18-fluoro-deoxyglucose standardised uptake value in hybrid PET/MR

AIM

To investigate the ability of apparent diffusion coefficient (ADC) heterogeneity index to discriminate liver metastases (LM) from normal-appearing liver (NAL) tissue as compared to common magnetic resonance imaging (MRI) metrics, and to investigate its correlation with 2-[¹⁸F]-fluoro-2-deoxy-d-glucose (¹⁸F-FDG) positron-emission tomography (PET) standardised uptake value (SUV).

MATERIALS AND METHODS

Thirty-nine liver metastases in 24 oncology patients (13 women, 11 men; mean age 56±13 years) with proven LM from heterogeneous sources were evaluated on a PET/MRI system. Abdominal sequences included Dixon and diffusion-weighted imaging (DWI) protocols with simultaneous PET. Tissue heterogeneity was calculated using the coefficient of variance (CV) of the ADC, and compared in LM and in NAL tissue of the same volume in an adjacent portion of the liver. The correlations between various ADC measures and PET SUV in distinguishing LM from NAL were evaluated.

RESULTS

A good correlation was found between ADCcv and SUVpeak (r=0.712). Moderate inverse correlation was found between ADCmin and SUVpeak (r=-0.536), and a weak inverse correlation between ADCmean and SUVpeak (r=-0.273). There was a significant difference between LM and NAL when ADCcv (p<0.0001) and ADCmin (p=0.001) were used. Receiver operating characteristic (ROC) analysis of SUV, ADCcv, ADCmin, and ADCmean produced an AUC of 0.989, 0.900, 0.742, and 0.623 respectively.

CONCLUSIONS

The ADCcv index is a potential biomarker of LM with better correlation to ¹⁸F-FDG PET SUVpeak than conventional MRI metrics, and may serve to quantitatively discriminate between LM and NAL.

Autoradiography and biopsy measurements of a resected hepatocellular carcinoma treated with 90 yttrium radioembolization demonstrate large absorbed dose heterogeneities

Purpose

Radioembolization is an alternative palliative treatment for hepatocellular carcinoma. Here, we examine the uptake differences between tumor tissue phenotypes and present a cross-section of the absorbed dose throughout a liver tissue specimen.

Methods and materials

A patient with hepatocellular carcinoma was treated with ⁹⁰Y radioembolization followed by liver tissue resection. Gamma camera images and autoradiographs were collected and biopsy tissue samples were analyzed using a gamma well counter and light microscopy.

Results

An analysis of 25 punched biopsy tissue samples identified 4 tissue regions: Normal tissue, viable tumor tissue with and without infarcted areas, and tumor areas with postnecrotic scar tissue. Autoradiography and biopsy tissue sample measurements showed large dose differences between viable and postnecrotic tumor tissue (159 Gy vs 23 Gy).

Conclusions

Radioembolization of 90 yttrium with resin microspheres produces heterogeneous-absorbed dose distributions in the treatment of unifocal hepatic malignancies that could not be accurately determined with current gamma camera imaging techniques.

Clinical and imaging-based prognostic factors in radioembolisation of liver metastases from colorectal cancer: a retrospective exploratory analysis

Background

The aim of this study was to investigate the relationship between absorbed dose and response of colorectal cancer liver metastases treated with [⁹⁰Y]-resin microspheres and to explore possible clinical and imaging derived prognostic factors.

Methods

FDG PET/CT was used to measure response of individual lesions to a measured absorbed dose, derived from post-treatment ⁹⁰Y PET imaging. Predicted dose was also derived from planning [^99mTc]-MAA SPECT data. Peak standardised uptake value and total lesion glycolysis (TLG) were explored as response measures, and compared to dose metrics including average dose (D_avg), biologically effective dose, minimum dose to 70% of lesion volume and volume receiving at least 50 Gy. Prognostic factors examined included baseline TLG, RAS mutation status, FDG heterogeneity and dose heterogeneity. In an exploratory analysis, response and clinico-pathological variables were evaluated and compared to overall survival.

Results

Sixty-three lesions were analysed from 22 patients. Poor agreement was seen between predicted and measured dose values. TLG was a superior measure of response, and all dose metrics were significant prognostic factors, with a D_avg of ~50 Gy derived as the critical threshold for a significant response (>50% reduction in TLG). No significant correlation was found between baseline TLG or RAS mutation status and response. Measured dose heterogeneity was a significant prognostic factor and when combined with D_avg had a positive predictive value for response >80%. In the exploratory analysis for prognostic factors of survival, low hepatic tumour burden and mean reduction in TLG >65% were independently associated with improved overall survival.

Conclusions

Lesions receiving an average dose greater than 50 Gy are likely to have a significant response. For lesions receiving less than 50 Gy, dose heterogeneity is a significant prognostic factor. Lesions receiving an average dose less than 20 Gy are unlikely to respond. A reduction in TLG may be associated with improved overall survival.

Electronic supplementary material

The online version of this article (doi:10.1186/s13550-017-0292-1) contains supplementary material, which is available to authorized users.

Computational simulation of the predicted dosimetric impact of adjuvant yttrium-90 PET/CT-guided percutaneous ablation following radioembolization

BACKGROUND

90Y PET/CT post-radioembolization imaging has demonstrated that the distribution of 90Y in a tumor can be non-uniform. Using computational modeling, we predicted the dosimetric impact of post-treatment 90Y PET/CT-guided percutaneous ablation of the portions of a tumor receiving the lowest absorbed dose. A cohort of fourteen patients with non-resectable liver cancer previously treated using 90Y radioembolization were included in this retrospective study. Each patient exhibited potentially under-treated areas of tumor following treatment based on quantitative 90Y PET/CT. 90Y PET/CT was used to guide electrode placement for simulated adjuvant radiofrequency ablation in areas of tumor receiving the lowest dose. The finite element method was used to solve Penne's bioheat transport equation, coupled with the Arrhenius thermal cell-death model to determine 3D thermal ablation zones. Tumor and unablated tumor absorbed-dose metrics (average dose, D50, D70, D90, V100) following ablation were compared, where D70 is the minimum dose to 70% of tumor and V100 is the fractional tumor volume receiving more than 100 Gy.

RESULTS

Compared to radioembolization alone, 90Y radioembolization with adjuvant ablation was associated with predicted increases in all tumor dose metrics evaluated. The mean average absorbed dose increased by 11.2 ± 6.9 Gy. Increases in D50, D70, and D90 were 11.0 ± 6.9 Gy, 13.3 ± 10.9 Gy, and 11.8 ± 10.8 Gy, respectively. The mean increase in V100 was 7.2 ± 4.2%. All changes were statistically significant (P < 0.01). A negative correlation between pre-ablation tumor volume and D50, average dose, and V100 was identified (ρ < - 0.5, P < 0.05) suggesting that adjuvant radiofrequency ablation may be less beneficial to patients with large tumor burdens.

CONCLUSIONS

This study has demonstrated that adjuvant 90Y PET/CT-guided radiofrequency ablation may improve tumor absorbed-dose metrics. These data may justify a prospective clinical trial to further evaluate this hybrid approach.

A Microdosimetric Analysis of Absorbed Dose to Tumor as a Function of Number of Microspheres per Unit Volume in 90Y Radioembolization

Differences in maximum tolerable absorbed dose to normal liver between (90)Y radioembolization and external-beam radiation therapy have been explained by citing differences in absorbed-dose heterogeneity at the microscopic level. We investigated microscopic absorbed-dose heterogeneity in radioembolization as a function of the number of microspheres per unit volume in tumor. The goal was to determine what effect the number of microspheres may have, if any, on tumor control in (90)Y radioembolization.

METHODS

(90)Y PET/CT data were combined with microscopic probability-density functions describing microsphere clustering to provide realistic simulation using Monte Carlo modeling on both a macroscopic and a microscopic level. A complete microdosimetric analysis using 100-μm voxels was performed on the basis of (90)Y PET/CT data from 19 tumors treated using radioembolization. Simulations were performed with average tumor microsphere-number densities from 200 to 70,000 spheres/mL. Monte Carlo simulations of each tumor and number density were repeated 20 times to establish SE. A 2-way balanced ANOVA was used to determine whether differences in microsphere-number density affected common tumor-dose metrics.

RESULTS

Decreasing the microsphere-number density resulted in a decrease in D70, the minimum dose to 70% of the tumor. The slope of the dose-volume histogram also decreased with decreasing microsphere-number density in all tumors. Compared with a density of 50,000 spheres/mL, decreases in D70 were statistically significant below 20,000 spheres/mL. However, these differences are unlikely to have clinical significance until the density decreases to below 5,000 spheres/mL. Although D70 was decreased at a low microsphere-number density, one can compensate for decreases by an increase in the average tumor-absorbed dose, that is, by increasing the radioembolization treatment dose.

CONCLUSION

Differences in microsphere-number density may have an effect on microscopic tumor absorbed-dose inhomogeneity. These results begin to explain differences in treatment planning strategies between glass and resin radioembolization devices.

A review of 3D image-based dosimetry, technical considerations and emerging perspectives in 90Y microsphere therapy

Yttrium-90 radioembolization (⁹⁰Y-RE) is a well-established therapy for the treatment of hepatocellular carcinoma (HCC) and also of metastatic liver deposits from other malignancies. Nuclear Medicine and Cath Lab diagnostic imaging takes a pivotal role in the success of the treatment, and in order to fully exploit the efficacy of the technique and provide reliable quantitative dosimetry that are related to clinical endpoints in the era of personalized medicine, technical challenges in imaging need to be overcome. In this paper, the extensive literature of current ⁹⁰Y-RE techniques and challenges facing it in terms of quantification and dosimetry are reviewed, with a focus on the current generation of 3D dosimetry techniques. Finally, new emerging techniques are reviewed which seek to overcome these challenges, such as high-resolution imaging, novel surgical procedures and the use of other radiopharmaceuticals for therapy and pre-therapeutic planning.

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

Background

The higher tolerated mean absorbed dose for selective internal radiation therapy (SIRT) with intra-arterially infused ⁹⁰Y microspheres compared to external beam therapy is speculated to be caused by absorbed dose inhomogeneity, which allows for liver regeneration. However, the complex liver microanatomy and rheology makes modelling less valuable if the tolerance doses are not based on the actual microsphere distribution. The present study demonstrates the sphere distribution and small-scale absorbed dose inhomogeneity and its correlation with the mean absorbed dose in liver tissue resected after SIRT.

Methods

A patient with marginally resectable cholangiocarcinoma underwent SIRT 9 days prior to resection including adjacent normal liver tissue. The resected specimen was formalin-fixed and sliced into 1 to 2-mm sections. Forty-one normal liver biopsies 6-8 mm in diameter were punched from these sections and the radioactivity measured. Sixteen biopsies were further processed for detailed analyses by consecutive serial sectioning of 15 30-μm sections per biopsy, mounted and stained with haematoxylin-eosin. All sections were scrutinised for isolated or conglomerate spheres. Small-scale dose distributions were obtained by applying a ⁹⁰Y-dose point kernel to the microsphere distributions.

Results

A total of 3888 spheres were found in the 240 sections. Clusters were frequently found as strings in the arterioles and as conglomerates in small arteries, with the largest cluster comprising 453 spheres. An increased mean absorbed dose in the punch biopsies correlated with large clusters and a greater coefficient of variation. In simulations the absorbed dose was 5–1240 Gy; 90% were 10-97 Gy and 45% were <30 Gy, the assumed tolerance in external beam therapy.

Conclusions

Sphere clusters were located in both arterioles and small arteries and increased in size with increasing sphere concentration, resulting in increased absorbed dose inhomogeneity, which contradicts earlier modelling studies.