SPECT/CT Image-Derived Absorbed Dose to Red Marrow Correlates with Hematologic Toxicity in Patients Treated with [177Lu]Lu-DOTATATE.

Hematologic toxicity, although often transient, is the most common limiting adverse effect during somatostatin peptide receptor radionuclide therapy. This study investigated the association between Monte Carlo-derived absorbed dose to the red marrow (RM) and hematologic toxicity in patients being treated for their neuroendocrine tumors. Methods: Twenty patients each receiving 4 treatment cycles of [177Lu]Lu-DOTATATE were included. Multiple-time-point 177Lu SPECT/CT imaging-based RM dosimetry was performed using an artificial intelligence-driven workflow to segment vertebral spongiosa within the field of view (FOV). This workflow was coupled with an in-house macroscale/microscale Monte Carlo code that incorporates a spongiosa microstructure model. Absorbed dose estimates to RM in lumbar and thoracic vertebrae within the FOV, considered as representations of the whole-body RM absorbed dose, were correlated with hematologic toxicity markers at about 8 wk after each cycle and at 3- and 6-mo follow-up after completion of all cycles. Results: The median of absorbed dose to RM in lumbar and thoracic vertebrae within the FOV (D median,vertebrae) ranged from 0.019 to 0.11 Gy/GBq. The median of cumulative absorbed dose across all 4 cycles was 1.3 Gy (range, 0.6-2.5 Gy). Hematologic toxicity was generally mild, with no grade 2 or higher toxicity for platelets, neutrophils, or hemoglobin. However, there was a decline in blood counts over time, with a fractional value relative to baseline at 6 mo of 74%, 97%, 57%, and 97%, for platelets, neutrophils, lymphocytes, and hemoglobin, respectively. Statistically significant correlations were found between a subset of hematologic toxicity markers and RM absorbed doses, both during treatment and at 3- and 6-mo follow-up. This included a correlation between the platelet count relative to baseline at 6-mo follow up: D median,vertebrae (r = -0.64, P = 0.015), D median,lumbar (r = -0.72, P = 0.0038), D median,thoracic (r = -0.58, P = 0.029), and D average,vertebrae (r = -0.66, P = 0.010), where D median,lumbar and D median,thoracic are median absorbed dose to the RM in the lumbar and thoracic vertebrae, respectively, within the FOV and D average,vertebrae is the mass-weighted average absorbed dose of all vertebrae. Conclusion: This study found a significant correlation between image-derived absorbed dose to the RM and hematologic toxicity, including a relative reduction of platelets at 6-mo follow up. These findings indicate that absorbed dose to the RM can potentially be used to understand and manage hematologic toxicity in peptide receptor radionuclide therapy.

Single-Time-Point Imaging for Dosimetry After [177Lu]Lu-DOTATATE: Accuracy of Existing Methods and Novel Data-Driven Models for Reducing Sensitivity to Time-Point Selection.

Estimation of the time-integrated activity (TIA) for dosimetry from imaging at a single time point (STP) facilitates the clinical translation of dosimetry-guided radiopharmaceutical therapy. However, the accuracy of the STP methods for TIA estimation varies on the basis of time-point selection. We constructed patient data-driven regression models to reduce the sensitivity to time-point selection and to compare these new models with commonly used STP methods. Methods: SPECT/CT performed at time period (TP) 1 (3-5 h), TP2 (days 1-2), TP3 (days 3-5), and TP4 (days 6-8) after cycle 1 of [177Lu]Lu-DOTATATE therapy involved 27 patients with 100 segmented tumors and 54 kidneys. Influenced by the previous physics-based STP models of Madsen et al. and Hänscheid et al., we constructed an STP prediction expression, TIA = A(t) × g(t), in a SPECT data-driven way (model 1), in which A(t) is the observed activity at imaging time t, and the curve, g(t), is estimated with a nonparametric generalized additive model by minimizing the normalized mean square error relative to the TIA derived from 4-time-point SPECT (reference TIA). Furthermore, we fit a generalized additive model that incorporates baseline biomarkers as auxiliary data in addition to the single activity measurement (model 2). Leave-one-out cross validation was performed to evaluate STP models using mean absolute error (MAE) and mean square error between the predicted and reference TIA. Results: At days 3-5, all evaluated STP methods performed very well, with an MAE of less than 7% (between-patient SD of <10%) for both kidneys and tumors. At other TPs, the Madsen method and data-driven models 1 and 2 performed reasonably well (MAEs < 17% for kidneys and < 32% for tumors), whereas the error with the Hänscheid method was substantially higher. The proof of concept of adding baseline biomarkers to the prediction model was demonstrated and showed a moderate enhancement at TP1, especially for estimating kidney TIA (MAE ± SD from 15.6% ± 1.3% to 11.8% ± 1.0%). Evaluations on 500 virtual patients using clinically relevant time-activity simulations showed a similar performance. Conclusion: The performance of the Madsen method and proposed data-driven models is less sensitive to TP selection than is the Hänscheid method. At the earliest TP, which is the most practical, the model incorporating baseline biomarkers outperforms other methods that rely only on the single activity measurement.

Accuracy and uncertainty analysis of reduced time point imaging effect on time-integrated activity for 177Lu-DOTATATE PRRT in patients and clinically realistic simulations.

BACKGROUND: Dosimetry promises many advantages for radiopharmaceutical therapies but repeat post-therapy imaging for dosimetry can burden both patients and clinics. Recent applications of reduced time point imaging for time-integrated activity (TIA) determination for internal dosimetry following 177Lu-DOTATATE peptide receptor radionuclide therapy have shown promising results that allow for the simplification of patient-specific dosimetry. However, factors such as scheduling can lead to sub-optimal imaging time points, but the resulting impact on dosimetry accuracy is still under investigation. We use four-time point 177Lu SPECT/CT data for a cohort of patients treated at our clinic to perform a comprehensive analysis of the error and variability in time-integrated activity when reduced time point methods with various combinations of sampling points are employed. METHODS: The study includes 28 patients with gastroenteropancreatic neuroendocrine tumors who underwent post-therapy SPECT/CT imaging at approximately 4, 24, 96, and 168 h post-therapy (p.t.) following the first cycle of 177Lu-DOTATATE. The healthy liver, left/right kidney, spleen and up to 5 index tumors were delineated for each patient. Time-activity curves were fit with either monoexponential or biexponential functions for each structure, based on the Akaike information criterion. This fitting was performed using all 4 time points as a reference and various combinations of 2 and 3 time points to determine optimal imaging schedules and associated errors. 2 commonly used methods of single time point (STP) TIA estimation are also evaluated. A simulation study was also performed with data generated by sampling curve fit parameters from log-normal distributions derived from the clinical data and adding realistic measurement noise to sampled activities. For both clinical and simulation studies, error and variability in TIA estimates were estimated with various sampling schedules. RESULTS: The optimal post-therapy imaging time period for STP estimates of TIA was found to be 3-5 days (71-126 h) p.t. for tumor and organs, with one exception of 6-8 days (144-194 h) p.t. for spleen with one STP approach. At the optimal time point, STP estimates give mean percent errors (MPE) within ± 5% and SD < 9% across all structures with largest magnitude error for kidney TIA (MPE = - 4.1%) and highest variability also for kidney TIA (SD = 8.4%). The optimal sampling schedule for 2TP estimates of TIA is 1-2 days (21-52 h) p.t. followed by 3-5 days (71-126 h) p.t. for kidney, tumor, and spleen. Using the optimal sampling schedule, the largest magnitude MPE for 2TP estimates is 1.2% for spleen and highest variability is in tumor with SD = 5.8%. The optimal sampling schedule for 3TP estimates of TIA is 1-2 days (21-52 h) p.t. followed by 3-5 days (71-126 h) p.t. and 6-8 days (144-194 h) p.t. for all structures. Using the optimal sampling schedule, the largest magnitude MPE for 3TP estimates is 2.5% for spleen and highest variability is in tumor with SD = 2.1%. Simulated patient results corroborate these findings with similar optimal sampling schedules and errors. Many sub-optimal reduced time point sampling schedules also exhibit low error and variability. CONCLUSIONS: We show that reduced time point methods can be used to achieve acceptable average TIA errors over a wide range of imaging time points and sampling schedules while maintaining low uncertainty. This information can improve the feasibility of dosimetry for 177Lu-DOTATATE and elucidate the uncertainty associated with non-ideal conditions.

The predictive value of pretherapy [68Ga]Ga-DOTA-TATE PET and biomarkers in [177Lu]Lu-PRRT tumor dosimetry.

PURPOSE: Metastatic neuroendocrine tumors (NETs) overexpressing type 2 somatostatin receptors are the target for peptide receptor radionuclide therapy (PRRT) through the theragnostic pair of 68Ga/177Lu-DOTATATE. The main purpose of this study was to develop machine learning models to predict therapeutic tumor dose using pre therapy 68Ga -PET and clinicopathological biomarkers. METHODS: We retrospectively analyzed 90 segmented metastatic NETs from 25 patients (M14/F11, age 63.7 ± 9.5, range 38-76) treated by 177Lu-DOTATATE at our institute. Patients underwent both pretherapy [68Ga]Ga-DOTA-TATE PET/CT and four timepoints SPECT/CT at ~ 4, 24, 96, and 168 h post-177Lu-DOTATATE infusion. Tumors were segmented by a radiologist on baseline CT or MRI and transferred to co-registered PET/CT and SPECT/CT, and normal organs were segmented by deep learning-based method on CT of the PET and SPECT. The SUV metrics and tumor-to-normal tissue SUV ratios (SUV_TNRs) were calculated from 68Ga -PET at the contour-level. Posttherapy dosimetry was performed based on the co-registration of SPECT/CTs to generate time-integrated-activity, followed by an in-house Monte Carlo-based absorbed dose estimation. The correlation between delivered 177Lu Tumor absorbed dose and PET-derived metrics along with baseline clinicopathological biomarkers (such as Creatinine, Chromogranin A and prior therapies) were evaluated. Multiple interpretable machine-learning algorithms were developed to predict tumor dose using these pretherapy information. Model performance on a nested tenfold cross-validation was evaluated in terms of coefficient of determination (R2), mean-absolute-error (MAE), and mean-relative-absolute-error (MRAE). RESULTS: SUVmean showed a significant correlation (q-value < 0.05) with absorbed dose (Spearman ρ = 0.64), followed by TLSUVmean (SUVmean of total-lesion-burden) and SUVpeak (ρ = 0.45 and 0.41, respectively). The predictive value of PET-SUVmean in estimation of posttherapy absorbed dose was stronger compared to PET-SUVpeak, and SUV_TNRs in terms of univariate analysis (R2 = 0.28 vs. R2 ≤ 0.12). An optimal trivariate random forest model composed of SUVmean, TLSUVmean, and total liver SUVmean (normal and tumoral liver) provided the best performance in tumor dose prediction with R2 = 0.64, MAE = 0.73 Gy/GBq, and MRAE = 0.2. CONCLUSION: Our preliminary results demonstrate the feasibility of using baseline PET images for prediction of absorbed dose prior to 177Lu-PRRT. Machine learning models combining multiple PET-based metrics performed better than using a single SUV value and using other investigated clinicopathological biomarkers. Developing such quantitative models forms the groundwork for the role of 68Ga -PET not only for the implementation of personalized treatment planning but also for patient stratification in the era of precision medicine.

Accuracy and uncertainty analysis of reduced time point imaging effect on time-integrated activity for 177Lu-DOTATATE PRRT in clinical patients and realistic simulations.

Background. Dosimetry promises many advantages for radiopharmaceutical therapies but repeat post-therapy imaging for dosimetry can burden both patients and clinics. Recent applications of reduced time point imaging for time-integrated activity (TIA) determination for internal dosimetry following 177 Lu-DOTATATE peptide receptor radionuclide therapy have shown promising results that allow for the simplification of patient-specific dosimetry. However, factors such as scheduling can lead to undesirable imaging time points, but the resulting impact on dosimetry accuracy is unknown. We use four-time point 177 Lu SPECT/CT data for a cohort of patients treated at our clinic to perform a comprehensive analysis of the error and variability in time-integrated activity when reduced time point methods with various combination of sampling points are employed. Methods. The study includes 28 patients with gastroenteropancreatic neuroendocrine tumors who underwent post-therapy SPECT/CT imaging at approximately 4, 24, 96, and 168 hours post-therapy (p.t.) following the first cycle of 177 Lu-DOTATATE. The healthy liver, left/right kidney, spleen and up to 5 index tumors were delineated for each patient. Time-activity curves were fit with either monoexponential or biexponential functions for each structure, based on the Akaike information criterion. This fitting was performed using all 4 time points as a reference and various combinations of 2 and 3 time points to determine optimal imaging schedules and associated errors. 2 commonly used methods of single time point (STP) TIA estimation are also evaluated. A simulation study was also performed with data generated by sampling curve fit parameters from log-normal distributions derived from the clinical data and adding realistic measurement noise to sampled activities. For both clinical and simulation studies, error and variability in TIA estimates were estimated with various sampling schedules. Results . The optimal post-therapy imaging time period for STP estimates of TIA was found to be 3-5 days (71-126 h) p.t. for tumor and organs, with one exception of 6-8 days (144-194 h) p.t. for spleen with one STP approach. At the optimal time point, STP estimates give mean percent errors (MPE) within +/-5% and SD < 9% across all structures with largest magnitude error for kidney TIA (MPE=-4.1%) and highest variability also for kidney TIA (SD=8.4%). The optimal sampling schedule for 2TP estimates of TIA is 1-2 days (21-52 h) p.t. followed by 3-5 days (71-126 h) p.t. for kidney, tumor, and spleen. Using the optimal sampling schedule, the largest magnitude MPE for 2TP estimates is 1.2% for spleen and highest variability is in tumor with SD=5.8%. The optimal sampling schedule for 3TP estimates of TIA is 1-2 days (21-52 h) p.t. followed by 3-5 days (71-126 h) p.t. and 6-8 days (144-194 h) p.t. for all structures. Using the optimal sampling schedule, the largest magnitude MPE for 3TP estimates is 2.5% for spleen and highest variability is in tumor with SD=2.1%. Simulated patient results corroborate these findings with similar optimal sampling schedules and errors. Many sub-optimal reduced time point sampling schedules also exhibit low error and variability. Conclusions. We show that reduced time point methods can be used to achieve acceptable average TIA errors over a wide range of imaging time points and sampling schedules while maintaining low uncertainty. This information can improve the feasibility of dosimetry for 177 Lu-DOTATATE and elucidate the uncertainty associated with non-ideal conditions.

177Lu-DOTATATE Theranostics: Predicting Renal Dosimetry From Pretherapy 68Ga-DOTATATE PET and Clinical Biomarkers.

PURPOSE: Pretreatment predictions of absorbed doses can be especially valuable for patient selection and dosimetry-guided individualization of radiopharmaceutical therapy. Our goal was to build regression models using pretherapy 68Ga-DOTATATE PET uptake data and other baseline clinical factors/biomarkers to predict renal absorbed dose delivered by 177Lu-DOTATATE peptide receptor radionuclide therapy (177Lu-PRRT) for neuroendocrine tumors. We explore the combination of biomarkers and 68Ga PET uptake metrics, hypothesizing that they will improve predictive power over univariable regression. PATIENTS AND METHODS: Pretherapy 68Ga-DOTATATE PET/CTs were analyzed for 25 patients (50 kidneys) who also underwent quantitative 177Lu SPECT/CT imaging at approximately 4, 24, 96, and 168 hours after cycle 1 of 177Lu-PRRT. Kidneys were contoured on the CT of the PET/CT and SPECT/CT using validated deep learning-based tools. Dosimetry was performed by coupling the multi-time point SPECT/CT images with an in-house Monte Carlo code. Pretherapy renal PET SUV metrics, activity concentration per injected activity (Bq/mL/MBq), and other baseline clinical factors/biomarkers were investigated as predictors of the 177Lu SPECT/CT-derived mean absorbed dose per injected activity to the kidneys using univariable and bivariable models. Leave-one-out cross-validation (LOOCV) was used to estimate model performance using root mean squared error and absolute percent error in predicted renal absorbed dose including mean absolute percent error (MAPE) and associated standard deviation (SD). RESULTS: The median therapy-delivered renal dose was 0.5 Gy/GBq (range, 0.2-1.0 Gy/GBq). In LOOCV of univariable models, PET uptake (Bq/mL/MBq) performs best with MAPE of 18.0% (SD = 13.3%), and estimated glomerular filtration rate (eGFR) gives an MAPE of 28.5% (SD = 19.2%). Bivariable regression with both PET uptake and eGFR gives LOOCV MAPE of 17.3% (SD = 11.8%), indicating minimal improvement over univariable models. CONCLUSIONS: Pretherapy 68Ga-DOTATATE PET renal uptake can be used to predict post-177Lu-PRRT SPECT-derived mean absorbed dose to the kidneys with accuracy within 18%, on average. Compared with PET uptake alone, including eGFR in the same model to account for patient-specific kinetics did not improve predictive power. Following further validation of these preliminary findings in an independent cohort, predictions using renal PET uptake can be used in the clinic for patient selection and individualization of treatment before initiating the first cycle of PRRT.

A Pipeline for Automated Voxel Dosimetry: Application in Patients with Multi-SPECT/CT Imaging After 177Lu-Peptide Receptor Radionuclide Therapy.

Patient-specific dosimetry in radiopharmaceutical therapy (RPT) is impeded by the lack of tools that are accurate and practical for the clinic. Our aims were to construct and test an integrated voxel-level pipeline that automates key components (organ segmentation, registration, dose-rate estimation, and curve fitting) of the RPT dosimetry process and then to use it to report patient-specific dosimetry in 177Lu-DOTATATE therapy. Methods: An integrated workflow that automates the entire dosimetry process, except tumor segmentation, was constructed. First, convolutional neural networks (CNNs) are used to automatically segment organs on the CT portion of one post-therapy SPECT/CT scan. Second, local contour intensity-based SPECT--SPECT alignment results in volume-of-interest propagation to other time points. Third, dose rate is estimated by explicit Monte Carlo (MC) radiation transport using the fast, Dose Planning Method code. Fourth, the optimal function for dose-rate fitting is automatically selected for each voxel. When reporting mean dose, we apply partial-volume correction, and uncertainty is estimated by an empiric approach of perturbing segmentations. Results: The workflow was used with 4-time-point 177Lu SPECT/CT imaging data from 20 patients with 77 neuroendocrine tumors, segmented by a radiologist. CNN-defined kidneys resulted in high Dice values (0.91-0.94) and only small differences (2%-5%) in mean dose when compared with manual segmentation. Contour intensity-based registration led to visually enhanced alignment, and the voxel-level fitting had high R 2 values. Across patients, dosimetry results were highly variable; for example, the average of the mean absorbed dose (Gy/GBq) was 3.2 (range, 0.2-10.4) for lesions, 0.49 (range, 0.24-1.02) for left kidney, 0.54 (range, 0.31-1.07) for right kidney, and 0.51 (range, 0.27-1.04) for healthy liver. Patient results further demonstrated the high variability in the number of cycles needed to deliver hypothetical threshold absorbed doses of 23 Gy to kidney and 100 Gy to tumor. The uncertainty in mean dose, attributable to variability in segmentation, averaged 6% (range, 3%-17%) for organs and 10% (range, 3%-37%) for lesions. For a typical patient, the time for the entire process was about 25 min (∼2 min manual time) on a desktop computer, including time for CNN organ segmentation, coregistration, MC dosimetry, and voxel curve fitting. Conclusion: A pipeline integrating novel tools that are fast and automated provides the capacity for clinical translation of dosimetry-guided RPT.

Synthetic pulmonary perfusion images from 4DCT for functional avoidance using deep learning.

Purpose.To develop and evaluate the performance of a deep learning model to generate synthetic pulmonary perfusion images from clinical 4DCT images for patients undergoing radiotherapy for lung cancer.Methods. A clinical data set of 58 pre- and post-radiotherapy99mTc-labeled MAA-SPECT perfusion studies (32 patients) each with contemporaneous 4DCT studies was collected. Using the inhale and exhale phases of the 4DCT, a 3D-residual network was trained to create synthetic perfusion images utilizing the MAA-SPECT as ground truth. The training process was repeated for a 50-imaging study, five-fold validation with twenty model instances trained per fold. The highest performing model instance from each fold was selected for inference upon the eight-study test set. A manual lung segmentation was used to compute correlation metrics constrained to the voxels within the lungs. From the pre-treatment test cases (N = 5), 50th percentile contours of well-perfused lung were generated from both the clinical and synthetic perfusion images and the agreement was quantified.Results. Across the hold-out test set, our deep learning model predicted perfusion with a Spearman correlation coefficient of 0.70 (IQR: 0.61-0.76) and a Pearson correlation coefficient of 0.66 (IQR: 0.49-0.73). The agreement of the functional avoidance contour pairs was Dice of 0.803 (IQR: 0.750-0.810) and average surface distance of 5.92 mm (IQR: 5.68-7.55).Conclusion. We demonstrate that from 4DCT alone, a deep learning model can generate synthetic perfusion images with potential application in functional avoidance treatment planning.