Proof-of-concept of DosiTest: A virtual multicentric clinical trial for assessing uncertainties in molecular radiotherapy dosimetry

Quality Assurance Considerations in Radiopharmaceutical Therapy Dosimetry Using PLANETDose: An International Atomic Energy Agency Study

Implementation of radiopharmaceutical therapy dosimetry varies depending on the clinical application, dosimetry protocol, software, and ultimately the operator. Assessing clinical dosimetry accuracy and precision is therefore a challenging task. This work emphasizes some pitfalls encountered during a structured analysis, performed on a single-patient dataset consisting of SPECT/CT images by various participants using a standard protocol and clinically approved commercial software. Methods: The clinical dataset consisted of the dosimetric study of a patient administered with [177Lu]Lu-DOTATATE at Tygerberg Hospital, South Africa, as a part of International Atomic Energy Agency-coordinated research project E23005. SPECT/CT images were acquired at 5 time points postinjection. Patient and calibration images were reconstructed on a workstation, and a calibration factor of 122.6 Bq/count was derived independently and provided to the participants. A standard dosimetric protocol was defined, and PLANETDose (version 3.1.1) software was installed at 9 centers to perform the dosimetry of 3 treatment cycles. The protocol included rigid image registration, segmentation (semimanual for organs, activity threshold for tumors), and dose voxel kernel convolution of activity followed by absorbed dose (AD) rate integration to obtain the ADs. Iterations of the protocol were performed by participants individually and within collective training, the results of which were analyzed for dosimetric variability, as well as for quality assurance and error analysis. Intermediary checkpoints were developed to understand possible sources of variation and to differentiate user error from legitimate user variability. Results: Initial dosimetric results for organs (liver and kidneys) and lesions showed considerable interoperator variability. Not only was the generation of intermediate checkpoints such as total counts, volumes, and activity required, but also activity-to-count ratio, activity concentration, and AD rate-to-activity concentration ratio to determine the source of variability. Conclusion: When the same patient dataset was analyzed using the same dosimetry procedure and software, significant disparities were observed in the results despite multiple sessions of training and feedback. Variations due to human error could be minimized or avoided by performing intensive training sessions, establishing intermediate checkpoints, conducting sanity checks, and cross-validating results across physicists or with standardized datasets. This finding promotes the development of quality assurance in clinical dosimetry.

The contest between internal and external-beam dosimetry: The Zeno's paradox of Achilles and the tortoise

Radionuclide therapy, also called molecular radiotherapy (MRT), has come of age, with several novel radiopharmaceuticals being approved for clinical use or under development in the last decade. External beam radiotherapy (EBRT) is a well-established treatment modality, with about half of all oncologic patients expected to receive at least one external radiation treatment over their disease course. The efficacy and the toxicity of both types of treatment rely on the interaction of radiation with biological tissues. Dosimetry played a fundamental role in the scientific and technological evolution of EBRT, and absorbed doses to the target and to the organs at risk are calculated on a routine basis. In contrast, in MRT the usefulness of internal dosimetry has long been questioned, and a structured path to include absorbed dose calculation is missing. However, following a similar route of development as EBRT, MRT treatments could probably be optimized in a significant proportion of patients, likely based on dosimetry and radiobiology. In the present paper we describe the differences and the similarities between internal and external-beam dosimetry in the context of radiation treatments, and we retrace the main stages of their development over the last decades.

Time-Activity data fitting in molecular Radiotherapy: Methodology and pitfalls

Absorbed radiation doses are essential in assessing the effects, e.g. safety and efficacy, of radiopharmaceutical therapy (RPT). Patient-specific absorbed dose calculations in the target or the organ at risk require multiple inputs. These include the number of disintegrations in the organ, i.e. the time-integrated activities (TIAs) of the organs, as well as other parameters describing the process of radiation energy deposition in the target tissue (i.e. mean energy per disintegration, radiation dose constants, etc). TIAs are then estimated by incorporating the area under the radiopharmaceutical's time-activity curve (TAC), which can be obtained by quantitative measurements of the biokinetics in the patient (typically based on imaging data such as planar scintigraphy, SPECT/CT, PET/CT, or blood and urine samples). The process of TAC determination/calculation for RPT generally depends on the user, e.g., the chosen number and schedule of measured time points, the selection of the fit function, the error model for the data and the fit algorithm. These decisions can strongly affect the final TIA values and thus the accuracy of calculated absorbed doses. Despite the high clinical importance of the TIA values, there is currently no consensus on processing time-activity data or even a clear understanding of the influence of uncertainties and variations in personalised RPT dosimetry related to user-dependent TAC calculation. As a first step towards minimising site-dependent variability in RPT dosimetry, this work provides an overview of quality assurance and uncertainty management considerations of the TIA estimation.

Patient-specific biokinetics and hybrid 2D/3D approach integration in OEDIPE software: Application to radioiodine therapy

BACKGROUND

The progression of targeted radionuclide therapy requires the development of dosimetry software accounting for patient-specific biokinetics. New functionalities were thus developed in the OEDIPE software, to deal with multiple 3D images or multiple planar images and a SPECT image.

MATERIEL & METHOD

Methods were implemented to recover patient biokinetics in volumes of interest. If several 3D SPECT images are available, they are registered to a reference CT scan. When several planar images and a single SPECT are available, the planar images are registered to the SPECT and counts of the planar images converted to activity. To validate these developments, six SPECT/CT and planar images of a Jaszczak phantom containing I-131 were acquired at different dates. Cumulated activity was estimated in each sphere using the SPECT/CT images only or the planar series associated to one SPECT/CT. Biokinetics and doses in lesions and in the lungs of a patient treated with I-131 for differentiated thyroid cancer were then estimated using four planar images and a SPECT/CT scan. Whole-body retention data were used to compare the biokinetics obtained from the planar and SPECT data.

RESULTS

Activities and cumulated activities estimated using OEDIPE in the phantom spheres agreed well with the reference values for both approaches. Results obtained for the patient compared well with those derived from whole-body retention data.

CONCLUSION

The implemented features allow automatic evaluation of patient-specific biokinetics from different series of patient images, enabling patient-specific dosimetry without the need for external software to estimate the cumulated activities in different VOIs.

Development of a validation imaging dataset for Molecular Radiotherapy dosimetry multicenter intercomparison exercises based on anthropomorphic phantoms

Validation of a Molecular Radiotherapy (MRT) dosimetry system requires imaging data for which an accompanying "ground truth" pharmacokinetic model and absorbed dose calculation are known.

METHODS

We present a methodology for production of a validation dataset for image based ¹⁷⁷Lu dotatate dosimetry calculations. A pharmacokinetic model is presented with activity concentrations corresponding to common imaging timepoints. Anthropomorphic 3D printed phantoms, corresponding to the organs at risk, have been developed to provide SPECT/CT and Whole Body imaging with known organ activities corresponding to common clinical timepoints.

RESULTS

Results for the accuracy of phantom filling reproduce the activity concentrations from the pharmacokinetic model for all timepoints and organs within measurement uncertainties, with a mean deviation of 0.6(8)%. The imaging dataset, ancillary data and phantoms designs are provided as a source of well characterized input data for the validation of clinical MRT dosimetry systems.

CONCLUSIONS

The combination of pharmacokinetic modelling with the use of anthropomorphic 3D printed phantoms are a promising procedure to provide data for the validation of Molecular Radiotherapy Dosimetry systems, allowing multicentre comparisons.

Predicting subacute hematological toxicity of 177Lu-DOTATATE therapy using healthy organs' uptake on post-treatment quantitative SPECT: A pilot study

The aim is to investigate the usefulness of 177Lu-DOTA-0-Tyr3-Octreotate (DOTATATE) healthy organs' (spleen, kidneys, bone marrow) standard uptake value for the prediction of subacute hematological toxicity in patients undergoing 177Lu-DOTATATE treatment. All patients referred from January 2021 to May 2022 for 177Lu-DOTATATE treatment were retrospectively screened. For each treatment session, baseline clinical data including age, sex, weight, delay between 177Lu-DOTATATE treatment and last cold somatostatin analogue intake were collected. Mean standardized uptake value (SUVmean) of healthy organs was measured and analyzed by generalized linear mixed effect models. Outcomes (significant decrease of platelets, hemoglobin levels and neutrophils) were assessed 1 month later, considering their within-subject biological coefficient of variation, published by the European Federation of Clinical Chemistry and Laboratory Medicine. A total of 9 patients (33 treatment sessions) were included. No predictive factors were identified for platelet and neutrophil decrease. Splenic SUVmean was found to be a significant predictor of hemoglobin levels decrease. Using an optimal threshold of ≥6.22, derived sensitivity and specificity to predict hemoglobin decrease were 85.7% [46.4; 99.0] and 76.9% [57.5; 89.2] respectively with an accuracy of 82.4%. Although not significantly predictive of hematological toxicity, bone marrow SUVmean and renal SUVmean were correlated with splenic SUVmean. Quantitative single photon emission computed tomography and healthy organs analysis might help to foresee hematological subacute toxicity in patients undergoing 177Lu-DOTATATE treatment and improve patient management.