Assessment of acquisition protocols for routine imaging of Y-90 using PET/CT

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.

Reduction of emission time for [68Ga]Ga-PSMA PET/CT using the digital biograph vision: a phantom study

BACKGROUND

The aim of this phantom study was to optimize the [⁶⁸Ga]Ga-PSMA PET/CT examination in terms of scan time duration and image reconstruction parameters, in combination with PSF and TOF modelling, in a digital Biograph Vision PET/CT scanner.

METHODS

Three types of phantoms were used: 1) soft-tissue tumor phantom consisting of six spheres mounted in a torso phantom; 2) bone-lung tumor phantom; 3) resolution phantom. Phantom inserts were filled with activity concentrations (ACs) that were derived from clinical data. Phantom data were acquired in list-mode at one bed position. Images with emission data ranging from 30 to 210 s in 30-s increments were reconstructed from a reference image acquired with 3.5-min emission. Iterative image reconstruction (OSEM), point-spread-function (PSF) and time-of-flight (TOF) options were applied using different iterations, Gaussian filters, and voxel sizes. The criteria for image quality was lesion detectability and lesion quantification, evaluated as contrast-to-noise ratio (CNR) and maximum AC (peak AC), respectively. A threshold value of CNR above 6 and percentage maximum AC (peak AC) deviation range of ±20% of the reference image were considered acceptable. The proposed single-bed scan time reduction was projected to a whole-body examination (patient validation scan) using the continuous-bed-motion mode.

RESULTS

Sphere and background ACs of 20 kBq/mL and 1 kBq/mL were selected, respectively. The optimized single-bed scan time was approximately 60 s using OSEM-TOF or OSEM-TOF+PSF (four iterations, 4.0-mm Gaussian filter and almost isotropic voxel size of 3.0-mm side length), resulting in a PET spatial resolution of 6.3 mm for OSEM-TOF and 5.5 mm for OSEM-TOF+PSF. In the patient validation, the maximum percentage difference in lesion quantification between standard and optimized protocol (whole-body scan time of 15 vs. 5 min) was below 19%.

CONCLUSIONS

A reduction of single-bed and whole-body scan time for [⁶⁸Ga]Ga-PSMA PET/CT compared to current recommended clinical acquisition protocols is postulated. Clinical studies are warranted to validate the applicability of this protocol.

Phantom-based acquisition time and image reconstruction parameter optimisation for oncologic FDG PET/CT examinations using a digital system

Background

New-generation silicon-photomultiplier (SiPM)-based PET/CT systems exhibit an improved lesion detectability and image quality due to a higher detector sensitivity. Consequently, the acquisition time can be reduced while maintaining diagnostic quality. The aim of this study was to determine the lowest ¹⁸F-FDG PET acquisition time without loss of diagnostic information and to optimise image reconstruction parameters (image reconstruction algorithm, number of iterations, voxel size, Gaussian filter) by phantom imaging. Moreover, patient data are evaluated to confirm the phantom results.

Methods

Three phantoms were used: a soft-tissue tumour phantom, a bone-lung tumour phantom, and a resolution phantom. Phantom conditions (lesion sizes from 6.5 mm to 28.8 mm in diameter, lesion activity concentration of 15 kBq/mL, and signal-to-background ratio of 5:1) were derived from patient data. PET data were acquired on an SiPM-based Biograph Vision PET/CT system for 10 min in list-mode format and resampled into time frames from 30 to 300 s in 30-s increments to simulate different acquisition times. Different image reconstructions with varying iterations, voxel sizes, and Gaussian filters were probed. Contrast-to-noise-ratio (CNR), maximum, and peak signal were evaluated using the 10-min acquisition time image as reference. A threshold CNR value ≥ 5 and a maximum (peak) deviation of ± 20% were considered acceptable. 20 patient data sets were evaluated regarding lesion quantification as well as agreement and correlation between reduced and full acquisition time standard uptake values (assessed by Pearson correlation coefficient, intraclass correlation coefficient, Bland–Altman analyses, and Krippendorff’s alpha).

Results

An acquisition time of 60 s per bed position yielded acceptable detectability and quantification results for clinically relevant phantom lesions ≥ 9.7 mm in diameter using OSEM-TOF or OSEM-TOF+PSF image reconstruction, a 4-mm Gaussian filter, and a 1.65 × 1.65 x 2.00-mm³ or 3.30 × 3.30 x 3.00-mm³ voxel size. Correlation and agreement of patient lesion quantification between full and reduced acquisition times were excellent.

Conclusion

A threefold reduction in acquisition time is possible. Patients might benefit from more comfortable examinations or reduced radiation exposure, if instead of the acquisition time the applied activity is reduced.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09993-4.

Hybrid kernelised expectation maximisation for Bremsstrahlung SPECT reconstruction in SIRT with 90Y micro-spheres

Background

Selective internal radiation therapy with Yttrium-90 microspheres is an effective therapy for liver cancer and liver metastases. Yttrium-90 is mainly a high-energy beta particle emitter. These beta particles emit Bremsstrahlung radiation during their interaction with tissue making post-therapy imaging of the radioactivity distribution feasible. Nevertheless, image quality and quantification is difficult due to the continuous energy spectrum which makes resolution modelling, attenuation and scatter estimation challenging and therefore the dosimetry quantification is inaccurate. As a consequence a reconstruction algorithm able to improve resolution could be beneficial.

Methods

In this study, the hybrid kernelised expectation maximisation (HKEM) is used to improve resolution and contrast and reduce noise, in addition a modified HKEM called frozen HKEM (FHKEM) is investigated to further reduce noise. The iterative part of the FHKEM kernel was frozen at the 72nd sub-iteration. When using ordered subsets algorithms the data is divided in smaller subsets and the smallest algorithm iterative step is called sub-iteration. A NEMA phantom with spherical inserts was used for the optimisation and validation of the algorithm, and data from 5 patients treated with Selective internal radiation therapy were used as proof of clinical relevance of the method.

Results

The results suggest a maximum improvement of 56% for region of interest mean recovery coefficient at fixed coefficient of variation and better identification of the hot volumes in the NEMA phantom. Similar improvements were achieved with patient data, showing 47% mean value improvement over the gold standard used in hospitals.

Conclusions

Such quantitative improvements could facilitate improved dosimetry calculations with SPECT when treating patients with Selective internal radiation therapy, as well as provide a more visible position of the cancerous lesions in the liver.

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.

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.

Automatic quantification package (Hone Graph) for phantom-based image quality assessment in bone SPECT: computerized automatic classification of detectability

OBJECTIVE

We previously developed a custom-design thoracic bone scintigraphy-specific phantom ("SIM² bone phantom") to assess image quality in bone single-photon emission computed tomography (SPECT). We aimed to develop an automatic assessment system for imaging technology in bone SPECT and demonstrate the validity of this system.

METHODS

Four spherical lesions of 13-, 17-, 22-, and 28-mm diameters in the vertebrae of SIM² bone phantom simulating the thorax were filled with radioactivity (target-to-background ratio: 4). Dynamic SPECT acquisitions were performed for 15 min; reconstructions were performed using ordered subset expectation maximization at 3-15-min timepoints. Consequently, 216 lesions (54 SPECT images) were obtained: 120 and 96 lesions were used for software development and validation, respectively. The developed software used statistical parametric mapping to rigidly register and automatically calculate quantitative indexes (contrast-to-noise ratio, % coefficient of variance, % detectability equivalence volume, recovery coefficient, target-to-normal bone ratio, and full width at half maximum). A detectability score (DS) was used to define the four observation types (4, excellent; 3, adequate; 2, average; 1, poor) to score hot spherical lesions. The gold standard for DSs was independently classified by three experienced board-certified nuclear medicine technologists using the four observation types; thereafter, a consensus regarding the gold standard for DSs was reached. Using 120 lesions for development, decision tree analysis was performed to determine DS based on the quantitative indexes. We verified the validation of the quantitative indexes and their threshold values for automatic classification using 96 lesions for validation.

RESULTS

The trends in the automatically calculated quantitative indices were consistent. Decision tree analysis produced four terminal groups; two quantitative indexes (% detectability equivalence volume and contrast-to-noise ratio) were used to classify DS. The automatically classified DSs exhibited an almost perfect agreement with the gold standard. The percentage agreement and kappa coefficient were 91.7% and 0.93, respectively, in 96 lesions for validation.

CONCLUSIONS

The developed software automatically classified the detectability of hot lesions in the SIM² bone phantom using the automatically calculated quantitative indexes, suggesting that this software could provide a means to automatically perform detectability analysis after data input that is excellent in reproducibility and accuracy.

Impact of contouring methods on pre-treatment and post-treatment dosimetry for the prediction of tumor control and survival in HCC patients treated with selective internal radiation therapy

INTRODUCTION

The aim of this study was to evaluate the impact of the contouring methods on dose metrics and their predictive value on tumor control and survival, in both situations of pre-treatment and post-treatment dosimetry, for patients with advanced HCC treated with SIRT.

METHODS

Forty-eight patients who underwent SIRT between 2012 and 2020 were retrospectively included in this study. Target volumes were delineated using two methods: MRI-based contours manually drawn by a radiologist and then registered on SPECT/CT and PET/CT via deformable registration (Pre-C_MRI and Post-C_MRI), ^99mTc-MAA-SPECT and ⁹⁰Y-microspheres-PET 10% threshold contouring (Pre-C_SPECT and Post-C_PET). The mean absorbed dose (Dm) and the minimal absorbed dose delivered to 70% of the tumor volume (D70) were evaluated with both contouring methods; the tumor-to-normal liver uptake ratio (TNR) was evaluated with MRI-based contours only. Tumor response was assessed using the mRECIST criteria on the follow-up MRIs.

RESULTS

No significant differences were found for Dm and TNR between pre- and post-treatment. TNR evaluated with radiologic contours (Pre-C_MRI and Post-C_MRI) were predictive of tumor control at 6 months on pre- and post-treatment dosimetry (OR 5.9 and 7.1, respectively; p = 0.02 and 0.01). All dose metrics determined with both methods were predictive of overall survival (OS) on pre-treatment dosimetry, but only Dm with MRI-based contours was predictive of OS on post-treatment images with a median of 23 months for patients with a supramedian Dm versus 14 months for the others (p = 0.04).

CONCLUSION

In advanced HCC treated with SIRT, Dm and TNR determined with radiologic contours were predictive of tumor control and OS. This study shows that a rigorous clinical workflow (radiologic contours + registration on scintigraphic images) is feasible and should be prospectively considered for improving therapeutic strategy.

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

PET images acquired after liver ⁹⁰Y radioembolization therapies are typically very noisy, which significantly challenges both visualization and quantification of activity distributions. To improve their noise characteristics, regularized iterative reconstruction algorithms such as block sequential regularized expectation maximization (Q.Clear for GE Healthcare, USA) have been proposed. In this study, we aimed to investigate the effects which different reconstruction algorithms may have on patient images, with reconstruction parameters initially narrowed down using phantom studies. Moreover, we evaluated the impact of these reconstruction methods on voxel-based dose distribution in phantom and patient studies (lesions and healthy livers). The International Electrotechnical Commission (IEC)/NEMA phantom, containing six spheres, was filled with ⁹⁰Y and imaged using a GE Discovery 690 PET/CT scanner with time-of-flight enabled. The images were reconstructed using Q.Clear (with β parameter ranging from 0 to 8000) and ordered subsets expectation maximization. The image quality and quantification accuracy were evaluated by computing the hot ([Formula: see text]) and cold ([Formula: see text]) contrast recovery coefficients, background variability (BV) and activity bias. Next, dose distributions and dose volume histograms were generated using MIM® software's SurePlan LiverY90 toolbox. Subsequently, parameters optimized in these phantom studies were applied to five patient datasets. Dose parameters, such as D_max, D_mean, D₇₀, and V_100Gy, were estimated, and their variability for different reconstruction methods was investigated. Based on phantom studies, the β parameter values optimized for image quality and quantification accuracy were 2500 and 300, respectively. When all investigated reconstructions were applied to patient studies, D_mean, D₅₀, D₇₀, and V_100Gy showed coefficients of variation below 8%; whereas the variability of D_max was up to 30% for both phantom and patient images. Although β = 300-1000 would provide accurate activity quantification for a region of interest, when considering activity/dose voxelized distribution, higher β value (e.g. 4000-5000) would provide the greatest accuracy for dose distributions. In this ⁹⁰Y radioembolization PET/CT study, the β parameter in regularized iterative (Q.Clear) reconstruction was investigated for image quality, accurate quantification and dose distributions based on phantom experiments and then applied to patient studies. Our results indicate that more accurate dose distribution can be achieved from smoother PET images, reconstructed with larger β values than those yielding the best activity quantifications but noisy images. Most importantly, these results suggest that quantitative measures, which are commonly used in clinics, such as SUV_max or SUV_peak( equivalent of D_max), should not be employed for ⁹⁰Y PET images, since their values would highly depend on the image reconstruction.

Local experience in Saudi Arabia of selective internal radiation therapy with yttrium-90 microspheres for the treatment of unresectable liver tumours

Aim:

To assess local experience in Saudi Arabia of the effectiveness and tolerability of selective internal radiation therapy (SIRT) in patients with liver tumours.

Materials and methods:

In a retrospective study, patients with unresectable liver tumours treated with yttrium-90-labelled resin microspheres (SIR-Spheres®, Sirtex, Australia) at the King Fahad Specialist Hospital Dammam (KFSHD) were followed up for at least 12 months, or until death. Data were extracted from medical records. The primary outcome measure was overall survival (OS). Tumour response was recorded using the modified Response Evaluation Criteria in Solid Tumors (mRECIST): tumour control (TC) was defined as the proportion of patients with stable disease (SD), partial response (PR) and complete response (CR).

Results:

A total of 30 patients (21 males, 9 females) with hepatocellular carcinoma (76⋅7%), metastatic colorectal cancer (13⋅3%), cholangiocarcinoma (6⋅7%) or carcinoid tumour (3⋅3%) were evaluated. Mean OS was 16⋅98 months. Eighty-five percent of patients showed TC, 35% had SD, 15% a PR and 35% a CR. No severe complications were observed. Four deaths were considered unrelated to treatment.

Findings:

In this cohort, SIRT showed similar efficacy to that in other studies, with an acceptable tolerability profile. SIRT appears a feasible procedure for liver tumour treatment in the KFSHD.

Dose volume histogram-based optimization of image reconstruction parameters for quantitative 90 Y-PET imaging

PURPOSE

⁹⁰ Y-microsphere radioembolization or selective internal radiation therapy is increasingly being used as a treatment option for tumors that are not candidates for surgery and external beam radiation therapy. Recently, volumetric ⁹⁰ Y-dosimetry techniques have been implemented to explore tumor dose-response on the basis of 3D ⁹⁰ Y-activity distribution from PET imaging. Despite being a theranostic study, the optimization of quantitative ⁹⁰ Y-PET image reconstruction still uses the mean activity concentration recovery coefficient (RC) as the objective function, which is more relevant to diagnostic and detection tasks than is to dosimetry. The aim of this study was to optimize ⁹⁰ Y-PET image reconstruction by minimizing errors in volumetric dosimetry via the dose volume histogram (DVH). We propose a joint optimization of the number of equivalent iterations (the product of the iterations and subsets) and the postreconstruction filtration (FWHM) to improve the accuracy of voxel-level ⁹⁰ Y dosimetry.

METHODS

A modified NEMA IEC phantom was used to emulate clinically relevant ⁹⁰ Y-PET imaging conditions through various combinations of acquisition durations, activity concentrations, sphere-to-background ratios, and sphere diameters. PET data were acquired in list mode for 300 min in a single-bed position; we then rebinned the list mode PET data to 60, 45, 30, 15, and 5 min per bed, with 10 different realizations. Errors in the DVH were calculated as root mean square errors (RMSE) of the differences in the image-based DVH and the expected DVH. The new optimization approach was tested in a phantom study, and the results were compared with the more commonly used objective function of the mean activity concentration RC.

RESULTS

In a wide range of clinically relevant imaging conditions, using 36 equivalent iterations with a 5.2-mm filtration resulted in decreased systematic errors in volumetric ⁹⁰ Y dosimetry, quantified as image-based DVH, in ⁹⁰ Y-PET images reconstructed using the ordered subset expectation maximization (OSEM) iterative reconstruction algorithm with time of flight (TOF) and point spread function (PSF) modeling. Our proposed objective function of minimizing errors in DVH, which allows for joint optimization of ⁹⁰ Y-PET iterations and filtration for volumetric quantification of the ⁹⁰ Y dose, was shown to be superior to conventional RC-based optimization approaches for image-based absorbed dose quantification.

CONCLUSION

Our proposed objective function of minimizing errors in DVH, which allows for joint optimization of iterations and filtration to reduce errors in the PET-based volumetric quantification ⁹⁰ Y dose, is relevant to dosimetry in therapy procedures. The proposed optimization method using DVH as the objective function could be applied to any imaging modality used to assess voxel-level quantitative information.

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.

Phantom validation of quantitative Y-90 PET/CT-based dosimetry in liver radioembolization

Background

PET/CT has recently been shown to be a viable alternative to traditional post-infusion imaging methods providing good quality images of ⁹⁰Y-laden microspheres after selective internal radiation therapy (SIRT). In the present paper, first we assessed the quantitative accuracy of ⁹⁰Y-PET using an anthropomorphic phantom provided with lungs, liver, spine, and a cylindrical homemade lesion located into the hepatic compartment. Then, we explored the accuracy of different computational approaches on dose calculation, including (I) direct Monte Carlo radiation transport using Raydose, (II) Kernel convolution using Philips Stratos, (III) local deposition algorithm, (IV) Monte Carlo technique (MCNP) considering a uniform activity distribution, and (V) MIRD (Medical Internal Radiation Dose) analytical approach. Finally, calculated absorbed doses were compared with those obtained performing measurements with LiF:Mg,Cu,P TLD chips in a liquid environment.

Results

Our results indicate that despite ⁹⁰Y-PET being likely to provide high-resolution images, the ⁹⁰Y low branch ratio, along with other image-degrading factors, may produce non-uniform activity maps, even in the presence of uniform activity. A systematic underestimation of the recovered activity, both for the tumor insert and for the liver background, was found. This is particularly true if no partial volume correction is applied through recovery coefficients. All dose algorithms performed well, the worst case scenario providing an agreement between absorbed dose evaluations within 20%. Average absorbed doses determined with the local deposition method are in excellent agreement with those obtained using the MIRD and the kernel-convolution dose calculation approach.

Finally, absorbed dose assessed with MC codes are in good agreement with those obtained using TLD in liquid solution, thus confirming the soundness of both calculation approaches. This is especially true for Raydose, which provided an absorbed dose value within 3% of the measured dose, well within the stated uncertainties.

Conclusions

Patient-specific dosimetry is possible even in a scenario with low true coincidences and high random fraction, as in ⁹⁰Y–PET imaging, granted that accurate absolute PET calibration is performed and acquisition times are sufficiently long. Despite Monte Carlo calculations seeming to outperform all dose estimation algorithms, our data provide a strong argument for encouraging the use of the local deposition algorithm for routine ⁹⁰Y dosimetry based on PET/CT imaging, due to its simplicity of implementation.

Improved quantitative 90 Y bremsstrahlung SPECT/CT reconstruction with Monte Carlo scatter modeling

PURPOSE

In ⁹⁰ Y microsphere radioembolization (RE), accurate post-therapy imaging-based dosimetry is important for establishing absorbed dose versus outcome relationships for developing future treatment planning strategies. Additionally, accurately assessing microsphere distributions is important because of concerns for unexpected activity deposition outside the liver. Quantitative ⁹⁰ Y imaging by either SPECT or PET is challenging. In ⁹⁰ Y SPECT model based methods are necessary for scatter correction because energy window-based methods are not feasible with the continuous bremsstrahlung energy spectrum. The objective of this work was to implement and evaluate a scatter estimation method for accurate ⁹⁰ Y bremsstrahlung SPECT/CT imaging.

METHODS

Since a fully Monte Carlo (MC) approach to ⁹⁰ Y SPECT reconstruction is computationally very demanding, in the present study the scatter estimate generated by a MC simulator was combined with an analytical projector in the 3D OS-EM reconstruction model. A single window (105 to 195-keV) was used for both the acquisition and the projector modeling. A liver/lung torso phantom with intrahepatic lesions and low-uptake extrahepatic objects was imaged to evaluate SPECT/CT reconstruction without and with scatter correction. Clinical application was demonstrated by applying the reconstruction approach to five patients treated with RE to determine lesion and normal liver activity concentrations using a (liver) relative calibration.

RESULTS

There was convergence of the scatter estimate after just two updates, greatly reducing computational requirements. In the phantom study, compared with reconstruction without scatter correction, with MC scatter modeling there was substantial improvement in activity recovery in intrahepatic lesions (from > 55% to > 86%), normal liver (from 113% to 104%), and lungs (from 227% to 104%) with only a small degradation in noise (13% vs. 17%). Similarly, with scatter modeling contrast improved substantially both visually and in terms of a detectability index, which was especially relevant for the low uptake extrahepatic objects. The trends observed for the phantom were also seen in the patient studies where lesion activity concentrations and lesion-to-liver concentration ratios were lower for SPECT without scatter correction compared with reconstruction with just two MC scatter updates: in eleven lesions the mean uptake was 4.9 vs. 7.1 MBq/mL (P = 0.0547), the mean normal liver uptake was 1.6 vs. 1.5 MBq/mL (P = 0.056) and the mean lesion-to-liver uptake ratio was 2.7 vs. 4.3 (P = 0.0402) for reconstruction without and with scatter correction respectively.

CONCLUSIONS

Quantitative accuracy of ⁹⁰ Y bremsstrahlung imaging can be substantially improved with MC scatter modeling without significant degradation in image noise or intensive computational requirements.

Comparison of Step-and-Shoot and Continuous-Bed-Motion PET Modes of Acquisition for Limited-View Organ Scans

Continuous-bed-motion (CBM) acquisition mode has been made commercially available in PET/CT scanners. CBM mode is designed for whole-body imaging, with a long scan length (multiple axial fields of view [aFOVs]) and short acquisition duration (2-3 min/aFOV). PET/CT has recently been used after ⁹⁰Y-microsphere therapy to quantify ⁹⁰Y activity distribution in the liver. Here we compared counting efficiencies along the bed-motion direction (z-axis) between CBM and step-and-shoot (SS) acquisition modes for limited-view organ scans, such as ⁹⁰Y PET/CT liver studies, that have short scan lengths (≤2 aFOVs) and long acquisition durations (10-30 min/aFOV). Methods: The counting efficiencies, that is, analytic sensitivities, in SS mode (single-aFOV and multiple-aFOV scans) and CBM mode were theoretically derived and experimentally validated using a cylindric ⁶⁸Ge phantom. The sensitivities along the z-axis were compared between the SS and CBM modes. Results: The analytic and experimental count profiles were in good agreement, validating the analytic models. For fixed scan durations, the overall coincidence counting efficiency in CBM mode was lower (∼60%) than those in SS modes, and the maximum sensitivity in CBM mode was 50% or less of that in 1-aFOV SS mode and 100% or less of that in 2-aFOV SS mode. Conclusion: The ability of CBM mode to tailor-fit the PET/CT scan length and local scan duration is not realized in studies with a short scan length (≤30 cm) and long scan duration (20 min/aFOV for the scanner). SS acquisition mode is preferable to CBM mode for limited-view organ and count-starved scans, such as ⁹⁰Y PET/CT liver scans, because of the higher counting efficiency of SS mode, which leads to better image quality and quantification precision.

Optimisation of reconstruction, volumetry and dosimetry for 99mTc-SPECT and 90Y-PET images: Towards reliable dose-volume histograms for selective internal radiation therapy with 90Y-microspheres

PURPOSE

In Selective Internal Radiation Therapy (SIRT), ^99mTc-MAA SPECT images are commonly used to predict microspheres distribution but recent works used ⁹⁰Y-microspheres PET images. Nevertheless, evaluation of the predictive power of ^99mTc-MAA has been hampered by the lack of reliable comparisons between ^99mTc-SPECT and ⁹⁰Y-PET images. Our aim was to determine the "in situ" optimisation procedure in order to reliably compare ^99mTc-SPECT and ⁹⁰Y-PET images and achieve optimal personal dosimetry.

METHODS

We acquired ^99mTc-SPECT/CT and ⁹⁰Y-PET/CT images of NEMA and Jaszczak phantoms. We found the best reconstruction parameters for quantification and for volume estimations. We determined adaptive threshold curves on the volumetric reconstruction. We copied the optimised volumes on the quantitative reconstruction, named here the "cross volumes" technique. Finally, we compared ^99mTc-SPECT and ⁹⁰Y-PET Dose Volume Histograms.

RESULTS

Our "in situ" optimisation procedure decreased errors on volumes and quantification (from -44.2% and -15.8% to -3.4% and -3.28%, respectively, for the 26.5mL PET phantom sphere). Moreover, ^99mTc-SPECT and ⁹⁰Y-PET DVHs were equivalent only after the optimisation procedure (difference in mean dose <5% for the three biggest spheres).

CONCLUSIONS

This work showed that a preliminary "in situ" phantom study was necessary to optimise volumes and quantification of ^99mTc-SPECT and ⁹⁰Y-PET images and allowed to achieve a reliable comparison between patient treatment planning and post implant dosimetry, notably by the use of the "cross volumes" technique. Methodology developed in this work will enable robust evaluations of the predictive power of ^99mTc-SPECT, as well as dose-response relationship and side effects in SIRT treatments.

A phantom study: Should 124 I-mIBG PET/CT replace 123 I-mIBG SPECT/CT?

PURPOSE

The isotope ¹²³ I is commonly labeled with meta-iodobenzylguanidine (mIBG) for imaging of neuroendocrine tumors, such as pheochromocytomas and neuroblastomas. ¹²³ I-mIBG SPECT/CT imaging is performed for staging, follow-up and selection of patients for treatment with ¹³¹ I mIBG. As an alternative to ¹²³ I, ¹²⁴ I-mIBG PET/CT may be used, potentially taking advantage of the superior PET image quality. The purpose of this study was to investigate whether ¹²⁴ I PET/CT improves image quality as compared with ¹²³ I SPECT/CT for equal patient effective radiation dose (in mSv).

METHODS

Phantom measurements were performed using the NEMA-2007 image quality phantom. SPECT and PET reconstruction settings were used with and without time-of-flight (TOF) and point-spread-function (PSF) modeling. As a measure of image quality, the contrast-to-noise ratio (CNR) was calculated. The ratio of the ¹²³ I to ¹²⁴ I activity concentration was determined at which the contrast-to-noise ratio was equal for both modalities. This metric was defined as the contrast equivalent activity ratio (CEAR).

RESULTS

CEARs of 47.7, 25.6, 23.1, 14.6, 10.0, and 9.1 were obtained for a TOF and PSF modeled ¹²⁴ I reconstruction method and an attenuation and scatter-corrected ¹²³ I reconstruction method for sphere sizes of 10 to 37 mm, respectively. As the effective radiation dose of ¹²⁴ I-mIBG is higher than of ¹²³ I-mIBG (in mSv/MBq), an equal effective dose corresponds to a CEAR of 5 to 10. Therefore, CEARs higher than 5 to 10 indicate that ¹²⁴ I PET/CT outperforms ¹²³ I SPECT/CT in the sense of image quality for equal patient effective radiation dose.

CONCLUSION

The CEAR is much larger than a factor of 5 to 10 (needed for equal patient effective radiation dose) for most of the reconstruction methods and sphere sizes. Therefore, ¹²⁴ I-mIBG PET/CT is expected to improve image quality and/or may be used to reduce effective patient dose as compared with ¹²³ I-mIBG SPECT/CT.

Practical reconstruction protocol for quantitative (90)Y bremsstrahlung SPECT/CT

PURPOSE

To develop a practical background compensation (BC) technique to improve quantitative (90)Y-bremsstrahlung single-photon emission computed tomography (SPECT)/computed tomography (CT) using a commercially available imaging system.

METHODS

All images were acquired using medium-energy collimation in six energy windows (EWs), ranging from 70 to 410 keV. The EWs were determined based on the signal-to-background ratio in planar images of an acrylic phantom of different thicknesses (2-16 cm) positioned below a (90)Y source and set at different distances (15-35 cm) from a gamma camera. The authors adapted the widely used EW-based scatter-correction technique by modeling the BC as scaled images. The BC EW was determined empirically in SPECT/CT studies using an IEC phantom based on the sphere activity recovery and residual activity in the cold lung insert. The scaling factor was calculated from 20 clinical planar (90)Y images. Reconstruction parameters were optimized in the same SPECT images for improved image quantification and contrast. A count-to-activity calibration factor was calculated from 30 clinical (90)Y images.

RESULTS

The authors found that the most appropriate imaging EW range was 90-125 keV. BC was modeled as 0.53× images in the EW of 310-410 keV. The background-compensated clinical images had higher image contrast than uncompensated images. The maximum deviation of their SPECT calibration in clinical studies was lowest (<10%) for SPECT with attenuation correction (AC) and SPECT with AC + BC. Using the proposed SPECT-with-AC + BC reconstruction protocol, the authors found that the recovery coefficient of a 37-mm sphere (in a 10-mm volume of interest) increased from 39% to 90% and that the residual activity in the lung insert decreased from 44% to 14% over that of SPECT images with AC alone.

CONCLUSIONS

The proposed EW-based BC model was developed for (90)Y bremsstrahlung imaging. SPECT with AC + BC gave improved lesion detectability and activity quantification compared to SPECT with AC only. The proposed methodology can readily be used to tailor (90)Y SPECT/CT acquisition and reconstruction protocols with different SPECT/CT systems for quantification and improved image quality in clinical settings.

A gate evaluation of the sources of error in quantitative90 Y PET

PURPOSE

Accurate reconstruction of the dose delivered by ⁹⁰ Y microspheres using a postembolization PET scan would permit the establishment of more accurate dose-response relationships for treatment of hepatocellular carcinoma with ⁹⁰ Y. However, the quality of the PET data obtained is compromised by several factors, including poor count statistics and a very high random fraction. This work uses Monte Carlo simulations to investigate what impact factors other than low count statistics have on the quantification of⁹⁰ Y PET.

METHODS

PET acquisitions of two phantoms-a NEMA PET phantom and the NEMA IEC PET body phantom-containing either ⁹⁰ Y or ¹⁸ F were simulated using gate. Simulated projections were created with subsets of the simulation data allowing the contributions of random, scatter, and LSO background to be independently evaluated. The simulated projections were reconstructed using the commercial software for the simulated scanner, and the quantitative accuracy of the reconstruction and the contrast recovery of the reconstructed images were evaluated.

RESULTS

The quantitative accuracy of the ⁹⁰ Y reconstructions were not strongly influenced by the high random fraction present in the projection data, and the activity concentration was recovered to within 5% of the known value. The contrast recovery measured for simulated ⁹⁰ Y data was slightly poorer than that for simulated ¹⁸ F data with similar count statistics. However, the degradation was not strongly linked to any particular factor. Using a more restricted energy range to reduce the random fraction in the projections had no significant effect.

CONCLUSIONS

Simulations of ⁹⁰ Y PET confirm that quantitative ⁹⁰ Y is achievable with the same approach as that used for ¹⁸ F, and that there is likely very little margin for improvement by attempting to model aspects unique to ⁹⁰ Y, such as the much higher random fraction or the presence of bremsstrahlung in the singles data.

Evaluating the Application of Tissue-Specific Dose Kernels Instead of Water Dose Kernels in Internal Dosimetry: A Monte Carlo Study

PURPOSE

The aim of this work is to evaluate the application of tissue-specific dose kernels instead of water dose kernels to improve the accuracy of patient-specific dosimetry by taking tissue heterogeneities into consideration.

MATERIALS AND METHODS

Tissue-specific dose point kernels (DPKs) and dose voxel kernels (DVKs) for yttrium-90 (⁹⁰Y), lutetium-177 (¹⁷⁷Lu), and phosphorus-32 (³²P) are calculated using the Monte Carlo (MC) simulation code GATE (version 7). The calculated DPKs for bone, lung, adipose, breast, heart, intestine, kidney, liver, and spleen are compared with those of water. The dose distribution in normal and tumorous tissues in lung, liver, and bone of a Zubal phantom is calculated using tissue-specific DVKs instead of those of water in conventional methods. For a tumor defined in a heterogeneous region in the Zubal phantom, the absorbed dose is calculated using a proposed algorithm, taking tissue heterogeneity into account. The algorithm is validated against full MC simulations.

RESULTS

The simulation results indicate that the highest differences between water and other tissue DPKs occur in bone for ⁹⁰Y (12.2% ± 0.6%), ³²P (18.8% ± 1.3%), and ¹⁷⁷Lu (16.9% ± 1.3%). The second highest discrepancy corresponds to the lung for ⁹⁰Y (6.3% ± 0.2%), ³²P (8.9% ± 0.4%), and ¹⁷⁷Lu (7.7% ± 0.3%). For ⁹⁰Y, the mean absorbed dose in tumorous and normal tissues is calculated using tissue-specific DVKs in lung, liver, and bone. The results are compared with doses calculated considering the Zubal phantom water equivalent and the relative differences are 4.50%, 0.73%, and 12.23%, respectively. For the tumor in the heterogeneous region of the Zubal phantom that includes lung, liver, and bone, the relative difference between mean calculated dose in tumorous and normal tissues based on the proposed algorithm and the values obtained from full MC dosimetry is 5.18%.

CONCLUSIONS

A novel technique is proposed considering tissue-specific dose kernels in the dose calculation algorithm. This algorithm potentially enables patient-specific dosimetry and improves estimation of the average absorbed dose of ⁹⁰Y in a tumor located in lung, bone, and soft tissue interface by 6.98% compared with the conventional methods.

Surefire infusion system versus standard microcatheter use during holmium-166 radioembolization: study protocol for a randomized controlled trial

Background

An anti-reflux catheter (ARC) may increase the tumor absorbed dose during radioembolization (RE) by elimination of particle reflux and its effects on hemodynamics. Since the catheter is fixed in a centro-luminal position, it may also increase the predictive accuracy of a scout dose administration before treatment. The purpose of the SIM trial is to compare the effects of ARC use during RE with holmium-166 (¹⁶⁶Ho) microspheres in patients with colorectal liver metastases (CRLM), with the use of a standard end-hole microcatheter.

Methods/Design

A within-patient randomized controlled trial (RCT) will be conducted in 25 patients with unresectable chemorefractory liver-dominant CRLM. Study participants will undergo a ¹⁶⁶Ho scout dose procedure in the morning and a therapeutic procedure in the afternoon. The ARC will be randomly allocated to the left/right hepatic artery, and a standard microcatheter will be used in the contralateral artery. SPECT/CT imaging will be performed for quantitative analyses of the microsphere distribution directly after the scout and treatment procedure. Baseline and follow-up investigations include ¹⁸F-FDG-PET + liver CT, clinical and laboratory examinations. The primary endpoint is the comparison of tumor to non-tumor (T/N) activity ratio in both groups. Secondary endpoints include comparisons of mean absorbed dose in tumors and healthy liver tissue, infusion efficiency, the predictive value of ¹⁶⁶Ho scout dose for tumor response. In the entire cohort, a dose-response relationship, clinical toxicity, and overall survival will be assessed. The sample was determined for the expectation that the ARC will increase the T/N ratio by 25 % (mean T/N ratio 2.0 vs. 1.6).

Discussion

The SIM trial is a within-patient RCT that will assess whether ¹⁶⁶Ho RE treatment can be optimized by using an ARC.

Trial registration

The SIM trial is registered at clinicaltrials.gov (NCT02208804). Registered on 31 July 2014.

Electronic supplementary material

The online version of this article (doi:10.1186/s13063-016-1643-3) contains supplementary material, which is available to authorized users.

Optimization of Image Reconstruction for 90Y Selective Internal Radiotherapy on a Lutetium Yttrium Orthosilicate PET/CT System Using a Bayesian Penalized Likelihood Reconstruction Algorithm

Imaging on a γ-camera with ⁹⁰Y after selective internal radiotherapy (SIRT) may allow for verification of treatment delivery but suffers relatively poor spatial resolution and imprecise dosimetry calculation. ⁹⁰Y PET/CT imaging is possible on 3-dimensional, time-of-flight machines; however, images are usually poor because of low count statistics and noise. A new PET reconstruction software using a Bayesian penalized likelihood (BPL) reconstruction algorithm (termed Q.Clear) was investigated using phantom and patient scans to optimize the reconstruction for post-SIRT imaging and clarify whether BPL leads to an improvement in clinical image quality using ⁹⁰Y. Methods: Phantom studies over an activity range of 0.5-4.2 GBq were performed to assess the contrast recovery, background variability, and contrast-to-noise ratio for a range of BPL and ordered-subset expectation maximization (OSEM) reconstructions on a PET/CT scanner. Patient images after SIRT were reconstructed using the same parameters and were scored and ranked on the basis of image quality, as assessed by visual evaluation, with the corresponding SPECT/CT Bremsstrahlung images by 2 experienced radiologists. Results: Contrast-to-noise ratio was significantly better in BPL reconstructions when compared with OSEM in phantom studies. The patient-derived BPL and matching Bremsstrahlung images scored higher than OSEM reconstructions when scored by radiologists. BPL with a β value of 4,000 was ranked the highest of all images. Deadtime was apparent in the system above a total phantom activity of 3.3 GBq. Conclusion: BPL with a β value of 4,000 is the optimal image reconstruction in PET/CT for confident radiologic reading when compared with other reconstruction parameters for ⁹⁰Y imaging after SIRT imaging. Activity in the field of view should be below 3.3 GBq at the time of PET imaging to avoid deadtime losses for this scanner.

Positron Emission Tomography: Current Challenges and Opportunities for Technological Advances in Clinical and Preclinical Imaging Systems

Positron emission tomography (PET) imaging is based on detecting two time-coincident high-energy photons from the emission of a positron-emitting radioisotope. The physics of the emission, and the detection of the coincident photons, give PET imaging unique capabilities for both very high sensitivity and accurate estimation of the in vivo concentration of the radiotracer. PET imaging has been widely adopted as an important clinical modality for oncological, cardiovascular, and neurological applications. PET imaging has also become an important tool in preclinical studies, particularly for investigating murine models of disease and other small-animal models. However, there are several challenges to using PET imaging systems. These include the fundamental trade-offs between resolution and noise, the quantitative accuracy of the measurements, and integration with X-ray computed tomography and magnetic resonance imaging. In this article, we review how researchers and industry are addressing these challenges.

Analysis of the influence of 111In on 90Y-bremsstrahlung SPECT based on Monte Carlo simulation

OBJECTIVE

⁹⁰Y-ibritumomab tiuxetan (Zevalin) which is used for the treatment of malignant lymphomas can be used for SPECT imaging based on bremsstrahlung from ⁹⁰Y beta particles. However, gamma rays emitted by ¹¹¹In, which is administered to evaluate the indication for the treatment, contaminate the ⁹⁰Y bremsstrahlung images. Our objective is to investigate the influence of ¹¹¹In on the ⁹⁰Y SPECT images using Monte Carlo simulation.

METHODS

We used an in-house developed simulation code for the Monte Carlo simulation of electrons and photons (MCEP). Two hot spheres with diameters of 40 mm were put in an elliptical phantom. Both spheres ("sphere 1" and "sphere 2") were filled with ⁹⁰Y and ¹¹¹In mixed solutions. The activities of ⁹⁰Y in sphere 1 and sphere 2 were 241 and 394 kBq/mL, respectively, and the ones of ¹¹¹In were 8.14 and 13.3 kBq/mL, respectively. The background activity of ⁹⁰Y was 38.6 kBq/mL, whereas that of ¹¹¹In was 1.30 kBq/mL; moreover, the acquisition time was 30 min. Two energy windows were used: one is 90-190 keV included the ¹¹¹In photopeak; the other is 90-160 keV. To evaluate the quality of the SPECT images, the contrast recovery coefficient (CRC) and the constant noise ratio (CNR) of the SPECT images were derived.

RESULTS

For the energy window between 90 and 160 keV, the ¹¹¹In count was 74 % of the total. In that case, the CRC values were 30.1 and 30.7 % for "sphere 1" and "sphere 2", respectively, whereas the CNR values were 6.8 and 12.1, respectively. For the energy window between 90 and 190 keV, the ¹¹¹In count reached 85 % of the total count. The CRC and CNR values were 38.6 and 40.0 % and 10.6 and 19.4, respectively.

CONCLUSIONS

Our simulation study revealed that the cross talk between ¹¹¹In and ⁹⁰Y in SPECT imaging is rather serious. Even for the energy window excluding the ¹¹¹In photopeak, the count ratio of ⁹⁰Y was less than 30 % of the total. However, the influence of ¹¹¹In on ⁹⁰Y-SPECT imaging cannot be ignored, and the count ratio because of ¹¹¹In is important to estimate the density of ⁹⁰Y.

Physics of pure and non-pure positron emitters for PET: a review and a discussion

With the increased interest in new PET tracers, gene-targeted therapy, immunoPET, and theranostics, other radioisotopes will be increasingly used in clinical PET scanners, in addition to ¹⁸F. Some of the most interesting radioisotopes with prospective use in the new fields are not pure short-range β⁺ emitters but can be associated with gamma emissions in coincidence with the annihilation radiation (prompt gamma), gamma-gamma cascades, intense Bremsstrahlung radiation, high-energy positrons that may escape out of the patient skin, and high-energy gamma rays that result in some e⁺/e⁻ pair production. The high level of sophistication in data correction and excellent quantitative accuracy that has been reached for ¹⁸F in recent years can be questioned by these effects. In this work, we review the physics and the scientific literature and evaluate the effect of these additional phenomena on the PET data for each of a series of radioisotopes: ¹¹C, ¹³N, ¹⁵O, ¹⁸F, ⁶⁴Cu, ⁶⁸Ga, ⁷⁶Br, ⁸²Rb, ⁸⁶Y, ⁸⁹Zr, ⁹⁰Y, and ¹²⁴I. In particular, we discuss the present complications arising from the prompt gammas, and we review the scientific literature on prompt gamma correction. For some of the radioisotopes considered in this work, prompt gamma correction is definitely needed to assure acceptable image quality, and several approaches have been proposed in recent years. Bremsstrahlung photons and ¹⁷⁶Lu background were also evaluated.

Characteristics of Bremsstrahlung emissions of (177)Lu, (188)Re, and (90)Y for SPECT/CT quantification in radionuclide therapy

PURPOSE

Beta particles emitted by radioisotopes used in targeted radionuclide therapies (TRT) create Bremsstrahlung (BRS) which may affect SPECT quantification when imaging these isotopes. The purpose of the current study was to investigate the characteristics of Bremsstrahlung produced in tissue by three β-emitting radioisotopes used in TRT.

METHODS

Monte Carlo simulations of (177)Lu, (188)Re, and (90)Y sources placed in water filled cylinders were performed. BRS yields, mean energies and energy spectra for (a) all photons generated in the decays, (b) photons that were not absorbed and leave the cylinder, and (c) photons detected by the camera were analyzed. Next, the results of simulations were compared with those from experiments performed on a clinical SPECT camera using same acquisition conditions and phantom configurations as in simulations.

RESULTS

Simulations reproduced relatively well the shapes of the measured spectra, except for (90)Y which showed an overestimation in the low energy range. Detailed analysis of the results allowed us to suggest best collimators and imaging conditions for each of the investigated isotopes. Finally, our simulations confirmed that the BRS contribution to the energy spectra in quantitative imaging of (177)Lu and (188)Re could be ignored.

CONCLUSIONS

For (177)Lu and (188)Re, BRS contributes only marginally to the total spectra recorded by the camera. Our analysis shows that MELP and HE collimators are the best for imaging these two isotopes. For (90)Y, HE collimator should be used.

PET/MRI of Hepatic 90Y Microsphere Deposition Determines Individual Tumor Response

Purpose

The purpose of our study is to determine if there is a relationship between dose deposition measured by PET/MRI and individual lesion response to yttrium-90 (⁹⁰Y) microsphere radioembolization.

Materials and Methods

26 patients undergoing lobar treatment with ⁹⁰Y microspheres underwent PET/MRI within 66 h of treatment and had follow-up imaging available. Adequate visualization of tumor was available in 24 patients, and contours were drawn on simultaneously acquired PET/MRI data. Dose volume histograms (DVHs) were extracted from dose maps, which were generated using a voxelized dose kernel. Similar contours to capture dimensional and volumetric change of tumors were drawn on follow-up imaging. Response was analyzed using both RECIST and volumetric RECIST (vRECIST) criteria.

Results

A total of 8 hepatocellular carcinoma (HCC), 4 neuroendocrine tumor (NET), 9 colorectal metastases (CRC) patients, and 3 patients with other metastatic disease met inclusion criteria. Average dose was useful in predicting response between responders and non-responders for all lesion types and for CRC lesions alone using both response criteria (p < 0.05). D70 (minimum dose to 70 % of volume) was also useful in predicting response when using vRECIST. No significant trend was seen in the other tumor types. For CRC lesions, an average dose of 29.8 Gy offered 76.9 % sensitivity and 75.9 % specificity for response.

Conclusions

PET/MRI of ⁹⁰Y microsphere distribution showed significantly higher DVH values for responders than non-responders in patients with CRC. DVH analysis of ⁹⁰Y microsphere distribution following treatment may be an important predictor of response and could be used to guide future adaptive therapy trials.

Monte Carlo simulation of PET and SPECT imaging of 90Y

PURPOSE

Yittrium-90 ((90)Y) is traditionally thought of as a pure beta emitter, and is used in targeted radionuclide therapy, with imaging performed using bremsstrahlung single-photon emission computed tomography (SPECT). However, because (90)Y also emits positrons through internal pair production with a very small branching ratio, positron emission tomography (PET) imaging is also available. Because of the insufficient image quality of (90)Y bremsstrahlung SPECT, PET imaging has been suggested as an alternative. In this paper, the authors present the Monte Carlo-based simulation-reconstruction framework for (90)Y to comprehensively analyze the PET and SPECT imaging techniques and to quantitatively consider the disadvantages associated with them.

METHODS

Our PET and SPECT simulation modules were developed using Monte Carlo simulation of Electrons and Photons (MCEP), developed by Dr. S. Uehara. PET code (MCEP-PET) generates a sinogram, and reconstructs the tomography image using a time-of-flight ordered subset expectation maximization (TOF-OSEM) algorithm with attenuation compensation. To evaluate MCEP-PET, simulated results of (18)F PET imaging were compared with the experimental results. The results confirmed that MCEP-PET can simulate the experimental results very well. The SPECT code (MCEP-SPECT) models the collimator and NaI detector system, and generates the projection images and projection data. To save the computational time, the authors adopt the prerecorded (90)Y bremsstrahlung photon data calculated by MCEP. The projection data are also reconstructed using the OSEM algorithm. The authors simulated PET and SPECT images of a water phantom containing six hot spheres filled with different concentrations of (90)Y without background activity. The amount of activity was 163 MBq, with an acquisition time of 40 min.

RESULTS

The simulated (90)Y-PET image accurately simulated the experimental results. PET image is visually superior to SPECT image because of the low background noise. The simulation reveals that the detected photon number in SPECT is comparable to that of PET, but the large fraction (approximately 75%) of scattered and penetration photons contaminates SPECT image. The lower limit of (90)Y detection in SPECT image was approximately 200 kBq/ml, while that in PET image was approximately 100 kBq/ml.

CONCLUSIONS

By comparing the background noise level and the image concentration profile of both the techniques, PET image quality was determined to be superior to that of bremsstrahlung SPECT. The developed simulation codes will be very useful in the future investigations of PET and bremsstrahlung SPECT imaging of (90)Y.

Quantitative Comparison of 124I PET/CT and 131I SPECT/CT Detectability

Radioiodine therapy with (131)I is used for treatment of suspected recurrence of differentiated thyroid carcinoma. Pretherapeutic (124)I PET/CT with a low activity (~1% of (131)I activity) can be performed to determine whether uptake of (131)I, and thereby the desired therapeutic effect, may be expected. However, false-negative (124)I PET/CT results as compared with posttherapeutic (131)I SPECT/CT have been reported by several groups. The purpose of this study was to investigate whether the reported discrepancies may be ascribed to a difference in lesion detectability between (124)I PET/CT and (131)I SPECT/CT and, hence, whether the administered (124)I activity is sufficient to achieve equal detectability.

METHODS

Phantom measurements were performed using the National Electrical Manufacturers Association 2007 image-quality phantom. As a measure of detectability, the contrast-to-noise ratio was calculated. The (124)I activity was expressed as the percentage of (131)I activity required to achieve the same contrast-to-noise ratio. This metric was defined as the detectability equivalence percentage (DEP).

RESULTS

Because lower DEPs were obtained for smaller spheres, a relatively low (124)I activity was sufficient to achieve similar lesion detectability between (124)I PET/CT and (131)I SPECT/CT. DEP was 1.5%, 1.9%, 1.9%, 4.4%, 9.0%, and 16.2% for spheres with diameters of 10, 13, 17, 18, 25, and 37 mm, respectively, for attenuation- and scatter-corrected SPECT versus point-spread function (PSF) model-based and time-of-flight (TOF) PET. For no-PSF no-TOF PET, DEP was 3.6%, 2.1%, 3.5%, 7.8%, 15.1%, and 23.3%, respectively.

CONCLUSION

A relatively low (124)I activity of 74 MBq (~1% of (131)I activity) is sufficient to achieve similar lesion detectability between (124)I PSF TOF PET/CT and (131)I SPECT/CT for small spheres (≤10 mm), since the reported DEPs are close to 1%. False-negative (124)I PET/CT results as compared with posttherapeutic (131)I SPECT/CT may be ascribed to differences in detectability for large lesions (>10 mm) and for no-PSF no-TOF PET, since DEPs are greater than 1%. On the basis of DEPs of 3.5% for lesion diameters of up to 17 mm on no-PSF no-TOF PET, (124)I activities as high as 170 MBq may be warranted to obtain equal detectability.

Beyond 18F-FDG: Characterization of PET/CT and PET/MR Scanners for a Comprehensive Set of Positron Emitters of Growing Application--18F, 11C, 89Zr, 124I, 68Ga, and 90Y

This study aimed to investigate image quality for a comprehensive set of isotopes ((18)F, (11)C, (89)Zr, (124)I, (68)Ga, and (90)Y) on 2 clinical scanners: a PET/CT scanner and a PET/MR scanner.

METHODS

Image quality and spatial resolution were tested according to NU 2-2007 of the National Electrical Manufacturers Association. An image-quality phantom was used to measure contrast recovery, residual bias in a cold area, and background variability. Reconstruction methods available on the 2 scanners were compared, including point-spread-function correction for both scanners and time of flight for the PET/CT scanner. Spatial resolution was measured using point sources and filtered backprojection reconstruction.

RESULTS

With the exception of (90)Y, small differences were seen in the hot-sphere contrast recovery of the different isotopes. Cold-sphere contrast recovery was similar across isotopes for all reconstructions, with an improvement seen with time of flight on the PET/CT scanner. The lower-statistic (90)Y scans yielded substantially lower contrast recovery than the other isotopes. When isotopes were compared, there was no difference in measured spatial resolution except for PET/MR axial spatial resolution, which was significantly higher for (124)I and (68)Ga.

CONCLUSION

Overall, both scanners produced good images with (18)F, (11)C, (89)Zr, (124)I, (68)Ga, and (90)Y.

Theranostic Imaging of Yttrium-90

This paper overviews Yttrium-90 ((90)Y) as a theranostic and nuclear medicine imaging of (90)Y radioactivity with bremsstrahlung imaging and positron emission tomography. In addition, detection and optical imaging of (90)Y radioactivity using Cerenkov luminescence will also be reviewed. Methods and approaches for qualitative and quantitative (90)Y imaging will be briefly discussed. Although challenges remain for (90)Y imaging, continued clinical demand for predictive imaging response assessment and target/nontarget dosimetry will drive research and technical innovation to provide greater clinical utility of (90)Y as a theranostic agent.

A multicentre comparison of quantitative (90)Y PET/CT for dosimetric purposes after radioembolization with resin microspheres : The QUEST Phantom Study

Purpose

To investigate and compare the quantitative accuracy of ⁹⁰Y imaging across different generation PET/CT scanners, for the purpose of dosimetry after radioembolization with resin microspheres.

Methods

A strict experimental and imaging protocol was followed by 47 international sites using the NEMA 2007/IEC 2008 PET body phantom with an 8-to-1 sphere-to-background ratio of ⁹⁰Y solution. The phantom was imaged over a 7-day period (activity ranging from 0.5 to 3.0 GBq) and all reconstructed data were analysed at a core laboratory for consistent processing. Quantitative accuracy was assessed through measures of total phantom activity, activity concentration in background and hot spheres, misplaced counts in a nonradioactive insert, and background variability.

Results

Of the 69 scanners assessed, 37 had both time-of-flight (ToF) and resolution recovery (RR) capability. These current generation scanners from GE, Philips and Siemens could reconstruct background concentration measures to within 10 % of true values over the evaluated range, with greater deviations on the Philips systems at low count rates, and demonstrated typical partial volume effects on hot sphere recovery, which dominated spheres of diameter <20 mm. For spheres >20 mm in diameter, activity concentrations were consistently underestimated by about 20 %. Non-ToF scanners from GE Healthcare and Siemens were capable of producing accurate measures, but with inferior quantitative recovery compared with ToF systems.

Conclusion

Current generation ToF scanners can consistently reconstruct ⁹⁰Y activity concentrations, but they underestimate activity concentrations in small structures (≤37 mm diameter) within a warm background due to partial volume effects and constraints of the reconstruction algorithm. At the highest count rates investigated, measures of background concentration (about 300 kBq/ml) could be estimated on average to within 1 %, 5 % and 2 % for GE Healthcare (all-pass filter, RR + ToF), Philips (4i8s ToF) and Siemens (2i21s all-pass filter, RR + ToF) ToF systems, respectively. Over the range of activities investigated, comparable performance between GE Healthcare and Siemens ToF systems suggests suitability for quantitative analysis in a scenario analogous to that of postradioembolization imaging for treatment of liver cancer.

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.

(90)Y-PET/CT Imaging Quantification for Dosimetry in Peptide Receptor Radionuclide Therapy: Analysis and Corrections of the Impairing Factors

PURPOSE

We evaluated the possibility to assess (90)Y-PET/CT imaging quantification for dosimetry in (90)Y-peptide receptor radionuclide therapy.

METHODS

Tests were performed by Discovery 710 Elite (GE) PET/CT equipment. A body-phantom containing radioactive-coplanar-spheres was filled with (90)Y water solution to reproduce different signal-to-background-activity-ratios (S/N). We studied minimum detectable activity (MDA) concentration, contrast-to-noise ratio (CNR), and full-width-at-half-maximum (FWHM). Subsequently, three recovery coefficients (RC)-based correction approaches were evaluated: maximum-RC, resolution-RC, and isovolume-RC. The analysis of the volume segmentation thresholding method was also assessed to derive a relationship between the true volume of the targets and the threshold to be applied to the PET images. (90)Y-PET/CT imaging quantification was then achieved on some patients and related with preclinical tests. Moreover, the dosimetric evaluation was obtained on the target regions.

RESULTS

CNR value was greater than 5 if the MDA was greater than 0.2 MBq/mL with no background activity and 0.5-0.7 MBq/mL with S/N ranging from 3 to 6. FWHM was equal to 7 mm. An exponential fitting of isovolume RCs-based correction technique was adopted for activity quantification. Adaptive segmentation thresholding exponential curves were obtained and applied for target volume identification in three signal-to-background-activity-ratios. The imaging quantification study and dosimetric evaluations in clinical cases was feasible and the results were coherent with those obtained in preclinical tests.

CONCLUSIONS

(90)Y-PET/CT imaging quantification is possible both in phantoms and in patients. Absorbed dose evaluations in clinical applications are strongly related to targets activity concentration.

Dosimetry of a (90)Y-hydroxide liquid brachytherapy treatment approach to canine osteosarcoma using PET/CT

A new treatment strategy based on direct injections of (90)Y-hydroxide into the tumor bed in dogs with osteosarcoma was studied. Direct injections of the radiopharmaceutical into the tumor bed were made according to a pretreatment plan established using (18)F-FDG images. Using a special drill, cannulas were inserted going through tissue, tumor and bone. Using these cannulas, direct injections of the radiopharmaceutical were made. The in vivo biodistribution of (90)Y-hydroxide and the anatomical tumor bed were imaged using a time-of-flight (TOF) PET/CT scanner. The material properties of the tissues were estimated from corresponding CT numbers using an electron-density calibration. Radiation absorbed dose estimates were calculated using Monte Carlo methods where the biodistribution of the pharmaceutical from PET images was sampled using a collapsing 3-D rejection technique. Dose distributions in the tumor bed and surrounding tissues were calculated, showing significant heterogeneity with multiple hot spots at injection sites. Dose volume histograms showed that approximately 33.9% of bone and tumor and 70.2% of bone marrow and trabecular bone received an absorbed dose over 200Gy; approximately 3.2% of bone and tumor and 31.0% of bone marrow and trabecular bone received a total dose of over 1000Gy.

Quantitative and qualitative assessment of Yttrium-90 PET/CT imaging

Yttrium-90 is known to have a low positron emission decay of 32 ppm that may allow for personalized dosimetry of liver cancer therapy with ⁹⁰Y labeled microspheres. The aim of this work was to image and quantify ⁹⁰Y so that accurate predictions of the absorbed dose can be made. The measurements were performed within the QUEST study (University of Sydney, and Sirtex Medical, Australia). A NEMA IEC body phantom containing 6 fillable spheres (10–37 mm ∅) was used to measure the ⁹⁰Y distribution with a Biograph mCT PET/CT (Siemens, Erlangen, Germany) with time-of-flight (TOF) acquisition. A sphere to background ratio of 8∶1, with a total ⁹⁰Y activity of 3 GBq was used. Measurements were performed for one week (0, 3, 5 and 7 d). he acquisition protocol consisted of 30 min-2 bed positions and 120 min-single bed position. mages were reconstructed with 3D ordered subset expectation maximization (OSEM) and point spread function (PSF) for iteration numbers of 1–12 with 21 (TOF) and 24 (non-TOF) subsets and CT based attenuation and scatter correction. Convergence of algorithms and activity recovery was assessed based on regions-of-interest (ROI) analysis of the background (100 voxels), spheres (4 voxels) and the central low density insert (25 voxels). For the largest sphere, the recovery coefficient (RC) values for the 30 min –2-bed position, 30 min-single bed and 120 min-single bed were 1.12±0.20, 1.14±0.13, 0.97±0.07 respectively. For the smaller diameter spheres, the PSF algorithm with TOF and single bed acquisition provided a comparatively better activity recovery. Quantification of Y-90 using Biograph mCT PET/CT is possible with a reasonable accuracy, the limitations being the size of the lesion and the activity concentration present. At this stage, based on our study, it seems advantageous to use different protocols depending on the size of the lesion.