Considerations of anthropometric, tissue volume, and tissue mass scaling for improved patient specificity of skeletal S values

Improving the dose-myelotoxicity correlation in radiometabolic therapy of bone metastases with 153Sm-EDTMP

PURPOSE

(153)Sm-ethylene diamine tetramethylene phosphonic acid ((153)Sm-EDTMP) is widely used to palliate pain from bone metastases, and is being studied for combination therapy beyond palliation. Conceptually, red marrow (RM) dosimetry allows myelotoxicity to be predicted, but the correlation is poor due to dosimetric uncertainty, individual sensitivity and biological effects from previous treatments. According to EANM guidelines, basic dosimetric procedures have been studied to improve the correlation between dosimetry and myelotoxicity in (153)Sm-EDTMP therapy.

METHODS

RM dosimetry for 33 treatments of bone metastases from breast, prostate and lung tumours was performed prospectively (with (99m)Tc-MDP) and retrospectively, acquiring whole-body scans early and late after injection. The (153)Sm-EDTMP activity was calculated by prospective dosimetry based on measured skeletal uptake and full physical retention, with the RM absorbed dose not exceeding 3.8 Gy. Patient-specific RM mass was evaluated by scaling in terms of body weight (BW), lean body mass (LBM) and trabecular volume (TV) estimated from CT scans of the L2–L4 vertebrae. Correlations with toxicity were determined in a selected subgroup of 27 patients, in which a better correlation between dosimetry and myelotoxicity was expected.

RESULTS

Skeletal uptakes of (99m)Tc and (153)Sm (Tc% and Sm%) were well correlated. The median Sm% was higher in prostate cancer (75.3 %) than in lung (60.5%, p = 0.005) or breast (60.8%, p = 0.008). PLT and WBC nadirs were not correlated with administered activity, but were weakly correlated with uncorrected RM absorbed doses, and the correlation improved after rescaling in terms of BW, LBM and TV. Most patients showed transient toxicity (grade 1–3), which completely and spontaneously recovered over a few days. Using TV, RM absorbed dose was in the range 2–5 Gy, with a median of 312 cGy for PLT in patients with toxicity and 247 cGy in those with no toxicity (p = 0.019), and 312 cGy for WBC in those with toxicity and 232 cGy in those with no toxicity (p = 0.019). ROC curves confirmed the correlations, yielding toxicity absorbed dose thresholds of 265 cGy for PLT and 232 cGy for WBC.

CONCLUSION

The best predictor of myelotoxicity and blood cells nadir was obtained scaling the RM absorbed dose in terms of the estimated TV. It seems clear that the increase in skeletal uptake due to the presence of bone metastases and the assumption of full physical retention cause an overestimation of the RM absorbed dose. Nevertheless, an improvement of the dose–toxicity correlation is easily achievable by simple methods, also leading to possible improvement in multifactorial analyses of myelotoxicity.

A method to correct for radioactivity in large vessels that overlap the spine in imaging-based marrow dosimetry of lumbar vertebrae

Accurate marrow dosimetry for radionuclide therapy based on imaging methods has been challenging because of a variety of factors. One of the uncertainties in image quantification of lumbar vertebrae is correction for radioactivity in large blood vessels anterior to the vertebrae. We developed a method to correct for background radioactivity contributed from blood in large vessels and tested it in a pilot study.

METHODS

CT images of 26 patients receiving (111)In- or (131)I-labeled conjugates were used to measure the inside diameters of the aorta and inferior vena cava (IVC) at the top of L2 and the bottom of L4 and to measure the length of this vessel segment. The volume was calculated for this vessel segment, and then the radioactivity in that volume at each imaging time was determined using a time-variant blood radioactivity concentration as established by serial blood samples. This vessel segment typically overlapped with lumbar vertebrae in anterior and posterior whole-body images. The contribution of this background radioactivity to the cumulated activity of the lumbar spine region of interest (ROI) from serial gamma-camera images was determined, taking into account differences in attenuation between vessel segments and lumbar vertebrae.

RESULTS

The total blood volumes varied from 25 to 94 mL, with a mean of 51 mL. This mean is 76% of the mean marrow volume of 3 lumbar vertebrae measured in some of these patients. Thirteen of the 14 patients evaluated for aortic position had the aortic segment completely within the L2-L4 ROI. For the IVC, a mean of 72% was in the L2-L4 ROI. Adjusting for radioactivity in major blood vessels that were in the ROI led to lower marrow dose estimates.

CONCLUSION

To improve the accuracy of lumbar spine imaging-based marrow dosimetry, one can adjust radioactivity in the large vessels by methods that measure the volume, position, and depth of vessels in the ROI.

CT volumetry of the skeletal tissues

Computed tomography (CT) is an important and widely used modality in the diagnosis and treatment of various cancers. In the field of molecular radiotherapy, the use of spongiosa volume (combined tissues of the bone marrow and bone trabeculae) has been suggested as a means to improve the patient-specificity of bone marrow dose estimates. The noninvasive estimation of an organ volume comes with some degree of error or variation from the true organ volume. The present study explores the ability to obtain estimates of spongiosa volume or its surrogate via manual image segmentation. The variation among different segmentation raters was explored and found not to be statistically significant (p value >0.05). Accuracy was assessed by having several raters manually segment a polyvinyl chloride (PVC) pipe with known volumes. Segmentation of the outer region of the PVC pipe resulted in mean percent errors as great as 15% while segmentation of the pipe's inner region resulted in mean percent errors within approximately 5%. Differences between volumes estimated with the high-resolution CT data set (typical of ex vivo skeletal scans) and the low-resolution CT data set (typical of in vivo skeletal scans) were also explored using both patient CT images and a PVC pipe phantom. While a statistically significant difference (p value <0.002) between the high-resolution and low-resolution data sets was observed with excised femoral heads obtained following total hip arthroplasty, the mean difference between high-resolution and low-resolution data sets was found to be only 1.24 and 2.18 cm3 for spongiosa and cortical bone, respectively. With respect to differences observed with the PVC pipe, the variation between the high-resolution and low-resolution mean percent errors was a high as approximately 20% for the outer region volume estimates and only as high as approximately 6% for the inner region volume estimates. The findings from this study suggest that manual segmentation is a reasonably accurate and reliable means for the in vivo estimation of spongiosa volume. This work also provides a foundation for future studies where spongiosa volumes are estimated by various raters in more comprehensive CT data

Image segmentation of trabecular spongiosa by visual inspection of the gradient magnitude

Recent advances in physical models of skeletal dosimetry utilize high-resolution 3-dimensional microscopic computed tomography images of trabecular spongiosa. These images are coupled to radiation transport codes to assess energy deposition within active bone marrow and trabecular endosteum. These transport codes rely primarily on the segmentation of the spongiosa images into bone and marrow voxels. Image thresholding has been the segmentation of choice for bone sample images because of its extreme simplicity. However, the ability of the segmentation to reproduce the physical boundary between bone and marrow depends on the selection of the threshold value. Statistical models, as well as visual inspection of the image, have been employed extensively to determine the correct threshold. Both techniques are affected by partial volume effect and can provide unexpected results if performed without care. In this study, we propose a new technique to threshold trabecular spongiosa images based on visual inspection of the image gradient magnitude. We first show that the gradient magnitude of the image reaches a maximum along a surface that remains almost independent of partial volume effect and that is a good representation of the physical boundary between bone and marrow. A computer program was then developed to allow a user to compare the position of the iso-surface produced by a threshold with the gradient magnitude. The threshold that produces the iso-surface that best coincides with the maximum gradient is chosen. The technique was finally tested with a set of images of a true bone sample with different resolutions, as well as with three images of a cube of Duocell aluminium foam of known mass and density. Both tests demonstrate the ability of the gradient magnitude technique to retrieve sample volumes or media volume fractions with 1% accuracy at 30 microm voxel size.

Correlations of Total Pelvic Spongiosa Volume With Both Anthropometric Parameters and Computed Tomography–Based Skeletal Size Measurements

Patient-specific dosimetry within the field of molecular radiotherapy continues to pose a challenge owing to the difficulty in predicting marrow toxicity. This study examined the correlation between total pelvic spongiosa volume (TPSV) and independent variables, which include both readily measured or calculated anthropometric parameters (AP), and image-based skeletal measurements requiring computed tomography (CT) images or skeletal radiographs. Fourteen (14) patients (5 male and 9 female) undergoing total hip arthroplasty (THA) were subjected to modified pelvic CT scans. These scans were utilized to estimate TPSV, which was comprised of the volumes of spongiosa within the L5 vertebra, os coxae, sacrum, and both proximal femurs. The APs investigated included total body height (TBH), total body mass (TBM), body mass index (BMI), body surface area (BSA), maximum effective mass (MEM), lean body mass (LBM), and fat-free mass (FFM). Skeletal measurements were also obtained from the CT images of the pelvic region. Correlation coefficients (r) were obtained for TPSV and each set of APs as well as each set of skeletal measurements. Total body height (r = 0. 80) and os coxae height (r = 0.83) had the highest correlation coefficients of all the APS and skeletal measurements, respectively. FFM (r = 0.50), LBM (r = 0.42), TBM (r = 0.11), and BSA (r = 0.11) did not correlate well with TPSV, which accounts for approximately 45% of total spongiosa seen throughout the skeleton at sites associated with active bone marrow. Skeletal height measurements appear to have a much higher correlation with TPSV than either their corresponding skeletal width measurements or parameters that are a function of an individual's TBM.

Accounting for beta‐particle energy loss to cortical bone via paired‐image radiation transport (PIRT)

Current methods of skeletal dose assessment in both medical physics (radionuclide therapy) and health physics (dose reconstruction and risk assessment) rely heavily on a single set of bone and marrow cavity chord-length distributions in which particle energy deposition is tracked within an infinite extent of trabecular spongiosa, with no allowance for particle escape to cortical bone. In the present study, we introduce a paired-image radiation transport (PIRT) model which provides a more realistic three-dimensional (3D) geometry for particle transport in the skeletal site at both microscopic and macroscopic levels of its histology. Ex vivo CT scans were acquired of the pelvis, cranial cap, and individual ribs excised from a 66-year male cadaver (BMI of 22.7 kg m(-2)). For the three skeletal sites, regions of trabecular spongiosa and cortical bone were identified and segmented. Physical sections of interior spongiosa were taken and subjected to microCT imaging. Voxels within the resulting microCT images were then segmented and labeled as regions of bone trabeculae, endosteum, active marrow, and inactive marrow through application of image processing algorithms. The PIRT methodology was then implemented within the EGSNRC radiation transport code whereby electrons of various initial energies are simultaneously tracked within both the ex vivo CT macroimage and the CT microimage of the skeletal site. At initial electron energies greater than 50-200 keV, a divergence in absorbed fractions to active marrow are noted between PIRT model simulations and those estimated under existing techniques of infinite spongiosa transport. Calculations of radionuclide S values under both methodologies imply that current chord-based models may overestimate the absorbed dose to active bone marrow in these skeletal sites by 0% to 27% for low-energy beta emitters (33P, 169Er, and 177Lu), by approximately 4% to 49% for intermediate-energy beta emitters (153Sm, 186Re, and 89Sr), and by approximately 14% to 76% for high-energy beta emitters (32p, 188Re, and 90Y). The PIRT methodology allows for detailed modeling of the 3D macrostructure of individual marrow-containing bones within the skeleton thus permitting improved estimates of absorbed fractions and radionuclide S values for intermediate-to-high energy beta emitters.

Backscatter and Dose Perturbations for Low- to Medium-Energy Electron Point Sources at the Interface between Materials with Different Atomic Numbers

Electron backscatter at interfaces between dissimilar media can affect dosimetry and should be taken into consideration in radiotherapy and in radiobiology experiments. Backscatter dose perturbations depend upon factors such as electron energy, medium atomic number (Z), and distance from the interface. This study quantifies the backscatter dose factor (BSDF) for electron point sources of energy between 0.1 to 3 MeV in water at the interface with scattering materials ranging in Z from (13)Al to (79)Au. A Monte Carlo code that performs dose calculations for monoenergetic and continuous-spectrum electron sources was developed using EGSnrc transport routines. The BSDF was quantified in a parallel layers geometry (BSDF(1D)) and three-dimensional voxel geometry (BSDF(3D)). The BSDF(1D) near the interface increased up to 52% with decreasing energy from 3 to 0.1 MeV and increasing Z from 13 to 79. The analysis of the BSDF(3D) showed a significant dependence of the scattered electron angular distribution on Z and energy, with a decrease in isotropy going from high to low Z. This effect proves the importance of considering the correct geometry when quantifying the BSDF for electron sources, especially when the dimensions of the relevant dose-collecting volume are comparable with the CSDA range of the source.

A hyperboliod representation of the bone-marrow interface within 3D NMR images of trabecular bone: applications to skeletal dosimetry

Recent advances in physical models of skeletal dosimetry utilize high-resolution NMR microscopy images of trabecular bone. These images are coupled to radiation transport codes to assess energy deposition within active bone marrow irradiated by bone- or marrow-incorporated radionuclides. Recent studies have demonstrated that the rectangular shape of image voxels is responsible for cross-region (bone-to-marrow) absorbed fraction errors of up to 50% for very low-energy electrons (<50 keV). In this study, a new hyperboloid adaptation of the marching cube (MC) image-visualization algorithm is implemented within 3D digital images of trabecular bone to better define the bone-marrow interface, and thus reduce voxel effects in the assessment of cross-region absorbed fractions. To test the method, a mathematical sample of trabecular bone was constructed, composed of a random distribution of spherical marrow cavities, and subsequently coupled to the EGSnrc radiation code to generate reference values for the energy deposition in marrow or bone. Next, digital images of the bone model were constructed over a range of simulated image resolutions, and coupled to EGSnrc using the hyperboloid MC (HMC) algorithm. For the radionuclides 33P, 117mSn, 131I and 153Sm, values of S(marrow<--bone) estimated using voxel models of trabecular bone were shown to have relative errors of 10%, 9%, <1% and <1% at a voxel size of 150 microm. At a voxel size of 60 microm, these errors were 6%, 5%, <1% and <1%, respectively. When the HMC model was applied during particle transport, the relative errors on S(marrow<--bone) for these same radionuclides were reduced to 7%, 6%, <1% and <1% at a voxel size of 150 microm, and to 2%, 2%, <1% and <1% at a voxel size of 60 microm. The technique was also applied to a real NMR image of human trabecular bone with a similar demonstration of reductions in dosimetry errors.