Skip to main content

Main menu

  • Home
  • Content
    • Current
    • Ahead of print
    • Past Issues
    • JNM Supplement
    • SNMMI Annual Meeting Abstracts
    • Continuing Education
    • JNM Podcasts
  • Subscriptions
    • Subscribers
    • Institutional and Non-member
    • Rates
    • Journal Claims
    • Corporate & Special Sales
  • Authors
    • Submit to JNM
    • Information for Authors
    • Assignment of Copyright
    • AQARA requirements
  • Info
    • Reviewers
    • Permissions
    • Advertisers
  • About
    • About Us
    • Editorial Board
    • Contact Information
  • More
    • Alerts
    • Feedback
    • Help
    • SNMMI Journals
  • SNMMI
    • JNM
    • JNMT
    • SNMMI Journals
    • SNMMI

User menu

  • Subscribe
  • My alerts
  • Log in
  • Log out
  • My Cart

Search

  • Advanced search
Journal of Nuclear Medicine
  • SNMMI
    • JNM
    • JNMT
    • SNMMI Journals
    • SNMMI
  • Subscribe
  • My alerts
  • Log in
  • Log out
  • My Cart
Journal of Nuclear Medicine

Advanced Search

  • Home
  • Content
    • Current
    • Ahead of print
    • Past Issues
    • JNM Supplement
    • SNMMI Annual Meeting Abstracts
    • Continuing Education
    • JNM Podcasts
  • Subscriptions
    • Subscribers
    • Institutional and Non-member
    • Rates
    • Journal Claims
    • Corporate & Special Sales
  • Authors
    • Submit to JNM
    • Information for Authors
    • Assignment of Copyright
    • AQARA requirements
  • Info
    • Reviewers
    • Permissions
    • Advertisers
  • About
    • About Us
    • Editorial Board
    • Contact Information
  • More
    • Alerts
    • Feedback
    • Help
    • SNMMI Journals
  • View or Listen to JNM Podcast
  • Visit JNM on Facebook
  • Join JNM on LinkedIn
  • Follow JNM on Twitter
  • Subscribe to our RSS feeds
Research ArticleClinical Investigation
Open Access

An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 2: Overall Variabilities in Absorbed Dose

Julia Brosch-Lenz, Suqi Ke, Hao Wang, Eric Frey, Yuni K. Dewaraja, John Sunderland and Carlos Uribe
Journal of Nuclear Medicine July 2023, 64 (7) 1109-1116; DOI: https://doi.org/10.2967/jnumed.122.265094
Julia Brosch-Lenz
1Department of Integrative Oncology, BC Cancer Research Institute, Vancouver, British Columbia, Canada;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Suqi Ke
2Division of Quantitative Sciences, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Hao Wang
2Division of Quantitative Sciences, Sidney Kimmel Comprehensive Cancer Center, Johns Hopkins University School of Medicine, Baltimore, Maryland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric Frey
3Rapid, LLC, Baltimore, Maryland;
4Department of Radiology, Johns Hopkins University, Baltimore, Maryland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yuni K. Dewaraja
5Department of Radiology, University of Michigan, Ann Arbor, Michigan;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
John Sunderland
6Department of Radiology, University of Iowa, Iowa City, Iowa;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Carlos Uribe
1Department of Integrative Oncology, BC Cancer Research Institute, Vancouver, British Columbia, Canada;
7Department of Functional Imaging, BC Cancer, Vancouver, British Columbia, Canada; and
8Department of Radiology, University of British Columbia, Vancouver, British Columbia, Canada
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • PDF
Loading

Visual Abstract

Figure
  • Download figure
  • Open in new tab
  • Download powerpoint

Abstract

Dosimetry for personalized radiopharmaceutical therapy has gained considerable attention. Many methods, tools, and workflows have been developed to estimate absorbed dose (AD). However, standardization is still required to reduce variability of AD estimates across centers. One effort for standardization is the Society of Nuclear Medicine and Molecular Imaging 177Lu Dosimetry Challenge, which comprised 5 tasks (T1–T5) designed to assess dose estimate variability associated with the imaging protocol (T1 vs. T2 vs. T3), segmentation (T1 vs. T4), time integration (T4 vs. T5), and dose calculation (T5) steps of the dosimetry workflow. The aim of this work was to assess the overall variability in AD calculations for the different tasks. Methods: Anonymized datasets consisting of serial planar and quantitative SPECT/CT scans, organ and lesion contours, and time-integrated activity maps of 2 patients treated with 177Lu-DOTATATE were made available globally for participants to perform dosimetry calculations and submit their results in standardized submission spreadsheets. The data were carefully curated for formal mistakes and methodologic errors. General descriptive statistics for ADs were calculated, and statistical analysis was performed to compare the results of different tasks. Variability in ADs was measured using the quartile coefficient of dispersion. Results: ADs to organs estimated from planar imaging protocols (T2) were lower by about 60% than those from pure SPECT/CT (T1), and the differences were statistically significant. Importantly, the average differences in dose estimates when at least 1 SPECT/CT acquisition was available (T1, T3, T4, T5) were within ±10%, and the differences with respect to T1 were not statistically significant for most organs and lesions. When serial SPECT/CT images were used, the quartile coefficients of dispersion of ADs for organs and lesions were on average less than 20% and 26%, respectively, for T1; 20% and 18%, respectively, for T4 (segmentations provided); and 10% and 5%, respectively, for T5 (segmentation and time-integrated activity images provided). Conclusion: Variability in ADs was reduced as segmentation and time-integration data were provided to participants. Our results suggest that SPECT/CT-based imaging protocols generate more consistent and less variable results than planar imaging methods. Effort at standardizing segmentation and fitting should be made, as this may substantially reduce variability in ADs.

  • radiopharmaceutical therapy
  • dosimetry
  • variability
  • standardization
  • dosimetry challenge

Recent clinical trials have demonstrated favorable patient outcomes and led to U.S. Food and Drug Administration approvals of 177Lu-based radiopharmaceuticals for the treatment of neuroendocrine tumors (177Lu-DOTATATE (1) in 2018) and metastatic castration-resistant prostate cancer (177Lu-PSMA-617 (2) in 2022). These approvals have motivated research for new targets and development of new radiopharmaceuticals by both academia and industry (3–7). Despite the initial promising results for radiopharmaceutical therapies, recurrence has also been reported (8,9). Patient-specific dosimetry (10) may allow personalization of administered activity to deliver maximized absorbed doses (ADs) to lesions while keeping normal-organ ADs below toxic levels. There is evidence that dosimetry-guided therapy increases the survival of patients who undergo liver radioembolization (11). However, dosimetry calculations are still not routinely performed for radiopharmaceutical therapies, partly because of the lack of standardized dosimetry tools, methods, and protocols. Dosimetry-based therapy-planning approaches are fundamentally limited by the precision of the AD estimates. However, relatively little is known about variability in ADs and the extent to which the different steps of the dosimetry workflow contribute to it.

Dosimetry calculations involve multiple steps (12–14): quantitative imaging of the distribution of the radiopharmaceutical over time, segmentation of lesions and organs of interest, estimation of the total number of disintegrations (time-integrated activity (TIA)) in each target region (e.g., organs and lesions), and conversion of TIA to AD using either organ-level or voxelized dosimetry methods. Alternatively, serial dose-rate images can be calculated first, followed by fitting and integration over time.

To better understand the relative contribution to variability of the various steps of the dosimetry workflow, the Society of Nuclear Medicine and Molecular Imaging (SNMMI) Dosimetry Task Force launched the 177Lu Dosimetry Challenge in 2021 (15). The challenge included 5 tasks (T1–T5). Three tasks investigated variability caused by different imaging protocols: serial SPECT/CT (T1), serial planar images (T2), or a hybrid approach (serial planar and 1 SPECT/CT image) (T3). Two additional tasks provided participants with volumes of interest (VOIs) (T4) and TIA maps (T5) with the aim of removing variability in segmentation and integration by removing sources of variability in the serial SPECT/CT workflow. The challenge did not address the impact of variability in ADs caused by image acquisition, calibration, or reconstruction.

The aim of the analysis presented in this work was to assess the source and magnitude of variability in AD estimates, for both organs and lesions, for the different tasks of the 177Lu dosimetry challenge and to inform standardization efforts.

MATERIALS AND METHODS

Patient Images and Data Collection

Datasets of 2 patients who underwent planar and SPECT/CT imaging at 4 time points after administration of 177Lu-DOTATATE therapy (16) were shared via the Deep Blue Data repository (https://deepblue.lib.umich.edu/data/collections/hm50ts030?locale=en) (17–21) of the University of Michigan. Sharing of patient images and data was approved by the University of Michigan Institutional Review Board, and both patients gave written informed consent. Table 1 summarizes the data provided. Maximum-intensity-projection images showing the provided VOIs are displayed in Figure 1.

View this table:
  • View inline
  • View popup
TABLE 1.

Patient Characteristics, Administered Activities, and Pre- and Post-therapeutic Imaging Information

FIGURE 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 1.

Maximum-intensity-projection images at 24 h after injection of 177Lu-DOTATATE for the 2 patients included in the dosimetry challenge. Contours of VOIs provided in T4 and T5 are shown. For patient A, average VOIs as measured from RTstructure files were 1,959 cm3 for liver, 247 cm3 for spleen, 467 cm3 for total kidney, 107 cm3 for lesion 1, and 3 cm3 for lesion 2; for patient B, they were 1,693 cm3 for liver, 229 cm3 for total kidney, 11 cm3 for lesion 1, 3 cm3 for lesion 2, 68 cm3 for lesion 3, and 22 cm3 for lesion 4. K = kidney.

Participants reported results on standardized spreadsheets tailored to each task. No lesion AD results were requested in T2 because of the overlap of lesions with organs on planar images. Separate spreadsheets were submitted for each patient. The data requested included information about methods, software used, intermediate values (e.g., VOI volumes, activities, and TIAs), and final ADs. Details of the methodology used for the challenge were given in our previous publication by Uribe et al. (15), including a full list of collected variables. Participants were encouraged to submit multiple calculations for the same patient using different types of software ormethodology (e.g., organ-based vs. voxel-based calculations). Participants were asked to briefly describe their dosimetry workflow in addition to the spreadsheet; unfortunately, however, only 1 participant mentioned partial-volume correction.

Data Collation

Data from all received submissions were extracted using the Python data analysis library (Pandas, version 1.3.5) and Python (version 3.9.5). All data were concatenated into a single data frame with columns corresponding to the specific variables collected. A full glossary of variables and column descriptions can be found in Supplemental Tables 1 and 2 (supplemental materials are available at http://jnm.snmjournals.org), as well as in the GitHub repository of the 177Lu Dosimetry Challenge (https://github.com/carluri/snmmi_dosimetry_challenge).

We curated the data, including identifying typographical, orientation (left/right ambiguity), and unit conversion issues and errors and evaluating data completeness. Participants were contacted for confirmation and clarification as needed. Results that were identified as containing mistakes in the calculations were removed from the analysis. A detailed description of the data curation process can be found in the supplemental material.

Statistical Analysis

Descriptive statistics, such as quartiles, means, and standard deviations (SDs) of ADs, were calculated separately for each task, patient, organ, and lesion. The quartile coefficient of dispersion (QCD) was calculated as the ratio of the difference between the 75th and 25th quartiles and the sum of the 75th and 25th quartiles of the data. The QCD was chosen since it is less sensitive to outliers in the data than is the coefficient of variation. For a normal distribution, the coefficient of variation is 1.4826 times the QCD.

A mixed-effects model was used to compare ADs among T1, T2, and T3 and, separately, among T1, T4, and T5. The model included the task as a fixed effect and the participant as a random effect. The analysis was performed separately for each patient and organ or lesion. Restricted maximum likelihood was applied to estimate the parameters in the models. All tests were 2-sided, and P values of 0.05 or less were considered to indicate statistically significance differences. The analysis was performed with Python (version 3.9.5) and R (version 4.1.0).

RESULTS

General Observations

We received 321 submissions from 51 institutions (Asia, 8; Australia, 2; Europe, 18; North America, 22; and South America, 1). Table 2 summarizes the submissions per task and patient. Some participants submitted multiple spreadsheets using different dosimetry methods (e.g. organ-level and voxelized).

View this table:
  • View inline
  • View popup
TABLE 2.

Number of Submissions per Task and per Patient

Both open-source and commercial dosimetry softwares were used, with details provided in Supplemental Figures 1 and 2. Approximately 27% of submissions used an in-house dose calculation approach.

Absorbed Doses

Figure 2 shows the distribution of mean ADs in organs and lesions from all submissions (after data curation). The numeric values of the descriptive statistics underlying Figure 2 can be found in Supplemental Table 1.

FIGURE 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 2.

Mean organ and lesion AD in Gray per task and per patient. Patient B had been splenectomized. T2 was based on planar images, T3 used hybrid imaging protocol of multiple planar images and 1 SPECT/CT image, T1 used multiple SPECT/CT images, T4 used multiple SPECT/CT images and provided VOIs, and T5 was based on provided TIA image and VOIs.

Table 3 shows the percentage difference calculated as the difference between the median of all submissions per task taking the median of all submissions from T1 as a reference, as well as the median of all submissions from T4 and T5 as a reference. Percentage differences were averaged separately for patients and for organs and lesions.

View this table:
  • View inline
  • View popup
TABLE 3.

Percentage Difference Between Medians of All Submissions per Task Relative to T1 and T4 Relative to T5

In general, ADs calculated from planar imaging protocols (T2) were lower than pure SPECT or hybrid protocols (Fig. 2). ADs for organs from T2 underestimated those from the pure SPECT protocol of T1 on average by 60%, ranging from −81% to −31% (Table 3). In contrast, the ADs of the hybrid protocol (T3) were similar to those of the pure SPECT/CT protocol (T1) for organs (on average 8% lower, ranging from −17% to 4%, Table 3). On average, the ADs for all organs were within ±10% of one another for T1, T4, and T5 (Table 3), suggesting no substantial bias between SPECT-based tasks. For lesions, a larger spread of values was observed (Fig. 2), but average percentage differences were within ±12%, ranging from −25% to +40% (Table 3). The provision of TIA maps in T5 yielded slightly smaller ranges in percentage difference of the medians of organ and lesion doses of T4 relative to T5 (Table 3).

Quantification of Variability

Figure 3 and Table 4 show the QCDs of ADs for all tasks averaged over organs and lesions. QCDs per organ and patient are provided in Supplemental Table 1. The QCDs for T2 were large, with an average of 69% for patient A and 46% for patient B (Table 4). The hybrid (T3) and pure SPECT (T1) protocols had similar QCDs for organs: on average 20% for patient A and 17% and 14% for patient B. Overall, the variability was reduced as segmentation and TIA activity data were provided (i.e., T1 to T4 to T5) (Figs. 2 and 3; Table 4). For the pure SPECT protocols (T1, T4, and T5), there was an overall reduction by a factor of about 1.5 in lesion QCD when VOIs were given to participants (T1 vs. T4, Table 4); the change for organs was smaller. Larger QCDs were observed for the smaller lesions (lesion 2 of patient A and lesion 2 of patient B) in T1, when segmentation and time integration were performed by participants. There was a substantial reduction in QCD for both organs and lesions when both VOIs and TIA maps were given to participants (T5, Table 4), resulting in QCDs of less than 7% for organs (excluding the liver of patient A) and less than 6% for lesions, independent of lesion size.

FIGURE 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 3.

QCD per task, organ, and patient.

View this table:
  • View inline
  • View popup
TABLE 4.

QCD per Task, Organ, and Patient

Statistical Analysis

The results of the statistical comparisons between ADs calculated for the different tasks are given in Table 5. In general, statistically significant differences were observed between T1 (pure SPECT) and T2 (planar imaging) but not between tasks that involved a SPECT scan, although there were a few exceptions.

View this table:
  • View inline
  • View popup
TABLE 5.

P Values for Comparisons of ADs Between Various Tasks and T1

DISCUSSION

All of the different steps in the dosimetry workflow potentially contribute to variability in AD estimates. First, the choice of imaging protocol (i.e., planar, SPECT, or hybrid imaging) can affect the measurement of the activity estimates that are the basis for dosimetry calculations. Second, the segmentation of organs and lesions can also affect the ADs. For voxel-based dosimetry, segmentation defines the spatial extent of VOIs where the AD is averaged, whereas for organ-based dosimetry it defines the organ mass and activity. The effect of segmentation on the AD is complicated since both the numerator (energy) and the denominator (mass) in the definition of dose (i.e., joules per kilogram) are affected. The estimation of total number of decays performed by curve fitting and time integration is influenced by the choice of fit function and the temporal limits of integration, neither of which is currently standardized. Finally, the choice of dosimetry method, software, and source of S-values or dose kernel can affect the final AD. The SNMMI 177Lu Dosimetry Challenge was designed to assess the variability in ADs caused by imaging protocol (T1, T2, and T3), segmentation (T4), and time-integration and dosimetry method (T5) on the final dosimetry results.

The planar protocol (T2) resulted in lower ADs by a factor of 2 and higher average QCDs than SPECT-based protocols, and these differences were statistically significant (Table 5). Lesion doses were not requested for T2 (planar protocols) because of the overlap with organs. Of note, T2 had the smallest number of submissions, with only 13 and 11 submissions for patients A and B, respectively.

The differences between ADs from the hybrid (T3) and pure SPECT (T1, T4, and T5) protocols were generally not statistically significant (Table 5). With the hybrid approach, the issues with overlapping structures are substantially reduced by use of the SPECT/CT image but can still affect the shape of the time–activity curve, especially for objects in high-uptake regions such as the lesions (Fig. 1). These results are consistent with previous reported data about the accuracy and precision of SPECT and hybrid protocols as compared with planar protocols (22–26).

The differences in ADs for the purely SPECT-based T1, T4, and T5 were generally small (Fig. 2; Table 3) and were not statistically significant for most organs and lesions (Table 5). However, statistically significant differences for the liver were observed, as can be explained by the presence of lesions in the livers of both patients (Fig. 1)—lesions might not have been excluded in the segmentation of healthy liver by all participants. This possibility suggests that standardization of segmentation methodologies should be considered to reduce variability. As demonstrated in Figure 2, the ranges of dose results were reduced when VOIs and TIA maps were provided, that is, comparing T1 with T4 and T5.

The largest average QCDs were found for the planar protocol, that is, T2 (Fig. 3; Table 4). We observed larger QCDs for the smaller lesions (Fig. 3, lesion 2 of patients A and B) than for the larger lesions. This result is expected given the difficulty and subjectivity associated with lesion segmentation and mass definition for those structures. This difficulty can further be related to the partial-volume effect, which is more pronounced for smaller lesions such as lesion 2 of patients A and B. In general, the QCDs, reflecting variation in ADs, were reduced as more information was provided to participants (i.e., T1 to T4 to T5), even for the challenging small lesions. Average organ QCDs changed little when VOIs were provided to participants (T1 to T4, Table 4), whereas lesion QCDs decreased by a factor of approximately 1.5. The QCDs of T5 were as low as 10% and 5% on average for all organs and 4% and 5% for all lesions of patients A and B, respectively (Table 4). The vast decrease in QCD between T1 and T5 from 30% to 4% for the small lesion 2 (∼3 mL) of patient A and from 39% to 6% for lesion 2 (∼3 mL) of patient B, compared with the moderate QCD decrease from 22% to 4% for lesion 1 (∼107 mL) of patient A and from 13% to 5% for lesion 3 (∼68 mL) of patient B, indicates that both segmentation and integration represent large sources of variability, especially for smaller objects. Segmentation can be further complicated when lesions are within an organ with little difference in contrast. The largest organ QCDs from all submissions for T5 were 32% and 7% for the liver of patients A and B, respectively (Fig. 3; Table 4; Supplemental Table 1). The size of this variation was unexpected given the data provided to the participants in T5. We attribute this variation to the presence of liver lesions and different decisions made by participants about what to include in the liver VOI (i.e., removing all lesions or only the lesions indicated by the challenge, Fig. 1). In general, the segmentation and TIA data provided in T4 and T5 substantially reduced variability as assessed by QCD (Fig. 3) and with respect to Figure 2, strongly suggesting that efforts to standardize segmentation (e.g., whether to include suspected lesions in normal tissues and whether to include the medulla and pelvis in kidney VOIs) may substantially reduce variability. Furthermore, providing TIA maps (i.e., standardization of fitting and integration) strongly reduced the variability in ADs and points to the integration approach as a source of substantial variability and a target for standardization.

The remaining variability in T5 can be attributed to several sources. Since this step included solely the conversion from TIA to AD, it may be related to differences in S values, dose kernels, or Monte Carlo simulations. Generally, these differences have been found to be small (<5%) (27–29), as is consistent with our independent findings for this dataset (30). Another potential source of variability is application of mass scaling to the S values (31) or density weighting to the dose kernels (27), both of which options are usually available in dosimetry software. Future analysis will focus on the effect of these factors on T5.

A limitation of this work is that it is based on only 2 patient datasets; this number was selected as a compromise between gaining more information on interpatient variability and the desire to attract a larger number of voluntary participants. The 2 patients chosen, however, illustrate some important characteristics and common challenges in the dosimetry workflow related to imaging protocol (planar vs. SPECT), segmentation, and integration. Specific characteristics of interest included a large tumor burden in the liver, significant differences between right and left kidney volumes, significant differences from standard phantom organ volumes, and lesion size and proximity to other high-uptake structures. In addition, some unintentional sources of variability are inherent in the design and implementation of the SNMMI 177Lu Dosimetry Challenge. For example, VOIs were provided both in the radiotherapy structure set (RTSTRUCT) of the DICOM standard and as voxelized masks to accommodate different capabilities in software available to participants. The process of voxelizing the RTSTRUCT resulted in different volume and activity estimates between the contours and the masks depending on the software used for the analysis. For example, one software package used by the organizers allowed the contours to include subvoxels, but the masks always contained complete voxels. Generally, the contour interpolation into subvoxels should be disabled for dosimetry purposes. This disabling resulted in average differences in volumes of 8% (range, 4%–14%) for organs and of 20% (range, 8%–37%) for lesions. Thus, the use of RTstructure by some participants and masks by others added to the variability in ADs. Similarly, VOIs were provided at each imaging time point, whereas the TIA map for T5 was generated using the first imaging time point as a reference. Applying to the TIA map the VOIs from a time point other than the first will consequently lead to different results. Although reflective of differences that may be observed at different sites in a clinical environment, the magnitude of the contribution of these sources to overall variability is likely specific to the design and implementation of the challenge. We did not expect that the provision of VOIs in RTstructure and mask format at multiple time points would affect the ADs and thus did not act to minimize this source of variability. Nevertheless, this difference when saving VOIs does highlight the potential for variability due to differences in software implementations and settings. Finally, this challenge did not address the impact of image acquisition, reconstruction, and quantification, which are also considered to be major contributors to variability of dose estimates.

The analysis of the SNMMI 177Lu Dosimetry Challenge data on overall dose variability elucidated several areas in which standardization or harmonization may be important to reduce variability across sites and methods. Our initial recommendations to reduce the variability of dose calculations based on the results of our analysis are as follows.

First, pure SPECT or hybrid SPECT/planar imaging protocols should be used for dosimetry instead of planar imaging–only protocols. The results presented here indicate that this recommendation would reduce variability and suggest a reduced bias.

Second, the development of segmentation guidelines for organs and lesions can help standardize the process and reduce the variability observed in this study. For example, guidance can be issued on how to deal with overlapping regions such as lesions in the liver and which regions of the kidney such as the renal pelvis should be included in the segmentation. In the meantime, we believe that publications should explicitly detail how these procedures are being performed to ensure that a correct comparison between results is being made.

Third, standardizing the use of fitting functions and integration methods would achieve some meaningful reduction in dose variability based on the reduction in variability from T4 to T5.

Fourth, more detailed and standardized reporting (32,33) of such details as dosimetry method (voxelized vs. organ-based), software used, appropriate use of mass scaling, and user decisions regarding the inclusion of kidney substructures should be developed to enable comparison of results from different centers and in different trials. A standardized nomenclature as emphasized in MIRD pamphlet 21 (34) would facilitate this reporting.

CONCLUSION

This analysis of all tasks of the SNMMI 177Lu Dosimetry Challenge highlights the need to move toward SPECT-based imaging protocols for dosimetry of radiopharmaceutical therapies. Standardizing segmentation and fitting methods and decisions is essential to reducing variability in AD. Removing these sources of variation from the dosimetry workflow reduced the variability to below 10% for organs and lesions.

DISCLOSURE

This work was partly supported by the SNMMI Value Initiative. Yuni Dewaraja acknowledges funding from grant R01CA240706 awarded by the National Cancer Institute for patient studies and resources made available by the University of Michigan Deep Blue Data Repository for data sharing. Eric Frey is a cofounder and part-owner of Rapid, LLC; receives royalty income from GE Healthcare; and acknowledges support from grants R44CA213782 and R01CA240779 awarded by the National Cancer Institute. Carlos Uribe acknowledges funding from Natural Science and Engineer Research Council (NSERC) discovery grant RGPIN-2021-02965. No potential conflict of interest relevant to this article was reported.

KEY POINTS

QUESTION: What is the overall variability in AD estimates in radiopharmaceutical therapy calculated using different methods at different centers?

PERTINENT FINDINGS: Planar imaging protocols resulted in substantially lower estimates of AD and the largest variability when compared with protocols that used SPECT/CT images. Overall variability using serial SPECT/CT as measured by QCD was less than 26% for organs and lesions. The results suggest that standardization of integration and segmentation may further reduce variability.

IMPLICATIONS FOR PATIENT CARE: Standardization of the dosimetry workflow will reduce variability in AD estimates and ultimately improve reliability—essential for dosimetry-based personalization of radiopharmaceutical therapies.

ACKNOWLEDGMENTS

We thank Bonnie Clarke, the director of research and discovery at the SNMMI, for her excellent support and her continuous help with data collection and communication of the challenge. Special thanks to all the participants: Aaron Scott, Adam Kesner, Albert Bartrés, Alessandro Desy, Alessia Milano, Anne-Laurène Wenger, Antoine Legrand, Arda Könik, Ashok Tiwari, Avery Peterson, Azadeh Akhavan, Benjamin Van, Carlos Montes Fuentes, Chae Moon Hong, Daniel McGowan, Daniele Pistone, David Adam, Diana McCrumb, Domenico Finocchiaro, Edoardo d’Andrea, Eric Brunner, Erin McKay, George Andl, Greta Mok, Heying Duan, Hina J. Shah, Ivan Yeung, Jacob Hesterman, Jeanette Libera-Körner, Joe Grudzinski, Johan Blakkisrud, Joshua Scheuermann, Juan Camilo Ocampo Ramos, Julia Brosch-Lenz, Keon Min Kim, Laetitia Imbert, Lara Bonney, Lawrence Hutton, Lorena Sandoval, Lukas Carter, Naoyuki Ukon, Natalie Cole, Nathaly Barbosa, Nuria Carrasco Vela, O Joo Hyun, Paulo Ferreira, Price Jackson, Rachele Danieli, Raquel Barquero, Richard Laforest, Sean McGurk, Seunggyun Ha, Shalini Subramanian, Sophia Wunder, Stéphanie Lamart, Stephen A. Graves, Su Bin Kim, Tay Young Soon, Teresa Pérez, Tiffany Beaumont, Valentina Ferri, Vikram Adhikarla, William Erwin, William Turner, and Yazdan Salimi

Footnotes

  • Published online Apr. 6, 2023.

  • © 2023 by the Society of Nuclear Medicine and Molecular Imaging.

Immediate Open Access: Creative Commons Attribution 4.0 International License (CC BY) allows users to share and adapt with attribution, excluding materials credited to previous publications. License: https://creativecommons.org/licenses/by/4.0/. Details: http://jnm.snmjournals.org/site/misc/permission.xhtml.

REFERENCES

  1. 1.↵
    1. Strosberg J,
    2. El-Haddad G,
    3. Wolin E,
    4. et al
    . Phase 3 trial of 177Lu-Dotatate for midgut neuroendocrine tumors. N Engl J Med. 2017;376:125–135.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Sartor O,
    2. de Bono J,
    3. Chi KN,
    4. et al
    . Lutetium-177-PSMA-617 for metastatic castration-resistant prostate cancer. N Engl J Med. 2021;385:1091–1103.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Fendler WP,
    2. Pabst KM,
    3. Kessler L,
    4. et al
    . Safety and efficacy of 90Y-FAPI-46 radioligand therapy in patients with advanced sarcoma and other cancer entities. Clin Cancer Res. 2022;28:4346–4353.
    OpenUrl
  4. 4.
    1. Fiedler L,
    2. Kellner M,
    3. Gosewisch A,
    4. et al
    . Evaluation of 177Lu[Lu]-CHX-A″-DTPA-6A10 Fab as a radioimmunotherapy agent targeting carbonic anhydrase XII. Nucl Med Biol. 2018;60:55–62.
    OpenUrl
  5. 5.
    1. Pooja D,
    2. Gunukula A,
    3. Gupta N,
    4. Adams DJ,
    5. Kulhari H
    . Bombesin receptors as potential targets for anticancer drug delivery and imaging. Int J Biochem Cell Biol. 2019;114:105567.
    OpenUrl
  6. 6.
    1. Wickstroem K,
    2. Karlsson J,
    3. Ellingsen C,
    4. et al
    . Synergistic effect of a HER2 targeted thorium-227 conjugate in combination with olaparib in a BRCA2 deficient xenograft model. Pharmaceuticals (Basel). 2019;12:155.
    OpenUrl
  7. 7.↵
    1. Grob NM,
    2. Schibli R,
    3. Béhé M,
    4. Mindt TL
    . Improved tumor-targeting with peptidomimetic analogs of minigastrin 177Lu-PP-F11N. Cancers (Basel). 2021;13:2629.
    OpenUrl
  8. 8.↵
    1. Ilhan H,
    2. Gosewisch A,
    3. Böning G,
    4. et al
    . Response to 225Ac-PSMA-I&T after failure of long-term 177Lu-PSMA RLT in mCRPC. Eur J Nucl Med Mol Imaging. 2021;48:1262–1263.
    OpenUrl
  9. 9.↵
    1. Kratochwil C,
    2. Bruchertseifer F,
    3. Rathke H,
    4. et al
    . Targeted α-therapy of metastatic castration-resistant prostate cancer with 225Ac-PSMA-617: swimmer-plot analysis suggests efficacy regarding duration of tumor control. J Nucl Med. 2018;59:795–802.
    OpenUrlAbstract/FREE Full Text
  10. 10.↵
    1. Danieli R,
    2. Milano A,
    3. Gallo S,
    4. et al
    . Personalized dosimetry in targeted radiation therapy: a look to methods, tools and critical aspects. J Pers Med. 2022;12:205.
    OpenUrl
  11. 11.↵
    1. Garin E,
    2. Tselikas L,
    3. Guiu B,
    4. et al
    . Personalised versus standard dosimetry approach of selective internal radiation therapy in patients with locally advanced hepatocellular carcinoma (DOSISPHERE-01): a randomised, multicentre, open-label phase 2 trial. Lancet Gastroenterol Hepatol. 2021;6:17–29.
    OpenUrl
  12. 12.↵
    1. Dewaraja YK,
    2. Frey EC,
    3. Sgouros G,
    4. et al
    . MIRD pamphlet no. 23: quantitative SPECT for patient-specific 3-dimensional dosimetry in internal radionuclide therapy. J Nucl Med. 2012;53:1310–1325.
    OpenUrlAbstract/FREE Full Text
  13. 13.
    1. Capala J,
    2. Graves SA,
    3. Scott A,
    4. et al
    . Dosimetry for radiopharmaceutical therapy: current practices and commercial resources. J Nucl Med. 2021;62(suppl 3):3S–11S.
    OpenUrlAbstract/FREE Full Text
  14. 14.↵
    1. Mora-Ramirez E,
    2. Santoro L,
    3. Cassol E,
    4. et al
    . Comparison of commercial dosimetric software platforms in patients treated with 177Lu-DOTATATE for peptide receptor radionuclide therapy. Med Phys. 2020;47:4602–4615.
    OpenUrl
  15. 15.↵
    1. Uribe C,
    2. Peterson A,
    3. Van B,
    4. et al
    . An international study of factors affecting variability of dosimetry calculations, part 1: design and early results of the SNMMI dosimetry challenge. J Nucl Med. 2021;62(suppl 3):36S–47S.
    OpenUrlAbstract/FREE Full Text
  16. 16.↵
    1. Dewaraja YK,
    2. Mirando DM,
    3. Peterson A,
    4. et al
    . A pipeline for automated voxel dosimetry: application in patients with multi-SPECT/CT imaging following 177Lu peptide receptor radionuclide therapy. J Nucl Med. 2022;63:1665–1672.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Dewaraja YK,
    2. Van BJ
    . Lu-177 DOTATATE anonymized patient datasets: pre-therapy diagnostic images. Deep Blue Data website. https://doi.org/10.7302/vqmy-g059. Published February 10, 2021. Modified November 18, 2022. Accessed March 15, 2023.
  18. 18.
    1. Dewaraja YK,
    2. Van BJ
    . Lu-177 DOTATATE anonymized patient datasets: multi-time point Lu-177 planar whole body scans. Deep Blue Data website. https://doi.org/10.7302/y4xd-s758. Published February 10, 2021. Modified November 18, 2022. Accessed March 15, 2023
  19. 19.
    1. Dewaraja YK,
    2. Van BJ
    . Lu-177 DOTATATE anonymized patient datasets: multi-time point Lu-177 SPECT/CT scans. Deep Blue Data website. https://doi.org/107302/0n8e-rz46. Published February 10, 2021. Modified November 18, 2022. Accessed March 15, 2023.
  20. 20.
    1. Dewaraja YK,
    2. Van BJ
    . Lu-177 DOTATATE anonymized patient datasets: lesion and organ volumes of interest. Deep Blue Data website. https://doi.org/10.7302/vhrh-qg23. Published May 19, 2021. Modified November 17, 2022. Accessed March 15, 2023.
  21. 21.↵
    1. Dewaraja YK,
    2. Van BJ
    . Lu-177 DOTATATE anonymized patient datasets: Lu-177 SPECT time integrated activity maps. Deep Blue Data website. https://doi.org/10.7302/9nct-bk44. Published August 23, 2021. Modified November 18, 2022. Accessed March 15, 2023.
  22. 22.↵
    1. He B,
    2. Du Y,
    3. Song X,
    4. Segars WP,
    5. Frey EC
    . A Monte Carlo and physical phantom evaluation of quantitative In-111 SPECT. Phys Med Biol. 2005;50:4169–4185.
    OpenUrlCrossRefPubMed
  23. 23.
    1. He B,
    2. Frey EC
    . Comparison of conventional, model-based quantitative planar, and quantitative SPECT image processing methods for organ activity estimation using In-111 agents. Phys Med Biol. 2006;51:3967–3981.
    OpenUrlCrossRefPubMed
  24. 24.
    1. He B,
    2. Wahl RL,
    3. Du Y,
    4. et al
    . Comparison of residence time estimation methods for radioimmunotherapy dosimetry and treatment planning: Monte Carlo simulation studies. IEEE Trans Med Imaging. 2008;27:521–530.
    OpenUrlCrossRefPubMed
  25. 25.
    1. Belli ML,
    2. Mezzenga E,
    3. Di Iorio V,
    4. et al
    . A whole body dosimetry protocol for peptide-receptor radionuclide therapy (PRRT): 2D planar image and hybrid 2D+3D SPECT/CT image methods. J Vis Exp. 2020;158:e60477.
    OpenUrl
  26. 26.↵
    1. Garkavij M,
    2. Nickel M,
    3. Sjögreen-Gleisner K,
    4. et al
    . 177Lu-[DOTA0,Tyr3] octreotate therapy in patients with disseminated neuroendocrine tumors: analysis of dosimetry with impact on future therapeutic strategy. Cancer. 2010;116:1084–1092.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Dieudonné A,
    2. Hobbs RF,
    3. Lebtahi R,
    4. et al
    . Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct Monte Carlo methods. J Nucl Med. 2013;54:236–243.
    OpenUrlAbstract/FREE Full Text
  28. 28.
    1. Brosch-Lenz J,
    2. Uribe C,
    3. Gosewisch A,
    4. et al
    . Influence of dosimetry method on bone lesion absorbed dose estimates in PSMA therapy: application to mCRPC patients receiving Lu-177-PSMA-I&T. EJNMMI Phys. 2021;8:26.
    OpenUrl
  29. 29.↵
    1. Grimes J,
    2. Celler A
    . Comparison of internal dose estimates obtained using organ-level, voxel S value, and Monte Carlo techniques. Med Phys. 2014;41:092501.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Brosch-Lenz J,
    2. Rahmim A,
    3. Sunderland J,
    4. Frey E,
    5. Dewaraja YK,
    6. Uribe C
    . Impact of dosimetry method on healthy organ and tumor absorbed dose estimates for lutetium-177-DOTATATE therapy of neuroendocrine tumors. Eur J Nucl Med Mol Imaging. 2022;49(suppl 1):S393–S394.
    OpenUrl
  31. 31.↵
    1. Howell RW,
    2. Wessels BW,
    3. Loevinger R,
    4. et al
    . The MIRD perspective 1999: medical internal radiation dose committee. J Nucl Med. 1999;40:3S–10S.
    OpenUrlFREE Full Text
  32. 32.↵
    1. Lassmann M,
    2. Chiesa C,
    3. Flux G,
    4. Bardiès M
    ; EANM Dosimetry Committee. EANM Dosimetry Committee guidance document: good practice of clinical dosimetry reporting. Eur J Nucl Med Mol Imaging. 2011;38:192–200.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Sgouros G,
    2. Bolch WE,
    3. Chiti A,
    4. et al
    . ICRU REPORT 96, dosimetry-guided radiopharmaceutical therapy. J ICRU. 2021;21:1–212.
    OpenUrl
  34. 34.↵
    1. Bolch WE,
    2. Eckerman KF,
    3. Sgouros G,
    4. Thomas SR
    . MIRD pamphlet no. 21: a generalized schema for radiopharmaceutical dosimetry—standardization of nomenclature. J Nucl Med. 2009;50:477–484.
    OpenUrlAbstract/FREE Full Text
  • Received for publication October 25, 2022.
  • Revision received February 7, 2023.
PreviousNext
Back to top

In this issue

Journal of Nuclear Medicine: 64 (7)
Journal of Nuclear Medicine
Vol. 64, Issue 7
July 1, 2023
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Complete Issue (PDF)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Journal of Nuclear Medicine.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 2: Overall Variabilities in Absorbed Dose
(Your Name) has sent you a message from Journal of Nuclear Medicine
(Your Name) thought you would like to see the Journal of Nuclear Medicine web site.
Citation Tools
An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 2: Overall Variabilities in Absorbed Dose
Julia Brosch-Lenz, Suqi Ke, Hao Wang, Eric Frey, Yuni K. Dewaraja, John Sunderland, Carlos Uribe
Journal of Nuclear Medicine Jul 2023, 64 (7) 1109-1116; DOI: 10.2967/jnumed.122.265094

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 2: Overall Variabilities in Absorbed Dose
Julia Brosch-Lenz, Suqi Ke, Hao Wang, Eric Frey, Yuni K. Dewaraja, John Sunderland, Carlos Uribe
Journal of Nuclear Medicine Jul 2023, 64 (7) 1109-1116; DOI: 10.2967/jnumed.122.265094
Twitter logo Facebook logo LinkedIn logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One
Bookmark this article

Jump to section

  • Article
    • Visual Abstract
    • Abstract
    • MATERIALS AND METHODS
    • RESULTS
    • DISCUSSION
    • CONCLUSION
    • DISCLOSURE
    • ACKNOWLEDGMENTS
    • Footnotes
    • REFERENCES
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • PDF

Related Articles

  • PubMed
  • Google Scholar

Cited By...

  • Computational Nuclear Oncology Toward Precision Radiopharmaceutical Therapies: Current Tools, Techniques, and Uncharted Territories
  • An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 4: Impact of Fitting Functions in Estimated Absorbed Doses
  • Dosimetry Software for Theranostic Applications: Current Capabilities and Future Prospects
  • An International Study of Factors Affecting Variability of Dosimetry Calculations, Part 3: Contribution from Calculating Absorbed Dose from Time-Integrated Activity
  • Google Scholar

More in this TOC Section

  • First-in-Human Study of 18F-Labeled PET Tracer for Glutamate AMPA Receptor [18F]K-40: A Derivative of [11C]K-2
  • Detection of HER2-Low Lesions Using HER2-Targeted PET Imaging in Patients with Metastatic Breast Cancer: A Paired HER2 PET and Tumor Biopsy Analysis
  • [11C]Carfentanil PET Whole-Body Imaging of μ-Opioid Receptors: A First in-Human Study
Show more Clinical Investigation

Similar Articles

Keywords

  • radiopharmaceutical therapy
  • dosimetry
  • variability
  • standardization
  • dosimetry challenge
SNMMI

© 2025 SNMMI

Powered by HighWire