Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2001263

Formats

Article sections

Authors

Related links

Neuroimage. Author manuscript; available in PMC 2008 June 1.

Published in final edited form as:

Published online 2007 March 16. doi: 10.1016/j.neuroimage.2007.03.004

PMCID: PMC2001263

NIHMSID: NIHMS24112

Yun Zhou,^{1} Susan M. Resnick,^{2} Weiguo Ye,^{1} Hong Fan,^{1} Daniel P. Holt,^{1} William E. Klunk,^{3} Chester A. Mathis,^{3} Robert Dannals,^{1} and Dean F. Wong^{1}

Author for correspondence: Yun Zhou, Ph.D., The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, 601 N. Caroline St., JHOC room 3245, Baltimore, MD 21287-0807, Phone: (410) 955-9798, Fax: (410) 955-0696, Email: ude.imhj@uohznuy

The publisher's final edited version of this article is available at Neuroimage

See other articles in PMC that cite the published article.

Reference tissue model (RTM) is a compartmental modeling approach that uses reference tissue time activity curve (TAC) as input for quantification of ligand-receptor dynamic PET without blood sampling. There are limitations in applying the RTM for kinetic analysis of PET studies using [^{11}C]Pittsburgh compound B ([^{11}C]PIB). For region of interest (ROI) based kinetic modeling, the low specific binding of [^{11}C]PIB in a target ROI can result in a high linear relationship between the output and input. This condition may result in amplification of errors in estimates using RTM. For pixel-wise quantification, due to the high noise level of pixel kinetics, the parametric images generated by RTM with conventional linear or nonlinear regression may be too noisy for use in clinical studies.

We applied RTM with parameter coupling and a simultaneous fitting method as a spatial constraint for ROI kinetic analysis. Three RTMs with parameter coupling were derived from a classical compartment model with plasma input: a RTM of 4 parameters (R_{1}, k′_{2R}, k_{4}, BP) (RTM4P); a RTM of 5 parameters (R_{1}, k_{2R}, NS, k_{6}, BP) (RTM5P); and a simplified RTM (SRTM) of 3 parameters (R_{1}, k′_{2R}, BP) (RTM3P). The parameter sets [k′_{2R}, k_{4}], [k_{2R}, NS, k_{6}], and k′_{2R} are coupled among ROIs for RTM4P, RTM5P, and RTM3P, respectively. A linear regression with spatial constraint (LRSC) algorithm was applied to the SRTM for parametric imaging. Logan plots were used to estimate the distribution volume ratio (DVR) (= 1 + BP (binding potential)) in ROI and pixel levels. Ninety-minute [^{11}C]PIB dynamic PET was performed in 28 controls and 6 individuals with mild cognitive impairment (MCI) on a GE Advance scanner. ROIs of cerebellum (reference tissue) and 15 other regions were defined on coregistered MRI’s.

The coefficients of variation of DVR estimates from RTM3P obtained by the simultaneous fitting method were lower by 77 - 89% (in striatum, frontal, occipital, parietal, and cingulate cortex) as compared to that by conventional single ROI TAC fitting method. There were no significant differences in both TAC fitting and DVR estimates between the RTM3P and the RTM4P or RTM5P. The DVR in striatum, lateral temporal, frontal and cingulate cortex for MCI group was 25% to 38% higher compared to the control group (p ≤ 0.05), even in this group of individuals with generally low PIB retention. The DVR images generated by the SRTM with LRSC algorithm had high linear correlations with those from the Logan plot (R^{2} = 0.99). In conclusion, the RTM3P with simultaneous fitting method is shown to be a robust compartmental modeling approach that may be useful in [^{11}C]PIB PET studies to detect early markers of Alzheimer’s disease where specific ROIs have been hypothesized. In addition, the SRTM with LRSC algorithm may be useful in generating R_{1} and DVR images for pixel-wise quantification of [^{11}C]PIB dynamic PET.

Positron emission tomography (PET) with [^{11}C]Pittsburgh compound B ([^{11}C]PIB) has been used for in vivo imaging of amyloid-β (Aβ) in Alzheimer’s disease (AD), mild cognitive impairment (MCI), and aging in humans (Mathis et al., 2004; Buckner et al., 2005; Klunk et al., 2004, 2005). The full chemical name for [^{11}C]PIB is [*N*-methyl-^{11}C]2-(4′-methylaminophenyl)-6-hydroxybenzothiazole (or [^{11}C]6-OH-BTA-1) that has binding affinity K_{D} = 1.4 nM for homogenates of post-mortem AD frontal cortex and K_{D} = 4.7 nM for synthetic Aβ (Mathis et al., 2003). Human studies using [^{11}C]PIB PET have indicated greater retention of [^{11}C]PIB in the brains of AD patients and subjects with MCI as compared to the healthy controls (Price et al., 2005; Lopresti et al., 2005; Mintun et al., 2006), as well as an inverse association between PIB retention and CSF Aβ (Fagan et al., 2006). [^{11}C]PIB is the most widely used PET imaging agent (the other is [^{18}F]FDDNP) in research studies aimed at improving early detection of AD, monitoring progression of Aβ deposition in the brain, and evaluating anti-amyloid and other therapies to stop progression of AD (Shoghi-Jadid et al., 2002; Mathis et al., 2004; Mintun 2005; Nichols et al., 2006; Nordberg 2004; Small et al., 2006; Wu et al., 2005).

The standard compartmental model with plasma input as well as model-independent spectral analysis and graphical analysis (Logan plot) with plasma input were used for [^{11}C]PIB kinetic analysis (Price at al., 2005; Verhoeff et al., 2004). Consistent results from the two previous studies demonstrated that 1) a 2-tissue compartmental model provided better cure fitting than 1-tissue compartmental model; 2) the Logan plot with plasma input is a robust approach to estimate [^{11}C]PIB distribution volume (DV) as compared to the 2-tissue compartmental model; and 3) there was no significant difference in the DV estimates for reference tissue (cerebellum) between controls and AD patients. The plasma input is usually obtained by arterial blood sampling during the PET study period, and tracer metabolism in plasma is corrected using the HPLC technique. This procedure is laborious and is associated with experimental errors and risks to subjects, particularly in the context of frequent longitudinal follow-up. Thus, the ability to conduct accurate studies without arterial sampling will increase the feasibility of [^{11}C]PIB for clinical practice and increase recruitment and retention of participants within the context of large, longitudinal studies. To quantify [^{11}C]PIB dynamic PET without arterial blood sampling, compartmental model with the plasma input derived from dynamic image data, graphical analysis (Logan plot) and a simplified reference tissue model (SRTM) with reference tissue input, and standardized uptake value ratio (SUVR) or target to reference tissue concentration ratio were evaluated (Edison et al., 2007; Fagan et al., 2006; Kemppainen et al., 2006; Lopresti et al., 2005; Price et al., 2005).

Reference tissue model is a compartmental modeling approach that uses the reference tissue time activity curve (TAC) as input (Cunningham et al., 1991; Gunn et al., 2000; Lammertsma et al., 1996; Lammertsma and Hume 1996; Morris et al., 2005; Watabe et al., 2000). In contrast to graphical analysis, the parameters of reference tissue model are estimated by fitting the model to the full time course of tissue TAC measured by PET. Analogous to the classical compartmental model with plasma input, the reference tissue model can be used to predict and simulate tissue tracer kinetics with given model parameters and reference tissue input. Compared to graphical analysis, the reference tissue model is commonly used to extract more physiological information from measured tracer kinetics, such as the relative tracer transport rate constant from vascular space to tissue. In addition, reference tissue models have been extended for kinetic analysis of dynamic PET with pharmacological challenges or cognitive activation during PET (Alpert et al., 2003; Votaw et al., 2002; Watabe et al., 1998; Zhou et al., 2006a).

On the other hand, in theory, there is a limitation in using reference tissue models for the tissue tracer kinetics of low or even negligible specific binding. This is because the low and negligible specific binding of [^{11}C]PIB in target tissue can result in a high linear relationship between the output and input for the reference tissue model. This situation may amplify the errors noticeably in estimates obtained through the reference tissue model. For example, a convergence problem in nonlinear fitting of SRTM to [^{11}C]doxepin ROI TAC occurred in tissues of low H_{1} receptor BP (Suzuki et al., 2005). It was also reported that the estimates of low BP obtained by conventional nonlinear SRTM fitting were not reliable in [^{11}C]PIB and [^{11}C]SB-13 PET studies (Verhoeff et al., 2004; Zhou et al., 2006b). Previous [^{11}C]PIB studies have reported that BP was close to 0, or the distribution volume ratio DVR (= 1+BP) was close to 1, in most cortical regions in controls. Estimates of PIB retention were also low in brain regions with negligible Aβ load in MCI and AD patients (Lopresti et al., 2005; Price et al., 2005). It has been demonstrated that the accuracy of estimates can be improved by simultaneously fitting a compartmental model with plasma input or reference tissue input to multiple ROI TACs (Buck et al., 1995; Cunningham et al., 2004, Endres et al., 2003; Ginovart et al., 2001; Raylman et al., 1994; Zhou et al., 2006b). For this approach, the coupled parameter or parameters can be estimated simultaneously for all ROIs. In contrast, the coupled parameter estimated by fitting the model to each single ROI TAC, the conventional ROI TAC model fitting procedure, usually varies among ROIs and this is not consistent with model assumption.

The use of higher order reference tissue models, i.e., more compartments used for target and reference tissues, has been proposed to reduce the bias of BP or DVR estimates from the SRTM (Endres et al., 2003; Kropholler et al., 2006; Wu and Carson, 2002; Zhou et al., 2006c). In this study, three reference tissue models were used for ROI kinetic modeling and their estimates were compared to investigate if there are any significant improvements using higher order reference tissue models. To obtain reliable estimates of model parameters, reference tissue models with parameter coupling were derived and implemented by simultaneous fitting for ROI based quantification. As a comparison, a reference tissue model with the conventional single ROI TAC fitting method was also applied to same ROI data set.

The parametric images that represent both the spatial distribution and quantification of the physiological parameters are generated by fitting a tracer kinetic model to the measured individual pixel time activity curves. However, due to the inherent high noise level of pixel kinetics of PET, the parametric images generated by conventional linear or nonlinear fitting are usually less accurate than those obtained by model fitting with spatial-temporal analysis (Gunn et al., 2002; Kimura et al., 1999; Turkheimer et al., 2003; Zhou et al., 2002, 2003). In this study, a linear regression with spatial constraint algorithm (LRSC) we previously developed (Zhou et al., 2003) was applied to the SRTM model for pixel-wise quantification of [^{11}C]PIB kinetics. The R_{1} and DVR images generated by the SRTM with the LRSC algorithm were compared to the estimates from ROI kinetic analysis and pixel-wise Logan plot. The [^{11}C]PIB dynamic PET data for 28 controls and 6 individuals with MCI were used in the study.

The reference tissue model is derived from classical compartmental model theory by eliminating plasma input with reference tissue TAC. In clinical ligand-receptor PET studies, a 2-tissue compartmental model with plasma input (Fig. 1) is commonly used to fit the measured reversible tracer kinetics for both target and reference tissues (Huang, et al., 1986; Koeppe et al., 1991; Lammertsma et al., 1996; Mintun et al., 1984; Price et al., 2005). The tracer kinetics described by Fig. 1 are based on the following assumptions: 1) rapid equilibrium between free and nonspecific binding in target tissue is attained; 2) the concentrations of tracer are homogenous in vascular space (C_{P}), free plus nonspecific binding compartment (C_{F+NS}), and specific binding compartment (C_{S}) for target tissue, free and nonspecific binding compartments (C_{F}, and C_{F+NS}) for reference tissue; and 3) the transport of tracer between compartments has first order kinetics. Based on above assumptions, the tracer kinetics in target and reference tissues are described by the following differential equations:

A representative 2-tissue compartmental model used in ligand-receptor PET studies for target and reference tissues. The concentrations in vascular space (C_{P}), free and nonspecific binding compartment (C_{F+NS}), specific binding compartment (C_{S}) for target **...**

$$\frac{d{C}_{F+NS}\mathit{\left(}t\mathit{\right)}}{dt}={K}_{1}{C}_{P}(t)-({k}_{2}+{k}_{3}){C}_{F+NS}(t)+{k}_{4}{C}_{S}$$

(1)

$$\frac{d{C}_{S}\mathit{\left(}t\mathit{\right)}}{dt}={k}_{3}{C}_{F+NS}-{k}_{4}{C}_{S}(t)$$

(2)

$$\frac{d{C}_{F}\mathit{\left(}t\mathit{\right)}}{dt}={K}_{1R}{C}_{P}(t)-({k}_{2R}+{k}_{5}){C}_{F}(t)+{k}_{6}{C}_{NS}$$

(3)

$$\frac{d{C}_{NS}\mathit{\left(}t\mathit{\right)}}{dt}={k}_{5}{C}_{F}-{k}_{6}{C}_{NS}(t)$$

(4)

Reference tissue model assumes that the tissue tracer activity contributed from vascular space is negligible, i.e., C_{T} = C_{F+NS} + C_{S}, and C_{R} = C_{F} + C_{NS}, where the C_{T} and C_{R} are the tracer concentrations in target and reference tissues measured by the PET scanner, C_{F+NS}(0) = C_{S}(0) = C_{F}(0) = C_{NS}(0)=0, *K*_{1} (ml/min/ml) is the transport rate constant from vascular space to target tissue, *k*_{2} (1/min) is the efflux rate constant from free plus nonspecific compartment to blood, *k*_{3} (1/min) is the rate of specific receptor binding, and *k*_{4} (1/min) is rate of dissociation from receptors, *K*_{1R} (ml/min/ml) is the transport rate constant from vascular space to reference tissue, *k*_{2R} (1/min) is the efflux rate from free compartment in reference tissue to blood, *k*_{5} (1/min) is the rate constant of nonspecific receptor binding, and *k*_{6} (1/min) is rate of dissociation from nonspecific binding. One common measure of tracer binding kinetics is the distribution volume (DV). The tracer DV in tissue or compartment is defined as the ratio of the tracer concentration in tissue or compartment to the tracer concentration in plasma at equilibrium condition. The primary measure for quantification of ligand-receptor dynamic PET is BP that is defined as BP = f_{2}B′_{max}/K_{D}, where f_{2} is the free fraction of tracer in the free and nonspecific binding compartment, B′_{max} (nM) is the available receptor density for tracer binding, and K_{D} (nM) is the tracer equilibrium dissociation constant. BP is an index of tracer specific binding to receptor (Huang et al., 1986; Mintun et al., 1984). In terms of model micro-parameters, BP in target tissue, DV in free plus nonspecific binding compartment (DV_{F+NS}), in target tissue (DV_{T}), and in reference tissue (DV_{REF}) can be expressed as BP = *k _{3}*/

The DVR or BP can be estimated directly by reference tissue models using the reference tissue TAC as input. The reference tissue model derived from a 2 compartments for both target and reference tissues (Fig. 1) has 7 parameters, and it usually results in model identity problems in clinical situations (Kropholler et al., 2006, Wu and Carson 2002). The following three reference tissue models with lower orders of model configuration described in Fig. 2 were rederived for models incorporating parameter coupling and were compared in the present study.

A conventionally employed full reference tissue model with 4 parameters (R_{1}, k_{2}, k_{3}, k_{4}) is shown in Fig. 2 A (Lammertsma et al., 1996). Under the assumption that rapid equilibrium is attained between free and nonspecific binding in reference tissue, the tracer concentration in reference tissue (C_{R}) is modeled with a single compartment. The reference tissue model with four parameters (R_{1}, k′_{2R}, k_{3}, k_{4}) (RTM4P) is then derived from 2 compartments for target tissue and 1 compartment for the reference tissue as below.

The tracer kinetics in target and reference tissue described by Fig. 2 A follow Eqs. (1)-(2) and Eq. (5) as below.

$$\frac{d{C}_{R}\mathit{\left(}t\mathit{\right)}}{dt}={K}_{1R}{C}_{p}(t)-{k\text{'}}_{2R}{C}_{R}(t)$$

(5)

where k′_{2R} (1/min) is the efflux rate from reference tissue to blood. Based on the assumption on DV_{F+NS}, i.e., K_{1}/k_{2} = K_{1R}/k′_{2R}. Let R_{1}=K_{1}/K_{1R}, we have k_{2} = (K_{1}/ K_{1R}) k′_{2R} = R_{1} k′_{2R}, By applying a Laplace transform to Eqs (1)-(2) and Eq. (5) with initial conditions C_{F+NS}(0) = C_{S}(0) = C_{R}(0) = 0, the operational equation for RTM4P can be expressed by parameters of R_{1}, k′_{2R}, k_{4} and BP:

$$\begin{array}{ll}{C}_{T}(t)& ={f}_{\mathit{RTM}\mathrm{4P}}& ={R}_{1}{C}_{R}+{C}_{R}(A\mathrm{\mathrm{exp}(\alpha t)+B\mathrm{\mathrm{exp}(\beta t))}}\end{array}$$

(6)

where
$\alpha =(-({k}_{2}+{k}_{3}+{k}_{4})+\sqrt{{({k}_{2}+{k}_{3}+{k}_{4})}^{2}-4{k}_{2}{k}_{4}})/2$,
$$\beta =(-({k}_{2}+{k}_{3}+{k}_{4})-\sqrt{{({k}_{2}+{k}_{3}+{k}_{4})}^{2}-4{k}_{2}{k}_{4}})/2$$,
$$A=\frac{{k}_{2}{k}_{3}+{k}_{2}(1-{R}_{1})(\alpha +{k}_{4})}{\alpha -\beta},\mathrm{B=\frac{-{k}_{2}{k}_{3}-{k}_{2}(1-{R}_{1})(\beta +{k}_{4})}{\alpha -\beta}}$$, k_{2} = R_{1} k′_{2R}, k_{3} = BP k_{4}, and is the mathematical operation of convolution. The k′_{2R} and k_{4} are the coupled parameters among all ROIs for RTM4P.

Symmetrically, the configuration for the second reference tissue model is usually referred to as the “Watabe” reference tissue model with parameters (R_{1}, k′_{2}, k_{2R}, k_{5}, k_{6}) (see Fig. 2 B (Endres et al., 2003; Gunn et al., 2001; Kropholler et al., 2006; Millet et al., 2002; Watabe et al., 2000). With the assumption of rapid equilibrium between free plus nonspecific binding and specific binding, tracer concentration C_{T} for target tissue can be modeled by a single compartment. Therefore, in addition to the Eq. (3)-(4) for tracer kinetics in reference tissue, the tracer kinetics in target tissue follows Eq. (7).

$$\frac{d{C}_{T}\mathit{\left(}t\mathit{\right)}}{dt}={K}_{1}{C}_{P}(t)-{k\text{'}}_{2}{C}_{T}(t)$$

(7)

where k′_{2} is the transport rate constant from target tissue to blood. Note that k′_{2} = k_{2}/(1 + BP), and k′_{2R} = k_{2R}/(1 + NS) (NS = k_{5}/k_{6}). The operational equation for RTM5P is then expressed by parameters of R_{1}, k_{2R}, NS, k_{6}, and BP as below.

$$\begin{array}{ll}{C}_{T}(t)& ={f}_{\mathit{\text{RTM}}5P}(t{R}_{1},{k}_{2R},NS,{k}_{6},BP)& & ={R}_{1}{C}_{R}+{C}_{R}(P\mathrm{\mathrm{exp}(-({k}_{5}+{k}_{6})t)+Q\mathrm{\mathrm{exp}(-{k\text{'}}_{2}t)})}\end{array}$$

(8)

where *k*_{5}=*k*_{6}*NS*,
${k\text{'}}_{2}=\frac{{R}_{1}{k}_{2R}}{(1+BP)(1+NS)}$,
$P=\frac{{R}_{1}{k}_{2R}{k}_{5}}{{k}_{5}+{k}_{6}-{k\text{'}}_{2}}$, and
$Q={R}_{1}({k}_{2R}-{k\text{'}}_{2})-P$. For RTM5P, the parameters k_{2R}, NS, k_{6} are coupled among all ROIs.

The last reference tissue model with configuration demonstrated by Fig. 2 C is commonly referred as the SRTM with parameters (R_{1}, k_{2}, and BP) (Lammertsma and Hume, 1996). The model with configuration Fig. 2 C assumes that the tracer concentrations in target and reference tissues can be modeled by a single compartment. The tracer kinetics in target and reference tissues for the SRTM are therefore determined by Eqs. (5) and (7). Based on DV_{F+NS} = DV_{REF}, we have k_{2} = (K_{1}/ K_{1R}) k′_{2R} = R_{1} k′_{2R}, and the SRTM with (R_{1}, k′_{2R}, BP) is referred to as RTM3P for consistency in this study. The operational equation for the RTM3P is then written as Eq. (3).

$$\begin{array}{ll}{C}_{T}(t)& ={f}_{\mathit{\text{RTM}}3P}(t{R}_{1},{k\text{'}}_{2R},BP)& & ={R}_{1}({C}_{R}+{k\text{'}}_{2R}(1-\frac{{R}_{1}}{1+BP}){C}_{R}\mathrm{\mathrm{exp}(-\frac{{R}_{1}{k\text{'}}_{2R}}{1+BP}t)})\end{array}$$

(9)

The k′_{2R} implies only one parameter coupled among all ROIs for RTM3P.

To embody the physiological assumptions into the model fitting process for the reference tissue model with parameter coupling, and to reduce the variability of estimates, simultaneously fitting a model to all ROI TACs has been used in ligand-receptor dynamic PET (Buck et al., 1995; Cunningham et al., 2004, Endres et al., 2003; Raylman et al., 1994; Zhou et al., 2006b, 2006c). For the above three reference tissue models with parameter coupling, the cost function to be minimized for RTM3P, RTM4P, and RTM5P are: $\sum _{j=1}^{N}\sum _{i=1}^{M}{w}_{i}{({C}^{j}({t}_{i})-{f}_{\mathit{\text{RTM}}3P}({t}_{i}{R}_{1}^{j},{k\text{'}}_{2R},B{P}^{j}))2}^{}$, $\sum _{j=1}^{N}\sum _{i=1}^{M}{w}_{i}{({C}^{j}({t}_{i})-{f}_{\mathit{\text{RTM}}4P}({t}_{i}{R}_{1}^{j},{k\text{'}}_{2R},{k}_{4,}B{P}^{j}))2}^{}$, and $\sum _{j=1}^{N}\sum _{i=1}^{M}{w}_{i}{({C}^{j}({t}_{i})-{f}_{\mathit{\text{RTM}}5P}({t}_{i}{R}_{1}^{j},{k}_{2R},NS,{k}_{6},B{P}^{j}))2}^{}$, respectively

Where M is the number of time frames for dynamic PET scans and N is number of ROI TACs, t_{i} is the mid time of ith frame of dynamic PET scanning, w_{i} is the duration of ith frame of dynamic PET scanning, *C ^{j}* (

For comparison, the parameters of the reference tissue models were also estimated by fitting RTM3P, RTM4P and RTM5P to each single ROI TAC. In contrast, the cost function for fitting model to individual jth ROI TAC is: $\sum _{i=1}^{M}{w}_{i}{({\mathrm{C}}^{\mathrm{j}}({t}_{i})-{f}_{\mathit{\text{RTM}}3P}({t}_{i}{R}_{1}^{j},{k\text{'}}_{2R}^{j},B{P}^{j}))2}^{}$, $\sum _{i=1}^{M}{w}_{i}{({\mathrm{C}}^{\mathrm{j}}({t}_{i})-{f}_{\mathit{\text{RTM}}4P}({t}_{i}{R}_{1}^{j},{k\text{'}}_{2R}^{j},B{P}^{j},{k}^{j}))2}^{}$, and $\sum _{i=1}^{M}{w}_{i}{({\mathrm{C}}^{\mathrm{j}}({t}_{i})-{f}_{\mathit{\text{RTM}}5P}({t}_{i}{R}_{1}^{j},{k}_{2R}^{j},N{S}^{j},{k}_{6}^{j},B{P}^{j}))2}^{}$ for RTM3P, RTM4P, and RTM5P, respectively.

To compare with previous results (Price et al., 2005; Lopresti et al., 2005), DVR was calculated as BP + 1 after model fitting. The Akaike information criterion (AIC) (Akaike 1976, Carson et al., 1993, Turkheimer et al., 2003) was calculated for simultaneous fitting. The AIC for the reference tissue model with conventional single ROI fitting method was calculated as $\sum _{j=1}^{N}\mathit{\text{AIC}}({\mathit{\text{ROI}}}^{j})$.

Based on the results from ROI kinetic analysis (the estimates of R_{1}, and DVR from RTM3P, RTM4P and RTM5P were almost same, see Comparison of reference tissue models subsection of Results section), the SRTM model was used to generate parametric images. To improve the accuracy of pixel-wise estimates, a linear regression with spatial constraint (LRSC) algorithm was applied to Eq. (10) and Eq. (11) to generate R_{1}, and DVR images, respectively (Zhou et al., 2003).

$${C}_{T}(t)={R}_{1}{C}_{\mathit{\text{REF}}}(t)+{k}_{2}\underset{0}{\overset{t}{\int}}{C}_{\mathit{\text{REF}}}(s)ds-{k}_{2}^{\text{'}}\underset{0}{\overset{t}{\int}}{C}_{T}(s)ds$$

(10)

$$\underset{0}{\overset{t}{\int}}{C}_{T}(s)ds=\mathit{\text{DVR}}\underset{0}{\overset{t}{\int}}{C}_{\mathit{\text{REF}}}(s)ds+{P}_{1}{C}_{\mathit{\text{REF}}}(t)-{P}_{2}{C}_{T}(t)$$

(11)

To perform pixel-wise statistical analysis, all the parametric images were spatially normalized to the standard space (pixel size: 2×2 mm^{2}, slice thickness 2 mm) using SPM2 (statistical parametric mapping software; Wellcome Department of Cognitive Neurology, London, UK). Because the *R*_{1} images contain greater structural information, the *R*_{1} images generated by LRSC were used to determine the parameters of spatial normalization. These transformation parameters were applied to all generated parametric images for each subject. Two iterations of the spatial normalization process were performed: 1) the parameters obtained by normalizing *R*_{1} images to the R_{1} template generated in our previous [^{11}C]raclopride study (Zhou et al., 2003), and 2) mean R_{1} images obtained by first iteration were used as a template for the second iteration.

The standard Logan plot with reference tissue input was used for estimation of the DVR of [^{11}C]PIB binding in previous studies (Lopresti et al., 2005, Mintun et al., 2006). As compared to the standard Logan plot with reference tissue input, one advantage of the simplified Logan plot with reference tissue input is that it eliminates the need to estimate the mean of k′_{2R} (Logan et al., 1996). The simplified Logan plot given by Eq. (12) was proposed for ROI based [^{11}C]PIB kinetic analysis in this study.

$$\frac{\underset{0}{\overset{\mathrm{t}}{\int}}{\mathrm{C}}_{\mathrm{T}}(\mathrm{s})\mathrm{ds}}{{\mathrm{C}}_{\mathrm{T}}(\mathrm{t})}=\mathit{\text{DVR}}\frac{\underset{0}{\overset{\mathrm{t}}{\int}}{\mathrm{C}}_{\text{REF}}(\mathrm{s})\mathrm{ds}}{{\mathrm{C}}_{\mathrm{T}}(\mathrm{t})}+\delta \mathrm{\text{for}\mathrm{\mathrm{t}\mathrm{t}}}$$

(12)

To obtain robust DVR estimates for pixel TACs of high noise levels, the simplified Logan plot in bilinear form as below was used to generate DVR images.

$$\underset{0}{\overset{\mathrm{t}}{\int}}{\mathrm{C}}_{\mathrm{T}}(\mathrm{s})\mathrm{ds}=\mathit{\text{DVR}}\underset{0}{\overset{\mathrm{t}}{\int}}{\mathrm{C}}_{\text{REF}}(\mathrm{s})\mathrm{ds}\mathrm{+\mathrm{\delta \mathrm{{\mathrm{C}}_{\mathrm{T}}(\mathrm{t})\mathrm{\mathrm{\text{for}\mathrm{\mathrm{t}\mathrm{t}}}}}}}$$

(13)

where t* = 40 min.

For comparison, the standard Logan plot using Eq. (14) was also applied to the ROI TACs to estimate ROI DVR in the study.

$$\frac{\underset{0}{\overset{\mathrm{t}}{\int}}{\mathrm{C}}_{\mathrm{T}}(\mathrm{s})\mathrm{ds}}{{\mathrm{C}}_{\mathrm{T}}(\mathrm{t})}=\mathit{\text{DVR}}\frac{\underset{0}{\overset{\mathrm{t}}{\int}}{\mathrm{C}}_{\text{REF}}(\mathrm{s})\mathrm{ds}+\frac{{\mathrm{C}}_{\text{REF}}(\mathrm{t})}{\overline{{k\text{'}}_{2\text{R}}}}}{{\mathrm{C}}_{\mathrm{T}}(\mathrm{t})}+\delta \mathrm{\mathrm{\text{for}\mathrm{\mathrm{t}\mathrm{t}}}}$$

(14)

where
$\overline{{k\text{'}}_{2R}}$ is a population average value of k′_{2R}. The
$\overline{{k\text{'}}_{2R}}$ estimated by a 2-tisue compartment, 4-parameter model with plasma input was 0.08/min for controls (n = 5), 0.07/min for MCI (n = 5) and AD (n = 4) (Price et al., 2005). The
$\overline{{k\text{'}}_{2R}}$ was also fixed at 0.149/min and 0.2/min for Logan plot with reference tissue input in [^{11}C]PIB PET studies by Lopresti and Mintun (2006), respectively. Note that the
$\overline{{k\text{'}}_{2R}}$ estimated by RTM3P for controls (n = 28) and MCIs (n = 6) was 0.05/min in our current study (see Results section). To investigate the effects of
$\overline{{k\text{'}}_{2R}}$ on the DVR estimates, we estimated DVR with a fixed series of
$\overline{{k\text{'}}_{2R}}$ values from 0.03/min to 0.5/min and compared these to the simplified Logan plot approach given by Eq. (12).

Subjects were 34 of the first 36 participants (excluding one with a clinical stroke and one with missing PET time frames due to scanner error) evaluated with [^{11}C]PIB as part of the ongoing neuroimaging substudy of the Baltimore Longitudinal Study of Aging (BLSA) (Resnick et al., 2000, 2003). At initial enrollment, BLSA neuroimaging participants were free of dementia and other central nervous system disorders, severe cardiovascular disease, and metastatic cancer (detailed in Resnick et al., 2000). [^{11}C]PIB studies were initiated in June 2005, and participants had been followed for up to 12 years with structural and functional imaging studies. Evaluation of diagnostic status followed established BLSA procedures, using prospective follow-up information (detailed in Gamaldo et al 2006). All participants received a detailed physical examination, including medical history updates and laboratory screening, neuropsychological testing, and assessment by the Clinical Dementia Rating (CDR) (Morris et al., 2001) scale in conjunction with the [^{11}C]PIB study. The CDR scores were typically based on informant interviews (spouse, child, or close friend) conducted by a certified examiner. Participant data was reviewed at a consensus diagnostic conference if the Blessed-Information-Memory-Concentration (BIMC; Blessed et al., 1968) test score was 3 or above or if the informant or subject CDR score was 0.5 or above. Diagnoses were made at consensus diagnostic conferences using DSM-III-R (APA, 1987) criteria for dementia and the NINCDS-ADRDA criteria (McKhann et al., 1984) using neuropsychological diagnostic tests and clinical data. A diagnosis of mild cognitive impairment not meeting criteria for dementia was made for participants who had cognitive impairment (typically memory) but did not have functional loss in activities of daily living. One participant in the present study met stringent diagnostic criteria for mild cognitive impairment and five additional participants scored 0.5 on the CDR scale, reflecting very mild cognitive impairment. Normal controls had scores of 0 on the CDR and were considered cognitively normal by our diagnostic procedures. It is important to emphasize that the MCI participants in this study are identified within the context of prospective longitudinal follow-ups and represent very mild cases of cognitive impairment in contrast to those followed in other studies who typically present with memory complaints.

Structural magnetic resonance imaging (MRI) scans and [^{11}C]PIB dynamic PET were acquired for each participant. MRI and PET imaging are typically performed on the same visit, but a major renovation of the MRI research scanners coincided with the initial PIB imaging studies. MRI scans were acquired within 3 days of the PET scan for 16 participants and within 1 to 2 years for 15 participants. Excluding 3 individuals, where it was necessary to use an MRI obtained 4, 5, and 10.5 years prior, respectively, the mean (SD) interval between MRI and PET was 0.6 (0.8) years.

MRI scans for anatomic localization were performed on a 1.5 Tesla GE Signa system using a spoiled gradient recalled acquisition sequence (124 slices with image matrix 256×256, pixel size 0.93×0.93 mm, slice thickness 1.5 mm). Dynamic PET [^{11}C]PIB studies were performed on a GE Advance scanner. The PET scanning was started immediately after intravenous bolus injection of 14.45 ± 1.01 (n = 34, range 11.02 to 16.36) (mean ± SD hereafter for “±”) mCi [^{11}C]PIB with a specific activity of 4.29 ± 1.49 (n = 34, range 0.96 to 6.83) Ci/μmol. Dynamic scans were acquired in 3D mode with acquisition protocol of 4×0.25, 8×0.5, 9×1, 2×3, 14×5 min (total 90 min, 37 frames). To minimize motion during PET scanning, all participants are fitted with thermoplastic face masks for the PET imaging. Ten-minute ^{68}Ge transmission scans acquired in 2-D mode were used for attenuation correction of the emission scans. Dynamic images were reconstructed using filtered back projection with a ramp filter (image size 128×128, pixel size 2×2 mm, slice thickness 4.25 mm), which resulted in a spatial resolution of about 4.5 mm FWHM at the center of the field of view. MRIs were coregistered to the mean of the first 20 min dynamic PET images using SPM2 with a mutual information method. In addition to the reference region (cerebellum), fifteen ROIs (1: caudate, 2: putamen, 3: thalamus, 4: lateral temporal, 5: mesial temporal, 6: orbital frontal, 7: prefrontal, 8: occipital, 9: superior frontal, 10: parietal, 11: anterior cingulate, 12: posterior cingulate, 13: pons, 14: midbrain, 15: white matter) were manually drawn on the coregistered MRIs (Price et al., 2005, Lopresti et al., 2005) and copied to the dynamic PET images to obtain ROI TACs for kinetic analysis. The ROIs estimates were also obtained by applying ROIs to parametric images.

There were 28 individuals (17 males, 11 females, age range 55-92 years, 78.6 ± 8.1) with CDR = 0 that were classified as the normal control group, and 6 individuals (2 males, 4 females, age range 77-89, 83.0 ± 4.2) in the MCI group. The difference in age between the control and MCI groups was not statistically significant.

The variation in DVR from RTM3P estimates from conventional single ROI TAC fitting method (cost function determined by single ROI TAC) was reduced remarkably by the simultaneous fitting method in caudate, putamen, and most cortical regions in the control group. As demonstrated by Fig. 3, the coefficients of variation (CV), defined as 100mean/(standard deviation) of DVR estimates, from the conventional method were reduced from 77% to 89% in caudate, putamen, orbital frontal cortex, prefrontal cortex, occipital cortex, superior frontal cortex, parietal cortex, anterior cingulate cortex, and posterior cingulate cortex in control group. The DVR estimated by the simultaneous fitting method were comparable to those estimated by the conventional method for the ROIs of thalamus, lateral temporal, mesial temporal cortex, pons, midbrain, and white matter in control group. For the MCI group, an improvement in the DVR estimates for the simultaneous fitting method relative to the conventional method was only found in the occipital cortex.

The coefficient of variation of DVR estimates (=100* (mean/(standard deviation))) for the control group (n = 28) and MCI group (n = 6). The DVR estimates were obtained from the RTM3P model with simultaneous fitting and single ROI TAC fitting methods. **...**

The mean plus standard deviation of k′_{2R} of RTM3P is shown in Fig. 4. k′_{2R}, the efflux rate constant from reference tissue to blood, is expected to be the same for all ROIs in the RTM3P model but shows high variation in estimates from the conventional single ROI TAC fitting method. The estimates of k′_{2R} vary from 0.01 ± 0.02 for prefrontal cortex to 0.15 ± 0.05 for thalamus in the control group, and from 0.01 ± 0.01 for occipital cortex to 0.14 ± 0.05 for thalamus in the MCI group. If k′_{2R} is estimated by the mean over all 15 ROIs, i.e., k′_{2R} (mean) =
$(\sum _{i=1}^{15}{k\text{'}}_{2R}(RO{I}_{i}))/15$, the k′_{2R}(mean) is significantly lower than that estimated by the simultaneous fitting approach (paired t test, p < 0.01) for the control group. The k′_{2R} (mean) in the control group tends to be lower than that in the MCI group (0.04 ± 0.01 versus 0.05 ± 0.01, p = 0.09), while the coupling method yields greater similarity between k′_{2R} for the control and MCI groups (0.05 ± 0.01 versus 0.05 ± 0.01, p = 0.42). R_{1} estimates from the conventional and simultaneous fitting methods do not differ significantly (p = 0.77). It was also found that R^{1} estimates from the conventional fitting method were high linearly correlated with those from the simultaneous fitting method as R_{1}(simultaneous fitting) = 0.96R_{1}(conventional) + 0.03 (R^{2} = 0.95).

The mean plus standard deviation of k′_{2R} estimates from the RTM3P model with simultaneous fitting and single ROI TAC fitting methods for the control group (n = 28) and MCI group (n = 6). k′_{2R} is the efflux rate constant from reference **...**

The reduction of variation for k′_{2R} and DVR estimates in parametric space for the simultaneous fitting approach is at the cost of higher residual sum of squares in kinetic space (cost function values) or AIC as compared to the conventional single ROI fitting approach. For the RTM3P model, the AIC of model fitting with the simultaneous fitting approach is higher (6 ± 3)% than that from the conventional fitting method. Note that the better fit or lower AIC mostly occurs in ROIs of lower DVR (< 1.5). Representative TACs are shown for a control in Fig 5 panel A1 and A2 and an MCI individual in Fig. 5 panel B1 and B2. As demonstrated by Fig. 5, the fitted ROI TACs from the simultaneous fitting and single ROI TAC fitting methods were quite comparable visually for (6 ± 3)% difference in AIC. Fig. 5 also illustrates that, as DVR increases, for example in posterior cingulate cortex, from controls to MCI, the difference in the fitted curves between the simultaneous fitting approach and conventional method tends to be smaller.

The results presented above suggest that estimates of the reference tissue model from the conventional single ROI TAC fitting method show high variability under conditions of low binding. Thus, in this section we base our comparison of parameter estimates from different reference tissue models on the simultaneous fitting approach. The AICs from SRTM3P, SRTM4P, and SRTM5P were -3184.39 ± 205.68, -3194.23 ± 195.01, and -3194.08 ± 201.67, respectively. Compared to RTM3P, there were no significant reductions in AIC by RTM4P (paired t test, p = 0.84) or RTM5P (paired t test, p = 0.85). Thus, there was no significant improvement in model fitting in kinetic space from RTM4P and RTM5P, as compared to RTM3P. Consistently, Fig. 6 shows that the estimates of R_{1}, and DVR from RTM3P were almost identical to those estimated from RTM4P and RTM5P. The k′_{2R} estimates from all three models were 0.05 ± 0.01 (n = 34), and there were no significant differences for k′_{2R} among the three models (p = 0.52 for RTM3P versus RTM4P and p = 0.71 for RTM3P versus RTM5P). Note that the k′_{2R} was calculated as k_{2R} /(1+NS) after fitting for RTM5P. The coupled parameter k_{4} estimates of RTM4P were 0.56 ± 0.24, and 0.40 ± 0.22 for controls and MCI subjects, respectively, and not significantly different (p = 0.16). For RTM5P, the estimates of NS (= k_{5}/k_{6}) were 0.17 ± 0.30, and 0.13 ± 0.24 for control and MCI, respectively (p = 0.76, not significant). The estimates from RTM5P for the coupled parameter k_{6} were 0.80 ± 0.23 and 0.63 ± 0.45 for control and MCI, respectively (p = 0.40, not significant).

The linear correlations between the DVR from Logan plot (simplified version) and the standard Logan plot with pre-determined $\overline{{k\text{'}}_{2R}}$ were:

- $\text{DVR}(\text{Logon plot}\overline{{k\text{'}}_{2R}}=0.03)=0.99\mathrm{DVR}(\text{Logon plot})+0.01({\mathrm{R}}^{2}=0.95)$,
- $\text{DVR}(\text{Logon plot}\overline{{k\text{'}}_{2R}}=0.05)=0.99\text{DVR}(\text{Logon plot})+0.01({\mathrm{R}}^{2}=0.98)$,
- $\text{DVR}(\text{Logon plot}\overline{{k\text{'}}_{2R}}=0.10)=1.00\text{DVR}(\text{Logon plot})+0.00({\mathrm{R}}^{2}=1.00)$, and
- $\text{DVR}(\text{Logon plot}\overline{{k\text{'}}_{2R}}=0.50)=1.00\text{DVR}(\text{Logon plot})+0.00({\mathrm{R}}^{2}=1.00)$.

Thus, the DVR estimates from the Logan plot using the simplified version employed in this study are almost the same as those from the standard Logan plot with the $\overline{{k\text{'}}_{2R}}$in the reported range. This observation is consistent with previous findings that the $\overline{{k\text{'}}_{2R}}$effect on DVR estimates for the Logan plot is negligible (Lopresti et al., 2006, Mintun et al., 2006).

The estimates of DVR from RTM3P had high linear correlations with those from the simplified Logan plot as DVR(Logan plot) = 0.87 DVR(RTM3P) +0.15 with R^{2} = 0.91. Paired t-test showed no significant differences between the DVR estimates from Logan plot and those from RTM3P (p = 0.37).

The DVR images generated by Logan plot were visually comparable with those generated by SRTM with LRSC approach. The linear correlations between the ROI values calculated on the DVR images generated by Logan plot and those calculated on DVR images generated by SRTM with LRSC were:

DVR(SRTM with LRSC) = 1.03DVR(Logan plot) - 0.07 with R^{2} = 0.99 (n = 34*15=510). The ROI estimates of DVR from ROI TAC based Logan plot were identical to those calculated directly on DVR images generated by Logan plot:

DVR(ROI TAC) = 1.00DVR (ROI on parametric image) - 0.00 with R^{2} = 1.00.

The statistics of DVR images in standard space showed that the standard deviation of DVR pixel estimators from SRTM with LRSC method was about 9% on average lower than those from the Logan plot in control group, and 1% on average lower for the MCI group. The mean images of R_{1} and DVR for the control (n = 28) and MCI (n = 6) groups with the mean MRI (n = 34) images are shown by Fig. 7. The R_{1} images of the control group are visually similar to the R_{1} images of the MCI group. The R_{1} estimates from ROI TAC fitting by RTM3P had the following linear relationship with those directly from R_{1} images:

Pixelwise mean of R_{1} and DVR for control group (n = 28) and MCI group (n = 6). The simplified reference tissue model with linear regression with spatial constraint parametric imaging algorithm was used for generation of R_{1} and DVR images. The MRI is the **...**

R_{1}(ROI TAC fitting) = 0.99R_{1}(ROI on parametric image) + 0.01 with R^{2} = 0.98.

Compared to DVR images, R_{1} images provided gray-white matter contrast consistent with that shown on the MRIs. This suggests that 1) R_{1} can be used for MRI-PET coregistration; 2) R_{1} can be used to determine spatial normalization parameters and provide a template for spatial normalization, for both control and MCI groups. Fig. 6 demonstrates that 1) In controls, DVR shows the highest values in thalamus, brain stem, and white matter; 2) The DVRs of mesial temporal cortex, thalamus, occipital cortex, brain stem, and white matter were similar in the control and MCI groups, and 3) The DVRs of caudate, putamen, and cortical regions, including frontal, lateral temporal, parietal and cingulate, in the MCI group were higher than those in the control group. The quantitative comparison of estimates between control and MCI groups is given in the following section.

The comparison of ROI estimates between control and MCI groups is summarized in Table. 1. The R_{1} and k′_{2R} estimates from the control group were similar to those from the MCI group. The DVR estimates from RTM3P were significantly greater than 1 (or BP > 0) in control and MCI groups for all 15 ROIs. However, the DVR estimates from Logan plots were not significant greater than 1 (or BP > 0) in mesial temporal cortex for both control and MCI groups (p = 0.052 and 0.24 for control and MCI, respectively), and not significantly greater than 1 (or BP > 0) in prefrontal cortex for the control group (p = 0.057). The difference in DVR between control and MCI groups was consistent across different estimates. Based on RTM3P, the DVRs in frontal and cingulate cortex for the MCI group were 38% higher on average than those for the control group. The DVRs in caudate, putamen, and lateral temporal cortex for the MCI group was 25% higher on average than those for the control group. There were no significant differences in DVRs between the control and MCI groups for thalamus, mesial temporal cortex, occipital cortex, pons, midbrain and white matter, although there was a trend for the DVR for thalamus to be higher in the MCI compared to control group (p = 0.11). Similar statistical inferences were also obtained for Logan plot and parametric image approaches.

Three reference tissue models, RTM3P, RTM4P, and RTM5P, were compared for 28 controls and 6 individuals with MCI, who had been studied using [^{11}C]PIB dynamic PET studies. Compared to the RTM3P model, RTM4P and RTM5P models did not yield significant improvements in AIC or estimates. It is notable that these results were based on a sample with relatively low cerebral [^{11}C]PIB specific binding. To generalize this comparison of the three models, it will be important to evaluate the three reference tissue models for [^{11}C]PIB PET with AD patients in future studies. In general, the performance of reference tissue models is dependent on the tracer kinetics in target tissue and reference tissue. For example, the underestimation of BP from SRTM was significantly reduced by a reference tissue model derived from 1 compartment for target tissue and 2 compartments for reference tissue (equivalent to RTM5P) in [^{11}C]carfentanil (Endres et al., 2003) and [^{11}C]diprenorphine dynamic PET studies (Zhou et al., 2006c). As compared to a 2-tissue compartment model with plasma input, the underestimation of [^{11}C]flumazenil BP estimates was (15 ± 0.6)%, (1 ± 1)%, and (15 ± 0.5)%, for RTM3P, RTM4P, and RTM5P respectively (Zhou et al., 2006c).

Variability in the estimates of DVR and k′_{2R} was noticeably reduced by the simultaneous fitting approach. Simultaneous fitting of multiple ROI TACs to the model of coupled parameters can be viewed as an approach that applies spatial constraints in parametric space to the model fitting in kinetic space at the ROI level (Zhou et al., 2002, 2003). However, as the DVR increases, for example in tissues with high density of Aβ, these improvements tend to be smaller. In other words, the DVR estimates from the conventional SRTM (R_{1}, k_{2}, BP) model fitting tend to be similar to those from RTM3P(R_{1}, k′_{2R}, BP) with simultaneous fitting methods if DVR 1 (e.g., DVR ≥ 2 or BP ≥ 1).

To compare the data from our study to previous results, a concentration ratio (CR) based semi-quantitative method was applied to the measured dynamic PET data for DVR estimates. As demonstrated by Fig. 8, the DVR estimated by CR was more sensitive to the pre-determined time frame. Consistent with TACs in Fig. 6, the mean of first 20 min scan [^{11}C]PIB images appears to be an appropriate time frame in trade off between contrast (gray matter versus white matter) or structural information and image noise level. The spatial distribution of CR images tends to be stable in time frames after 40 min post tracer injection. The CR over frames from 40 to 90 minutes is higher (13 ± 9)% than the DVR estimated from the Logan plot. However, the CR and Logan plot are highly correlated:

The mean images of concentration ratio images (= target(pixel)/cerebellum (ROI)) for the time periods [0 20] to [60 90] for the control (n = 28) and MCI (n = 6) groups.

CR([40 90]) = 1.39 DVR(Logan plot) − 0.33 (R^{2} = 0.94).

The ROI based kinetic analysis showed that the statistical power to distinguish between control and MCI group for CR([40 90]) method is the same as that for the Logan plot and the RTM3P or SRTM model. These results were quite consistent with previous results (Lopresti et al., 2006).

It is important to note that [^{11}C]PIB kinetics are different or even opposite to those for FDG in terms image contrast and noise. As shown by Fig. 8 and TACs in Fig. 6, the images obtained with time frame [40 90] are of high noise levels, and the contrast between gray matter and white matter is not consistent with that in early phase [0 20]. Thus, images obtained in later time phase (such as [40 60] or [40 90]) are not appropriate for multi-modality image coregistration. Due to heterogeneity and uncertainty in [^{11}C]PIB in spatial accumulation over brain tissues in MCI or AD patients, and less structural information, the CR images from the later phase are not recommended for use in spatial normalization or as a template. Fig. 9 showed that the statistical power to distinguish between controls and MCI individuals by DVR estimates from RTM3P fitting tends to be stable if PET study time is more than 60 min. The difference in DVR from RTM3P between 60 min and 90 min study is less than 4% in the ROIs including caudate, putamen, thalamus, cortex, pons and midbrain for both controls and MCI. Compared to a 90 min study, the DVR in white matter from a 60 min study is higher by 12% for controls and 10% for MCI group on average. In addition to DVR estimates from RTM3P, the R_{1} estimates from SRTM model with linear regression provide relative tracer transport rate (from blood) information, and R_{1} images appear reliable for MRI-PET coregistration and spatial normalization. Taking these factors into consideration, the reference tissue model RTM3P is suggested for kinetic analysis, especially where shorter (< 90 min) PET imaging times are required.

The p values for t-tests between control and MCI groups as a function of PET study time. The t-tests are based on the DVR estimates from the RTM3P model with the simultaneous fitting method. The regions of interest (ROIs) are numbered as: 1: caudate, **...**

In summary, reference tissue models with parameter coupling were derived and implemented by simultaneous model fitting for ROI kinetic analysis. A previously developed parametric imaging algorithm, linear regression with spatial constraint for the SRTM model, was evaluated. For comparison, the Logan plot with reference tissue input was applied to both ROI kinetic analysis and parametric imaging. Twenty eight controls and six MCI participants, imaged with [^{11}C]PIB dynamic PET, were evaluated in this study. Compared to conventional individual ROI TAC fitting, the variation of DVR estimates was reduced by the simultaneous ROI fitting approach, especially in tissues of low or negligible specific binding in this group of individuals with generally low PIB retention. There were no significant differences in both ROI TAC fitting and DVR estimates between the RTM3P and the RTM4P or RTM5P with simultaneous fitting for parameter estimation. As it produces similar results to the more complex RTM4P and RTM5P models, the simpler RTM3P model is proposed for ROI TAC based kinetic analysis in studies using [^{11}C]PIB PET. Thus, we compared the R_{1} and DVR images generated by the SRTM with LRSC algorithm to those from the RTM3P ROI kinetic analysis and the DVR images from Logan plot. The RTM3P with simultaneous fitting method is shown to be a robust compartmental modeling approach that may be useful in [^{11}C]PIB PET studies to detect early markers of Alzheimer’s disease where specific ROIs have been hypothesized, especially in situtations where PIB retention is not high. In addition, the SRTM with LRSC algorithm is recommended for generation of R_{1} and DVR images for pixel-wise quantification of [^{11}C]PIB dynamic PET.

We thank the cyclotron, PET, and MRI imaging staff of the Johns Hopkins Medical Institutions; Andrew H. Crabb for data transfer and computer administration. This study was supported in part by the Intramural Research Program, National Institute on Aging, NIH and by N01-AG-3-2124. This work was presented in part at The 53nd Society of Nuclear Medicine Annual Conference, June 3-7, 2006, San Diego, California, U.S.A.

**Publisher's Disclaimer: **This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

- Akaike H. An information criteria (AIC) Math Sci. 1976;14:5–9.
- Alpert NM, Badgaiyan RD, Livni E, Fischman AJ. A novel method for noninvasive detection of neuromodulatory changes in specific neurotransmitter systems. Neuroimage. 2003;19(3):1049–1060. [PubMed]
- Blessed G, Tomlinson BE, Roth M. The association between quantitative measures of dementia and senile change in the cerebral gray matter of elderly subjects. Br J Psychiatry. 1968;114:797–811. [PubMed]
- Buck A, Westera G, vonSchulthess GK, Burger C. Modeling alternatives for cerebral carbon-11-iomazenil kinetics. J Nucl Med. 1995 1996 Apr;37(4):699–705. [PubMed]
- Buckner RL, Snyder AZ, Shannon BJ, LaRossa G, Sachs R, Fotenos AF, Sheline YI, Klunk WE, Mathis CA, Morris JC, Mintun MA. Molecular, structural, and functional characterization of Alzheimer’s disease: evidence for a relationship between default activity, amyloid, and memory. J Neurosci. 2005;25(34):7709–17. [PubMed]
- Cunningham VJ, Hume SP, Price GR, Ahier RG, Cremer JE, Jones AK. Compartmental analysis of diprenorphine binding to opiate receptors in the rat in vivo and its comparison with equilibrium data in vitro. J Cereb Blood Flow Metab. 1991;11(1):1–9. [PubMed]
- Edison P, Archer HA, Hinz R, Hammers A, Pavese N, Tai YF, Hotton G, Cutler D, Fox N, Kennedy A, Rossor M, Brooks DJ. Amyloid, hypometabolism, and cognition in Alzheimer disease: an [11C]PIB and [18F]FDG PET study. Neurology. 2007;68(7):501–8. [PubMed]
- Cunningham VJ, Matthews JC, Gunn RN, Rabiner EA, Gee AD. Identification and interpretation of microparameters in neuroreceptor compartmental models. Neuroimage. 2004;22(suppl 2):T13. [PubMed]
- Endres CJ, Slifstein M, Frankle G, Talbot PS, Laruelle M. Simultaneous modeling of multiple TACs with parameter coupling can be used to enforce the assumption of uniform non-specific binding when applying reference-tissue models, and potentially improve parameter estimation. 2003. Brain03 and BrainPET′03. J Cereb Blood Flow Metab. 2003;23(suppl)
- Fagan AM, Mintun MA, Mach RH, Lee SY, Dence CS, Shah AR, Larossa GN, Spinner ML, Klunk WE, Mathis CA, Dekosky ST, Morris JC, Holtzman DM. Inverse relation between in vivo amyloid imaging load and cerebrospinal fluid Abeta(42) in humans. Ann Neurol. 2006;59(3):512–9. [PubMed]
- Gamaldo A, Moghekar A, Kilada S, Resnick SM, Zonderman AB, O’Brien R. Effect of a clinical stroke on the risk of dementia in a prospective cohort. Neurology. 2006;67(8):1363–1369. [PubMed]
- Ginovart N, Wilson AA, Meyer JH, Hussey D, Houle S. Positron emission tomography quantification of [(11)C]-DASB binding to the human serotonin transporter: modeling strategies. J Cereb Blood Flow Metab. 2001;21(11):1342–1253. [PubMed]
- Gunn RN, Gunn SR, Cunningham VJ. Positron emission tomography compartmental models. J Cereb blood Flow Metab. 2001;21:635–652. [PubMed]
- Gunn RN, Gunn SR, Turkheimer FE, Aston JA, Cunningham VJ. Positron emission tomography compartmental models: a basis pursuit strategy for kinetic modeling. J Cereb Blood Flow Metab. 2002;22(12):1425–39. [PubMed]
- Herholz K, Wienhard K, Heiss W-D. Validity of PET studies in brain tumors. Cerebrovascular and Brain Metabolism Review. 1990;2:240–265. [PubMed]
- Huang SC, Barrio JR, Phelps ME. Neuroreceptor assay with positron emission tomography: equilibrium versus dynamic approaches. J Cereb Blood Flow Metab. 1986;6(5):515–521. [PubMed]
- Kemppainen NM, Aalto S, Wilson IA, Nagren K, Helin S, Bruck A, Oikonen V, Kailajarvi M, Scheinin M, Viitanen M, Parkkola R, Rinne JO. Voxel-based analysis of PET amyloid ligand [11C]PIB uptake in Alzheimer disease. Neurology. 2006;67(9):1575–80. [PubMed]
- Kimura Y, Hsu H, Toyama H, Senda M, Alpert NM. Improved signal-to-noise ratio in parametric images by cluster analysis. Neuroimage. 1999;9(5):554–61. [PubMed]
- Klunk WE, Engler H, Nordberg A, Wang Y, Blomqvist G, Holt DP, Bergstrom M, Savitcheva I, Huang GF, Estrada S, Ausen B, Debnath ML, Barletta J, Price JC, Sandell J, Lopresti BJ, Wall A, Koivisto P, Antoni G, Mathis CA, Langstrom B. Imaging brain amyloid in Alzheimer’s disease with Pittsburgh Compound-B. Ann Neurol. 2004;55(3):306–19. [PubMed]
- Klunk WE, Lopresti BJ, Ikonomovic MD, Lefterov IM, Koldamova RP, Abrahamson EE, Debnath ML, Holt DP, Huang GF, Shao L, DeKosky ST, Price JC, Mathis CA. Binding of the positron emission tomography tracer Pittsburgh compound-B reflects the amount of amyloid-beta in Alzheimer’s disease brain but not in transgenic mouse brain. J Neurosci. 2005;25(46):10598–606. [PubMed]
- Kropholler MA, Boellaard R, Schuitemaker A, Folkersma H, van Berckel BN, Lammertsma AA. Evaluation of reference tissue models for the analysis of [(11)C](R)-PK11195 studies. J Cereb Blood Flow Metab. 2006 in press. [PubMed]
- Koeppe RA, Holthoff VA, Frey KA, Kilbourn MR, Kuhl DE. Compartmental analysis of [11C]flumazenil kinetics for the estimation of ligand transport rate and receptor distribution using positron emission tomography. J Cereb Blood Flow Metab. 1991;11:735–744. [PubMed]
- Lammertsma AA, Bench CJ, Hume SP, Osman S, Gunn K, Brooks DJ, Frackowiak RSJ. Comparison of Methods for Analysis of Clinical [
^{11}C]Raclopride Studies. J Cereb blood Flow Metab. 1996a;16:42–52. [PubMed] - Lammertsma AA, Hume SP. Simplified reference tissue model for PET receptor studies. NeuroImage. 1996;4:153–158. [PubMed]
- Logan J, Fowler JS, Volkow ND, et al. Graphical analysis of reversible radioligand binding from time-activity measurements applied to [N-11C-methuyl]-(-)-cocaine PET studies in human subjects. J Cereb blood Flow Metab. 1990;10:740–747. [PubMed]
- Logan J, Fowler JS, Volkow ND, Wang G-J, Ding Y-S, Alexoff DL. Distribution volume ratios without blood sampling from graphic analysis of PET data. J Cereb Blood Flow Metab. 1996;16:834–840. [PubMed]
- Lopresti BJ, Klunk WE, Mathis CA, Hoge JA, Ziolko SK, Lu X, Meltzer CC, Schimmel K, Tsopelas ND, DeKosky ST, Price JC. Simplified quantification of Pittsburgh Compound B amyloid imaging PET studies: a comparative analysis. J Nucl Med. 2005;46(12):1959–72. [PubMed]
- McKhann G, Drachman D, Folstein M, Katzman R, Price D, Stadlan EM. Clinical diagnosis of Alzheimer’s disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer’s Disease. Neurology. 1984;34(7):939–44. [PubMed]
- Mathis CA, Wang Y, Holt DP, Huang GF, Debnath ML, Klunk WE. Synthesis and evaluation of 11C-labeled 6-substituted 2-arylbenzothiazoles as amyloid imaging agents. Med Chem. 2003;46(13):2740–54. [PubMed]
- Mathis CA, Wang Y, Klunk WE. Imaging beta-amyloid plaques and neurofibrillary tangles in the aging human brain. Curr Pharm Des. 2004;10(13):1469–92. [PubMed]
- Mathis CA, Klunk WE, Price JC, DeKosky ST. Imaging technology for neurodegenerative diseases: progress toward detection of specific pathologies. Arch Neurol. 2005;62(2):196–200. [PubMed]
- Marquardt DW. An algorithm for least-squares estimations of nonlinear parameters. J Soc Ind App Math. 1963;11:431–441.
- Millet P, Graf C, Buck A, Walder B, Ibanez V. Evaluation of the reference tissue models for PET and SPECT benzodiazepine binding parameters. Neuroimage. 2002;17(2):928–942. [PubMed]
- Mintun MA. Utilizing advanced imaging and surrogate markers across the spectrum of Alzheimer’s disease. CNS Spectr. 2005;10(11 Suppl 18):13–6. [PubMed]
- Mintun MA, Raichle ME, Kilbourn MR, Wooten GF, Welch MJ. A quantitative model for the in vivo assessment of drug binding sites with positron emission tomography. Ann Neurol. 1984;15:217–227. [PubMed]
- Mintun MA, Larossa GN, Sheline YI, Dence CS, Lee SY, Mach RH, Klunk WE, Mathis CA, DeKosky ST, Morris JC. [
^{11}C]PIB in a nondemented population: potential antecedent marker of Alzheimer disease. Neurology. 2006;67(3):446–52. [PubMed] - Morris ED, Yoder KK, Wang C, Normandin MD, Zheng QH, Mock B, Muzic RF, Jr, Froehlich JC. ntPET: a new application of PET imaging for characterizing the kinetics of endogenous neurotransmitter release. Mol Imaging. 2005;4(4):473–89. [PubMed]
- Morris JC, Storandt M, Miller JP, McKeel DW, Price JL, Rubin EH, Berg L. Mild cognitive impairment represents early-stage Alzheimer disease. Arch Neurol. 2001;58(3):397–405. [PubMed]
- Nichols L, Pike VW, Cai L, Innis RB. Imaging and in vivo quantitation of β-amyloid: an exemplary biomarker for Alzheimer’s disease? Biol Psychiatry. 2006;59(10):940–947. [PubMed]
- Nordberg A. PET imaging of amyloid in Alzheimer’s disease. Lancet Neurol. 2004;3(9):519–527. [PubMed]
- Price JC, Klunk WE, Lopresti BJ, Lu X, Hoge JA, Ziolko SK, Holt DP, Meltzer CC, DeKosky ST, Mathis CA. Kinetic modeling of amyloid binding in humans using PET imaging and Pittsburgh Compound-B. J Cereb Blood Flow Metab. 2001;2005;25(11):1528–47. [PubMed]
- Raylman RR, Hutchins GD, Beanlands RS, Schwaiger M. Modeling of carbon-11-acetate kinetics by simultaneously fitting data from multiple ROIs coupled by common parameters. J Nucl Med. 1994 May;35(5):909–13. [PubMed]
- Reimold M, Mueller-Schauenburg W, Becker GA, Reischl G, Dohmen BM, Bares R. Non-invasive assessment of distribution volume ratios and binding potential: tissue heterogeneity and interindividually averaged time-activity curves. Eur J Nucl Med Mol Imaging. 2004;31(4):564–77. [PubMed]
- Resnick SM, Goldszal AF, Davatzikos C, Golski S, Kraut MA, Metter EJ, Bryan RN, Zonderman AB. One-year age changes in MRI brain volumes in older adults. Cereb Cortex. 2000;10(5):464–72. [PubMed]
- Resnick SM, Pham DL, Kraut MA, Zonderman AB, Davatzikos C. Longitudinal magnetic resonance imaging studies of older adults: a shrinking brain. J Neurosci. 2003;23(8):3295–301. [PubMed]
- Shoghi-Jadid K, Small GW, Agdeppa ED, Kepe V, Ercoli LM, Siddarth P, Read S, Satyamurthy N, Petric A, Huang SC, Barrio JR. Localization of neurofibrillary tangles and beta-amyloid plaques in the brains of living patients with Alzheimer disease. Am J Geriatr Psychiatry. 2002;10(1):24–35. [PubMed]
- Small GW, Kepe V, Ercoli LM, Siddarth P, Bookheimer SY, Miller KJ, Lavretsky H, Burggren AC, Cole GM, Vinters HV, Thompson PM, Huang SC, Satyamurthy N, Phelps ME, Barrio JR. PET of brain amyloid and tau in mild cognitive impairment. N Engl J Med. 2006;355(25):2652–63. [PubMed]
- Suzuki A, Tashiro M, Kimura Y, Mochizuki H, Ishii K, Watabe H, Yanai K, Ishiwata K, Ishii K. Use of reference tissue models for quantification of histamine H1 receptors in human brain by using positron emission tomography and [11C]doxepin. Ann Nucl Med. 2005;19(6):425–433. [PubMed]
- Turkheimer FE, Hinz R, Gunn RN, Aston JA, Gunn SR, Cunningham VJ. Rank-shaping regularization of exponential spectral analysis for application to functional parametric mapping. Phys Med Biol. 2003;48(23):3819–41. [PubMed]
- Turkheimer FE, Hinz R, Cunningham VJ. On the undecidability among kinetic models: from model selection to model averaging. J Cereb blood Flow Metab. 2003;23:490–498. [PubMed]
- Verhoeff NP, Wilson AA, Takeshita S, Trop L, Hussey D, Singh K, Kung HF, Kung MP, Houle S. In-vivo imaging of Alzheimer disease beta-amyloid with [11C]SB-13 PET. Am J Geriatr Psychiatry. 2004;12(6):584–95. [PubMed]
- Votaw JR, Howell LL, Martarello L, Hoffman JM, Kilts CD, Lindsey KP, Goodman MM. Measurement of dopamine transporter occupancy for multiple injections of cocaine using a single injection of [F-18]FECNT. Synapse. 2002;44(4):203–210. [PubMed]
- Watabe H, Endres CJ, Carson RE. Modeling methods for the determination of dopamine release with [
^{11}C]raclopride and constant infusion. Neuroimage. 1998;7:A57. - Wu C, Pike VW, Wang Y. Amyloid imaging: from benchtop to bedside. Curr Top Dev Biol. 2005;70:171–213. [PubMed]
- Wu Y, Carson RE. Noise reduction in the simplified reference tissue model for neuroreceptor functional imaging. J Cereb Blood Flow Metab. 2002;22(12):1440–52. [PubMed]
- Zhou Y, Cloughesy T, Hoh CK, Black K, Phelps ME. A modeling-based factor extraction method for determining spatial heterogeneity of Ga-68 EDTA kinetics in brain tumors. IEEE Transactions on Nuclear Science. 1997;44(6):2522–2527.
- Zhou Y, Huang SC, Bergsneider M, Wong DF. Improved parametric image generation using spatial-temporal analysis of dynamic PET studies. NeuroImage. 2002;15:697–707. [PubMed]
- Zhou Y, Endres CJ, Brasic JR, Huang SC, Wong DF. Linear regression with spatial constraint to generate parametric images of ligand-receptor dynamic PET studies with a simplified reference tissue model. Neuroimage. 2003;18(4):975–989. [PubMed]
- Zhou Y, Chen MK, Endres CJ, Ye W, Brasic JR, Alexander M, Crabb AH, Guilarte TR, Wong DF. An extended simplified reference tissue model for the quantification of dynamic PET with amphetamine challenge. Neuroimage. 2006a;33(2):550–563. [PubMed]
- Zhou Y, Resnick SM, Ye W, Fan H, Holt D, Klunk WE, Mathis CA, Dannals R, Wong DF. Evaluation of reference tissue models for quantification of [
^{11}C] Pittsburgh compound B dynamic PET studies. J Nucl Med. 2006b;47(S1):117P. - Zhou Y, Ye W, Brasic JR, Hilton J, Wong DF. Compartmental model with plasma input versus reference tissue model for quantification of liganreceptor dynamic PET studies. J Nucl Med. 2006c;47(S1):117P.

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |