A method for tumor dosimetry based on hybrid planar-SPECT/CT images and semiautomatic segmentation

Pareto optimization of SPECT acquisition and reconstruction settings for 177Lu activity quantification

BACKGROUND

The aim was to investigate the noise and bias properties of quantitative 177Lu-SPECT with respect to the number of projection angles, and the number of subsets and iterations in the OS-EM reconstruction, for different total acquisition times.

METHODS

Experimental SPECT acquisition of six spheres in a NEMA body phantom filled with 177Lu was performed, using medium-energy collimators and 120 projections with 180 s per projection. Bootstrapping was applied to generate data sets representing acquisitions with 20 to 120 projections for 10 min, 20 min, and 40 min, with 32 noise realizations per setting. Monte Carlo simulations were performed of 177Lu-DOTA-TATE in an anthropomorphic computer phantom with three tumours (2.8 mL to 40.0 mL). Projections representing 24 h and 168 h post administration were simulated, each with 32 noise realizations. Images were reconstructed using OS-EM with compensation for attenuation, scatter, and distance-dependent resolution. The number of subsets and iterations were varied within a constrained range of the product number of iterations × number of projections ≤ 2400 . Volumes-of-interest were defined following the physical size of the spheres and tumours, the mean activity-concentrations estimated, and the absolute mean relative error and coefficient of variation (CV) over noise realizations calculated. Pareto fronts were established by analysis of CV versus mean relative error.

RESULTS

Points at the Pareto fronts with low CV and high mean error resulted from using a low number of subsets, whilst points at the Pareto fronts associated with high CV but low mean error resulted from reconstructions with a high number of subsets. The number of projection angles had limited impact.

CONCLUSIONS

For accurate estimation of the 177Lu activity-concentration from SPECT images, the number of projection angles has limited importance, whilst the total acquisition time and the number of subsets and iterations are parameters of importance.

Relationship Between Absorbed Dose and Response in Neuroendocrine Tumors Treated with [177Lu]Lu-DOTATATE

Peptide receptor radionuclide therapy presents the possibility of tracing and quantifying the uptake of the drug in the body and performing dosimetry, potentially allowing individualization of treatment schemes. However, the details of how neuroendocrine tumors (NETs) respond to different absorbed doses are insufficiently known. Here, we investigated the relationship between tumor-absorbed dose and tumor response in a cohort of patients with NETs treated with [177Lu]Lu-DOTATATE. Methods: This was a retrospective study based on 69 tumors in 32 patients treated within a clinical trial. Dosimetry was performed at each cycle of [177Lu]Lu-DOTATATE, rendering 366 individual absorbed dose assessments. Hybrid planar-SPECT/CT imaging using [177Lu]Lu-DOTATATE was used, including quantitative SPECT reconstruction, voxel-based absorbed dose rate calculation, semiautomatic image segmentation, and partial-volume correction. Changes in tumor volume were used to determine tumor response. The volume for each tumor was manually delineated on consecutive CT scans, giving a total of 712 individual tumor volume assessments. Tumors were stratified according to grade. The relationship between absorbed dose and response was investigated using mixed-effects models and logistic regression. Tumors smaller than 4 cm3 were excluded. Results: In grade 2 NETs, a clear relationship between absorbed dose and volume reduction was observed. Our observations suggest a 90% probability of partial tumor response for an accumulated tumor-absorbed dose of at least 135 Gy. Conclusion: Our findings are in accordance with previous observations regarding the relationship between tumor shrinkage and absorbed dose. Moreover, our data suggest an absorbed dose threshold for partial response in grade 2 NETs. These observations provide valuable insights for the design of dosimetry-guided peptide receptor radionuclide therapy schemes.

Accuracy of patient-specific I-131 dosimetry using hybrid whole-body planar-SPECT/CT I-123 and I-131 imaging

PURPOSE

This study aimed to assess the accuracy of patient-specific absorbed dose calculations for tumours and organs at risk in radiopharmaceutical therapy planning, utilizing hybrid planar-SPECT/CT imaging.

METHODS

Three Monte Carlo (MC) simulated digital patient phantoms were created, with time-activity data for mIBG labelled to I-123 (LEHR and ME collimators) and I-131 (HE collimator). The study assessed the accuracy of the mean absorbed doses for I-131-mIBG therapy treatment planning. Multiple planar whole-body (WB) images were simulated (between 1 to 72 h post-injection (p.i)). The geometric-mean image of the anterior and posterior WB images was calculated, with scatter and attenuation corrections applied. Time-activity curves were created for regions of interest over the liver and two tumours (diameters: 3.0 cm and 5.0 cm) in the WB images. A corresponding SPECT study was simulated at 24 h p.i and reconstructed using the OS-EM algorithm, incorporating scatter, attenuation, collimator-detector response, septal scatter and penetration corrections. MC voxel-based absorbed dose rate calculations used two image sets, (i) the activity distribution represented by the SPECT images and (ii) the activity distribution from the SPECT images distributed uniformly within the volume of interest. Mean absorbed doses were calculated considering photon and charged particle emissions, and beta emissions only. True absorbed doses were calculated by MC voxel-based dosimetry of the known activity distributions for reference.

RESULTS

Considering photon and charged particle emissions, mean absorbed dose accuracies across all three radionuclide-collimator combinations of 3.8 ± 5.5% and 0.1 ± 0.9% (liver), 5.2 ± 10.0% and 4.3 ± 1.7% (3.0 cm tumour) and 15.0 ± 5.8% and 2.6 ± 0.6% (5.0 cm tumour) were obtained for image set (i) and (ii) respectively. Considering charged particle emissions, accuracies of 2.7 ± 4.1% and 5.7 ± 0.7% (liver), 3.2 ± 10.2% and 9.1 ± 1.7% (3.0 cm tumour) and 13.6 ± 5.7% and 7.0 ± 0.6% (5.0 cm tumour) were obtained for image set (i) and (ii) respectively.

CONCLUSION

The hybrid WB planar-SPECT/CT method proved accurate for I-131-mIBG dosimetry, suggesting its potential for personalized treatment planning.

Kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy based on multiple small VOIs

PURPOSE

The aim was to investigate the use of multiple small VOIs for kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy.

METHOD

The study was based on patient and simulated SPECT images in anthropomorphic geometries. Images were reconstructed using two reconstruction programs (local LundaDose and commercial Hermia) using OS-EM with and without resolution recovery (RR). Five small VOIs were placed to determine the average activity concentration (AC) in each kidney. The study consisted of three steps: (i) determination of the number of iterations for AC convergence based on simulated images; (ii) determination of recovery-coefficients (RCs) for 2 mL VOIs using a separate set of simulated images; (iii) assessment of operator variability in AC estimates for simulated and patient images. Five operators placed the VOIs, using for guidance: a) SPECT/CT with RR, b) SPECT/CT without RR, and c) CT only. For simulated images, time-integrated ACs (TIACs) were evaluated. For patient images, estimated ACs were compared with results of a previous method based on whole-kidney VOIs.

RESULTS

Eight iterations and ten subsets were sufficient for both programs and reconstruction settings. Mean RCs (mean ± SD) with RR were 1.03 ± 0.02 (LundaDose) and 1.10 ± 0.03 (Hermia), and without RR 0.91 ± 0.03 (LundaDose) and 0.94 ± 0.03 (Hermia). Most stable and accurate estimates of the AC were obtained using five 2-mL VOIs guided by SPECT/CT with RR, applying them to images without RR, and including an explicit RC for recovery correction.

CONCLUSION

The small VOI method based on five 2-mL VOIs was found efficient and sufficiently accurate for kidney dosimetry in [177Lu]Lu-DOTA-TATE therapy.

Future trends for patient-specific dosimetry methodology in molecular radiotherapy

Molecular radiotherapy is rapidly expanding, and new radiotherapeutics are emerging. The majority of treatments is still performed using empirical fixed activities and not tailored for individual patients. Molecular radiotherapy dosimetry is often seen as a promising candidate that would allow personalisation of treatments as outcome should ultimately depend on the absorbed doses delivered and not the activities administered. The field of molecular radiotherapy dosimetry has made considerable progress towards the feasibility of routine clinical dosimetry with reasonably accurate absorbed-dose estimates for a range of molecular radiotherapy dosimetry applications. A range of challenges remain with respect to the accurate quantification, assessment of time-integrated activity and absorbed dose estimation. In this review, we summarise a range of technological and methodological advancements, mainly focussed on beta-emitting molecular radiotherapeutics, that aim to improve molecular radiotherapy dosimetry to achieve accurate, reproducible, and streamlined dosimetry. We describe how these new technologies can potentially improve the often time-consuming considered process of dosimetry and provide suggestions as to what further developments might be required.

Tumor dosimetry using 177Lu: influence of background activity, measurement method and reconstruction algorithm

BACKGROUND

Image-based tumor dosimetry after radionuclide therapy, using the isotope ¹⁷⁷Lu, finds application e.g., for tumor-to-organ dose comparison and for dose response evaluation. When the tumor extent is not much larger than the image resolution, and when ¹⁷⁷Lu is found in nearby organs or other tumors, an accurate determination of tumor dose is particularly challenging. Here a quantitative evaluation of three different methods for determining the ¹⁷⁷Lu activity concentration in a phantom is performed, and the dependence on a variety of parameters is described. The phantom (NEMA IEC body phantom) has spheres of different size in a background volume, and sphere-to-background ¹⁷⁷Lu activity concentration ratios of infinity, 9.5, 5.0 and 2.7 are applied. The methods are simple to implement and well-known from the literature. They are based on (1) a large VOI encompassing the whole sphere, without background activity and with volume information from other sources, (2) a small VOI located in the sphere center, and (3) a VOI consisting of voxels with voxel value above a certain percentage of the maximum voxel value.

RESULTS

The determined activity concentration varies significantly with sphere size, sphere-to-background ratio, SPECT reconstruction method and method for determining the concentration. Based on the phantom study, criteria are identified under which the activity concentration can be determined with a maximal error of 40% even in the presence of background activity.

CONCLUSIONS

Tumor dosimetry is feasible in the presence of background activity using the above-mentioned methods, provided appropriate SPECT reconstructions are applied and tumors are selected for dosimetry analysis according to the following criteria for the three methods: (1) solitary tumor with diameter > 15 mm, (2) tumor diameter > 30 mm and tumor-to-background ratio > 2, and (3) tumor diameter > 30 mm and tumor-to-background ratio > 3.

Relationships between uptake of [68Ga]Ga-DOTA-TATE and absorbed dose in [177Lu]Lu-DOTA-TATE therapy

BACKGROUND

Somatostatin receptor ⁶⁸Ga PET imaging is standard for evaluation of a patient's suitability for ¹⁷⁷Lu peptide receptor radionuclide therapy of neuroendocrine tumours (NETs). The ⁶⁸Ga PET serves to ensure sufficient somatostatin receptor expression, commonly evaluated qualitatively. The aim of this study is to investigate the quantitative relationships between uptake in ⁶⁸Ga PET and absorbed doses in ¹⁷⁷Lu therapy.

METHOD

Eighteen patients underwent [⁶⁸Ga]Ga-DOTA-TATE PET imaging within 20 weeks prior to their first cycle of [¹⁷⁷Lu]Lu-DOTA-TATE. Absorbed doses for therapy were estimated for tumours, kidney, spleen, and normal liver parenchyma using a hybrid SPECT/CT-planar method. Gallium-68 activity concentrations were retrieved from PET images and also used to calculate SUVs and normalized SUVs, using blood and tissue for normalization. The ⁶⁸Ga activity concentrations per injected activity, SUVs, and normalized SUVs were compared with ¹⁷⁷Lu activity concentrations 1 d post-injection and ¹⁷⁷Lu absorbed doses. For tumours, for which there was a variable number per patient, both inter- and intra-patient correlations were analysed. Furthermore, the prediction of ¹⁷⁷Lu tumour absorbed doses based on a combination of tumour-specific ⁶⁸Ga activity concentrations and group-based estimates of the effective half-lives for grade 1 and 2 NETs was explored.

RESULTS

For normal organs, only spleen showed a significant correlation between the ⁶⁸Ga activity concentration and ¹⁷⁷Lu absorbed dose (r = 0.6). For tumours, significant, but moderate, correlations were obtained, with respect to both inter-patient (r = 0.7) and intra-patient (r = 0.45) analyses. The correlations to absorbed doses did not improve when using ⁶⁸Ga SUVs or normalized SUVs. The relationship between activity uptakes for ⁶⁸Ga PET and ¹⁷⁷Lu SPECT was stronger, with correlation coefficients r = 0.8 for both inter- and intra-patient analyses. The ¹⁷⁷Lu absorbed dose to tumour could be predicted from the ⁶⁸Ga activity concentrations with a 95% coverage interval of - 65% to 248%.

CONCLUSIONS

On a group level, a high uptake of [⁶⁸Ga]Ga-DOTA-TATE is associated with high absorbed doses at ¹⁷⁷Lu-DOTA-TATE therapy, but the relationship has a limited potential with respect to individual absorbed dose planning. Using SUV or SUV normalized to reference tissues do not improve correlations compared with using activity concentration per injected activity.

Milestones in dosimetry for nuclear medicine therapy

Nuclear Medicine therapy has reached a critical juncture with an unprecedented number of patients being treated and an extensive list of new radiopharmaceuticals under development. Since the early applications of these treatments dosimetry has played a vital role in their development, in both aiding optimisation and enhancing safety and efficacy. To inform the future direction of this field, it is useful to reflect on the scientific and technological advances that have occurred since those early uses. In this review, we explore how dosimetry has evolved over the years and discuss why such initiatives were conceived and the importance of maintaining standards within our practise. Specific milestones and landmark publications are highlighted and a thematic review and significant outcomes during each decade are presented.

EANM dosimetry committee recommendations for dosimetry of 177Lu-labelled somatostatin-receptor- and PSMA-targeting ligands

The purpose of the EANM Dosimetry Committee is to provide recommendations and guidance to scientists and clinicians on patient-specific dosimetry. Radiopharmaceuticals labelled with lutetium-177 (177Lu) are increasingly used for therapeutic applications, in particular for the treatment of metastatic neuroendocrine tumours using ligands for somatostatin receptors and prostate adenocarcinoma with small-molecule PSMA-targeting ligands. This paper provides an overview of reported dosimetry data for these therapies and summarises current knowledge about radiation-induced side effects on normal tissues and dose-effect relationships for tumours. Dosimetry methods and data are summarised for kidneys, bone marrow, salivary glands, lacrimal glands, pituitary glands, tumours, and the skin in case of radiopharmaceutical extravasation. Where applicable, taking into account the present status of the field and recent evidence in the literature, guidance is provided. The purpose of these recommendations is to encourage the practice of patient-specific dosimetry in therapy with 177Lu-labelled compounds. The proposed methods should be within the scope of centres offering therapy with 177Lu-labelled ligands for somatostatin receptors or small-molecule PSMA.

Personalized Dosimetry in Targeted Radiation Therapy: A Look to Methods, Tools and Critical Aspects

Targeted radiation therapy (TRT) is a strategy increasingly adopted for the treatment of different types of cancer. The urge for optimization, as stated by the European Council Directive (2013/59/EURATOM), requires the implementation of a personalized dosimetric approach, similar to what already happens in external beam radiation therapy (EBRT). The purpose of this paper is to provide a thorough introduction to the field of personalized dosimetry in TRT, explaining its rationale in the context of optimization and describing the currently available methodologies. After listing the main therapies currently employed, the clinical workflow for the absorbed dose calculation is described, based on works of the most experienced authors in the literature and recent guidelines. Moreover, the widespread software packages for internal dosimetry are presented and critical aspects discussed. Overall, a selection of the most important and recent articles about this topic is provided.

Quantitative SPECT/CT for dosimetry of peptide receptor radionuclide therapy

Neuroendocrine tumors (NETs) are uncommon malignancies of increasing incidence and prevalence. As these slow growing tumors usually overexpress somatostatin receptors (SSTRs), the use of ⁶⁸Ga-DOTA-peptides (gallium-68 chelated with dodecane tetra-acetic acid to somatostatin), which bind to the SSTRs, allows for PET based imaging and selection of patients for peptide receptor radionuclide therapy (PRRT). PRRT with radiolabeled somatostatin analogues such as ¹⁷⁷Lu-DOTATATE (lutetium-177-[DOTA,Tyr³]-octreotate), is mainly used for the treatment of metastatic or inoperable NETs. However, PRRT is generally administered at a fixed injected activity in order not to exceed dose limits in critical organs, which is suboptimal given the variability in radiopharmaceutical uptake among patients. Advances in SPECT (single photon emission computed tomography) imaging enable the absolute quantitative measure of the true radiopharmaceutical distribution providing for PRRT dosimetry in each patient. Personalized PRRT based on patient-specific dosimetry could improve therapeutic efficacy by optimizing effective tumor absorbed dose while limiting treatment related radiotoxicity.

Analysis of absorbed dose in radioimmunotherapy with 177Lu-trastuzumab using two different imaging scenarios: a pilot study

OBJECTIVES

Internal organ dosimetry is an important procedure to demonstrate the reliable application of 177Lu-trastuzumab radioimmunotherapy for human epidermal growth factor receptor-positive metastatic breast cancers. We are reporting the first human dosimetry study for 177Lu-trastuzumab. Another objective of our study was to calculate and compare the absorbed doses for normal organs and tumor lesions in patients before radioimmunotherapy with 177Lu-trastuzumab using two different imaging scenarios.

METHODS

Eleven patients (48.27 ± 8.95 years) with a history of metastatic breast cancer were included in the study. Postadministration of 177Lu-trastuzumab (351.09 ± 23.89 MBq/2 mg), acquisition was performed using planar and hybrid imaging scenarios at 4, 24, 72 and 168 h. Single-photon emission computed tomography/computed tomography imaging was performed at 72 h postinjection. Acquired images were processed using Dosimetry Toolkit software for the estimation of normalized cumulated activity in organs and tumor lesions. OLINDA/EXM 2.0 software was used for absorbed dose calculation in both scenarios.

RESULTS

Significant difference in normalized cumulated activity and the absorbed dose is noted between two imaging scenarios for the organs and tumor lesions (P < 0.05). Mean absorbed dose (mGy/MBq) estimated from heart, lungs, liver, spleen, kidney, adrenal, pancreas and colon using planar and hybrid scenarios were 0.81 ± 0.19 and 0.63 ± 0.17; 0.75 ± 0.13 and 0.32 ± 0.06; 1.26 ± 0.25 and 1.01 ± 0.17; 0.68 ± 0.22 and 0.53 ± 0.16; 0.91 ± 0.3 and 0.69 ± 0.24; 0.18 ± 0.04 and 0.11 ± 0.02; 0.25 ± 0.22 and 0.09 ± 0.02 and 0.75 ± 0.61 and 0.44 ± 0.28, respectively.

CONCLUSIONS

On the basis of our dosimetric evaluation, we concluded that radioimmunotherapy with 177Lu-trastuzumab is well tolerated to be implemented in routine clinical practice against HER2 positive metastatic breast cancer. Liver is the main critical organ at risk. Hybrid scenario demonstrated significantly lower absorbed doses in organs and tumors compared to the multiplanar method.

Dosimetric quantities of neuroendocrine tumors over treatment cycles with 177Lu-DOTA-TATE

Tumor dosimetry was performed for ¹⁷⁷Lu-DOTA-TATE with the aims of better understanding i) the range and variation of the tumor absorbed doses (ADs), ii) how different dosimetric quantities evolve over the treatment cycles, and iii) whether this evolution differs depending on the tumor grade. Such information is important for radiobiological interpretation and may inform the design of alternative administration schemes. Methods: Data come from 41 patients with neuroendocrine tumors (NETs) of grade 1 (n = 23) or 2 (n = 18), that had received between 2 and 9 treatment cycles. Dosimetry was performed for 182 individual lesions, giving in total 880 individual AD assessments across all cycles. Hybrid planar-SPECT/CT imaging was used, including quantitative SPECT reconstruction, voxel-based absorbed-dose-rate calculation, semi-automatic image segmentation, and partial-volume correction. Linear mixed-effect models were used to analyze changes over cycles in tumor ADs, absorbed-dose rates and activity concentrations at day-1, effective half-times, and tumor volumes. Tumors smaller than 8 ml were excluded from analyses. Results: Tumor ADs ranged between 2 and 77 Gy per cycle. On average the AD decreased over the cycles, with significantly different rates (P < 0.05) for grade 1 and 2 NETs of 6% and 14% per cycle, respectively. The absorbed-dose rates and activity concentrations at day-1 decreased by similar amounts. The effective half-times were less variable but shorter for grade 2 than grade 1 (P < 0.001). For grade 2 NETS the tumor volumes decreased, with a similar tendency in grade 1. Conclusion: The tumor AD, absorbed-dose rate and activity uptake decrease, in parallel with tumor volumes, between ¹⁷⁷Lu-DOTA-TATE treatment cycles, particularly for grade 2 NETs. The effective half-times vary less but are lower for grade 2 than grade 1 NETs. These results may indicate the development of radiation-induced fibrosis and could have implications for the design of future treatment and dosimetry protocols.

Modelling of energy-dependent spectral resolution for SPECT Monte Carlo simulations using SIMIND

PURPOSE

Monte Carlo (MC) modelling techniques have been used extensively in Nuclear Medicine (NM). The theoretical energy resolution relationship ( ∝ 1 / E ), does not accurately predict the gamma camera detector response across all energies. This study aimed to validate the accuracy of an energy resolution model for the SIMIND MC simulation code emulating the Siemens Symbia T16 dual-head gamma camera.

METHODS

Measured intrinsic energy resolution data (full width half maximum (FWHM) values), for Ba-133, Lu-177, Am-241, Ga-67, Tc-99m, I-123, I-131 and F-18 sources in air, were used to create a fitted model of the energy response of the gamma camera. Both the fitted and theoretical models were used to simulate intrinsic and extrinsic energy spectra using three different scenarios (source in air; source in simple scatter phantom and a clinical voxel-based digital patient phantom).

RESULTS

The results showed the theoretical model underestimated the FWHM values at energies above 160.0 keV up to 23.5 keV. In contrast, the fitted model better predicted the measured FWHM values with differences less than 3.3 keV. The I-131 in-scatter energy spectrum simulated with the fitted model better matched the measured energy spectrum. Higher energy photopeaks, (I-123: 528.9 keV and I-131: 636.9 keV) simulated with the fitted model, more accurately resembled the measured photopeaks. The voxel-based digital patient phantom energy spectra, simulated with the fitted and theoretical models, showed the potential impact of an incorrect energy resolution model when simulating isotopes with multiple photopeaks.

CONCLUSION

Modelling of energy resolution with the proposed fitted model enables the SIMIND user to accurately simulate NM images. A great improvement was seen for high-energy photon emitting isotopes (e.g. I-131), as well as isotopes with multiple photopeaks (e.g. Lu-177, I-131 and Ga-67) in comparison to the theoretical model. This will result in accurate evaluation of radioactivity quantification, which is vital for dosimetric purposes.

Influence of sampling schedules on [177Lu]Lu-PSMA dosimetry

BACKGROUND

Individualized dosimetry is recommended for [¹⁷⁷Lu]Lu-PSMA radioligand therapy (RLT) which is resource-intensive and protocols are often not optimized. Therefore, a simulation study was performed focusing on the determination of efficient optimal sampling schedules (OSS) for renal and tumour dosimetry by investigating different numbers of time points (TPs).

METHODS

Sampling schedules with 1-4 TPs were investigated. Time-activity curves of the kidneys and two tumour lesions were generated based on a physiologically based pharmacokinetic (PBPK) model and biokinetic data of 13 patients who have undergone [¹⁷⁷Lu]Lu-PSMA I&T therapy. Systematic and stochastic noise of different ratios was considered when modelling time-activity data sets. Time-integrated activity coefficients (TIACs) were estimated by simulating the hybrid planar/SPECT method for schedules comprising at least two TPs. TIACs based on one single SPECT/CT measurement were estimated using an approximation for reducing the number of fitted parameters. For each sampling schedule, the root-mean-squared error (RMSE) of the deviations of the simulated TIACs from the ground truths for 1000 replications was used as a measure for accuracy and precision.

RESULTS

All determined OSS included a late measurement at 192 h p.i., which was necessary for accurate and precise tumour TIACs. OSS with three TPs were identified to be 3-4, 96-100 and 192 h with an additional SPECT/CT measurement at the penultimate TP. Kidney and tumour RMSE of 6.4 to 7.7% and 6.3 to 7.8% were obtained, respectively. Shortening the total time for dosimetry to e.g. 96 h resulted in kidney and tumour RMSE of 6.8 to 8.3% and 9.1 to 11%, respectively. OSS with four TPs showed similar results as with three TPs. Planar images at 4 and 68 h and a SPECT/CT shortly after the 68 h measurement led to kidney and tumour RMSE of 8.4 to 12% and 12 to 16%, respectively. One single SPECT/CT measurement at 52 h yielded good approximations for the kidney TIACs (RMSE of 7.0%), but led to biased tumour TIACs.

CONCLUSION

OSS allow improvements in accuracy and precision of renal and tumour dosimetry for [¹⁷⁷Lu]Lu-PSMA therapy with potentially less effort. A late TP is important regarding accurate tumour TIACs.

Feasibility and limitations of quantitative SPECT for 223Ra

The aim of this paper is to investigate the feasibility and limitations of activity-concentration estimation for ²²³Ra using SPECT. Phantom measurements are performed using spheres (volumes 5.5 mL to 26.4 mL, concentrations 1.6 kBq mL^-1 to 4.5 kBq mL^-1). Furthermore, SPECT projections are simulated using the SIMIND Monte Carlo program for two geometries, one similar to the physical phantom and the other being an anthropomorphic phantom with added lesions (volumes 34 mL to 100 mL, concentrations 0.5 kBq mL^-1 to 4 kBq mL^-1). Medium-energy and high-energy collimators, 60 projections with 55 s per projection and a 20% energy window at 82 keV are employed. For the Monte Carlo simulated images, Poisson-distributed noise is added in ten noise realizations. Reconstruction is performed (OS-EM, 40 iterations, 6 subsets) employing compensation for attenuation, scatter, and collimator-detector response. The estimated concentrations in the anthropomorphic phantom are also corrected using recovery coefficients. Errors for the largest sphere in the physical phantom range from -25% to -34% for the medium-energy collimator and larger deviations for smaller spheres. Corresponding results for the high-energy collimator are -15% to -31%. The corresponding Monte Carlo simulations show standard deviations of a few percentage points. For the anthropomorphic phantom, before application of recovery coefficients the bias ranges from -16% to -46% (medium-energy collimator) and -10% to -28% (high-energy collimator), with standard deviations of 2% to 14% and 1% to 16%. After the application of recovery coefficients, the biases range from -3% to -35% (medium energy collimator) and from 0% to -18%. The errors decrease with increasing concentrations. Activity-concentration estimation of ²²³Ra with SPECT is feasible, but problems with repeatability need to be further studied.

Technical Note: Optimal sampling schedules for kidney dosimetry based on the hybrid planar/SPECT method in 177 Lu-PSMA therapy

PURPOSE

Accurate and precise renal dosimetry during ¹⁷⁷ Lu-labeled prostate-specific membrane antigen (PSMA) radioligand therapy is crucial for therapy decisions. Sampling schedules for estimating the necessary time-integrated activity coefficients (TIACs) are not optimized and standardized for clinical practice. Therefore, a simulation study to determine optimal sampling schedules (OSSs) was performed on 13 virtual ¹⁷⁷ Lu-PSMA I&T therapy patients.

METHOD

A total of 880 clinically feasible sampling schedules for planar imaging (three time points) were investigated. To simulate the hybrid planar/SPECT method, an additional quantitative SPECT/CT measurement following one planar image was considered. For each sampling schedule and patient, the activity values were generated separately. Measurement noise was modeled by drawing random numbers of log-normal distributions. The used fractional standard deviations (FSD) differed depending on the imaging modality. For activity values assigned to planar imaging, systematic noise between 25% and 75% of the total noise was simulated. After fitting with a mono-exponential function, the root-mean-squared errors of the deviations of the simulated TIACs from the ground truth for 1000 replications were used to determine the OSS. The uncertainties of the TIACs and renal dose coefficients were estimated.

RESULTS

For the hybrid planar/SPECT method, OSSs were determined to be (3-4, 72-76, 124-144)  h post injection (p.i.) with the quantitative SPECT/CT scan shortly after the second measurement. The accuracy and precision of the determined TIACs were in the range of (-3.0 ±  6.2)% and (-1.0  ±  6.5)%. This precision was improved by a factor 2-3 compared to dosimetry based on planar images only. Similar results were obtained for the renal dose coefficients. The virtual patients' renal dose coefficients were (0.68  ± 0.24)  Gy/GBq indicating that a population-based method yields an uncertainty of 35%.

CONCLUSIONS

Dosimetry based on the hybrid planar/SPECT method with OSS outperforms dosimetry based on planar images. The high variability in dose coefficients between the virtual patients demonstrates the need for individualized dosimetry.