Abstract
123I-CLINDE is a radiotracer developed for SPECT and targets the 18-kDa translocator protein (TSPO). TSPO is upregulated in glial cells and used as a measure of neuroinflammation in a variety of central nervous system diseases. The aim of this study was to examine the test–retest variability of 123I-CLINDE binding in healthy subjects. Methods: SPECT scans were acquired over 90 min in 16 healthy controls (9 women, 8 mixed-affinity binders [MABs] and 8 high-affinity binders [HABs] twice with an interval of 35 ± 15 d). Arterial input functions were based on individual blood measurements in 8 subjects and a population-based approach in combination with individual whole-blood time–activity curves in the other 8 subjects. Seven brain volumes of interest were extracted and quantified by SUVs and by 2-tissue-compartment modeling for calculation of distribution volumes (VT). Test–retest variability was measured by percentage difference (PD), the absolute PD, intraclass correlation coefficient (ICC), and coefficient of variation. Results: The absolute PD of brain SUV and the VT had similar values. The ICC values were higher for VTs than for brain SUVs, which were both moderate to high; however, lower ICC values were observed when calculated separately for HABs and MABs. Test–retest reproducibility was higher in subjects with immediate centrifugation of blood samples. The population-based method efficiently recovered data with delayed centrifugation. The VT of a 49-y-old male HAB was 7.5 ± 1.4 mL/cm3 compared with 4.6 ± 1.4 mL/cm3 of a sex- and age-matched MAB. The SUVs of a 49-y-old male HAB and MAB were 1.03 ± 0.14 and 0.88 ± 0.15 g/mL, respectively. Conclusion: The test–retest reproducibility of 123I-CLINDE is comparable or better than that reported for commonly used PET TSPO tracers. Because of the binding of 123I-CLINDE to blood cells and peripheral tissues, SUV is not a sufficient surrogate of VT from 2-tissue-compartment modeling. The population-adjusted method has the potential to reduce the complexity of blood analyses of TSPO tracers.
In brain parenchyma, the translocator protein (TSPO) is predominantly expressed by glial cells (astrocytes and microglia). Subcellularly, TSPO is mainly located in the outer mitochondrial membrane and functioning as a transporter of cholesterol and porphyrins and importer of proteins (1). In the blood, TSPO is expressed mainly in monocytes and polymorphnuclear neutrofils and to a lesser extent in other white blood cells, platelets, and erythrocytes (2).
After brain injury, TSPO levels are increased within activated microglia and astrocytes, inspiring the investigation of TSPO as a biomarker of brain injury and repair in brain diseases such as stroke, glioma, and neurodegenerative disorders (3–6). TSPO molecular imaging is currently exclusively performed in research protocols but may have future clinical roles in diagnosis and therapeutic drug monitoring in neuroinflammatory disorders.
Test–retest variability is a way to assess the usefulness of radioligands; a radioligand with high reproducibility is better suited for longitudinal studies and requires smaller sample sizes in group comparisons. So far, TSPO imaging test–retest studies have demonstrated moderate to good reproducibility for 11C-PK11195 (7) and better reproducibility for the second-generation tracers 11C-DPA-713 (8) and 11C-PBR28 (9,10).
123I-CLINDE is a SPECT TSPO ligand (11). 123I-CLINDE SPECT imaging has previously been applied in patients with stroke, glioblastoma multiforme, and anti N-methyl-d-aspartate receptor encephalitis (12–14). 123I-CLINDE is susceptible to the rs6971 polymorphism, which affects the affinity of the radioligand binding to the TSPO and separates individuals into high-affinity binders (HABs), mixed-affinity binders (MABs), and low-affinity binders (15).
The aim of the current study was to examine the test–retest reproducibility measured as percentage difference (PD) and reliability measured as intraclass correlation coefficient (ICC) and coefficient of variation (COV) of 123I-CLINDE binding in healthy subjects based on the brain SUV and 2-tissue-compartment modeling (2TCM) with arterial blood as the input function.
MATERIALS AND METHODS
Subjects and Genotyping
This study was conducted in accordance with the Declaration of Helsinki at the Copenhagen University Hospital, Rigshospitalet, Denmark. The ethical committee of the Copenhagen Capital Region (H-2-2010-086) approved the study protocol. All subjects provided written informed consent.
Because of low TSPO signal in glioblastoma tissue from a patient who was a low-affinity binder (12), we did not expect to obtain reliable TSPO signals in healthy low-affinity binders. Consequently, after screening for the rs6971 polymorphism (12), 16 healthy subjects (MAB-to-HAB, 1:1; 9 women) were included and SPECT-scanned twice with an interval of 35 ± 15 d. All healthy subjects had normal physical and neurologic examinations and blood tests. Table 1 shows a summary of age, sex, TSPO genotype, and scan information.
Arterial Blood Sampling and Centrifugation
From a radial arterial cannula, blood samples for measuring radioactivity in plasma and whole blood were drawn manually at 0.25, 0.5, 0.75, 1, 1.5, 2, 2.5, 3, 4, 6, 8, 10, 12, 15, 20, 25, 30, 35, 45, 65, and 85 min after injection of 123I-CLINDE. Additionally, samples were drawn at 0.5, 4, 10, 20, 30, 45, 65, and 85 min for metabolite analysis by radio–high-performance liquid chromatography (HPLC) (16). Fractions of HPLC eluent were collected using a fraction collector device (Foxy Jr FC144; Teledyne Isco) and counted off-line in a well-counter (2480 Wizard2 γ-Counter; Perkin Elmer) for accurate measurements of parent tracer. For the first 8 subjects, the blood was stored on ice for an estimated 30–90 min until centrifugation (Eppendorf centrifuge 5800R with 3,500 rpm at 4°C in 7 min). During the course of the study, we discovered that 123I-CLINDE continued to bind to TSPO in formed blood cells, presumably the monocytes, until separating plasma from the sample by centrifugation. For the remaining 8 subjects, all blood samples were stored on ice and centrifuged within a maximum of 5 min. The total radioactivity in plasma and whole blood was measured in a well-counter (Cobra 5003; Packard Instruments), and data were decay-corrected to the time of injection. Further, the radioligand purity was measured for each batch of 123I-CLINDE, and trapping efficiency of the HPLC column (plasma control) was determined for each scan by spiking water and blank plasma with 123I-CLINDE. All equipment was cross-calibrated to the SPECT scanner. Blood screening including white blood cell differential count was performed at each scan.
Image Acquisition
To block thyroidal uptake of free radioiodine, 200 mg of potassium perchlorate were administered intravenously 20 min before injection of 123I-CLINDE. A headband was applied to minimize head movements during scanning. Dynamic SPECT scanning with a triple-head IRIX camera (Philips Medical) was started simultaneously with a bolus injection of 123I-CLINDE (MAP Medical Technologies) and lasted for 90 min. The protocol consisted of 10 × 2 min frames, followed by 7 × 10 min frames.
T1-weighted MRI was acquired once for each subject using a 3-T Prisma MR scanner (Siemens) 3 ± 49 d (mean ± SD; median, 0 d) after the corresponding SPECT scan. All brain scans were normal as interpreted by a neuroradiologist.
Image Analysis
Image Preprocessing
For each subject, the weighted mean SPECT image was coregistered to the T1-weighted MR image by interactive image overlay (17). Automatic delineation of volumes of interest (VOIs) was performed using probability maps and a data processing pipeline (18). Because of the limited spatial resolution of SPECT images, no segmentation of gray matter, white matter, and cerebrospinal fluid was attempted. Seven VOIs were investigated in this study: midbrain, thalamus, pons, striatum, hippocampus, cerebellum (without vermis), and neocortex. These VOIs were defined as volume-weighted means of either the left- and the right-hemisphere or of smaller VOIs: striatum as the mean of caudate nucleus and putamen and neocortex as the mean of frontal cortex, superior temporal gyrus, medial inferior temporal gyrus, parietal cortex, sensory motor cortex, and occipital cortex.
Brain and Blood Uptake
The radioactive concentrations in blood and brain were studied independently to investigate the test and retest variability. Therefore, PD of the SUVmean was calculated, where SUV (g/mL) was defined as time–activity curve (kBq/mL) in the blood or brain VOIs normalized by the injected dose per body weight (in MBq/kg) (19), and SUVmean was defined as the area under the SUV curve normalized by time length.
Kinetic Modeling
A 2TCM with arterial plasma function as input (K1 to k4 with a fixed blood volume in brain tissue of 5%) was used to quantify 123I-CLINDE binding as distribution volume (VT), as this has previously been demonstrated to be the optimal method (12).
The arterial input function was generated in 2 ways, depending on the blood sample centrifugation.
Method 1: For Studies with Immediate Centrifugation of Blood Samples
To account for tracer purity and trapping efficiency, the measured parent fraction (PF) from the fraction collector was first fitted by a multiexponential function and adjusted as below:Eq. 1The difference between 123I-CLINDE purity and plasma control indicates the part failed to be trapped by the extraction column of the radio-HPLC.
Then the metabolite-corrected plasma input function was calculated as the product of the plasma time–activity curve Cp(t) and the adjusted parent fraction PF_adj(t):Eq. 2
Method 2: For Studies with Delayed Centrifugation of Blood Samples
The plasma–to–whole-blood ratio (PoB) and adjusted parent fraction as a function of time were generated for all scans that were accompanied by blood samples undergoing immediate centrifugation. The mean of these functions, respectively, were generated as population-based PoB (PB_PoB) and population-based parent fraction (PB_PF):Eq. 3where N is the number of scans. To generate the metabolite-corrected plasma input function (Cmcp), these 2 functions were then applied on the measured whole-blood time–activity curve (Cw) of the studies with delayed centrifugation as below:Eq. 4Image and data preprocessing were performed using MATLAB 8.1 (R2013a; Mathworks Inc.), and kinetic modeling was done using PMOD (version 3.0; PMOD Technologies Inc.).
Statistical Analysis
To evaluate the test–retest reproducibility and variability of 123I-CLINDE, the following 4 matrices were used: PD, absolute PD, ICC, and COV.
PD
Test–retest reproducibility was calculated as the difference between test and retest outcome measures divided by their average, outcome measures being SUVmean of blood or brain uptake, and VT. The absolute PD was also calculated for comparison to other studies.
ICC
ICC evaluates the reliability by measuring within-subject variability relative to between-subject variability as below:Eq. 5whereEq. 6MSSbtw is the mean sum of squares between subjects and MSSwithin is the mean sum of squares within subjects, is the mean of test and retest measures of subject i, and is the mean of all observations (test and retest of all subjects). is one of the measures of subject i. N is the number of subjects, and k is the number of repeated measures (k = 2). If MSSbtw = 0 there is no reliability, and ICC = −1; if MSSwithin = 0 there is maximum reliability, and ICC = 1.
COV
The COV of outcome measures was calculated as the SD divided by the mean across all observations for each VOI.
To test effects of centrifugation, genotype, or the test–retest, a linear model fitting using generalized least squares (GLS) was performed in R (version 3.2.2). In general, a GLS model performs similar to a standard linear model but can furthermore account for a complex covariance structure between observations. In our case, the GLS model allows for different variances between residuals in repeated measures and for correlation between them, unless otherwise stated.
Details about each model applied will be given in the “Results” section. For GLS models, mean and SE are given. In other cases, mean and SD are given. The paired-sample t test, Wilcoxon rank-sum test, and Pearson linear correlation were performed using MATLAB 8.1 (R2013a; Mathworks Inc.).
RESULTS
The injected radioactivity in test and retest scans was not statistically different, and neither was 123I-CLINDE purity or plasma control (Table 1).
Delayed Centrifugation and Discrepancy Between Blood and Brain PD
In the first 8 subjects with delayed centrifugation of blood samples, we found a large test–retest PD of SUVmean in the plasma parent compound (1.0% ± 31.3%), but much lower PD in the brain (0.7% ± 8.0% - 2.2% ± 18.8% in 7 VOIs).
Effect of Genotype
In subjects with immediately centrifuged samples, the parent fraction SUVmean was higher in MABs (0.37 ± 0.07, mean ± SE) than in HABs (0.22 ± 0.05, mean ± SE; P = 0.0006, GLS model, age and sex adjusted), with a MAB-to-HAB ratio of 1.65. A genotype-related difference was also observed in PoB (Supplemental Fig. 1; supplemental materials are available at http://jnm.snmjournals.org). The mean PoB time–activity curve of HABs (n = 2 × 3) and that of MABs (n = 2 × 5) had different kinetics, especially during the first 20 min; however, the difference was not significant. The SUVmean 0–20 min of MABs correlated negatively with monocyte counts (P = 0.0341, R2 = 0.4486). A negative trend was also observed for HABs (P = 0.1938).
Arterial Input Function
Immediately Centrifuged Blood Samples
After parent fraction adjustment using Equation 1, the parent fraction between MAB and HAB was still significantly different (P = 0.0002, GLS model, age and sex adjusted). Furthermore, the parent fraction curves had different starting points: 94% intact 123I-CLINDE left for MABs and 77% for HABs; at 20 min after injection, 53% for MABs and 36% for HABs; and at 90 min, 38% left for MABs and 27% for HABs.
The arterial input function of a MAB and a HAB (Figs. 1C and 1D) showed faster washout in the HAB. Because of the observed differences between genotypes, population-based PoB functions and population-based PF functions were generated separately for these 2 genetic groups. Before applying the population-based functions to adjust data with delayed centrifugation, method validation was done on the data with immediate centrifugation using a leave-one-out procedure, where population-based functions were generated excluding the scan that the functions were applied on.
These PB_PoB functions were applied to the whole-blood time–activity curve Cw first to estimate the plasma time–activity curve Cp. The differences between estimated Cp and the measured data were 0.2% ± 4.3% for MABs (n = 2 × 5) and 0.6% ± 10.8% for HABs (n = 2 × 3); none of the differences was significant. Thereafter, PB_PF was applied to Cp to estimate the metabolite-corrected plasma input function Cmcp. Differences between population-adjusted estimates and measured Cmcp were 0.6% ± 9.1% for MABs and 0.5% ± 8.0% for HABs; none of the differences was significant.
Delayed Centrifuged Blood Samples
Immediate centrifugation increased the SUVmean of PoB from 0.95 to 1.07 (P = 0.0003, GLS model, genotype adjusted), which was 12.4% more than the SUVmean of PoB from subjects whose samples underwent delayed centrifugation (descriptive data of PoB are given in Fig. 2). Thus, the arterial input functions were underestimated in scans with delayed centrifugation of blood samples (Figs. 1A and 1B). Therefore, population-adjusted input functions were calculated in these cases.
After adjustment using PB_PoB and PB_PF on Cw, test–retest PD of Cp decreased from −22.7% ± 27.6% to −9.4% ± 15.2% (absolute PD of Cp decreased from 27.7% to 14.4%), and PD of Cmcp changed from −2.5% ± 34.6% to −6.4% ± 15.5% (absolute PD of Cmcp decreased from 28.0% to 14.1%) (Supplemental Fig. 2).
Test–Retest Variability
Brain Uptake
The brain SUVmean of 7 VOIs was studied. The SUVmean of a male HAB and a male MAB, at age 49 y, were 1.03 ± 0.14 and 0.88 ± 0.15 g/mL (mean ± SE), respectively, and the difference was borderline significant (P = 0.0535, GLS with different residual variances). The higher binding in HABs than in MABs was also shown in Figure 1, illustrated as SUVs of neocortex with slower washout in HABs. The SUVmean of HABs was 16.8% higher than that of MABs (descriptive data are shown in Fig. 3A). Ten of 16 subjects had lower brain SUVmean in test than retest scans.
Because the beginning of a brain time–activity curve can be influenced by blood flow, SUVmean 60–90 min was also studied. The SUVmean 60–90 min of a male HAB and a male MAB aged 49 y were 0.55 ± 0.09 and 0.43 ± 0.09 g/mL (mean ± SE), respectively (P = 0.0105, GLS with different residual variances), and the SUVmean 60–90 min of HABs was 27.6% higher than MABs (descriptive data are provided in Fig. 3B).
The PDs of brain SUVmean across all subjects in 7 VOIs are shown in (Table 2). Test and retest SUVmean of the same subject correlated to a high degree, with correlation coefficients of 0.86 ± 0.07 (P < 0.05, n = 16). The PD of MABs was on average −3.6% ± 8.1% whereas that of HABs was −7.0% ± 7.8%; however, the difference was not significant (GLS model, age and sex adjusted), and the absolute PD was not sensitive to genotype. ICC values for MABs and HABs and the total samples are given in Table 2. The COV was 20.0% ± 2.3% in MABs and 9.2% ± 2.0% in HABs. Furthermore, ICC and COV differed between MABs and HABs (P = 2.29e-05, P = 1.38e-05, respectively), with a MAB-to-HAB ratio of 1.8 for ICC and 2.2 for COV.
VTs
The VTs of 7 VOIs derived from 2TCM of scans with immediately and delayed centrifuged blood samples, with and without population-based adjustment of the input functions, are given in Supplemental Tables 1–3. Figure 1 illustrates the 2TCM fitting using the measured or population-based adjusted input functions.
The VTs of scans with delayed centrifuged blood samples (16 scans, 7 VOIs) before and after population-based adjustment was significantly different (P = 0.0176, GLS with different residual variances), with a VT value of 5.3 ± 1.1 mL/cm3 before the adjustment and 4.7 ± 1.1 mL/cm3 (mean ± SE) after the adjustment. Delayed centrifugation of samples led to an underestimation of input functions and therefore caused an overestimation of the VTs by 12.4% on average.
The VTs (32 scans, 7 VOIs) of a male HAB and a male MAB aged 49 y (mean age) were 7.5 ± 1.4 and 4.6 ± 1.4 mL/cm3 (mean ± SE), respectively (P = 0.0001, GLS with different residual variances) using population-based adjusted input functions in scans with delayed blood centrifugation and measured input functions in scans with immediate blood centrifugation. VTs of HABs were on average 65% higher than MABs, and there was no test–retest difference (Fig. 3C). Test–retest variability measures were calculated for the immediately centrifuged scans and the population-adjusted delayed centrifuged scans both separately and combined (Table 2). For subject 16, the measured purity and plasma standard were not available, thus no adjustment was done.
For the immediately centrifuged scans, test VT was 0.28 mL/cm3 lower than retest VT across all 7 VOIs (P = 0.0122, GLS model, age and sex adjusted). This GLS model assumes a homogeneous test–retest difference across VOIs, because there were no interactions between the 7 VOIs and test–retest scans.
The PD between test–retest of all subjects (n = 16) across all 7 VOIs evaluated as repeated measures was on average −11.7% ± 8.6% for MABs and −12.1% ± 8.3% for HABs, and the difference was not significant (GLS, age and sex adjusted). The COV was 16.7% ± 3.4% for MABs and 26.5% ± 7.7% for HABs. Opposite to the brain SUVmean, ICC and COV estimated by the VTs were borderline significantly higher for HABs than MABs (P = 0.0533, P = 0.0492, respectively), with a HAB-to-MAB ratio of 1.6 for ICC and 1.7 for COV.
DISCUSSION
In this study, we investigated the test–retest variability of the second-generation TSPO SPECT ligand 123I-CLINDE in a group of 16 healthy volunteers.
Population-Based Approach
We used a population-based approach to adjust data from subjects with delayed blood centrifugation. The feasibility of the approach was validated on the immediately centrifuged data in a leave-one-out manner to avoid bias. After adjustment, the test–retest difference decreased in blood components. If time-dependent binding to blood cells is also observed in other TSPO ligands, this approach could be applied to data if centrifugation is delayed. Furthermore, this approach has the potential to reduce the complexity and time consumption associated with blood handling and determination of metabolites, because only whole-blood counts are needed.
123I-CLINDE SPECT Variability
ICC values of VTs were higher for 123I-CLINDE SPECT than previously reported in healthy controls for the first-generation TSPO tracer 11C-PK11195 PET (7) and comparable to second-generation TSPO tracers such as 11C-PBR28 (9) and 11C-DPA-713 (8). ICC is a statistical analysis to compare the between-subject variation with the within-subject variation. If HAB and MAB subjects are pooled and not measured separately for each genetic group, ICC values will be higher and reflect not only the variability of the method but also a known biologic difference in affinity. In contrast to previous studies in second-generation TSPO tracers, we also calculated ICC and COV separately for each genetic group, revealing higher ICC and COV values for HABs than MABs. A plausible explanation for a generally modest reliability when subdividing into genotypes is that the TSPO signal, because of its nonprominent presence under healthy conditions, has a relatively low signal-to-background ratio compared with a disease situation with an amplified TSPO signal driven by disease-related activated glial cells. For MABs, the signal-to-background ratio is lower than for HABs, consequently leading to a comparably lower reliability. PD of brain SUVmean and VT and the absolute PD for 123I-CLINDE SPECT were slightly smaller than previously reported for 11C-PBR28 (9).
SUV as Semiquantitative Measure
In a clinical setting, SUV is often used as a semiquantitative surrogate measure to VT. This is due to several complications related to the complexity of the blood sampling and analysis. In this paper, we studied the genotype difference and test–retest variability using both SUV and VT to investigate the feasibility of using SUV as an outcome measure for a TSPO ligand. From the brain SUVmean, it is possible to differentiate MAB and HAB (Fig. 4), with HABs having an approximately 17% higher binding than MABs when quantified using SUVmean of the whole time–activity curve and approximately 28% higher when quantified using SUVmean 60–90 min. However, HABs had on average a 65% higher binding than MABs using VT estimated by 2TCM modeling. This can be explained by higher affinity of radioligand to TSPO in HABs and the binding of radiotracer to the abundant TSPO in both blood and peripheral tissues. Thus, less free radioligand would be available in the plasma of HABs compared with MABs, which would remain undetected using brain SUV as an outcome. In a clinical setting, patients can have a longitudinally varying disease-related peripheral immune response and thereby varying peripheral TSPO binding sites. This may cause changing concentrations of parent TSPO tracer in the plasma, which as a consequence influences the brain uptake. In Yoder et al. (20), SUV was compared with VT in quantifying TSPO changes in baboons’ brain caused by lipopolysaccharide using 11C-PBR28, and SUV failed to detect neuroinflammation as indicated by VT changes. It might be ascribed to the violation of assumptions: the arterial input time–activity curve has a consistent pattern from study to study, and the area under the curve of time–activity curve is proportional to the injected dose per body weight. Our study clearly indicated the violation of these assumptions in the case of TSPO imaging. Therefore, using only brain SUV to quantify TSPO binding has its limitations.
Effect of Blood Centrifugation
Switching to immediate centrifugation improved the test–retest reliability in blood significantly. Because TSPO is also expressed in blood, for example, in monocytes (2), we hypothesized that 123I-CLINDE continues to bind to TSPO in blood cells and that intact 123I-CLINDE in a blood sample will gradually bind to blood cells causing an underestimation of plasma radioactivity. Our hypothesis is supported by the increase in PoB ratio by immediate centrifugation. The negative correlation between PoB in MABs and monocyte counts also indicate that CLINDE binding is related to TSPO expression in blood. Additionally, the lower starting point of PoB time–activity curve of HABs than time–activity curve of MABs (Supplemental Fig. 1) indicates that binding of CLINDE to blood cells is dependent on the affinity/TSPO polymorphism (21). Previously, in Rizzo et al. (22) an extra compartment was proposed to describe TSPO ligand binding to the vascular walls. On the basis of our observations, the endothelial compartment may additionally include binding TSPO in blood cells. However, fast and systematic centrifugation of blood samples will still be required to produce consistent results.
CONCLUSION
The variability of 123I-CLINDE binding was investigated in a test–retest setting with 16 healthy subjects. The test–retest variability estimated by ICC and PD was favorable compared with the first-generation TSPO ligand 11C-PK11195 and comparable to or slightly better than second-generation TSPO ligands 11C-DPA-713 and 11C-PBR28. SUV as a surrogate of the output measure VT from gold standard 2TCM has its limitations. Immediate centrifugation of blood samples is essential and can prevent underestimation of the plasma input function. A population-based method could efficiently recover data with delayed centrifugation and has the potential to be applied to studies using other TSPO ligands.
DISCLOSURE
This work was financially supported by the European Union's Seventh Framework Program (FP7/2007-2013) under grant agreement no. HEALTH-F2-2011-278850 (INMiND) and the Danish Council for Independent Research (grant ID 0602-0228-8B). No other potential conflict of interest relevant to this article was reported.
Acknowledgments
We acknowledge Svitlana Olsen and Glenna Skouboe for technical assistance and Brice Ozenne for his assistance in statistical analysis.
Footnotes
↵* Contributed equally to this work.
Published online Nov. 10, 2016.
- © 2017 by the Society of Nuclear Medicine and Molecular Imaging.
REFERENCES
- Received for publication August 27, 2016.
- Accepted for publication October 21, 2016.