Accuracy of quantitative reconstructions in SPECT/CT imaging

SPECT/CT imaging: quantifying 99mTc-MDP concentration in the spine and pelvis

OBJECTIVE

This study aimed to identify a relatively robust SUV for guiding clinical practice through quantitative measurement and comparison of various normalization methods based on the SUV of 99mTc-MDP in the normal spine and pelvis using an integrated SPECT/CT scanner.

METHODS

Between June 2017 and September 2019, a total of 500 oncology patients (mean age, 60.9; men, 66.0%) who underwent bone SPECT/CT scans with 99mTc-MDP were enrolled in this retrospective study. The mean SUV (SUVmean) of 4962 spinal and pelvic bones was calculated based on the patients' body weight (BW), lean body mass (LBM), bone mineral content (BMC), body surface area (BSA), and body mass index (BMI), defined as SUVbw, SUVlbm, SUVbmc, SUVbsa, and SUVbmi, respectively. The coefficients of variation (CoVs) of the aforementioned parameters were compared, and the correlation and multiple linear regression analyses were used to compare the extent to which these parameters were affected by sex, age, height, weight, BMI, and CT values.

RESULTS

The average SUVs in the normal spine and pelvis displayed a relatively wide variability: 4.573 ± 1.972 for SUVbw, 3.555 ± 1.517 for SUVlbm, 0.163 ± 0.071 for SUVbmc, 0.124 ± 0.052 for SUVbsa, and 1.668 ± 0.732 for SUVbmi. In general, SUVbsa had relatively lowest CoV (42.1%) in all vertebrae and pelvis compared with other SUVs. For correlation analyses, all SUVs displayed weak but significant correlations with age and CT values. For regression analyses, SUVbsa was influenced only by age, BMI, and CT values independently. The effects of these variables on SUVbsa were all smaller than those on conventional SUVbw.

CONCLUSIONS

The SUVs of 99mTc-MDP in normal bone derived from quantitative bone SPECT/CT could serve as a reference for evaluating tumor bone metastasis, but it should be assessed on a site-specific basis. SUVbsa exhibited superior robustness among all the SUV normalization variations, indicating potential clinical applications.

Quantifying partial volume effect in SPECT and planar imaging: optimizing region of interest for activity concentration estimation in different sphere sizes

INTRODUCTION

To quantify the partial volume effect in single photon emission tomography (SPECT) and planar images of Carlson phantom as well as providing an optimum region of interest (ROI) required to more accurately estimate the activity concentration for different sphere sizes.

METHODS

131 I solution with the 161.16 kBq/ml concentration was uniformly filled into the different spheres of Carlson phantom (cold background condition) with the diameters of 7.3, 9.2, 11.4, 14.3, 17.9, 22.4 and 29.9 mm, and there was no background activity. In the hot background condition, the spheres were filled with the solution of 131 I with the 1276.5 kBq/ml addition to the background activity concentration of 161.16 kBq/ml in all the phantoms. The spheres were mounted inside the phantom and underwent SPECT and planar images. ROI was drawn closely on the boundary of each sphere image and it was extended to extract the true count.

RESULTS

In the cold background condition, the recovery coefficient (RC) value for SPECT images ranged between 0.8 and 1.03. However, in planar imaging, the RC value was 0.72 for the smallest sphere size and it increased for larger spheres until 0.98 for 29.9 mm. In the hot background condition, the RC value for sphere diameters larger than 20 mm was overestimated more than in the cold background condition. The ROI/size required to more accurately determine activity concentration for the cold background ranged from 1.18 to 2.7. However, in the hot background condition, this ratio varied from 1.34 to 4.05.

CONCLUSION

In the quantification of partial volume effects, the spill-out effect seems to play a crucial role in the distribution of the image counts beyond the boundaries of the image pixels. However, more investigations are needed to accurately characterize limitations regarding the object size, background levels, and other factors.

Quantification and dosimetry of small volumes including associated uncertainty estimation

BACKGROUND

Accurate quantification of radioactivity in a source of interest relies on accurate registration between SPECT and anatomical images, and appropriate correction of partial volume effects (PVEs). For small volumes, exact registration between the two imaging modalities and recovery factors used to correct for PVE are unreliable. There is currently no guidance relating to quantification or the associated uncertainty estimation for small volumes.

MATERIAL AND METHODS

A method for quantification of small sources of interest is proposed, which uses multiple oversized volumes of interest. The method was applied to three Na[131I]I activity distributions where a Na[131I]I capsule was situated within a cylindrical phantom containing either zero background, uniform background or non-uniform background and to a scenario with small lesions placed in an anthropomorphic phantom. The Na[131I]I capsule and lesions were quantified using the proposed method and compared with measurements made using two alternative quantification methods. The proposed method was also applied to assess the absorbed dose delivered to a bone metastasis following [131I]mIBG therapy for neuroblastoma including the associated uncertainty estimation.

RESULTS

The method is accurate across a range of activities and in varied radioactivity distributions. Median percentage errors using the proposed method in no background, uniform backgrounds and non-uniform backgrounds were - 0.4%, - 0.3% and 1.7% with median associated uncertainties of 1.4%, 1.4% and 1.6%, respectively. The technique is more accurate and robust when compared to currently available alternative methods.

CONCLUSIONS

The proposed method provides a reliable and accurate method for quantification of sources of interest, which are less than three times the spatial resolution of the imaging system. The method may be of use in absorbed dose calculation in cases of bone metastasis, lung metastasis or thyroid remnants.

Reference values of mandibular condylar growth activity: a study of SUVmax with quantitative bone SPECT/CT

OBJECTIVE

To assess the condylar growth activity (CGA) with quantitative bone single photon emission computed tomography/computed tomography (SPECT/CT) to establish reference values of maximum standardized uptake value (SUVmax) and cutoff values for identifying active unilateral condylar hyperplasia (UCH) in different ages.

METHODS

We analyzed the CGA of 58 UCH patients and that of 125 volunteers as a control group by SUVmax of quantitative bone SPECT/CT imaging. The SUVmax and the uptake difference between bilateral condyles among different age groups were analyzed. SUVmax cutoff values for detecting active condyle were calculated by receiver operating characteristic curve analysis.

RESULTS

The condylar SUVmax in 10-19, 20-29, 30-39, 40-49 and 50-59 years old groups of volunteers were 6.24 ± 1.39, 4.76 ± 0.98, 3.23 ± 0.64, 3.00 ± 0.61 and 2.90 ± 0.53, respectively. The uptake difference between bilateral condyles in the control group was 3.84% ± 1.71%. The affected condylar SUVmax was significantly higher than that of the contralateral condyle in active UCH patients (6.03 ± 2.85 vs. 3.96 ± 1.07; Z = -5.264; P = 0.000). SUVmax of the affected condyles in active UCH patients was not statistically higher than condylar SUVmax in the corresponding age group (6.03 ± 2.85 vs. 5.50 ± 1.41; Z = -0.173; P = 0.863). SUVmax of the unaffected condyles was significantly lower than condylar SUVmax in the corresponding age group (3.96 ± 1.07 vs. 5.50 ± 1.41; Z = -5.833; P = 0.000). SUVmax cutoff values for identifying active condyle were 6.26 and 4.53 in patients of 13-19 and 20-29 years old, respectively.

CONCLUSIONS

The condylar SUVmax varied with age. Different cutoff values of condylar SUVmax should be employed for diagnosing active UCH for patients in different ages.

Quantitative evaluation of single-photon emission computed tomography findings in lower extremity possible without computed tomography-based attenuation correction

OBJECTIVE

For performing accurate quantitative analysis of single-photon emission computed tomography (SPECT)/computed tomography (CT) images, CT-based attenuation correction (CTAC) is considered to be necessary. However, the effect on quantitative values for an examined area close to the body surface, such as in the lower extremity, has yet to be elucidated. We performed the present investigation to determine the possibility of quantitative evaluation using a SPECT standalone device without CT.

METHODS

Validation was performed using clinical data of patients who underwent a lower extremity SPECT/CT examination, with grouping based on presence or absence of CTAC, scatter correction and resolution recovery. Using a reference group in which all types of correction were applied, standardized uptake values (SUVs), including maximum (SUVmax) and peak (SUVpeak), were examined in each group and compared.

RESULTS

As compared to the reference group, the difference in quantitative values became smaller in the order of the applied scatter correction and resolution recovery, applied resolution recovery, applied scatter correction, and neither scatter correction or resolution recovery applied groups, with no significant difference between the reference group and that with neither scatter correction or resolution recovery applied. A similar tendency was seen for both SUVmax and SUVpeak.

CONCLUSIONS

In bone SPECT quantitative examinations of the lower extremity, quantitative evaluation without CTAC is possible without the use of scatter correction or resolution recovery. Thus, quantitative evaluation can be performed with use of a standalone SPECT device without CT.

Evaluation of Iodine-123 and Iodine-131 SPECT activity quantification: a Monte Carlo study

Purpose

The quantitative accuracy of Nuclear Medicine images, acquired for both planar and SPECT studies, is influenced by the isotope-collimator combination as well as image corrections incorporated in the iterative reconstruction process. These factors can be investigated and optimised using Monte Carlo simulations. This study aimed to evaluate SPECT quantification accuracy for ¹²³I with both the low-energy high resolution (LEHR) and medium-energy (ME) collimators and ¹³¹I with the high-energy (HE) collimator.

Methods

Simulated SPECT projection images were reconstructed using the OS-EM iterative algorithm, which was optimised for the number of updates, with appropriate corrections for scatter, attenuation and collimator detector response (CDR), including septal scatter and penetration compensation. An appropriate calibration factor (CF) was determined from four different source geometries (activity-filled: water-filled cylindrical phantom, sphere in water-filled (cold) cylindrical phantom, sphere in air and point-like source), investigated with different volume of interest (VOI) diameters. Recovery curves were constructed from recovery coefficients to correct for partial volume effects (PVEs). The quantitative method was evaluated for spheres in voxel-based digital cylindrical and patient phantoms.

Results

The optimal number of OS-EM updates was 60 for all isotope-collimator combinations. The CF_point with a VOI diameter equal to the physical size plus a 3.0-cm margin was selected, for all isotope-collimator geometries. The spheres’ quantification errors in the voxel-based digital cylindrical and patient phantoms were less than 3.2% and 5.4%, respectively, for all isotope-collimator combinations.

Conclusion

The study showed that quantification errors of less than 6.0% could be attained, for all isotope-collimator combinations, if corrections for; scatter, attenuation, CDR (including septal scatter and penetration) and PVEs are performed. ¹²³I LEHR and ¹²³I ME quantification accuracies compared well when appropriate corrections for septal scatter and penetration were applied. This can be useful in departments that perform ¹²³I studies and may not have access to ME collimators.

Comparing calculated and experimental activity and dose values obtained from image-based quantification of 90Y SPECT/CT Data

PURPOSE

Selective internal radiation therapy (SIRT) is a treatment for various kinds of liver tumours by injecting ⁹⁰Y bearing microspheres into the liver vessels. To perform meaningful post-treatment dosimetry, quantitative imaging is performed.

METHODS

This work uses a Monte-Carlo based reconstruction software with scatter and attenuation correction and collimator modelling that allows the quantification of ⁹⁰Y bremsstrahlung SPECT/CT data. A dataset comprising 17 patients and measurements on a Jaszczak phantom, a NEMA IEC Body phantom and an anthropomorphic liver phantom are analysed and activities and dose values are acquired. These measured values are compared with applied activities and pre-treatment calculations, allowing to assess the quality of the SPECT reconstruction. A detailed uncertainty budget is presented, including uncertainties of the dose calibrator, the count rate, non-included interactions and other factors.

RESULTS

The applied method is validated by finding measurements repeatable within the given uncertainty, and it is shown the influence of various parameters on the reconstruction process is negligible. Furthermore, activities and doses measured in the phantoms show good agreement with calculated values, if they are corrected for partial volume effects.

CONCLUSIONS

The strict observation of metrological requirements and the creation of an uncertainty budget increase the reliability and traceability of this novel approach to ⁹⁰Y dosimetry. It gives an example of successful voxel-based dosimetry based on quantitative ⁹⁰Y SPECT/CT image data.

The effect of calibration factors and recovery coefficients on 177Lu SPECT activity quantification accuracy: a Monte Carlo study

BACKGROUND

Different gamma camera calibration factor (CF) geometries have been proposed to convert SPECT data into units of activity concentration. However, no consensus has been reached on a standardised geometry. The CF is dependent on the selected geometry and is further affected by partial volume effects. This study investigated the effect of two CF geometries and their corresponding recovery coefficients (RCs) on the quantification accuracy of 177Lu SPECT images using Monte Carlo simulations.

METHODS

The CF geometries investigated were (i) a radioactive-sphere surrounded by non-radioactive water (sphere-CF) and (ii) a cylindrical phantom uniformly filled with radioactive water (cylinder-CF). Recovery coefficients were obtained using the sphere-CF and cylinder-CF, yielding the sphere-RC and cylinder-RC values, respectively, for partial volume correction (PVC). The quantification accuracy was evaluated using four different-sized spheres (15.6-65.4 ml) and a kidney model with known activity concentrations inside a cylindrical, torso and patient phantom. Images were reconstructed with the 3D OS-EM algorithm incorporating attenuation, scatter and detector-response corrections. Segmentation was performed using the physical size and a small cylindrical volume inside the cylinder for the sphere-CF and cylinder-CF, respectively.

RESULTS

The sphere quantification error (without PVC) was better for the sphere-CF (≤ - 5.54%) compared to the cylinder-CF (≤ - 20.90%), attributed to the similar geometry of the quantified and CF spheres. Partial volume correction yielded comparable results for the sphere-CF-RC (≤ 3.47%) and cylinder-CF-RC (≤ 3.53%). The accuracy of the kidney quantification was poorer (≤ 22.34%) for the sphere-CF without PVC compared to the cylinder-CF (≤ 2.44%). With PVC, the kidney quantification results improved and compared well for the sphere-CF-RC (≤ 3.50%) and the cylinder-CF-RC (≤ 3.45%).

CONCLUSION

The study demonstrated that upon careful selection of CF-RC combinations, comparable quantification errors (≤ 3.53%) were obtained between the sphere-CF-RC and cylinder-CF-RC, when all corrections were applied.

Validation of a Monte Carlo modelled gamma camera for Lutetium-177 imaging

This study validated a model of the Siemens Symbia T16 dual-head SPECT/CT gamma camera created using the Monte Carlo program SIMIND for ¹⁷⁷Lu. The validation was done by comparing experimental and simulated gamma camera performance criteria tests for the ¹⁷⁷Lu 208 keV photopeak with a medium-energy collimator. Results showed good agreement between the experimental and simulated values. These results illustrated that SIMIND could emulate the Symbia T16 successfully and therefore, can be used with confidence to model ¹⁷⁷Lu images.

Quantitative SPECT/CT-Technique and Clinical Applications

The continuous development of SPECT over the past 50 years has led to improved image quality and increased diagnostic confidence. The most influential developments include the realization of hybrid SPECT/CT devices, as well as the implementation of attenuation correction and iterative image reconstruction techniques. These developments have led to a preference for SPECT/CT devices over SPECT-only systems and to the widespread adoption of the former, strengthening the role of SPECT/CT as the workhorse of Nuclear Medicine imaging. New trends in the ongoing development of SPECT/CT are diverse. For example, whole-body SPECT/CT images, consisting of acquisitions from multiple consecutive bed positions in the manner of PET/CT, are increasingly performed. Additionally, in recent years, some interesting approaches in detector technology have found their way into commercial products. For example, some SPECT cameras dedicated to specific organs employ semiconductor detectors made of cadmium telluride or cadmium zinc telluride, which have been shown to increase the obtainable image quality by offering a higher sensitivity and energy resolution. However, the advent of quantitative SPECT/CT which, like PET, can quantify the amount of tracer in terms of Bq/mL or as a standardized uptake value could be regarded as most important development. It is a major innovation that will lead to increased diagnostic accuracy and confidence, especially in longitudinal studies and in the monitoring of treatment response. The current work comprises two main aspects. At first, physical and technical fundamentals of SPECT image formation are described and necessary prerequisites of quantitative SPECT/CT are reviewed. Additionally, the typically achievable quantitative accuracy based on reports from the literature is given. Second, an extensive list of studies reporting on clinical applications of quantitative SPECT/CT is provided and reviewed.

Analysis of quantitative [I-123] mIBG SPECT/CT in a phantom and in patients with neuroblastoma

Purpose

To determine the accuracy of quantitative SPECT, intersystem and interpatient standardized uptake value (SUV) calculation consistency for a manufacturer-independent quantitative SPECT/CT reconstruction algorithm, and the range of SUVs of normal and neoplastic tissue.

Methods

A NEMA body phantom with 6 spheres (ranging 10–37 mm) was filled with a known activity-to-volume ratio and used to determine the contrast recovery coefficient (CRC) for each visible sphere, and the measured SUV accuracy of those spheres and background water solution. One hundred eleven 123I-metaiodobenzylguanidine ([I-123] mIBG) SPECT/CT examinations from 43 patients were reconstructed using SUV SPECT® (HERMES Medical Solutions Inc.); 42 examinations were acquired using a GE Infinia Hawkeye 4 SPECT/CT, and 69 were acquired on a Siemens Symbia Intevo SPECT/CT. Inter scanner SUV analysis of 9 regions of normal [I-123] mIBG tissue uptake was conducted. Intrapatient mean SUV variability was calculated by measuring normal liver uptake within patients scanned on both cameras. The intensity of uptake by neoplastic tissue in the images was quantified using maximum SUV and, if present, compared over time.

Results

The phantom results of the visible spheres and background resulted in accuracy calculations better than 5–10% with CRC correction. Interscanner SUV variability showed no statistical difference (average p value 0.559; range 0.066–1.0) among the 9 normal tissues analyzed. Intrapatient liver mean SUV varied ≤ 16% as calculated for 28 patients (87 examinations) studied on both scanners. In one patient, a thoracic tumor evaluated over 10 time points (18 months) underwent a 74% (3.1/12.0) reduction in maximum SUV with treatment.

Conclusion

The results demonstrate quantitative accuracy to better than 10%, and both consistent SUV calculation between 2 different SPECT/CT scanners for 9 tissues, and low intrapatient measurement variability for quantitative SPECT/CT analysis in a pediatric population with neuroblastoma. Quantitative SPECT/CT offers the opportunity for objective analysis of tumor response using [I-123] mIBG by normalizing the uptake to injected dose and patient weight, as is done for PET.

Quantification evaluation of 99mTc-MDP concentration in the lumbar spine with SPECT/CT: compare with bone mineral density

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 of osteoporosis. The aim of this study is to evaluate the feasibility of quantitative bone SPECT/CT for measuring lumbar standard uptake value (SUV) in patients with different bone-mineral density (BMD), and investigate the correlation between SUV measured with ^99mTc-methylene diphosphonate (MDP) SPECT/CT and BMD assessment by dual-energy X-ray absorptiometry (DXA).

METHODS

A retrospective analysis of 62 cases ^99mTc-MDP whole-body bone imaging and local lumbar SPECT/CT tomography were performed. According to the results of dual-energy X-ray bone density examination, they were divided into normal group, osteopenic group, and osteoporosis group. The raw SPECT data were reconstructed using flash3D which includes attenuation correction, scatter compensation, and collimator resolution recovery, SPECT images from this algorithm were calibrated for SUV analysis. Comparing difference of lumbar SUV in different BMD subjects, and investigating the correlation between lumbar SUV and BMD. Data were analyzed by one-way ANOVA and Pearson regression analysis using SPSS 17.0 software.

RESULTS

The maximum SUV (SUV_max) and mean SUV (SUV_mean) of L1-L4 vertebral in 62 subjects were 7.39 ± 1.84 and 4.90 ± 1.27, respectively. The average BMD was 0.85 ± 0.15 (g/cm²), and the average CT value was 145.88 ± 53.99 (HU). The SUV_max, SUV_mean, BMD, and CT values of the lumbar spine were statistically significantly different among the three groups (F = 24.089, 30.501, 94.847, 30.241, all p < 0.001), and the osteopenic group was significantly lower than the normal group (all p < 0.001), the osteoporosis group was significantly lower than the normal group and the osteopenic group (all p < 0.001). Lumbar SUV_max, SUV_mean, and BMD were significantly negatively correlated with age (r = - 0.328 to - 0.442, all p < 0.05), and positively correlated with body weight and CT value (r = 0.299-0.737, all p < 0.05), but no significant correlation with height (r = 0.006-0.175, all p > 0.05). Lumbar SUV_max and SUV_mean increased significantly with the increase of BMD (r = 0.638, 0.632, p < 0.001).

CONCLUSION

The SUV of lumbar spine in ^99mTc-MDP bone SPECT/CT was significantly different among subjects with different BMD, and the SUV was positively correlated with BMD. These findings justify that quantitative bone SPECT/CT is an applicable tool for clinical quantification of bone metabolism in osteoporosis patients.

A simulation-based method to determine optimal sampling schedules for dosimetry in radioligand therapy

AIM

For dosimetry in radioligand therapy, the time-integrated activity coefficients (TIACs) for organs at risk and for tumour lesions have to be determined. The used sampling scheme affects the TIACs and therefore the calculated absorbed doses. The aim of this work was to develop a general and flexible method, which analyses numerous clinically applicable sampling schedules using true time-activity curves (TACs) of virtual patients.

METHOD

Nine virtual patients with true TACs of the tumours were created using a physiologically-based pharmacokinetic (PBPK) model and individual biokinetic data of five patients with neuroendocrine tumours and four with meningioma. ¹¹¹In-DOTATATE was used for pre-therapeutic dosimetry. In total, 15,120 sampling schemes, each consisting of 4 time points, were investigated. Gaussian noise of different levels was added to the corresponding true time-activity points. A bi-exponential function was used to fit the simulated time-activity data. For each investigated sampling schedule, 1000 replications were performed. Patient-specific and population-specific optimal sampling schedules were determined using the relative root-mean-square error (rRMSE). Furthermore, the fractions of TIACs a˜ deviating >5% (f_Δa˜>5%) and >10% (f_Δa˜>10%) from the true TIAC a˜_true were used for additional evaluations e.g. to investigate the effect of varying single time points.

RESULTS

Almost all patient-specific and all population-specific optimal sampling schedules have t₄≥96h for all noise levels. Changing the latest time point from the population-specific optimal time to e.g. 48h leads to a median increase of f_Δa˜>10% from 0.1% to 88% for the lowest investigated noise level. Using the determined population-specific optimal schedules, results in more accurate and precise results than established schedules from the literature.

CONCLUSION

A method of determining the optimal sampling schedule for dosimetry, which considers clinical working hours and measurement uncertainties, has been developed and applied. The simulation study shows that optimised sampling schedules result in high accuracy and precision of the determined TIACs.

Influence of the SPECT calibration source position on the absorbed dose calculation for 131I-NaI therapy using GATE simulations

Many research groups have studied nuclear medicine image quantification to improve its accuracy in dose estimation. This work aims to evaluate the influence of the source calibration position for absorbed dose calculation for a ¹³¹I-NaI therapy using Monte Carlo (MC) simulations. The calibration approach consisted of a cylindrical phantom filled with water. A cylindrical ¹³¹I source with 361.1 ± 3.6 kBq ml^-1 was positioned at the center of the phantom and its outer part. Images were acquired with 150 00 counts per projection image acquired with SPECT detector (high counts density-HCD) and 3000 counts per projection (low counts density-LCD). MC simulations, performed with GATE code, were validated by comparing the S values of a water sphere uniformly filled with ¹³¹I, as from the sphere model of OLINDA/EXM 1.1. Calibration factors deviation between central and peripheral calibrations is more significant for HCD (18.3%) than for LCD images (3.7%). The 3D dose distribution map obtained from GATE resulted in a dose factor equal to 1.5 × 10^-3 mGy/(MBq.s). For both HCD and LCD images, the commonly used approach, which employs the central source calibration to obtain the dose from a peripheral source, resulted in dose overestimation. Results suggest that organ dose calculation can be improved considering the organ position in the field of view. Finally, patients' radiation protection in dosimetry studies could be improved considering the calibration source position, due to the superior accuracy in dose calculation.

A new pharmacokinetic model for 90Y-ibritumomab tiuxetan based on 3-dimensional dosimetry

Monoclonal antibodies (mAbs) are key components in several therapies for cancer and inflammatory diseases but current knowledge of their clinical pharmacokinetics and distribution in human tissues remains incomplete. Consequently, optimal dosing and scheduling in clinics are affected. With sequential radiolabeled mAb-based imaging, radiation dosing in tissues/organs can be calculated to provide a better assessment of mAb concentrations in tissues. This is the first pharmacokinetic model of ⁹⁰Y-Ibritumomab tiuxetan (⁹⁰Y-IT) in humans to be described, based on three-dimensional (3D) dosimetry using single-photon emission computed-tomography coupled with computed-tomography. 19 patients with follicular lymphoma were treated initially with ⁹⁰Y-IT in the FIZZ trial. Based on a compartmental approach individualising the vascular compartment within studied organs, this study proposes a reliable pharmacokinetic (PK) five-compartment model replacing the currently used two-compartment model and constitutes a new direction for further research. This model provides exchange constants between the different tissues, Area Under the Curve of ¹¹¹In-IT in blood (AUC) and Mean Residence Time (MRT) that have not been reported so far for IT. Finally, the elimination process appears to occur in a compartment other than the liver or the spleen and suggests the metabolism of mAbs may take place mainly on the vascular compartment level.

Accuracy of 177Lu activity quantification in SPECT imaging: a phantom study

BACKGROUND

The aim of the study is to assess accuracy of activity quantification of 177Lu studies performed according to recommendations provided by the committee on Medical Internal Radiation Dose (MIRD) pamphlets 23 and 26. The performances of two scatter correction and three segmentation methods were compared. Additionally, the accuracy of tomographic and planar methods for determination of the camera normalization factor (CNF) was evaluated. Eight phantoms containing inserts of different sizes and shapes placed in air, water, and radioactive background were scanned using a Siemens SymbiaT SPECT/CT camera. Planar and tomographic scans with 177Lu sources were used to measure CNF. Images were reconstructed with our SPEQToR software using resolution recovery, attenuation, and two scatter correction methods (analytical photon distribution interpolated (APDI) and triple energy window (TEW)). Segmentation was performed using a fixed threshold method for both air and cold water scans. For hot water experiments three segmentation methods were compared as folows: a 40% fixed threshold, segmentation based on CT images, and our iterative adaptive dual thresholding (IADT). Quantification error, defined as the percent difference between experimental and true activities, was evaluated.

RESULTS

Quantification error for scans in air was better for TEW scatter correction (<6%) than for APDI (<11%). This trend was reversed for scans in water (<10% for APDI and <14% for TEW). For hot water, the best results (<18% for small objects and <5% for objects >100 ml) were obtained when APDI and IADT were used for scatter correction and segmentation, respectively. Additionally, we showed that planar acquisitions with scatter correction and tomographic scans provide similar CNF values. This is an important finding because planar acquisitions are easier to perform than tomographic scans. TEW and APDI resulted in similar quantification errors with APDI showing a small advantage for objects placed in medium with non-uniform density.

CONCLUSIONS

Following the MIRD recommendations for data acquisition and reconstruction resulted in accurate activity quantification (errors <5% for large objects). However, techniques for better organ/tumor segmentation must still be developed.

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.

Quantitative Imaging for Targeted Radionuclide Therapy Dosimetry - Technical Review

Targeted radionuclide therapy (TRT) is a promising technique for cancer therapy. However, in order to deliver the required dose to the tumor, minimize potential toxicity in normal organs, as well as monitor therapeutic effects, it is important to assess the individualized internal dosimetry based on patient-specific data. Advanced imaging techniques, especially radionuclide imaging, can be used to determine the spatial distribution of administered tracers for calculating the organ-absorbed dose. While planar scintigraphy is still the mainstream imaging method, SPECT, PET and bremsstrahlung imaging have promising properties to improve accuracy in quantification. This article reviews the basic principles of TRT and discusses the latest development in radionuclide imaging techniques for different theranostic agents, with emphasis on their potential to improve personalized TRT dosimetry.

Accuracy of Rhenium-188 SPECT/CT activity quantification for applications in radionuclide therapy using clinical reconstruction methods

The main applications of ¹⁸⁸Re in radionuclide therapies include trans-arterial liver radioembolization and palliation of painful bone-metastases. In order to optimize ¹⁸⁸Re therapies, the accurate determination of radiation dose delivered to tumors and organs at risk is required. Single photon emission computed tomography (SPECT) can be used to perform such dosimetry calculations. However, the accuracy of dosimetry estimates strongly depends on the accuracy of activity quantification in ¹⁸⁸Re images. In this study, we performed a series of phantom experiments aiming to investigate the accuracy of activity quantification for ¹⁸⁸Re SPECT using high-energy and medium-energy collimators. Objects of different shapes and sizes were scanned in Air, non-radioactive water (Cold-water) and water with activity (Hot-water). The ordered subset expectation maximization algorithm with clinically available corrections (CT-based attenuation, triple-energy window (TEW) scatter and resolution recovery was used). For high activities, the dead-time corrections were applied. The accuracy of activity quantification was evaluated using the ratio of the reconstructed activity in each object to this object's true activity. Each object's activity was determined with three segmentation methods: a 1% fixed threshold (for cold background), a 40% fixed threshold and a CT-based segmentation. Additionally, the activity recovered in the entire phantom, as well as the average activity concentration of the phantom background were compared to their true values. Finally, Monte-Carlo simulations of a commercial [Formula: see text]-camera were performed to investigate the accuracy of the TEW method. Good quantification accuracy (errors  <10%) was achieved for the entire phantom, the hot-background activity concentration and for objects in cold background segmented with a 1% threshold. However, the accuracy of activity quantification for objects segmented with 40% threshold or CT-based methods decreased (errors  >15%), mostly due to partial-volume effects. The Monte-Carlo simulations confirmed that TEW-scatter correction applied to ¹⁸⁸Re, although practical, yields only approximate estimates of the true scatter.

Analytically based photon scatter modeling for a multipinhole cardiac SPECT camera

PURPOSE

Dedicated cardiac SPECT scanners have improved performance over standard gamma cameras allowing reductions in acquisition times and/or injected activity. One approach to improving performance has been to use pinhole collimators, but this can cause position-dependent variations in attenuation, sensitivity, and spatial resolution. CT attenuation correction (AC) and an accurate system model can compensate for many of these effects; however, scatter correction (SC) remains an outstanding issue. In addition, in cameras using cadmium-zinc-telluride-based detectors, a large portion of unscattered photons is detected with reduced energy (low-energy tail). Consequently, application of energy-based SC approaches in these cameras leads to a higher increase in noise than with standard cameras due to the subtraction of true counts detected in the low-energy tail. Model-based approaches with parallel-hole collimator systems accurately calculate scatter based on the physics of photon interactions in the patient and camera and generate lower-noise estimates of scatter than energy-based SC. In this study, the accuracy of a model-based SC method was assessed using physical phantom studies on the GE-Discovery NM530c and its performance was compared to a dual energy window (DEW)-SC method.

METHODS

The analytical photon distribution (APD) method was used to calculate the distribution of probabilities that emitted photons will scatter in the surrounding scattering medium and be subsequently detected. APD scatter calculations for ^99mTc-SPECT (140 ± 14 keV) were validated with point-source measurements and 15 anthropomorphic cardiac-torso phantom experiments and varying levels of extra-cardiac activity causing scatter inside the heart. The activity inserted into the myocardial compartment of the phantom was first measured using a dose calibrator. CT images were acquired on an Infinia Hawkeye (GE Healthcare) SPECT/CT and coregistered with emission data for AC. For comparison, DEW scatter projections (120 ± 6 keV ) were also extracted from the acquired list-mode SPECT data. Either APD or DEW scatter projections were subtracted from corresponding 140 keV measured projections and then reconstructed with AC (APD-SC and DEW-SC). Quantitative accuracy of the activity measured in the heart for the APD-SC and DEW-SC images was assessed against dose calibrator measurements. The difference between modeled and acquired projections was measured as the root-mean-squared-error (RMSE). APD-modeled projections for a clinical cardiac study were also evaluated.

RESULTS

APD-modeled projections showed good agreement with SPECT measurements and had reduced noise compared to DEW scatter estimates. APD-SC reduced mean error in activity measurement compared to DEW-SC in images and the reduction was statistically significant where the scatter fraction (SF) was large (mean SF = 28.5%, T-test p = 0.007). APD-SC reduced measurement uncertainties as well; however, the difference was not found to be statistically significant (F-test p > 0.5). RMSE comparisons showed that elevated levels of scatter did not significantly contribute to a change in RMSE (p > 0.2).

CONCLUSIONS

Model-based APD scatter estimation is feasible for dedicated cardiac SPECT scanners with pinhole collimators. APD-SC images performed better than DEW-SC images and improved the accuracy of activity measurement in high-scatter scenarios.

Pericervical Injection of 99mTc-Nanocolloid Is Superior to Peritumoral Injection for Sentinel Lymph Node Detection of Endometrial Cancer in SPECT/CT

PURPOSE

Scintigraphic mapping of sentinel lymph node (SLN) is increasingly performed in patients with endometrial carcinoma although its routine clinical use is still under investigation. The purpose of this study was to compare preoperative SLN detection by means of SPECT/CT using pericervical (PC) versus hysteroscopic peritumoral (PT) injection.

PATIENTS AND METHODS

One hundred forty consecutive patients with endometrial carcinoma who underwent surgery and preoperative SLN SPECT/CT with Tc-nanocolloid were included. Seventy women received hysteroscopic injection at 3 PT sites, and 70 women received PC injection at 3-, 6-, 9-, and 12-o'clock positions. Each patient underwent SPECT/CT followed by modified radical hysterectomy with lymphadenectomy on the day after. Histopathological results were collected for validation.

RESULTS

Three hundred thirty-four SLNs were detected by SPECT/CT in 106 patients (mean, 3.15; range, 1-9). The detection rate after PC nanocolloid injection was 83% versus 69% after PT injection (Pearson χ test, P = 0.049). However, PT application resulted in a higher rate of para-aortic SLNs (PC: 60% vs PT: 38% of positive scans, P = 0.02). SPECT/CT yielded an overall sensitivity of 70% for the SLN detection in women with lymph node metastases with 3 false-negative cases. Failure to detect SLN was mostly associated with uptake in the reticuloendothelial system (liver, spleen, and bone marrow) or peritoneal diffusion in both cohorts. Negative scans after PT application often showed a minor to even failing injection depot.

CONCLUSIONS

Pericervical injection leads to a significantly better detection rate of SLN on SPECT/CT while reducing invasiveness of the injection procedure. Failure to detect SLN seems to be associated with major venous drainage.

Organ-specific SPECT activity calibration using 3D printed phantoms for molecular radiotherapy dosimetry

Background

Patient-specific absorbed dose calculations for molecular radiotherapy require accurate activity quantification. This is commonly derived from Single-Photon Emission Computed Tomography (SPECT) imaging using a calibration factor relating detected counts to known activity in a phantom insert.

Methods

A series of phantom inserts, based on the mathematical models underlying many clinical dosimetry calculations, have been produced using 3D printing techniques. SPECT/CT data for the phantom inserts has been used to calculate new organ-specific calibration factors for ^99mTc and ¹⁷⁷Lu. The measured calibration factors are compared to predicted values from calculations using a Gaussian kernel.

Results

Measured SPECT calibration factors for 3D printed organs display a clear dependence on organ shape for ^99mTc and ¹⁷⁷Lu. The observed variation in calibration factor is reproduced using Gaussian kernel-based calculation over two orders of magnitude change in insert volume for ^99mTc and ¹⁷⁷Lu. These new organ-specific calibration factors show a 24, 11 and 8 % reduction in absorbed dose for the liver, spleen and kidneys, respectively.

Conclusions

Non-spherical calibration factors from 3D printed phantom inserts can significantly improve the accuracy of whole organ activity quantification for molecular radiotherapy, providing a crucial step towards individualised activity quantification and patient-specific dosimetry. 3D printed inserts are found to provide a cost effective and efficient way for clinical centres to access more realistic phantom data.

Measuring total liver function on sulfur colloid SPECT/CT for improved risk stratification and outcome prediction of hepatocellular carcinoma patients

Background

Assessment of liver function is critical in hepatocellular carcinoma (HCC) patient management. We evaluated parameters of [^99mTc] sulfur colloid (SC) SPECT/CT liver uptake for association with clinical measures of liver function and outcome in HCC patients.

Methods

Thirty patients with HCC and variable Child-Turcotte-Pugh scores (CTP A5-C10) underwent [^99mTc]SC SPECT/CT scans for radiotherapy planning. Gross tumor volume (GTV), anatomic liver volume (ALV), and spleen were contoured on CT. SC SPECT image parameters include threshold-based functional liver volumes (FLV) relative to ALV, mean liver-to-spleen uptake ratio (L/S_mean), and total liver function (TLF) ratio derived from the product of FLV and L/S_mean. Optimal SC uptake thresholds were determined by ROC analysis for maximizing CTP classification accuracy. Image metrics were tested for rank correlation to composite scores and clinical liver function parameters. Image parameters of liver function were tested for association to overall survival with Cox proportional hazard regression.

Results

Optimized thresholds on SC SPECT were 58 % of maximum uptake for FLV, 38 % for L/S_mean, and 58 % for TLF. TLF produced the highest CTP classification accuracy (AUC = 0.93) at threshold of 0.35 (sensitivity = 0.88, specificity = 0.86). Higher TLF was associated with lower CTP score: TLF_A = 0.6 (0.4–0.8) versus TLF_B = 0.2 (0.1–0.3), p < 10⁻⁴. TLF was rank correlated to albumin and bilirubin (|R| > 0.63). Only TLF >0.30 was independently associated with overall survival when adjusting for CTP class (HR = 0.12, 95 % CI = 0.02–0.58, p = 0.008).

Conclusions

SC SPECT/CT liver uptake correlated with differential liver function. TLF was associated with improved overall survival and may aid in personalized oncologic management of HCC patients.

Multicenter evaluation of single-photon emission computed tomography quantification with third-party reconstruction software

Reliable and reproducible quantification is essential in many clinical situations. Previously, single-photon emission computed tomography (SPECT) has not been considered a quantitative imaging modality, but recent advances in reconstruction algorithm development have made SPECT quantitative. In this study, we investigate the reproducibility of SPECT quantification with phantoms in a multicenter setting using novel third-party reconstruction software. A total of five hospitals and eight scanners (three GE scanners and five Siemens scanners) participated in the study. A Jaszczak phantom without inserts was used to calculate counts to activity concentration conversion factors. The quantitative accuracy was tested using the NEMA-IEC phantom with six spherical inserts (diameters from 10 to 37 mm) filled to an 8 : 1 insert-background concentration ratio. Phantom studies were reconstructed at one central location using HERMES HybridRecon applying corrections for attenuation, collimator-detector response, and scatter. Spherical volumes of interest with the same diameter as the inserts were drawn on the images and recovery coefficients for the spheres were calculated. The coefficient of variation (CoV) of the NEMA-IEC phantom recovery coefficients ranged from ∼19 to 5% depending on the insert diameter so that the lowest CoV was obtained with the largest spheres. The intersite CoV was almost equal to intrasite CoV. In conclusion, quantitative SPECT is reproducible in a multicenter setting with third-party reconstruction software.

Quantitative accuracy of (177)Lu SPECT reconstruction using different compensation methods: phantom and patient studies

Background

In targeted radionuclide therapy (TRT), accurate quantification using SPECT/CT images is important for optimizing radiation dose delivered to both the tumour and healthy tissue. Quantitative SPECT images are regularly reconstructed using the ordered subset expectation maximization (OSEM) algorithm with various compensation methods such as attenuation (A), scatter (S) and detector and collimator response (R). In this study, different combinations of the compensation methods are applied during OSEM reconstruction and the effect on the ¹⁷⁷Lu quantification accuracy is studied in an anthropomorphic torso phantom. In addition, the phantom results are reflected to (177)Lu-DOTA-Tyr3-octreotate (¹⁷⁷Lu-DOTATATE)-treated patient data and kidney absorbed dose estimates.

Methods

The torso phantom was imaged with nine various sized (0.4–104.4 cm³) spherical inserts, filled with known ¹⁷⁷Lu activity ranging from 0.5 to 105.5 MBq. Images were reconstructed using OSEM algorithm using A, AR and ARS compensation method combinations. The compensation method combinations were compared by calculating the concentration recovery coefficient (cRC) for each insert. In addition, ten ¹⁷⁷Lu-DOTATATE-treated patient’s post-therapy dosimetry acquisitions were reconstructed, and the absorbed dose to kidneys was estimated.

Results

cRC values depend on the insert size for all compensation methods. AR and ARS produced significantly higher cRC values than attenuation correction alone. There were no cRC value differences between the methods for the smallest 1-cm-diameter insert, cRC being 0.18. However, the collimator and detector response compensation method (R) made the 1.3-cm-diameter insert clearly visible and improved cRC estimate from 0.19 to 0.43. ARS produced slightly higher cRC values for small- and medium-sized inserts than AR. On the patient data, a similar trend could be seen. AR and ARS produced higher kidney activities than using attenuation correction alone; the total absorbed doses to the right and left kidneys were on average 15 and 20 % higher for AR and 19 and 25 % higher for ARS, respectively. The effective half-life decay estimated from time-activity curves however showed no notable difference between the compensation methods.

Conclusions

The highest cRC values were achieved by applying ARS compensation during reconstruction. The results were notably higher than those using attenuation correction alone. Similarly, higher activity estimates and thus higher absorbed dose estimates were found in patient data when all compensation methods were applied. ARS improved cRC especially in small-sized sources, and it thus might aid tumour dosimetry for ¹⁷⁷Lu PRRT treatments.

Intravenous contrast-enhanced CT can be used for CT-based attenuation correction in clinical (111)In-octreotide SPECT/CT

BACKGROUND

CT-based attenuation correction (CT-AC) using contrast-enhancement CT impacts (111)In-SPECT image quality and quantification. In this study we assessed and evaluated the effect.

METHODS

A phantom (5.15 L) was filled with an aqueous solution of In-111. Three SPECT/CT scans were performed: (A) no IV contrast, (B) with 100-mL IV contrast, and (C) with 200-mL IV contrast added. Scan protocol included a localization CT, a low-dose CT (LD), and a full-dose CT (FD). Phantom, LD and FD scan series were performed at 90, 120, and 140 kVp. Phantom data were evaluated looking at mean counts in a central volume. Ten patients referred for (111)In-octreotide scintigraphy were scanned according to our clinical (111)In-SPECT/CT protocol including a topogram, a LD (140 kVp), and a FD (120 kVp). The FD/contrast-enhanced CT was acquired in both arterial (FDAP) and venous phase (FDVP) following a mono-phasic IV injection of 125-mL Optiray (4.5 mL/s). For patient data, we report image quality, Krenning scores, and mean/max values for liver and tumor regions.

RESULTS

Phantoms: in uncorrected emission data, mean counts (average ± SD) decreased with increasing IV concentration: (A) 119 ± 9, (B) 113 ± 8, and (C) 110 ± 9. For all attenuation correction (AC) scans, the mean values increased with increasing iodine concentration.

PATIENTS

there were no visible artifacts in single photon emission computed tomography (SPECT) following CT-AC with contrast-enhanced CT. The average score of image quality was 4.1 ± 0.3, 3.8 ± 0.4, and 4.2 ± 0.4 for LD, arterial phase, and venous phase, respectively. A total of 16 lesions were detected. The Krenning scores of 13/16 lesions were identical across all scan series. The max pixel values for the 16 lesions showed generally lower values for LD than for contrast-enhanced CT.

CONCLUSIONS

In (111)In-SPECT/CT imaging of phantoms and patients, the use of IV CT contrast did neither degrade the SPECT image quality nor affect the clinical Krenning score. Reconstructed counts in healthy liver tissues were unaffected, and there was a generally lower count value in lesions following CT-AC based on the LD non-enhanced images. Overall, for clinical interpretation, no separate low-dose CT is required for CT-AC in (111)In-SPECT/CT.

In vivo quantification of (177)Lu with planar whole-body and SPECT/CT gamma camera imaging

Background

Advances in gamma camera technology and the emergence of a number of new theranostic radiopharmaceutical pairings have re-awakened interest in in vivo quantification with single-photon-emitting radionuclides. We have implemented and validated methodology to provide quantitative imaging of ¹⁷⁷Lu for 2D whole-body planar studies and for 3D tomographic imaging with single-photon emission computed tomography (SPECT)/CT.

Methods

Whole-body planar scans were performed on subjects to whom a known amount of [¹⁷⁷Lu]-DOTA-octreotate had been administered for therapy. The total radioactivity estimated from the images was compared with the known amount of the radionuclide therapy administered. In separate studies, venous blood samples were withdrawn from subjects after administration of [¹⁷⁷Lu]-DOTA-octreotate while a SPECT acquisition was in progress and the concentration of the radionuclide in the venous blood sample compared with that estimated from large blood pool structures in the SPECT reconstruction. The total radioactivity contained within an internal SPECT calibration standard was also assessed.

Results

In the whole-body planar scans (n = 28), the estimated total body radioactivity was accurate to within +4.6 ± 5.9 % (range −17.1 to +11.2 %) of the correct value. In the SPECT reconstructions (n = 12), the radioactivity concentration in the cardiac blood pool was accurate to within −4.0 ± 7.8 % (range −16.1 to +7.5 %) of the true value and the internal standard measurements (n = 89) were within 2.0 ± 8.5 % (range −16.3 to +24.2 %) of the known amount of radioactivity contained.

Conclusions

In our hands, state-of-the-art hybrid SPECT/CT gamma cameras were able to provide accurate estimates of in vivo radioactivity to better than, on average, ±10 % for use in biodistribution and radionuclide dosimetry calculations.

(99m)Tc-Annexin A5 quantification of apoptotic tumor response: a systematic review and meta-analysis of clinical imaging trials

PURPOSE

(99m)Tc-Annexin A5 has been used as a molecular imaging probe for the visualization, characterization and measurement of apoptosis. In an effort to define the quantitative (99m)Tc-annexin A5 uptake criteria that best predict tumor response to treatment, we performed a systematic review and meta-analysis of the results of all clinical imaging trials found in the literature or publicly available databases.

METHODS

Included in this review were 17 clinical trials investigating quantitative (99m)Tc-annexin A5 (qAnx5) imaging using different parameters in cancer patients before and after the first course of chemotherapy and/or radiation therapy. Qualitative assessment of the clinical studies for diagnostic accuracy was performed using the QUADAS-2 criteria. Of these studies, five prospective single-center clinical trials (92 patients in total) were included in the meta-analysis after exclusion of one multicenter clinical trial due to heterogeneity. Pooled positive predictive values (PPV) and pooled negative predictive values (NPV) (with 95% CI) were calculated using Meta-Disc software version 1.4.

RESULTS

Absolute quantification and/or relative quantification of (99m)Tc-annexin A5 uptake were performed at baseline and after the start of treatment. Various quantitative parameters have been used for the calculation of (99m)Tc-annexin A5 tumor uptake and delta (Δ) tumor changes post-treatment compared to baseline including: tumor-to-background ratio (TBR), ΔTBR, tumor-to-noise ratio, relative tumor ratio (TR), ΔTR, standardized tumor uptake ratio (STU), ΔSTU, maximum count per pixel within the tumor volume (Cmax), Cmax%, absolute ΔU and percentage (ΔU%), maximum ΔU counts, semiquantitative visual scoring, percent injected dose (%ID) and %ID/cm(3). Clinical trials investigating qAnx5 imaging have included patients with lung cancer, lymphoma, breast cancer, head and neck cancer and other less common tumor types. In two phase I/II single-center clinical trials, an increase of ≥25% in uptake following treatment was considered a significant threshold for an apoptotic tumor response (partial response, complete response). In three other phase I/II clinical trials, increases of ≥28%, ≥42% and ≥47% in uptake following treatment were found to be the mean cut-off levels in responders. In a phase II/III multicenter clinical trial, an increase of ≥23% in uptake following treatment was found to be the minimum cut-off level for a tumor response. In one clinical trial, no significant difference in (99m)Tc-annexin A5 uptake in terms of %ID was found in healthy tissues after chemotherapy compared to baseline. In two other clinical trials, intraobserver and interobserver measurements of (99m)Tc-annexin A5 tumor uptake were found to be reproducible (mean difference <5%, kappa =  0.90 and 0.82, respectively) and to be highly correlated with treatment outcome (Spearman r = 0.99, p < 0.0001). The meta-analysis demonstrated a pooled positive PPV of 100% (95% CI 92 - 100%) and a pooled NPV of 70% (95% CI 55 - 82%) for prediction of a tumor response after the first course of chemotherapy and/or radiotherapy in terms of ΔU%. In a symmetric sROC analysis, the AUC was 0.919 and the Q* index was 85.21 %.

CONCLUSION

Quantitative (99m)Tc-annexin A5 imaging has been investigated in clinical trials for the assessment of apoptotic tumor responses. This meta-analysis showed a high pooled PPV and a moderate pooled NPV with ΔU cut-off values ranging between 20% and 30%. Standardization of quantification and harmonization of results are required for high-quality clinical research. A standardized uptake value score (SUV, ΔSUV) using quantitative SPECT/CT imaging may be a promising approach to the simple, reproducible and semiquantitative assessment of apoptotic tumor changes.

Comparison of internal dose estimates obtained using organ-level, voxel S value, and Monte Carlo techniques

PURPOSE

The authors' objective was to compare internal dose estimates obtained using the Organ Level Dose Assessment with Exponential Modeling (OLINDA/EXM) software, the voxel S value technique, and Monte Carlo simulation. Monte Carlo dose estimates were used as the reference standard to assess the impact of patient-specific anatomy on the final dose estimate.

METHODS

Six patients injected with 99mTc-hydrazinonicotinamide-Tyr3-octreotide were included in this study. A hybrid planar/SPECT imaging protocol was used to estimate 99mTc time-integrated activity coefficients (TIACs) for kidneys, liver, spleen, and tumors. Additionally, TIACs were predicted for 131I, 177Lu, and 90Y assuming the same biological half-lives as the 99mTc labeled tracer. The TIACs were used as input for OLINDA/EXM for organ-level dose calculation and voxel level dosimetry was performed using the voxel S value method and Monte Carlo simulation. Dose estimates for 99mTc, 131I, 177Lu, and 90Y distributions were evaluated by comparing (i) organ-level S values corresponding to each method, (ii) total tumor and organ doses, (iii) differences in right and left kidney doses, and (iv) voxelized dose distributions calculated by Monte Carlo and the voxel S value technique.

RESULTS

The S values for all investigated radionuclides used by OLINDA/EXM and the corresponding patient-specific S values calculated by Monte Carlo agreed within 2.3% on average for self-irradiation, and differed by as much as 105% for cross-organ irradiation. Total organ doses calculated by OLINDA/EXM and the voxel S value technique agreed with Monte Carlo results within approximately ±7%. Differences between right and left kidney doses determined by Monte Carlo were as high as 73%. Comparison of the Monte Carlo and voxel S value dose distributions showed that each method produced similar dose volume histograms with a minimum dose covering 90% of the volume (D90) agreeing within ±3%, on average.

CONCLUSIONS

Several aspects of OLINDA/EXM dose calculation were compared with patient-specific dose estimates obtained using Monte Carlo. Differences in patient anatomy led to large differences in cross-organ doses. However, total organ doses were still in good agreement since most of the deposited dose is due to self-irradiation. Comparison of voxelized doses calculated by Monte Carlo and the voxel S value technique showed that the 3D dose distributions produced by the respective methods are nearly identical.

Estimation of radioactivity in single-photon emission computed tomography for sentinel lymph node biopsy in a torso phantom study

OBJECTIVE

The number of lymph nodes to be removed is determined from residual counts. Advance estimation of residual radioactivity in lymphatic nodes before a biopsy is useful for reducing surgical operation time. The purpose of this study was to estimate the total radioactivity of a small hotspot in single-photon emission computed tomography (SPECT) of a torso phantom.

METHODS

A cross-calibration study was performed to convert counts in SPECT images to radioactivity. A simulation study was performed to estimate the size of the volume of interest (VOI) covering a hotspot corrupted with full-width at half-maximum between 8 and 16 mm. The estimation of total radioactivity was validated in a torso phantom study using small sources.

RESULTS

True radioactivity was approximately equal to integrated values of hotspots using the VOI with a diameter of 40 mm in our simulation study. The difference was less than 18% in cases of more than 9.4 kBq.

CONCLUSION

The total radioactivity in small sources simulating a typical sentinel node was estimated from SPECT images using a VOI of 40 mm in a torso phantom study. Because the difference from actual values was less than 10% on average when radioactivities were more than 9.4 kBq, the total radioactivity of a lymph node can be estimated in a clinical examination.

Initial study of radiological and clinical efficacy radioembolization using 188Re-human serum albumin (HSA) microspheres in patients with progressive, unresectable primary or secondary liver cancers

Background

The aim of this initial study was to evaluate the clinical and radiological effectiveness of radioembolization (RE) using ¹⁸⁸Re-Human Serum Albumin (HSA) microspheres in patients with advanced, progressive, unresectable primary or secondary liver cancers, not suitable to any other form of therapy.

Material/Methods

Overall, we included 13 patients with 20 therapy sessions. Clinical and radiological responses were assessed at 6 weeks after therapy, and then every 3 months. The objective radiological response was classified according to Response Evaluation Criteria in Solid Tumors (RECIST) v.1.0 by sequential MRI. Adverse events were evaluated using NCI CTCAE v.4.03.

Results

There were 4 patients with hepatocellular carcinoma (HCC), 6 with metastatic colorectal cancer (mCRC), 2 with neuroendocrine carcinoma (NEC), and 1 patient with ovarian carcinoma. Mean administered activity of ¹⁸⁸Re HSA was 7.24 GBq (range 3.8–12.4) A high microspheres labeling efficacy of over 97±2.1% and low urinary excretion of ¹⁸⁸Re (6.5±2.3%) during first 48-h follow-up. Median overall survival (OS) for all patients was 7.1 months (CI 6.2–13.3) and progression-free survival (PFS) was 5.1 months (CI 2.4–9.9). In those patients who had a clinical partial response (PR), stable disease (SD), and disease progression (DP) as assessed 6 weeks after therapy, the median OS was 9/5/4 months, respectively, and PFS was 5/2/0 months, respectively. The treatment adverse events (toxicity) were at an acceptable level. Initially and after 6 weeks, the CTC AE was grade 2, while after 3 months it increased to grade 3 in 4 subjects. This effect was mostly related to rapid cancer progression in this patient subgroup.

Conclusions

The results of this preliminary study indicate that RE using ¹⁸⁸Re HSA is feasible and a viable option for palliative therapy in patients with extensive progressive liver cancer. It was well tolerated by most patients, with a low level of toxicity during the 3 months of follow-up.

JADA: a graphical user interface for comprehensive internal dose assessment in nuclear medicine

PURPOSE

The main objective of this work was to design a comprehensive dosimetry package that would keep all aspects of internal dose calculation within the framework of a single software environment and that would be applicable for a variety of dose calculation approaches.

METHODS

Our MATLAB-based graphical user interface (GUI) can be used for processing data obtained using pure planar, pure SPECT, or hybrid planar/SPECT imaging. Time-activity data for source regions are obtained using a set of tools that allow the user to reconstruct SPECT images, load images, coregister a series of planar images, and to perform two-dimensional and three-dimensional image segmentation. Curve fits are applied to the acquired time-activity data to construct time-activity curves, which are then integrated to obtain time-integrated activity coefficients. Subsequently, dose estimates are made using one of three methods.

RESULTS

The organ level dose calculation subGUI calculates mean organ doses that are equivalent to dose assessment performed by OLINDA/EXM. Voxelized dose calculation options, which include the voxel S value approach and Monte Carlo simulation using the EGSnrc user code DOSXYZnrc, are available within the process 3D image data subGUI.

CONCLUSIONS

The developed internal dosimetry software package provides an assortment of tools for every step in the dose calculation process, eliminating the need for manual data transfer between programs. This saves times and minimizes user errors, while offering a versatility that can be used to efficiently perform patient-specific internal dose calculations in a variety of clinical situations.

Three-dimensional personalized dosimetry for 188Re liver selective internal radiation therapy based on quantitative post-treatment SPECT studies

We demonstrate that accurate patient-specific distributions of microspheres labeled with ¹⁸⁸Re and resulting absorbed doses can be obtained from single-photon emission computed tomography (SPECT) studies performed after ¹⁸⁸Re selective internal radiation therapy when accurate correction methods are employed in image reconstruction. Our quantitative image reconstruction algorithm includes corrections for attenuation, resolution degradations and scatter as well as a window-based compensation for contamination. The procedure has been validated using four phantom experiments containing an 18 ml cylindrical source (82-93 MBq of ¹⁸⁸Re activity) simulating a liver tumor. In addition, we applied our approach to post-therapy SPECT studies of ten patients with progressive primary or metastatic liver carcinomas. Our quantitative algorithm accurately (within 9%) recovered ¹⁸⁸Re activity from four phantom experiments. In addition, for two patients that received three scans, deviations remained consistent between the measured and the reconstructed activities that were determined from studies with differing severity of the dead-time effect. The analysis of absorbed doses for patient studies allowed us to hypothesize that D90 (the minimum dose received by 90% of the tumor volume) may be a reliable metric relating therapy outcomes to the calculated doses. Among several considered metrics, only D90 showed statistically significant correlation with the overall survival.

MIRD pamphlet No. 24: Guidelines for quantitative 131I SPECT in dosimetry applications

The reliability of radiation dose estimates in internal radionuclide therapy is directly related to the accuracy of activity estimates obtained at each imaging time point. The recently published MIRD pamphlet no. 23 provided a general overview of quantitative SPECT imaging for dosimetry. The present document is the first in a series of isotope-specific guidelines that will follow MIRD 23 and focuses on one of the most commonly used therapeutic radionuclides, (131)I. The purpose of this document is to provide guidance on the development of protocols for quantitative (131)I SPECT in radionuclide therapy applications that require regional (normal organs, lesions) and 3-dimensional dosimetry.

Quantitative SPECT/CT: SPECT joins PET as a quantitative imaging modality

The introduction of combined modality single photon emission computed tomography (SPECT)/CT cameras has revived interest in quantitative SPECT. Schemes to mitigate the deleterious effects of photon attenuation and scattering in SPECT imaging have been developed over the last 30 years but have been held back by lack of ready access to data concerning the density of the body and photon transport, which we see as key to producing quantitative data. With X-ray CT data now routinely available, validations of techniques to produce quantitative SPECT reconstructions have been undertaken. While still suffering from inferior spatial resolution and sensitivity compared to positron emission tomography (PET) imaging, SPECT scans nevertheless can be produced that are as quantitative as PET scans. Routine corrections are applied for photon attenuation and scattering, resolution recovery, instrumental dead time, radioactive decay and cross-calibration to produce SPECT images in units of kBq.ml(-1). Though clinical applications of quantitative SPECT imaging are lacking due to the previous non-availability of accurately calibrated SPECT reconstructions, these are beginning to emerge as the community and industry focus on producing SPECT/CT systems that are intrinsically quantitative.

The accuracy and reproducibility of SPECT target volumes and activities estimated using an iterative adaptive thresholding technique

OBJECTIVE

Our aim was to design a practical and reproducible image segmentation method for calculations of total absorbed doses in organs and tumours for internally delivered radioisotopes. We have built upon our previously proposed use of two separate thresholds and employed an iterative technique for semiautomatic selection of background regions for segmenting an object of interest using thresholds that depend on the source-to-background ratio of activity concentrations.

METHODS

The parameters of curves relating volume and activity thresholds to source-to-background ratio were established using phantoms with 20 different inserts. The accuracy of our technique was validated using a second phantom experiment, whereas the reproducibility of volume, activity and dose estimates of organs and tumours was investigated using 13 patient studies. The accuracy and reproducibility of segmentations achieved were assessed using images reconstructed with three different methods that ranged from a standard clinical reconstruction to an advanced quantitative reconstruction approach.

RESULTS

In the validation phantom experiment, bottle volumes and activities measured using iterative adaptive thresholding agreed on average with the true values to within 4%, regardless of the reconstruction method used. In the patient studies, volumes and activities estimated from the single-photon emission computed tomography images reconstructed with clinical software agreed with the volumes and activities estimated using the advanced reconstruction approach to within 6%, whereas the corresponding doses agreed to within 4%.

CONCLUSION

The proposed iterative adaptive thresholding technique can accurately determine object volume and activity, which allows standard clinical reconstructions to generate absorbed dose estimates that are similar to those values obtained using more advanced reconstruction methods.

Personalized image-based radiation dosimetry for routine clinical use in peptide receptor radionuclide therapy: pretherapy experience

Patient-specific dose calculations are not routinely performed for targeted radionuclide therapy procedures, partly because they are time consuming and challenging to perform. However, it is becoming widely recognized that a personalized dosimetry approach can help plan treatment and improve understanding of the dose-response relationship. In this chapter, we review the procedures and essential elements of an accurate internal dose calculation and propose a simplified approach that is aimed to be practical for use in a busy nuclear medicine department.

An evidence-based review of quantitative SPECT imaging and potential clinical applications

SPECT has traditionally been regarded as nonquantitative. Advances in multimodality γ-cameras (SPECT/CT), algorithms for image reconstruction, and sophisticated compensation techniques to correct for photon attenuation and scattering have, however, now made quantitative SPECT viable in a manner similar to quantitative PET (i.e., kBq cm(-3), standardized uptake value). This review examines the evidence for quantitative SPECT and demonstrates clinical studies in which the accuracy of the reconstructed SPECT data has been assessed in vivo. SPECT reconstructions using CT-based compensation corrections readily achieve accuracy for (99m)Tc to within ± 10% of the known concentration of the radiotracer in vivo. Quantification with other radionuclides is also being introduced. SPECT continues to suffer from poorer photon detection efficiency (sensitivity) and spatial resolution than PET; however, it has the benefit in some situations of longer radionuclide half-lives, which may better suit the biologic process under examination, as well as the ability to perform multitracer studies using pulse height spectroscopy to separate different radiolabels.

Two methods to generate templates for template-based partial volume effect correction: SPECT phantom experiments

In this paper, we explore the applicability of template-based compensation for the partial volume effect (PVE) for situations where (i) the image has multiple uptake sites (tumors and organs) but only one of them is treated as a region of interest (ROI) with the boundaries available from a high-resolution modality and (ii) no information regarding activities inside or outside this ROI is a priori available. We modeled this situation by performing SPECT acquisitions of phantoms containing 21 containers, which had different shapes and sizes and were filled with different levels of activity. In our analysis, each of these containers was treated as an individual ROI. We compared the performance of two methods of template construction. In method 1, the ROI template value was obtained from a conventionally reconstructed (without PVEC) image. In method 2, the ROI template value was directly (bypassing the PVE-affected conventional image) calculated from projections using region-based reconstruction. Our processing shows that method 1 resulted in consistent (activities for all 21 ROIs were improved) but relatively weak PVE compensation (errors of recovered total activities were equal to or lower than 10% for 5 ROIs only). Application of method 2 resulted in a selective (activities for 19 ROIs were improved) but considerably better compensation when compared to method 1 (errors of recovered total activities were equal to or lower than 10% for 10 ROIs).

Scanning linear estimation: improvements over region of interest (ROI) methods

In tomographic medical imaging, a signal activity is typically estimated by summing voxels from a reconstructed image. We introduce an alternative estimation scheme that operates on the raw projection data and offers a substantial improvement, as measured by the ensemble mean-square error (EMSE), when compared to using voxel values from a maximum-likelihood expectation-maximization (MLEM) reconstruction. The scanning-linear (SL) estimator operates on the raw projection data and is derived as a special case of maximum-likelihood estimation with a series of approximations to make the calculation tractable. The approximated likelihood accounts for background randomness, measurement noise and variability in the parameters to be estimated. When signal size and location are known, the SL estimate of signal activity is unbiased, i.e. the average estimate equals the true value. By contrast, unpredictable bias arising from the null functions of the imaging system affect standard algorithms that operate on reconstructed data. The SL method is demonstrated for two different tasks: (1) simultaneously estimating a signal's size, location and activity; (2) for a fixed signal size and location, estimating activity. Noisy projection data are realistically simulated using measured calibration data from the multi-module multi-resolution small-animal SPECT imaging system. For both tasks, the same set of images is reconstructed using the MLEM algorithm (80 iterations), and the average and maximum values within the region of interest (ROI) are calculated for comparison. This comparison shows dramatic improvements in EMSE for the SL estimates. To show that the bias in ROI estimates affects not only absolute values but also relative differences, such as those used to monitor the response to therapy, the activity estimation task is repeated for three different signal sizes.