TARE is an effective treatment option for intermediate- or advanced-stage HCC or in cases of portal vein thrombosis. In addition to the low complication rate, the reported outcomes of TARE are comparable to or even better than those of other treatment options such as TACE or systemic chemotherapy [17–20]. In TARE, radiopharmaceuticals are directly injected into the tumor-feeding arteries and do not show redistribution after initial embolization. Thus, calculation of its treatment dose is easier than that of other systemically administered radioactive drugs. Pretreatment planning angiography is currently performed using 99mTc-MAA to simulate radiation doses to the tumor, normal liver tissues, and lungs.
In addition to pretreatment scans, many institutions perform a post-treatment scan to assess the results of TARE and identify possible complications occurring because of exposure of unexpected organs to radioactivity. Because 90Y is a pure beta-emitter, the post-treatment scans are performed using highly sensitive scanners to acquire bremsstrahlung gamma images or PET images. The radiation dosimetry methods using post-treatment images are based on the MIRD schema. Several studies have applied the MIRD schema to bremsstrahlung imaging [21, 22], and Gulec et al. reported the use of this schema for the dosimetry of 90Y-microspheres confined to the liver [23]. The usual MIRD schema assumes a uniform distribution of radioactivity in a certain organ or tissues and its use is very simple.
However, in the real world, therapeutic radiopharmaceuticals show heterogeneous distributions in-target tissues or organs because of vasculature heterogeneity, anatomical variation, and tissue necrosis. This study adopted voxel-based dosimetry, which considered intra-tumoral heterogeneity of 90Y-microsphere distribution [24]. In this voxel-based dosimetry study, the VSV was calculated for each voxel to create an absorbed dose map. This is an easier approach than that involving patient-specific Monte Carlo simulations and showed accurate dosimetry results in uniform-density organs such as the liver [24]. We applied a VSV from a 90Y dataset published previously [13].
Based on the voxel-based dose map, the TDv of each lesion was successfully calculated. The considerable variation in TDv, from 40 Gy to 177 Gy, may be the main reason behind the difference in treatment outcomes. In survival analysis, small tumor size, low AST level, and low TDv (median, ≥81 Gy) were significant prognostic factors for successful treatment. The multivariate analysis confirmed TDv as the only independent prognostic factor. Previous studies reported clinical factors related to hepatic function (total bilirubin, albumin) and tumor aggressiveness (alpha-fetoprotein level, portal vein thrombosis, and tumor size) as significant prognostic factors in patients treated with TARE [25, 26]. In our study, these factors were not significant, probably because of the small sample size. Further investigation in a larger cohort is needed to determine the role of these factors.
When applying a partition model of dosimetry to pretreatment simulation scans, 120 Gy is deemed a target dose for effective treatment. In our study, some lesions were outside this target range. Kao et al. reported that HCC tumors receiving > 105 Gy showed a higher response rate when the tumor dose was calculated using a partition model-based dosimetry model on 99mTc-MAA SPECT/CT [27]. Kokabi et al. estimated the tumor dose using 90Y bremsstrahlung SPECT/CT and reporting a significant correlation between a tumor dose of > 105 Gy and prolonged survival [28]. We observed an obvious difference in PFS between groups showing high and low absorbed doses. In Fig. 4, the cutoff of the TDv was set at 80–120 Gy. We also observed significant correlations between the absorbed dose and PFS in correlation analysis; this finding is consistent with that of a previous study [29].
Because this study performed voxel-based dosimetry, it was possible to calculate the voxel-wise intra-tumoral heterogeneity. High intra-tumoral heterogeneity was assumed to be related to LCF, despite the high average tumor dose. The SD and CV of dose distributions in a tumor were measured as simple indices for heterogeneity [30]. Although there was no significant difference in these values between the LCS and LCF groups in our study, further studies including more cases are warranted to investigate the effects of intra-tumoral dose distributions.
This study has several limitations. First, a small number of patients with inoperable HCC who underwent a single TARE session with 90Y-resin microspheres were included in this study. Further studies with large cohorts are required to investigate the role of voxel-based dosimetry in clinical practice. Second, we used only simple indices for intra-tumoral heterogeneity. However, the results of this study may be used as preliminary data for determining the role of intra-tumoral heterogeneity of absorbed doses in treatment response.