Abstract
Several methods are in use for analyzing 11C-Pittsburgh compound-B (11C-PiB) data. The objective of this study was to identify the method of choice for measuring longitudinal changes in specific 11C-PiB binding. Methods: Dynamic 90-min 11C-PiB baseline and follow-up scans (interval, 30 ± 5 mo) were obtained for 7 Alzheimer disease (AD) patients, 11 patients with mild cognitive impairment (MCI), and 11 healthy controls. Parametric images were generated using reference parametric mapping (RPM2), reference Logan values, and standardized uptake value volume ratios (SUVr), the latter for intervals between 60 and 90 (SUVr60–90) and 40 and 60 (SUVr40–60) minutes after injection. In all analyses, cerebellar gray matter was used as a reference region. A global cortical volume of interest was defined using a probability map–based template. Percentage change between baseline and follow-up was derived for all analytic methods. Results: SUVr60–90 and SUVr40–60 overestimated binding with 13% and 10%, respectively, compared with RPM2. Reference Logan values were on average 6% lower than RPM2. Both SUVr measures showed high intersubject variability. Over time, R1, the delivery of tracer to the cortex relative to that to the cerebellum, decreased in AD patients (P < 0.05) but not in MCI patients and controls. Simulations showed that SUVr, but not RPM2 and reference Logan values, was highly dependent on uptake period and that changes in SUVr over time were sensitive to changes in flow. Conclusion: To reliably assess amyloid binding over time—for example, in drug intervention studies—it is essential to use fully quantitative methods for data acquisition and analysis.
- Alzheimer disease
- positron emission tomography
- 11C-Pittsburgh compound-B
- receptor parametric mapping
- reference Logan
- SUVr
Neuropathologically, Alzheimer disease (AD) is characterized by the presence of senile plaques that consist mainly of amyloid-β (Aβ) (1). Amyloid burden can be measured in vivo using PET and the ligand N-methyl-11C-2-4(4′-methylaminophenyl)-6-hydroxy-benzothiazole, also known as Pittsburgh compound-B (PiB) (2). Previous studies using 11C-PiB have shown high diagnostic accuracy for the detection of AD (2,3). Longitudinal studies, however, have provided inconsistent findings with either no (4–6) or modest (7–9) changes in 11C-PiB binding over time in AD patients. Apart from other methodologic considerations, one possible explanation for these inconsistent results could be the method used to quantify specific 11C-PiB binding.
In a previous study, receptor parametric mapping (RPM2, a basis function implementation of the simplified reference tissue model) (10) was identified as the best parametric method for analyzing 11C-PiB data (11). Compared with other methods, RPM2 was least sensitive to noise and showed the highest contrast. The difference in quantitative performance between RPM2 and reference Logan (12) was small, but RPM2 showed slightly better image quality. SUVr (standardized uptake value ratio; cortical uptake normalized to cerebellar gray matter uptake), however, showed higher variability and poorer image quality.
In most clinical 11C-PiB studies, SUVr has been used to measure amyloid load (2,6,13). This is understandable because SUVr has numerous advantages, such as computational simplicity, shorter scan duration, and less vulnerability to patient movement. Nevertheless, this method is sensitive to differences in both washin and washout of the tracer between subjects (14). Although SUVr may be acceptable for diagnostic purposes—for example, PiB-positive versus PiB-negative—more accurate quantification methods may be needed for longitudinal studies aimed to measure changes in 11C-PiB binding over time.
Therefore, the aim of this longitudinal study was to compare the most frequently used methods for the analysis of 11C-PiB scans (RPM2, reference Logan, and SUVr) in relation to changes in binding over time.
MATERIALS AND METHODS
Subjects
Data for 7 AD patients, 11 patients with mild cognitive impairment (MCI), and 11 healthy controls were used (5). All subjects underwent baseline and follow-up 11C-PiB scans, with an interval of 30 ± 5 mo (range, 24–48 mo). Subjects underwent 11C-PiB scans within an interval of 2 ± 1 (at baseline) and 1 ± 1 (at follow-up) months of the clinical evaluation. Written informed consent was obtained from all subjects after a complete written and verbal description of the study. The study was approved by the Medical Ethics Review Committee of the VU University Medical Center.
PET
PET scans were obtained using an ECAT EXACT HR+ scanner (Siemens/CTI), equipped with a neuroinsert to reduce the contribution of scattered photons from outside the field of view of the scanner. This scanner enables the acquisition of 63 transaxial planes over a 15.5-cm axial field of view, thus allowing the whole brain to be imaged in 1 bed position. The properties of this scanner have been reported elsewhere (15). All subjects received a venous cannula for tracer injection. First, a 10-min transmission scan was obtained in 2-dimensional acquisition mode using 3 retractable rotating 68Ge line sources. This scan was used to correct the subsequent emission scan for tissue attenuation. Next, a dynamic emission scan in 3-dimensional acquisition mode was started simultaneously with the intravenous injection of 351 ± 82 (baseline) or 377 ± 91 (follow-up) MBq of 11C-PiB, with specific activities of 41 ± 22 and 88 ± 40 GBq·μmol−1, respectively. 11C-PiB was injected using an infusion pump (Med-Rad; Beek) at a rate of 0.8 mL·s−1, followed by a flush of 42 mL of saline at 2.0 mL·s−1. The emission scan consisted of 23 frames with progressive increases in frame duration (1 × 15, 3 × 5, 3 × 10, 2 × 30, 3 × 60, 2 × 150, 2 × 300, and 7 × 600 s) for a total scan duration of 90 min. Patient motion was restricted with a head holder and monitored by checking the position of the head using laser beams.
MR Imaging
All subjects underwent a structural MR imaging scan using a 1.5-T Sonata scanner (Siemens Medical Solutions) at baseline and at follow-up (mean interval between PET and MR imaging, 2 ± 1 and 1 ± 1 mo, respectively). The scan protocol included a coronal T1-weighted 3-dimensional MPRAGE (magnetization-prepared rapid acquisition gradient echo) (slice thickness, 1.5 mm; 160 slices; matrix size, 256 × 256; voxel size, 1 × 1×1.5 mm3; echo time, 3.97 ms; repetition time, 2,700 ms; inversion time, 950 ms; flip angle, 8°), which was used for coregistration, segmentation, and region-of-interest (ROI) definition.
Image and Data Analysis
All PET sinograms were corrected for dead time, tissue attenuation using the transmission scan, decay, scatter, and randoms and were reconstructed using a standard filtered backprojection algorithm and a Hanning filter with a cutoff at 0.5 times the Nyquist frequency. A zoom factor of 2 and a matrix size of 256 × 256 × 63 were used, resulting in a voxel size of 1.2 × 1.2 × 2.4 mm3 and a spatial resolution of approximately 7 mm full width at half maximum at the center of the field of view. MR images were aligned to corresponding PET images using a mutual information algorithm. Data were further analyzed using PVElab, a software package that uses a probability map based on 35 delineated ROIs that have been validated previously (16). For the evaluation of the different analytic methods, a global cortical ROI was used. This ROI was composed from the volume-weighted average of the orbital frontal, medial inferior frontal, superior frontal, parietal, superior temporal, medial inferior temporal, and entorhinal cortices together with the hippocampus and posterior cingulate.
Kinetic Analysis
Data were analyzed on a voxel-by-voxel level using RPM2 (17,18), reference Logan values (12), and SUVrs (19). Cerebellar gray matter was used as reference tissue because of its low levels of fibrillar amyloid in AD patients (20).
We applied a simplified reference tissue model with basis lookup function (21), henceforth referred to as RPM2 (11), to the full 90-min dynamic PET data (10). The outcome measure of RPM2, BPND (nondisplaceable binding potential), is a quantitative measure of specific binding. For RPM2, the dynamic scan was first processed using RPM/simplified reference tissue model. This first step provided parametric images of R1 (the delivery of tracer to the cortex [K1] relative to that to the cerebellum [K1′]), BPND, and k2′. Next, the median value of k2′ was determined and subsequently fixed in the second run of RPM2 processing of the dynamic scan. Consequently, k2′ is fixed on an individual scan basis. The parametrically obtained BPND reflects the concentration of specifically bound tracer relative to that of free and nonspecifically bound tracer in tissue under equilibrium conditions. Furthermore, parametric maps of relative delivery (R1) were also generated using RPM2. Reference Logan is based on integration of the differential model equations for target and reference regions. The outcome measure DVR (distribution volume ratio) represents the ratio of distribution volumes of the cortex and cerebellum. For reference Logan, we used the implementation published by Logan et al. (12). In this implementation, it is not required to fix k2′, and DVR is derived directly as the slope of the linear part of the graphical plot. More information can be found in Yaqub et al. (11) where we described and evaluated the performance of the various parametric methods in detail. SUVr60–90 and SUVr40–60 are the ratios of tissue concentrations in the cortex and cerebellum, measured in the time frame from 60 to 90 min and from 40 to 60 min after injection, respectively. The global cortical ROI was projected onto the various parametric images. For the present comparison, results obtained using RPM2 were expressed as BPND + 1, which corresponds to the outcome measures obtained using reference Logan and SUVr. Percentage changes over time within methods for all groups were calculated using percentage change (%) = 100 × (follow-up value – baseline value)/(baseline value). Next, relative differences between methods for both baseline and follow-up conditions were calculated using relative difference (%) = 100 × (method A – method B)/method B.
Simulations
Simulations were performed to assess effects of flow variations on accuracy of 11C-PiB binding parameters. Parameters, derived from clinical studies (22), were used in combination with a typical plasma input function. For the reference region the following parameters were used: blood volume fraction VB (proportion of tissue volume occupied by intravascular blood) = 0.05, together with K1 = 0.32 mL·cm−3·min−1, k2 = 0.16 min−1, k3 = 0.025 min−1, k4 = 0.033 min−1, BPND = 0.76 (=k3/k4), and VT (volume of distribution) = 3.5 (=K1/k2·(1 + k3/k4)). Parameters for a typical AD region were set at VB = 0.05, K1 = 0.32 mL·cm−3·min−1, k2 = 0.16 min−1, k3 = 0.075 min−1, k4 = 0.033 min−1, BPND = 2.25, and VT = 6.5.
Flow changes were simulated by proportionally changing R1, defined as the K1 ratio between target and reference regions, while keeping the K1/k2 ratio constant. In the simulation, a change in flow is simulated by altering K1 (and k2), thereby assuming no change in first-pass extraction. This simulation reflects a change in flow between reference and cerebral AD regions (i.e., a heterogeneous flow change) at follow-up. R1 was varied from 0.6 to 1.4. In addition, a second simulation was performed by keeping K1′ in the reference region constant and only changing K1 in cerebral AD regions. Both K1 in cerebral AD and K1′ in cerebellar reference regions (and proportionally k2, keeping K1/k2 constant) were varied from K1 = K1′ = 0.19 to 0.48 mL·cm−3·min−1.
For all simulations, SUVrs were calculated for several uptake times and with several simulated flow variations. For comparison, BPND + 1 (=DVR) was obtained using RPM2 and reference Logan applied to the entire simulated 90-min reference and cortical time–activity curves. For all measures (i.e., SUVr, RPM2, and reference Logan-based BPND + 1), percentage bias compared with true or simulated DVR values were determined. In addition, the percentage change in all parameters as a result of both global and heterogeneous K1 changes was calculated by comparing these values with true or simulated BPND + 1 (that was kept constant at follow-up).
Simulations were performed both with and without adding noise to the time–activity curve, using a noise model according to Yaqub et al. (23). In the case of noisy simulations, 100 time–activity curves per simulation were generated. The average results from SUVr, RPM2, and reference Logan then were evaluated to study the effect of changes in R1 and K1. Noisy simulations showed results near-identical to those obtained without noise. Therefore, results from the simulations without noise will be reported, such that only flow or K1 effects are illustrated.
Statistics
Demographic and clinical differences between groups were assessed using ANOVA with post hoc least significant difference tests and age as covariate. At both baseline and follow-up, mean parameter values of 11C-PiB binding for the different methods were compared using ANOVA with post hoc least significant difference tests and age as covariate. Finally, group differences in R1 values at baseline and follow-up were examined using ANOVA with adjustment for age. R1 changes over time within groups were assessed using paired-samples t tests. Data are presented as mean ± SD, unless otherwise stated.
RESULTS
Demographic and clinical characteristics of the 3 diagnostic groups are presented in Table 1. Groups did not differ with respect to age or mean interval between baseline and follow-up.
Demographic and Clinical Characteristics According to Diagnostic Group
Baseline and Follow-up Binding Measures
Significant differences in 11C-PiB binding between groups, both at baseline and at follow-up, were found with all methods (Table 2). For all subjects together, SUVr60–90 was on average 14% and 13% higher at baseline and follow-up, respectively, than corresponding (BPND + 1) values obtained with RPM2. For SUVr40–60, these overestimations were 9% and 10% at baseline and follow-up, respectively. Reference Logan values were on average 6% lower than RPM2 values both at baseline and at follow-up.
Binding Values and Percentage Change of 11C-PIB According to Analytic Method
Longitudinal Changes in Binding Measures
Percentage changes between baseline and follow-up differed between methods, especially in AD patients (Table 2). Although BPND + 1 (RPM2), DVR (reference Logan), and SUVr40–60 were relatively stable (0% ± 6%, −1% ± 5%, and 0% ± 6%, respectively), SUVr60–90 decreased with 4% ± 8% in AD patients. Both SUVr measures showed larger variability than RMP2 values, especially at follow-up (Fig. 1). Differences were less pronounced for MCI patients and controls (RPM2: 6 ± 7, 2% ± 3%; reference Logan: 5 ± 6, 2% ± 3%; SUVr60–90: 8 ± 9, 3% ± 4%; and SURr40–60: 5 ± 7, 6% ± 6% for MCI and controls, respectively). Parametric images of an AD patient at both baseline and follow-up are presented in Figure 2 for all analytic models.
Global cortical binding of 11C-PiB in AD patients using RPM2 (A), reference Logan (B), SUVr60–90 (C), and SUVr40–60 (D).
Parametric images of AD patient with 27 mo of follow-up, showing increased RPM2 BPND + 1 (4%), reference Logan DVR (2%), and SUVr40–60 (2%), whereas SUVr60–90 decreased with 10%. R1 was 2% lower at follow-up.
Relative Tracer Delivery
At baseline, mean R1 values of 11C-PiB were 0.87 ± 0.05 in patients with AD, 0.89 ± 0.03 in MCI patients, and 0.87 ± 0.04 in controls. At follow-up, R1 values were 0.83 ± 0.06 in AD patients, 0.88 ± 0.03 in MCI patients, and 0.88 ± 0.04 in controls (Fig. 3). ANOVA adjusted for age revealed no differences in R1 between groups at baseline (F(2,28) = 0.65, P = 0.59). At follow-up, however, AD patients showed lower R1 values than MCI patients and controls at trend level (F(2,28) = 2.82, P = 0.06) (Fig. 3). Paired-samples t tests revealed a significant decrease in R1 over time in AD patients (t(6) = 2.85, P < 0.05), whereas no changes in R1 were found for MCI patients (t(10) = 1.73, P = 0.12) or controls (t(10) = 1.03, P = 0.33).
R1 values at baseline and follow-up for AD patients (▲ and △), MCI patients ( and ○), and controls (▪ and □). Significant decrease in R1 was found in AD group only. *P < 0.05.
Simulations
Simulations revealed that SUVr was dependent on uptake period (Figs. 4A and 5A). In general, SUVr overestimated BPND + 1 (=DVR) for all simulated K1 variations from 60 min after injection onward. Moreover, changes in both K1 (global flow changes) and R1 (heterogeneous flow changes), as shown in Figures 4B and 5B, induced both positive and negative bias in SUVr changes when compared with reference (i.e., baseline) conditions (K1 = 0.32 mL·cm−3·min−1 and R1 = 1, respectively). Bias and flow dependence of SUVr were larger than those of BPND + 1 (RPM2) and DVR (reference Logan), whereas the latter 2 methods showed comparable results.
(A) Percentage bias in SUVr (relative to BPND + 1) as function of time for various K1 values with R1 = 1 (i.e., K1 = K1′). For comparison, BPND + 1 obtained with RPM2 and reference Logan are indicated at 90 min after injection. (B) Percentage bias in change in SUVr (relative to change in BPND + 1) as function of time for various follow-up K1 values, baseline K1 = 0.32 mL·cm−3·min−1 and with R1 = 1 (i.e., K1 = K1′) both at baseline and at follow-up. For comparison, BPND + 1 obtained with RPM2 and reference Logan are indicated at 90 min after injection. RPM2 and reference Logan results for all simulated K1 values are plotted at 90 min after injection. x-axis represents mid-time of 10-min period for calculating SUVr measures.
Regional flow (i.e., K1 in simulations) variation effects (R1). (A) Percentage bias in SUVr, compared with true (simulated) DVR = BPND + 1 as function of uptake time. For comparison, results obtained using RPM2 and reference Logan are indicated at 90 min after injection. (B) Percentage change in SUVr as function of regional flow variation indicating artificial variation in SUVr as result of variations in R1. For comparison, results obtained using RPM2 and reference Logan are indicated at 90 min after injection. Small biases are now also observed for RPM2 and reference Logan but with smaller range than those seen with SUVr. RPM2 and reference Logan results for all simulated R1 values are plotted at 90 min after injection. For both methods, direction of bias (positive or negative) showed same trend as function of R1 as those seen with SUVr but with smaller amplitudes. x-axis represents mid-time of 10-min period for calculating SUVr measures.
DISCUSSION
This study directly compared changes in 11C-PiB binding parameters using 4 different analytic methods. It revealed marked differences between methods. Although both kinetic methods (RPM2 and reference Logan) showed relatively stable estimates of 11C-PiB binding in AD patients over time, SUVr60–90 demonstrated a decrease in 11C-PiB uptake. Although this was not observed with SUVr40–60, compared with the quantitative methods, both SUVr measures showed larger variability between subjects. This variability could be related to changes over time in relative tracer delivery (R1) to the region of interest, because this decreased over time in AD patients.
Both SUVr measures overestimated RPM2 values with 9%–14%. Carson et al. (14), using 18F-cyclofoxy, showed that tissue ratios such as SUVr—due to its sensitivity to differences in clearance rate—can overestimate specific binding substantially. In addition, it was shown that this bias was different for high- and low-binding areas. Using 11C-PiB, Lopresti et al. (19) were the first to describe that SUVr40–60 and SUVr40–90 showed a large positive bias, compared with quantitative methods, but that the percentage bias was fairly similar between low- and high-binding areas. Overestimation of SUVr, compared with quantitative methods, has also been observed using the novel amyloid tracers 18F-flutametamol (24) and 18F-florbetapir (25). Reference Logan estimates were about 6% lower than RPM2 estimates, likely due to the sensitivity of reference Logan graphical analysis to statistical noise (26). In the present study, there was a difference of 17%–19% between SUVr and reference Logan values. In a recent study in MCI patients (27), SUVrs were up to 31% higher than reference Logan DVRs. These discrepancies between methods have important implications for determination of quantitative thresholds for PiB positivity, because these thresholds depend critically on the method that is used to analyze the data.
Both SUVr measures showed larger variability than RPM2 values, especially at follow-up, which can be explained by 2 processes. First, at later time points the SUVr becomes constant over time. Yet, SUVr will overestimate the true DVR as the tissue curves follow the clearance of plasma input curves. In other words, there is no true equilibrium between plasma concentration and tissue concentration (because the concentration in plasma is lower than in tissue, there will be a net transport of tracer from tissue to plasma). The violation of equilibrium between reference and target regions will be different because of differences in specific binding between these regions and subjects. Because of these differences in nonequilibrium between subjects, there will more variability of SUVr between subjects. With the performance of dynamic scans and kinetic analysis, the change and variability of the input function (that affects both reference and target region) is considered, and consequently BPND estimated with RPM2 can be expected to be more reproducible. A second reason may be the variations in flow. To see to what degree the various parametric methods (SUVr, reference Logan, and RPM2) are affected by flow, we performed several simulations. In general, the simulation data were in good agreement with the clinical findings. First, SUVr generally overestimated simulated BPND + 1 (=DVR) values, whereas RPM2 and reference Logan–based values showed minimal bias. Second, these simulations showed that both global and regional (i.e., heterogeneous) K1 changes over time could result in artificial changes in SUVr over time, a phenomenon that was less prominent for RPM2 or reference Logan. Yet, some small biases and dependence on changes in K1 were also observed for RPM2 and reference Logan, which can be explained by the contribution of signal from the blood volume fraction. Blood volume fractions were included in the simulations to generate realistic time–activity curves but are, by definition, not considered by both RPM2 and reference Logan. Finally, results obtained using SUVr strongly depended on the specific uptake interval, with slightly less bias and flow dependence for earlier time intervals. This is consistent with clinical data where SUVr40–60 seemed to show less variability over time than SUVr60–90. It is already known that, in AD patients, flow changes over time occur because of disease progression (28), and, consequently, it may be expected that this will result in more variability in SUVr over time. Changes in R1 seen in AD patients indicate that these flow changes are indeed present. Therefore, it is highly likely that changes in SUVr observed in the present series of AD patients do not reflect changes in specific 11C-PiB binding but are rather due to changes in perfusion during the course of the disease. Apart from heterogeneous flow changes (reflected by changes in R1), also relatively large global flow changes are likely to occur both in healthy subjects and in patients. A recent study by Bremmer et al. (29) showed that day-to-day variations in global cerebral blood flow were about 30% under normal conditions. As shown by the simulations, these global flow variations could add to the clinically observed variability in SUVr. Similar effects of flow on SUVrs may be expected when using radiotracers other than 11C-PiB. The flow dependence is caused by the lack of equilibrium of tracer distributions between blood and tissue and the tissue compartments. Therefore this flow dependency occurs for any tracer, although the degree of this effect differs between tracers depending on their kinetic behavior. As such, these findings imply that any study (but particularly when using 11C-PiB or tracers with similar kinetic behavior) where variations in blood flow can be expected should not be analyzed using SUVr. Consequently, for accurate quantification of longitudinal amyloid imaging studies, dynamic scanning protocols and fully quantitative data analysis methods are essential. This is especially true for longitudinal studies with disease-modifying agents aiming to lower amyloid load in the brain. Studies that use suboptimal methods such as SUVr carry the inherent risk that ineffective drugs are not identified appropriately or, more importantly, that potential effective drugs are dismissed, especially when effect sizes are small. Recently, effects of bapineuzumab on fibrillar amyloid load in AD patients, as measured using 11C-PiB and PET, were reported (8). This study found a significant reduction in mean 11C-PiB uptake across 6 targeted ROIs in patients in the treatment arm, compared with those in the placebo group. This is a landmark study, because it was the first, to our knowledge, to show a central effect of a therapeutic approach, aimed at lowering cerebral amyloid load in patients with AD. For the analysis of 11C-PiB PET scans, SUVr60–90 was used. Results of the present study, however, indicate that SUVr60–90 is susceptible to flow changes, which may be different in bapineuzumab and placebo groups. As such, it is impossible to differentiate between decreases in amyloid load due to treatment and decreases due to flow artifacts. This differentiation is important because potentially, with more stringent methodology, the ineffectiveness of bapineuzumab (30) in the treatment of AD could have been identified in a much earlier phase of development.
One could argue that repeated 90-min dynamic scans cause a selection bias, because only relatively few patients can undergo such a demanding procedure. This is indeed debatable for patients with moderate to severe AD. It is, however, most likely that amyloid-lowering drugs are most effective in the early stages of the disease, justifying inclusion of patients with mild AD or, preferably, prodromal AD or individuals with autosomal dominant AD in a preclinical stage. Dynamic scanning protocols and fully quantitative data analysis methods are necessary in these patients, as a solid baseline measurement to monitor their disease course and treatment response is needed because it is these patients who will progress to advanced stages with associated blood flow changes. These subjects are perfectly capable of undergoing dynamic scans when carefully prepared and monitored during scanning.
In the case of longitudinal 11C-PiB studies, a steady-state approach could be an alternative for dynamic scanning (31). Steady-state levels of 11C-PiB can be achieved using a bolus-with-continuous-infusion protocol, which can be performed outside the PET camera. When a reliable steady state has been achieved, usually a short scanning period is sufficient. This method would combine advantages of a good quantitative measure that is independent of (relative) flow changes with short scan duration. However, this method still needs to be tested and validated for 11C-PiB.
CONCLUSION
SUVr should not be used for longitudinal 11C-PIB studies, especially when only small changes in specific binding can be expected.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. This work was financially supported by the Internationale Stichting Alzheimer Onderzoek (ISAO, grant 05512) and the American Health Assistance Foundation (AHAF, grant A2005-026). No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online Aug. 12, 2013.
- © 2013 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication September 10, 2012.
- Accepted for publication April 19, 2013.