Original Research ARTICLE
Washout allometric reference method (WARM) for parametric analysis of [11C]PIB in human brains
- 1Department of Nuclear Medicine and PET Centre, Aarhus University Hospital, Aarhus, Denmark
- 2Department of Neurology, Aarhus University Hospital, Aarhus, Denmark
- 3Department of Neuroscience and Pharmacology, University of Copenhagen, Copenhagen, Denmark
Rapid clearance and disappearance of a tracer from the circulation challenges the determination of the tracer's binding potentials in brain (BPND) by positron emission tomography (PET). This is the case for the analysis of the binding of radiolabeled [11C]Pittsburgh Compound B ([11C]PIB) to amyloid-β (Aβ) plaques in brain of patients with Alzheimer's disease (AD). To resolve the issue of rapid clearance from the circulation, we here introduce the flow-independent Washout Allometric Reference Method (WARM) for the analysis of washout and binding of [11C]PIB in two groups of human subjects, healthy aged control subjects (HC), and patients suffering from AD, and we compare the results to the outcome of two conventional analysis methods. We also use the rapid initial clearance to obtain a surrogate measure of the rate of cerebral blood flow (CBF), as well as a method of identifying a suitable reference region directly from the [11C]PIB signal. The difference of average absolute CBF values between the AD and HC groups was highly significant (P < 0.003). The CBF measures were not significantly different between the groups when normalized to cerebellar gray matter flow. Thus, when flow differences confound conventional measures of [11C]PIB binding, the separate estimates of CBF and BPND provide additional information about possible AD. The results demonstrate the importance of data-driven estimation of CBF and BPND, as well as reference region detection from the [11C]PIB signal. We conclude that the WARM method yields stable measures of BPND with relative ease, using only integration for noise reduction and no model regression. The method accounts for relative flow differences in the brain tissue and yields a calibrated measure of absolute CBF directly from the [11C]PIB signal. Compared to conventional methods, WARM optimizes the Aβ plaque load discrimination between patients with AD and healthy controls (P = 0.009).
The marker [11C]PIB is a radiotracer of amyloid-β (Aβ) plaque load, used clinically to diagnose Alzheimer's disease (AD) (Cohen et al., 2012). Despite its wide use, several challenges face the actual quantification of this and other Aβ probes (Ducharme et al., 2013), as recently reviewed by Kepe et al. (2013). Unfortunately, no dose-escalation and competitive displacement studies have been performed for any Aβ imaging agent. Thus, the in vivo binding properties of these tracers have not been fully characterized (Villemagne et al., 2012), including the differences of binding to specific tissue types (Fodero-Tavoletti et al., 2009) such as estrogen-related receptors. Functionally, AD syndromes are associated with degeneration of specific functional networks, and amyloid deposition as measured with [11C]PIB explains at most a small amount of the clinico-anatomic heterogeneity in AD (Lehmann et al., 2013). Even if [11C]PIB binds uniquely to extracellular Aβ plaque formations, the relevance of such binding may be coupled to autophagy as recent evidence suggests that Aβ secretion and plaque formation depend on autophagy (Nilsson et al., 2013). Given its widespread clinical use, enigmatic binding properties, and the relative high cost of positron emission tomographic procedures, better methods for the quantification of this tracer in relation to other clinical parameters adds to the clinical relevance of the tracer.
The rapid initial clearance from the circulation suggests that [11C]PIB is subject to flow-limited uptake, such that the parametric analysis of the binding must take blood flow and washout kinetics into account. There is therefore a need to couple Aβ measures to both metabolic function and CBF deficits (Johannsen et al., 2000), as well as the loss of flow variability observed in AD (Rodell et al., 2012).
Three methodological problems complicate the correct quantification of binding of flow limited and rapidly metabolising tracers like [11C]PIB in human brain (Møller et al., 2009).
The first is the rapid disappearance of the tracer from the circulation and the consequent brief exchange with brain tissue. Conventional binding studies by equilibrium methods are affected by this absence of a continuing source of tracer in the circulation. Reference and binding regions independently clear tracer from the respective volumes of distribution, and the tracer in the reference region no longer is a proper surrogate for the tracer in the circulation. The time-activity functions of different regions now depend on regional properties of binding, blood flow, and blood-brain barrier permeability, rather than on a common source of tracer in the circulation.
The second is the uncertain choice of a region of reference of specific flow and no specific binding. The cerebellum is held to be little affected by amyloid deposition in AD, and the cerebellar gray matter often serves as reference region for measures of [11C]PIB retention (Landau et al., 2013), but other regions have been considered as well, including whole cerebellum (Joshi et al., 2012), and the Pons (Klunk et al., 2007; Knight et al., 2011).
The third potential pitfall is the influence of cerebral blood flow and blood-brain barrier permeability differences on regionally specific binding of [11C]PIB. If the tracer clears from multiple compartments in a single region of interest with different quantities of exchangeable and bound tracer, it is possible that both flow and permeability differences can mimic or mask changes of binding.
The problems facing conventional [11C]PIB quantification are manifest when [11C]PIB retention is evaluated in healthy subjects and patients with AD with different degrees of Aβ deposition. For example, it has been reported that PIB may bind differentially to polymorphic Aβ aggregates in some humans as well as in animals (Rosen et al., 2010; Ikonomovic et al., 2012). Additionally, apparent retention of PIB is evident in cerebral white matter both in vivo by PET (Fodero-Tavoletti et al., 2009) and in vitro by postmortem auto-radiography (Klunk et al., 2004; Svedberg et al., 2009).
In the present study we aimed to establish a method of global parametric mapping of the binding potential (BPND) of [11C]PIB that would take the particular kinetic properties of [11C]PIB exchange with brain tissue into account. For evaluation of these properties we compared tree non-invasive assays of [11C]PIB binding in brain of healthy subjects and a group of patients suffering from AD.
First, we applied a new model of tracer clearance, Washout allometric reference method (WARM), which we designed to map the washout of the tracer from regions of specific and non-specific binding. Previous approaches to the analysis of wash-out of tracers included the early “Height-Over-Area” method of Zierler (1965) and the recent “Hypotime” method of Møller et al. (2009), from both of which the current approach borrows, taking into account the methodological weaknesses discussed by Kanno and Uemura (1975).
Second, measures of regional tracer clearances were used to identify a reference region of negligible specific binding.
Third, the simplified reference region method (SRTM) (Lammertsma and Hume, 1996), also accounts for flow-dependent differences in tracer delivery, but contrary to the WARM method it assumes that differential equations modeling the radioactivity in a region of interest and a reference region are continuously coupled by exchange with a well defined arterial contribution of tracer (Lammertsma and Hume, 1996).
Fourth, the reference region subsequently served to obtain parametric maps of BPND by means of a well established and clinically popular and simple ratio measure (SUVR) of the area under the specific retention curve (AUC), at a presumed optimal time range 40–60 min after i.v. administration of the tracer, relative to the AUC of the retention in the reference region.
Altogether, we evaluated the methods that are used to establish binding potentials, as well as the effect of flow-dependent correction on the binding potential of [11C]PIB.
Six patients with AD (four women and two men) with an average age of 65 (SD = 7) years and moderately reduced Mini-Mental State Examination (MMSE) scores of 22–25 volunteered to complete the tomography. The patients were recruited from the local Dementia Clinic and screened by an experienced neurologist to fulfill the criteria for probable Alzheimer s disease.
Eight healthy age-matched HC volunteers with a mean age of 68 (SD = 5) recruited by public advertisement served as controls. They all had a normal physical and neurological examination and had a MMSE between 28 and 30. To exclude cognitive impairment they were furthermore examined with the Danish version of CAMCOG (Lolk et al., 2000).
We obtained written informed consent from all subjects to the protocols approved by the Regional Science Ethics Committee in accordance with the Declarations of Helsinki. We previously reported some PET results from the same subjects (Rodell et al., 2012; Gjedde et al., 2013).
Positron Emission Tomography
All subjects had positron emission recordings, one or two with [15O]water and one with [11C]PIB, in the 3D mode of the ECAT High Resolution Research Tomograph (HRRT, CTI/Siemens, Knoxville, TN, USA) in a quiet room with the subjects resting in a supine position with eyes open. One of the male AD patients only completed the [11C]PIB recording. The images were reconstructed with 3D-OP-OSEM point spread function reconstruction (Varrone et al., 2009) using 10 iterations and 16 subsets with FWHM at approximately 1.5 mm. The reconstructed images were corrected for random and scatter events, detector efficiency variations, and dead time. Tissue attenuation scans were performed using a rotating 68Ge source. Dynamic emission recordings lasting 3 min (21 frames) were initiated upon bolus intravenous injection of [15O]water (500 MBq) or injection of [11C]PIB (500 MBq). Catheters (Artflon and Venflon, Becton Dickinson, Swindon, UK) were inserted in the right radial artery and left cubital vein and arterial blood radioactivity was measured every half second for the duration of the PET scan by an automated blood sampling system (Allogg AB, Mariefred, Sweden), cross-calibrated with the tomograph, and then corrected for external delay and dispersion. For anatomical orientation, high-resolution T1-weighted MR images were obtained at 1.5 or 3 T (GE Sigma Systems).
Image registration and segmentation
The summed emission recordings of [15O]water and [11C]PIB were automatically co-registered to the individual MRI scans using a six parameter affine transformation. Individual MR Images were co-registered to a locally generated version of the common stereotactic space (ICBM, Montreal Neurologic Institute) (Mazziotta et al., 2001) using a combination of linear and non-linear registrations (Collins et al., 1994; Grabner et al., 2006). After the calculation of the final 16 mm non-linear PET-Talairach transformation grid, dynamic emission recordings were re-sampled into common coordinates. Regional BPND, R1, and CBF measures were obtained from parametric PET image maps using standard model based segmentation (Collins et al., 1994; Grabner et al., 2006). The regions analyzed were cerebral cortex excluding cerebellum (CORT), putamen (PU), caudate nucleus (CN) frontal (FL), occipital (OL), parietal (PL), and temporal (TL) lobes, as well as white matter (WM), and the cerebellar gray matter (CERB).
Quantification of [11C]PIB Retention
Reference region ratio measure (SUVR)
The [11C]PIB retention can be calculated by determining the accumulation relative to a reference tissue to obtain a ratio measure (SUVR). The ratio measure is the fraction of the region-of-interest integral of [11C]PIB accumulation at steady-state, assumed to have been established no later than this time after injection (ts = 40 min), extended to the end (te = 60 min), relative to the integral of the [11C]PIB accumulation observed in the same period in the reference region. In the reference region, we assume the accumulated tracer as function of time, mND(t) to represent non-specific binding after delivery of the tracer by homogeneous flow to all voxels of the reference region. In this context, we further assumed the interval from ts = 40 min to te = 60 min to be sufficient to establish steady-state or secular equilibrium in all regions, as the basis for the definition of the volume of distribution of the tracer, VT, as the sum of the volumes of distribution of non-displaceable tracer (VND) and an additional volume of distribution of displaceably bound tracer.
Washout allometric reference method (WARM)
In the case of negligible input from the circulation after the initial brief uptake, the tissue time-activity curves of the radio ligand are established by the radioactivity initially persed to the tissue and the subsequent washout from the brain regions of uptake. The WARM method specifically takes this condition into account and uses only the differences among washout rates from regions with different properties of binding, blood flow rates, and blood-brain barrier permeability. The condition means that the differential Equations (1) and (2)
are linked only while the tracer is dispersed from well-defined ca, i.e., during the brief uptake period until maximum peak (within 2–10 min timeframe) when washout is assumed to be negligible.
Equations (1) and (2) also form the base of the two compartment Simplified Reference Tissue Model (SRTM) (Lammertsma and Hume, 1996). The term K1 is the unidirectional clearance of the tracer ca by the tissue, KND1 is the clearance of the tracer ca by the reference region, m* and m*ND are the measured PET signal in the tissue (with displaceable binding) and reference, respectively. The term k2a defines the apparent measurable washout rate constant for the ROI. The term k2 is the unknown washout rate for non-specifically bound tracer of the same region of interest, and kND2 defines the measurable washout rate of non-specifically bound tracer in the reference tissue into the plasma. The uncoupling of the first and the second term on the right hand side of the equations means that elimination of the first K1 and KND1 terms yields the equations for the washout part of the signal.
The rate constants k2 and kND2 of washout of non-specifically bound tracer from tissue to plasma are linked by R1,
The distribution volume ratio(DVR) can be expressed as an allometric relationship between the logarithmic of the fraction of remaining tracer in a ROI relative to the deposited amount before washout, and the reference region where is the scaling exponent, found in this formulation by log–log linearization as,
where the ratio m*(t)/m*(0) is defined by the amount of tracer m*(t) remaining relative to the amount of tracer m*(0) initially deposited before the washout. Similarly for m*ND(t)/m*ND(0).
For direct calculation without linearization, BPND can be found using the operational equation
Intuitively, when log transformed, the fraction part (i.e., DVR) of this equation states that the nominator is the accumulated log-signal for the reference tissue relative to how much was present before washout, this difference is scaled by the initial tracer amount of the ROI. The denominator describes the accumulated log-signal for a ROI or voxel relative to how much was present before washout. This difference is scaled by the start amount of the reference region. Thus the fraction is corrected both for flow, i.e., initially deposited tracer, and the exponential behavior of the washout. For a more detailed derivation please refer to the Theory section.
Figure 1 illustrates the dynamics of the nominator, the denominator, and the BPND(T) terms of Equations (43) and (40) (in the Theory section) for simulated ROI and reference curves with 20% Gaussian noise added. As seen, the nominator and the denominator from Equation (40) (in panel 2 from the left) are stabilized by the integration in Equation (43) (panel 3 from the left). Panel 4 illustrates the convergence of the BPND(T) estimates toward the theoretical result.
Figure 1. Behavior of the parts of the warm method. The first leftmost panel shows two initial simulated mono-exponential washout curves, a reference curve (kND2 = 0.92) in red, and a binding curve (k2a = 0.98) in blue. The black curve shows the flow normalized version of the reference curve normalized to the ROI curve, with 20% Gaussian noise added. The second panel shows the denominators and nominators from Equation (40) as function of time, i.e., the fraction between these two ideally linear factors form the binding potential. The third panel shows the similar integrated nominator and denominator from Equation (43). The integration stabilizes the result. The fourth panel shows the PBND from Equation (40) in blue, and from Equation (43) in black. The dotted line is the theoretical true value for PBND. The PBND from both Equations (40) and (43) converge toward the true value with increasing tomographic duration.
Simplified reference tissue method
The simplified reference tissue method (SRTM) (Lammertsma and Hume, 1996) yields binding potential when a single tissue compartment model fits the data. SRTM solves differential equations similar to (25) and (26) (Equations (1) and (5) in Lammertsma and Hume (1996)). The method applies standard non-linear regression analysis to establish the relationship between the tissue concentrations of the region of interest and a reference region tissue such as typically the cerebellar gray matter in the case of [11C]PIB. The parameters estimated are R1, k2, and BPND. R1 accounts for differences of delivery to the regions of interest and reference.
Flow dependence of specific binding measure
With the regional flow ratio measure R1 derived from [15O]water or [11C]PIB analysis for each voxel with a signal n(t), we simulated flow-adjusted reference curves (nND(t)) for the corresponding voxel of the image, determining the tracer washout by the actual flow measured in the voxel. Here, the term nND(t) represents the simulated dynamic wash-out of the tracer that would have been recorded from an individual voxel in the absence of binding. The result is a simulated image of the dynamic passages of the tracer through every voxel as functions of time in the absence of any binding in any voxel. The simulation is based on the dynamic behavior in a fixed reference region (cerebellar gray matter) given by mND(t), which is the dynamic time-activity curve of the reference region. Figure 2 shows an example of the real and simulated non-displaceable (unbound) tracer time-activity curves for a small white matter region and a small putamen region, as well as the real time-activity curves for the ROIs. The obvious variation in the flow corrected reference curves illustrate that a single uncorrected reference curve may bias the result significantly. Further details can be found in the Theory section.
Figure 2. Time activity curves of PIB from a single subject in different regions of interest (ROI). In blue is shown the measured cerebellar reference curve. The purple and yellow curves are the putamen curves as well as the estimated reference curve. The red and green curves are the curves for a small white matter region and the corresponding estimated reference curve.
Quantification of [15O]Water Accumulation as CBF (KH2O1)
The use of positron emission tomography with [15O]water in subjects included in the present analysis for the purpose of correction of measures of CBF for the effect of CO2 tension in blood were reported previously (Rodell et al., 2012). We quantified the CBF as the unidirectional blood-brain [15O]water clearance (KH2O1) in units of ml hg−1 min−1 with the linearized two-compartment model (Blomquist, 1984) modification of Ohta et al. (1996) and the Lawson–Hanson non-negative least squares solution to general least squares functions (Lawson and Hanson, 1974).
Automatic Detection of Reference Region
In order to validate the cerebellar gray matter as a non-specifically binding region, we mapped the distribution of the washout index Θ image, which has unit of time, as defined for the Hypotime method (Møller et al., 2009). Distinctive regions of high Θ values (i.e., close to T = Tomography duration) are indicative of a potential reference region, as previously demonstrated for the tracer [11C]WAY-100635 ([11C]WAY) (Hirvonen et al., 2007; Møller et al., 2009) where the cerebellar white matter served as the reference region.
Wash-In and Wash-Out Phases
Due to the rapid removal of the tracer from the circulation, analysis of the dynamic [11C]PIB record revealed an initial high frequency signal of arterial origin in the first 2 min of recording, followed first by a maximum peak (within the timeframe 2–10 min) from which there was only by washout. For the WARM method, we split the signal into the three time frames of, first, arterial phase, second, peak uptake, and third, wash-out. We used the maximum peak within (the 2–10 min peak uptake phase) for estimation of the relative uptake coefficient RPIB1 and the calibrated surrogate CBF index directly from the [11C]PIB signal. We used the wash-out phase from maximum peak within 2–60 min of the dynamic record for estimation of the binding potential BPND.
Signal-To-Noise Ratio and Stability
In order to establish the temporal dependence and stability of the WARM method in relation to tomography duration, we analyzed the performance of the WARM and SRTM methods at different scan-time. After 10 min, the signal-to-noise ratio increased for [11C]PIB because of the decay of the radioactivity and washout of the tracer. After 60 min, the standard deviation reached almost half the signal, as illustrated in Figure 3 (left panel). Figure 3 (right panel) shows the stability of the SRTM and WARM method results with respect to time in the tomograph. The WARM method results converge after 60 min, while the SRTM results retain some dependence on time. Considering the scale used for BPND, the results of the two methods are in relative good agreement. Considering the SNR and stability together, we regard a maximum time of 60 min to be an acceptable compromise. For this reason, we confined the calculations of binding potential with the WARM method to the 2–60 min time window. The SRTM calculations were made with a maximum duration of 60 min.
Figure 3. The left panel shows the relative noise (SNR) of a typical PIB scan as a function of tomography time duration for a region covering the cerebral cortex without cerebellum. After 60 min the noise level is more than half the signal, probably due to the loss of the signal due to decay of [11C] and washout of the tracer. the right panel shows the stability of the BPND measure for the WARM method and SRTM, for each method is shown the result of applying the analysis to images with different end frame times. The curves represent the parametric BPND values for a typical AD patient for a single voxel in frontal lobe and a striatal region mean value. As seen the methods are in relatively good agreement, but the WARM method converge with more and more data considered.
In order to test the methods' power to distinguish the AD from the HC subject groups by quantification of the amount of retention, we extracted regional parametric values for the PIB retention from the parametric images with the three different methods WARM, SRTM and SUVR. For the WARM, SRTM methods, we report mean absolute BPND values and for the SUVR method we extracted DVR values.
The mean absolute BPND values from the WARM and SRTM methods were in very good agreement for cerebral cortex values for the AD patients, while for the HC subjects the WARM method values are lower for most of cortex. The WARM method therefore yielded more significant differences between the AD and HC groups than the SRTM method for cortex (CORT), putamen (PU), frontal lobe (FL), occipital lobe (OL), temporal lobe (TL), and parietal lobe (PL). Only caudate nucleus (CN), white matter (WM) and cerebellum (CERB) reference values were less significant than with the SRTM method. Compared with the SUVR method, both the WARM and SRTM methods yielded greater between-group differences of retention. Table 1 lists the significance levels, and Table 2 lists the mean parametric estimates for each region and group.
Table 1. Regional statistical significance (student t-test) between AD vs. HC subject values for both [11C]PIB retention and flow estimates (*p < 0.05) (**p < 0.001).
Table 2. BPND and DVR values measured for AD vs. HC subjects using three different methods (*indicates statistical significance between subject groups p < 0.05).
Closer examination of the WM binding revealed a considerable amount of white matter retention with the WARM method, although the flow normalization accounted for some of this. The SRTM method yielded negative binding values for these WM areas (not visible within the specified range). We judge the negative values to be due to difficulties of regression to the WM dynamics.
To illustrate the regional distribution of the binding, Figure 4 shows the group mean binding images for the three different methods. The WARM and SRTM method results are reported as BPND values, scaled similarly, while the SUVR values are reported as DVR values in the 1–3 range.
Figure 4. Images of [11C]PIB binding in brain. Upper three rows show mean binding images of sections of the brains of a group of AD patients (N = 6), and lower three rows show similar sections for age-matched healthy controls (N = 8). For each group is shown the BPND values derived from the WARM and SRTM methods. For the SUVR method, the mean DVR measures are illustrated in the range 1–3.
Cerebral Blood Flow
Surprisingly, the absolute CBF values derived from [15O]water (Table 3) most significantly distinguished the AD from the HC groups (Table 1), especially after correction for the within-group variation of arterial CO2 tensions. The low standard deviations compared to the retention values probably accounted for the more significant differences among the CBF estimates. The CBF values of the reference region (cerebellar gray matter) differed among the groups. This difference disappeared when we considered the relative flow indices, i.e., the R1 values (Table 4). Thus, normalization may remove important differences and should be used with caution (Mayr, 1982; Borghammer et al., 2009a,b). The R1 estimates obtained with [15O]water and [11C]PIB were similar (Table 4), and for both methods the temporal and parietal lobes were the only regions with significantly different values when we compared AD and HC subjects. Hence, it is feasible to use R1 estimates from the [11C]PIB sequence to correct for intra-subject flow normalization, also when inter-individual flow differences may be of clinical value as a marker of AD pathology.
Table 3. Regional absolute CBF measures before and after correction for PCO2 (Rodell et al., 2012).
We used a site-specific scale factor to obtain a simulated absolute flow estimates from the PIB peak arrival characteristics, but this factor was derived only from healthy subjects. We did not find the same statistical differences for the simulated CBF values (Table 5) as for the standard CBF estimates. As illustrated for select regions in Figure 5, this failure may be due to the high variability observed in the simulated CBF values, as generally there was good correlation for the group means seen in Figure 6, except for white matter and cerebellar gray matter. Here, the simulated CBF values did not yield the same CBF value decrease for the AD subjects.
Table 5. Regional surrogate CBF measures before and after correction for PCO2 (Rodell et al., 2012).
Figure 5. Distribution and mean relationships between the [11C]PIB derived surrogate CBF values (lower panels) and the standard [15O]water derived CBF measures (upper panels) for selected regions in healthy controls (left panels) and AD patients (right panels).
Figure 6. Left panel shows relationship between the [11C]PIB estimated surrogate CBF measures and the standard [15O]water-measured CBF values for averages of healthy controls and AD patients (ml hg−1 min−1). Right panel shows the linear relationship between the different regions averaged for the HC and AD groups.
A clear advantage of the Hypotime method was the production of washout indices in the form of Θ images with unit of time as an intermediate result (Møller et al., 2009). High Θ values (Θ > 1660 s relative to T = 1890 s) are indicative of the region's potential for service as reference. This was previously demonstrated in studies of the WAY tracer, where the reference region was selected as the cerebellar white matter. However, for the present tracer, the cerebellar gray matter consistently stood out as a reference region candidate. Figure 7 shows the average Θ washout index image for all subjects, for comparison with an insert of the reference region previously identified for WAY (Møller et al., 2009).
Figure 7. The figure shows the average washout index Θ image calculated over the time interval T = 1890 s for all subjects. Mean values >1660 s, superimposed on a non-linear registered average MRI brain, the cyan insert shows the similar reference region for the way compound as found in Møller et al. (2009).
For [11C]PIB, the mono-exponential washout kinetics with rapid disappearance of tracer from the circulation complicates the qualitative assessment of binding with existing methods. Here we demonstrate the use of the novel flow-independent WARM to calculate the binding potential of [11C]PIB. The calculation is based on the clearance of the tracer from the brain, relative to the initial deposits in brain tissue, compared to the same measure in a non-binding reference region with the same initial deposit and washout characteristics of the cerebellar gray matter. We measure the binding either regionally by allometric linearization using Equation (41), or by direct parametric mapping using Equation (43). Compared to the standard SRTM and SUVR methods, WARM yielded binding measures with greater statistical discrimination between groups of AD and HC subjects for cortical gray-matter (excluding cerebellum), putamen, frontal, occipital, temporal, and parietal lobes, and lower differences for caudate nucleus, white matter, and especially cerebellum, as expected. This has relevance for the enhanced use of quantified [11C]PIB retention for clinical discrimination of AD, regardless of the relationship between amyloid-β deposits and the disease.
The WARM method revealed a converging stability of the BPND measure as a function of tomography duration with acceptable stability at 60 min, at a time when the SNR is still favorable. The method also takes advantage of direct calculation without regression. The direct calculation is important to its use because it serves to reduce tomography duration without compromise of the quantification, without the loss of the uptake signal, which we have shown to be useful for flow estimation.
Based on the Θ washout index with unit of time, defined in Møller et al. (2009), we found that the cerebellar gray matter consistently is the best choice of reference region with no specific binding of [11C]PIB. The striking difference between the cerebellar reference region for tracer WAY (Hirvonen et al., 2007; Møller et al., 2009) and [11C]PIB, testifies to the value of directly confirming the presence of a true reference area from the PET sequence.
We also showed the feasibility of using the initial K1 signal as a surrogate measure of absolute CBF, directly from the [11C]PIB image. However, compared to absolute CBF measures with [15O]water, the CBF index derived from [11C]PIB had greater variability. As [11C]PIB is comparatively lipophilic, there are probable permeability differences from [15O]water images (Gjedde et al., 2013), which may account for some of the variability. The absolute CBF estimates differed significantly for all regions, especially after correction for arterial PCO2 differences (Rodell et al., 2012), notably also for the cerebellar reference region. When we calculated the R1 values relative to the cerebellar reference, much of the difference was eliminated, and only parietal and temporal lobes had significant decline in the AD group. This observation shows that the normalization to relative flow measures may mask important disease specific information in AD, although it does enable the correction for the intra-individual flow bias (Borghammer et al., 2009a). The observation is of concern also to the MRI estimation of CBF indices which rarely are reported as calibrated measures. Clinically, the absolute CBF value should be regarded as an important imaging parameter for the diagnosis of AD, alongside the [11C]PIB retention, the fluorodeoxyglucose (FDG) uptake, and the accumulation of other markers.
Non-specific binding in white matter (WM) areas is found both when assessed in vivo by PET (Fodero-Tavoletti et al., 2009) and by postmortem autoradiography (Svedberg et al., 2009). This implies some non-specific binding to WM, also when tissue slices are superfused with the tracer, rather than accumulating the tracer from the circulation after transport across the blood-brain barrier. However, we find that the WM binding often observed with [11C]PIB is due both to slow washout kinetics and to non-specific retention that is not explicable by R1 differences.
In conclusion, we demonstrate that the WARM method yields a stable measure of PIB's binding potential with relative simplicity and reasonable tomography duration, employing only integration for noise reduction with no need for model regression. The washout index Θ is readily used to assess the suitability of a potential reference region. The method accounts for relative flow differences in the tissue and yields a calibrated measure of the absolute CBF, obtained directly from the [11C]PIB signal. Comparable with the SRTM and SUVR methods, the WARM method provides better discrimination between AD subjects and healthy controls. Taken together these characteristics (Table 6) merit further investigation of WARM for clinical use with [11C]PIB, and the method may be equally applicable to other washout tracers.
This section describes theoretical considerations regarding the three different approaches to the determination of the binding potential of [11C]PIB (BPND) or DVR.
Reference Region Ratio Measure (SUVR)
The [11C]PIB retention can be calculated by determining the accumulation relative to a reference tissue to obtain a ratio measure (SUVR). The ratio measure is the fraction of the region-of-interest integral of [11C]PIB accumulation at steady-state, assumed to have been established no later than this time after injection (ts = 40 min), extended to the end (te = 60 min). The volume of distribution of the tracer, VT and the volume of distribution of non-displaceable tracer (VND) define the binding potential BPND.
The signal measured by PET is the accumulated mass of the tracer (VT) relative to the integral of ca(t), the radioactivity of the tracer in the arterial circulation i.e.,
The variable of major interest to [11C]PIB binding to Aβ formations is the quantity of specific (i.e., displaceable) binding. This binding is expressed as the binding potential (BPND), the displaceable tracer relative to the amount of non-displaceable tracer in the tissue. In any given region of interest (ROI), an area of reference may be referred to as a region where equal amounts of non-displaceable tracer enter and leave the tissue. Such a reference fulfills the two requirements that,
(1) tracer enters the tissue in proportion to blood flow,
(2) unbound tracer (i.e., tracer dissolved in VND) clears the tissue at the rate determined by the blood flow, as defined for a tracer subject to flow-limited exchange across the blood-brain barrier.
For the reference region, VT ≡ VND, and Equation (10) defines the steady-state volume of distribution,
As the total amount of radioactivity in the arterial blood are the same in the two regions, i.e., the ROI (Equation 10), and the reference region (Equation 11), the binding potential is obtained from the Equation (9).
Equation (12) is valid when requirements (1) and (2) are fulfilled. However, for a tracer that disappears rapidly from the blood stream, (1) is invalid after the initial wave of tracer has passed the tissue. When the concentrations ca(t) in the circulation primarily depend on the wash-out and hence mainly on the CBF. Requirement (2) depends on the regional CBF and is valid for tracers subject to flow-limited exchange with brain tissue.
It is improbable that we would find a reference region for all tissues and tracers that uphold these requirements, unless additional tomography is completed with a similar tracer of no specific binding and hence no displacement. However, when we identify a single homogeneous reference region for [11C]PIB, we can mimic the tracer's behavior as reference for any given ROI.
Initial [11C]PIB Distribution as Surrogate CBF Measure
The two factors that affect the accumulation of [11C]PIB after the initial distribution are the magnitude of blood flow, which mediates washout, and the degree of binding, which limits the rate of washout. In order to distinguish these factors and hence to assess the binding of [11C]PIB in any given region or voxel, we determine the relative flow ratios R1 as the magnitude of CBF in a region or voxel, relative to CBF in the reference region, defined as R1 = CBF/CBFND.
Cerebral blood flow can be assessed in two different ways, either directly by means of PET with [15O]water (CBF), or a surrogate measure (CBFPIB) approximated indirectly from the [11C]PIB signal at the peak of distribution. This approximation is made by first excluding the instantaneous arterial distribution in the first 2 min after injection. In the subsequent 8 min, the initial distribution of the [11C]PIB signal depends largely on the wash-in. Until maximum peak values are reached at Δtp, the wash-in can be expressed as the unidirectional clearance (K1) from the blood into the brain when no tracer has left the brain yet.
We recently determined the permeability-surface product (PS) for [11C]PIB (Gjedde et al., 2013). There the K1 measure from the initial distribution of the [11C]PIB signal is related to the CBF by the Renkin–Crone formula (Crone, 1963; Renkin, 1964),
where by definition the extraction fraction E from blood into tissue is,
From Equation (13), it follows that the wash-in ratio
relative to a reference region is linked to the similar ratio for blood flow R1 = CBF/CBFND by the relative extraction fraction relative to the extraction fraction of the reference region END
This relation indicates that the RPIB1 approximates R1 for extraction fractions similar to the reference region, which is the case for compounds with sufficiently high PS products and for regions where the CBF is close to the CBFND of the reference region. With these limitations, and for want of known CBF values, RPIB1 may be estimated from the maximum signal intensity value m(Δtp), such that RPIB1 = m(Δtp)/mND(Δtrefp), relative to the maximum signal intensity of the reference region at time Δtrefp.
A tentative absolute measure for K1 can be estimated by additionally accounting for differences in the delay of the maximum peak and by normalizing for weight and dose, (m(Δtp)/Δtp) * Δtrefp * weight/dose provides an delay normalized estimate of the magnitude of K1. The estimate of K1 is scaled by a site- and tomograph-specific constant Ksite for cerebral cortex values in order to derive tentative absolute flow estimates CBFPIB from the [11C]PIB images. This derivation requires calibration to a normal material of estimates of cortex K1 and CBF values, with tight adherence to the protocol used in the calibration.
Flow Dependence of Specific Binding Measure
With the regional flow ratio measure R1 derived from [15O]water or [11C]PIB analysis for each voxel with a signal n(t), we simulated flow-adjusted reference curves (nND(t)) for the corresponding voxel of the image, determining the tracer washout by the actual flow measured in the voxel. Here, the term nND(t) represents the simulated dynamic wash-out of the tracer that would have been recorded from an individual voxel in the absence of binding. The result is a simulated image of the dynamic passages of the tracer through every voxel as functions of time in the absence of any binding in any voxel. The simulation is based on the dynamic behavior in a fixed reference region (cerebellar gray matter) given by mND(t), which is the dynamic time-activity curve of the reference region.
In order to determine the simulated reference curve nND(t) from the measured reference curve mND(t), let a non-binding voxel or region have a flow CBFn, and let it refer to a measured region with flow CBFND. Each region has washout rates given by k2 = CBFn/VND and kND2 = CBFND/VND, respectively. Then R1 relates to the variables as indicated by the following equations,
where each voxel is the site of mono-exponential washout given by rate constants of the magnitudes,
where R1 relates the k2 terms as,
such that Equations (19) and (20) yield nND(t) as,
where nND(t) is the simulated flow-adjusted reference curve corresponding to the measured curve n(t) when equal amounts of tracer enter, i.e., n(0) = nND(0). When expressed for every voxel, the procedure yields a new dynamic reference image nND(t) that enables the comparisons among the flow-involved SUVR method and the flow-devolved WARM and SRTM methods, as the flow-devolved version of Equation (12) can now be expressed as,
Figure 2 shows an example of the real and simulated non-displaceable (unbound) tracer time-activity curves for a small white matter region and a small putamen region, as well as the real time-activity curves for the ROIs. The variation in the simulated flow corrected reference curves, illustrate that a single uncorrected reference curve may bias the result significantly and may only be valid for regions where the flow is equal to in this case the cerebellum.
Washout Allometric Reference Method (WARM)
In the case of negligible input from the circulation after the initial brief uptake, the differential equations (25) and (26)
are linked only while the tracer is dispersed from well-defined ca, i.e., during the brief uptake period until maximum peak (within 2–10 min timeframe) when washout is assumed to be negligible.
The term K1 is the unidirectional clearance of the tracer ca by the tissue, KND1 is the clearance of the tracer ca by the reference region, m* and m*ND are the measured PET signal in the tissue (with displaceable binding) and reference, respectively. The term k2a defines the apparent measurable washout rate constant for the ROI. The term k2 is the unknown washout rate for non-specifically bound tracer of the same region of interest, and kND2 defines the measurable washout rate of non-specifically bound tracer in the reference tissue into the plasma. The uncoupling of the first and the second term on the right hand side of the equations means that elimination of the first K1 and KND1 terms yields the equations.
This defines the first-order decay for 1-compartment first order kinetics, assuming that association and dissociation is sufficiently rapid in the tissue compartments. The Equations (27) and (28) each predict a mono-exponential washout from the time t0 where the [11C]PIB signal depends only on the washout rate. The total volume of distribution in the binding region then is simplified to (Lammertsma and Hume, 1996; Møller et al., 2009)
As the differential Equations (25) and (26) are linked in the initial phase (described by K1), the ratio R1 = K1/KND1 accounts for the difference of delivery to the regions of interest and reference. For the non-specifically bound tracer, we assume that the volume of distribution is the same in all regions, i.e.,
Consequently the rate constants k2 and kND2 of washout of non-specifically bound tracer from tissue to plasma are similarly linked by R1,
by combining which Equations (30) and (32) BPND is found as
The late uncoupled mono-exponential differential equations (27) and (28) can be rearranged to
Integrating both sides the gives
where m*(0) = m*(t0) and m*ND(0) = m*ND(t0) are the initial value at start of the washout t0 = 2 min. By dividing Equation (37) by (36)
The right hand side of this equation can be expressed in terms of the binding potential by use of Equation (33).
where (1 + BPND) is the distribution volume ratio(DVR). Alternatively written
where (1 + BPND) is the distribution volume ratio(DVR).
Equation (39) can also be expressed as an allometric relationship between the logarithmic of the fraction of remaining f tracer in a ROI and the reference region where is the scaling exponent, found in this formulation by log–log linearization as,
For direct calculation without linearization, Equation (39) was expressed in terms of integrations and differences of the logarithms, assuming constant BPND and R1, which greatly reduces noise in BPND estimates.
That we further simplified with Equation (17) to the operational equation,
When log transformed, the fraction (i.e., DVR) part of this equations states that the nominator is the accumulated log-signal for the reference tissue relative to how much was present before washout, this difference is scaled by the initial tracer amount of the ROI. The denominator describes the accumulated log-signal for a ROI or voxel relative to how much was present before washout. This difference is scaled by the start amount of the reference region. Thus the fraction is corrected both for flow, i.e., initially deposited tracer, and the exponential behavior of the washout. Figure 1 illustrates the behavior of the nominator, denominator and BPND(T) of Equations (43) and (40) for simulated ROI and reference curves with 20% added Gaussian noise. As seen the nominator and denominator from Equation (40) (in panel 2 from the left) are stabilized by the integration in Equation (43) (panel 3 from the left). Panel 4 illustrates how the BPND(T) values converge towards the theoretical result.
Simplified Reference Tissue Method
The simplified reference tissue method (SRTM) (Lammertsma and Hume, 1996) yields binding potential when a single tissue compartment model fits the data. SRTM solves differential equations similar to (25) and (26) (equations (1) and (5) in the paper of Lammertsma and Hume (1996)). The method assumes that these differential equations are coupled by tracer distribution in the circulation throughout the entire duration of the tomography.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Borghammer, P., Aanerud, J., and Gjedde, A. (2009a). Data-driven intensity normalization of PET group comparison studies is superior to global mean normalization. Neuroimage 46, 981–988. doi: 10.1016/j.neuroimage.2009.03.021
Borghammer, P., Cumming, P., Aanerud, J., Förster, S., and Gjedde, A. (2009b). Subcortical elevation of metabolism in Parkinson's disease–a critical reappraisal in the context of global mean normalization. Neuroimage 47, 1514–1521. doi: 10.1016/j.neuroimage.2009.05.040
Cohen, A. D., Rabinovici, G. D., Mathis, C. A., Jagust, W. J., Klunk, W. E., and Ikonomovic, M. D. (2012). Using Pittsburgh compound B for in vivo PET imaging of fibrillar amyloid-beta. Adv. Pharmacol. 64, 27–81. doi: 10.1016/B978-0-12-394816-8.00002-7
Collins, D. L., Neelin, P., Peters, T. M., and Evans, A. C. (1994). Automatic 3D inter-subject registration of MR volumetric data in standardized talairach space. J. Comput. Assist. Tomogr. 18, 192–205. doi: 10.1097/00004728-199403000-00005
Ducharme, S., Guiot, M. C., Nikelski, J., and Chertkow, H. (2013). Does a positive Pittsburgh compound B scan in a patient with dementia equal Alzheimer disease? JAMA Neurol. 70, 912–914. doi: 10.1001/jamaneurol.2013.420
Fodero-Tavoletti, M. T., Rowe, C. C., McLean, C. A., Leone, L., Li, Q. X., Masters, C. L., et al. (2009). Characterization of PiB binding to white matter in Alzheimer disease and other dementias. J. Nucl. Med. 50, 198–204. Epub 2009 Jan 21. doi: 10.2967/jnumed.108.057984
Grabner, G., Janke, A. L., Budge, M. M., Smith, D., Pruessner, J., and Collins, D. L. (2006). Symetric atlasing and model based segmentation: an application to the Hippocampus in older adults. Med. Image Comput. Comput. Assist. Interv. 9(Pt 2), 58–66. doi: 10.1007/11866763_8
Hirvonen, J., Kajander, J., Allonen, T., Oikonen, V., Någren, K., and Hietala, J. (2007). Measurement of serotonin 5-HT1A receptor binding using positron emission tomography and [carbonyl-(11)C]WAY-100635-considerations on the validity of cerebellum as a reference region. J. Cereb. Blood Flow Metab. 27, 185–195. doi: 10.1038/sj.jcbfm.9600326
Ikonomovic, M. D., Abrahamson, E. E., Price, J. C., Hamilton, R. L., Mathis, C. A., Paljug, W. R., et al. (2012). Early AD pathology in a [C-11]PiB-negative case: a PiB-amyloid imaging, biochemical, and immunohistochemical study. Acta Neuropathol. 123, 433–447. doi: 10.1007/s00401-012-0943-2
Joshi, A. D., Pontecorvo, M. J., Clark, C. M., Carpenter, A. P., Jennings, D. L., Sadowsky, C. H., et al. (2012). Performance characteristics of amyloid PET with florbetapir F 18 in patients with Alzheimer's disease and cognitively normal subjects. J. Nucl. Med. 53, 378–384. doi: 10.2967/jnumed.111.090340
Kanno, I., and Uemura, K. (1975). Some experimental errors in calculating regional cerebral blood flow from the intracarotid 133Xe clearance curve. A quantitative evaluation employing a digital model. Stroke 6, 370–375. doi: 10.1161/01.STR.6.4.370
Kepe, V., Moghbel, M. C., Långström, B., Zaidi, H., Vinters, H. V., Huang, S. C., et al. (2013). Amyloid-β positron emission tomography imaging probes: a critical review. J. Alzheimers Dis. 36, 613–631. doi: 10.3233/JAD-130485
Klunk, W. E., Engler, H., Nordberg, A., Wang, Y., Blomqvist, G., Holt, D. P., et al. (2004). Imaging brain amyloid in Alzheimer's disease with Pittsburgh compound-B. Ann. Neurol. 55, 306–319. doi: 10.1002/ana.20009
Klunk, W. E., Price, J. C., Mathis, C. A., Tsopelas, N. D., Lopresti, B. J., Ziolko, S. K., et al. (2007). Amyloid deposition begins in the striatum of presenilin-1 mutation carriers from two unrelated pedigrees. J. Neurosci. 27, 6174–6184. doi: 10.1523/JNEUROSCI.0730-07.2007
Knight, W. D., Okello, A. A., Ryan, N. S., Turkheimer, F. E., Rodríguez Martinez de Llano, S., Edison, P., et al. (2011). Carbon-11-Pittsburgh compound B positron emission tomography imaging of amyloid deposition in presenilin 1 mutation carriers. Brain 134(Pt 1), 293–300. Epub 2010 Nov 16. doi: 10.1093/brain/awq310
Landau, S. M., Breault, C., Joshi, A. D., Pontecorvo, M., Mathis, C. A., Jagust, W. J., et al. (2013). Alzheimer's disease neuroimaging initiative. Amyloid imaging with Pittsburgh compound B and florbetapir: comparing radiotracers and quantification methods. J. Nucl. Med. 54, 70–77. doi: 10.2967/jnumed.112.109009
Lehmann, M., Ghosh, P. M., Madison, C., Laforce, R. Jr., Corbetta-Rastelli, C., Weiner, M. W., et al. (2013). Diverging patterns of amyloid deposition and hypometabolism in clinical variants of probable Alzheimer's disease. Brain 136(Pt 3), 844–858. doi: 10.1093/brain/aws327.
Lolk, A., Nielsen, H., Andersen, K., Andersen, J., and Kragh-Srensen, P. (2000). CAMCOG as a screening instrument for dementia: the Odense study. Cambridge cognitive examination. Acta Psychiatr. Scand. 102, 331–335. doi: 10.1034/j.1600-0447.2000.102005331.x
Møller, M., Rodell, A., and Gjedde, A. (2009). Parametric mapping of 5HT1A receptor sites in the human brain with the Hypotime method: theory and normal values. J. Nucl. Med. 50, 1229–1236. doi: 10.2967/jnumed.108.053322
Mazziotta, J. A., Toga, A. W., Evans, A. C., Fox, P. T., Lancaster, J., Zilles, K., et al. (2001). A probabilistic atlas and reference system for the human brain: international consortium for brain mapping (ICBM). Philos. Trans. R. Soc. Lond. B Biol. Sci. 356, 1293–1322. doi: 10.1098/rstb.2001.0915
Nilsson, P., Loganathan, K., Sekiguchi, M., Matsuba, Y., Hui, K., Tsubuki, S., et al. (2013). Aβ Secretion and plaque formation depend on autophagy. Cell Rep. 5, 61–69. doi: 10.1016/j.celrep.2013.08.042
Ohta, S., Meyer, E., Fujita, H., Reutens, D. C., Evans, A., and Gjedde, A. (1996). Cerebral [15O]water clearance in humansdetermined by PET: I. Theory and normal values. J. Cereb. Blood Flow Metab. 16, 765–780. doi: 10.1097/00004647-199609000-00002
Rodell, A. B., Aanerud, J., Braendgaard, H., and Gjedde, A. (2012). Low residual CBF variability in Alzheimer's disease after correction for CO(2) effect. Front. Neuroenergetics. 4:8. doi: 10.3389/fnene.2012.00008
Rosen, R. F., Ciliax, B. J., Wingo, T. S., Gearing, M., Dooyema, J., Lah, J. J., et al. (2010). Deficient high-affinity binding of Pittsburgh compound B in a case of Alzheimer's disease. Acta Neuropathol. 119, 221–233. doi: 10.1007/s00401-009-0583-3
Svedberg, M. M., Hall, H., Hellström-Lindahl, E., Estrada, S., Guan, Z., Nordberg, A., et al. (2009). [(11)C]PIB-amyloid binding and levels of Abeta40 and Abeta42 in postmortem brain tissue from Alzheimer patients. Neurochem. Int. 54, 347–357. doi: 10.1016/j.neuint.2008.12.016
Varrone, A., Sjöholm, N., Eriksson, L., Gulyás, B., Halldin, C., and Farde, L. (2009). Advancement in PET quantification using 3D-OP-OSEM point spread function reconstruction with the HRRT. Eur. J. Nucl. Med. Mol. Imag. 36, 1639–1650. doi: 10.1007/s00259-009-1156-3
Villemagne, V. L., Klunk, W. E., Mathis, C. A., Rowe, C. C., Brooks, D. J., Hyman, B. T., et al. (2012). Aβ imaging: feasible, pertinent, and vital to progress in Alzheimer's disease. Eur. J. Nucl. Med. Mol. Imag. 39, 209–219. doi: 10.1007/s00259-011-2045-0
Keywords: Alzheimer's disease, CBF, Aβ, PIB, flow normalization, parametric imaging
Citation: Rodell A, Aanerud J, Braendgaard H and Gjedde A (2013) Washout allometric reference method (WARM) for parametric analysis of [11C]PIB in human brains. Front. Aging Neurosci. 5:45. doi: 10.3389/fnagi.2013.00045
Received: 05 August 2013; Paper pending published: 10 September 2013;
Accepted: 07 November 2013; Published online: 27 November 2013.
Edited by:Hari S. Sharma, Uppsala University, Sweden
Reviewed by:J. Arturo García-Horsman, University of Helsinki, Finland
Aruna Sharma, Uppsala University Hospital, Sweden
Copyright © 2013 Rodell, Aanerud, Braendgaard and Gjedde. 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) or licensor 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: Anders Rodell, Department of Nuclear Medicine and PET Centre, Aarhus University Hospital, Norrebrogade 44, Aarhus 8000, Denmark e-mail: email@example.com