ORIGINAL RESEARCH article
Noninvasive Relative Quantification of [11C]ABP688 PET Imaging in Mice Versus an Input Function Measured Over an Arteriovenous Shunt
- 1Molecular Imaging Center Antwerp, University of Antwerp, Wilrijk, Belgium
- 2Department of Nuclear Medicine, Antwerp University Hospital, Edegem, Belgium
- 3Bio-Imaging Lab, University of Antwerp, Wilrijk, Belgium
- 4CHDI Foundation, Princeton, NJ, United States
Impairment of the metabotropic glutamate receptor 5 (mGluR5) has been implicated with various neurologic disorders. Although mGluR5 density can be quantified with the PET radiotracer [11C]ABP688, the methods for reproducible quantification of [11C]ABP688 PET imaging in mice have not been thoroughly investigated yet. Thus, this study aimed to assess and validate cerebellum as reference region for simplified reference tissue model (SRTM), investigate the feasibility of a noninvasive cardiac image-derived input function (IDIF) for relative quantification, to validate the use of a PET template instead of an MRI template for spatial normalization, and to determine the reproducibility and within-subject variability of [11C]ABP688 PET imaging in mice. Blocking with the mGluR5 antagonist MPEP resulted in a reduction of [11C]ABP688 binding of 41% in striatum (p < 0.0001), while no significant effect could be found in cerebellum (−4.8%, p > 0.99) indicating cerebellum as suitable reference region for mice. DVR-1 calculated using a noninvasive IDIF and an arteriovenous input function correlated significantly when considering the cerebellum as the reference region (striatum: DVR-1, r = 0.978, p < 0.0001). Additionally, strong correlations between binding potential calculated from SRTM (BPND) with DVR-1 based on IDIF (striatum: r = 0.980, p < 0.0001) and AV shunt (striatum: r = 0.987, p < 0.0001). BPND displayed higher discrimination power than VT values in determining differences between wild-types and heterozygous Q175 mice, an animal model of Huntington's disease. Furthermore, we showed high agreement between PET- and MRI-based spatial normalization approaches (striatum: r = 0.989, p < 0.0001). Finally, both spatial normalization approaches did not reveal any significant bias between test-retest scans, with a relative difference below 5%. This study indicates that noninvasive quantification of [11C]ABP688 PET imaging is reproducible and cerebellum can be used as reference region in mice.
Glutamate is the most prominent neurotransmitter in the brain. The metabotropic glutamate receptors (mGluRs) are G-protein coupled receptors which modulate synaptic transmission and neuronal excitability (1). Impairment of the mGluR Group I (mGluR1 and mGluR5) has been implicated with various neurologic disorders, including Huntington's disease (HD) (2). HD is an autosomal dominant neurodegenerative disorder (3) caused by an expanded CAG repeat in exon 1 of the gene encoding the protein huntingtin (HTT) (4). Subjects with HD exhibit as main neuropathological feature a progressive neuronal cell loss in the caudate-putamen (5), which receives input from different areas of the basal ganglia as well as glutamatergic inputs from thalamus and cortex (6). Thus, mGluR5 is hypothesized to play an important role in the pathogenesis of HD (7, 2) and it represents an interesting target to image in vivo by means of positron emission tomography (PET). Among the developed radiotracers to image mGluR5, there is [11C]ABP688 (3-(6-methyl-pyridin-2-ylethynyl)-cyclohex-2-enone-O-11C-methyl-oxime), which binds to the mGluR5 allosteric-binding site with high affinity and selectivity (8). Although [11C]ABP688 PET imaging has successfully been used to quantify mGluR5 at preclinical and clinical level (8–11), the methods for reproducible quantification of [11C]ABP688 PET imaging in mice have not been thoroughly investigated yet.
Volume of distribution (VT) quantification requires the knowledge of the arterial input function (IPF), which can be measured by serial arterial blood sampling. However, in small animals, the amount of blood is very limited and therefore we applied an arteriovenous (AV) shunt coupled with a γ coincidence counter (12). The AV shunt surgery in small animals is typically an end of life procedures and does not allow longitudinal imaging over months in the same animals. A noninvasive approach is the extraction of the whole blood input function from the PET images. This image-derived input function (IDIF) is typically obtained from the lumen of the left ventricle of the heart (13) end-of-diastole in ECG monitored animals. In order to assess [11C]ABP688 quantification for longitudinal studies, first we determined the reproducibility of the IDIF by comparison with blood sampling via AV shunt and γ coincidence counter.
Nevertheless, the need of an IPF could be overcome altogether in the presence of a reference region in the brain devoid of the targeted receptors by using reference tissue models (14). Previous studies with [11C]ABP688 have shown cerebellum or cerebellar gray matter to be an optimal reference region in rats (10), baboons (11), and in humans (15) although mGluR5 is known to be present in small quantities in the cerebellum (16, 17). The valid use of a reference region for [11C]ABP688 quantification has not been investigated in mice. For this reason, the second aim of this study was to evaluate and validate the reference tissue model with the cerebellum as reference region for [11C]ABP688 PET imaging in mice.
Additionally, as HD is characterized by brain atrophy, primarily in striatum and cortex (5), we investigated the potential benefit of individual magnetic resonance (MR) images compared to the use of a [11C]ABP688 PET template for spatial normalization for quantification of [11C]ABP688 PET binding in the mouse brain.
Previous preclinical test-retest studies in rats and baboons indicated that [11C]ABP688 provides reproducible outcome measures with an average percentage difference below 10% (11, 18). On the other hand, same day test-retest studies in rhesus monkeys and in humans reported larger intra-individual variability, with an average increased in [11C]ABP688 uptake during the retest scan (15, 19, 20). As test-retest stability of [11C]ABP688 has never been investigated in mice, the final aim of this study was to determine the reproducibility and within-subject variability.
Materials and Methods
Six months old (n = 29 per genotype) and 9 months old (n = 4 per genotype) male heterozygous (HET) Q175 mice (21, 22) and wild-type (WT) Q175 littermates were obtained from Jackson Laboratories (Bar Harbour, Maine, USA) were included in this study. Animals were single-housed in individually ventilated cages under a 12 h light/dark cycle in a temperature- and humidity-controlled environment with food and water ad libitum. The animals were acclimatized to the facility for at least 1 week before the start of procedures, which were performed according to the European Committee Guidelines (decree 2010/63/CEE) and the Animal Welfare Act (7 USC 2131). All experiments were in compliance with the ARRIVE guidelines, and they were approved by the Ethical Committee for Animal Testing (ECD 2014-92) at the University of Antwerp (Belgium) and all applicable institutional and European guidelines for the care and use of animals were followed.
To assess noninvasive quantification of [11C]ABP688, HET and WT Q175 mice (n = 6 per genotype, 6 months old) were scanned with an AV shunt in order to determine VT invasively and compare it to noninvasive approaches, namely VT with IDIF and non-displaceable binding potential (BPND) calculated from distribution volume ratio-1 (DVR-1) as well as using the simplified reference tissue model (SRTM). Comparison of discrimination power between genotypes was performed for the different quantification methods. To validate cerebellum as reference region, a blocking experiment with an mGluR5 antagonist was performed and compared to baseline scans (n = 4 per condition, 9 months old). In order to determine the accuracy of a [11C]ABP688 PET template for quantifying mGluR5 the V in different regions was quantified in a total of 36 mice (HET and WT Q175, n = 18 per genotype, 6 months old) following two different spatial normalization approaches based respectively on individual MR images or the aforementioned [11C]ABP688 PET template. Finally, for test-retest stability of [11C]ABP688, a total of 10 mice were included (HET and WT Q715, n = 5 for each genotype, 6 months old).
[11C]ABP688 was prepared using an automated synthesis module (Carbosynthon I, Comecer, The Netherlands). Synthesis of [11C]ABP688 was accomplished by reacting of 0.5 mg desmethyl-ABP-688 (E/Z) with [11C]CH3SO3CF3 in 400 μl of acetone in presence of 10 μl of NaOH, followed by purification and filtration as previously described (8). Average radiochemical purity was 98.6 ± 1.2%, while the mean specific radioactivity was 73 ± 16 GBq/μmol.
[11C]ABP688 Dynamic microPET Scan
MicroPET/Computed tomography (CT) imaging was performed on two Siemens Inveon PET-CT scanners (Siemens Preclinical Solution, USA). The animals were anesthetized using isoflurane (Forene, Belgium) in medical oxygen (induction 5%, maintenance 1.5%), catheterized in the tail vein for intravenous (i.v) bolus injection of the tracer and positioned onto the scanner. Respiration and heart rate of the animal were constantly monitored using the Monitoring Acquisition Module (Minerve, France), with body temperature of the animals maintained at 37 ± 1°C using a feedback-controlled warm air flow (Minerve, France) during the entire scanning period. A full body image was acquired in a single PET bed position, thus including the lumen of the left ventricle of the heart in the field of view (FOV) for the calculation of the IDIF.
To measure the arterial input function, an AV shunt was surgically inserted into the femoral vein and artery prior to PET scan. After positioning the animal onto the scanner, the shunt was connected to a peristaltic pump: tubing from the artery was led through the Twilite detector (23) and ran through the pump. The tubing coming from the vein was connected on the output line of the pump, together with a second line for tracer injection with a bolus of [11C]ABP688. During injection, the peristaltic pump was stopped to prevent backflow. To perform [11C]ABP688 dynamic microPET scan with cardiac gating, cardiac electrocardiogram (ECG) signal were obtained from electrode tubes covered with ECG gel that were placed over the animal's front legs and one hind leg.
At the onset of the 60 min dynamic microPET scan, mice were injected with a bolus of [11C]ABP688 over a 12 s interval (1 ml/min) using an automated pump (Pump 11 Elite, Harvard Apparatus, USA). Tracer activity was injected keeping the cold dose within tracer conditions (< 1.50 μg/kg). PET data were acquired in list mode and the cardiac gate trigger signals were inserted into the list mode stream. Following the microPET scan, a 10 min 80 kV/500 μA CT scan was performed for attenuation and scatter correction. The AV shunt surgery was performed in the animals allocated for the first study (invasive quantification of [11C]ABP688), while the cardiac gating was included in the test-retest study. For the blocking experiment, the highly selective non-competitive mGluR5 antagonist 1,2-methyl-6-(phenylethynyl)-pyridine (MPEP) was used. MPEP was dissolved in saline and administered with an i.v. bolus injection (6 mg/kg) 10 minutes before the injection of the radiotracer. To limit the possible effect of a circadian variation in mGluR5 availability, test-retest scans were acquired at similar time of the day, while the scans for the blocking experiment were performed at the exact same time of the day.
A total of 3 WT and 2 HET Q175 mice were excluded from the [11C]ABP688 PET template validation study due to issues related to either tracer injection or image acquisition. Data on the body weight of the animals, injected radioactivity, injected mass, and number of animals for each study are reported in Supplementary Table 1.
Image Processing and Analysis
Acquired PET data were histogrammed and reconstructed into 33 frames of increasing duration (12 × 10 s, 3 × 20 s, 3 × 30 s, 3 × 60 s, 3 × 150 s, and 9 × 300 s). Iterative PET image reconstruction of the images was performed using 4 iterations and 16 subsets of the 2-dimensional ordered-subset expectation maximization (2D-OSEM) algorithm (24) following Fourier rebinning. Normalization, dead time, CT-based attenuation and single-scatter simulation scatter corrections were applied. PET image frames were reconstructed on a 128 × 128 × 159 grid with 0.776 × 0.776 × 0.776 mm3 ignoring the cardiac gating trigger signals. These reconstructions were used for quantifying the brain uptake. Additional cardiac gated reconstruction was also obtained by dividing each heart cycle (defined by the gating trigger signals in the list-mode stream) into 4 bins and the bin containing the diastole images were used to derive the IDIF.
The arterial IPF was obtained at a 1 s sampling interval from the whole blood activity derived from the Twilite count detection coupled with the AV shunt. To reduce the noise in the Twilite data, a three-exponential function was fitted to the decaying part of the IPF. The delay between the IPF measured at the shunt and the true cerebral IPF was estimated trough a two-tissue compartmental model fit with an extra free parameter for the delay to the striatal and cortical time activity curves (TACs). The same delay was shared between the different regions. The IDIF was obtained from the whole blood activity derived from the PET images by delineating a region-of-interest (threshold set to 50% of max) in the lumen of the left ventricle of the heart. The ventricular region was delineated on an early time frame exhibiting maximal activity in the lumen of the left ventricle.
PET images were processed and analyzed using PMOD 3.6 software (Pmod Technologies, Zurich, Switzerland). For spatial normalization of the PET/CT images, a [11C]ABP688 PET template was generated using only the data of the WT animals (n = 16). First, individual static PET images covering the whole scan duration (i.e., 60 min) were generated and spatially aligned to their individual MRI by applying the non-linear warping CT to MRI transformation. The transformation was calculated by co-registering the individual animals CT image to its corresponding MR image. The same transformation could be used to align the PET images as the PET/CT images were intrinsically co-registered as acquired on the same gantry. Then, all individual MR images were spatially transformed through a rigid body registration to the MRI of the first animal. The MR to MR transformations obtained during creation of the WT MR template were then also applied to the corresponding static PET images. The static PET MR images (n = 16) were averaged resulting in a [11C]ABP688 WT template ([11C]ABP688 PET template). To validate the [11C]ABP688 PET template, spatial normalization of the PET/CT images was performed through rigid body image co-registration of the PET images to both (i) individual MR images (based on the CT to MR transformation) and (ii) [11C]ABP688 PET template (based on a transformation calculated from PET to PET image registration). Following validation, the [11C]ABP688 PET template was applied to all other PET images for quantification (i.e., test-retest study, validation IDIF, and DVR-1 vs. SRTM comparison). The volumes-of-interest (VOIs) were determined based on an existing MRI template with predefined VOIs (25). Using the predefined VOIs of the template, standardized uptake value (SUV) TACs of different regions (striatum, cerebral cortex, hippocampus, thalamus, and cerebellum) were extracted from the images. Following kinetic modeling, TACs were fitted by a standard two-tissue compartment model (2TCM) with blood volume fixed at 0.036 mL/cm3 and by a Logan model (26) to calculate the total volume of distribution VT using the measured arterial IPF and IDIF (V and V). As there was a high correlation between 2TCM and Logan (e.g., in striatum: r = 0.998, r2 = 0.997), the regression line was close to identity line, and the Bland-Altman plot showing limited bias (4.09%) as well as negligible 95% confidence intervals (from 1.81 to 6.38%), only VT calculated with Logan is reported. From the VT values, DVR-1 was calculated with the cerebellum as reference region. In addition, the binding potential BPND for these regions was calculated using the SRTM (14) with the cerebellum as reference tissue.
VT, DVR-1 and BPND values were calculated to validate IDIF as input function for kinetic modeling. VT values were calculated for the blocking experiment to validate the presence of reference region. VT values obtained from either the individual MRI or PET template approach, and they were compared to evaluate the accuracy of the [11C]ABP688 PET template for spatial normalization. Finally, VT and BPND values were calculated to determine test-retest stability of [11C]ABP688 PET imaging.
Additionally, voxel-based parametric VT images were generated using Logan model with the IDIF as input function, while BPND images were calculated using the SRTM with the cerebellum as reference region. Parametric images are represented as averages over the group (HET and WT) in stack coronal slices selected from a 3D coronal/sagittal/transversal mouse brain view.
To validate the [11C]ABP688 template for PET quantification, individual MR images were obtained in the same week of the microPET/CT scan. The animals were anesthetized using isoflurane in a mixture of N2/O2 (induction 5%, maintenance 1.5%) and placed in prone position onto the scanner (7T Biospec, Bruker, Germany). Body temperature was maintained at 37 ± 1°C by means of rectal thermistor with a feedback-controlled warm air circuitry (MR-compatible Small Animal Heating System, SA Instruments, Inc. USA). Three-dimensional turbo rapid acquisition with relaxation enhancement (turboRARE) images were acquired with repetition time 3,185 ms, echo time 44 ms, echo train length 8, and matrix size 128 × 64 × 40. FOV was 25.6 × 13 × 10 mm3 and resolution of 0.2 × 0.2 × 0.25 mm3. The MR image acquisition procedure lasted 21 min. Data were acquired using ParaVision 5.1 (Bruker, Germany).
All data were assessed for normality (Shapiro-Wilk test). Since no evidence against normality was found, parametric tests were performed. Differences between VT quantification based on either AV shunt or IDIF were evaluated separately with the paired-t test for each region. Two-way ANOVA was applied to the blocking experiment to compare baseline and blockade scans in the different regions. Pearson's correlation tests were used to determine the relationship between VT and DVR-1 quantified with the different input functions (namely, AV shunt and IDIF), between DVR-1 and BPND, and to determine correlations between VT values obtained with both the individual MRI and PET template. Additionally, the agreement between individual VT measurements based on MRI and PET templates was visualized by plotting the percentage difference between the two parameters against their averages in a Bland-Altman plot (27). For the test-retest study, reproducibility of the data was determined by the intraclass correlation coefficient (ICC). Pearson's correlation tests as well as Bland-Altman plots were used to compare VT and BPND values for test-retest scans. In addition, the percentage relative change between test and retest were calculated as relative difference = |retest – test|/retest x 100%. Finally, the mean ± standard deviation (SD) of the intra-animal coefficient of variation (COV) was calculated as follow:
where G represents the group, N is the number of animals in the group, and are respectively the mean and standard deviation of the test and retest values for animal i. Paired t-test was performed to investigate any methodological difference between test and retest scans. All analyses were performed with GraphPad Prism (v 6.0) statistical software, with the exception of the ICC, which were calculated in JMP Pro 13 software (SAS Institute Inc., USA) and the power analysis, which was calculated with G*Power software (http://www.gpower.hhu.de/). The data are represented as mean ± standard deviation (SD), unless specified otherwise. All tests were two-tailed and significance was set at p < 0.05.
Quantification Based on IDIF and IPF Correlate Significantly When Using Cerebellum as the Reference Region
To evaluate a noninvasive approach for quantification of [11C]ABP688, arterial IPF based on AV shunt and IDIF were compared. Average input function TACs for AV shunt and IDIF approaches are shown in Figure 1A. The IDIF tail values were higher than the AV shunt values at corresponding time points, indicating that IDIF overestimates the radioactivity present in the blood (Figure 1A). Accordingly, VT values using IDIF are significantly lower (p < 0.0001) than with an AV shunt IPF in both WT and HET Q175 mice (n = 6 per genotype; Figure 1B). As shown in Figure 2A and reported in Supplementary Table 2, moderate correlations were found between VT values obtained from the 2 input functions in striatum (r = 0.629, p = 0.028), cortex (r = 0.557, p = 0.059), hippocampus (r = 0.558, p = 0.059), and thalamus (r = 0.536, p = 0.072). However, when considering the cerebellum as reference region, correlations between DVR-1 based on AV shunt and IDIF were strong and highly significant in all regions (r = 0.978, p < 0.0001 in striatum; r = 0.967, p < 0.0001 in cortex; r = 0.955, p < 0.0001 in hippocampus; r = 0.934, p < 0.0001 in thalamus) (Figure 2B and Supplementary Table 2).
Figure 1. Comparison of invasive and noninvasive input functions for [11C]ABP688 quantification. (A) Average SUV TACs of the IPF based on arteriovenous (AV) shunt and image-derived input function (IDIF) (n = 12) (note the logarithmic scale on the y-axis). (B) Comparison of volume of distribution (VT) quantification of [11C]ABP688 based on Logan plot in striatum (left) and cortex (right). As a consequence of the overestimation of the tail of the IPF by the IDIF, the VT is significantly lower compared to VT based on AV shunt IPF. Paired-t test (n = 12). ****p < 0.0001. IPF, input function; WT, wild type; HET, heterozygous.
Figure 2. Correlations between invasive and noninvasive [11C]ABP688 quantification. (A) Volume of distribution (VT) values showed a moderate correlation when comparing the 2 input function approaches in both striatum (left) and cortex (right). (B) Distribution volume ratio (DVR-1) of [11C]ABP688 displayed strong correlations between input function based on arteriovenous shunt and image-derived input function in both striatum (left) and cortex (right). Pearson's correlation test.
A blocking experiment was performed to investigate cerebellum as a possible reference region in mice (Figure 3). Average standardized uptake value (SUV) time-activity curves of animals injected with MPEP displayed a clear reduction of [11C]ABP688 binding compared to baseline in striatum, cortex, and hippocampus, while in cerebellum no reduction could be observed (Figure 3B). Accordingly, MPEP administration resulted in a statistically significant decrease of [11C]ABP688 binding in receptor-rich regions (striatum: −42%, p < 0.0001; cortex: −40%, p < 0.0001; hippocampus: −33%, p < 0.0001). On the contrary, in the cerebellum only a negligible no significant reduction of V (Logan) could be observed (−4.8%, p > 0.99) (Figure 3C).
Figure 3. Validation of cerebellum as reference region for [11C]ABP688 in mice. (A) Average parametric VT images during baseline scan (BL) and following administration of MPEP (6 mg/kg) 10 min before tracer injection in WT mice. (B) SUV time-activity curves for different brain region during baseline and blocking scans. (C) VT (Logan) quantification showed a statistically significant reduction in mGluR5-rich regions, while no significant changes were found in cerebellum. n = 4 per condition. ****p < 0.0001. BL, baseline; MPEP, 1,2-methyl-6-(phenylethynyl)-pyridine; STR, striatum; CTX, cortex; HIP, hippocampus; CB, cerebellum.
As the administration of mGluR5 antagonist MPEP did not modify V values in cerebellum, we selected this region as reference region for relative quantification. Thus, we quantified BPND with SRTM [BPND(SRTM)] and compared it to the aforementioned DVR-1 values based on AV shunt and IDIF. Correlations between BPND(SRTM) and DVR-1 based on AV shunt were strong and highly significant in all the investigated regions (r = 0.987, p < 0.0001 in striatum; r = 0.979, p < 0.0001 in cortex; r = 0.976, p < 0.0001 in hippocampus; r = 0.953, p < 0.0001 in thalamus) (Figure 4 and Supplementary Table 3). Similarly, highly significant correlations were found between BPND(SRTM) and DVR-1 based on IDIF in all the investigated regions (r = 0.980, p < 0.0001 in striatum; r = 0.953, p < 0.0001 in cortex; r = 0.959, p < 0.0001 in hippocampus; r = 0.953, p < 0.0001 in thalamus) (Figure 4 and Supplementary Table 3).
Figure 4. Correlations between distribution volume ratio (DVR-1) and binding potential (BPND) values for [11C]ABP688. In striatum (A) and cortex (B), DVR-1 calculated with AV shunt as well as IDIF strongly correlated with BPND determined using the simplified reference tissue model (SRTM). Pearson's correlation test. AV, arteriovenous; IDIF, image-derived input function.
Since both V, V, and BPND(SRTM) were obtained from the same animals, we investigated whether the use of an IDIF or BPND would affect the quantification of disease-related changes in WT and HET Q175 mice. As shown in Supplementary Figure 1, the genotypic difference in striatal V and V displayed was comparable (−7.0%, p = 0.40 and −5.9%, p = 0.35, respectively), while striatal BPND(SRTM) could detect a larger difference between genotypes (−10.3%, p = 0.16), however none of the comparison reached statistical significance due to the limited sample size (n = 6 per genotype). In accordance, the power analysis performed on these values confirmed the higher discrimination power for BPND(SRTM) compared to V and V (Supplementary Table 4).
MRI and PET Template-Based Spatial Normalization Approaches Showed High Agreement in [11C]ABP688 Quantification
To evaluate the necessity of individual MR images for [11C]ABP688 quantification, we compared VT values based on individual MR images and the [11C]ABP688 PET template. The results of the different spatial normalization approaches for [11C]ABP688 are summarized in Figure 5 and Table 1. A strong and highly significant correlation was found between spatial normalization approaches in all investigated regions (r = 0.989, p < 0.0001 in striatum; r = 0.983, p < 0.0001 in cortex; r = 0.985, p < 0.0001 in hippocampus; r = 0.986, p < 0.0001 in thalamus) (Table 1). The correlations between approaches coincided with the identity line (Figure 5A). In addition, the Bland-Altman plot showed high agreement between the two approaches as visible by the low bias obtained (red dashed line; Figure 5B), which was ≤ 1% for all the investigated regions (Table 1).
Figure 5. Comparison of spatial normalization approaches for [11C]ABP688 VT regional quantification. (A) Correlation between VT values in striatum (left) and cortex (right) based on individual MR images- and PET template-based spatial normalizations. Pearson's correlation test. Dashed line represents identity line. (B) Bland-Altman plot to compare the spatial normalization approaches in striatum (left) and cortex (right). The bias between the two approaches corresponds to the difference between the mean (red dashed line) and X axis (solid line). The dotted lines represent the 95% limits of agreement (mean difference ± 1.96 × SD of the differences), respectively. WT, wild type; HET, heterozygous.
Noninvasive Quantification of [11C]ABP688 Imaging Is Reproducible
Methodologically no significant difference between test and retest scans was observed in the injected dose (5.24 ± 0.91 MBq and 4.89 ± 1.06 MBq, respectively; p = 0.39), specific activity (74.3 ± 18.5 GBq/μmol and 68.7 ± 14.2 GBq/μmol, respectively; p = 0.27), body weight (28.2 ± 1.8 g and 27.8 ± 1.6 g, respectively; p = 0.59) or injected mass (1.23 ± 0.08 μg/kg and 1.25 ± 0.05 μg/kg, respectively; p = 0.39). Finally, there was no significant difference between the start time for the test and retest scans (Δt = 59 ± 62 min; p = 0.62). Average TACs for striatum and cerebellum were consistent between test and retest scans and genotype as shown in Supplementary Figure 2.
Values and differences between test and retest scans are summarized in Table 2. For [11C]ABP688 V, the mean relative difference between test and retest was lower than 4.1% in all the investigated regions, with a SD ranging from 7.4 to 17.9%. ICC values were low in thalamus (0.35) and between 0.46 and 0.61 in the other regions (hippocampus, striatum and cortex, respectively). Intra-animal COV was lower than 10% in all investigated regions, with a SD ranging from 6.8 to 7.4%. In addition, the Bland-Altman plot showed a negligible bias between test and retest (0.30%), although the 95% confidence intervals were relatively large (−25.4%; 24.7%) (Figure 6).
Table 2. Reproducibility of test-retest parameters for [11C]ABP688 volume of distribution (VT) and binding potential (BPND).
Figure 6. Reproducibility of test-retest quantification for [11C]ABP688. Correlations of the volume of distribution (VT) (A) and binding potential (BPND) (B) for the regions of interest in both WT and HET Q175 mice (n = 5 per genotype). Pearson's correlation test. Bland-Altman plot to compare test-retest quantification of VT (C) and BPND (D). The bias between the two approaches corresponds to the difference between the mean (red dashed line) and X axis (solid line). The dotted lines represent the 95% limits of agreement (mean difference ± 1.96 x SD of the differences), respectively. (E) Average parametric images (n = 5) representing the VT (left) and BPND (right) for both WT and HET Q175 mice. STR, striatum; HIP, hippocampus; THAL, thalamus; CTX, cortex; WT, wild type; HET, heterozygous.
For [11C]ABP688 BPND(SRTM), the mean relative difference between test and retest was lower than 3.5% in all the investigated regions, with a SD ranging from 5.7 to 9.6%. ICC value was low in thalamus (0.32) and between 0.53 and 0.62 in the other regions (hippocampus, striatum and cortex, respectively). Intra-animal COV was lower than 10.2% in all investigated regions, with a SD ranging from 6.1 to 7.7%. Finally, the Bland-Altman plot confirmed the negligible bias between test and retest (−0.04%), despite the relatively large 95% confidence intervals (−29.6%; 29.5%) (Figure 6).
Scatter plots comparing the individual outcome values for both VT and BPND and average parametric images for both VT and BPND are shown in Figure 6.
The present study investigated several imaging parameters for reproducible quantification of [11C]ABP688 PET imaging in mice.
Image-Derived Input Function Reproducibility and Reference Region Validation for Noninvasive Quantification in Mice
The employment of a noninvasive input function is extremely attractive for longitudinal PET imaging as it circumvents all challenges and limitations related to the arterial blood sampling. Thus, in the present study, we compared the invasive arterial IPF measured with AV shunt to a proposed noninvasive IDIF measured in the left ventricle of the heart. The IDIF resulted in a broader peak than the AV shunt. This can be attributed to the coarser time sampling and the contributions from both the left and right ventricular blood pools as well as the ventricular walls, three regions for which the peak activity concentration does not occur at the same time. The proposed IDIF resulted in an activity peak comparable to the AV shunt IPF, however, the IDIF overestimated the blood activity at the tails of the curve, resulting in higher values than the AV shunt IPF at corresponding times. This is likely due to spill over activity from nearby regions (e.g., myocardium) which is negligible at the time of the peak but becomes significant at later times. Since the difference between IDIF and AV shunt IPF increased during the scan, it was not possible to simply scale the IDIF e.g., using a single blood sample. As a consequence of the overestimation of the tail, VT values calculated with IDIF were reduced compared to the values obtained with the AV shunt IPF (p < 0.0001). As a result, moderate correlations could be established between VT values determined with IDIF and AV shunt (i.e., r = 0.629 in striatum). Nonetheless, when comparing the VT values obtained from the 2 IPFs between WT and HET Q175 mice, where a reduction of mGluR5 binding is expected (28), the absolute genotypic VT difference was comparable (e.g., in striatum, AV shunt: −7.0%, p = 0.40; IDIF: −5.9%, p = 0.36). This suggests that values determined with the IDIF, although they do not directly match those obtained using the AV shunt, might be representative of the quantification based on AV shunt IPF. However, future studies are necessary to investigate whether noninvasive IDIF could be employed as alternative to evaluate noninvasively group differences in longitudinal or interventional studies.
A limitation of the present study was the lack of metabolite-corrected IPF. This was not performed because of the limited amount of blood can be obtained from a specific mouse, which does not allow to collect multiple samples from the same animal when using standard techniques for metabolite analysis. One way to circumvent this limitation in future studies could be the use of different cohorts of animals in order to generate a population-based correction curve, possibly taking into account a potential age-dependent metabolism of the radiotracer.
Previous clinical and preclinical studies have reported the application of the cerebellum or cerebellar gray matter as reference region for [11C]ABP688 quantification (10, 15, 19). Therefore, we calculated DVR-1 values based on VT from both IDIF and AV shunt using the cerebellum as reference region and we found strong and highly significant correlations between the 2 measurements (e.g., in striatum, r = 0.978; p < 0.0001). This finding indicates that DVR-1 quantification using IDIF provides highly comparable measurements to the arterial IPF and therefore it can be utilized for noninvasive quantification of [11C]ABP688 in the presence of a valid reference region. Nevertheless, the use of SRTM remains the preferable choice for noninvasive quantification of [11C]ABP688 if the reference region is validated. Validation of the reference region was performed with a blocking experiment by administration of the mGluR5 antagonist MPEP before the administration of the radiotracer. As the blocking of mGluR5 had a non-significant and negligible effect on cerebellar VT, we concluded that cerebellum can be used as reference region for [11C]ABP688 quantification. The 4.8% cerebellar reduction found in the present study is in line with the previously reported validation in rats where blocking with MPEP resulted in 7% reduction of the cerebellar VT (based on Logan graphical analysis) (10). Other studies in baboons and humans validated the use of cerebellar gray matter as reference region (15, 11). Although the use of cerebellar gray matter might be more accurate, in mice a clear distinction between these regions is not possible due to the limited resolution of the PET camera. Thus, we validated the reference tissue model with the cerebellum as reference region by comparing BPND values calculated using the SRTM to the DVR-1 calculated with IDIF as well as with the AV shunt. As this resulted in strong and highly significant correlations (e.g., in striatum, r = 0.980 and r = 0.987, respectively). Interestingly, BPND(SRTM) resulted in a higher discrimination between genotypes (e.g., in striatum, −10.3%, p = 0.16) than the VT quantification based on either AV shunt IPF or IDIF. Altogether these findings indicate that the SRTM using the cerebellum as reference region is suitable for quantification of [11C]ABP688 binding in mice.
[11C]ABP688 PET Template for Spatial Normalization
As a dedicated high resolution small animal MR scanner it is not always available to conduct preclinical PET studies and also in order to limit as much as possible the anesthesia sessions, we investigated the relevance of individual MRI-based spatial normalization to quantify [11C]ABP688. When comparing the MRI-based spatial normalization to the [11C]ABP688 PET template approach, we found that [11C]ABP688 VT values strongly correlated (e.g., in striatum, r = 0.989, p < 0.0001), with no statistical significant difference in the quantification of [11C]ABP688. This finding indicates that PET-based spatial normalization is comparable to the MRI-based spatial normalization and therefore the use of individual MR images is not essential to obtain reproducible [11C]ABP688 quantification as we previously also confirmed for rats (29). This could be related to the fairly high spatial information in [11C]ABP688 PET images by the abundancy of mGluR5 throughout the brain. However, it is important to note that a disease condition might have a detrimental effect on the accuracy of the PET template for spatial normalization due to an altered signal. Therefore, we included diseased HD mice as in HD the levels of mGluR5 are known to be altered (30, 31). Importantly, imaging HD mice did not invalidate the high agreement between the 2 spatial normalization approaches.
Test-Retest Stability of [11C]ABP688 PET Imaging in Mice
As reproducible quantification is fundamental in order to perform interventional and longitudinal studies, we examined the test-retest stability of [11C]ABP688 PET in healthy and diseased mice. Both BPND and VT values showed a very low relative group difference between test and retest (< 5%) as well as negligible bias with Bland-Altman plots. Even though the average values between test and retest were highly reproducible, standard deviations of these measurements were relatively large, indicating high reproducibility on a group-level but limitations for the single animal. Nonetheless, these findings are in line with previous test-retest studies in rats and baboons, where a percentage difference below 10% was reported (11, 32, 18), with no statistically significant changes between the 2 scans. On the contrary, in a test-retest study in humans, a significant increase in BPND values was reported in the retest measurement, with regional differences up to 73% (15, 20). Recent studies confirmed test-retest regional increased uptake of [11C]ABP688 during retest scan in rhesus monkeys (19). This increased [11C]ABP688 binding during the retest scan could be related to a circadian variation of the mGluR5 during the day. In line with this hypothesis, a previous study comparing sleep deprived subjects to controls reported increased [11C]ABP688 binding following sleep deprivation (33). Additionally, a recent report showed that [11C]ABP688 binding in rats is significantly different during distinct phases of the day, supporting the hypothesis that glutamate binding to mGluR5 undergoes a circadian variation (34).
In conclusion, we proposed a cardiac noninvasive input function (IDIF) to quantify the [11C]ABP688 volume of distribution in mice as well as the use of a cerebellar reference region and the SRTM method for the relative quantification of mGluR5 in mice.
Additionally, the good agreement between spatial normalization approaches indicates that a [11C]ABP688 PET template can be used for reproducible regional quantification during interventional or longitudinal studies in the mouse brain. Taken together, our results indicate that noninvasive relative quantification of [11C]ABP688 PET imaging can be perform in mice.
JV, SSta, and LM conceived and designed the study. JV, DB, LK, DT, LW, AV, SStr, JW, CD, LM, and SSta were involved in execution of the experimental design, data acquisition and data interpretation for the study. DB, JV, and SSta were involved in drafting and editing the manuscript and figures. All authors approved the final manuscript and they are accountable for the content of the work.
This work was funded by CHDI Foundation, Inc., a nonprofit biomedical research organization exclusively dedicated to developing therapeutics that will substantially improve the lives of HD-affected individuals.
Conflict of Interest Statement
This work was funded by CHDI Foundation, Inc., a nonprofit biomedical research organization exclusively dedicated to developing therapeutics that will substantially improve the lives of HD-affected individuals.
DB has a PhD fellowship from the Research Foundation Flanders (FWO, 11W2516N/11W2518N). The authors thank Philippe Joye and Caroline Berghmans of the Molecular Imaging Center Antwerp (MICA) for their expertise and valuable technical assistance.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fneur.2018.00516/full#supplementary-material
Supplementary Figure 1. Comparison of [11C]ABP688 striatal quantification in WT and HET Q175 mice (n = 6 per genotype) using volume of distribution (VT) based on Logan plot with AV shunt (inv) (A) and IDIF (B) as well as binding potential (BPND) using SRTM (C). WT, wild type; HET, heterozygous; AV, arteriovenous; IDIF, image-derived input function; SRTM, simplified reference tissue model; p, p-value. Data are represented as mean ± standard deviation.
Supplementary Figure 2. Average SUV TACs for the test and retest scans of WT and HET Q175 mice (n = 5 per genotype). STR, striatum; CB, cerebellum; WT, wild type; HET, heterozygous. Data are represented as mean ± standard error mean.
Supplementary Table 1. Scan parameters for each of the performed studies. AV, arteriovenous; IDIF, image-derived input function; WT, wild type; HET, heterozygous.
Supplementary Table 2. Relationship between VT and DVR-1 [11C]ABP688 quantification based on invasive (AV shunt) and noninvasive input function (IDIF). AV, arteriovenous; IDIF, image-derived input function; DVR, distribution volume ratio; r, Pearson's correlation; r2, coefficient of determination.
Supplementary Table 3. Correlations between DVR-1 and BPND for [11C]ABP688. AV, arteriovenous; IDIF, image-derived input function; DVR, distribution volume ratio; BPND, binding potential; SRTM, simplified reference tissue model; r, Pearson's correlation; r2, coefficient of determination.
Supplementary Table 4. Representative power analysis to detect striatal difference in [11C]ABP688 quantification between WT and HET Q175 mice at 6 months of age.
2. Ribeiro FM, Paquet M, Cregan SP, Ferguson SS. Group I metabotropic glutamate receptor signalling and its implication in neurological disease. CNS Neurol Disord Drug Targets (2010) 9:574–95. doi: 10.2174/187152710793361612
5. Vonsattel JP, Myers RH, Stevens TJ, Ferrante RJ, Bird ED, Richardson EP. Jr. Neuropathological classification of Huntington's disease. J Neuropathol Exp Neurol. (1985) 44:559–77. doi: 10.1097/00005072-198511000-00003
6. Joel D, Weiner I. The organization of the basal ganglia-thalamocortical circuits: open interconnected rather than closed segregated. Neuroscience (1994) 63:363–79. doi: 10.1016/0306-4522(94)90536-3
8. Ametamey SM, Kessler LJ, Honer M, Wyss MT, Buck A, Hintermann S, et al. Radiosynthesis and preclinical evaluation of 11C-ABP688 as a probe for imaging the metabotropic glutamate receptor subtype 5. J Nucl Med. (2006) 47:698–705. Available online at: http://jnm.snmjournals.org/content/47/4/698.full.pdf+html
9. Treyer V, Streffer J, Wyss MT, Bettio A, Ametamey SM, Fischer U, et al. Evaluation of the metabotropic glutamate receptor subtype 5 using PET and 611C-ABP688: assessment of methods. J Nucl Med. (2007) 48:1207–15. doi: 10.2967/jnumed.107.039578
10. Elmenhorst D, Minuzzi L, Aliaga A, Rowley J, Massarweh G, Diksic M, et al. In vivo and in vitro validation of reference tissue models for the mGluR(5) ligand [11C]ABP688. J Cereb Blood Flow Metab. (2010) 30:1538–49. doi: 10.1038/jcbfm.2010.65
11. Delorenzo C, Milak MS, Brennan KG, Kumar JS, Mann JJ, Parsey RV. In vivo positron emission tomography imaging with [11C]ABP688: binding variability and specificity for the metabotropic glutamate receptor subtype 5 in baboons. Eur J Nucl Med Mol Imaging (2011) 38:1083–94. doi: 10.1007/s00259-010-1723-7
13. Locke LW, Berr SS, Kundu BK. Image-derived input function from cardiac gated maximum a posteriori reconstructed PET images in mice. Mol Imaging Biol. (2011) 13:342–7. doi: 10.1007/s11307-010-0347-4
15. Delorenzo C, Kumar JS, Mann JJ, Parsey RV. In vivo variation in metabotropic glutamate receptor subtype 5 binding using positron emission tomography and [11C]ABP688. J Cereb Blood Flow Metab. (2011) 31:2169–80. doi: 10.1038/jcbfm.2011.105
16. Romano C, Sesma MA, Mcdonald CT, O'malley K, Van Den Pol AN Olney JW. Distribution of metabotropic glutamate receptor mGluR5 immunoreactivity in rat brain. J Comp Neurol. (1995) 355:455–69. doi: 10.1002/cne.903550310
17. Neki A, Ohishi H, Kaneko T, Shigemoto R, Nakanishi S, Mizuno N. Metabotropic glutamate receptors mGluR2 and mGluR5 are expressed in two non-overlapping populations of Golgi cells in the rat cerebellum. Neuroscience (1996) 75:815–26. doi: 10.1016/0306-4522(96)00316-8
18. Elmenhorst D, Aliaga A, Bauer A, Rosa-Neto P. Test-retest stability of cerebral mGluR(5) quantification using [11C]ABP688 and positron emission tomography in rats. Synapse (2012) 66:552–60. doi: 10.1002/syn.21542
19. Sandiego CM, Nabulsi N, Lin SF, Labaree D, Najafzadeh S, Huang Y, et al. Studies of the metabotropic glutamate receptor 5 radioligand [11C]ABP688 with N-acetylcysteine challenge in rhesus monkeys. Synapse (2013) 67:489–501. doi: 10.1002/syn.21656
20. Delorenzo C, Gallezot JD, Gardus J, Yang J, Planeta B, Nabulsi N, et al. In vivo variation in same-day estimates of metabotropic glutamate receptor subtype 5 binding using [11C]ABP688 and [18F]FPEB. J Cereb Blood Flow Metab. (2016) 37:2716–27. doi: 10.1177/0271678X16673646
21. Heikkinen T, Lehtimaki K, Vartiainen N, Puolivali J, Hendricks SJ, Glaser JR, et al. Characterization of neurophysiological and behavioral changes, MRI brain volumetry and 1H MRS in zQ175 knock-in mouse model of Huntington's disease. PLoS ONE (2012) 7:e50717. doi: 10.1371/journal.pone.0050717
22. Menalled LB, Kudwa AE, Miller S, Fitzpatrick J, Watson-Johnson J, Keating N, et al. Comprehensive behavioral and molecular characterization of a new knock-in mouse model of Huntington's disease: zQ175. PLoS ONE (2012) 7:e49838. doi: 10.1371/journal.pone.0049838
23. Alf MF, Wyss MT, Buck A, Weber B, Schibli R, Kramer SD. Quantification of brain glucose metabolism by 18F-FDG PET with real-time arterial and image-derived input function in mice. J Nucl Med. (2013) 54:132–8. doi: 10.2967/jnumed.112.107474
25. Mirrione MM, Schiffer WK, Fowler JS, Alexoff DL, Dewey SL, Tsirka SE. A novel approach for imaging brain-behavior relationships in mice reveals unexpected metabolic patterns during seizures in the absence of tissue plasminogen activator. Neuroimage (2007) 38:34–42. doi: 10.1016/j.neuroimage.2007.06.032
26. Logan J, Fowler JS, Volkow ND, Wolf AP, Dewey SL, Schlyer DJ, et al. Graphical analysis of reversible radioligand binding from time-activity measurements applied to [N-11C-methyl]-(-)-cocaine PET studies in human subjects. J Cereb Blood Flow Metab. (1990) 10:740–7. doi: 10.1038/jcbfm.1990.127
28. Abd-Elrahman KS, Hamilton A, Hutchinson SR, Liu F, Russell RC, Ferguson SSG. mGluR5 antagonism increases autophagy and prevents disease progression in the zQ175 mouse model of Huntington's disease. Sci Signal (2017) 10:eaan6387. doi: 10.1126/scisignal.aan6387
29. Verhaeghe J, Wyffels L, Wyckhuys T, Stroobants S, Staelens S. Rat brain normalization templates for robust regional analysis of [11C]ABP688 positron emission tomography/computed tomography. Mol Imaging (2014) 13:1–14. doi: 10.2310/7290.2014.00037
31. Gulyas B, Sovago J, Gomez-Mancilla B, Jia Z, Szigeti C, Gulya K, et al. Decrease of mGluR5 receptor density goes parallel with changes in enkephalin and substance P immunoreactivity in Huntington's disease: a preliminary investigation in the postmortem human brain. Brain Struct Funct. (2015) 220:3043–51. doi: 10.1007/s00429-014-0812-y
32. Miyake N, Skinbjerg M, Easwaramoorthy B, Kumar D, Girgis RR, Xu X, et al. Imaging changes in glutamate transmission in vivo with the metabotropic glutamate receptor 5 tracer [11C] ABP688 and N-acetylcysteine challenge. Biol Psychiatry (2011) 69:822–4. doi: 10.1016/j.biopsych.2010.12.023
33. Hefti K, Holst SC, Sovago J, Bachmann V, Buck A, Ametamey SM, et al. Increased metabotropic glutamate receptor subtype 5 availability in human brain after one night without sleep. Biol Psychiatry (2013) 73:161–8. doi: 10.1016/j.biopsych.2012.07.030
Keywords: [11C]ABP688, input function, mGluR5, PET imaging, test-retest
Citation: Verhaeghe J, Bertoglio D, Kosten L, Thomae D, Verhoye M, Van Der Linden A, Wyffels L, Stroobants S, Wityak J, Dominguez C, Mrzljak L and Staelens S (2018) Noninvasive Relative Quantification of [11C]ABP688 PET Imaging in Mice Versus an Input Function Measured Over an Arteriovenous Shunt. Front. Neurol. 9:516. doi: 10.3389/fneur.2018.00516
Received: 26 March 2018; Accepted: 11 June 2018;
Published: 29 June 2018.
Edited by:Freimut Dankwart Juengling, St. Claraspital Basel, Switzerland
Reviewed by:Christine DeLorenzo, Stony Brook University, United States
Adriaan Anthonius Lammertsma, VU University Medical Center, Netherlands
Copyright © 2018 Verhaeghe, Bertoglio, Kosten, Thomae, Verhoye, Van Der Linden, Wyffels, Stroobants, Wityak, Dominguez, Mrzljak and Staelens. 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 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: Steven Staelens, email@example.com