Simulating the effect of cerebral blood flow changes on regional quantification of [ 18 F]flutemetamol and [ 18 F]florbetaben studies

Abstract

Global and regional changes in cerebral blood flow (CBF) can result in biased quantitative estimates of amyloid load by PET imaging. Therefore, the current simulation study assessed effects of these changes on amyloid quantification using a reference tissue approach for [¹⁸F]flutemetamol and [¹⁸F]florbetaben. Previously validated pharmacokinetic rate constants were used to simulate time-activity curves (TACs) corresponding to full dynamic and dual-time-window acquisition protocols. CBF changes were simulated by varying the tracer delivery (K₁) from +25 to −25%. The standardized uptake value ratio (SUVr) was computed and TACs were fitted using reference Logan (RLogan) and the simplified reference tissue model (SRTM) to obtain the relative delivery rate (R₁) and volume of distribution ratio (DVR). RLogan was least affected by CBF changes (χ² = 583 p < 0.001, χ² = 81 p < 0.001, for [¹⁸F]flutemetamol and [¹⁸F]florbetaben, respectively) and the extent of CBF sensitivity generally increased for higher levels of amyloid. Further, SRTM-derived R₁ changes correlated well with simulated CBF changes (R² > 0.95) and SUVr’s sensitivity to CBF changes improved for later uptake-times, with the exception of [¹⁸F]flutemetamol cortical changes. In conclusion, RLogan is the preferred method for amyloid quantification of [¹⁸F]flutemetamol and [¹⁸F]florbetaben studies and SRTM could be additionally used for obtaining a CBF proxy.

Introduction

Amyloid-beta accumulation (Aβ) in the brain is one of the hallmarks of Alzheimer’s disease (AD). It can be visualised and quantified using positron emission tomography (PET) and both static or dynamic scanning protocols can be used.^1,2 The static protocol has the advantage of a short scan duration together with relatively simple processing and analytical steps, while the dynamic protocol provides higher accuracy at the cost of a much longer scan, more complex processing and advanced kinetic analysis. The semi-quantitative parameter obtained from a static scan, the standardized uptake value ratio (SUVr), depends on post-injection starting time and duration of the acquisition, and may be affected by changes in cerebral blood flow (CBF) as demonstrated for [¹¹C]PiB.^3–6 In contrast, the non-displaceable binding potential (BP_ND), derived from a dynamic PET scan, is less sensitive to noise, more robust against CBF changes, and for [¹¹C]PiB it has been shown to be the parameter of choice when measuring longitudinal changes in amyloid burden.^5,6 This characteristic may be especially important in situations where changes in CBF can occur, such as when measuring disease progression or treatment response in clinical trials. As a compromise between the protocols mentioned above, dynamic data acquisitions from a dual-time-window protocol have gained attention, in which data are acquired separately for early and late phases of tracer uptake to reduce overall scanning time, maintain high quantitative accuracy and provide tracer delivery information.^7,8

Cerebral blood flow declines with age and differs per brain region. ⁹ Compared with young adults (25 years), elderly (late 70s) can present with up to 25% CBF reductions, with an average annual CBF decline in grey matter of approximately 0.5%.^9,10 On top of this global decline, day-to-day whole brain CBF fluctuations of around 30% have been reported in test–retest studies; however, this percentage was assumed to consist of both physiological as well as measurement error. ¹¹ Furthermore, drugs may also exert an effect on CBF,^12,13 thus potentially compromising the measurement of amyloid changes associated to pharmacological interventions. In the context of AD, additional relative CBF changes may be present, as focal reductions in CBF have been observed in several brain regions.^14–16 A study on non-demented older adults scanned with [¹¹C]PiB showed that subjects with elevated [¹¹C]PiB signals experienced greater relative CBF variations, mainly in those regions showing increased amyloid deposition. ¹⁷ Given that CBF differs per region, declines with age, and that additional regional CBF changes occur during the course of Alzheimer’s disease, it is important to assess whether and, if so, to which extent these changes affect quantification of amyloid load. Furthermore, the additional effects of acquisition start-time and duration on quantitative accuracy need to be understood, due to its relevance in clinical and research practice. ³

The effect of CBF changes on amyloid load quantification has been assessed previously for [¹¹C]PiB, [¹⁸F]florbetaben and semi-quantitative measures of [¹⁸F]florbetapir.^5,7,18 These studies showed that large cortical CBF reductions resulted in a maximum change of ±10% in SUVr. For large global CBF reductions, a larger bias was observed for [¹¹C]PIB as compared with [¹⁸F]florbetaben and [¹⁸F]florbetapir.^5,7,18 These between tracer differences are to be expected due to differences in tracer kinetics and corresponding equilibrium times. However, the effect of CBF changes on quantification of [¹⁸F]flutemetamol scans remains unknown, just as the potential effects of such changes on novel dual-time-window protocols. Therefore, the present simulation study aimed to assess the effects of regional (i.e. target and reference tissues) and global CBF changes on quantitative amyloid measures derived from static, dynamic and dual-time window scanning protocols using reference tissue approaches, ⁸ thereby focusing on [¹⁸F]flutemetamol and [¹⁸F]florbetaben, the two tracers used within the AMYPAD consortium. A second aim was to assess whether the relative delivery rate (R₁) of both ligands could be used to accurately monitor changes in CBF.^16,19

Materials and methods

Subjects and PET data

Previously reported clinical PET data from 6 [¹⁸F]flutemetamol and 20 [¹⁸F]florbetaben subjects (both control and AD subjects), consisting of regional time-activity curves (TACs) and whole blood and/or plasma input curves, were used for this simulation study.^20,21 All participants provided written informed consent in accordance with the Declaration of Helsinki. The Ethical Committee of the University Hospitals Leuven approved the study protocol for the [¹⁸F]flutemetamol study (EudraCT 2007-000784-19, Registered 8 February 2007). For [¹⁸F]florbetaben, the local Institutional Review Board of University of Leipzig, the National Radiation Safety Committee, and the German Federal Institute for Drugs and Medical Devices approved the study protocol (EudraCT 2006-003882-15, Registered 2006). The first group received a bolus injection (<40 s) of 181 ± 5 MBq [¹⁸F]flutemetamol and were scanned on a Siemens HiRez Biograph 16 PET/CT scanner.^3,20 The second group received an injection of 300 ± 60 MBq [¹⁸F]florbetaben (over 90 s) and were scanned on a ECAT HR+ Siemens/CTI scanner. ²¹ As described previously, whole blood curves and metabolite-corrected plasma input curves were available for [¹⁸F]flutemetamol, and metabolite-corrected plasma input curves and discrete whole blood samples for [¹⁸F]florbetaben. ⁸ Subsequently, regional TACs were fitted using the reversible two-tissue compartment model (four rate constants) with additional blood volume fraction parameter (2T4k_V_b) to obtain pharmacokinetic rate constants for the target (composite cortical region consisting of anterior and posterior cingulate, frontal, parietal, and lateral and medialtemporal cortex) and reference tissue (grey matter cerebellum). ⁸

TAC simulations

Based on the pharmacokinetic rate constants estimated from the clinical data mentioned above, realistic target and reference tissue TACs of 130 min duration were simulated using the 2T4k_V_b model (see Figure 1 and Table 1 for the rate constants). ⁸ Target tissue TACs were simulated for a clinically observed range of BP_ND values (captured by five simulated BP_ND values, Figure 1), and noise (0, 1 and 2%) corresponding to regions of interest (ROI) of various sizes was added to the target tissue TACs. ²² This resulted in a total of 50 TACs per noise level (identical TACs in case of 0% noise) and for each simulated BP_ND (see Figure 1 for a flowchart of the method). No noise was added to reference tissue TACs due to its relatively large volume.

Figure 1.

Schematic overview of the applied methods. 1. [¹⁸F]flutemetamol DVR_SIM (simulated volume of distribution ratio) ranged from 1.022 to 1.778, [¹⁸F]florbetaben DVR_SIM ranged from 1.026 to 2.051. 2. All steps were repeated for 1 and 2% noise. 3. All values are minutes post injection (p.i.). *Here DVR refers to BP_ND+1. Abbreviations: V_b: blood volume; BL: baseline (no cerebral blood flow change); SRTM: simplified reference tissue model; RLOGAN: reference Logan, SUVr: standardised uptake value ratio; FLOW: with changes in cerebral blood flow.

Table 1.

Rate constants used for CBF simulations.

	Global		Local cortical		Local cerebellar
	Target	Reference	Target	Reference	Target	Reference
[18F]flutemetamol
K1	0.186–0.310	0.240–0.400	0.186–0.310	0.320 a	0.248 a	0.240–0.400
k2	0.060–0.100	0.077–0.129	0.060–0.010	0.103 a	0.080 a	0.077–0.129
k3	0.008–0.028	0.018	0.008–0.028	0.018	0.008–0.028	0.018
k4	0.020	0.050	0.020	0.050	0.020	0.050
[18F]florbetaben
K1	0.170–0.283	0.188–0.313	0.170–0.283	0.250 a	0.226 a	0.188–0.313
k2	0.052–0.087	0.057–0.095	0.052–0.087	0.076 a	0.069 a	0.057–0.095
k3	0.010–0.030	0.007	0.010–0.030	0.007	0.010–0.030	0.007
k4	0.010	0.007	0.010	0.007	0.010	0.007

^aBaseline parameter, corresponding to a 0% CBF change. Units: K₁ in ml × g⁻¹ × min⁻¹ and k₂ in min⁻¹.

Note: For both tracers, the blood volume fraction parameter (V_b) was set to 0.05.

Simulating CBF changes

Given that the rate constant for tracer delivery (K₁) can be considered a proxy for blood flow as long as the extraction remains constant ( $K_1=\mathit{Flow}\times \mathit{Extraction}$ ), changes in CBF were simulated by varying both K₁ and tissue clearance (k₂), in order to maintain the non-displaceable volume of distribution constant $(V_{\mathrm{ND}}=K_1/k_2)$ . ²³ Both increases and decreases were simulated (from +25 to −25%) in three different scenarios (see Table 1 and Figure 2 for the simulation parameters and a schematic overview of the simulations):

Figure 2.

Schematic overview of the simulated CBF changes. C _P : tracer plasma radioactivity concentration, C _F+NS : free and non-specifically bound concentration in tissue, C _{S+NS_SLOW} : specifically bound and slow non-specifically bound concentration in tissue, all units Bq*ml⁻¹. Units of the rate constants: K₁ in ml × g⁻¹ × min⁻¹ and k₂, k₃, k₄ in min⁻¹.

1. Global CBF changes: equal variations in K₁ and k₂ for target and reference tissues.

2. Cortical CBF changes: variations in K₁ and k₂ only for the target tissue, reference tissue parameters were kept constant.

3. Cerebellar CBF changes: variations in K₁ and k₂ only for the reference tissue, target tissue parameters were kept constant.

Dual-time-window TACs

In addition to the full 130 min TACs, dual-time-window TACs (here called ‘interval TACs’) were created by removing data points from target and reference tissue TACs according to the dual-time-window protocol. ⁸ In short, the first 110 min were used (corresponding to a dynamic scanning protocol) and data points between the early (0–30 min post injection (p.i.)) and late acquisition phase (90–110 min p.i.) were deleted.

Estimating parameters of interest

SUVr was calculated from the 130 min full TACs, for starting times (t*) ranging from 70 to 110 min p.i. (interval t*: 3× 5 min, 4 × 2.5 min, 3 × 5min), and durations of 10, 15 and 20 min. Next, 110 min of the full TACs (corresponding to a dynamic scanning protocol) and interval TACs were fitted using reference Logan (RLogan) to estimate the volume of distribution ratio (DVR). ²⁴ The implementation of this model did not require fixing k₂ʹ (as per equation (7)). ²⁴ The linearization start times (t*) evaluated ranged from 50 to 80 min p.i. (in steps of 10 min) and interpolation of interval target tissue TACs was performed using cubic interpolation, as the method is routinely used for parametric imaging. Finally, 110 min of the full TACs and the interval TACs were fitted using the simplified reference tissue model (SRTM) to estimate R₁, k₂ and BP_ND. ²⁵ Given that SRTM requires a continuous input to fit the target tissue TACs, missing data points of the reference tissue interval TACs were interpolated using the 2T4k_V_b model together with a typical, tracer specific input function as validated previously. ⁸ In addition, for both tracers, boundary values (optimized for full, noiseless TACs) were set for all kinetic parameters to prevent physiologically implausible results (Supplementary Figure 1). Optimal boundaries were defined based on the simulated parameter range and lower boundaries were fine-tuned (based on parameter histograms) in case of k₂. This procedure selected the boundary value that resulted in least bias (calculated from the simulated BP_ND+1) across CBF scenarios and amyloid binding levels

Statistical analysis

For all analyses (i.e. all amyloid levels, starting times, durations, scanning protocols and for TACs both with and without noise), parameter estimates obtained in the absence of CBF variations were used as baseline parameters (SUVr_{_BL}, R_{1_BL}, DVR_{BL_SRTM} and DVR_{BL_RLOGAN}, Table 1). Next, percentage change as a result of CBF variations was calculated relative to each baseline for ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ , R_{1_FLOW}, ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ and RLogan ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}:$

$${\mathrm{\Delta }\mathrm{PAR}}_{\mathrm{FLOW}}\mathrm{ }\left(\mathrm{\%}\right)\hspace{0.17em}=\hspace{0.17em}\frac{{\mathrm{ }{\mathrm{PAR}}_{\mathrm{FLOW}}-\mathrm{ }\mathrm{PAR}}_{\mathrm{BL}\mathrm{ }}^{\mathrm{ }}}{{\mathrm{PAR}}_{\mathrm{BL}}}·100\mathrm{\%}$$

where PAR corresponds to the parameter of interest (i.e. SUVr, DVR or R₁). These parameters were then used to assess the sensitivity of the methods to variations in CBF. This was done for different levels of amyloid load, different scanning protocols, noise levels and uptake times.

Non-parametric Kruskal–Wallis tests were used to assess differences in sensitivity to CBF variations between methods and simulated BP_ND values. ²⁶ In addition, the maximum percentage change due to CBF variations was compared between methods and for the whole spectrum of simulated BP_ND values. Furthermore, to assess the relationship between $\mathrm{\Delta }$ R_{1_FLOW} and simulated $\mathrm{\Delta }$ CBF, the coefficient of determination (R²) was calculated based on Pearson’s correlation coefficient. Next, differences in sensitivity were compared between full and interval TACs fitted with SRTM and RLogan using the post hoc Mann–Whitney U test and boxplot analyses, where outliers were defined as points outside the whiskers (created using the Tukey method ²⁷ ). Subsequently, effects of noise on sensitivity to CBF changes and its relationship with all parameters were evaluated based on boxplot analyses. Finally, using Kruskal–Wallis tests, effects of altering starting times (t*) on CBF sensitivity were evaluated for RLogan and altering starting times and acquisition duration on CBF sensitivity for SUVr.

Results

Evaluation of outcome parameter(s)

The resulting boundary values for [¹⁸F]flutemetamol were: R₁: lower boundary (LB) = 1e-6 upper boundary (UB) = 2.00, k₂: LB = 0.018 UB = 2.00, BP_ND: LB = −0.50 UB = 5.00. For [¹⁸F]florbetaben, boundary values were: R₁: LB = 1e-6 UB = 2.00, k₂: LB = 0.028 UB = 2.00, BP_ND: LB = −0.50 UB = 5.00.

Sensitivity to CBF changes – Comparison between methods

[¹⁸F]flutemetamol

Overall, ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ was least sensitive to changes in CBF (χ² = 583, p < 0.001). With respect to global CBF changes, all models showed less than 6% change. More specifically, ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ showed least overall sensitivity to global CBF changes (χ² = 72, p < 0.001), while ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ showed the smallest maximum change (3.7%) followed by ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ (5.1%) and ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ (max. 5.8%). For cortical and cerebellar changes, $\mathrm{ }{\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ was least sensitive to CBF changes (χ² = 995 p < 0.001, maximum regional change: 2.0% compared with ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ 13.4% and ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ 4.8%), as shown in Figure 3. Sensitivity to simulated CBF changes increased for higher simulated levels of amyloid load for global ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : χ² = 280, p < 0.001; ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : χ² = 104, p < 0.001; ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : χ² = 239, p < 0.001), cortical ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : χ² = 63, p < 0.001; ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : χ² = 118, p < 0.001; ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : χ² = 142, p <0.001) and cerebellar changes ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : χ²  = 161, p <0.001; ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : χ²  = 92, p < 0.001; ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : no significant difference).

Figure 3.

[¹⁸F]flutemetamol: sensitivity to CBF changes across methods for three levels of amyloid load. Low amyloid (blue): DVR = 1.022, intermediate amyloid (pink): DVR = 1.400, high amyloid (red): DVR = 1.778. Four values (between 10 and 10.6%) for SRTM cortical CBF changes high amyloid are not shown. RLogan linearization time was 50–110 min p.i., SUVr uptake time was 90–110 min p.i.

[¹⁸F]florbetaben

Overall, ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ was least sensitive to CBF changes compared with the other methods (χ²  = 81, p < 0.001), in particular for cerebellar changes (χ²  = 457, p < 0.001), and showed smallest maximum change ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}:$ 1.2%, ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}:$ 6.1%, ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : 2.4%). With respect to global CBF changes, all methods showed a very similar sensitivity pattern with CBF decreases resulting in lower and CBF increases in higher changes (maximum change: ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}:\mathrm{ }5.2\mathrm{\%},\mathrm{ }{\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : 7.5% and ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : 6.1%), as shown in Figure 4. With respect to cortical CBF changes, ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ was least sensitive (χ²  = 68. p < 0.001) and showed the smallest maximum change ( ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : 3.7%, ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : 4.7% and ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : 7.3%). Furthermore, in general, the extent of CBF sensitivity increased for higher levels of simulated amyloid load, for global ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : χ² = 290, p < 0.001; ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : χ² = 171, p < 0.001; ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : χ² = 274, p < 0.001), cortical ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : χ² = 131, p < 0.001; ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : χ² = 105, p < 0.001; SUVr: χ² = 184, p < 0.001) and cerebellar changes ( ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ : χ² = 13, p =  0.011: ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ : χ² = 201, p < 0.001; ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ : no significant difference).

Figure 4.

[¹⁸F]florbetaben: sensitivity to CBF changes across methods for three levels of amyloid load. Low amyloid (blue): DVR = 1.026, intermediate amyloid (pink): DVR = 1.538, high amyloid (red): DVR = 2.051. RLogan linearization time was 50–110 min p.i., SUVr uptake time was 90–110 min p.i.

Relative delivery and CBF changes

For [¹⁸F]flutemetamol, correlations between simulated ΔCBF and ΔR_{1_FLOW} were high for all amyloid levels both for cortical (R² = 0.98 to 1.00, p < 0.001) as well as cerebellar CBF changes (R² = 0.97 to 1.00, p < 0.001). For [¹⁸F]florbetaben, high correlations were also observed between simulated ΔCBF and ΔR_{1_FLOW} across amyloid levels both for cortical (R² = 0.98 to 0.99, p < 0.001) and cerebellar CBF changes (R² = 0.95 to 0.99, p < 0.001). Very similar results were obtained for R₁ estimates derived from interval TACs (data not shown).

Sensitivity to CBF changes – Full versus interval TACs

For [¹⁸F]florbetaben, full TACs were found to be less sensitive to CBF changes than interval TACs (p < 0.001), while for [¹⁸F]flutemetamol, the opposite effect was found (p < 0.001). Post hoc tests revealed that for [¹⁸F]florbetaben this effect was present for both models, while for [¹⁸F]flutemetamol the effect was only present for SRTM. Furthermore, for both tracers, most outliers were present for interval compared to full TACs (outliers interval TACs: [¹⁸F]flutemetamol: 91.67%, [¹⁸F]florbetaben: 66.67%). These outliers were only observed for TACs fitted with SRTM, predominantly corresponding to low simulated amyloid levels (Figure 5). For [¹⁸F]florbetaben, interval TAC-derived ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ was more sensitive to cerebellar CBF changes compared to full TAC-derived ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ , in case of intermediate (Figure 5) to high (data not shown) amyloid levels.

Figure 5.

Comparison of full and interval TACs and their sensitivity to CBF changes. Coloured boxes represent full- and uncoloured boxes interval TACs. Upper row [¹⁸F]flutemetamol, bottom row [¹⁸F]florbetaben. Note: whiskers were defined according to the Tukey method and outliers are depicted as red crosses.

Sensitivity to CBF changes – The effect of noise

Overall, effects of noise on the sensitivity of the different methods to CBF changes were minimal for full TACs (maximum difference: 1 percentage point (p.p.)), see Supplementary Figure 2). Furthermore, for interval TACs, more outliers and a larger sensitivity were observed at higher compared to lower noise levels.

Sensitivity to CBF changes – Acquisition start time and duration

Overall, later uptake times showed a decreased sensitivity to CBF changes for both tracers ([¹⁸F]flutemetamol: χ² = 54, p < 0.001, [¹⁸F]florbetaben: χ² = 127, p < 0.001). Post hoc tests revealed the effect was not present for [¹⁸F]flutemetamol cortical changes. Furthermore, CBF increases resulted in a higher ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ , while CBF decreases caused the opposite effect, except for [¹⁸F]flutemetamol cortical changes (Figure 6). Varying the duration of the uptake period from 20 to 15 or 10 min provided essentially identical results for [¹⁸F]flutemetamol. With respect to [¹⁸F]florbetaben, ${\mathrm{SUVr}}_{\mathrm{FLOW}}\mathrm{ }$ calculated over 20 min showed significantly less CBF sensitivity compared with ${\mathrm{SUVr}}_{\mathrm{FLOW}}\mathrm{ }$ calculated over 10 min (p =  0.005). With respect to $\mathrm{ }{\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ , there was a small effect of linearization start time with later start times (corresponding to shorter acquisitions) being less sensitive to CBF changes ([¹⁸F]flutemetamol: χ² = 78, p < 0.001, [¹⁸F]florbetaben: χ² =  8, p =  0.04). More specifically, global and cortical CBF changes resulted in a small decrease in CBF sensitivity for later linearization times ([¹⁸F]flutemetamol: χ² = 66, p < 0.001 and χ² =  49, p < 0.001, respectively; [¹⁸F]florbetaben: χ² = 36 and p < 0.001, χ² = 14 p =  0.003, respectively). On the other hand, for [¹⁸F]florbetaben, cerebellar CBF changes resulted in a small increase in the models sensitivity (χ² = 44, p < 0.001; Supplementary Figure 3).

Figure 6.

The effect of uptake time on sensitivity of SUVr to CBF changes. Upper row is [¹⁸F]flutemetamol and bottom row is [¹⁸F]florbetaben, depicted for intermediate amyloid (DVR = 1.400 and DVR = 1.538 for both tracers respectively) Red dots resemble CBF increases and blue dots CBF decreases.

Discussion

The current simulation study assessed the effect of global and regional CBF changes on regional amyloid quantification based on static, dynamic and dual-time-window scanning-protocols using reference tissue approaches for [¹⁸F]flutemetamol and [¹⁸F]florbetaben. The results of the present simulation study indicate that, compared with ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}\mathrm{ }$ and $\mathrm{ }{\mathrm{SUVr}}_{\mathrm{FLOW}}$ , ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}\mathrm{ }$ was less affected by changes in CBF for both [¹⁸F]flutemetamol and [¹⁸F]florbetaben. Furthermore, this sensitivity to CBF changes increased with increasing levels of amyloid. Finally, changes in the relative delivery rate R₁, obtained with SRTM, were highly correlated with simulated CBF changes.

The finding that ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ was robust against changes in CBF is in line with a previous study that found the same result for [¹¹C]PiB ⁵ and this robustness may be due to the linearity of the model. ²⁴ In contrast, SRTM was found to be somewhat more sensitive to changes in CBF for both tracers, which is in agreement with work of Bullich et al. ⁷ Their findings differ from the results reported here, in that they reported a slightly higher sensitivity of SRTM to cortical and a lower sensitivity to global CBF changes. This discrepancy may be attributed to differences in simulation design, the different cortical composite region and differences in fitter settings (such as parameter boundary values) between the studies.^7,22 In the present study, parameter boundary values were optimised to prevent fit parameters from becoming undetermined (i.e. k₂ and BP_ND), in particular for lower levels of amyloid.

The finding that ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ is sensitive to changes in CBF is probably due to the fact that the assumptions underlying SRTM are violated for both tracers, i.e. the requirement that tracer kinetics are well described by a one-tissue compartment model in both the target and reference tissue. ²⁵ More specifically, for both tracers the target and reference tissue kinetics have been shown to be better described by a two-tissue compartment model.^3,21 In contrast, RLogan does not assume any specific number of target tissue compartments, robustly estimating ${\mathrm{DVR}}_{\mathrm{FLOW}}$ independent of the underlying compartmental separation. Furthermore, most outliers were present for interval compared to full TACs (outliers interval TACs: [¹⁸F]flutemetamol: 91.67%, [¹⁸F]florbetaben: 66.67%). These outliers were only observed for ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ and in particular for TACs corresponding to low amyloid levels.

With respect to ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ , effects of CBF changes were relatively small for both tracers (see Figures 3 and 4, maximum 5.1 and 6.1% for [¹⁸F]flutemetamol and [¹⁸F]florbetaben, respectively). However, SUVr’s sensitivity to global CBF changes increased for higher amyloid levels, as also seen for the other models. A comparable finding has been reported for [¹⁸F]florbetapir, where amnestic MCI and AD subjects showed higher CBF sensitivity compared to controls, for large CBF decreases (−40%). ¹⁸ The maximum change in ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ was within the maximum reported change for other amyloid tracers that assessed slightly more extreme flow variations.^5,18 In addition, the sensitivity of ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ to CBF changes was, as expected, dependent on acquisition start time, with later start times being more robust against CBF fluctuations, except for [¹⁸F]flutemetamol cortical changes. For [¹⁸F]flutemetamol, the direction of the relation between sensitivity to CBF changes and acquisition start time was, although less pronounced, comparable with [¹¹C]PiB findings. More specifically, for cortical CBF changes, an increased sensitivity was observed for later uptake-times, while this relationship was inverted for global changes. ⁵ These results suggest that imaging data acquired at expedited instead of delayed uptake times should be interpreted carefully, as a greater impact of CBF changes on SUVr is expected.

It is of interest to note that, in some cases, the confounding effects of perfusion changes were different between the two tracers. TACs were simulated using parameters derived from existing studies.^20,21 Although sample sizes of those studies were different, it is unlikely that this had an effect on the final results, as both datasets consisted of 50% AD patients and 50% healthy controls. The most likely explanation for these differences in results is the difference in tracer kinetics of the two tracers, which can be deduced from the rate constants used for the simulations (Table 1). As mentioned above, reference tissue kinetics of both tracers were best described by including a second compartment, thereby violating the assumptions of a reference tissue model to some degree. More specifically, individual k₃ and k₄ parameters, as well as their ratio, differed between both tracers. This, in combination with differences in target tissue kinetics, could explain the observed differences between both tracers, which were most pronounced for SRTM, a model that assumes single tissue kinetics in both target and reference tissues. In other words, the difference in sensitivity to perfusion changes between both tracers is most likely due some degree of violation of underlying model assumptions, which may be different for the two tracers.

For both tracers, strong relationships (R² ≥ 0.95) were observed between SRTM-derived R₁ changes and simulated changes in CBF. This finding was in agreement with the high correlation observed between SRTM-derived R₁ and ¹⁵ O-H₂O PET, reported by a combined ¹⁵ O–H₂O PET and [¹⁸F]florbetapir study. ¹⁶ More specifically, it suggests that for both [¹⁸F]flutemetamol and [¹⁸F]florbetaben, R₁ could be used as a proxy for measuring CBF.

Finally, the effect of noise (resembling regions of interest) on sensitivity to CBF changes for ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{RLOGAN}}$ , ${\mathrm{DVR}}_{\mathrm{FLOW}\_\mathrm{SRTM}}$ and ${\mathrm{SUVr}}_{\mathrm{FLOW}}\mathrm{ }$ was limited across amyloid levels. The effect of noise corresponding to a voxel level needs further validation, preferably using parametric imaging data.

As mentioned before, the main goal of this work was to assess the effect of global and regional CBF changes on regional amyloid quantification for [¹⁸F]flutemetamol and [¹⁸F]florbetaben studies. Overall, a maximum CBF-induced change in amyloid outcome measure of 13.4% (SRTM) was observed across methods. However, it must be noted that the evaluated range of CBF changes was more extreme than the change one would typically expect during the timespan of a clinical trial or longitudinal study from an ageing or early AD perspective. On the other hand, in late AD stages, the CBF changes may be larger and drugs may also exert effects on CBF.^12,13,15 Although the exact size of the effect is unknown for many drugs, the effect of regularly used substances such as caffeine and nicotine ranges from −22 to +25%.^28,29 Therefore, besides the effect of drugs or stimulants, effects of CBF changes on ${\mathrm{DVR}}_{\mathrm{FLOW}}$ and ${\mathrm{SUVr}}_{\mathrm{FLOW}}$ will likely be smaller than the effects reported in the present study.

The present results suggest that if large changes in amyloid load are to be expected, the possible confounding effects due to changes in perfusion are likely to be insignificant, obviating the need for a dynamic scanning protocol. However, in drug development studies, the drug may have unknown effects on perfusion, which must be well understood prior to implementation of simplified methods based on a static scanning protocol. This is especially relevant in secondary prevention trials, where measured changes in amyloid load may be small. In those cases, a dynamic or dual-time-window protocol should be used to assess whether changes in perfusion have occurred and to estimate their magnitude. Based upon this information, it can be decided whether a static protocol is sufficient or a dynamic protocol is needed to address a specific clinical or research question. Finally, it should be noted that a dynamic or dual-time-window scanning protocol also can provide a measurement of relative perfusion, which may be an additional relevant source of information. It should be noted that an exhaustive assessment of additional sources of bias was outside the scope of the present study. However, it is known from previous studies that SUVr typically shows a positive bias, while RLogan tends to underestimate the gold standard, and these intrinsic biases challenge the pooling of data from different methods.^5,7,20,21 Moreover, other factors such as image contrast and test-retest reliability may play a role when deciding on the method of choice for analysing data. ⁶ This implies that additional research, in particular test-retest and parametric imaging studies, is warranted.

Finally, given that the linear model (RLogan) outperformed the non-linear model (SRTM) with respect to robustness against CBF changes, one could conclude that SRTM is not the optimal model for describing the tracer’s kinetics. Alternatively, a linearised form, basis function implementation of SRTM (RPM) ³⁰ could be evaluated, which has shown improved performance, especially for noisy data, compared to the original implementation of SRTM for [¹⁸F]flutemetamol. ²⁰ The present study focussed on data corresponding to a region of interest level. As a next step, future studies could validate the application of parametric methods, such as RPM and SRTM2, ³¹ against the present results. The implementation of these methods (for example RPM’s basis function settings) would ideally be validated using imaging data.

Conclusion

RLogan was least affected by changes in cerebral blood flow and is therefore, at least within this context, the preferred method for regional amyloid quantification of [¹⁸F]flutemetamol and [¹⁸F]florbetaben. The same data could be analysed again using SRTM in order to determine the relative rate of delivery R₁, which showed good correlation with CBF. Finally, in most cases, effects of CBF changes on SUVr were relatively small, in particular for later uptake-times.

Supplemental Material

JCB918029 Supplemental Material - Supplemental material for Simulating the effect of cerebral blood flow changes on regional quantification of [¹⁸F]flutemetamol and [¹⁸F]florbetaben studies

Supplemental material, JCB918029 Supplemental Material for Simulating the effect of cerebral blood flow changes on regional quantification of [¹⁸F]flutemetamol and [¹⁸F]florbetaben studies by Fiona Heeman, Maqsood Yaqub, Isadora Lopes Alves, Kerstin Heurling, Santiago Bullich, Juan D Gispert, Ronald Boellaard, Adriaan A Lammertsma and on behalf of the AMYPAD Consortium in Journal of Cerebral Blood Flow & Metabolism