|
|
||||||||
Basic Science Investigations |
1 Division of Nuclear Medicine, Department of Radiology, University of Michigan Medical Center, Ann Arbor, Michigan
2 Department of Radiology, University of Michigan Medical Center, Ann Arbor, Michigan
3 Division of Hematology and Oncology, Department of Internal Medicine, University of Michigan Medical Center, Ann Arbor, Michigan
| ABSTRACT |
|---|
|
|
|---|
Key Words: tumor dosimetry 131I radioimmunotherapy lymphoma SPECT registration
| INTRODUCTION |
|---|
|
|
|---|
To date, clinical correlates to the degree of response in 131I-tositumomab therapy have not been easy to find. A logistic regression analysis that included patients from 6 clinical trials (phase I, II, or III) covering the full range of grades of NHL found that elevated lactate dehydrogenase (LDH), relapse after prior radiotherapy, and high tumor burden correlated with a lesser degree of response (4). Tumor radiation absorbed dose was not tested as an explanatory variable in that study.
Initial results for estimates of tumor radiation absorbed dose by a hybrid pretherapy conjugate-view intratherapy SPECT method were reported previously for 33 of the 52 patients scanned by SPECT (5). In addition, correlation of the absorbed dose for each individual tumor with the degree of response was reported for 20 of these 33 patients (6). However, because most of these 20 patients had >1 tumor and each tumor was related to a response, the statistical test of the significance of the relation between tumor radiation dose and response relied on the correctness of a complicated correction for interdependence in the software.
The correlation of tumor volume reduction 12 wk after therapy with estimated radiation dose was also performed for partial response (PR) patients and previously published (7). A negative slope was found in the plot of volume reduction of composite tumor versus dose estimate by pretherapy conjugate views, which is counterintuitive. In contrast, the experimental points with dose estimates from the hybrid method were well fit by sigmoid shapes, with volume reduction increasing with dose, as expected. The P value representing the significance of the relationship relative to no doseresponse relationship approached statistical significance for all individual tumors (P = 0.06) and reached statistical significance for the subset of individual tumors that had small pretherapy volume (P = 0.03) (7).
In this article, we compare activity estimates from the hybrid method with those from the conjugate-view method, report dosimetric results from the hybrid method for 19 additional patients, and expand our examination of the results with both methods. Importantly, we analyze the same tumor CT volumes by both methods, even though the volumes are grouped into fewer, larger (composite) tumors by the conjugate-view method. In the statistical analysis, we use the patients average tumor radiation dose to provide more robust testing of statistical significance of the doseresponse relationship. In addition, because dependence of tumor activity on time is an important dosimetric consideration, we intercompare SPECT timeactivity curves during evaluation for 6 individual tumors of 1 patient and compare their sum with the curve for the corresponding composite tumor measured by conjugate views.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Baseline CT, Administered Activity, and Informed Consent
Patients underwent a baseline CT scan, usually about 1 wk before the start of any nuclear medicine procedures and usually after receiving oral and intravenous contrast. They underwent pretherapy conjugate-view imaging, after a predose of tositumomab followed by a tracer dose of 131I-tositumomab antibody. The tracer evaluation and the therapy were as described (5). Patients gave their separate written informed consent for all SPECT imaging that was not part of the normal 131I-tositumomab protocol. This imaging received separate University of Michigan Internal Review Board approval.
Hybrid-SPECT Dosimetry
Except as described below, the procedure for the imaging and the methods for estimation of tumor radiation dose with conjugate views and with the hybrid method were as described (5,6). The hybrid procedure assumed that the shape of the timeactivity curve for each individual tumor during therapy was the same as that for the associated composite tumor during evaluation. Our previous measurements of timeactivity curves during therapy for 9 individual tumors in a patient showed 8 of the 9 curves to be similar, so the measurement of individual shapes did not seem essential (5). Moreover, in another patient, the shape of the pretherapy conjugate-view curve for the composite tumor was very similar to that for the intratherapy SPECT curve for the sum of 4 individual tumors associated with that composite tumor (5). Therefore, we chose to continue using the pretherapy conjugate-view curve shape to relieve us of the necessity of acquiring and analyzing a set of SPECT measurements after therapy.
The hybrid-SPECT procedure uses a total activity recovery coefficient, defined as the activity measured divided by the true activity. The coefficient is determined as a function of volume from a phantom study with spheres of various volumes (8). The total-activity recovery-coefficient correction factor, defined as the inverse of the coefficient, multiplies the initial estimate of tumor activity. Previously (5,6,9), when a tumor had a large volume, a value of 1 was used for the correction factor. In this study, the value read off the curve fitted to the phantom data is used regardless of the tumor volume. Monte Carlo simulation has indicated the new procedure to be correct (Y. Dewaraja, unpublished data, September 1999). Also, the new procedure yielded better agreement for an experimental liver activity measurement in an anthropomorphic phantom: With the factor fixed at 1, the activity estimate was 59.5% high, whereas, with the fitted value, the estimate was correct within 5%. The fitted curve for the correction factor for head 3 of the Prism 3000XP camera (Marconi Medical Systems, Cleveland, OH) (8) is shown in Figure 1. The dose estimates for all large individual tumors in the patients who had been evaluated previously were revised to accommodate the new correction.
|
![]() |
Patient Response
For each patient, a posttherapy CT scan assessed the response to therapy at 7, 13, and 26 wk after tracer dose and every 3 mo thereafter, until disease progression. A complete response (CR) was complete disappearance of all measurable and evaluable disease for at least 4 wk, and a PR was
50% reduction in the sum of the products of the longest perpendicular diameters of measurable lesions for at least 4 wk, with no new lesions.
Pretherapy SPECT Time Series
Patient 64 was scanned with SPECT 6 times after the tositumomab tracer infusion that was administered for evaluation purposes. Normal conjugate-view scans were obtained either immediately before or immediately after the SPECT scans. The time period covered was 163 h. The time series provided an opportunity to compare the measured SPECT timeactivity curve for the sum of 6 individual tumors to the measured conjugate-view timeactivity curve for the composite tumor and also to compare the SPECT timeactivity curves for the individual tumors among themselves. Unlike in the generation of the dose estimates, scatter correction was not used to generate the values for the timeactivity plots. This was because the individual tumor dose estimates during evaluation were quite noisy and that noise was worse with scatter compensation, making the plots harder to interpret.
Data Analysis
Before calculating statistical significance, the dose values for PR patients and for CR patients with both conjugate-view dose estimation and with hybrid-SPECT dose estimation were checked to see if they were normally distributed. They were not, and so all of the dose data was transformed by taking the logarithm to the base 10. This change shifted the dose values toward a normal distribution and then significance was calculated. A 2-sample Student t test available within the application Excel (Microsoft, Redmond, WA) or from the Statistical Analysis System ([SAS] SAS, Cary, NC) was used to compare the dose means from different populations. Significance for the difference between the mean value for the patients average dose when including only axillary tumors compared with that when including only chest, abdominal, and pelvic tumors was calculated by SAS with a 2-tailed test because there was no reason to expect a difference. No declaration of variance equality or inequality was required by the software.
Significance for the difference between the mean value for the patients average dose over evaluated tumors for CR patients versus that for PR patients was calculated by Excel with a 1-tailed test because it was expected that the mean value for CR patients would be greater than that for PR patients. Equal variance was assumed for the 2 types of responders with conjugate views because the measured variances differed by a factor of <2. Unequal variance was assumed with the hybrid method because the measured values differed by a factor of >2.
Total tumor burden and elevated LDH level were 2-valued variables that were associated with degree of response over a large number of patients from multiple trials with tositumomab (4). A multiple logistic regression analysis was performed for the 47 patients above to see if combining radiation dose with those 2 variables produced a correlation with response (10). SAS was used.
In addition, a pseudopredictive analysis was performed for the patients average tumor dose and response. It served 2 purposes. The first was to characterize the amount of dose overlap between PR and CR patientsthat is, if there were no overlap, the predictive values would be 100%. The second was to provide a set of pseudopredictive values for previously untreated patients. A justification is given in Koral et al. (6). A similar analysis was performed for each of the two 2-valued variables mentioned above for comparison to the results with dose.
To perform the pseudopredictive analysis, we defined the positive predictive value (PPV) to be the fraction of all patients who achieved a CR who were predicted to do so by a high average tumor dose, low tumor burden, or low LDH level. Also, we defined the negative predictive value (NPV) to be the fraction of all patients who achieved a PR who were predicted to do so by a low average tumor dose, high tumor burden, or high LDH level. Finally, we defined the accuracy (A) to be the fraction of all patients who achieved the response predicted. Also, for prediction by dose, we set the dividing line between a high and low dose to be the average value for the mean of the dose for CR patients and the mean of the dose for PR patients, as is standard.
We also made the following auxiliary definitions. We defined TP to be the number of patients predicted to achieve a CR who did; FP to be the number of patients predicted to achieve a CR who did not; TN to be the number of patients predicted to achieve a PR who did; and FN to be the number of patients predicted to achieve a PR who did not. Then, PPV, NPV, and A were evaluated as follows:
![]() |
![]() |
![]() |
| RESULTS |
|---|
|
|
|---|
To confirm an earlier finding that axillary tumors had low dose values with the hybrid method (5), the mean value over 8 patients for the patients average tumor dose using only axillary tumors was compared with the mean value over 47 patients using only chest, abdominal, and pelvic tumors. The mean for axillary tumors (162 cGy) was less than one fourth the mean for nonaxillary tumors (624 cGy). Variance about the mean was 774 cGy2 for the axillary tumors and 183,612 cGy2 for the nonaxillary tumors. The difference in the 2 means was highly significant statistically (P < 0.0001).
Because dose estimates for axillary tumors were uniformly low compared with other tumors, and because all patients with an analyzed axillary tumor had a CR, the axillary tumors were eliminated from the analysis of radiation dose versus response to have the best chance for obtaining a statistically significant difference between CRs and PRs based on radiation dose. If the axillary tumors had been included, they would have decreased the patients average tumor dose for 3 CRs, added 5 new CRs with low values for patients average tumor dose, and so the separation of all CRs versus all PRs would have been less.
The axillary tumors were not included as follows: The patients for whom we had evaluated only 1 or more axillary tumors by the hybrid method were excluded. The number available for the response analysis thus dropped from 52 to 47. Also, for the 3 patients for whom we had evaluated both 1 or more axillary tumors and 1 or more abdominal or pelvic tumors, the patients average tumor dose was computed from only his or her abdominal or pelvic tumors.
For the 47 patients, total tumor burden and LDH level are listed in columns 2 and 3 of Table 1. Each patients average tumor dose, as evaluated by the hybrid method, is shown in column 5. The number of individual tumors that were averaged is shown in column 4. The dose values range from a minimum of 125 cGy to a maximum of 25.4 Gy. The average dose estimate for each patient using pretherapy conjugate views also was calculated excluding his or her axillary tumors and is presented in column 7, with the number of composite tumors that were averaged shown in column 6. The dose values range from a minimum of 222 cGy to a maximum of 76.6 Gy. The patients response to therapy is listed in column 8.
|
|
The calculated negative and positive pseudopredictive values and accuracy for pretherapy conjugate views and for hybrid SPECT, as well as for LDH and for tumor burden, are compared in Table 2. All 3 values were greater with the hybrid method than with the conjugate-view method. For accuracy alone, in order, from the poorest to the best, were conjugate-view dose, hybrid-SPECT dose, tumor burden, and LDH.
|
|
The 3 variables in the multiple logistic regression model did not predict the degree of response in a statistically significant manner with either dose-estimation method. The P value was 0.51 with the hybrid-SPECT method and 0.58 with the conjugate-view method.
For patient 64, the plot of the percentage injected dose versus time after infusion for the sum over 6 individual tumors from pretherapy SPECT and for pretherapy conjugate views for the corresponding composite tumor each had about the same area under the curve. This result is the same as the published result for patient 43 (5).
The shapes of the timeactivity curve for the 6 individual tumors of patient 64 were similar. This result is also the same as the published result for 8 of 9 tumors of patient 43 (5). Looking at the fine differences between the curves for individual tumors for patient 64, 2 of the tumors showed at least some rise at 19 h compared with 0 h. The other 4 were characterized by constant or decreasing values at 19 h compared with 0 h. The shape for 1 of the 4 tumors was close to that for a decaying exponential. The 2 tumors that had an activity rise at 19 h had slightly higher values at the last time point (163 h) than the values for the other 4 tumors.
| DISCUSSION |
|---|
|
|
|---|
The difference in radiation dose estimated by the hybrid method between axillary and other tumors is not caused by a consistent error in technique as far as our Monte Carlo investigations indicate. Therefore, we conjecture that it is due to lower blood flow to the axillary tumors or to some other host factor.
We infer that early prediction of response may be possible for a subset of all patients with evaluated chest, abdominal, or pelvic tumors (in this sample, it would be for 16 of a total of 47 patients, or 34%)that is, if the mean radiation dose is large, it could be predicted with confidence that the patient will have a CR. If the dose is smaller, then a prediction cannot, at present, be made (in this sample, 8 had a PR but even more, 23, achieved a CR).
Both dose methods yielded a higher mean tumor radiation dose, on average, for a patient who went on to a CR compared with a patient who went on to a PR. This result is in agreement with the trend for a correspondence between higher radiation dose and better response of refractory patients that was observed when the larger doses in myeloablative therapy (12) are combined with the smaller doses in nonmyeloablative 131I-tositumomab therapy (13). In both studies, the dose-estimation method was based on pretherapy conjugate views. The patients had a CR, PR, minor response, or disease progression. No statistical test was applied to the combined data (K.F. Koral and R.L. Wahl, unpublished data, June 1995).
Patients with a low average tumor dose proceeded to either a PR or a CR. The difference in response for these patients may be attributable to a difference in the effect of the unlabeled antibody, even though the same amount was infused for all patients. However, in comparable patients, there is as yet no experimental evidence for such a difference in unlabeled antibody effect with equal predose.
A different possibility is that the accuracy of our method for dose estimation is not high enough or that the sampling of tumors in the chest, abdomen, and pelvis was somehow insufficient for our purpose. In the hybrid method, the conversion factor to produce the activity estimate from total reconstructed counts varies with the ratio of the background to the tumor activity concentration. The estimate for this ratio in patients varies widely between different tumors and is cause for concern that some of the individual dose values may be inaccurate. Use of a super-high-resolution reconstruction algorithm in the future could possibly eliminate this problem by making the counts-to-activity conversion factor independent of the level of background activity. The sampling problem could be reduced by the use of a SPECT camera with a large field of view.
Another possibility is that the follicular class of NHL combines patients whose disease is diverse. This has been shown to be the case in the diffuse large B-cell class of NHL. In that class, gene expression profiling has identified 2 subgroups with significantly different overall survival (14). If such profiling could have separated our patients into subgroups, the dose values might have been more predictive.
The lack of higher P values is in contrast to recent results with radiopharmaceutical therapy of neuroblastoma with 131I-methyliodobenzylguanidine (MIBG). In that disease-therapy system, with 16 PRs, 3 mixed responses, 14 cases of stable disease, and 9 cases of progressive disease, tumor self-absorbed radiation dose was a significant predictor of degree of response (P < 0.01) (15). Intratherapy conjugate views were used and a single tumor was selected for each patient on the basis of MIBG uptake.
The lack of a bigger difference between the hybrid method and conjugate views in predicting response is in contrast to the large difference in the results for the 2 methods in the study of volume-reduction versus dose (7) that was cited previously. The reason may be that conglomeration of the dose and the volume reduction of several individual tumors into a composite tumor, as was done in the conjugate-view part of that study, may be very suboptimal in predicting the volume reduction from dose. This, in turn, would result because the dose and the volume from 2 or more tumors, which are perhaps unlike, are first both combined and then the ratio is effectively taken to get the dose for the composite tumor, leading to a value that is perhaps not representative of the true dose for any of the constituent tumors. Similarly, the volumes before and after therapy are first both combined and then the ratio is taken to get the volume reduction for the composite tumor. This procedure could lead to calculated values of volume reduction that are unrealistic. Because both variables can have large errors, an incorrect correspondence is not surprising.
It appears that taking into account the shape differences among the timeactivity curves for individual tumors would affect the cumulated activity for the individual tumors to only a small extent (on the order of 10%). Therefore, from the results we have on only 2 patients, there is no evidence that measuring such shape differences for individual tumors would greatly affect the observed degree of separation for patient average radiation dose between CRs and PRs.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Kenneth F. Koral, PhD, Division of Nuclear Medicine, University of Michigan, B1G412 UH, Ann Arbor, MI 48109-0028.
E-mail: kenkoral{at}umich.edu
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E. Mah and K. M. Spicer Comparison of Medium- and High-Energy Collimators for 131I-Tositumomab Dosimetry J. Nucl. Med. Technol., September 1, 2007; 35(3): 148 - 153. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. K. Tan, R. W. Wassenaar, and W. Zeng Collimator Selection, Acquisition Speed, and Visual Assessment of 131I-Tositumomab Biodistribution in a Phantom Model J. Nucl. Med. Technol., December 1, 2006; 34(4): 224 - 227. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. K. Dewaraja, S. J. Wilderman, M. Ljungberg, K. F. Koral, K. Zasadny, and M. S. Kaminiski Accurate Dosimetry in 131I Radionuclide Therapy Using Patient-Specific, 3-Dimensional Methods for SPECT Reconstruction and Absorbed Dose Calculation J. Nucl. Med., May 1, 2005; 46(5): 840 - 849. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Goldenberg and R. M. Sharkey Radioimmunotherapy of Non-Hodgkin's Lymphoma Revisited J. Nucl. Med., February 1, 2005; 46(2): 383 - 384. [Full Text] [PDF] |
||||
![]() |
R. L. Wahl Tositumomab and 131I Therapy in Non-Hodgkin's Lymphoma J. Nucl. Med., January 1, 2005; 46(1_suppl): 128S - 140S. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Sgouros, K. S. Kolbert, A. Sheikh, K. S. Pentlow, E. F. Mun, A. Barth, R. J. Robbins, and S. M. Larson Patient-Specific Dosimetry for 131I Thyroid Cancer Therapy Using 124I PET and 3-Dimensional-Internal Dosimetry (3D-ID) Software J. Nucl. Med., August 1, 2004; 45(8): 1366 - 1372. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. F. Koral, K. Zasadny, R. L. Wahl, and R. M. Sharkey Bias Reduction in Correlation of Radiation-Absorbed Dose with Response J. Nucl. Med., July 1, 2004; 45(7): 1271 - 1272. [Full Text] [PDF] |
||||
![]() |
K. E. Britton Radioimmunotherapy of Non-Hodgkin's Lymphoma J. Nucl. Med., May 1, 2004; 45(5): 924 - 925. [Full Text] [PDF] |
||||
![]() |
R. M. Sharkey, A. Brenner, J. Burton, G. Hajjar, S. P. Toder, A. Alavi, A. Matthies, D. E. Tsai, S. J. Schuster, E. A. Stadtmauer, et al. Radioimmunotherapy of Non-Hodgkin's Lymphoma with 90Y-DOTA Humanized Anti-CD22 IgG (90Y-Epratuzumab): Do Tumor Targeting and Dosimetry Predict Therapeutic Response? J. Nucl. Med., December 1, 2003; 44(12): 2000 - 2018. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. F. Koral, M. S. Kaminski, R. L. Wahl, G. Sgouros, S. Squeri, A. M. Ballanrud, K. Kolbert, J. B. Teitcher, K. S. Panageas, R. D. Finn, et al. Correlation of Tumor Radiation-Absorbed Dose with Response Is Easier to Find in Previously Untreated Patients J. Nucl. Med., September 1, 2003; 44(9): 1541 - 1543. [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |