Optimization of Y-90 Radioembolization Imaging for Post-Treatment Dosimetry on a Long Axial Field-of-View PET/CT Scanner

BACKGROUND

PET imaging after yttrium-90 (Y-90) radioembolization is challenging because of the low positron fraction of Y-90 (32 × 10-6). The resulting low number of events can be compensated by the high sensitivity of long axial field-of-view (LAFOV) PET/CT scanners. Nevertheless, the reduced event statistics require optimization of the imaging protocol to achieve high image quality (IQ) and quantification accuracy sufficient for post-treatment dosimetry.

METHODS

Two phantoms (NEMA IEC and AbdoMan phantoms, mimicking human liver) filled with Y-90 and a 4:1 sphere (tumor)-to-background ratio were scanned for 24 h with the Biograph Vision Quadra (Siemens Healthineers). Eight patients were scanned after Y-90 radioembolization (1.3-4.7 GBq) using the optimized protocol (obtained by phantom studies). The IQ, contrast recovery coefficients (CRCs) and noise were evaluated for their limited and full acceptance angles, different rebinned scan durations, numbers of iterations and post-reconstruction filters. The s-value-based absorbed doses were calculated to assess their suitability for dosimetry.

RESULTS

The phantom studies demonstrate that two iterations, five subsets and a 4 mm Gaussian filter provide a reasonable compromise between a high CRC and low noise. For a 20 min scan duration, an adequate CRC of 56% (vs. 24 h: 62%, 20 mm sphere) was obtained, and the noise was reduced by a factor of 1.4, from 40% to 29%, using the full acceptance angle. The patient scan results were consistent with those from the phantom studies, and the impacts on the absorbed doses were negligible for all of the studied parameter sets, as the maximum percentage difference was -3.89%.

CONCLUSIONS

With 2i5s, a 4 mm filter and a scan duration of 20 min, IQ and quantification accuracy that are suitable for post-treatment dosimetry of Y-90 radioembolization can be achieved.

The American Brachytherapy Society consensus statement for permanent implant brachytherapy using Yttrium-90 microsphere radioembolization for liver tumors

PURPOSE

To develop a multidisciplinary consensus for high quality multidisciplinary implementation of brachytherapy using Yttrium-90 (90Y) microspheres transarterial radioembolization (90Y TARE) for primary and metastatic cancers in the liver.

METHODS AND MATERIALS

Members of the American Brachytherapy Society (ABS) and colleagues with multidisciplinary expertise in liver tumor therapy formulated guidelines for 90Y TARE for unresectable primary liver malignancies and unresectable metastatic cancer to the liver. The consensus is provided on the most recent literature and clinical experience.

RESULTS

The ABS strongly recommends the use of 90Y microsphere brachytherapy for the definitive/palliative treatment of unresectable liver cancer when recommended by the multidisciplinary team. A quality management program must be implemented at the start of 90Y TARE program development and follow-up data should be tracked for efficacy and toxicity. Patient-specific dosimetry optimized for treatment intent is recommended when conducting 90Y TARE. Implementation in patients on systemic therapy should account for factors that may enhance treatment related toxicity without delaying treatment inappropriately. Further management and salvage therapy options including retreatment with 90Y TARE should be carefully considered.

CONCLUSIONS

ABS consensus for implementing a safe 90Y TARE program for liver cancer in the multidisciplinary setting is presented. It builds on previous guidelines to include recommendations for appropriate implementation based on current literature and practices in experienced centers. Practitioners and cooperative groups are encouraged to use this document as a guide to formulate their clinical practices and to adopt the most recent dose reporting policies that are critical for a unified outcome analysis of future effectiveness studies.

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.

Characterization of point-spread function specification error on Geometric Transfer Matrix partial volume correction in [11C]PiB amyloid imaging

Purpose

Partial-volume correction (PVC) using the Geometric Transfer Matrix (GTM) method is used in positron emission tomography (PET) to compensate for the effects of spatial resolution on quantitation. We evaluate the effect of misspecification of scanner point-spread function (PSF) on GTM results in amyloid imaging, including the effect on amyloid status classification (positive or negative).

Methods

Twenty-nine subjects with Pittsburgh Compound B ([¹¹C]PiB) PET and structural T1 MR imaging were analyzed. FreeSurfer 5.3 (FS) was used to parcellate MR images into regions-of-interest (ROIs) that were used to extract radioactivity concentration values from the PET images. GTM PVC was performed using our “standard” PSF parameterization [3D Gaussian, full-width at half-maximum (w) of approximately 5 mm]. Additional GTM PVC was performed with “incorrect” parameterizations, taken around the correct value. The result is a set of regional activity values for each of the GTM applications. For each case, activity values from various ROIs were combined and normalized to produce standardized uptake value ratios (SUVRs) for nine standard [¹¹C]PiB quantitation ROIs and a global region. GTM operating-point characteristics were determined from the slope of apparent SUVR versus w curves.

Results

Errors in specification of w on the order of 1 mm (3D) mainly produce only modest errors of up to a few percent. An exception was the anterior ventral striatum in which fractional errors of up to 0.29 per millimeter (3D) of error in w were observed.

Conclusion

While this study does not address all the issues regarding the quantitative strengths and weakness of GTM PVC, we find that with reasonable caution, the unavoidable inaccuracies associated with PSF specification do not preclude its use in amyloid quantitation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00403-5.

3D voxel-based dosimetry to predict contralateral hypertrophy and an adequate future liver remnant after lobar radioembolization

INTRODUCTION

Volume changes induced by selective internal radiation therapy (SIRT) may increase the possibility of tumor resection in patients with insufficient future liver remnant (FLR). The aim was to identify dosimetric and clinical parameters associated with contralateral hepatic hypertrophy after lobar/extended lobar SIRT with ⁹⁰Y-resin microspheres.

MATERIALS AND METHODS

Patients underwent ⁹⁰Y PET/CT after lobar or extended lobar (right + segment IV) SIRT. ⁹⁰Y voxel dosimetry was retrospectively performed (PLANET Dose; DOSIsoft SA). Mean absorbed doses to tumoral/non-tumoral-treated volumes (NTL) and dose-volume histograms were extracted. Clinical variables were collected. Patients were stratified by FLR at baseline (T0-FLR): < 30% (would require hypertrophy) and ≥ 30%. Changes in volume of the treated, non-treated liver, and FLR were calculated at < 2 (T1), 2-5 (T2), and 6-12 months (T3) post-SIRT. Univariable and multivariable regression analyses were performed to identify predictors of atrophy, hypertrophy, and increase in FLR. The best cut-off value to predict an increase of FLR to ≥ 40% was defined using ROC analysis.

RESULTS

Fifty-six patients were studied; most had primary liver tumors (71.4%), 40.4% had cirrhosis, and 39.3% had been previously treated with chemotherapy. FLR in patients with T0-FLR < 30% increased progressively (T0: 25.2%; T1: 32.7%; T2: 38.1%; T3: 44.7%). No dosimetric parameter predicted atrophy. Both NTL-Dmean and NTL-V30 (fraction of NTL exposed to ≥ 30 Gy) were predictive of increase in FLR in patients with T0 FLR < 30%, the latter also in the total cohort of patients. Hypertrophy was not significantly associated with tumor dose or tumor size. When ≥ 49% of NTL received ≥ 30 Gy, FLR increased to ≥ 40% (accuracy: 76.4% in all patients and 80.95% in T0-FLR < 30% patients).

CONCLUSION

NTL-Dmean and NTL exposed to ≥ 30 Gy (NTL-V30) were most significantly associated with increase in FLR (particularly among patients with T0-FLR < 30%). When half of NTL received ≥ 30 Gy, FLR increased to ≥ 40%, with higher accuracy among patients with T0-FLR < 30%.

The joint use of 99mTc-MAA-SPECT/CT and cone-beam CT optimizes radioembolization planning

Purpose

To determine which imaging method used during radioembolization (RE) work-up: contrast-enhanced computed tomography (CECT), ^99mTc-MAA-SPECT/CT or cone beam-CT (CBCT), more accurately predicts the final target volume (TgV) as well as the influence that each modality has in the dosimetric calculation.

Methods

TgVs from ^99mTc-MAA-SPECT/CT, CECT and CBCT were consecutively obtained in 24 patients treated with RE and compared with ⁹⁰Y PET/CT TgV. Using the TgVs estimated by each imaging modality and a fictitious activity of 1 GBq, the corresponding absorbed doses by tumor and non-tumoral parenchyma were calculated for each patient. The absorbed doses for each modality were compared with the ones obtained using ⁹⁰Y PET/CT TgV.

Results

^99mTc-MAA-SPECT/CT predicted ⁹⁰Y PET/CT TgV better than CBCT or CECT, even for selective or superselective administrations. Likewise, ^99mTc-MAA-SPECT/CT showed dosimetric values more similar to those obtained with ⁹⁰Y PET/CT. Nevertheless, CBCT provided essential information for RE planning, such as ensuring the total coverage of the tumor and, in cases with more than one feeding artery, splitting the activity according to the volume of tumor perfused by each artery.

Conclusion

The joint use of ^99mTc-MAA-SPECT/CT and CBCT optimizes dosimetric planning for RE procedures, enabling a more accurate personalized approach.

A proof-of-concept study of the in-vivo validation of a computational fluid dynamics model of personalized radioembolization

Radioembolization (RE) with yttrium-90 (90Y) microspheres, a transcatheter intraarterial therapy for patients with liver cancer, can be modeled computationally. The purpose of this work was to correlate the results obtained with this methodology using in vivo data, so that this computational tool could be used for the optimization of the RE procedure. The hepatic artery three-dimensional (3D) hemodynamics and microsphere distribution during RE were modeled for six 90Y-loaded microsphere infusions in three patients with hepatocellular carcinoma using a commercially available computational fluid dynamics (CFD) software package. The model was built based on in vivo data acquired during the pretreatment stage. The results of the simulations were compared with the in vivo distribution assessed by 90Y PET/CT. Specifically, the microsphere distribution predicted was compared with the actual 90Y activity per liver segment with a commercially available 3D-voxel dosimetry software (PLANET Dose, DOSIsoft). The average difference between the CFD-based and the PET/CT-based activity distribution was 2.36 percentage points for Patient 1, 3.51 percentage points for Patient 2 and 2.02 percentage points for Patient 3. These results suggest that CFD simulations may help to predict 90Y-microsphere distribution after RE and could be used to optimize the RE procedure on a patient-specific basis.

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.

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

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

Phantom, clinical, and texture indices evaluation and optimization of a penalized-likelihood image reconstruction method (Q.Clear) on a BGO PET/CT scanner

INTRODUCTION

The aim of this study was to evaluate the behavior of a penalized-likelihood image reconstruction method (Q.Clear) under different count statistics and lesion-to-background ratios (LBR) on a BGO scanner, in order to obtain an optimum penalization factor (β value) to study and optimize for different acquisition protocols and clinical goals.

METHODS

Both phantom and patient images were evaluated. Data from an image quality phantom were acquired using different Lesion-to-Background ratios and acquisition times. Then, each series of the phantom was reconstructed using β values between 50 and 500, at intervals of 50. Hot and cold contrasts were obtained, as well as background variability and contrast-to-noise ratio (CNR). Fifteen ¹⁸ F-FDG patients (five brain scans and 10 torso acquisitions) were acquired and reconstructed using the same β values as in the phantom reconstructions. From each lesion in the torso acquisition, noise, contrast, and signal-to-noise ratio (SNR) were computed. Image quality was assessed by two different nuclear medicine physicians. Additionally, the behaviors of 12 different textural indices were studied over 20 different lesions.

RESULTS

Q.Clear quantification and optimization in patient studies depends on the activity concentration as well as on the lesion size. In the studied range, an increase on β is translated in a decrease in lesion contrast and noise. The net product is an overall increase in the SNR, presenting a tendency to a steady value similar to the CNR in phantom data. As the activity concentration or the sphere size increase the optimal β increases, similar results are obtained from clinical data. From the subjective quality assessment, the optimal β value for torso scans is in a range between 300 and 400, and from 100 to 200 for brain scans. For the recommended torso β values, texture indices present coefficients of variation below 10%.

CONCLUSIONS

Our phantom and patients demonstrate that improvement of CNR and SNR of Q.Clear algorithm which depends on the studied conditions and the penalization factor. Using the Q.Clear reconstruction algorithm in a BGO scanner, a β value of 350 and 200 appears to be the optimal value for 18F-FDG oncology and brain PET/CT, respectively.

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.

On the assessment of spatial resolution of PET systems with iterative image reconstruction

Spatial resolution is an important metric for performance characterization in PET systems. Measuring spatial resolution is straightforward with a linear reconstruction algorithm, such as filtered backprojection, and can be performed by reconstructing a point source scan and calculating the full-width-at-half-maximum (FWHM) along the principal directions. With the widespread adoption of iterative reconstruction methods, it is desirable to quantify the spatial resolution using an iterative reconstruction algorithm. However, the task can be difficult because the reconstruction algorithms are nonlinear and the non-negativity constraint can artificially enhance the apparent spatial resolution if a point source image is reconstructed without any background. Thus, it was recommended that a background should be added to the point source data before reconstruction for resolution measurement. However, there has been no detailed study on the effect of the point source contrast on the measured spatial resolution. Here we use point source scans from a preclinical PET scanner to investigate the relationship between measured spatial resolution and the point source contrast. We also evaluate whether the reconstruction of an isolated point source is predictive of the ability of the system to resolve two adjacent point sources. Our results indicate that when the point source contrast is below a certain threshold, the measured FWHM remains stable. Once the contrast is above the threshold, the measured FWHM monotonically decreases with increasing point source contrast. In addition, the measured FWHM also monotonically decreases with iteration number for maximum likelihood estimate. Therefore, when measuring system resolution with an iterative reconstruction algorithm, we recommend using a low-contrast point source and a fixed number of iterations.

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.

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.