2.1. SRS imaging protocol
All images were acquired on a Symbia Intevo 16 (Siemens Medical Solution, Erlangen, Germany) with low middle-energy general purpose (LMEGP) collimators. The energy windows for 111In-SPECT were selected as 172 keV ± 7.5% and 247 keV ± 7.5% for the main window, 15% and 8% for the lower and upper windows of 172 keV, respectively, and 10% for the lower window of 247 keV (Fig. 1). SPECT scan was performed for 20 min (60 steps, 40 s/step, 128×128 matrix; magnification, 1.0; pixel size, 4.8 mm). SPECT images were obtained using the ordered subset expectation maximization (OSEM) method (Flash3D; Siemens Healthcare, Erlangen, Germany) with 12 iterations and 6 subsets. Images were smoothed using a Gaussian filter (9.6 mm at full width half maximum). Scattering correction was performed using the triple energy window (TEW) method at 171 keV and dual energy window (DEW) method at 245 keV. Attenuation correction was performed using computed tomography attenuation correction (CTAC). CT acquisition for attenuation correction was performed using adaptive dose modulation (CARE Dose 4D; Siemens Healthcare, Erlangen, Germany) at a tube voltage of 130 kV and tube current of 80 mA. The CT data were reconstructed to a slice thickness of 3 mm, and medium sharp and attenuation kernels (B50s and B31s, respectively) were used. All image analyses were performed using the Daemon Research Image Processor (DRIP; Fujifilm Toyama Chemical Co., Ltd., Tokyo, Japan) version 3.0.2.0, and OsiriX software (Pixmeo, Bernex, Switzerland) v. 12.0.3.
2.2. Phantom studies
2.2.1. Phantom design
The characteristics of each energy peak were confirmed by phantom experiments. A NEMA IEC body phantom (Data Spectrum, NC, USA) was used to simulate the upper abdomen. Six hot spheres (φ37, 28, 22, 17, 13, and 10 mm) and lung insert (φ44 mm) were placed in the phantom. The renal excretion rate of 111In-pentetreotide in the human body after 24 h is 85%, and the uptake in the upper abdomen relative to the whole body is approximately 60 %–70% [19]. The residual rate of 111In-pentetreotide in the body after 24 h is approximately 20 MBq. Therefore, the entire phantom contained 13 MBq of the 111In solution (Fig. 2).
The cross-calibration factor (CCF), which is necessary to calculate the quantitative value, was measured. A polyethylene bottle (φ95 mm, volume = 1170 mL) was used for the CCF; The bottle was filled with 3.47 kBq/mL of 111In solution for the measurement.
2.2.2. Data analysis
Three types of SPECT images were used: images acquired at 172 keV ± 7.5% and 247 keV ± 7.5% in the main window (171 keV image and 245 keV image, respectively), and images acquired by summing the two windows (sum image). To clarify the PVE, the RC was determined for each window. A new quantitative value, the indium uptake index (IUI), was calculated by correcting the RC.
Regions of interest (ROIs) were placed at the hot spheres and lung insert in the SPECT images for each window, adjusting the position and size with reference to the CT images. The shape of the ROI was a circle that matched the shape of the hot sphere and lung insert. The hot sphere was placed in the slice where it was most depicted in the three slices, before and after the slice corresponding to the hot sphere of the lung insert (Fig. 3). The maximum counts in the three slices before and after each hot sphere and lung insert were measured from the ROIs. The RCs were calculated as follows:
$$RC = \frac{{C}_{max,i }}{{C}_{max,44}}$$
1
C max,44 is the maximum count in the lung insert, and Cmax,i is the maximum count in each hot sphere.
Based on the relationship between the calculated RCs and the diameters of the hot spheres and lung insert, spline interpolation was performed between the measurements to observe a total of 60 RCs using the original program written in Python. Wilcoxon's signed rank sum test was used to examine the differences between the three measured RCs. P < 0.05 was considered statistically significant.
To correct the PVE for accumulations of 44 mm or less in diameter, the maximum count of the hot spheres and lung insert must be divided by the RC corresponding to their diameter. In calculating the quantitative value, the dosage and weight of the patient must be considered. The standardized uptake value (SUV) method is commonly used to calculate quantitative values independent of body weight and dose [22]. In this study, the IUI was calculated with reference to this method. The IUI of each hot sphere and lung insert was calculated as follows:
$$IUI = \frac{\frac{{C}_{max i,energy}}{{RC}_{i,energy}} }{Dose radioactivity \left(Bq\right) / Weight \left(g\right)}\times CCF$$
2
C max i,energy is the maximum count of each hot sphere in each window, RCi,energy is the RC corresponding to the diameter in each window, Dose radioactivity is the amount of radioactivity actually contained in the phantom (Bq), and Weight is the weight of the phantom (g). CCF is calculated from the radioactivity concentration in the bottle and the counts from the SPECT images.
In this study, 17 mm spheres were measured to confirm the quantitative accuracy for accumulations of less than 20 mm in diameter. The quantitative values for the 17 mm sphere were calculated using two patterns: IUI and SUV. In addition, the percentage (%) difference was calculated to confirm the difference between the ideal quantitative value and calculated quantitative value. The %difference was calculated as follows:
$$\% difference = \frac{{Index}_{ref}-{Index}_{cal}}{{Index}_{ref}} \times 100$$
3
Index ref is the ideal quantitative value of IUI and SUV (IUI, SUV = 8.8), and Indexcal is the quantitative value of IUI and SUV by calculation.
2.3. Case study
2.3.1. Patient protocol
All patients underwent SRS at the Cancer Institute Hospital of JFCR between April 2016 and March 2020. Fourteen patients (17 sites) were included in the study. Four patients had normal liver metastases, seven had liver metastases (two pancreatic tumors), and five had pancreatic head and body tumors (two liver metastases). In all patients, 111In-pentetreotide (Octreoscan; Fujifilm Toyama Chemical Co., Ltd., Tokyo, Japan) was administered at a dose of 180.1 ± 14.8 MBq. The selection criteria for patients were as follows: at least one SPECT/CT session performed 24 hours after 111In-pentetreotide administration, no liver lesions on modalities other than SPECT/CT in patients with normal liver, and no accumulation of lesions on SPECT/CT in patients with liver metastases or pancreatic head or body tumors. In cases of liver metastases and pancreatic head and body tumors, SPECT/CT showed accumulation in the lesions.
2.3.2. Data analysis
The IUI was calculated for all 17 sites using the RC and CCF calculated by the phantom experiments.
For the calculation of IUI in normal liver cases, ROIs of a size that enclosed the liver were placed visually in the liver area of the case images (Fig. 4). The IUI in the liver region was calculated for each window using the obtained maximum count, dose for each case, and body weight.
For the accumulation of liver metastases and tumors in the head and body of the pancreas, ROIs were set according to the size of the tumor, and the maximum counts in the tumor were measured (Fig. 5). Tumor size was measured using the previous CT and MRI. Based on the measured length and diameter of each tumor, the RC was selected, and the IUI was calculated using the maximum count.
Table 1
Definition for Krenning score (0–4)
Score
|
Intensity
|
0
|
None (no uptake)
|
1
|
Very low
|
2
|
Less than or equal to that of the liver
|
3
|
Greater than that of the liver
|
4
|
Greater than that of the spleen
|
Krenning scores for all sites were determined according to Table 1. For cases with normal liver, the score was set to 2 based on the criteria of the Krenning score. The correlation coefficient was used to clarify the relationship between the Krenning score and IUI.
2.4.3. Statistical analysis
To check for differences between the groups, Friedman's test was used for the calculated IUI for each energy window. To clarify the relationship between the Krenning score and IUI, box plot analysis was used, and Spearman's rank correlation coefficient was calculated. Statistical significance was set at p < 0.05. All statistical analyses were performed using Easy R (Saitama Medical Center, Jichi Medical University, Saitama, Japan) version 1.54, and the graphical user interface of R (The R Foundation for Statistical Computing, Vienna, Austria) version 3.6.2 [23].