Impact of image reconstruction method on dose distributions derived from 90Y PET images: phantom and liver radioembolization patient studies

Voxel-Based Dosimetry as a Means for Treatment Personalisation in Radioembolization: A Systematic Review

INTRODUCTION

Radionuclide therapy including 90Y radioembolization is an established form of brachytherapy for treatment of malignancy including hepatocellular carcinoma. Currently, there are several methods available to estimate patient absorbed dose, including voxel-based dosimetry, that can achieve a level of personalisation in the planning and outcome assessments of radioembolization. Despite the advantages of voxel-based dosimetry, it remains a relatively new concept in radioembolization. This study evaluates if voxel-based dosimetry was associated with improved treatment efficacy in radioembolization planning.

METHODS

A systematic review was conducted by searching relevant databases (Medline Ovid, PubMed, Embase Ovid, CINAHL Complete, Cochrane Library, CENTRAL, Australian New Zealand Clinical Trials Registry, ClinicalTrials.gov, WHO International Trials Registry, Google Scholar) for literature regarding voxel-based dosimetry in radioembolization.

RESULTS

A total of 41 papers were included for this systematic review. Review of these studies revealed that voxel-based dosimetry can benefit numerous aspects of radioembolization in radionuclide therapy including predicting tumour response, toxicity and patient survival. Numerous studies also indicated that voxel-based dosimetry in radioembolization is a more accurate approach in establishing a dose-effect relationship in targeted radionuclide therapy when compared to other methods. Despite these promising findings, these studies did not investigate or comment on the accuracy of voxel-based dosimetry.

CONCLUSION

The evidence from this review highlights that voxel-based dosimetry can improve treatment efficacy in radioembolization planning. However, further studies are required to validate the accuracy and feasibility of voxel-based dosimetry in clinical practice.

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

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

Contralateral Hypertrophy Post Yttrium-90 Transarterial Radioembolization in Patients With Hepatocellular Carcinoma and Portal Vein Tumor Thrombus

Objectives Contralateral hypertrophy of non-irradiated liver following Yttrium-90 (90Y) transarterial radioembolization (TARE) is increasingly recognized as an option to facilitate curative surgical resection in patients that would otherwise not be surgical candidates due to a small future liver remnant (FLR). This study aimed to investigate the correlation between patient features and liver hypertrophy and identify potential predictors for liver growth in patients with hepatocellular carcinoma (HCC) and portal vein tumor thrombus (PVTT) undergoing TARE. Methodology Twenty-three patients with HCC and PVTT were included. Contralateral liver hypertrophy was assessed at six months posttreatment based on CT or MRI imaging. Thirteen patient features were selected for statistical and prediction analysis. Univariate Spearman correlation and analysis of variance (ANOVA) tests were performed. Subsequently, four feature-selection methods based on multivariate analysis were used to improve model generalization performance. The selected features were applied to train linear regression models, with fivefold cross-validation to assess the performance of the predicted models. Results The ratio of disease-free target liver volume to spared liver volume and total liver volume showed the highest correlations with contralateral hypertrophy (P-values = 0.03 and 0.05, respectively). In three out of four feature-selection methods, the feature of disease-free target liver volume to total liver volume ratio was selected, having positive correlations with the outcome and suggesting that more hypertrophy may be expected when more volume of disease-free liver is irradiated. Conclusions Contralateral hypertrophy post-90Y TARE can be an option for facilitating surgical resection in patients with otherwise small FLR.

A single centre intercomparison between commercial treatment planning systems for 90Y radioembolization using virtual and experimental phantoms

INTRODUCTION

Dedicated Treatment Planning Systems (TPSs) were developed to personalize 90Y-transarterial radioembolization. This study evaluated the agreement among four commercial TPSs assessing volumes of interest (VOIs) volumes and dose metrics.

METHODS

A homogeneous (EH) and an anthropomorphic phantom with hot and cold inserts (EA) filled with 99mTc-pertechnetate were acquired with a SPECT/CT scanner. Their virtual versions (VH and VA, respectively) and a phantom with activity inside a single voxel (VK) were generated by an in-house MATLAB script. Images and delineated VOIs were imported into the TPSs to compute voxel-based absorbed dose distributions with various dose deposition approaches: local deposition method (LDM) and dose kernel convolution (DKC) with/without local density correction (LDC). VOI volumes and mean absorbed doses were assessed against their median value across TPSs. Dose-volume histograms (DVHs) and VK-derived dose profiles were evaluated.

RESULTS

Small (<2.1 %) and large (up to 42.4 %) relative volume differences were observed on large (>500 ml) and small VOIs, respectively. Mean absorbed doses relative differences were < 3 % except for small VOIs with steep dose gradients (up to 89.1 % in the VA Cold Sphere VOI). Within the same TPS, LDC negligibly affected the mean absorbed dose, while DKC and LDM showed differences up to 63 %. DHVs were mostly overlapped in experimental phantoms, with some differences in the virtual versions. Dose profiles agreed within 1 %.

CONCLUSION

TPSs showed an overall good agreement except for small VOI volumes and mean absorbed doses of VOIs with steep dose gradients. These discrepancies should be considered in the dosimetry uncertainty assessment, thus requiring an appropriate harmonization.

A novel figure of merit to investigate 68Ga PET/CT image quality based on patient weight and lesion size using Q.Clear reconstruction algorithm: A phantom study

INTRODUCTION

Q.Clear is a Bayesian penalised-likelihood algorithm that uses a β-value for positron emission tomography(PET)/computed tomography(CT) image reconstruction(IR). Our study proposes a novel figure of merit, named CRBV, to compare the Q.Clear performances using 68Ga PET/CT image with the ordered-subset-expectation-maximization(OSEM) algorithm and to identify the optimal β-values for these images using two phantoms mimicking normal and overweight patients.

METHODS

NEMA IQ phantom with or without a ring of water-filled plastic bags (NEMAstd and NEMAow, respectively) was acquired and reconstructed with OSEM and Q.Clear at various β-values and minutes/bed position(min/bp). Contrast recovery(CR), background variability(BV) and CRBV were calculated. Highest CRBV values were used to identify optimal β-value ranges.

RESULTS

Q.Clear with 250 ≤ β ≤ 800 improved CRBV compared to OSEM for all the investigated spheres and acquisition setups. Outside of this range, Q.Clear still outperformed OSEM with few exceptions depending on spheres diameters and phantoms(e.g.,β-value = 1600 for diameters ≤ 17 mm using the NEMAow phantom). Regarding the CRBV performance for IR optimization, for the 4 min/bp NEMAstd IR, β-values = 300 ÷ 350 allowed to simultaneously optimize all diameters(except for the 10 mm); for the NEMAow IR, β-values = 350 ÷ 500 were needed for diameters > 20 mm, while β-values = 200 ÷ 250 were selected for the remaining diameters. For the 2 min/bp, β-value = 500 was suitable for diameters > 17 mm in both NEMAstd and NEMAow IR, while for smaller diameters β-value = 200 and β-values = 250 ÷ 350 were obtained for NEMAstd and NEMAow, respectively.

CONCLUSION

Almost all tested β-values of Q.Clear improved the CRBV compared to OSEM. In both phantoms, simulating normal and over-weight patients, optimal β-values were found according to lesion sizes and investigated acquisition times.

Effects of Respiratory Motion on Y-90 PET Dosimetry for SIRT

Respiratory motion degrades the quantification accuracy of PET imaging by blurring the radioactivity distribution. In the case of post-SIRT PET-CT verification imaging, respiratory motion can lead to inaccuracies in dosimetric measures. Using an anthropomorphic phantom filled with ⁹⁰Y at a range of clinically relevant activities, together with a respiratory motion platform performing realistic motions (10–15 mm amplitude), we assessed the impact of respiratory motion on PET-derived post-SIRT dosimetry. Two PET scanners at two sites were included in the assessment. The phantom experiments showed that device-driven quiescent period respiratory motion correction improved the accuracy of the quantification with statistically significant increases in both the mean contrast recovery (+5%, p = 0.003) and the threshold activities corresponding to the dose to 80% of the volume of interest (+6%, p < 0.001). Although quiescent period gating also reduces the number of counts and hence increases the noise in the PET image, its use is encouraged where accurate quantification of the above metrics is desired.

Overview of commercial treatment planning systems for targeted radionuclide therapy

INTRODUCTION

Targeted Radionuclide Therapy (TRT) is a branch of cancer medicine dealing with the therapeutic use of radioisotopes associated with biological vectors accumulating in the tumors/targets, indicated as Molecular Radiotherapy (MRT), or directly injected into the arteries that supply blood to liver tumour vasculature, indicated as Selective RT (SRT). The aim of this work is to offer a panoramic view on the increasing number of commercially-available TRT treatment planning systems (TPSs).

MATERIALS AND METHODS

A questionnaire was sent to manufacturers' representatives. Academic software were not considered. Questions were grouped as follows: general information, clinical workflow, calibration procedure, image processing/reconstruction, image registration and segmentation tools, time-activity curve (TAC) fitting and absorbed dose calculation.

RESULTS

All software reported have CE-marking. TPSs were divided between SRT-dedicated software [4] and MRT [5] dosimetry software. In SRT, since no kinetic process is involved, absorbed dose calculation does not require TAC fitting, and image registration is not fully developed in all TPS. All software requires a radionuclide-specific calibration. In SRT, a relative image calibration can be obtained by scaling the counts to a known activity. Automated VOI contouring and rigid/deformable propagation between different acquisitions time-points is implemented in most TPSs, although DICOM export is rare. Different TAC fits are available depending on the number of time-points. Voxel S-value and Local deposition methods are the most frequent dosimetric approaches; dose-voxel kernel convolution and semi-Monte Carlo method are also available.

CONCLUSIONS

Available TPSs allows performing personalized dosimetry in clinical practice. Individual variations in methodology/algorithms must be considered in the standardisation/harmonization processes.

Yttrium-90 quantitative phantom study using digital photon counting PET

BACKGROUND

PET imaging of 90Y-microsphere distribution following radioembolisation is challenging due to the count-starved statistics from the low branching ratio of e+/e- pair production during 90Y decay. PET systems using silicon photo-multipliers have shown better 90Y image quality compared to conventional photo-multiplier tubes. The main goal of the present study was to evaluate reconstruction parameters for different phantom configurations and varying listmode acquisition lengths to improve quantitative accuracy in 90Y dosimetry, using digital photon counting PET/CT.

METHODS

Quantitative PET and dosimetry accuracy were evaluated using two uniform cylindrical phantoms specific for PET calibration validation. A third body phantom with a 9:1 hot sphere-to-background ratio was scanned at different activity concentrations of 90Y. Reconstructions were performed using OSEM algorithm with varying parameters. Time-of-flight and point-spread function modellings were included in all reconstructions. Absorbed dose calculations were carried out using voxel S-values convolution and were compared to reference Monte Carlo simulations. Dose-volume histograms and root-mean-square deviations were used to evaluate reconstruction parameter sets. Using listmode data, phantom and patient datasets were rebinned into various lengths of time to assess the influence of count statistics on the calculation of absorbed dose. Comparisons between the local energy deposition method and the absorbed dose calculations were performed.

RESULTS

Using a 2-mm full width at half maximum post-reconstruction Gaussian filter, the dosimetric accuracy was found to be similar to that found with no filter applied but also reduced noise. Larger filter sizes should not be used. An acquisition length of more than 10 min/bed reduces image noise but has no significant impact in the quantification of phantom or patient data for the digital photon counting PET. 3 iterations with 10 subsets were found suitable for large spheres whereas 1 iteration with 30 subsets could improve dosimetry for smaller spheres.

CONCLUSION

The best choice of the combination of iterations and subsets depends on the size of the spheres. However, one should be careful on this choice, depending on the imaging conditions and setup. This study can be useful in this choice for future studies for more accurate 90Y post-dosimetry using a digital photon counting PET/CT.