Abstract
In hepatic 90Y radioembolization, pretreatment 99mTc-macroaggregated albumin (99mTc-MAA) nuclear imaging is used for lung shunt analysis, evaluation of extrahepatic deposition, and sometimes for treatment planning, using a partition model. A high level of agreement between pretreatment 99mTc-MAA distribution and final 90Y-microsphere distribution is assumed. The aim of this study was to investigate the value of pretreatment 99mTc-MAA SPECT to predict intrahepatic posttreatment 90Y-microsphere distribution. Methods: Volumes of interest (VOIs) were delineated on pretreatment contrast-enhanced CT or MR images according to Couinaud liver segmentation. All VOIs were registered to the 99mTc-MAA SPECT and 90Y SPECT images. The 99mTc-MAA SPECT and 90Y SPECT activity counts were normalized to the total administered activity of 90Y. For each VOI, this practice resulted in a predictive amount of 90Y (MBq/cm3) based on 99mTc-MAA SPECT in comparison with an actual amount of 90Y based on 90Y SPECT. Bland–Altman analysis was used to investigate the agreement of the activity distribution between 99mTc-MAA SPECT and 90Y SPECT. Results: A total of 39 procedures (225 VOIs) in 31 patients were included for analysis. The overall mean difference between pretreatment and posttreatment distribution of activity concentration for all segments was −0.022 MBq/cm3 with 95% limits of agreement of −0.581 to 0.537 MBq/cm3 (−28.9 to 26.7 Gy absorbed dose). A difference of >10%, >20%, and >30% of the mean activity per milliliter was found in, respectively, 153 (68%), 97 (43%), and 72 (32%) of the 225 segments. In every 99mTc-MAA procedure, at least 1 segment showed an under- or overestimation of >10%. The position of the catheter tip during administrations, as well as the tumor load of the liver segments, significantly influenced the disagreement. Conclusion: In current clinical practice, 99mTc-MAA distribution does not accurately predict final 90Y activity distribution. Awareness of the importance of catheter positioning and adherence to specific recommendations may lead to optimization of individualized treatment planning based on pretreatment imaging.
Radioembolization with 90Y microspheres is widely used for treatment of primary or metastatic liver malignancies. Selective injection of these microspheres in the hepatic artery results in high absorbed tumor doses while largely sparing the surrounding normal liver parenchyma, which is dependent mainly on the portal vein for its blood supply (1–3).
In a pretreatment angiographic procedure the anatomy of the liver vasculature is evaluated, and hepaticoenteric anastomoses that may lead to extrahepatic deposition of activity are occluded by coil embolization (4). Thereafter 99mTc-macroaggregated albumin particles (99mTc-MAA) are injected in the liver artery supplying the target volume, and the distribution of 99mTc-MAA is visualized by scintigraphy. Most centers that perform radioembolization use the distribution of 99mTc-MAA to calculate the lung shunt fraction and to detect any extrahepatic deposition of activity (5,6). Furthermore, it is assumed that 99mTc-MAA can also be used to predict the intrahepatic distribution of 90Y-microspheres, and as such, 99mTc-MAA is sometimes used for individualized treatment planning by the so-called partition model (7–9).
In most of the patients, the prescribed activities for radioembolization are calculated with methods based on liver weight or on a combination of body surface area and tumor liver involvement (10). The absorbed dose to the tumor and any accompanying toxicity effects to the normal liver parenchyma may be observed only after the actual treatment. Dosimetry can be used for individualized treatment planning and aims to optimize treatment efficacy with acceptable toxicity. Pretreatment dosimetry, however, requires a scout or safety dose as a reference for the treatment, for example by using the 99mTc-MAA distribution as a reference for posttreatment dose distribution.
The partition model is suggested as an alternative means of activity calculation in patients with a limited number of hypervascular liver tumors, optimizing the administered activity in individual patients (5). Tumor-to-nontumor activity ratios on pretreatment 99mTc-MAA SPECT are used to calculate activities that better reflect the intrahepatic dose distribution. In clinical practice, this means that the partition-model–based activity may be much higher than prescribed activities that are based on the more conventional methods, especially in patients with hypervascular tumors, having high tumor-to-nontumor activity ratios (11).
The partition model relies on 99mTc-MAA as a predictor for 90Y-microsphere distribution. However, the predictive value of 99mTc-MAA for the distribution of 90Y-microspheres in the liver is still a matter of debate (12). Parameters that may influence distribution differences between 99mTc-MAA and 90Y-microspheres include interval differences in catheter position, physiologic variances in hepatic blood flow, size and morphology differences between 99mTc-MAA particles and 90Y-microspheres, tumor histopathology, and tumor load. These and other factors may all limit the agreement between 99mTc-MAA and 90Y-microsphere distribution.
The aim of this study was to investigate the value of 99mTc-MAA to predict 90Y-microsphere distribution. Insight on this matter is essential for further development of any dose calculation method based on pretreatment 99mTc-MAA distribution.
MATERIALS AND METHODS
Patients
All patients who were treated with radioembolization from the start of our program in February 2009 up to February 2012 were retrospectively analyzed. The institutional review board approved this study and waived the requirements for patient informed consent. Patients who had received both the pretreatment administration of 99mTc-MAA and the treatment with 90Y-microspheres were included. Exclusion criteria were missing data, malregistration of imaging data, and procedures with multiple administrations, rendering it impossible to relate individual administrations to specific target volumes. In some patients, 2 separate lobar procedures were included in the study analysis. Subsegmental administrations were not performed.
Radioembolization
All procedures were performed according to international consensus (13). In short, during a pretreatment angiographic procedure, a 5-French catheter was used to evaluate hepatic vascular anatomy and to identify nontarget vessels leading to organs other than the liver. In general, the gastroduodenal artery and the right gastric artery were coil-embolized to prevent extrahepatic deposition of activity. Any other vessels branching off near the injection site and leading to nontarget organs were embolized as well. The cystic artery was not prophylactically embolized. Consecutively, a scout dose of 99mTc-MAA (150 MBq, 0.8 mg in 3.0 mL, TechneScan LyoMaa; Mallinckrodt Medical B.V.) was injected via a slow-pulsed injection, followed by planar imaging and SPECT, to check for inadvertent extrahepatic deposition. 99mTc-MAA was prepared immediately before administration; imaging was performed immediately after administration. If the 99mTc-MAA was not distributed to any nontarget area (including a hepatopulmonary shunt ≤ 20%), the patient was scheduled for treatment (mean interval, 12 d; range, 0–23 d). However, in the case of inadvertent 99mTc-MAA distribution, a second pretreatment procedure was performed to embolize the culprit vessels. In those cases, only the most recent 99mTc-MAA SPECT study was used for analyses. No vessels were coil-embolized after final 99mTc-MAA administration. Resin microspheres (SIR-Spheres; Sirtex) were used for treatment. Activity calculations were based on the body surface area method (prescribed activity in GBq, body surface area – 0.2 + fractional tumor involvement). Patients with a lung shunt fraction of >20% were excluded from treatment, and in patients with a lung shunt fraction of 10%–15% or 15%–20% a reduction of 20% and 40% was applied, respectively. Per protocol, the interventional radiologist placed the catheter tip in the same position during both procedures. The microspheres were infused slowly, with intermittent contrast injection and digital subtraction angiography to check for stasis. All administered activity was corrected for any residual activity after treatment.
Imaging
Pre- and posttreatment imaging was performed on a dual-head γ-camera (Forte [Philips] for 7 procedures and Symbia [Siemens Health Care] for 32 procedures). Pretreatment 99mTc-MAA planar and SPECT images were acquired on a 128 × 128 matrix using a 129.5- to 150.5-keV energy window and a low-energy general-purpose collimator. For posttreatment 90Y bremsstrahlung SPECT imaging, the combination of a high-energy general-purpose collimator and a wide 50- to 250-keV energy window was used, which yields images with a favorable combination of sensitivity and contrast (14). SPECT imaging was performed with 120 projections over a noncircular orbit of 180° (Forte; 30 s/projection) or 360° (Symbia; 20 s/projection). Data were reconstructed using ordered-subsets expectation maximization (5 iterations, 8 subsets) including attenuation correction and a gaussian postreconstruction filter of 5 mm in full width at half maximum. The reconstructed voxel size was 4.7 × 4.7 × 4.7 mm, 3.9 × 3.9 × 3.9 mm, and 4.8 × 4.8 × 4.8 mm, for Forte images, Symbia 99mTc-MAA images, and Symbia 90Y images, respectively.
Analysis
Contrast-enhanced CT or MR pretreatment images were used for liver segmentation according to the Bismuth adaptation of the Couinaud classification of liver anatomy (15). Software was developed at our institution for this purpose (Research Volumetool, version 1.3.3) (16). A maximum of 8 segments (volumes of interest, or VOIs) were delineated per patient. Prior liver resection and lobar procedures resulted in fewer segments. All delineated VOIs were manually registered to the 99mTc-MAA SPECT and 90Y SPECT images (Fig. 1). Procedures for which coregistration was impossible or inaccurate because of differences in liver position between the different scans were excluded from analysis.
The pretreatment 99mTc-MAA SPECT and the posttreatment 90Y SPECT images were converted into units of 90Y activity concentration by normalization of the total number of reconstructed counts in the VOIs to the total administered activity of 90Y. For each VOI, this practice resulted in a predictive amount of 90Y (MBq/cm3) based on 99mTc-MAA SPECT in comparison with an actual amount of 90Y based on 90Y SPECT. A homogeneous distribution of the activity inside a VOI, no activity distribution outside the liver, and no interval change in liver morphology were assumed. To illustrate the clinical implications of the disagreement, a map of the absorbed dose in grays was also calculated, using a conversion factor of 49.7 Gy/(MBq/cm3) (17).
The injection positions of 99mTc-MAA and 90Y-microspheres were retrospectively analyzed. Three observers independently reviewed the agreement between the 2 injection positions (per procedure) on fluoroscopy images on a 4-point scale (1, very poor agreement, difference > 10 mm; 2, poor agreement, difference > 5–10 mm; 3, good agreement, difference > 3–5 mm; 4, very good agreement, difference < 3 mm). A subgroup of patients with suboptimal agreement between the injection positions (average score ≤ 2.5) was selected. In this subgroup, the average difference in the catheter positions as scored by 3 reviewers was 5 mm or more. The injection positions were also classified as close to a major bifurcation (<10 mm) or not close to a major bifurcation, and segments were classified as having >25% tumor involvement or ≤25%.
Statistical Analysis
A commercial statistical software package (SPSS for Windows, version 20.0; SPSS Inc.) was used for data analysis. Bland–Altman plots were used for evaluating agreement between pre- and posttreatment activity distributions (18,19). A Bland–Altman graph is the preferred method to test for agreement between 2 instruments that are intended to measure the same parameter, in our case 99mTc-MAA SPECT and 90Y SPECT to measure final 90Y-microsphere distribution. In a Bland–Altman plot, the difference between the 2 methods is plotted against the mean of the 2 methods. The error was estimated by the mean difference (dm) and the SD of the differences (s). The 95% limits of agreement were calculated by dm ± 2s. Because the 99mTc-MAA was normalized to the 90Y activity, the expected mean difference is zero. In Bland–Altman analysis, the width of the distribution (i.e., 95% limits of agreement) is a measurement of the agreement between the 2 methods. An absorbed dose map in grays was calculated for translation to clinical practice. Cutoff levels (10%, 20%, and 30%) for the difference from the mean were used to evaluate variability in activity distribution. The Fisher exact test was used to test differences between subgroups with differences in tumor involvement, tumor cell type, and catheter positions.
RESULTS
Eighty patients were evaluated for radioembolization (Fig. 2). In 18 patients, contraindications to therapy were found during or after the 99mTc-MAA procedure. The remaining 62 patients underwent 73 treatment procedures. Multiple injections of activity in the same procedure led to exclusion of 27 procedures, and missing data led to exclusion of 4 procedures. In another 3 procedures, technical difficulties were encountered during coregistration because of substantial differences in liver position between the different scans (n = 2) or because of segmentation problems (n = 1). In the latter patient, extensive disease in the liver made it impossible to delineate the individual liver segments. A total of 39 procedures in 31 patients were included for analysis (Table 1). The mean administered activity of 90Y-microspheres was 1,002 MBq per procedure (range, 207–1,912 MBq). In 8 patients, 2 separate procedures were included, one for the left liver lobe and one for the right lobe. A total of 225 liver segments were analyzed.
The overall mean difference between pretreatment and posttreatment distribution of activity concentration for all segments was −0.022 MBq/cm3, with an SD of the mean of 0.285 MBq/cm3. A Bland–Altman plot was constructed with the absolute differences against their mean (Fig. 3). The 95% limits of agreement of the differences were −0.581 and 0.537 MBq/cm3, which correspond to 95% limits of agreement of −28.9 and 26.7 Gy absorbed dose.
A difference of >10%, >20%, and >30% of the mean activity per milliliter was found in, respectively, 153 (68%), 97 (43%), and 72 (32%) of the 225 segments (Fig. 4). In every 99mTc-MAA procedure, at least 1 segment showed an under- or overestimation of >10%. A >20% and >30% difference in at least 1 segment was found for 35 (90%) and 32 (82%) of 39 procedures, respectively (Fig. 5). A substantial difference in agreement between 99mTc-MAA and 90Y activity distribution was found for every procedure (Fig. 6).
Interestingly, the distribution differences were found to be smaller for segments with greater tumor involvement. The mean difference and 95% limits of agreement were −0.027 ± 0.603 and 0.007 ± 0.381 MBq/cm3 for segments with ≤25% and >25% tumor involvement, respectively. In grays, the 95% limits of agreement (−18.6 to 19.3 Gy) for segments with >25% tumor involvement had a smaller width than the 95% limits of agreement for segments with ≤25% tumor involvement (−31.3 to 28.6 Gy). A relative difference of >10%, >20%, and >30% was found in, respectively, 18 (51%), 9 (26%), and 7 (20%) of the 35 segments with >25% tumor involvement and in 136 (72%), 88 (47%), and 64 (34%) of the 190 segments with ≤25% tumor involvement. This proved significant for the cutoff values of >10% (P = 0.028) and >20% (P = 0.026), but significance was not reached for a cutoff value of >30% (P = 0.118).
A suboptimal agreement on catheter tip position was found for 9 patients (11 procedures). The mean difference and 95% limits of agreement were −0.008 ± 0.622 MBq/cm3 for procedures with a suboptimal agreement on catheter tip position (Fig. 7) and −0.026 ± 0.556 MBq/cm3 for procedures with an optimal agreement on catheter tip position. A relative difference of >10%, >20%, and >30% was found in, respectively, 54 (79%), 33 (49%), and 24 (35%) of the 68 segments for procedures with a suboptimal agreement on catheter tip position and in 100 (64%), 64 (41%), and 47 (30%) of the 157 segments for procedures with an optimal agreement on catheter tip position (Fig.7). This reached significance for the cutoff value of >10% (P = 0.020) but not for >20% (P = 0.307) and >30% (P = 0.438).
Significant differences were not found between procedures with the catheter tip near a major bifurcation (<10 mm) and procedures without, as long as the position was the same during the 99mTc-MAA and the 90Y injection. In procedures with the catheter tip close to a major bifurcation, as well as suboptimal agreement in catheter position between the 99mTc-MAA and the 90Y injection, a significant difference in activity distribution was found. A relative difference of >10%, >20%, and >30% was found in, respectively, 19 (95%), 12 (60%), and 9 (45%) of the 20 segments for procedures with these 2 characteristics, and in 133 (65%), 85 (41%), and 62 (31%) of the 205 segments for procedures without. This reached significance for the cutoff value of >10% (P = 0.005) but not for >20% (P = 0.155) and >30% (P = 0.209).
Evaluation of other parameters, including other combinations of parameters, did not yield any significant results. Procedures in patients with colorectal metastases (relatively hypovascular tumors) did not show any difference from procedures in patients with more hypervascular tumors, nor did treatment approach with regard to left lobar versus right lobar treatments.
DISCUSSION
It is expected that patients may greatly benefit from individualized treatment planning (10). Most promising in this regard is the so-called partition model. This method was previously shown to accurately predict treatment response and survival (11,20). The expected absorbed dose to the tumor is calculated on 99mTc-MAA SPECT. The tumors are delineated on morphologic images, whereas the dose distribution is estimated by calculation of 99mTc-MAA in the tumors, the normal liver, and the lungs. The prescribed activity may be calculated such that it does not exceed the maximum safe absorbed dose to the normal liver and lungs. Agreement between 99mTc-MAA and subsequent 90Y is therefore crucial for the accuracy of the partition method.
The presented data, however, show a substantial disagreement between 99mTc-MAA and 90Y activity distribution. In 68% of all segments, a difference of >10% between 99mTc-MAA and 90Y activity distribution was found. In every procedure, at least 1 segment showed a >10% difference. These findings raise concern about the validity of the partition method. However, regardless of the found disagreement, early studies were able to show the accuracy of the partition method nevertheless, although dose–effect relationships with regard to toxicity on the normal liver parenchyma were ignored (11,20). Knowledge of the existence, the magnitude, and the etiology of disagreement between 99mTc-MAA and 90Y activity distribution should ultimately lead to the improved validity of these methods. Specific technical and methodologic recommendations may help to overcome this issue.
Knesaurek et al. visually assessed the correlation between 99mTc-MAA and 90Y-distribution and found that correlation could vary from poor to relatively good (voxel-based Spearman rank correlation varied from 0.451 to 0.818) (21). However, their methodology is questionable. Their correlations indicate that 90Y activity is higher when 99mTc-MAA is higher. This does not imply that for individual measurements, the distributions are equal to (or close approximations of) each other. Agreement as described by Bland and Altman (18) is a more appropriate method of comparing 2 measurements of the same variable (i.e., activity distribution).
There are several factors that may have caused the disagreement in 99mTc-MAA and 90Y distribution. First, difference in catheter position between the 2 procedures seems to be a key factor. Although the catheter tip was positioned at the same location, small deviations (approximately 5–10 mm) were still found. A significantly increased disagreement between 99mTc-MAA and 90Y distribution was found in these procedures. In another study with more substantial differences in catheter tip position between the 99mTc-MAA and subsequent 90Y procedure, investigators found that position differences (P < 0.001) and a catheter position close to important side branches or bifurcations (P < 0.01) led to significant visually assessed distribution differences between 99mTc-MAA and 90Y (12). In the current study, we confirmed this finding quantitatively. In particular, patients with a mismatch in catheter tip position and injections close to a bifurcation showed significant disagreement. Selective administration of microspheres distal to the proper hepatic artery may largely overcome this issue. It was already shown that selective administrations are beneficial to prevent extrahepatic deposition of microsphere activity (22). In addition, differences in the intraluminal cross-sectional position of the catheter tip may have influenced the disagreement (Fig. 7). These differences are known to cause substantial differences in preferential flow (23–25). Innovative catheter designs with spacers to keep the catheter tip centered in the arterial lumen during injection may lead to more comparable injection positions and subsequent improved agreement between 99mTc-MAA and 90Y activity distribution (26).
Second, differences in the number, density, size, and morphology of the radiopharmaceuticals may also have resulted in a different activity distribution. The number of 99mTc-MAA particles (1–2 × 105 particles) is significantly lower than the number of 90Y microspheres applied (resin: 40–80 × 106; glass: 1.2–8 × 106 particles), whereas the density of 99mTc-MAA particles (1.1 g/mL) is lower than that of 90Y microspheres (resin: 1.6 g/mL; glass: 3.3 g/mL) (27–29). The particle size distribution of 99mTc-MAA is such that over 90% are within 10–90 μm in size (mean, 15 μm). The mean size of 90Y microspheres is 32 ± 10 μm, and the morphology of the spheric 90Y microspheres is also considerably different from the macroaggregated random shape of 99mTc-MAA particles (27). The embolization effect of the much larger number of 90Y-microspheres may result in flow alterations that alter the distribution of the particles.
Third, the use of bremsstrahlung 90Y SPECT after treatment leads to a degree of blurring and quantitative uncertainty. To overcome this methodologic problem, we chose to evaluate larger segment-based volumes instead of using a voxel-based analysis. The uncertainty in stochastic effects that are responsible for the measurement accuracy is far less in larger volumes. Larger VOIs were also useful to overcome quantification errors caused by coregistration artifacts. As a consequence, we were able to study quantitative distribution differences on only a segmental level. Any existing disagreement will be underestimated, since differences in certain areas of the segment may level out differences in other areas. One way to improve part of the methodology may be the use of 90Y PET imaging instead of SPECT. PET facilitates more accurate quantification based on improved spatial resolution and may aid in the study of distribution on a subsegmental or tumor level (30–32).
And lastly, the studied population was a heterogeneous group with regard to the histopathology of the primary tumor. We could not differentiate whether histopathologic features had any influence on the results. Remarkable was the finding that a larger degree of tumor involvement was associated with a better agreement between 99mTc-MAA and 90Y. This finding may be due to a lesser degree of random distribution of activity.
On one hand, the results lead to concerns about the agreement between 99mTc-MAA and 90Y distribution and, consequently, the validity of activity calculation methods that are based on the assumption of agreement, such as the partition method. On the other hand, it has been shown that, although based on the false assumption of agreement, the partition method still offers huge advantages over existing methods with regard to the prediction of treatment outcome and individualized treatment planning. This may be explained by the fact that the partition method is prescribed only for patients with a limited number of hypervascular tumors. Our results show that the agreement between 99mTc-MAA and 90Y distribution in segments with a high tumor involvement is considerably better. Moreover, these patients are generally treated superselectively, beyond the major bifurcation of the proper hepatic artery, and also with fewer embolic glass microspheres. All these factors seem to contribute to a better agreement between 99mTc-MAA and 90Y distribution.
On the basis of the current study, the following 2 recommendations may lead to optimization of the predictive value of a pretreatment scout dose. First, the catheter tip should be placed in exactly the same position during both procedures, possibly augmented by the use of catheters that center the tip within the lumen; second, for each administration, the catheter tip should be placed distal to major bifurcations, with selective administration in each branch to prevent preferential flow. In patients with large hypervascular tumors limited to a single lobe, the agreement between 99mTc-MAA and 90Y distribution is best, but caution should be taken in patients with multiple small liver metastases in both lobes, especially when an embolic effect of the microspheres is anticipated. The latter could be the case in patients who are treated with high-dosage resin microspheres and have small livers and prior treatments with antiangiogenic drugs such as bevacizumab.
To overcome the limitations of 99mTc-MAA as a scout dose, our group has developed a new generation of microspheres for multimodality image-guided radioembolization: 166Ho-poly(l-lactic acid) microspheres (33). The radioisotope 166Ho is embedded in microspheres of poly(l-lactic acid). It emits β-radiation (half-life, 26.8 h; maximum energy, 1.77 and 1.85 MeV) and γ-radiation (γ-energy, 80.6 keV) and is paramagnetic, because the element holmium is chemically part of the lanthanide group, like gadolinium (34–36). The microspheres can be visualized in vivo with several clinical imaging modalities, including SPECT and MR imaging (14,37,38). The particles used for pretreatment evaluation and actual treatment are exactly the same, and SPECT and MR imaging are used to combine high sensitivity with high spatial-temporal resolution and superior soft-tissue contrast to optimize dosimetry before and after treatment. The performance of 166Ho-microspheres as a scout dose to predict distribution of the therapy dose is currently under investigation.
CONCLUSION
Individualized treatment planning methods may be used for optimized safety and efficacy of radioembolization treatments. By definition, these methods are based on predictive scout dose distribution within the target volume. The limited agreement between 99mTc-MAA and 90Y distribution in current clinical practice raises concern about the validity of these methods. Care should be taken to use proper administration techniques to overcome this limitation.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. No potential conflict of interest relevant to this article was reported.
Footnotes
Published online Jun. 7, 2013.
- © 2013 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication November 23, 2012.
- Accepted for publication March 12, 2013.