- 1Department of Nuclear Medicine, Gødstrup Hospital, Herning, Denmark
- 2Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
- 3University Clinic in Nephrology and Hypertension, Department of Medicine, Gødstrup Hospital and Aarhus University, Herning, Denmark
Introduction: The aim of this study was to investigate the impact of the use of varying input parameters on resulting bone plasma clearance (Ki) and other kinetic modelling parameters in a group of patients with chronic kidney disease–mineral and bone disorder (CKD-MBD).
Methods: Raw PET/CT data and blood data were systematically analyzed using five different VOIs for the input functions in the left ventricle and in the thoracic aorta. Standardized VOIs were placed in four thoracic vertebrae and the results pooled and averaged. The basic image-derived input functions (IDIFs) were corrected for partial volume effect and spill-over and modified by substitution of the terminal image exponential with the corresponding plasma-exponentials derived from blood samples. Ki was then calculated using both a non-linear regression (NLR) analysis and a graphical Patlak analysis and compared.
Results: Our original results were reproducible with an inter-observer difference of approximately 6%. The correction factors varied with the VOI volumes from 0.73 ± 0.17 for the largest LV-VOI (48.7 ± 25.3 cm3) to 0.99 ± 0.10 for the AO-VOI (3.4 ± 1.2 cm3). The mean NLR-Ki results varied between 0.0378 ± 0.0112 and 0.0432 ± 0.0095 ml/min ml−1 with a fixed vB and 0.0408 ± 0.0111 and 0.045 ± 0.0102 ml/min ml−1 with a free-fitted vB. The corresponding Patl-Ki-results varied between 0.0302 ± 0.0071 and 0.0325 ± 0.0070 ml/min ml−1, having lesser differences and variances. The input functions with least variance and mean differences compared with NLR results were derived from the left ventricle with a VOI volume of 19.2 ± 11.3 cm3 corrected for PVE and Bg with a mean Ki-difference: 0.0097 ± 0.0370 ml/min ml−1 and 95% confidence limits (−0.023 to 0.004).
Conclusions: Our results indicated that a VOI with a volume of approximately 20 cm3 with a correction factor of 0.83 ± 0.13 results in Patlak results with the least variance and difference compared with the NLR results. The use of free-fitted vB in the NLR analysis showed the most robust results in all input series. The Patlak results were in comparison generally lower than the NLR results (−17.3% to −23.4%) but very robust across the various input series and with results comparable to previously published data and are therefore recommended for future analysis.
1. Introduction
The practical implementation of correct and accurate quantification of dynamic positron emission and computed tomographic (PET/CT) tracer kinetic studies of metabolism is non-trivial, the success of which is dependent on the choice of many analysis-input parameters, both scanner and protocol related, all of which need to be understood and optimized when locally implementing a new examination method for either research or clinical purposes.
Historically, the use of dynamic PET/CT tracer kinetic studies has been restricted primarily to research purposes at larger, specialized university departments with many years of experience in the necessary local implementation and development of analysis protocols and software, with the clinical use of dynamic PET/CT techniques being restricted to a few commercially developed functional imaging packages for a limited range of organs, e.g., 18F-FDG-PET/CT (fluorodeoxyglucose) and 82Rb-PET/CT (82Rubidium) for myocardial perfusion examinations. Currently, with the more widespread availability of PET/CT scanner functional options (list-mode data acquisition, vendor supplied re-binning reconstruction software, simplified kinetic analysis software tools), the use of a kinetic analysis of dynamic PET/CT acquisitions for other organs, with various radiopharmaceutical tracers, is becoming an attractive and attainable possibility in the more general clinical setting. However, before deviating from the use of the commercially available application packages, a number of technical questions need to be addressed.
Within the field of studying bone metabolism, a recent paper by Puri et al. (1) answered in detail some of the questions we needed answers to when we, more than 5 years ago, implemented dynamic tracer kinetic studies of metabolism using 18F-NaF-PET/CT in a group of patients with chronic kidney disease–mineral bone disorder (CKD-MBD). The goals of our original study (2) were to implement methods for determining bone plasma clearance (Ki) in this patient population, based on the methods as previously applied to a population of patients with osteoporosis (3–5) and ultimately to derive a standardized input curve (“semi-population input function,” see Section 2.1.4) for future routine use in a site-specific analysis of patients with CKD-MBD (6).
To achieve these original goals, our volume-of-interest (VOI) definition used a contour in the range of 45%–65% of the maximum value within a box limiting the left ventricle of the heart (LV). This resulted in a VOI filling most of the LV, which was separated from the background VOI in the myocardial wall by at least 2 voxels (≥6.4 mm). The VOIs were subject to partial volume effect (PVE) and spill-over of background activity (Bg) to and from the surrounding myocardial wall and it was necessary to correct for PVE and Bg as described by Cook et al. (3) through the use of a measured recovery coefficient (RCß). The mean subject-specific measured RCβ was 0.69 ± 0.15 (2), which seemed a little low for a modern PET/CT scanner when using 18F-NaF as a tracer.
In addition, our original work used a fixed blood-volume fraction parameter (vB) of 0.05 instead of a free-fitted vB in the non-linear regression analysis (NLR). Our implemented method of substituting the final image exponentials of the input curves with exponentials derived from plasma samples using logarithmic transformation differed from that described by Frost et al. (4) and Blake et al. (6).
In the review process for publication of this original study, this selection of parameters and method implementation, which were based on the available literature at the time (3, 6), was questioned, raising discussion and debate regarding the following: (1) the choice of VOI definition and size, (2) the best/most correct way to correct for PVE and spill-over from the background, and (3) the most correct value for the fractional blood volume parameter (vB).
As these issues were considered limitations in our original study, this work presents an attempt to improve and clarify the following points: (1) the effects of VOI definition, size, and correction factor on the input functions and the dynamic results; (2) the difference caused by the use of fixed vB vs. free-fitted vB; (3) whether the results obtained with our “logarithmic multiplicative method” are different from the results obtained with the original “exponential additive method” as described by Blake et al. (6); (4) whether the inclusion of an additional blood sample at 90 mpi (minutes post injection) at the end of the whole-body scan results in a better fit between image and blood data; and (5) whether a semi-population function (SP-function) derived from the optimized analysis of CKD-MBD patients differs from the corresponding SP-function derived from patients with osteoporosis.
The latter point is important, in order to indicate the necessity, or not, of requiring separate SP-functions for individual disease populations, as differences in bone turn-over may affect the shape, and thus the area under the curve (AUC) of the plasma curve (7).
2. Methods
2.1. Patients
Raw data from 12 patients with CKD-MBD enrolled in our previous study of method implementation (2) were reanalyzed for reproducibility. For the optimized analyses in this study, one obese patient was excluded due to extremely poor counting statistics resulting in outlying data in all series with small VOI volumes ≤1 ml.
2.2. Blood samples
In the original study, 5-ml venous blood samples were collected at −5, 30, 40, 50, 60, and 90 min after injection and prepared for well counting. The well-counter and PET/CT scanner were cross-calibrated as previously described by Vrist et al. (2). Whole blood and plasma data from our original study were reused for this study.
To convert measured activity from image-derived whole blood to plasma activity curves, plasma to whole blood activity ratios (PWR) were calculated for each of the samples.
In addition, all plasma values were transformed using the natural logarithm function. The slope and intercept of the resulting line at 40–60 and 40–90 mpi (plasma exponentials) were then determined by linear regression analysis and extrapolated back to the time for the peak. The plasma exponentials were used as a substitution for the corresponding image exponentials for the construction of the various input functions as described below. The 90-min samples were included for two reasons:
1. Interpolation of plasma data for calculation of Ki at the time of whole-body (WB) data acquisition, instead of extrapolation.
2. Comparison of the Ki results obtained using the 40–60 mpi plasma-curve substitution with results using the 40–90 mpi plasma-curve substitution, as in theory, the later sampling should result in a better fit to the original curve.
2.3. Image acquisition
The original PET/CT images were acquired on a Siemens Biograph mCT-4R 64 slice PET/CT scanner with a 22-cm axial field of view (FOV). The participants were positioned with the heart and the thoracic vertebrae Th7–Th10 centered in the FOV. After an intravenous bolus injection of 150 MBq 18F-NaF flushed with 20 ml isotonic saline, a 60-min list-mode dynamic scan was acquired immediately followed by a WB scan from the middle of the femur to the vertex of the skull acquired in 6–7 FOVs of 3 min per bed position.
2.4. Image reconstruction
The original PET images for dynamic analysis were re-binned into 50-time frames: 20 × 3 s, 12 × 5 s, 4 × 30 s, and 14 × 240 s. The reconstruction of PET scans used filtered back-projection with a Gaussian filter of 5 mm and a matrix size of 256 × 256 (3.2 mm × 3.2 mm × 1.4 mm). All dynamic images were automatically decay-corrected to the study injection time (study reference time). Image data from the WB scan were automatically decay-corrected to the start of the WB scan requiring additional decay correction to the study injection time for comparison with dynamic data.
Low-dose CT scans were performed, and the images reconstructed in three utilization-dependent series: (1) attenuation correction, (2) localization and identification of the thoracic vertebrae in the dynamic scan, and (3) localization of the relevant bone regions in the WB scan.
2.5. Image analysis
PMOD® version 4.206 software (PMOD Technologies LLC, Switzerland) was used for the non-linear regression analysis of the dynamic data and analysis of the static WB data.
2.6. Bone VOIs
For the reanalysis of the original data, all VOIs were constructed as described in the original study (2), whereas for the new optimization studies the vertebral VOIs (Th7/8–Th10/11) were drawn using a circular region of interest (ROI) in a single slice of 1-voxel thickness (1.4 mm) with a diameter of 4 voxels (12.8 mm), which was centered in the spongious bone and propagated through six slices resulting in a cylinder with a volume of approximately 1 ml in each of the four vertebrae (Figure 1A). Care was taken to avoid cortical bone and areas of obvious abnormal bone turnover (e.g., in a compressed vertebrae or for near lying large osteophytes). The data from the four individual VOIs were then unified to one combined VOI using the Union function in the PMOD VOI tools to improve counting statistics.
Figure 1. The various VOIs. (A) Vertebrae Th7-Th10; (B) CT-image–My-ocardial Bg (blue); (C) AO-Bg (pink); (D) LV-Orig (yellow); (E) LV-New (yellow); (F) LV-Fix (yellow); (G) AO-Fix (yellow); (H) AO-Peak (yellow).
These combined VOIs were kept constant for each subject dataset and used throughout all studies of the various input functions described below.
2.7. Background VOIs
2.7.1. LV studies
The background (myocardial Bg-VOI) was a ROI drawn on a single slice using the brush tool with a width of 1–3 pixels and propagation of the ROI over at least six slices. Its placement in the myocardial wall and as far as possible from the LV contour was controlled using the CT scan in the optimization studies (Figure 1B).
2.7.2. Aorta studies
Two cylindrical VOIs were centered around the aorta at a level not lower than the middle of vertebrae Th7. The inner cylinder had a diameter (38–40 mm) approximately 2–3 voxels wider than the activity in the aorta and the outer cylinder had a diameter 4 mm wider than the inner cylinder. The lengths of the cylinders were 30–40 mm, depending on the curvature of the aorta and the necessity to avoid inclusion of, for example, the intercostal arteries, which might result in a false high background activity (Figure 1C). The background activity (CBg) was calculated as the ratio between the differences in activity and the volume between inner and outer cylinders:
2.8. Input VOIs for image-derived input functions (IDIFs)
The VOIs for the IDIFs were constructed in five different ways. In the two studies using hot contouring of activity in the LV, the LV was delimited using a box VOI (approximately 40 mm × 40 mm × 40 mm). Within this VOI, a hot contour was drawn using PMOD's contour tool with the cutoff values described below.
The VOI definitions for these basic (unmodified) input functions were:
1. For reanalysis with original VOI definition (LV-Orig): a hot contour delineating 50%–70% of the max value in the box (Figure 1D). VOI volume: 50.8 ± 25.2.All contours were visually inspected for overlap with myocardial background and having a distance of at least two voxels (∼6.4 mm) from the myocardial background in all planes.
2. For LV-VOIs less prone to PVE and spill-over (LV-New): a hot contour delineating 60%–70% of the max value in the box (Figure 1E). VOI volume: 19.2 ± 11.3.
3. For LV-Fixed VOI studies (LV-Fix): A 12.8-mm diameter spherical VOI (∼4 voxels) was placed over the volume with highest activity in the center of the left ventricle (Figure 1F). VOI volume: 1.0 ± 0.6.
4. For the Aorta-Fixed VOI studies (AO-Fix): a cylindrical VOI with a diameter of 4 voxels (12.8 mm) was centered between the background VOIs described above, with the same length as the background VOI (Figure 1G). VOI volume: 3.4 ± 1.2.
5. For the Aorta-Peak VOI studies (AO-Peak): In the aorta a spherical peak-VOI (∼1 ml) was constructed with the peak VOI contouring tool using the inner background cylindrical VOIs as the delimiter in all frames (Figure 1H). VOI volume: 1.0 ± 0.6.
2.9. Correction of IDIFs for partial volume effect and spill-over from background
The basic input functions were derived from the activity data in the various input VOIs and converted to plasma activity by multiplying the whole blood activity with the PWR in all frames.
All basic input functions – LV-Orig, LV-New, AO-Fix, and AO-Peak with the exception of the LV-Fix – were corrected for PVE and Bg using the recovery coefficient RCß as described previously (2, 3):
where RIDIF is activity measured in VOIs in either LV or AO, CBg is background activity, and CA is the “true” activity in arterial blood, which after 30 min equals the activity in venous blood (3).
For the comparison and estimation of the influence of background activity, all input data were also corrected using the simpler calibration factor RCCF (1) without correction for spill-over from background activity:
This equation can be used with VOIs with small volumes ≤1 ml, e.g., LV-Fix placed so far from surrounding background activity (>25 mm) that spill-over from background activity mostly can be ignored. An illustrative phantom measurement has shown that activity spill-over beyond this distance is rather constant and less than 8.2% (see Supplementary Section 1.5).
After correction, basic input functions were given identifiers by adding either -ß or -CF to the geometry nomenclature—e.g., “LV-VOI-Orig-ß” and so on.
2.10. Modification of the basic IDIFs
In all the basic image curves the terminal exponentials from 40 to 60 mpi were replaced by an exponential calculated from the plasma samples using either 40–60 mpi plasma exponentials or 40–90 mpi plasma exponentials.
All image values of the IDIFs were transformed using the natural logarithm function. The slope and intercept for the resulting line from 40 to 60 mpi (“terminal exponential”) was then determined by linear regression analysis and extrapolated back to the time for the peak.
The values of this line were subtracted from the values of the entire logarithmic curve to obtain the residual curve of the initial “fast image exponentials.”
The input curves were then reconstituted in two ways:
1. The logarithmically transformed plasma curve derived from the plasma samples (40–60 mpi and 40–90 mpi) were added to the logarithmically transformed residual curve and then exponentially retransformed as in our original study (2)—“the multiplicative logarithmic” (“Log”) method with the given identifiers “LV-Orig-Pl-40-90-Log (t60), and so on.
2. Both the logarithmically transformed residual curve and the plasma curve (40–60 mpi and 40–90 mpi) were retransformed using the exponential function before being added (6)—“the additive exponential” (“Exp”) method with the given identifiers “LV-Orig-Pl-40-90-Exp (t60),” and so on). The terminology “(t60)” indicates that the resulting Ki values were obtained using the 60 mpi data point.
All input curve combinations used for calculating bone plasma clearance are presented in Supplementary Section 1.1 (Table S1).
2.11. Bone plasma clearance
Ki ml/min−1·ml−1 was calculated as the mean value of four thoracic vertebrae (Th7–Th10 or Th8–Th11). No correction for delay was made as we found the delay to be of only a few seconds, and attempts at correcting the very noisy data in the first few acquisition frames failed to make the data more consistent.
2.12. Non-linear regression (NLR) analysis
The PMOD® version 4.206 software (PMOD Technologies LLC, Switzerland) was used to perform a two-tissue compartment dynamic NLR analysis of 18F-NaF-turnover as described by Hawkins et al. (8). The exchange of 18F-NaF between the compartments—plasma, extravascular, and bone—is described by the kinetic parameters K1-k4 and the parameter for regional bone plasma clearance Ki is defined as
This NLR method was further analyzed for two values of the fractional blood volume vB, which was either fixed at 0.05 (vB-Fix) or determined by PMOD as a free-fitted (vB-Free) parameter.
2.13. Patlak analysis of dynamic studies
Assuming the efflux rate constant k4 to be negligibly small (k4 ≈ 0 min−1), the Patlak graphical analysis (9) provides a simpler alternative analysis method for measuring Ki as described by equation 6 (5):
This equation approximates a straight-line fit with Ki as the slope where CB and CPl are the respective concentrations of tracer bound in bone and freely diffusible in plasma at each time point (t). V0 is the intercept of the ordinate and represents the apparent volume of distribution.
Ki was calculated from the 60-min dynamic PET/CT scan using a bone TAC and each of the various IDIF modifications or selected semi-population input functions as described in Section 2.14. All Patlak results generated by the various input functions were compared to the basic corrected IDIF within each series and to the corresponding PMOD results.
2.14. Semi-population input functions (SPIFs) for static scan analysis
For the future analysis of static WB scans with the “best” input function as defined below, we derived an optimized three-exponential SPIF, as previously described by Blake et al. (6). The SPIF was derived from a population residual curve (PopRes) and then added to the terminal exponential derived from the plasma samples (described above under “Blood samples”), where the PopRes was derived from the corrected basic IDIFs scaled to a reference dose of 100 MBq, as used in the study by Puri et al. (1):
The residual curve represents the sum of the early fast exponentials and was derived by subtracting the terminal exponential (all data values ≥40 mpi) from the entire image-derived curve.
All residual curves were adjusted so that the times of peak count rate for all curves were coincident with the most frequent unadjusted peak time (16.5 s). The residual curves were then averaged to define the PopRes for combination with the plasma exponential to make the SPIF.
A mathematical model for the PopRes was fitted and is described in Supplementary Section 1.4.
2.15. Static scan analysis
As in the previous study (2), the static scan analysis was performed using a modified Patlak analysis with only two data points (4, 6). The first (0,V0) was obtained from either the original dynamic function or the reconstructed SPIF using the individual patient's blood samples. The second time point was obtained as the start time for the WB scan and the activity in the vertebrae (CB(t), CPl(t)) at that time using the same VOIs as for the dynamic scan but adjusted for proper alignment. The values of the static scan data were corrected for decay to the time of injection/start of the dynamic scan for comparison with the dynamic data.
The Ki values were then calculated as the slope of the line between (0,V0) and (CB(t), CPl(t)) (4–6).
2.16. Statistics
The results are presented as mean ± standard deviation (SD) but, as the majority of the datasets obtained using the different input functions were not supposed to be normally distributed, visual representation of the data are presented in box-and-whisker plots showing Max, 75% quartile, Median, 25% quartile, and Min values as well as the difference (the box) between the 75% quartile and the 25% quartile (the interquartile range (IQR)).
Differences-between-method results were evaluated using Bland–Altman plots (10) showing the mean differences between the corresponding data points and the upper and lower 95% confidence limits.
Correlations between Ki values obtained using different analysis models were calculated using Pearson's correlation coefficient but were not used as selection criteria when choosing the best agreement between the methods/parameter choices.
The percentage coefficient of variation of the PopRes curves was obtained as the ratio between the SD and the mean of the PopRes curve.
The PMOD χ2 test was used to evaluate the model fit of the input functions to the applied model.
The original and reanalyzed data were compared using a paired, two-tailed t-test.
3. Results
3.1. Reproducibility of original data
Comparison of the reanalyzed and original data (LV-Orig β series) using the original method definitions are summarized in Supplementary Section 1.2. The reanalyzed and original data (PWR, RCß, AUCs of derived input curves, NLR-Ki and Patlak-Ki results) were all comparable.
The mean inter-observer difference for the NLR-ß-Ki values were 0.0006 ± 0.0052, r = 0.91, and p = 0.72 (NS) and −0.0018 ± 0.0042. r = 0.89, and p = 0.17 (NS) for the Patlak-ß-Ki values.
3.2. VOI volumes and corresponding recovery coefficients
The VOI volumes for input VOIs, Bg, and Bone with their corresponding RCs for the basic input IDIFs are shown in Table 1. The LV-Orig VOI has the lowest RCß (0.73 ± 0.17) compared with the AO-Fix VOI's RCß (0.99 ± 0.09). In contrast, the AO-Peak VOIs had RCs of 1.73 ± 0.62. The RCCFs were generally higher as spill-over was not included.
3.3. Comparison of the basic image derived input functions
All median and range data supporting the following observational results for the uncorrected and corrected Basic IDIF AUCs can be found in the Supplementary Section 1.3 (Figure S1A,B) box plots:
Correction of the basic IDIFs with either RCß or RCCF tended to shift the AUCs toward higher AUC values (Supplementary Section 1.3, Figure S1B) but with less change in AUC values for smaller VOIs, the highest values being 377.53 ± 66.88 kBq·min and the lowest 293.57 ± 66.09 kBq min. The AO-Fix curves were almost unchanged while the AO-Peak curves showed lower values (Δ-AUC-mean: −187.1 ± 122.21 kBq min) with an unacceptably wide 95% CI of −426.6 to 52.5 kBq min.
3.4. Input functions with plasma-exponential substitutions
The effect of substituting the final image exponential in the corrected basic IDIFs with plasma exponentials is shown in the box plots (Supplementary Figure S2A–D).
IDIFs that were modified using the exponential additive method generally had a lower IQR and range for the 40–90 mpi data analysis compared with the corresponding 40–60 mpi analysis curves (Δ-IQR: −33.6–1.8; range: −7.7 to −19.7). However, the mean AUCs for the LV input functions were not significantly different from the corresponding basic input functions as shown in Supplementary Section 1.3 (Table S3). This also applied to the SPIFs reconstructed from the LV-New-ß series. The small differences in AUCs between the SPIFs and their corresponding IDIFs were not significant.
In Figure 2A, the model of our PopRes was compared with data from the observed, optimized PopRes curve. The visual fit was very good and use of the curve results in comparable Ki results, as described in the following subsections. The mean difference was −0.40 ± 2.79 with a 95% CI of −5.87 to 5.06. The AUC1800sec was 4,078 kBq s.
Figure 2. Population residual curves. (A) Modeled Pop-Res curve (red) compared with observed data (blue squares). (B) Modeled Pop-Res curve (red) compared with Pop-Res model (green) constructed using parameters published by Puri et al. (1).
The mathematical best-fit curve for our PopRes is presented in Supplementary Section 1.4.
3.5. Ki results using the various input functions
3.5.1. Patlak analysis
All Patlak results for the various corrected, basic IDIFs were compared as illustrated in Figure 3A–D. The box plot in Figure 3A shows the distribution of data for all the input series and their values in Table 2.
Figure 3. Cross comparison of Patlak-Ki-results using the basic corrected IDIFs. (A) Box-plot shoving the distribution of data within the various series. (B–D) Bland-Altman plots of the differences between the basic LV-Orig-ß input function the basic input functions most comparable to LV-Orig-ß. “Upper” and “Lower” refers to 95%-confidence limits.
Results for the LV-New-ß and LV-Fix-CF (Figure 3A) show comparable median and mean values, whereas the LV-Orig-ß are a little lower and the AO-Fix-ß and AO-Fix-CF a little higher, but the observed differences are not significantly different from zero, as shown in the corresponding Bland–Altman plots (Figure 3B–D). The data in Table 3A show the LV-New-ß to have the smallest confidence interval and the least mean difference and variance compared with the LV-Orig-ß function. However, apart from the AO-Peak results, which show an unacceptable large variance (Figure 3A), the differences and variation were generally small.
In order to select the derived input function with the least difference and variance compared with the corresponding corrected basic input function, the Patlak data were compared serially as shown in Figure 4 and in Supplementary Section 2.1 (Figure S3). In all series, the Ki results at t = 60 mpi substitution with plasma exponentials 40–90 mpi using the exponential method showed the least difference and variation compared with the basic corrected IDIF, with the smallest values observed for the LV-New-ß series compared with the LV-Original method (Table 3B).
Figure 4. Patlak analysis. Serial comparison of the various modifications of the basic IDIFs. Box-plots (A–D) shows the distribution of values using the various input functions.
For all input functions, the mean differences between the use of the exponential (additive) method and the logarithmic (multiplicative) method are found to be very small, but the smallest variance is observed for the use of the exponential method (Table 3B).
3.5.2. Non-linear regression (NLR) analysis
For completeness, all NLR-fitted kinetic parameter results for the differing input analysis methods are found in Supplementary Section 2.2 (Table S4A,B).
The mean Ki values obtained for the various corrected, basic IDIFs were determined using either a constricted vB-Fix of 0.05 (Supplementary Table S4A) or calculated as part of the PMOD parametric fitting. The Ki values using free-fitted vB resulted in higher Ki values for all five input functions (Supplementary Table S4B).
The distribution of Ki values using the corrected five basic input functions are shown in Figure 5A,B. The results are comparable except for the AO-Peak-IDIF, which shows the largest difference and variance and hence was excluded from further analysis.
Figure 5. NLR-Ki-results using fixed vB vs. free-fitted vB. (A) NLR-corrected, vB = 0.05. (B) NLR-corrected, free-fitted vB.
The differences between Ki values using the various input functions and LV-Orig-ß (Supplementary Section 2.2, Figures S4A–F and Supplementary Table S5) were not significant but the confidence intervals were rather variable with the smallest variance found for the LV-New-ß with vB-free-fit (Supplementary Figure S4D). The LV-Fix input function showed the largest confidence interval in both modes.
The vB-Free had a mean value of 0.01 ± 0.008 across the input modes.
3.5.3. Comparison of Patlak results with NLR results
Based on the prior analyses (Sections 3.5.1, 3.5.2), the NLR-Ki results for the LV-Orig-ß series, LV-New-ß series, and AO-Fix-ß series were chosen for comparison with the corresponding Patlak-Ki results, with the distribution of data shown in Figure 6A–H and the corresponding Bland–Altman plots in Supplementary Section 2.3 (Figure S5A–H). Quantitative results are presented in Table 5 and the mean differences in Supplementary Section 2.3 (Table S6). The smallest difference and confidence interval are seen with the LV-New-ß series and, as such, this series is identified as our choice of optimum analysis method. All observed variances are comparable with Δ-Ki ≤ 0.013 ml/min· min−1 (Supplementary Section 2.3, Table S6).
Figure 6. Distribution of optimized NLR-Ki-values compared with the corresponding Patlak-values. (A,B) LV-Orig-ß. (C,D) LV-New-ß. (E,F) LV-Fix-CF. (G,H) AO-Fix-ß. All analyses with either fixed vB (left panels) or free-fitted vB (right panels).
Table 5. Comparison of NLR- (fixed and free vB) and Patlak- Ki-analysis for selected basic input functions and their derivatives.
3.5.4. Comparison of dynamic results with results using semi-population functions
For this analysis, we exclusively used data for the LV-New-ß-Input series (optimum analysis) and the semi-population functions derived from these (Figure 2A). In Figure 7A,B, Patlak-analysis plots from the same patient are shown using the LV-New-Pl-40-90-Exp and the corresponding SP-Pl-40-90-Exp input functions. In both plots, data from the static scan are shown as a red square lying close to the regression line. The distribution of data using dynamic input functions are shown in the box plots in Figure 7C,D with quantitative differences shown in the Bland–Altman plots of Figure 7E,F, for the Log (multiplicative) and Exp (additive) analysis methods, respectively. The observed differences within the dynamic series were smallest for the Exp series, as seen in Figure 7F.
Figure 7. Examples of static scan data (red square) plotted in combination with dynamic data (black triangles) for the same patient using (A): Patl-LV-New Pl-40-90-Exp input function and (B): SP-New-Pl-40-90-Exp input function. (C,D): Comparison of Ki-results obtained by the dynamic input functions and SP-input functions using the logarithmically substituted and the exponentially substituted IDIFs, respectively. (E,F): Bland-Altman analysis of the corresponding logarithmically and exponentially substituted IDIFs. “Upper” and “Lower” refers to 95%-Confidence limits.
3.5.5. Static scan analysis using two-point Patlak analysis
The results for the static scan analyses using the two-point Patlak analysis with either the LV-New-Pl-40-90 dynamic input functions or the derived SPIF are shown in Figure 8A–D and Table 6. As for the dynamic input functions, the series modified using exponential substitution of the final exponentials showed the lowest mean Ki-difference from zero with the 95% confidence interval of −0.0027 to 0.0029.
Figure 8. Static WB-scan analysis. Comparison of Ki-results obtained by dynamic LV-New-Pl-40-60 functions and the corresponding derived SP-input functions using Patlak Two-Point analysis. (A,B) Distribution of data and differences using logarithmically substituted IDIFs. (C,D) Distribution of data and differences using exponentially substituted IDIFs. “Upper” and “Lower” refers to 95%-confidence limits.
Compared to the dynamic results presented above, the Ki values measured using the static scan analysis were 0.0019 ± 0.0017 ml/min ml−1 higher than those for the corresponding dynamic scan analysis. However, all differences were well within the 95% confidence limits (−0.0015 to 0.0053 for both the logarithmic and exponential series) and thus were not significantly different from zero.
4. Discussion
4.1. Comparison with original results
As a test of inter-observer reproducibility, all original image data (2) were reanalyzed using the original method parameters (Section 3.1). The observed slight difference (NS) in the plasma-to-whole-blood ratio (PWR-reanalyzed = 1.16 ± 0.02. vs. PWR-Orig = 1.17 ± 0.03) lies within the range of values presented in Figure 3 in Vrist et al. (2) and is most likely due to the more selective exclusion of plasma outliers in the reanalysis in the present study.
The RCß correction factor used in the present study (0.72 ± 0.17) was comparable to that of the original study (0.69 ± 0.15).
The AUC of the basic input curves in this study is 374 ± 59 kBq min, which is comparable with the 353 ± 59 kBq min result from the original study, with the resulting Ki values of the NLR and Patlak analyses differing slightly, but with an acceptable inter-observer difference (NLR results: 1.3%; Patlak results: 6.3%).
Reproducible values between the original data analysis of this and the original study (2) are a prerequisite for being able to attribute subsequently observed differences in kinetic parameters as being due to this study's use of “optimized” input functions while keeping all other parameters (PWR, Bg-VOIs, and Bone-VOIs) almost constant in all series.
This assumption is supported by previous repeatability studies on the use of 18F-NaF-PET to evaluate SUV as a measure for uptake in cancer patients (11, 12), as well as 18F-NaF-PET being previously used successfully for studies on changes in Ki values before and after treatment of osteoporosis (4, 7, 13).
4.2. Effects on basic input functions by VOI size and correction factors
Comparing the VOI volumes and their corresponding RCs, it is obvious there is a correlation between these parameters, as the RCs are derived from the activity in the VOIs and compared with the actual blood sample activities (Table 1): the smaller the VOI, the closer to 1 the correction factor, except for the AO-Peak VOI in which the factors are much higher than 1. This correlation has been described previously in cardiac FDG-PET studies (14) and can be ascribed to a combination of PVE and an increasing mean activity (kBq/ml) in the VOIs as the smaller volumes are placed, more selectively, over areas with the highest activity.
The RCß corrects the activities measured in the input VOI for both PVE and Bg (3), whereas the RCCF only corrects for PVE. The RCß (0.99 ± 0.09) and RCCF (0.99 ± 0.09) are identical for the AO-Fix-ß VOIs, indicating the background in and around the aortic wall is negligible, at least at this level of the thoracic aorta. These factors are comparable to the coefficients reported in a study by Puri et al. (15). However, for the AO-Peak VOIs, the correction factors are unrealistically high (1.73 ± 0.63) and result in a much larger variance in the derived Ki results when compared with the other corrected input functions.
As the correction factors were primarily lower than 1, the general effect on the uncorrected input functions was to increase the AUCs as shown in Supplementary Figure S1. The closer the correction factor is to 1, as in the AO-Fix curves, the lesser the shift in values.
The exception to this was the AO-Peak curves which, due to the correction factors of 1.73 ± 0.63, resulted in smaller AUCs and larger variance. The most likely reason for the poor AO-Peak VOI performance is due to its inherent poor counting statistics caused by the small VOI volume of only a few voxels in combination with the amount of activity injected and the short time resolution in the initial dynamic acquisition bins.
4.3. Effect of substituting the final exponentials of basic input functions with plasma exponentials
The general effect of substituting the final image exponentials with plasma exponentials was smaller IQRs, as seen in the box plots in Supplementary Figure S2. However, the AUC mean differences were not significantly different across the series of input functions (Supplementary Section 1.3, Table S3).
The substitution methods producing AUC results (kBq/ml min) closest to the basic corrected IDIFs were the LV-New-Pl-40-90-Log (Δ-Mean: −1.17, 95% Cl −17.7 to 15.3) and the LV-New-40-90-Exp (Δ-Mean: 3.2, 95% Cl −9.7 to 16.2). For comparison, the corresponding values for the AO-Fix-40-90-Exp were much higher (Δ-Mean: −11.0, 95% Cl −44.4 to 22.5; p = 0.04).
4.4. Comparison of Ki results using the various analysis methods
4.4.1. Patlak analysis
The Patlak analysis was applied to all basic IDIFs in order to find the IDIFs with the least variances of Ki values, after which they were compared with the results using the basic LV-Orig-ß-IDIF. Our choice for comparing to LV-Orig-ß-IDIF was not because we believe our original analysis method to be “a reference standard” but simply to assess the performance of this method compared with the others (“comparative reference”).
From the box plot in Figure 3A, the LV-New-ß, LV-Fix-ß, and AO-Fix-ß are identified as all being qualified methods. From the Bland–Altman plots in Figure 3B–D and the data in Table 3, it is shown that no IDIF showed Ki results significantly different from those of the basic LV-Orig-ß-IDIF but the LV-New-ß IDIF resulted in the least difference when compared with the LV-Orig-ß-IDIF (0.0009 ± 0.0042 ml/min ml−1).
The modified IDIF with the least variance and least difference of Ki values compared with the basic corrected IDIF was selected from a serial comparison of the Ki values obtained by the modified IDIFs as shown in Figure 4, Tables 2 and 3, and in Supplementary Section 2.1 (Figure S3).
Of all series, the LV-New-IDIFs with plasma substitution 40-90-Exp showed the least variance, with the AO series showing the highest variances, but the mean differences between the modified IDIFs and the corresponding basic IDIFs were not significantly different from zero.
The Ki results of the input functions obtained using small VOIs tend to be a little higher than the LV-Orig series but this may simply be due to these functions having smaller AUCs, which, in turn, results in a higher ratio between activity uptake in the bones compared with the lower activity in the plasma. Generally, as can be seen in Figure 3A and Tables 2 and 3, the variations in Ki across the various basic and derived input functions are small, but even so, the Ki mean difference of 0.0023 ± 0.0046 ml/min ml−1 for the Patl-AO-Fix-ß-IDIF compared with Patl-LV-Orig-ß was significantly higher than that of Patl-LV-New-ß (p = 0.01).
Overall, our data are comparable with results reported in other studies. In the original study by Hawkins et al. (8) of thoracic vertebrae in normal individuals using image-derived input functions, the mean Ki value was 0.093 ± 0.0071. In a study by Frost et al. (16) of hemodialysis patients suspected for adynamic bone disease and osteoporosis, the mean Ki results in the lumbar spine were 0.028 ± 0.012 and 0.027 ± 0.005, respectively. Thus, for example, the mean basic Patlak results for our patients (0.0302 ± 0.0071 to 0.0325 ± 0.0070) are well within this rather wide range.
4.4.2. Non-linear regression (NLR) analysis
The NLR analysis based on the Hawkins model has often been used as a reference standard for the less complicated Patlak analysis (2, 8, 17, 18). In our original study, we used the PMOD Software® with the default vB of 0.05, which is a value previously reported by Messa et al. (19). However, it has been questioned why we did not allow the vB to be fitted freely; therefore, we have examined the possible differences using both NLR-Fix vB = 0.05 and NLR-Free vB as shown in Figure 5 and Table 4 and in Supplementary Section 2.2 (Figure S4 and Table S5).
The general effect on Ki values using the free-fitted vB tended to be 3.7%–8.0% higher values (5.3% on average). However, none of these differences were significantly different from zero, but the AO-Peak-ß showed a very high SD of 0.35 ml/min ml−1 with a CV of approximately 60% compared with 27% for the other IDIFs. Based on this observation, we excluded this variant from further studies, as it would probably be too insensitive for the detection of small changes in Ki.
The NLR-Free-fit vB analysis showed a smaller range of Ki values, especially for the LV-input series (Supplementary Section 2.2, Table S4), and this mode will be chosen for future dynamic study analyses. This is in accordance with the results reported by Puri et al. (20), showing the CV for Ki to be approximately 15% whereas the K1–k4 parameters had CVs of at least 30%, which is likely caused by the NLR method being sensitive to statistical noise in the input data.
In this study, we found the averaged free-fitted value for vB to be 0.01 ± 0.008. This is about one-fifth of the value of 0.05 originally reported by Messa et al. (19) and five times the value of 0.002 used in more recent studies by Puri et al. (21). We have no explanation for these differences, but differences between bone regions and disease populations seem plausible.
4.4.3. Comparison of Patlak results with NLR results
The variances of the NLR results were much higher than the variances of the Patlak results. As discussed in the previous subsection, this is probably due to the NLR method having a greater sensitivity to statistical noise in input data, whereas the Patlak analysis appears to be very robust, as indicated by the relatively small differences observed across the various input series (Table 3).
It is well known that Patlak results are generally lower than NLR results, which have previously been reported to be −7% by Brenner et al. (17), −13% by Installé et al. (22), and −23.7% by Puri et al. (23).
In the present study, the differences were in the range of −17.3% to −26.1% for the NLR-Fix vB and −22.8% to −29.2% for the NLR-Free fit vB, with the smallest difference being for the LV-New-ß series in both comparisons (Supplementary Section 2.3, Table S6). In our original study, we found a difference of −17.4 ± 10% (2), so the optimized data are within the previously reported range.
Consequently, our aim was to select the Patlak series showing the least variance and difference when compared with the PMOD-NLR analysis, both with fixed vB and free-fitted vB as shown in Figure 6 and Supplementary Section 2.3 (Table S6, Figure S5).
4.4.4. Comparison of dynamic results with results using semi-population functions
Based on the above results and discussion, we chose to derive our optimized population residual functions from the LV-New-Pl-40-90-Exp input function, which was least variable and least different from the corresponding basic input data. This population residual was used to construct the semi-population curves to be used as input functions for the static scan analysis as originally described by Frost et al. (4) and further used in, for example, studies by Vrist et al. (2) and Blake et al. (6).
As shown in Figure 2A, the parameters for our model CKD-MBD population residual curve result in a good fit with the curve derived as an average of the residuals obtained from our analyses. Compared with an osteoporosis population model curve calculated using the parameters recently published by Puri et al. (1) (Figure 2B), there are some differences in peak height and full-width-half-maximum (FWHM) of the peaks, even though our data were normalized to the injected dose of 100 MBq as used in the study by Puri et al. (1). As a result, the AUC1800s of our curve (4,078 kBq s) was approximately 12% lower than that of 4,560 kBq s reported by Puri et al. (1). The reason for the observed differences in the curve forms is probably due to “delay and dispersion,” which is to be expected as our input function were obtained from the left ventricle, whereas Puri et al. (1) used the abdominal aorta as input for the study of bone metabolism in lumbar vertebrae. In comparison, our aorta input curves at the level of Th7 showed no significant difference in shape or peak height, with a delay of only 3 s compared with those of the left ventricle.
However, apart from the technical issues regarding the curve differences, the interesting question of whether these differences are caused by differences in bone metabolism in different patient populations (1, 7) cannot be conclusively answered using the data from this study or our original study, as we were unable to use identical VOI definitions and placement in the arterial system. The data show that differences in the placement of input VOIs matters but, despite this, the AUCs of the two different population residual curves only differ by approximately 12% without corrections. Thus, the use of an “universal input model” for various anatomical regions and/or patient groups may be feasible but should be done with great care and possibly with the inclusion of a relevant correction for delay and dispersion (24) in order for the derived results to be comparable between patient groups. For studies of serial changes in bone plasma clearance within the same patient, this is probably less important.
The resulting Ki values for the semi-population functions showed no significant difference for either the Log series or the Exp series, but the variance and 95% confidence interval for the differences were much smaller for the Exp series. Thus, this input function derivation will be chosen for future studies as this should make the detection of smaller differences possible when comparing serial measurements of Ki, such as previously described by Frost et al. (4) before and after treatment with teriparatide, compared to use of a semi-population function with a wider variance and confidence interval.
4.4.5. Static scan analysis
In this study the LV-New input functions showed no significant difference between Ki values obtained using the dynamic input function or the corresponding SPIF-Exp series, but a little higher value was obtained using the SPIF-Log series (Table 6). The results are comparable to the original dynamic results but with a little lower Ki value compared to our original static scan Ki result of 0.0395 ± 0.011 ml/min ml−1. This may be explained by the AUC of our original SPIF (331 ± 70 kBq min−1) being on average 9% smaller than our optimized AUC of 362 ± 47 kBq min, thus resulting in a higher Ki value.
5. Conclusions
We have established that our original logarithmic/multiplicative method is valid, producing quantitative results comparable to the exponential/additive method, but with the results having a little more variance.
We have shown that recovery coefficients/correction factors are correlated to the volume of input VOIs and their location and that the correct choice of correction depends on the applied analysis geometry.
Within limits, the corrected input functions have comparable AUCs and yield correspondingly comparable Ki results but with differing variances.
The use of a blood sample taken at 90 mpi, corresponding to the end of the entire 18F-NaF-PET study protocol after the static data acquisition and not just at the end of the dynamic data collection, improves the fit between the original IDIF and the IDIF with the final exponential substituted by the plasma exponential using blood samples between 40 and 90 mpi.
Even though differences between the Ki results obtained using the various input functions are small and comparable, we recommend the use of the same analysis implementation technique for future comparative studies due to the possibility of the differing variances making serial changes in Ki more difficult to detect for the use of non-optimized analysis parameters/methods.
A comparison of two models for population residual curves for differing patient populations showed significant differences in peak shape (peak height and FWHM), but a difference of approximately 12% only in total AUC, which indicates that a comparison of data may be possible, but the input curves should ideally be obtained with comparable anatomical input VOIs or at least corrected for differences in delay and dispersion. To answer the question of whether the input curves vary due to changes in bone metabolism and/or between various patient population studies would require a comparative study utilizing the same analysis technique in order to avoid methodological pitfalls.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Regional Committee on Biomedical Research Ethics (case number: 1-10-72-340-16) and the Danish Health and Medicines Authority (EudraCT number: 2016-005160-34). The patients/participants provided written informed consent to participate in this study.
Author contributions
JT and CF were responsible for the study design, image and data analysis, and drafting the manuscript. MV and JB were responsible for the recruitment of patients and original study data. JT, MV, JB, and CF were responsible for the critical revision of and intellectual contribution to the final article. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnume.2023.1235800/full#supplementary-material
Abbreviations
AO, aorta (thoracalis); Bg, spill-over of activity from background; CI, confidence Interval; FWHM, full-width-half-maximum—width of a peak at 50% of its maximum value; IDIF, image-derived input function; IQR, interquartile range (difference between 75% quartile and 25% quartile); LV, left ventricle; NS, non-significant; PWR, plasma to whole blood ratio; RC, recovery coefficient; ROI, region of interest; vB-Fix, fraction of blood volume (in bone) with a fixed value; vB-free, fraction of blood volume (in bone) with a free-fitted value; VOI, volume of interest.
References
1. Puri T, Frost ML, Moore AEB, Cook GJR, Blake GM. Input function and modeling for determining bone metabolic flux [18F] sodium fluoride PET imaging: a step-by-step guide. Med Phys. (2023) 50(4):2017–88. doi: 10.1002/mp.16125
2. Vrist MH, Bech JN, Lauridsen TG, Theil J. Comparison of [18F] NaF PET/CT dynamic analysis methods and a static analysis method including derivation of a semi-population input function for site-specific measurements of bone formation in a population with chronic kidney disease-mineral and bone disorder. EJNMMI Res. (2021) 11:117–30. doi: 10.1186/s13550-021-00859-7
3. Cook GJR, Lodge MA, Marsden PK, Dynes A, Fogelman I. Non-invasive assessment of skeletal kinetics using fluorine-18 fluoride positron emission tomography: evaluation of image and population-derived arterial input functions. Eur J Nucl Med. (1999) 26(11):1424–9. doi: 10.1007/s002590050474
4. Frost ML, Siddique M, Blake GM, Moore AEB, Schleyer PJ, Dunn JT, et al. Differential effects of teriparatide on regional bone formation using (18)F-fluoride positron emission tomography. J Bone Miner Res. (2011) 26(5):1002–11. doi: 10.1002/jbmr.305
5. Siddique M, Blake GM, Frost ML, Moore AEB, Puri T, Marsden PK, et al. Estimation of regional bone metabolism from whole-body 18F-fluoride PET static images. Eur J Nucl Med Mol Imaging. (2012) 39(2):337–43. doi: 10.1007/s00259-011-1966-y
6. Blake GM, Siddique M, Puri T, Frost ML, Moore AE, Cook GJR, et al. A semipopulation input function for quantifying static and dynamic 18F-fluoride PET scans. Nucl Med Commun. (2012) 33:881–8. doi: 10.1097/MNM.0b013e3283550275
7. Blake GM, Siddique M, Frost ML, Moore AE, Fogelman I. Radionuclide studies of bone metabolism: do bone uptake and bone plasma clearance provide equivalent measurements of bone turnover? Bone. (2011) 49(3):537–42. doi: 10.1016/j.bone.2011.05.031
8. Hawkins RA, Choi Y, Huang SC, Hoh CK, Dahlbom M, Schiepers C, et al. Evaluation of the skeletal kinetics of fluorine-18-fluoride ion with PET. J Nucl Med. (1992) 33(5):633–42.1569473
9. Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood to brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab. (1983) 3(1):1–7. doi: 10.1038/jcbfm.1983.1
10. Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. (1986) 1:307–10. doi: 10.1016/S0140-6736(86)90837-8
11. Kurdziel KA, Shih JH, Apolo AB, Lindenberg L, Mena E, McKinney YY, et al. The kinetics and reproducibility of 18F-sodium fluoride for oncology using current PET camera technology. J Nucl Med. (2012) 53:1175–84. doi: 10.2967/jnumed.111.100883
12. Lin C, Bradshaw T, Perk T, Harmon S, Eickhoff J, Jallow N, et al. Repeatability of quantitative 18F-NaF PET: a multicenter study. J Nucl Med. (2016) 57:1872–9. doi: 10.2967/jnumed.116.177295
13. Frost ML, Moore AE, Siddique M, Blake GM, Laurent D, Borah B, et al. 18F-fluoride PET as a noninvasive imaging biomarker for determining treatment efficacy of bone active agents at the hip: a prospective, randomized, controlled clinical study. J Bone Miner Res. (2013) 28(6):1337–47. doi: 10.1002/jbmr.1862
14. Lin K-P, Huang S-C, Choi Y, Brunken RC, Schelbert HR, Phelps ME. Correction of spillover radioactivities for estimation of the blood time-activity curve from the imaged LV chamber in cardiac dynamic FDG PET studies. Phys Med Biol. (1995) 40:629–42. doi: 10.1088/0031-9155/40/4/009
15. Puri T, Blake GM, Siddique M, Frost ML, Cook GJR, Marsden PK, et al. Validation of new image-derived arterial input functions at the aorta using 18F-fluoride positron emission tomography. Nucl Med Commun. (2011) 32(6):486–95. doi: 10.1097/MNM.0b013e3283452918
16. Frost ML, Compston JE, Goldsmith D, Moore AE, Blake GM, Siddique M, et al. 18F-fluoride positron emission tomography measurements of regional bone formation in hemodialysis patients with suspected adynamic bone disease. Calcif Tissue Int. (2013) 93(5):436–47. doi: 10.1007/s00223-013-9778-7
17. Brenner W, Vernon C, Muzi M, Mankoff DA, Link JM, Conrad EU, et al. Comparison of different quantitative approaches to 18F-fluoride PET scans. J Nucl Med. (2004) 45(9):1493–500.15347716
18. Assiri R, Knapp K, Fulford J, Chen J. Correlation of the quantitative methods for the measurement of bone uptake and plasma clearance of 18F-NaF using positron emission tomography. Systematic review and meta-analysis. Eur J Radiol. (2022) 146-110081:1–15. doi: 10.1016/j.ejrad.2021.110081
19. Messa C, Goodman WG, Hoh CK, Choi Y, Nissenson AR, Salusky IB, et al. Bone metabolic activity measured with positron emission tomography and [18F]fluoride ion in renal osteodystrophy: correlation with bone histomorphometry. J Clin Endocrinol Metab. (1993) 77(4):949–55. doi: 10.1210/jcem.77.4.8408470
20. Puri T, Frost ML, Cook GJ, Blake GM. [18f] sodium fluoride PET kinetic parameters in bone imaging. Tomography. (2021) 7(4):843–54. doi: 10.3390/tomography7040071
21. Puri T, Siddique MM, Frost ML, Moore AEB, Blake GM. A short dynamic scan method of measuring bone metabolic flux using [18F]NaF PET. Tomography. (2021) 7(4):623–35. doi: 10.3390/tomography7040053
22. Installé J, Nzeusseu A, Bol A, Depresseux G, Devogelaer J-P, Lonneux M. 18F-fluoride PET for monitoring therapeutic response in Paget’s disease of bone. J Nucl Med. (2005) 46(10):1650–8.
23. Puri T, Blake GM, Frost ML, Siddique M, Moore AEB, Marsden PK, et al. Comparison of six quantitative methods for the measurement of bone turnover at the hip and lumbar spine using 18F-fluoride PET-CT. Nucl Med Commun. (2012) 33(6):597–606. doi: 10.1097/MNM.0b013e3283512adb
Keywords: [18F] sodium fluoride, PET imaging, bone metabolism, kinetic analysis, modeling, input function, comparative analysis, CKD-MBD
Citation: Theil J, Vrist MH, Bech JN and Fynbo CA (2023) A systematic evaluation of five different image-derived input functions for the clinical implementation of 18F-NaF bone PET/CT in patients with chronic kidney disease–mineral and bone disorder. Front. Nucl. Med. 3:1235800. doi: 10.3389/fnume.2023.1235800
Received: 6 June 2023; Accepted: 10 July 2023;
Published: 28 July 2023.
Edited by:
Mario Petretta, IRCCS SYNLAB SDN, ItalyReviewed by:
Tommaso Volpi, Yale University, United StatesMercy Iyabode Akerele, University of Florida, United States
© 2023 Theil, Vrist, Bech and Fynbo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jørn Theil, anRoZWlsQGRhZGxuZXQuZGs=