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.

Comparison of a 3D CZT and conventional SPECT/CT system for quantitative Lu-177 SPECT imaging

PURPOSE

Next-generation SPECT/CT systems with CdZnTe (CZT) digital detectors in a ring-like setup are emerging to perform quantitative Lu-177 SPECT imaging in clinical routine. It is essential to assess how the shorter acquisition time might affect the image quality and uncertainty on the mean absorbed dose of the tumors and organs at risk compared to a conventional system.

METHODS

A NEMA Image Quality phantom was scanned with a 3D CZT SPECT/CT system (Veriton, by Spectrum Dynamics) using 6 min per bed position and with a conventional SPECT/CT system (Symbia T16, by Siemens) using 16 min per bed position. The sphere-to-background ratio was 12:1 and the background activity concentration ranged from 0.52 to 0.06 MBq/mL. A clinical reconstruction protocol for dosimetry purposes was determined for both systems by maximizing the sphere-to-background ratio while keeping the coefficient of variation of the background as low as possible. The corresponding image resolution was determined by the matching filter method and used for a dose uncertainty assessment of both systems following an established uncertainty model..

RESULTS

The optimized iterative reconstruction protocol included scatter and attenuation correction for both systems and detector response modeling for the Siemens system. For the 3D CZT system, 6 iterations and 8 subsets were combined with a Gaussian post-filter of 3 mm Full Width Half Maximum (FWHM) for post-smoothing. For the conventional system, 16 iterations and 16 subsets were applied with a Gaussian post-smoothing filter of 1 mm FWHM. For these protocols, the sphere-to-background ratio was 18.5% closer to the true ratio for the conventional system compared to the 3D CZT system when considering the four largest spheres. Meanwhile, the background coefficient of variation was very similar for both systems. These protocols resulted in SPECT image resolution of 14.8 mm and 13.6 mm for the 3D CZT and conventional system respectively. Based on these resolution estimates, a 50% dose uncertainty corresponded to a lesion volume of 28 mL for the conventional system and a lesion volume of 33 mL for the 3D CZT system.

CONCLUSIONS

An optimized reconstruction protocol for a Veriton system with 6 min of acquisition time per bed position resulted in slightly higher dose uncertainties than a conventional Symbia system using 16 min of acquisition time per bed position. Therefore, a 3D CZT SPECT/CT allows to significantly reduce the acquisition times with only a very limited impact on dose uncertainties such that quantitative Lu-177 SPECT/CT imaging becomes much more accessible for treatment concurrent dosimetry. Nevertheless, the uncertainty of SPECT-based dose estimates remains high.

3D printed non-uniform anthropomorphic phantoms for quantitative SPECT

BACKGROUND

A 3D printing grid-based method was developed to construct anthropomorphic phantoms with non-uniform activity distributions, to be used for evaluation of quantitative SPECT images. The aims were to characterize the grid-based method and to evaluate its capability to provide realistically shaped phantoms with non-uniform activity distributions.

METHODS

Characterization of the grid structures was performed by printing grid-filled spheres. Evaluation was performed by micro-CT imaging to investigate the printing accuracy and by studying the modulation contrast ([Formula: see text]) in SPECT images for 177Lu and 99mTc as a function of the grid fillable-volume fraction (FVF) determined from weighing. The grid-based technique was applied for the construction of two kidney phantoms and two thyroid phantoms, designed using templates from the XCAT digital phantoms. The kidneys were constructed with a hollow outer container shaped as cortex, an inner grid-based structure representing medulla and a solid section representing pelvis. The thyroids consisted of two lobes printed as grid-based structures, with void hot spots within the lobes. The phantoms were filled with solutions of 177Lu (kidneys) or 99mTc (thyroids) and imaged with SPECT. For verification, Monte Carlo simulations of SPECT imaging were performed for activity distributions corresponding to those of the printed phantoms. Measured and simulated SPECT images were compared qualitatively and quantitatively.

RESULTS

Micro-CT images showed that printing inaccuracies were mainly uniform across the grid. The relationships between the FVF from weighing and [Formula: see text] were found to be linear (r = 0.9995 and r = 0.9993 for 177Lu and 99mTc, respectively). The FVF-deviations from the design were up to 15% for thyroids and 4% for kidneys, mainly related to possibilities of cleaning after printing. Measured and simulated SPECT images of kidneys and thyroids exhibited similar activity distributions and quantitative comparisons agreed well, thus verifying the grid-based method.

CONCLUSIONS

We find the grid-based technique useful for the provision of 3D printed, realistically shaped, phantoms with non-uniform activity distributions, which can be used for evaluation of different quantitative methods in SPECT imaging.

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.

A list-mode multi-energy window low-count SPECT reconstruction method for isotopes with multiple emission peaks

BACKGROUND

Single-photon emission computed tomography (SPECT) provides a mechanism to perform absorbed-dose quantification tasks for [Formula: see text]-particle radiopharmaceutical therapies ([Formula: see text]-RPTs). However, quantitative SPECT for [Formula: see text]-RPT is challenging due to the low number of detected counts, the complex emission spectrum, and other image-degrading artifacts. Towards addressing these challenges, we propose a low-count quantitative SPECT reconstruction method for isotopes with multiple emission peaks.

METHODS

Given the low-count setting, it is important that the reconstruction method extracts the maximal possible information from each detected photon. Processing data over multiple energy windows and in list-mode (LM) format provide mechanisms to achieve that objective. Towards this goal, we propose a list-mode multi energy window (LM-MEW) ordered-subsets expectation-maximization-based SPECT reconstruction method that uses data from multiple energy windows in LM format and include the energy attribute of each detected photon. For computational efficiency, we developed a multi-GPU-based implementation of this method. The method was evaluated using 2-D SPECT simulation studies in a single-scatter setting conducted in the context of imaging [[Formula: see text]Ra]RaCl[Formula: see text], an FDA-approved RPT for metastatic prostate cancer.

RESULTS

The proposed method yielded improved performance on the task of estimating activity uptake within known regions of interest in comparison to approaches that use a single energy window or use binned data. The improved performance was observed in terms of both accuracy and precision and for different sizes of the region of interest.

CONCLUSIONS

Results of our studies show that the use of multiple energy windows and processing data in LM format with the proposed LM-MEW method led to improved quantification performance in low-count SPECT of isotopes with multiple emission peaks. These results motivate further development and validation of the LM-MEW method for such imaging applications, including for [Formula: see text]-RPT SPECT.

Quantitative SPECT (QSPECT) at high count rates with contemporary SPECT/CT systems

Background

Accurate QSPECT is crucial in dosimetry-based, personalized radiopharmaceutical therapy with ¹⁷⁷Lu and other radionuclides. We compared the quantitative performance of three NaI(Tl)-crystal SPECT/CT systems equipped with low-energy high-resolution collimators from two vendors (Siemens Symbia T6; GE Discovery 670 and NM/CT 870 DR).

Methods

Using up to 14 GBq of ^99mTc in planar mode, we determined the calibration factor and dead-time constant under the assumption that these systems have a paralyzable behaviour. We monitored their response when one or both detectors were activated. QSPECT capability was validated by SPECT/CT imaging of a customized NEMA phantom containing up to 17 GBq of ^99mTc. Acquisitions were reconstructed with a third-party ordered subset expectation maximization algorithm.

Results

The Siemens system had a higher calibration factor (100.0 cps/MBq) and a lower dead-time constant (0.49 μs) than those from GE (75.4–87.5 cps/MBq; 1.74 μs). Activities of up to 3.3 vs. 2.3–2.7 GBq, respectively, were quantifiable by QSPECT before the observed count rate plateaued or decreased. When used in single-detector mode, the QSPECT capability of the former system increased to 5.1 GBq, whereas that of the latter two systems remained independent of the detectors activation mode.

Conclusion

Despite similar hardware, SPECT/CT systems’ response can significantly differ at high count rate, which impacts their QSPECT capability in a post-therapeutic setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00421-3.

Optimization of SPECT/CT imaging protocols for quantitative and qualitative 99mTc SPECT

Background

The introduction of hybrid SPECT/CT devices enables quantitative imaging in SPECT, providing a methodological setup for quantitation using SPECT tracers comparable to PET/CT. We evaluated a specific quantitative reconstruction algorithm for SPECT data using a ^99mTc-filled NEMA phantom. Quantitative and qualitative image parameters were evaluated for different parametrizations of the acquisition and reconstruction protocol to identify an optimized quantitative protocol.

Results

The reconstructed activity concentration (AC_rec) and the signal-to-noise ratio (SNR) of all examined protocols (n = 16) were significantly affected by the parametrization of the weighting factor k used in scatter correction, the total number of iterations and the sphere volume (all, p < 0.0001). The two examined SPECT acquisition protocols (with 60 or 120 projections) had a minor impact on the AC_rec and no significant impact on the SNR.

In comparison to the known AC, the use of default scatter correction (k = 0.47) or object-specific scatter correction (k = 0.18) resulted in an underestimation of AC_rec in the largest sphere volume (26.5 ml) by − 13.9 kBq/ml (− 16.3%) and − 7.1 kBq/ml (− 8.4%), respectively. An increase in total iterations leads to an increase in estimated AC and a decrease in SNR. The mean difference between AC_rec and known AC decreased with an increasing number of total iterations (e.g., for 20 iterations (2 iterations/10 subsets) = − 14.6 kBq/ml (− 17.1%), 240 iterations (24i/10s) = − 8.0 kBq/ml (− 9.4%), p < 0.0001). In parallel, the mean SNR decreased significantly from 2i/10s to 24i/10s by 76% (p < 0.0001).

Conclusion

Quantitative SPECT imaging is feasible with the used reconstruction algorithm and hybrid SPECT/CT, and its consistent implementation in diagnostics may provide perspectives for quantification in routine clinical practice (e.g., assessment of bone metabolism). When combining quantitative analysis and diagnostic imaging, we recommend using two different reconstruction protocols with task-specific optimized setups (quantitative vs. qualitative reconstruction). Furthermore, individual scatter correction significantly improves both quantitative and qualitative results.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-021-00405-3.

Experimental validation of absolute SPECT/CT quantification for response monitoring in patients with coronary artery disease

BACKGROUND

Quantitative SPECT enables absolute quantification of uptake in perfusion defects. The aim of this experimental study is to assess quantitative accuracy and precision of a novel iterative reconstruction technique (Evolution; GE Healthcare) for the potential application of response monitoring using ^99mTc-tetrofosmin SPECT/CT in patients with coronary artery disease (CAD).

METHODS

Acquisitions of an anthropomorphic torso phantom with cardiac insert containing defects (with varying sizes), filled with ^99mTc-pertechnetate, were performed on a SPECT/CT (Discovery 670 Pro, GE Healthcare). Subsequently, volumes of interest of the defects were manually drawn on CT to assess the recovery coefficient (RC). Bull's eye plots were composed to evaluate the uptake per segment. Finally, ^99mTc-tetrofosmin SPECT/CT scans of 10 CAD patients were used to illustrate clinical application.

RESULTS

The phantom study indicated that Evolution showed convergence after 7 iterations and 10 subsets. The average repeatability deviation of all configurations was 2.91% and 3.15% (%SD mean) for filtered (Butterworth) and unfiltered data, respectively. The accuracy after post-filtering was lower compared to the unfiltered data with a mean (SD) RC of 0.63 (0.05) and 0.70 (0.07), respectively (p < 0.05). More artificial defects were found on Bull's eye plots created with the unfiltered data compared to filtered data. Eight out of ten patients showed significant changes in uptake before and after treatment (p < 0.05).

CONCLUSION

Quantification of ^99mTc-tetrofosmin SPECT/CT seems feasible for CAD patients when 7 iterations (10 subsets), Butterworth post-filtering (cut off frequency 0.52 in cycles/cm, order of 5) and manual CT-delineation are applied. However, future prospective patient studies are required for clinical application.

Impact of dead time on quantitative 177Lu-SPECT (QSPECT) and kidney dosimetry during PRRT

BACKGROUND

Dead time may affect the accuracy of quantitative SPECT (QPSECT), and thus of dosimetry. The aim of this study was to quantify the effect of dead time on ¹⁷⁷Lu-QSPECT and renal dosimetry following peptide receptor radionuclide therapy (PRRT) of neuroendocrine tumours.

METHODS

QSPECT/CT was performed on days 1 and 3 during 564 personalized ¹⁷⁷Lu-octreotate cycles in 166 patients. The dead-time data for each scanning time point was compiled. The impact of not correcting QSPECT for the dead time was assessed for the kidney dosimetry. This was also estimated for empiric PRRT by simulating in our cohort a regime of 7.4 GBq/cycle.

RESULTS

The probability to observe a larger dead time increased with the injected activity. A dead-time loss greater than 5% affected 14.4% and 5.7% of QSPECT scans performed at days 1 and 3, respectively. This resulted in renal absorbed dose estimates that would have been underestimated by more than 5% in 5.7% of cycles if no dead-time correction was applied, with a maximum underestimation of 22.1%. In the simulated empiric regime, this potential dose underestimation would have been limited to 6.2%.

CONCLUSION

Dead-time correction improves the accuracy of dosimetry in ¹⁷⁷Lu radionuclide therapy and is warranted in personalized PRRT.

Comprehensive SPECT/CT system characterization and calibration for 177Lu quantitative SPECT (QSPECT) with dead-time correction

Background

Personalization of ¹⁷⁷Lu-based radionuclide therapy requires implementation of dosimetry methods that are both accurate and practical enough for routine clinical use. Quantitative single-photon emission computed tomography/computed tomography (QSPECT/CT) is the preferred scanning modality to achieve this and necessitates characterizing the response of the camera, and calibrating it, over the full range of therapeutic activities and system capacity. Various methods to determine the camera calibration factor (CF) and the deadtime constant (τ) were investigated, with the aim to design a simple and robust protocol for quantitative ¹⁷⁷Lu imaging.

Methods

The SPECT/CT camera was equipped with a medium energy collimator. Multiple phantoms were used to reproduce various attenuation conditions: rod sources in air or water-equivalent media, as well as a Jaszczak phantom with inserts. Planar and tomographic images of a wide range of activities were acquired, with multiple energy windows for scatter correction (double or triple energy window technique) as well as count rate monitoring over a large spectrum of energy. Dead time was modelled using the paralysable model. CF and τ were deduced by curve fitting either separately in two steps (CF determined first using a subset of low-activity acquisitions, then τ determined using the full range of activity) or at once (both CF and τ determined using the full range of activity). Total or segmented activity in the SPECT field of view was computed. Finally, these methods were compared in terms of accuracy to recover the known activity, in particular when planar-derived parameters were applied to the SPECT data.

Results

The SPECT camera was shown to operate as expected on a finite count rate range (up to ~ 350 kcps over the entire energy spectrum). CF and τ from planar (sources in air) and SPECT segmented Jaszczak data yielded a very good agreement (CF < 1% and τ < 3%). Determining CF and τ from a single curve fit made dead-time-corrected images less prone to overestimating recovered activity. Using triple-energy window scatter correction while acquiring one or more additional energy window(s) to enable wide-spectrum count rate monitoring (i.e. ranging 55–250 or 18–680 keV) yielded the most consistent results across the various geometries. The final, planar-derived calibration parameters for our system were a CF of 9.36 ± 0.01 cps/MBq and a τ of 0.550 ± 0.003 μs. Using the latter, the activity in a Jaszczak phantom could be quantified by QSPECT with an accuracy of 0.02 ± 1.10%.

Conclusions

Serial planar acquisitions of sources in air using an activity range covering the full operational capacity of the SPECT/CT system, with multiple energy windows for wide-spectrum count rate monitoring, and followed by simultaneous determination of CF and τ using a single equation derived from the paralysable model, constitutes a practical method to enable accurate dead-time-corrected QSPECT imaging in a post-¹⁷⁷Lu radionuclide therapy setting.

Quantitative bone SPECT/CT: high specificity for identification of prostate cancer bone metastases

Purpose

Bone scintigraphy with ^99mTc-labeled diphosphonates can identify prostate cancer bone metastases with high sensitivity, but relatively low specificity, because benign conditions such as osteoarthritis can also trigger osteoblastic reactions. We aimed to investigate the diagnostic performance of ^99mTc-2,3-dicarboxy propane-1,1-diphosphonate (^99mTc-DPD) uptake quantification by single-photon emission computed tomography coupled with computed tomography (SPECT/CT) for distinguishing prostate cancer bone metastases from spinal and pelvic osteoarthritic lesions.

Methods

We retrospectively assessed 26 bone scans from 26 patients with known prostate cancer bone metastases and 13 control patients with benign spinal and pelvic osteoarthritic changes without known neoplastic disease. Quantitative SPECT/CT (xSPECT, Siemens Symbia Intevo, Erlangen, Germany) was performed and standardized uptake values (SUVs) were quantified with measurements of SUV_max and SUV_mean (g/mL) in all bone metastases for the prostate cancer group and in spinal and pelvic osteoarthritic changes for the control group. We used receiver operating characteristics (ROC) curves to determine the optimum SUV_max cutoff value to distinguish between bone metastases and benign spinal and pelvic lesions.

Results

In total, 264 prostate cancer bone metastases were analyzed, showing a mean SUV_max and SUV_mean of 34.6 ± 24.6 and 20.8 ± 14.7 g/mL, respectively. In 24 spinal and pelvic osteoarthritic lesions, mean SUV_max and SUV_mean were 14.2 ± 3.8 and 8.9 ± 2.2 g/mL, respectively. SUV_max and SUV_mean were both significantly different between the bone metastases and osteoarthritic groups (p ≤ 0.0001). Using a SUV_max cutoff of 19.5 g/mL for prostate cancer bone metastases in the spine and pelvis, sensitivity, specificity, positive and negative predictive values were 87, 92, 99 and 49%, respectively.

Conclusion

This study showed significant differences in quantitative ^99mTc-DPD uptake on bone SPECT/CT between prostate cancer bone metastases and spinal and pelvic osteoarthritic changes, with higher SUV_max and SUV_mean in metastases. Using a SUV_max cutoff of 19.5 g/mL, high specificity and positive predictive value for metastases identification in the spine and pelvis were found, thus increasing accuracy of bone scintigraphy.

ROC study and SUV threshold using quantitative multi-modal SPECT for bone imaging

Background

We investigated the clinical performance of a quantitative multi-modal SPECT/CT reconstruction platform for yielding radioactivity concentrations of bone imaging with ^99mTc-methylene diphosphonate (MDP) or ^99mTc-dicarboxypropane diphosphonate (DPD). The novel reconstruction incorporates CT-derived tissue information while preserving the delineation of tissue boundaries. We assessed image-based reader concordance and confidence, and determined lesion classification and SUV thresholds from ROC analysis.

Methods

Seventy-two cancer patients were scanned at three US and two German clinical sites, each contributing two experienced board-certified nuclear medicine physicians as readers. We compared four variants of the reconstructed data resulting from the Flash3D (F3D) and the xSPECT Bone™ (xB) iterative reconstruction methods and presented images to the readers with and without a fused CT, resulting in four combinations. We used an all-or-none approach for inclusion, compiling results only when a reader completed all reads in a subset. After the final read, we conducted a “surrogate truth” reading, presenting all data to each reader. For any remaining discordant lesions, we conducted a consensus read. We next undertook ROC analysis to determine SUV thresholds for differentiating benign and lesional uptake.

Results

On a five-point rating scale of image quality, xB was deemed better by almost two points in resolution and one point better in overall acceptance compared to F3D. The absolute agreement of the rendered decision between the nine readers was significantly higher with CT information either inside the reconstruction (xB, xBCT) or simply through image fusion (F3DCT): 0.70 (xBCT), 0.67 (F3DCT), 0.64 (xB), and 0.46 (F3D). The confidence level to characterize the lesion was significantly higher (3.03x w/o CT, 1.32x w/CT) for xB than for F3D. There was high correlation between xB and F3D scores for lesion detection and classification, but lesion detection confidence was 41% higher w/o CT, and 21% higher w/CT for xB compared to F3D. Without CT, xB had 6.6% higher sensitivity, 7.1% higher specificity, and 6.9% greater AUC compared to F3D, and similarly with CT-fusion. The overall SUV-criterion (SUV_c) of xB (12) exceeded that for xSPECT Quant™ (xQ; 9), an approach not using the tissue delineation of xB. SUV critical numbers depended on lesion volume and location. For non-joint lesions > 6 ml, the AUC for xQ and xB was 94%, with SUV_c > 9.28 (xQ) or > 9.68 (xB); for non-joint lesions ≤ 6 ml, AUCs were 81% (xQ) and 88% (xB), and SUV_c > 8.2 (xQ) or > 9.1 (xB). For joint lesions, the AUC was 80% (xQ) and 83% (xB), with SUV_c > 8.61 (xQ) or > 13.4 (xB).

Conclusion

The incorporation of high-resolution CT-based tissue delineation in SPECT reconstruction (xSPECT Bone) provides better resolution and detects smaller lesions (6 ml), and the CT component facilitates lesion characterization. Our approach increases confidence, concordance, and accuracy for readers with a wide range of experience. The xB method retained high reading accuracy, despite the unfamiliar image presentation, having greatest impact for smaller lesions, and better localization of foci relative to bone anatomy. The quantitative assessment yielded an SUV-threshold for sensitively distinguishing benign and malignant lesions. Ongoing efforts shall establish clinically usable protocols and SUV thresholds for decision-making based on quantitative SPECT.

Does quantification have a role to play in the future of bone SPECT?

Routinely, there is a visual basis to nuclear medicine reporting: a reporter subjectively places a patient's condition into one of multiple discrete classes based on what they see. The addition of a quantitative result, such as a standardised uptake value (SUV), would provide a numerical insight into the nature of uptake, delivering greater objectivity, and perhaps improved patient management.For bone scintigraphy in particular quantification could increase the accuracy of diagnosis by helping to differentiate normal from abnormal uptake. Access to quantitative data might also enhance our ability to characterise lesions, stratify and monitor patients' conditions, and perform reliable dosimetry for radionuclide therapies. But is there enough evidence to suggest that we, as a community, should be making more effort to implement quantitative bone SPECT in routine clinical practice?We carried out multiple queries through the PubMed search engine to facilitate a cross-sectional review of the current status of bone SPECT quantification. Highly cited papers were assessed in more focus to scrutinise their conclusions.An increasing number of authors are reporting findings in terms of metrics such as SUVmax. Although interest in the field in general remains high, the rate of clinical implementation of quantitative bone SPECT remains slow and there is a significant amount of validation required before we get carried away.

Comparison of standardized uptake values between 99mTc-HDP SPECT/CT and 18F-NaF PET/CT in bone metastases of breast and prostate cancer

Background

Despite recent technological advances allowing for quantitative single-photon emission computed tomography (SPECT), quantitative SPECT has not been widely used in the clinical practice. The aim of this study is to evaluate the feasibility of quantitative SPECT for measuring metastatic bone uptake in breast and prostate cancer by comparing standard uptake values (SUVs) measured with ^99mTc-HDP SPECT/CT and ¹⁸F-NaF PET/CT.

Methods

Twenty-six breast and 27 prostate cancer patients at high risk of bone metastases underwent both ^99mTc-HDP SPECT/CT and ¹⁸F-NaF PET/CT within 14 days of each other. The SPECT and PET data were reconstructed using ordered-subset expectation-maximization algorithms achieving quantitative images. Metastatic and benign skeletal lesions visible in both data sets were identified, and their maximum, peak, and mean SUVs (SUV_max, SUV_peak, and SUV_mean) were determined. SUV ratios (SUVRs) between the lesions and adjacent normal appearing bone were also calculated. Linear regression was used to evaluate the correlations between the SUVs of SPECT and PET and Bland-Altman plots to evaluate the differences between the SUVs and SUVRs of SPECT and PET.

Results

A total of 231 skeletal lesions, 129 metastatic and 102 benign, were analyzed. All three SUV measures correlated very strongly between SPECT and PET (R² ≥ 0.80, p < 0.001) when all lesions were included, and the PET SUVs were significantly higher than SPECT SUVs (p < 0.001). The median differences were 21%, 12%, and 19% for SUV_max, SUV_peak, and SUV_mean, respectively. On the other hand, the SUVRs were similar between SPECT and PET with median differences of 2%, − 9%, and 2% for SUVR_max, SUVR_peak, and SUVR_mean, respectively.

Conclusion

The strong correlation between SUVs and similar SUVRs of ^99mTc-HDP SPECT/CT and ¹⁸F-NaF PET/CT demonstrate that SPECT is an applicable tool for clinical quantification of bone metabolism in osseous metastases in breast and prostate cancer patients.

Biodistribution, pharmacokinetics, and organ-level dosimetry for 188Re-AHDD-Lipiodol radioembolization based on quantitative post-treatment SPECT/CT scans

Background

Rhenium-188-labelled-Lipiodol radioembolization is a safe and cost-effective treatment for primary liver cancer. In order to determine correlations between treatment doses and patient response to therapy, accurate patient-specific dosimetry is required. Up to date, the reported dosimetry of ¹⁸⁸Re-Lipiodol has been based on whole-body (WB) planar imaging only, which has limited quantitative accuracy. The aim of the present study is to determine the in vivo pharmacokinetics, bio-distribution, and organ-level dosimetry of ¹⁸⁸Re-AHDD-Lipiodol radioembolization using a combination of post-treatment planar and quantitative SPECT/CT images. Furthermore, based on the analysis of the pharmacokinetic data, a practical and relatively simple imaging and dosimetry method that could be implemented in clinics for ¹⁸⁸Re-AHDD-Lipiodol radioembolization is proposed.

Thirteen patients with histologically proven hepatocellular carcinoma underwent ¹⁸⁸Re-AHDD-Lipiodol radioembolization. A series of 2–3 WB planar images and one SPECT/CT scan were acquired over 48 h after the treatment. The time-integrated activity coefficients (TIACs, also known as residence-times) and absorbed doses of tumors and organs at risk (OARs) were determined using a hybrid WB/SPECT imaging method.

Results

Whole-body imaging showed that ¹⁸⁸Re-AHDD-Lipiodol accumulated mostly in the tumor and liver tissue but a non-negligible amount of the pharmaceutical was also observed in the stomach, lungs, salivary glands, spleen, kidneys, and urinary bladder. On average, the measured effective half-life of ¹⁸⁸Re-AHDD-Lipiodol was 12.5 ± 1.9 h in tumor. The effective half-life in the liver and lungs (the two organs at risk) was 12.6 ± 1.7 h and 12.0 ± 1.9 h, respectively. The presence of ¹⁸⁸Re in other organs was probably due to the chemical separation and subsequent release of the free radionuclide from Lipiodol.

The average doses per injected activity in the tumor, liver, and lungs were 23.5 ± 40.8 mGy/MBq, 2.12 ± 1.78 mGy/MBq, and 0.11 ± 0.05 mGy/MBq, respectively. The proposed imaging and dosimetry method, consisting of a single SPECT/CT for activity determination followed by ¹⁸⁸Re-AHDD-Lipiodol clearance with the liver effective half-life of 12.6 h, resulted in TIACs estimates (and hence, doses) mostly within ± 20% from the reference TIACs (estimated using three WB images and one SPECT/CT).

Conclusions

The large inter-patient variability of the absorbed doses in tumors and normal tissue in ¹⁸⁸Re-HDD-Lipiodol radioembolization patients emphasizes the importance of patient-specific dosimetry calculations based on quantitative post-treatment SPECT/CT imaging.

Electronic supplementary material

The online version of this article (10.1186/s40658-018-0227-6) contains supplementary material, which is available to authorized users.

Accuracy and reproducibility of simplified QSPECT dosimetry for personalized 177Lu-octreotate PRRT

Background

Routine dosimetry is essential for personalized ¹⁷⁷Lu-octreotate peptide receptor radionuclide therapy (PRRT) of neuroendocrine tumors (NETs), but practical and robust dosimetry methods are needed for wide clinical adoption. The aim of this study was to assess the accuracy and inter-observer reproducibility of simplified dosimetry protocols based on quantitative single-photon emission computed tomography (QSPECT) with a limited number of scanning time points. We also updated our personalized injected activity (IA) prescription scheme.

Methods

Seventy-nine NET patients receiving ¹⁷⁷Lu-octreotate therapy (with a total of 279 therapy cycles) were included in our study. Three-time-point (3TP; days 0, 1, and 3) QSPECT scanning was performed following each therapy administration. Dosimetry was obtained using small volumes of interest activity concentration sampling for the kidney, the bone marrow and the tumor having the most intense uptake. Accuracy of the simplified dosimetry based on two-time-point (2TP; days 1 and 3, monoexponential fit) or a single-time-point (1TP_D3; day 3) scanning was assessed, as well as that of hybrid methods based on 2TP for the first cycle and 1TP (day 1 or 3; 2TP/1TP_D1 and 2TP/1TP_D3, respectively) or no imaging at all (based on IA only; 2TP/no imaging (NI)) for the subsequent induction cycles. The inter-observer agreement was evaluated for the 3TP, 2TP, and hybrid 2TP/1TP_D3 methods using a subset of 60 induction cycles (15 patients). The estimated glomerular filtration rate (eGFR), body size descriptors (weight, body surface area (BSA), lean body weight (LBW)), and products of both were assessed for their ability to predict IA per renal absorbed dose at the first cycle.

Results

The 2TP dosimetry estimates correlated highly with those from the 3TP data for all tissues (Spearman r > 0.99, P < 0.0001) with small relative errors between the methods, particularly for the kidney and the tumor, with median relative errors not exceeding 2% and interdecile ranges spanning over less than 6% and 4%, respectively, for the per-cycle and cumulative estimates. For the bone marrow, the errors were slightly greater (median errors < 6%, interdecile ranges < 14%). Overall, the strength of correlations of the absorbed dose estimates from the simplified methods with those from the 3TP scans tended to progressively decrease, and the relative errors to increase, in the following order: 2TP, 2TP/1TP_D3, 1TP_D3, 2TP/1TP_D1, and 2TP/NI. For the tumor, the 2TP/NI scenario was highly inaccurate due to the interference of the therapeutic response. There was an excellent inter-observer agreement between the three observers, in particular for the renal absorbed dose estimated using the 3TP and 2TP methods, with mean errors lesser than 1% and standard deviations of 5% or lower. The eGFR · LBW and eGFR · BSA products best predicted the ratio of IA to the renal dose (GBq/Gy) for the first cycle (Spearman r = 0.41 and 0.39, respectively; P < 0.001). For the first cycle, the personalized IA proportional to eGFR · LBW or eGFR · BSA decreased the range of delivered renal absorbed dose between patients as compared with the fixed IA. For the subsequent cycles, the optimal personalized IA could be determined based on the prior cycle renal GBq/Gy with an error of less than 21% in 90% of patients.

Conclusions

A simplified dosimetry protocol based on two-time-point QSPECT scanning on days 1 and 3 post-PRRT provides reproducible and more accurate dose estimates than the techniques relying on a single time point for non-initial or all cycles and results in limited patient inconvenience as compared to protocols involving scanning at later time points. Renal absorbed dose over the 4-cycle induction PRRT course can be standardized by personalizing IA based on the product of eGFR with LBW or BSA for the first cycle and on prior renal dosimetry for the subsequent cycles.

Tc-HDP quantitative SPECT/CT in transthyretin cardiac amyloid and the development of a reference interval for myocardial uptake in the non-affected population

Background

99mTechnetium-HDP (HDP) bone scans differentiate transthyretin (ATTR) cardiac amyloid from other infiltrative myocardial diseases. These scans are not quantitative and are assessed by comparing myocardial uptake to bone. This study examined whether quantitative HDP SPECT/CT can discriminate individuals with cardiac ATTR from the population without this disease.

Methods

HDP thoracic xSPECT/CT QUANT (xQUANT) was performed in 29 patients: ATTR cardiac amyloid (n = 6); AL cardiac amyloid (n = 1); other infiltrative myocardial disease (n = 4); no known infiltrative cardiac disease (n = 18). SUVmax measured within volumes of interest for whole heart, ascending aorta blood pool, and specific bones. HDP myocardial uptake calculated as whole heart minus blood pool.

Results

The cardiac ATTR group had greater HDP myocardial uptake than those with no known infiltrative disease (p = 0.002). AL and other myocardial diseases had uptake indistinguishable from the group with no known infiltrative cardiac disease. The SUVmaxima were sufficiently similar between individuals without cardiac ATTR that a 99% reference interval for HDP uptake could be calculated, providing an upper limit cut point of SUVmax 1.2. Individuals with cardiac ATTR had SUVmax well above this cut point.

Conclusion

Quantitative SPECT/CT can measure HDP myocardial uptake in individuals with normal hearts and those with cardiac ATTR without recourse to comparison with bone. It enables calculation of a reference interval for HDP myocardial uptake in the population without ATTR cardiac amyloid. Using this reference interval single individuals with cardiac ATTR can be accurately discriminated from the non-affected population. This technique uses a NIST traceable calibration source, potentially allowing development of multicentre clinical decision limits. Its role in disease management warrants further assessment.

Phantom and clinical evaluation of the effect of full Monte Carlo collimator modelling in post-SIRT yttrium-90 Bremsstrahlung SPECT imaging

Background

Post-therapy SPECT/CT imaging of ⁹⁰Y microspheres delivered to hepatic malignancies is difficult, owing to the continuous, high-energy Bremsstrahlung spectrum emitted by ⁹⁰Y. This study aimed to evaluate the utility of a commercially available software package (HybridRecon, Hermes Medical Solutions AB) which incorporates full Monte Carlo collimator modelling. Analysis of image quality was performed on both phantom and clinical images in order to ultimately provide a recommendation of an optimum reconstruction for post-therapy ⁹⁰Y microsphere SPECT/CT imaging.

A 3D-printed anthropomorphic liver phantom was filled with ⁹⁰Y with a sphere-to-background ratio of 4:1 and imaged on a GE Discovery 670 SPECT/CT camera. Datasets were reconstructed using ordered-subsets expectation maximization (OSEM) 1–7 iterations in order to identify the optimal OSEM reconstruction (5 iterations, 15 subsets). Quantitative analysis was subsequently carried out on phantom datasets obtained using four reconstruction algorithms: the default OSEM protocol (2 iterations, 10 subsets) and the optimised OSEM protocol, both with and without full Monte Carlo collimator modelling. The quantitative metrics contrast recovery (CR) and background variability (BV) were calculated.

The four algorithms were then used to retrospectively reconstruct 10 selective internal radiation therapy (SIRT) patient datasets which were subsequently blind scored for image quality by a consultant radiologist.

Results

The optimised OSEM reconstruction (5 iterations, 15 subsets with full MC collimator modelling) increased the CR by 42% (p < 0.001) compared to the default OSEM protocol (2 iterations, 10 subsets). The use of full Monte Carlo collimator modelling was shown to further improve CR by 14% (30 mm sphere, CR = 90%, p < 0.05).

The consultant radiologist had a significant preference for the optimised OSEM over the default OSEM protocol (p < 0.001), with the optimised OSEM being the favoured reconstruction in every one of the 10 clinical cases presented.

Conclusions

OSEM (5 iterations, 15 subsets) with full Monte Carlo collimator modelling is quantitatively the optimal image reconstruction for post-SIRT ⁹⁰Y Bremsstrahlung SPECT/CT imaging. The use of full Monte Carlo collimator modelling for correction of image-degrading effects significantly increases contrast recovery without degrading clinical image quality.

Quantitative 3D scintigraphy shows increased muscular uptake of pyrophosphate in idiopathic inflammatory myopathy

BACKGROUND

Nuclear imaging is increasingly being used in the diagnostic work-up of idiopathic inflammatory myopathy (IIM). Increased muscular uptake of technetium-99m-pyrophosphate (^99mTc-PYP) has hitherto been assessed qualitatively by planar scintigraphy. We set out to perform quantitative tomographic scintigraphy in IIM.

RESULTS

Ninety IIM patients and 48 control subjects underwent ^99mTc-PYP single-photon emission computed tomography (SPECT)/CT of the upper and lower body. Scans were evaluated visually by an intensity score (1-4) and quantitatively by the mean standardized uptake value (SUV_mean) in thigh muscles after semi-automated segmentation of these. Furthermore, a SUV_mean gradient down along the thighs was determined by linear regression of the slice-by-slice activity. Interobserver analyses were performed on qualitative evaluations. Compared to controls, patients more often had a high intensity score (p < 0.0001), but interobserver analyses revealed only moderate agreement. The thigh muscular ^99mTc-PYP activity (SUV_mean) was 60% higher in patients than in controls, p < 0.0001, albeit with a wide range. There was an activity gradient down the thigh muscle, the proximal tracer uptake being highest, and this gradient was steeper in patients than in controls; the activity decreased by 0.00024 and 0.00010 SUV_mean mm^-1, respectively, along the thighs.

CONCLUSIONS

The muscular uptake of ^99mTc-PYP was significantly higher in patients than in healthy controls by qualitative and quantitative assessment. The tracer uptake was higher in the proximal than in the distal part of the thigh muscle, and SUV_mean gradients differed between groups. Hence, tomographic nuclear imaging allowing for quantification of the ^99mTc-PYP uptake might contribute to the diagnosis of IIM, and SPECT/CT of the lower body might suffice.

Evaluation of sequential SPECT and CT for targeted radionuclide therapy dosimetry

PURPOSE

In targeted radionuclide therapy (TRT), a prior knowledge of the absorbed dose biodistribution is essential for pre-therapy treatment planning. Previously, we showed that non-rigid organ-by-organ registration in sequential quantitative SPECT images improved dose estimation. This study aims to investigate if sequential CT can further improve TRT dosimetric accuracy.

METHODS

We simulated SPECT/CT acquisitions at 1, 12, 24, 72 and 144 h In-111 Zevalin post-injection using an analytical MEGP projector, modeling attenuation, scatter and collimator-detector response. We later recruited a patient injected with 222 MBq In-111 DTPAOC imaged at 3 SPECT/CT sessions for clinical evaluations. Four registration schemes were evaluated: whole-body-based registration performed on sequential (1) SPECT (WB-SPECT) or (2) CT (WB-CT) images; organ-based registration applied on organs individually segmented from sequential (3) SPECT (O-SPECT) or (4) CT (O-CT) images. Voxel-by-voxel integration was performed followed by Y-90 voxel-S-kernel convolution. Organ-absorbed doses, iso-dose curves, dose-volume histograms (DVHs) were generated for targeted organs for analysis.

RESULTS

In simulation study, organ-absorbed dose errors were (- 8.66 ± 2.83)%, (- 2.51 ± 3.69)%, (- 9.23 ± 3.28)%, (- 7.17 ± 2.53)% for liver, (- 14.81 ± 4.91)%, (- 3.60 ± 4.37)%, (- 18.13 ± 4.44)%, (- 11.34 ± 4.22)% for spleen, for O-SPECT, O-CT, WB-SPECT and WB-CT registrations, respectively. For all organs, O-CT showed superior results. Results of iso-dose contour, DVHs were in accordance with the organ-absorbed doses. In clinical studies, the results were also consistent which showed O-CT method deviated the most from the result with no registration.

CONCLUSIONS

We conclude that if both sequential SPECT/CT scans are available, CT organ-based registration method can more effectively improve the 3D dose estimation. Sequential low-dose CT scans might be considered to be included in the standard TRT protocol.

Validation of a calibration method using the cross-calibration factor and system planar sensitivity in quantitative single-photon emission computed tomography imaging

The present study aimed to validate the absolute quantitative accuracy of a calibration method for single-photon emission computed tomography (SPECT) using cross-calibration factor (CCF)- and system sensitivity-based calibration methods. The CCF obtained with different reconstruction parameters was evaluated using a cylindrical phantom (diameter 20 cm, height 20 cm). SPECT images were acquired with a positron emission tomography/computed tomography (CT) phantom. Subsequently, they were reconstructed by using ordered subset expectation maximization with resolution recovery, scatter, and CT-based attenuation correction. All reconstructed SPECT counts were converted to activity concentrations based on the CCF and system planar sensitivity. We placed 12 circular regions of interest, 37 mm in diameter, on the phantom background, and the converted activity concentration and relative measurement error were assessed. The CCF obtained using a cylindrical phantom was affected by the iterative update number and post-smoothing filter function. The activity concentration calibrated using the CCF showed over- and underestimation. However, the activity concentration obtained from the system planar sensitivity was similar to that gained using the phantom. The values obtained using the system planar sensitivity were within 10% of the activity concentrations obtained with the phantom. These findings demonstrated that the calibration method using system planar sensitivity provides accurate quantification within 10% of the true activity concentration. Further clinical examination is required to validate the present results.

Objective evaluation of reconstruction methods for quantitative SPECT imaging in the absence of ground truth

Quantitative single-photon emission computed tomography (SPECT) imaging is emerging as an important tool in clinical studies and biomedical research. There is thus a need for optimization and evaluation of systems and algorithms that are being developed for quantitative SPECT imaging. An appropriate objective method to evaluate these systems is by comparing their performance in the end task that is required in quantitative SPECT imaging, such as estimating the mean activity concentration in a volume of interest (VOI) in a patient image. This objective evaluation can be performed if the true value of the estimated parameter is known, i.e. we have a gold standard. However, very rarely is this gold standard known in human studies. Thus, no-gold-standard techniques to optimize and evaluate systems and algorithms in the absence of gold standard are required. In this work, we developed a no-gold-standard technique to objectively evaluate reconstruction methods used in quantitative SPECT when the parameter to be estimated is the mean activity concentration in a VOI. We studied the performance of the technique with realistic simulated image data generated from an object database consisting of five phantom anatomies with all possible combinations of five sets of organ uptakes, where each anatomy consisted of eight different organ VOIs. Results indicate that the method provided accurate ranking of the reconstruction methods. We also demonstrated the application of consistency checks to test the no-gold-standard output.