Abstract
To determine whether the therapeutic effect of 90Y-ibritumomab might be enhanced by a full course of rituximab followed by single dose of 90Y-ibritumomab, the trial included pre- and post-rituximab treatment imaging with 111In-ibritumomab and blood pharmacokinetics. Comparison of the pre- and post-rituximab imaging and blood data allowed for the assessment of impact of rituximab on 90Y-ibritumomab dosimetry. Methods: Seventeen patients with relapsed B cell non-Hodgkin lymphoma first received 250 mg/m2 of rituximab plus 185 MBq of 111In-ibritumomab for initial dosimetry evaluation. In weeks 2–4, patients received 3 weekly 375 mg/m2 doses of rituximab. In week 6, patients received a 250 mg/m2 dose of rituximab plus 185 MBq of 111In-ibritumomab for a second dosimetry evaluation. Five sequential, whole-body γ-camera images were acquired after the 111In-ibritumomab injection. Calculated radiation doses were based on patient-specific organ masses. For each patient, paired comparison of radiation doses before and after rituximab treatment was performed. Paired comparison of residence times for spleen and tumor was also performed. Results: Before rituximab treatment, the median radiation dose (mGy/MBq) was 0.48 (range, 0.24–0.86) for total body, 3.7 (range, 2.1–11.6) for liver, 6.1 (range, 1.8–17.8) for spleen, 3.3 (range, 2.0–4.7) for kidneys, 2.4 (range, 1.3–3.7) for heart wall, 1.1 (range, 0.4–2.3) for lungs, 0.79 (range, 0.32–1.22) for marrow from blood, and 18.1 (range, 4.7–98.9) for tumor. Paired comparisons were performed in 16 patients only because human antimurine antibody developed in 1 patient. The median change was 0.007 mGy/MBq for body, −0.14 mGy/MBq for liver, −0.31 mGy/MBq for kidneys, 0.38 mGy/MBq for heart wall, −0.17 mGy/MBq for lungs, and 0.046 mGy/MBq for marrow from blood. The median change in residence time was −0.92 h for spleen and −0.24 h for tumor. The changes were statistically insignificant for total body, liver, kidneys, lungs, and marrow from blood. The median residence times, or mGy/MBq if there were no volume changes, decreased 24% for spleen (P = 0.0005) and 28% for tumor (P = 0.005). The median radiation dose to heart wall increased 16%, which was statistically significant (P = 0.002). Conclusion: Changes in 90Y-ibritumomab dosimetry after 4 wk of rituximab treatment were not significant for most organs, except for the heart wall. The reduction of spleen and tumor residence times is more likely to be due to the therapeutic effects of rituximab.
- monoclonal antibodies
- radiobiology/dosimetry
- radionuclide therapy
- 90Y-ibritumomab
- rituximab
- Zevalin
- dosimetry
The non-Hodgkin lymphomas (NHL) are a diverse group of lymphoid neoplasms that collectively rank fifth in cancer incidence and mortality in the United States (1). Chemotherapy treatment with or without radiotherapy can lead to remissions; however, virtually all treated patients with low-grade NHL repeatedly relapse, with shorter remissions after each subsequent course of chemotherapy (2). Additionally, no chemotherapeutic regimen has been shown to be superior to another in prolonging survival in low-grade NHL (3). Ibritumomab is a murine IgGl κ-monoclonal antibody that is linked to 111In for imaging and 90Y for therapy via the chelator tiuxetan (isothiocyanatobenzyl M–diethylenetriaminepentaacetic acid). Ibritumomab targets the CD20 antigen, which is present on 95% of B-cell lymphomas (4). 90Y-Ibritumomab (Zevalin; BIOGEN/IDEC Pharmaceuticals) yields an overall response rate of 73%, with a response duration of 12 mo when administered as a single nonmyeloablative fraction in relapsed follicular lymphoma (5).
The standard dosing of 90Y-ibritumomab is 11.1 or 14.8 MBq/kg (0.3 or 0.4 mCi/kg, respectively) for platelet count of 100,000–149,000/mm3 or greater than or equal to 150,000/mm3, respectively. The maximal administration dose is 1,184 MBq (32 mCi), regardless of patient body weight (6). The radiolabeled antibody is given after 250 mg/m2 of rituximab (Rituxan; BIOGEN/IDEC Pharmaceuticals). While the dose administration of 90Y-ibritumomab is mainly limited by myelotoxicity, unlabeled anti-CD20 mAb, rituximab, does not have a severe myelosuppressive effect (7,8). It was hypothesized that the therapeutic response might be further enhanced by a full course of rituximab, followed by single dose of 90Y-ibritumomab, because these 2 treatments do not have spatial additive adverse effects on the same normal organ. However, there was concern that an additional large amount of unlabeled anti-CD20 mAb, rituximab, might adversely alter the pharmacokinetics and dosimetry of 90Y-ibritumomab.
To determine whether the therapeutic effect of 90Y-ibritumomab might be enhanced by a full course of rituximab, followed by single dose of 90Y-ibritumomab, the trial included pre– and post–rituximab treatment imaging with 111In-ibritumomab and blood pharmacokinetics. Comparison of paired quantitative imaging and blood data before and after a full course of rituximab treatment allowed for the assessment of changes in 90Y-ibritumomab dosimetry as a result of rituximab treatment. To further improve accuracy, dosimetry results were determined using patient-specific organ masses.
MATERIALS AND METHODS
Patients and Study Design
Seventeen patients with histologically confirmed, follicular, or transformed CD20-positive B-cell NHL who had failed at least 1 chemotherapy regimen were enrolled in the present study. No anticancer treatment was allowed for 3 wk before study initiation (6 wk if rituximab, nitrosourea, fludarabine, or mitomycin C), and patients had to have fully recovered from all toxicities associated with prior treatment. Bone marrow involvement with lymphoma was less than 25% from bilateral bone marrow biopsy and aspirate evaluation. Exclusion criteria included prior external radiotherapy to greater than 25% of active bone marrow and the presence of human antimurine antibody (HAMA) reactivity with prior exposure to murine antibodies. Fifteen patients were enrolled at University of Alabama at Birmingham, 1 patient was enrolled at Stanford University, and 1 patient was enrolled at Mayo Clinic. All patients signed written informed consent approved by the corresponding institutional review board. All patients received a first dose of 250 mg of rituximab per square meter, plus 185 MBq (5 mCi) of 111In-ibritumomab for the initial dosimetric evaluation. In weeks 2–4, patients received 3 weekly 375 mg/m2 doses of rituximab. In week 6, patients received a 250 mg/m2 dose of rituximab, plus 185 MBq (5 mCi) of 111In-ibritumomab for a second dosimetric evaluation (Fig. 1). Including 1 preload rituximab dose, patients received a total of 1,375 mg of rituximab per square meter in 4 wk before the second 111In-ibritumomab dose injection, similar to a full course of rituximab (4 × 375 mg/m2). All patients had measurable lesions assessed with CT at baseline, about 2 wk before the first injection of 111In-ibritumomab. CT slices were 3 or 5 mm in thickness. 90Y-ibritumomab dosimetry was derived from 111In-ibritumomab imaging and paired dosimetric results before and after 3 weekly doses of 375 mg/m2.
Rituximab levels were compared in 16 of 17 patients. Paired dosimetric results could not be compared for patient 10 because HAMA developed in the patient after 1,375 mg/m2 doses of rituximab in 4 wk. Among the 16 patients, patient 7 had a large body size and weighed 159 kg. Although parts of this patient's arms were outside the γ-camera field of view, all organs were inside the field of view; thus, dose estimates for this patient were not affected.
Rituximab and 111In/90Y-Ibritumomab
Rituximab is a chimeric IgG1 κ-monoclonal antibody. It specifically recognizes the CD20 antigen expressed on normal B cells and most malignant B-cell lymphomas. Rituximab is licensed in the United States as therapy for patients with relapsed or refractory low-grade or follicular, CD20-positive, B-cell NHL. A full course of rituximab consists of 4 weekly doses of 375 mg/m2.
Ibritumomab is a murine IgGl κ-monoclonal antibody that is linked to 111In (for imaging) and 90Y (for therapy) via the chelator tiuxetan (isothiocyanatobenzyl M–diethylenetriaminepentaacetic acid). 90Y-ibritumomab is licensed for the treatment of patients with relapsed or refractory low-grade, follicular, or transformed B-cell NHL, including patients with rituximab-refractory follicular NHL. 90Y-ibritumomab is administered in 2 steps: step 1 includes a single 250 mg/m2 infusion of rituximab before a fixed dose of 185 MBq (5 mCi; 1.6 mg of antibody). Step 2 follows step 1 by 7–9 d and consists of a second 250 mg/m2 infusion of rituximab before 14.8 MBq (0.4 mCi/kg) of 90Y-ibritumomab per kilogram.
Quantification of 90Y/111In in Blood
Whole-blood samples were collected at 30 min and 4, 24, 72, and 144 h after the administration of 111In-ibritumomab. 111In activity was measured using a γ-well counter, and counts from the sample were converted to activity by co-counting samples with a reference source (9). The reference source containing 19 kBq (0.5 μCi) of 111In in 1 mL was pipetted from a calibrated stock solution with a known concentration. Activities of 111In-ibritumomab were corrected for decay of the radioisotope from the time of drug administration to the time of measurement. 90Y-ibritumomab activity in the whole blood was derived from 111In activity adjusted by a small difference in physical decay half-life. 90Y-ibritumomab concentration was expressed as percentage injected dose per milliliter (%ID/mL) and analyzed using a biexponential clearance (α and β) model. The cumulated activity or residence time was determined with the biexponential clearance parameters (9). Pharmacokinetic results were determined for the first dosimetry study in all 17 patients. Pharmacokinetic results were determined in 15 of 17 patients for the second dosimetry study because HAMA developed in 1 patient after 3 weekly 375 mg/m2 doses of rituximab, and blood samples were lost in 1 patient.
Quantitative Imaging
Quantitative imaging for 111In has been described previously (9,10). Planar conjugate whole-body images were acquired with a dual-detector γ-camera interfaced to a nuclear medicine computer system. Medium-energy collimators were used to image 111In with energy windows centered at 171 and 245 keV (15% width). Transmission scan images were obtained using a 57Co sheet source containing about 370 MBq (10 mCi). The same medium-energy collimators were used to image 57Co with an energy window centered at 122 keV (15% width). Planar conjugate views of the whole body were acquired immediately after the administration of 111In-ibritumomab and at 4, 24, 72, and 144–168 h. Whole-body images were acquired in a 256 × 1,024 word matrix. Each patient's position on the imaging table and vertical positions of the camera detectors were recorded in the first imaging session and were used throughout the sequential imaging studies for reproducible detector-to-patient positioning (11,12). A 10-mL 111In reference source containing 1.9 MBq (50 μCi) was placed on the imaging table at least 10 cm away from the patient's feet. This calibrated source was used to correct for potential system drift and to convert image region-of-interest (ROI) counts to units of radioactivity adjusted with experimentally measured volume effect.
Before processing images, the operator reviewed 111In-ibritumomab images and CT images with the physician who was familiar with the patient's medical record and CT images for determining organ and tumor ROIs. Major organs that were visible above body background tissue after clearance of the blood pool were contoured, except for the bladder. Although tumors were visible for each patient, tumors were quantified only if they met the following criteria for adequate accuracy: at least 1 g in mass, tumor-to-background pixel count ratio of at least 1.5, and clear tumor ROI boundary. Counts in selected ROIs were corrected for background. The thickness of the background region was equivalent, after accounting for photon attenuation, to the overlapped background tissues in the organ from the anterior and posterior views.
The importance of attenuation correction in quantitative imaging has long been recognized (13,14). We determined activities in the liver and lungs using the geometric-mean quantification.Eq. 1where ISource,anterior and ISource,posterior represent the background-corrected counts per second in the source ROI in anterior and posterior views, respectively. The attenuation-correction factor, ACF, was determined byEq. 2where Npt and Nno pt represent the liver or lung ROI counts in the 57Co transmission images with and without the patient, respectively. Equation 2 corrects for the energy differences between 57Co and 111In.
The values for μ57Co,liver or lungs and μ111In,liver or lungs were determined experimentally using large-volume phantoms (9). The value of f, a self-attenuation factor for source thickness, was also experimentally measured instead of calculated using the formula μst/sinh(μst/2) (13,15). c is the camera system calibration factor (counts/s/activity) for that source adjusted with experimentally measured volume effect.
For a source organ that can be clearly visualized only in a single view, previous study suggested that source quantification using a single view is more suitable than using conjugated views (14,16). Kidney, spleen, and tumor activity were quantified using an effective point-source method:Eq. 3where ISource is background counts corrected, counts per second in the source ROI of a single-view image. μeff is the effective linear attenuation coefficient determined using phantoms of 150, 50, and 10 mL. The source depth d is determined from CT images. c is the camera system calibration factor (counts/s/activity) for that source. The depth of kidneys, spleen, and tumors from the body surface were measured from CT images. Phantom sizes and depths were selected to mimic sizes and depths for the liver, spleen, kidneys, and tumors in patients (9). The self-attenuation factor for source thickness is not explicitly included in Equation 2 because our phantom data suggested that the self-attenuation factor was close to 1 for 150- and 50-mL sources.
The uptake of 90Y/111In in organs and tumors was expressed as percentage of injected dose at various imaging time points. Cumulated activity (and residence time) and biologic clearance half-life were determined by fitting the uptake data with a monoexponential curve. If fit of a monoexponential curve was not possible as uptake data continuously increased during the period of sequential imaging, cumulated activity was determined using the trapezoid method, and a conservative estimation of the tail was used with a clearance rate of the physical half-life. The cumulated activity or residence time of 90Y was determined from the biodistribution of 111In and adjusted for the small difference in physical half-life between 111In and 90Y.
Organ and Tumor Volumetry Using Patient-Specific Measurements
Baseline CT images before the first injection of 185 MBq (5 mCi) of 111In-ibritumomab were used to determine a patient-specific organ mass to increase the accuracy of the dose estimate. Volume measurements were performed on CT Digital Imaging and Communications in Medicine images using the Eclipse 3-dimensional radiotherapy treatment planning system (Varian Medical System). The tumor, liver, spleen, kidney, and heart volumes were manually defined by individual slice contour drawings on 3- to 5-mm-thick slices. Lung volumes were determined using a contrast-difference edge-detection method, with manual correction for air in the main bronchi. A nominal tissue density of 1 g/cm3 was used for soft tissues, except for the lungs, for which we assumed the lung density to be 0.26 g/cm3 (17). The heart region was defined as extending below the level in which the pulmonary trunk branches into the left and right pulmonary arteries. The activity measured in the heart was assumed to be associated with circulating radiolabeled antibody rather than with uptake in the heart muscle and the mass of the heart wall was calculated.
This retrospective dosimetry analysis was part of a clinical trial to assess the response of the combination of full-course rituximab plus 90Y-ibritumomab. Posttreatment CT images were obtained after the completion of all treatments. CT images after the last rituximab dose and before the second dosimetric imaging study were not available. Volumes in the liver, lungs, kidneys, and heart were assumed to be unchanged before and after a full course of rituximab. However, we were not able to assume that the volume of tumors or spleens remained unchanged after a full course of rituximab. It is known that the full course of rituximab, a 4-dose schema approved by the Food and Drug Administration, induced an overall response rate of 48% (18). To assess if spleen volume remained unchanged, the product of the spleen long and short axes (19) in 111In-ibritumomab images were measured before and after rituximab treatment.
Radiation Dosimetry with Patient-Specific Organ Mass
Residence times for the liver, spleen, kidneys, lungs, heart, tumor, marrow, and total body were determined from quantitative imaging or counting as described above. The residence time for the remainder of the body was defined as residence time of the total body minus residence times of the liver, spleen, kidneys, lungs, heart, and marrow. These residence times were input into an OLINDA-EXM program for radiation dose calculation (20). The voiding bladder module in the OLINDA-EXM program was used, assuming a voiding interval of 3 h during the day and an interval of 6 h at night (21). The masses of liver, spleen, kidneys, lungs, heart wall, and body of individual patients were used to modify phantom masses in the OLINDA-EXM program. The mass of the individual heart wall was calculated assuming a ratio of 316 to 770 between the heart wall and the heart (wall + content) (22).
To determine 90Y concentration in red marrow, the red marrow–to–blood concentration ratio was ascertained with a widely used method described by Sgouros (23). Although the residence time of red marrow for each patient was calculated assuming a reference mass of 1,120 g, the dosimetry result of red marrow from blood was still patient-specific for 90Y because the patient-specific marrow dose from β-emissions in the blood does not require an explicit estimate of marrow mass (24).
For the high-energy β-emitter 90Y, the absorbed fraction is less than 0.89 for a tumor mass that is less than 20 g (25). Tumor mass and residence time were input into the spheric model module of the OLINDA-EXM program for tumor dose calculation (20).
Statistics
Paired comparisons between tissue radiation dose estimates before and after a full course of rituximab were performed (26). Normal distribution assumption for the difference before and after was examined by skewness, kurtosis, and omnibus normality tests. If normality could not be rejected by each of these tests, results obtained from a paired t test were reported. If normality was rejected, results obtained from a nonparametric test (Wilcoxon signed rank) were reported. Differences in total body, kidneys, lungs, and marrow doses were analyzed with a t test. Differences in liver, spleen, heart wall, and tumor were analyzed with the Wilcoxon signed rank test. A comparison was considered statistically significant if the P value was less than 0.05.
RESULTS
Organ and Tumor Masses
The median and range for body weights and organ masses are listed in Table 1. Compared with the mass of the reference man's phantom, the observed mass could be 2.2 times higher for the body, 2.5 times lower for liver, 3.8 times higher for spleen, 1.8 times higher for kidneys, 1.9 times lower for lungs, and 2.8 times lower for heart wall (Table 1). Tumor masses ranged from 1 to 282 g (n = 28). In those 28 tumors, 9 tumors had masses between 1 and 10 g, 14 tumors had masses between 11 and 50 g, and 5 tumors had masses greater than 50 g.
Impact on 90Y-Ibritumomab Pharmacokinetics
Blood activities versus time fitted well with a biexponential model. The coefficient of determination (R2) for biexponential curve fitting was greater than or equal to 0.99 for all studies (n = 32), except for 1 (R2 = 0.97). Before rituximab treatment, the median peak concentration of 90Y-ibritumomab in the whole blood was 0.017 %ID/mL and ranged from 0.012 to 0.035 %ID/mL (n = 17). The median β- (slow) effective clearance phase was 31.9 h and ranged from 11.8 to 55.3 h. The median cumulated activity per unit injected dose or area under the curve was 0.0024 h/mL and ranged from 0.0007 to 0.0037 h/mL. The median radiation dose from blood to red marrow was 0.79 mGy/MBq and ranged from 0.32 to 1.22 mGy/MBq.
After 4 wk of rituximab treatment, the median difference in peak concentration of 90Y-ibritumomab in whole blood was −0.0015 %ID/mL and ranged from −0.005 to 0.0025 %ID/mL (n = 15). The median difference in the β- (slow) effective clearance phase was 1.0 h and ranged from −14.1 to 24.6 h. The median difference in cumulated activity per unit injected dose was 0.0001 h/mL and ranged from −0.0011 to 0.0015 h/mL. The median difference in marrow dose from blood was 0.046 mGy/MBq and ranged from −0.300 to 0.430 mGy/MBq (Fig. 2A). There was no statistically significant difference in red marrow dose before and after rituximab treatment (Table 2).
Impact on 90Y-Ibritumomab Dosimetry
The median radiation dose and range (n = 17) for total body, liver, spleen, kidneys, heart wall, lungs, and tumors are listed in Table 2. The paired difference (after rituximab minus before rituximab, n = 16) in radiation dose for the total body (Fig. 2B), liver (Fig. 2C), kidneys (Fig. 2E), heart wall (Fig. 2F), and lungs (Fig. 2G) are also listed in Table 2. The median and range for paired difference in residence times were calculated for the spleen (Fig. 2D) and tumor (Fig. 2H) (Table 2).
There was no significant difference in radiation dose to the total body, liver, kidneys, and lungs before and after rituximab treatment (Table 2). The median radiation dose to the heart wall increased 16% and was statistically significant.
However, the median residence time decreased 24% for the spleen and 28% for tumors after rituximab treatment. These decreases were statistically significant (Table 2). Dosimetry estimates for the tumor or spleen would decrease by the same percentage as that of residence time if there were no volume changes. However, 5 of 27 tumors (excluding 1 tumor of a HAMA patient) were no longer detectable from γ-camera images at the second dosimetry infusion, and posttreatment CT confirmed these tumor responses. For planar image measurements on spleen long and short axis product, 8 of 16 patients (excluding 1 HAMA patient) had a decrease in spleen size greater than 15% at the second dosimetry infusion. The median decrease in spleen size was 14.5%, and the range was 0%−69%.
DISCUSSION
The first 131I-rituximab study for the effect of prior rituximab treatment on subsequent radioimmunotherapy has been recently reported by Illidge et al. (27). For dosimetric evaluation, they assessed changes in the total body dose and marrow radiation dose by examining changes in the effective half-life of the total body and rituximab concentration level in the serum. The current analysis was the first 90Y-ibritumomab dosimetry report for paired comparison in radiation dose change in mGy/MBq for the total body, liver, kidneys, heart wall, lungs, and marrow from blood pharmacokinetics. Radiation dose directly reflects complete information on activity concentration in the tissue over integrated time.
Uptake and clearance can be significantly different when the same monoclonal antibody (Lym-1, CC49, ChL6) is radiolabeled with 131I and 111In (28–36). The difference between 131I-antibody and 111In-antibody was quite prominent in the blood clearance. 111In or another radiometal-labeled antibody generally has faster initial clearance (α-phase), followed by slower secondary clearance (β-phase). This faster initial clearance in the blood is somewhat related to greater liver uptake of the 111In-antibody than the 131I-antibody. Naturally, the impact of rituximab treatment on radiation dosimetry can be different between the 111In/90Y-rituximab (111In/90Y-ibritumomab) reported here and the 131In-rituximab reported by Illidge et al. (27). While Illidge et al. (27) noticed significant increases in radiation dose to total body and marrow based on the increased effective half-life for 131In-rituximab, there was no statistically significant difference in radiation dose to total body and marrow from blood for 111In/90Y-ibritumomab in the present study. We also observed no statistically significant difference in radiation dose to liver, kidneys, and lungs. However, we did observe a statistically significant decrease in residence time for spleen and tumors, as expected. These decreases were in agreement with the observation made by Illidge et al. (27) for 131In-rituximab. We noticed that an increase in radiation dose to the heart wall was statistically significant (P = 0.002), and the median radiation dose increase was 16%. The clinical impact of this increase in the heart wall is unknown, because the heart wall is generally not a dose-limiting factor for 111In/90Y-ibritumomab.
To further improve dosimetry accuracy, we used patient-specific organ mass for dose calculation. For example, patient 16 had a relatively higher liver dose (Fig. 2C) and a small liver mass of 719 g (Table 1). The error could be substantial if a reference liver mass of 1,910 g was used. The median organ masses for the liver and heart were close to the reference mass in OLINDA-EXM. The median mass for lungs and kidneys was approximately 25% higher in this group of patients. This group of patients had moderate splenomegaly; the largest spleen mass was 697 g. The organ radiation doses reported here were generally in the same range as previously reported by Cremonesi et al. (37) and Fisher et al. (10) using patient-specific organ masses.
This retrospective dosimetry analysis was part of a clinical trial to assess the response of the combination of full-course rituximab plus 90Y-ibritumomab. The posttreatment CT images were 13 wk after the first dose of 111In-ibritumomab; the tumor or spleen volume at the second 111In-ibritumomab imaging study was unknown. Five of 27 tumors were no longer detectable from γ-camera images at the second dosimetry infusion, and posttreatment CT confirmed these tumor responses. One logical interpretation of this finding is the treatment effect of the rituximab; another possibility could be complete blocking of all CD20 sites in these tumors by rituximab, preventing 111In- ibritumomab binding at the second dosimetry infusion. The median decrease in spleen size (long and short axis product) was 14.5%, and the range was 0%−69% between the two 111In-ibritumomab image sets. While the residence time or cumulated activity per dose injection accurately reflects the change of spleen and tumor in response to rituximab treatment, percentage changes in tumor and spleen dosimetry would be identical to the percentage changes in residence time if there were no volume changes.
CONCLUSION
Changes in 90Y-ibritumomab dosimetry were not significant for most organs after rituximab treatment, except for a 16% increase to the heart wall. The reduction of residence time for spleen and tumor is likely due to therapeutic effects of rituximab on malignant cells in the spleen and tumor.
Acknowledgments
We thank Gayle Hines for her excellent work in patient imaging and sample counting. The clinical trial was supported by BIOGEN/IDEC Inc.
Footnotes
-
COPYRIGHT © 2010 by the Society of Nuclear Medicine, Inc.
References
- Received for publication May 26, 2009.
- Accepted for publication October 7, 2009.