Abstract
The data that have been used in almost all calculations of MIRD S value absorbed dose and effective dose are based on stylized anatomic computational phantoms and tissue-weighting factors adopted by the International Commission on Radiological Protection (ICRP) in its publication 60. The more anatomically realistic phantoms that have recently become available are likely to provide more accurate effective doses for diagnostic agents. 68Ga-DOTATATE is a radiolabeled somatostatin analog that binds with high affinity to somatostatin receptors, which are overexpressed in neuroendocrine tumors and can be used for diagnostic PET/CT-based imaging. Several studies have reported effective doses for 68Ga-DOTATATE using the stylized Cristy–Eckerman (CE) phantoms from 1987; here, we present effective dose calculations using both the ICRP 60 and more updated formalisms. Methods: Whole-body PET/CT scans were acquired for 16 patients after 68Ga-DOTATATE administration. Contours were drawn on the CT images for spleen, liver, kidneys, adrenal glands, brain, heart, lungs, thyroid gland, salivary glands, testes, red marrow (L1–L5), muscle (right thigh), and whole body. Dosimetric calculations were based on the CE phantoms and the more recent ICRP 110 reference-voxel phantoms. Tissue-weighting factors from ICRP 60 and ICRP 103 were used in effective dose calculations for the CE phantoms and ICRP 110 phantoms, respectively. Results: The highest absorbed dose coefficients (absorbed dose per unit activity) were, in descending order, in the spleen, pituitary gland, kidneys, adrenal glands, and liver. For ICRP 110 phantoms with tissue-weighting factors from ICRP 103, the effective dose coefficient was 0.023 ± 0.003 mSv/MBq, which was significantly lower than the 0.027 ± 0.005 mSv/MBq calculated for CE phantoms with tissue-weighting factors from ICRP 60. One of the largest differences in estimated absorbed dose coefficients was for the urinary bladder wall, at 0.040 ± 0.011 mGy/MBq for ICRP 110 phantoms compared with 0.090 ± 0.032 mGy/MBq for CE phantoms. Conclusion: This study showed that the effective dose coefficient was slightly overestimated for CE phantoms, compared with ICRP 110 phantoms using the latest tissue-weighting factors from ICRP 103. The more detailed handling of electron transport in the latest phantom calculations gives significant differences in estimates of the absorbed dose to stem cells in the walled organs of the alimentary tract.
Dosimetry for nuclear medicine imaging agents is performed to ensure that the long-term radiation risks of the imaging procedure are minimal compared with the benefits of obtaining a correct image-based diagnosis. In the late 1960s and early 1970s, the MIRD Committee of the Society of Nuclear Medicine and Molecular Imaging published a formalism that standardized the process for computing tissue-absorbed doses through the introduction of 2 essential terms, cumulated activity (Ã) and radionuclide S value (S) (1). Cumulated activity, currently known by the term time-integrated activity (abbreviated TIA but still represented by Ã) (2), depends on the pharmacokinetic properties of the diagnostic agent. These are typically characterized by longitudinal quantitative imaging in an appropriate patient population. S value is defined as the absorbed dose to a target region per unit TIA in the source region, such that the product of TIA and S value gives the absorbed dose contribution from a source region to a target region. S values depend on the emission characteristics of the radionuclide used in the diagnostic agent and on the reference anatomic model used in the Monte Carlo radiation transport simulations of the emitted radiation particles and photons. The S values that have been used in effective dose calculations over the past several decades have been derived primarily from the stylized Cristy–Eckerman (CE) phantom series developed in 1987 (3) and used extensively by the International Commission on Radiological Protection (ICRP). Reflecting the computational and imaging capabilities of the time, the CE phantoms comprised stylized organs composed of simple geometries defined by mathematic surface equations to describe the inner anatomy and outer body contour. In 2009, the ICRP and the International Commission on Radiation Units and Measurements released adult male and female CT-based voxelized reference phantoms (ICRP publication 110) (4), with organ and tissue masses matched to reference values (ICRP publication 89) (5). Specific absorbed fractions for the ICRP 110 voxelized reference phantoms have been published (ICRP publication 133) (6), as well as the latest reference tissue–weighting factors (ICRP publication 103) (7).
In the present study, we used 68Ga-DOTATATE PET/CT scans to compare dosimetry based on the most recent ICRP standards (ICRP 110 phantoms and ICRP 103 tissue-weighting factors) with dosimetry based on the prior standards (CE phantoms and ICRP 60 tissue-weighting factors (8)) and with previously reported dosimetry based on these prior standards (9,10). Such comparisons are important to understanding the impact of transitioning from the older phantoms and tissue-weighting factors to the most recently recommended ones.
MATERIALS AND METHODS
Patients
68Ga-DOTATATE PET/CT scans from 16 Brazilian patients with somatostatin-avid tumors were used in this study (11 female and 5 male; mean age, 53.8 ± 13.0 y; range, 36–79 y) (Table 1). PET was indicated for staging, follow-up, or planning of peptide receptor radionuclide therapy. Five of the patients had gastroenteropancreatic neuroendocrine tumors, 5 had medullary thyroid cancer with increasing calcitonin levels, 3 had neuroectodermic tumors (malignant pheochromocytoma or paraganglioma), 2 had associated adrenal and thyroid tumors, and 1 had ectopic Cushing syndrome requiring tumor localization. The institutional review board at São Paulo University School of Medicine approved this study, and all subjects signed an informed-consent form.
PET/CT Imaging
Two to four whole-body PET/CT scans were acquired at 2–240 min after injection of 68Ga-DOTATATE (mean, 131.2 ± 26.3 MBq) (Table 1). Imaging was performed in 3-dimensional time-of-flight mode on a Discovery PET/CT 690 system (GE Healthcare) at 2 min per bed position. The acquisition matrix was 192 × 192 pixels, the in-slice pixel size was 3.27 mm, and the transaxial slice thickness was 3.27 mm. Ordered-subsets expectation maximization was used for PET reconstruction, with CT-based attenuation correction using automatic exposure control for dose reduction.
Normal-Tissue Dosimetry
Overview
Patient-tissue activity concentrations rather than organ-specific activities were taken from the patient data and transposed into the 2 reference phantom models. Normal-tissue absorbed dose coefficients (absorbed dose per unit activity) were calculated using the MIRD S-value–based methodology (1,2). OLINDA/EXM software, version 1 (11), was used for the CE phantom dosimetry calculations, and a Microsoft Excel spreadsheet with S values obtained from the University of Florida was used for the ICRP 110 phantom dosimetry calculations. TIA (i.e., total number of nuclear transformations between 2 time points) was calculated from time–activity data within the source regions, defined by drawing organ contours on the CT portion of the first PET/CT image for each patient. These contours were transposed to later time points after deformable registration of the scans had been performed using Velocity software (version 3.1; Varian Medical Systems Inc.). Contours were drawn for the spleen, liver, kidneys, adrenal glands, brain, heart, lungs, thyroid gland, salivary glands, testes, red marrow (L1–L5), muscle (right thigh), and whole body. Dosimetry was performed for each patient and inserted into the respective phantom, after which the average for all patients was calculated.
TIA Calculations
The average activity concentration is , where A is activity, V is volume, rS is source region, and t is time after injection. Volume was derived from each patient’s PET/CT images. The specific organ/tissue activity (Bq) for a phantom κ was calculated for each patient’s organ/tissue activity as follows:Eq. 1
where M is mass and ρ is tissue density. Mass was taken from OLINDA/EXM (11,12) for the CE phantom calculations and from ICRP 133 for the ICRP 110 phantom calculations. Tissue density was according to ICRP 110. Source region rS in Equation 1 is an index over all tissues and the whole body for which TIA was assigned.
Excel was used to fit a monoexponential expression to each patient’s specific organ/tissue time–activity curve. The expressions were analytically integrated from 0 to infinity, resulting in the TIA corresponding to each phantom’s source region. Preserving activity concentration rather than total activity can lead to a discrepancy in the whole-body activity between the phantom and the patient at the time of administration. The TIA coefficient (TIAC) is defined as the TIA divided by the administered activity. In this scenario, the administered activity is the total virtually administered activity to the phantom rather than the administered activity to the patient. Scaling of each TIAC by the phantom-to-patient whole-body mass ratio accounts for this difference. The TIAs (Bq⋅s) used for the calculations were as follows:Eq. 2The TIACs (s) whereEq. 3where ã is TIAC, wb is whole body, and A0 is administered activity.
Remainder of Body
The remainder-of-body TIAC (s) was calculated as the difference between the TIAC of the whole body, the TIAC of the target region, and the TIAC of the source regions summed:Eq. 4where rb is rest of body, wb is whole body, and rT is target region. The corresponding phantom masses (kg) wereEq. 5Source region rS in Equations 4 and 5 is an index over all source regions for which a TIAC is assigned, except for the target region. The remainder-of-body mass excluded the contents of the walled organs.
Urinary Bladder
The TIAC for urinary bladder contents was calculated using the MIRD bladder-voiding model (13) as implemented in OLINDA/EXM. A 2-h voiding interval was used, with the whole-body biologic clearance half-life obtained from the PET/CT images and a fraction of 1. This TIAC was used both for the CE phantom calculations and for the ICRP 110 phantom calculations.
Pituitary Gland
To account for partial-volume effects arising from the small size of the pituitary gland, a hybrid dosimetric method developed by Plyku et al. (14) for small tumors was used. Briefly, the pituitary gland was modeled as a 0.6-g unit-density sphere. A 10% threshold of the maximum PET value was used for contouring. The volume of the pituitary gland contoured on the PET images was systematically larger than the ICRP volume and presumably larger than the real volume because of activity spill-out. However, this volume contains the activity in the pituitary gland plus some background activity (Bq) in the PET volume outside the real volume, which needs to be subtracted:Eq. 6where pg is pituitary gland and bkg is background. The background activity concentration for the pituitary gland was determined by drawing a contour in the brain. A monoexponential expression was fitted to the time–activity data points and analytically integrated from zero to infinity, and the TIAC was calculated according to Equations 2 and 3 for the respective phantom.
Absorbed Dose Calculations
The MIRD Committee–derived equation was used for all calculations of absorbed dose coefficients (absorbed dose per unit activity (Gy/Bq)) (2):Eq. 7where d is absorbed dose coefficient. The TIACs for the CE phantoms were used as input in OLINDA/EXM to obtain estimates of the CE phantom–derived absorbed dose coefficient. Because the CE phantoms do not have S values for the pituitary gland, the self-dose absorbed dose coefficient to the pituitary gland was calculated using the unit-density sphere model provided by OLINDA/EXM. An Excel spreadsheet was used for the ICRP dosimetric calculations.
Remainder-of-Body Contribution to Absorbed Dose
The CE phantom–based dosimetry specifies 25 target regions and 28 source regions. In contrast, the ICRP 110 phantom dosimetry allows TIAC apportionment in up to 76 source regions and provides the absorbed dose coefficients to 41 target regions. More significantly, the CE phantom–based formalism includes whole-body–to–individual-tissue S values, whereas no such S values are tabulated for the ICRP 110 phantoms. The remainder TIAC that has not been otherwise allocated is given the term remainder of body in CE phantom calculations. Target absorbed dose coefficients (Gy/Bq) for remainder of body are given by (15)Eq. 8with the remainder-of-body to target S value (Gy/Bq⋅s) given byEq. 9Because almost all tissues in the body are accounted for by the ICRP 110 phantoms, whole-body–to–individual-tissue S values have not been calculated for the ICRP 110 phantoms. The contribution of the absorbed dose coefficient from remainder of body to each target tissue is obtained by apportioning the remainder-of-body TIAC (s) to each remaining source tissue not previously accounted for:Eq. 10where is the index over source regions identified in the ICRP 110 phantoms that have not been specifically assigned a TIAC.
The absorbed dose coefficient (Gy/Bq) to each target region of the ICRP 110 phantoms is given byEq. 11
Effective Dose Calculations
The effective dose coefficient (effective dose per unit activity (Sv/Bq)) for CE phantom calculations in an individual patient is given byEq. 12where e is effective dose coefficient, wT is ICRP 60 tissue-weighting factor, and h is equivalent dose coefficient (Sv/Bq).
The radiation-weighting factor for all 68Ga radiation emissions is equal to 1. Accordingly,Eq. 13where α denotes the sex of each patient.
The effective dose coefficient per unit activity (Sv/Bq) is obtained by averaging the male- and female-averaged patient-specific effective dose coefficient:Eq. 14where and are the number of males and females, respectively.
The effective dose coefficient (Sv/Bq) for ICRP isEq. 15where the tissue-weighting factor is from ICRP 103.
The percentage differences between the CE phantom values and the ICRP 110 phantom values were calculated as follows:Eq. 16
Statistical Analysis
Statistical analysis was performed using Prism software (version 7.03; GraphPad Software Inc.). All data are presented as the mean value ± SD. Groups were compared using a 2-tailed Student t test, and differences between groups were considered significant for P values of less than 0.05.
RESULTS
The absorbed dose coefficients calculated for the ICRP 110 phantoms are shown in Figure 1. The tissues with the highest 7 absorbed dose coefficients, as well as the effective dose coefficient, are listed and compared with previously published data (9,10) in Table 2. ICRP 110 phantom–based calculations for the 16 patients are shown in Figure 2. TIACs are listed in Table 3. The absorbed dose coefficients obtained in this study are compared with those in the literature (9) in Figure 3. The ICRP 110 phantom–derived effective dose coefficient (0.023 ± 0.003 mSv/MBq) obtained using tissue-weighting factors from ICRP 103 was significantly lower (P = 0.0114) than that (0.027 ± 0.005 mSv/MBq) derived from the CE phantoms using tissue-weighting factors from ICRP 60. Significant differences in calculated absorbed dose coefficients between the CE phantoms and the ICRP 110 phantoms were found for gallbladder wall, breast, thymus, pancreas, small intestine wall, stomach wall, uterus, ovaries, skin, heart wall, urinary bladder wall (P < 0.0001), and pituitary gland (P = 0.0252). The contributions from remainder of body, source tissues, and target tissue to the total absorbed dose coefficients for kidney and urinary bladder wall are shown in Figure 4.
DISCUSSION
The main objective of radiopharmaceutical dosimetry for diagnostic agents is to ensure that the risks of the imaging procedure are minimal compared with the benefits. Accordingly, dosimetry for diagnostic agents is not performed for any individual patient. Rather, agent-specific pharmacokinetic data from an appropriate patient population are collected and used in conjunction with an internationally recognized standard reference geometry that provides the S values required to perform the calculation. Standardization of all elements required to arrive at an absorbed dose estimate is important in maintaining consistency across centers and ensuring that the values reported in publications reflect actual differences arising from the properties of the agent rather than from different data acquisition protocols and corresponding S values. The release of a new set of reference geometries and associated tissue-weighting factors requires that physicists and nuclear medicine practitioners transition from the older, CE phantom, set of standards to the newer, ICRP 110 phantom, set. Using data from 68Ga-DOTATATE, we have investigated the effect on absorbed dose calculations arising from the switch in phantoms and corresponding S values.
In the ICRP 110 phantom–based calculations, the remainder activity was uniformly distributed in the remaining tissues using the respective organ/tissue weights, including blood, according to ICRP 133. The remainder activity in walled organs was uniformly distributed in the wall according to Equation 10 and not in the contents (urinary bladder contents were the exception because urinary bladder was a source region). This explanation accounts to some extent for the higher absorbed dose coefficients in these walled organs for the ICRP 110 phantoms than for the CE phantoms using the contents as the source. There are also differences in tissues definitions between the ICRP 110 phantoms and the CE phantoms; for example, in the ICRP 110 phantoms, the right colon wall includes the ascending colon wall and the right half of the transverse colon wall, and the left colon wall includes the descending colon wall and the left half of the transverse colon wall. In the CE phantoms, the lower large intestine corresponds to the descending colon, sigmoid colon, and rectum, and the upper large intestine corresponds to the ascending colon and transverse colon (16). In addition, the lungs are divided into 4 target regions in the ICRP 110 phantoms, compared with 1 target region in the CE phantoms.
When the calculated self-dose S values were compared between the CE and ICRP 110 phantoms, the CE phantom values were higher for both females (range, −7.7% to +65.3%) and males (range, −961% to +59.6%), with the largest difference being for the urinary bladder contents to the urinary bladder wall (+65.3% in females and +59.6% in males) and an outlier for the self-dose breast S value (−961%) in males. The standard mass for the male breast was 351 g (11,12) and 26.2 g (6) for the CE and ICRP 110 phantoms, respectively, and is probably the main explanation for the large difference observed in the respective S values. A comparison of the CE and the ICRP 110 phantoms using the averaged male and female self-dose S values for all organs and tissues gave a difference of +26.3% and +27.1%, respectively, with a respective tissue mass difference of −2.2% and −6.7%. However, a comparison using the specific TIACs and organ/tissue masses from the ICRP 110 method in the CE method and comparing with the CE phantom results gave an averaged male and female absorbed dose difference of 2.1% for the urinary bladder wall and 1.8% for the kidneys; thus, the effect of mass differences on the results is small because of the scaling of input activity to the phantom organ mass.
In general, the results showed lower absorbed doses using ICRP 110 phantoms than using CE phantoms. This finding is likely explained by 2 fundamental differences between these phantoms. The first is a difference in organ topology, with the ICRP 110 phantoms having less interorgan tissue space and a better accounting and placement of distributed tissue (e.g., marrow, adipose tissue, and muscle). The second difference is that S values derived from the ICRP 110 phantoms have a better handling of the fraction of electron energy absorbed. This characteristic explains to some extent the larger differences seen in thin-walled organs (e.g., urinary bladder wall, which had a 125% higher absorbed dose coefficient in the CE phantom calculations than in the ICRP 110 phantom calculations). This advantage leads to a different and better approach to handling TIACs not specifically assigned to source organs. It is now possible to deposit TIACs uniformly throughout the remaining non-source tissue and then use specific source-to-target S values to account for the absorbed dose due to a TIAC that is not specifically allocated. Other differences in the calculations (e.g., updated radionuclide decay data (17) and tissue densities) had negligible effects on the differences, as previously shown by Hadid et al. (18).
Dosimetric calculations for 68Ga-DOTATATE have previously been published by Walker et al. (9) and Sandström et al. (10). The percentage difference between results obtained in this study using the CE phantoms and those of Walker et al. ranged from −229% to +86.5% (brain and adrenal glands, respectively). The absorbed dose coefficient to the adrenal glands reported by Sandström et al. was 21% lower than that obtained in this study. The large adrenal-gland difference for Walker et al. could be explained by their not using the adrenal glands as a source organ (no contour was drawn), whereas our study and that of Sandström et al. did use drawn contours. In addition, the absorbed dose coefficient to the salivary glands calculated by Walker et al. (0.012 ± 0.008 mGy/MBq) was 52% lower than that obtained using the ICRP 110 phantom calculation (0.025 ± 0.010 mGy/MBq). Similarly, the salivary glands were not included in the CE phantoms as a source or target region. In addition, Walker et al. used data exclusively from male patients, whereas our study used the average between male and female. For both the salivary and the pituitary glands, there were significant differences (P < 0.035) between the male and female absorbed dose coefficients (Fig. 1). In addition, because the ICRP 110 phantoms use the salivary glands both as a source and as a target region, the absorbed dose considers not only the self-dose but also contributions from adjacent and surrounding tissues (unlike version 1 of the OLINDA/EXM sphere model, for example). Furthermore, Walker et al. did not perform whole-body PET/CT and made estimates for the activity in the nonimaged extremities to get the whole-body activity. Sandström et al. also did not perform whole-body PET/CT but, rather, scanned from the base of the skull to the proximal femur. Because we performed whole-body (head to toes) scans in this work, we were able to account for all the activity in the body.
Figure 4 depicts the averaged male and female contribution from different source tissues to the kidneys and urinary bladder wall. These 2 tissues were chosen to illustrate the differences in relative contribution to the total absorbed dose between a solid organ and a walled tissue, respectively. As expected, self-dose is the dominant contributor to the total absorbed dose for the kidneys. In contrast, urinary bladder contents provide a greater contribution to the bladder wall dose than does the wall itself, as can be seen in the ICRP 110 bladder wall dose calculation. Urinary bladder wall self-dose S values are available for the ICRP 110 phantom calculations but not for the CE phantom calculations. The absorbed dose coefficient for the kidneys was 5.9% higher for the CE phantoms than for the ICRP 110 phantoms, as can be explained mainly by a 29.7% higher self-dose S value, a 30.7% lower average kidney TIAC, and the difference in contributions from source tissues and remainder of body. For the urinary bladder wall, the absorbed dose coefficient was 125% higher for the CE phantoms than for the ICRP 110 phantoms, as can be explained mainly by a 62.5% higher urinary bladder contents-to-wall S value, no difference in TIAC (the TIACs for urinary bladder contents were calculated using the MIRD bladder-voiding model as implemented in OLINDA/EXM, version 1, and is not phantom-dependent), and the difference in contributions from source tissues and remainder of body (including urinary bladder wall for the ICRP 110 phantoms). This difference stems from a difference in methods, with the CE phantoms using only contents as a source for walled organs (e.g., stomach, gallbladder, intestine, and urinary bladder) but the ICRP 110 phantoms using organ wall as a source as well.
As shown in Figures 2 and 3, the calculated absorbed dose coefficients for the different organs/tissues varied widely among patients, with the differences being in many cases larger than those between the ICRP 110 and CE phantom calculations. About half the compared organs/tissues (mainly those that are small or walled) differed significantly between the ICRP 110 and CE phantom calculations. The calculated absorbed doses for larger solid organs (e.g., liver) did not significantly differ between the ICRP 110 and CE phantoms. The significant differences in the smaller or walled organs represent relatively low absorbed doses compared with larger solid organs such as the spleen. Tissue-absorbed doses are already low for diagnostic agents. The modest reduction in tissue-absorbed doses and in effective dose will have a minimal impact on routine clinical practice. The main impact of these updated calculations will be in the development of new imaging agents. In this context, lower absorbed doses would allow slightly greater administered activities and therefore higher-quality images.
Another area that would be impacted, although negatively so, is the therapeutic use of radionuclides. Even though dosimetry for radiopharmaceutical therapy should account for individual patient anatomy and the spatial distribution of activity within organs, standard phantom-based dosimetry methods are still used for initial assessments of likely toxicity and for reporting to regulatory authorities.
The tissue-weighting factors used with the ICRP 110 phantoms were from ICRP 103, whereas those used with the CE phantoms were from ICRP 60. The main change from the earlier tissue-weighting factors is a decrease for gonads, an increase for breast, the use of separate weighting factors for salivary glands and brain, and an increase for the remainder category, both in value and in number of tissues included (7,8). The calculated effective dose coefficient was significantly lower for the ICRP 110 phantoms than for the CE phantoms. The recommended administered activity of 68Ga-DOTATATE (range, 100–200 MBq) (19)—and the average administered in this study—was 131.2 MBq. According to the ICRP 110 phantoms, this amount of activity results in an effective dose of 3.0 mSv (2.3–4.6 mSv for the recommended activity range), which is lower than the 3.6 mSv calculated using the CE phantoms. Accordingly, effective doses have been overestimated using the CE phantoms, as compared with the ICRP 110 phantoms in combination with specific absorbed fractions from ICRP 133 and tissue-weighting factors from ICRP 103.
CONCLUSION
CE phantom–based dosimetry has been shown to overestimate effective dose. The differences that were found between CE phantom–based dosimetry and ICRP 110–based dosimetry for 68Ga-DOTATATE may require a reevaluation of the dosimetry for other diagnostic agents as well. These new standards should be adopted rapidly so as to minimize confusion in the dosimetry literature.
DISCLOSURE
This work was supported by FAPESP 13/03876-4 (“Avaliação do uso do 68Ga-peptídeo análogo de somatostatina PET/CT como ferramenta diagnóstica em tumores neuroendócrinos e sua correlação com marcadores moleculares”) and by NIH R01 CA116477. No other potential conflict of interest relevant to this article was reported.
Acknowledgments
Drs. George Barberio Coura Filho and Ana Amélia Fialho de Oliveira Hoff, from Instituto do Cancer do Estado de São Paulo, are acknowledged for helping us acquire the PET/CT images. We also thank Dr. William Goodwin, from the Department of Biomedical Engineering at the University of Florida, for valuable information and insight on the ICRP 110 voxelized phantoms.
Footnotes
Published online Feb. 9, 2018.
- © 2018 by the Society of Nuclear Medicine and Molecular Imaging.
REFERENCES
- Received for publication October 18, 2017.
- Accepted for publication December 14, 2017.