Longitudinal evaluation of surrogates of regional cerebral blood flow computed from dynamic amyloid PET imaging

Abstract

Surrogates of neuronal activity, typically measured by regional cerebral blood flow (rCBF) or glucose metabolism, can be estimated from dynamic amyloid PET imaging. Using data for 149 participants (345 visits) from the Baltimore Longitudinal Study of Aging, we assessed whether the average of early amyloid frames (EA) and R₁ computed from dynamic ¹¹C-Pittsburgh compound B (PiB) PET can serve as surrogates of rCBF computed from ¹⁵O-H₂O-PET. R₁ had the highest longitudinal test–retest reliability. Interquartile range (IQR) of cross-sectional Pearson correlations with rCBF was 0.60–0.72 for EA and 0.63–0.72 for R₁. Correlations between rates of change were lower (IQR 0.22–0.50 for EA, 0.25–0.55 for R₁). Values in the Alzheimer’s metabolic signature meta-ROI were negatively associated with age and exhibited longitudinal declines for each PET measure. In age-adjusted analyses, meta-ROI rCBF and R₁ were lower among amyloid+ individuals; EA and R₁ were lower among males. Regional PiB-based measures, in particular R₁, can be suitable surrogates of rCBF. Dynamic PiB-PET may obviate the need for a separate scan to measure neuronal activity, thereby reducing patient burden, radioactivity exposure, and cost.

Keywords

Amyloid cerebral blood flow longitudinal neurodegeneration neuronal activity

Introduction

Both structural and functional brain changes are important for the study of Alzheimer’s disease (AD). Fibrillar amyloid accumulation, a hallmark structural brain neuropathology of AD, begins about two decades prior to clinical diagnosis.¹ Ten percent of cognitively normal individuals aged 50 have amyloidosis, with prevalence increasing to over 40% by age 90.² ¹¹C-Pittsburgh compound B (PiB) was the first radiotracer used in human studies of amyloid PET.³ Several studies have demonstrated that ¹⁸F-fluorodeoxyglucose (FDG) PET measures of glucose metabolism, an indicator of brain function, are correlated with the mean of the early PiB time frames^4–8 as well as the R₁ parameter reflecting relative radiotracer delivery computed from dynamic PiB-PET.⁹ These results suggest that dynamic PiB imaging can be used to assess both the structural proteinopathies and functional properties of the brain. Similar results have been obtained with other amyloid radiotracers, including ¹⁸F-florbetapir (FBP)¹⁰ and ¹⁸F-florbetaben.^7,11 Additionally, regional cerebral blood flow (rCBF) as measured by dynamic ¹⁵O-H₂O-PET has been shown to correlate with PiB-K₁, the unidirectional influx rate of radiotracer from plasma to tissue in units of $ml / (ml \cdot min)$ ,^12,13 as well as with the unitless relative delivery parameter R₁ given by the ratio of tissue K₁ to reference region K₁.^14,15 Cerebral glucose metabolism and blood flow are reflective of neuronal, in particular synaptic activity.^16,17 Hence, ¹⁵O-H₂O and FDG-PET have been used as measures of functional neurodegeneration, which, as the aforementioned studies show, can also be assessed via dynamic amyloid imaging.

This dual-biomarker imaging of both fibrillar amyloid burden and neuronal activity via dynamic amyloid PET has important clinical implications: it may obviate the need for a separate PET scan to measure neuronal activity, thereby reducing patient burden and radioactivity exposure. The availability of information on neuronal activity in addition to amyloid burden can enable better differential diagnosis of dementias.¹⁸ Furthermore, understanding functional neurodegeneration is important for prognosis. For example, studies show that co-occurrence of amyloidosis and neurodegeneration is associated with steeper cognitive declines than the presence of either alone.^19–23 Therefore, characterization of neuronal activity, a proxy for neurodegeneration, alongside amyloidosis is key to a more complete understanding of cognitive decline. It is particularly important to assess longitudinal changes in neuronal activity to identify factors that can slow down or prevent cognitive decline and associated brain changes.

Previous reports that demonstrated the use of dynamic amyloid imaging for computing surrogates of neuronal activity have been limited to cross-sectional associations. It remains unclear whether surrogate measures computed from dynamic amyloid imaging are reliable longitudinally, and whether they exhibit longitudinal changes similar to those of rCBF. In this paper, we use rCBF as measured by static ¹⁵O-H₂O-PET and dynamic PiB-PET to address these gaps in knowledge by examining within-individual longitudinal associations and population-level longitudinal trends in rCBF, mean of early PiB time frames, and PiB-R₁. We hypothesized that both early PiB time frames and R₁ would be highly correlated with rCBF cross-sectionally and longitudinally, and that they would have similar associations with age, sex, and amyloid group.

Material and methods

Participants

We selected visits with both a ¹⁵O-H₂O and a dynamic ¹¹C-PiB-PET scan in the Baltimore Longitudinal Study of Aging (BLSA), and a structural magnetic resonance image (MRI) within six months of PET scanning. “Baseline” was defined as the earliest visit meeting this inclusion criterion. The resulting sample consisted of 159 individuals with a total of 355 visits acquired between June 2005 and September 2017. Ten individuals were used to define the optimal time frame for the early amyloid image. Subsequent statistical analyses were conducted using data for the remaining 149 individuals with 345 visits (Table 1).

Table 1.

Participant characteristics at baseline.

Characteristic	n = 149
Age, years, mean (SD)	76.7 (7.93)
Female, n (%)	75 (50.3%)
White, n (%)	114 (76.5%)
Education, years, mean (SD)	17.2 (2.37)
APOE ɛ4 + , n (%)^a	44 (29.7%)
Amyloid + , n (%)	47 (31.5%)
Clinical diagnosis, n (%)
Cognitively normal	143 (96.0%)
Mild cognitive impairment	4 (2.68%)
Dementia	2 (1.34%)
Follow-up duration, years, mean (SD)	2.50 (2.93)
Individuals with:
1 visit, n (%)	61 (40.9%)
2 visits	36 (24.2%)
3 visits	25 (16.8%)
4 visits	10 (6.71%)
5 visits	9 (6.04%)
6–8 visits	8 (5.37%)

APOE genotype was unavailable for one individual.

Research protocols were conducted in accordance with the United States federal policy for the protection of human research subjects contained in Title 45 Part 46 of the Code of Federal Regulations (45 CFR 46), approved by local Institutional Review Boards (IRB), and all participants gave written informed consent at each visit. BLSA PET study is governed by the IRB of the Johns Hopkins Medical Institutions (protocol number NA_00051793), and the BLSA study is overseen by the National Institute of Environmental Health Sciences (NIEHS) IRB. At enrollment into the neuroimaging substudy, which began in 1994, participants were free of CNS disease (dementia, stroke, bipolar illness, epilepsy), severe cardiac disease (myocardial infarction, coronary artery disease requiring angioplasty or coronary artery bypass surgery), severe pulmonary disease, and metastatic cancer.

Image acquisition and processing

Structural imaging

MRI acquisition is detailed in Supplement. Anatomical labels and regional brain volumes were obtained for each MRI using the MUSE labeling approach²⁴ with atlases that have been harmonized to account for differences in scanners and acquisition parameters.²⁵ This harmonization enabled the use of volumetric data acquired on the three different scanners (355 visits) together in statistical analyses. We used 118 MUSE regions of interest (ROI) covering the cerebral cortex, lobar white matter (WM), and subcortical gray matter (GM) structures (Supplementary Table). In our analyses, we focused on an AD metabolic signature composite region, the meta-ROI (comprised of bilateral angular, posterior cingulate, and inferior temporal gyri),²⁶ given that dynamic PiB-PET may be useful in assessing AD-related neurodegeneration in addition to gaining insight into a hallmark neuropathology of AD. We also evaluated a composite of all cortical GM regions.

Transmission scan

A transmission scan in two-dimensional mode using a ⁶⁸Ge rotating pin source was obtained prior to the emission scans and was used for attenuation correction during PET reconstruction.

¹⁵O-H₂O-PET imaging

After completion of the transmission scan, a static ¹⁵O-H₂O-PET scan was acquired in resting state on a GE Advance scanner in three-dimensional mode following an intravenous bolus injection of 370 MBq of ¹⁵O-H₂O. Image data were obtained for 60 s once the total radioactivity counts in the brain reached threshold levels. The image was reconstructed using filtered backprojection (image matrix = 128 × 128, 35 slices, pixel size = 2 × 2 mm, slice thickness = 4.25 mm). Each ¹⁵O-H₂O-PET scan was rigidly registered onto the corresponding MRI using normalized mutual information.²⁷ The MUSE label image was transformed from MRI to PET space using the resulting rigid registration matrix, and was resampled to the PET space voxel grid using nearest neighbor interpolation. The ¹⁵O-H₂O-PET images were normalized by cerebellar GM to obtain rCBF images consistent with PiB-R₁ (see below), which is based on a cerebellar GM reference region. The cerebellum has been used previously as a reference region for analysis of ¹⁵O-H₂O-PET.^14,28

PiB-PET imaging

Ten minutes after the ¹⁵O-H₂O-PET scan, a dynamic PiB-PET scan was obtained over 70 min (33 frames) on a GE Advance scanner in three-dimensional mode. Acquisition began immediately following an intravenous bolus injection of approximately 555 MBq of PiB, and time frames were binned according to the following protocol: 4 × 0.25, 8 × 0.5, 9 × 1, 2 × 3, and 10 × 5 min. The dynamic scan was reconstructed using filtered backprojection with a ramp filter, yielding a spatial resolution of approximately 4.5 mm FWHM at the center of the field of view (image matrix = 128 × 128, 35 slices, pixel size = 2 × 2 mm, slice thickness = 4.25 mm). Each frame was aligned to the mean of the first 2 min to correct for motion. Average of the first 20 min of PiB-PET scans was rigidly registered onto the corresponding MRI, and the MUSE label image was transformed from MRI to PET space. Radiotracer kinetics were modeled using a simplified reference tissue model with a cerebellar GM reference region, and distribution volume ratio (DVR) and relative delivery parameter (R₁) images were estimated using linear regression with spatial constraint.²⁹ R₁ is given by $K_{1} / K_{1}^{ref}$ , where K₁ is the unidirectional influx rate from plasma into tissue and $K_{1}^{ref}$ is the unidirectional influx rate into the reference region. We separated individuals into amyloid − and + groups based on a mean cortical PiB-DVR threshold of 1.062 derived from a Gaussian mixture model (Supplement). Early amyloid image (EA) was computed as the average of the time frames within the 0.75–2.5-min interval, normalized by cerebellar GM activity. This time interval was selected to maximize regional correlation with rCBF in a separate training data set (see Supplement).

Regional rCBF, PiB-EA, and PiB-R₁ mean values were computed in native PET space using transformed MUSE labels.

Data from the BLSA are available upon request from the BLSA website (https://www.blsa.nih.gov). All requests are reviewed by the BLSA Data Sharing Proposal Review Committee and are also subject to approval from the NIH IRB.

Statistical analyses

Intraclass correlation coefficients

To assess the reliability and stability of each PET measure, we computed the intraclass correlation coefficient (ICC) for each ROI. This entailed fitting a linear mixed effects model with a fixed and random intercept as well as a fixed slope to account for change over time. ICC was computed as the ratio of the variance of the random intercept to the total variance (i.e. sum of the random intercept and noise variances).

Within-ROI correlations

To determine whether EA or R₁ is a better surrogate of rCBF, we computed all pairwise Pearson correlations at baseline for each ROI. We repeated these correlations for rates of change, computed for 52 individuals with at least three PET measurements. For each ROI, rate of change in rCBF, denoted $\frac{Δ rCBF}{Δ t}$ , was defined as the coefficient of the age term in individual-specific linear regression models with rCBF as the outcome, and intercept and age as explanatory variables. We similarly computed the rates of change in EA and R₁.

Cross-sectional and longitudinal associations

We assessed cross-sectional and longitudinal associations of each meta-ROI PET measure with age, sex, and amyloid group using linear mixed effects models. We also explored models including APOE ɛ4 positivity as an additional independent variable. We first converted each PET measurement to a z-score by subtracting the mean and dividing by the standard deviation at baseline. For each PET measure, we fitted a linear mixed effects model with longitudinal PET measure as the dependent variable, and intercept, baseline age (centered at 75), sex, amyloid status, time from baseline, and interactions of baseline age, sex, and, amyloid status with time as fixed effects. We included a random intercept per individual. Statistical significance was defined as p < 0.05.

Power calculations

For each meta-ROI PET measure, we fitted a random slopes linear mixed effects model with age as the independent variable. Mean annualized change was computed as the fixed effect estimate for the age term, and standard deviation (SD) was computed as the standard deviation of the random slope. We computed the number of participants required per group to detect a groupwise difference equal to 50% of the estimated mean annualized rate of change given the estimated variance–covariance matrix based on a two-sided test with α = 0.05 and 80% statistical power, assuming that each individual is observed for four visits separated by two-year intervals.

Software

For statistical analyses, we used R version 3.4.3.³⁰ Linear mixed effects models were fitted using the nlme package.³¹ Plots were generated using ggplot.³² Power calculations were performed using the lmmpower function in the longpower package.³³

Results

Figure 1 shows rCBF, EA, and R₁ images for a participant with higher and a participant with lower levels of rCBF.

Figure 1.

Sample PET images for an individual with higher (top row) and lower (bottom row) levels of rCBF, superimposed on corresponding MRIs.

Intraclass correlation coefficients

PiB-R₁ ICCs were the highest in every ROI except for the right occipital pole, where EA had a higher ICC (Supplementary Figure 1). Median (interquartile range [IQR]) ICC across ROIs was 0.58 (0.50–0.64), 0.67 (0.57–0.73), and 0.81 (0.73–0.84) for rCBF, EA, and R₁, respectively. Meta-ROI ICCs fell in the highest quartile: 0.74, 0.81, and 0.86 for rCBF, EA, and R₁, respectively (Supplementary Table 1). Longitudinal trends in the meta-ROI for each PET measure are illustrated in Figure 2.

Figure 2.

Longitudinal meta-ROI PET measure trends.

Within-ROI correlations

Cross-sectionally, rCBF had higher correlations with R₁ compared to EA in 75 of the 118 ROIs (Supplementary Figure 2(a)). Median (IQR) Pearson correlations across ROIs was 0.65 (0.60–0.72), 0.67 (0.63–0.72), and 0.88 (0.85–0.92) for rCBF-EA, rCBF-R₁, and EA-R₁, respectively. Meta-ROI correlations fell in the highest quartile: 0.79, 0.76, and 0.93 for rCBF-EA, rCBF-R₁, and EA-R₁, respectively (Figure 3(a), Supplementary Table 2).

Figure 3.

(a) Cross-sectional and (b) longitudinal correlations among meta-ROI PET measures.

Rate of longitudinal change in rCBF was more strongly correlated with rate of change in R₁ than EA in 71 of the 118 ROIs (Supplementary Figure 2(b)). Median (IQR) Pearson correlations across ROIs was 0.39 (0.22–0.50), 0.42 (0.25–0.55), and 0.82 (0.73–0.89) for rCBF-EA, rCBF-R₁, and EA-R₁, respectively. Meta-ROI correlations fell in the highest quartile: 0.61, 0.65, and 0.94 for rCBF-EA, rCBF-R₁, and EA-R₁, respectively (Figure 3(b), Supplementary Table 2).

Cross-sectional and longitudinal associations

We did not observe any associations between the PET measures investigated and APOE ɛ4 positivity. Therefore, APOE was not included in the final linear mixed effects model. Associations with baseline age (β = −0.035, p = 0.0007 for rCBF; β = −0.021, p = 0.043 for EA; β = −0.022, p = 0.026 for R₁) and time (β = −0.062, p = 0.0012 for rCBF; β = −0.084, p < 0.0001 for EA; β = −0.061, p = 0.0001 for R₁) were statistically significant for each meta-ROI PET measure (Table 2).

Table 2.

Cross sectional and longitudinal associations of meta-ROI PET measures with age, sex, and amyloid group.

Independent/ Dependent variable	rCBF (z-score)	EA (z-score)	R₁ (z-score)
Intercept	0.21 (0.11) t = 1.86 p = 0.064	0.24 (0.12) t = 2.04 p = 0.043	0.37 (0.11) t =−3.36 p = 0.0009
Baseline age	−0.035 (0.010) t = −3.46 p = 0.0007	−0.021 (0.010) t = −2.04 p = 0.043	−0.022 (0.010) t = −2.25 p = 0.026
Male sex	−0.048 (0.15) t =−0.310 p = 0.76	−0.35 (0.16) t = −2.19 p = 0.030	−0.50 (0.15) t = −3.26 p = 0.0014
Amyloid + group	−0.39 (0.15) t = −2.56 p = 0.011	−0.12 (0.15) t = −0.779 p = 0.44	−0.31 (0.14) t = −2.22 p = 0.028
Time	−0.062 (0.019) t = −3.29 p = 0.0012	−0.084 (0.017) t = −4.94 p < 0.0001	−0.061 (0.015) t = −4.10 p = 0.0001
Baseline age × time	−0.0033 (0.0020) t = −1.60 p = 0.11	−0.0032 (0.0018) t = −1.74 p = 0.08	−0.0029 (0.0016) t = −1.81 p = 0.072
Sex × time	−0.018 (0.027) t = −0.648 p = 0.52	0.0091 (0.025) t = 0.368 p = 0.71	0.023 (0.021) t = 1.08 p = 0.28
Amyloid group × time	0.049 (0.032) t = 1.54 p = 0.13	0.0017 (0.029) t = 0.582 p = 0.56	−0.008 (0.025) t = −0.326 p = 0.74

Note: Statistically significant associations (p < 0.05) are in bold.

There were some differences across the different PET measures in relation to cross-sectional associations with sex and amyloid group. While the directions of the associations agreed, there were differences in the magnitude of the estimated regression coefficients and their statistical significance. rCBF and R₁ were lower in the amyloid + group (β = −0.39, p = 0.011 for rCBF; β = −0.31, p = 0.028 for R₁), whereas this association was not statistically significant for EA (β = −0.12, p = 0.44). EA and R₁ levels were lower among males (β = −0.35, p = 0.030 for EA; β = −0.50, p = 0.0014 for R₁), while this was not statistically significant for rCBF (β = −0.048, p = 0.76). We did not find statistically significant associations of baseline age, sex, or amyloid group with rate of longitudinal change in any of the PET measures. For a qualitative assessment of the differences between amyloid+ and amyloid−, we computed mean images per amyloid group (not adjusted for covariates) and subtracted them (amyloid− − amyloid+). The difference images (Figure 4) further demonstrated the similarity between rCBF and R₁ images, whereas the EA images exhibited subtler group differences.

Figure 4.

Amyloid group difference images (amyloid− – amyloid+) for rCBF, EA, and R₁ (not adjusted for covariates).

The lack of statistical significance for the amyloid status and EA association raised the possibility that EA may be contaminated by amyloid-binding PiB signal, nullifying our ability to detect lower neuronal activity among amyloid+ individuals. We assessed whether EA reflects amyloid burden by computing cross-sectional Pearson correlation between continuous PiB-DVR and EA at baseline among the 47 PiB+ individuals in our sample. This correlation (0.18, p = 0.22), while not statistically significant, was higher than the correlations of PiB-DVR with PiB-R₁ (−0.0065, p = 0.96) and rCBF (0.065, p = 0.66).

Power calculations

The mean (SD) annualized change in meta-ROI rCBF, EA, and R1 was −0.0031 (0.00070), −0.0025 (0.00035), and −0.0020 (0.0012), respectively. We found that to detect a groupwise difference equal to 50% of the estimated mean annualized change, 318, 243, and 232 participants would be needed per group for rCBF, EA, and R1, respectively.

Discussion

We evaluated whether two measures derived from dynamic PiB-PET scans, namely the average of the early frames (EA) and the relative delivery parameter R₁, can serve as surrogates of rCBF measured by static ¹⁵O-H₂O-PET. Compared to rCBF and EA, R₁ had the highest longitudinal test–retest reliability and required the smallest sample to detect groupwise differences in longitudinal rates of change in the meta-ROI, an AD metabolic signature composite region. IQR of cross-sectional Pearson correlations of PiB-based measures with rCBF was around 0.6–0.7, whereas longitudinal correlation IQR was around 0.2–0.5. The meta-ROI exhibited cross-sectional and longitudinal correlations in the highest quartiles. All meta-ROI PET measures had statistically significant associations cross-sectionally with age and showed statistically significant longitudinal declines. In age-adjusted analyses, meta-ROI rCBF and R₁ were lower among individuals with elevated amyloid levels (adjusting for sex), whereas EA and R₁ were lower among males (adjusting for amyloid status). The lack of statistical significance for lower EA among amyloid+ individuals can be explained by the weak correlation we observed between EA and amyloid burden, potentially compromising EA as a surrogate of rCBF. Given the higher ICC of R₁ and its lower correlation with amyloid burden compared to EA, R₁ is a more favorable surrogate measure for rCBF. Our findings confirm previous cross-sectional reports of early PiB time frames²⁸ and PiB-R₁^14,15 as surrogates of cerebral blood flow as measured by ¹⁵O-H₂O-PET. Our study adds to the existing body of knowledge by evaluating the longitudinal agreement between rCBF and these surrogate measures from dynamic PiB scans, demonstrating that PiB-based measures, in particular R₁, can be longitudinally reasonably reliable surrogates of rCBF in certain individual ROIs and the meta-ROI.

The recently published National Institute on Aging and Alzheimer’s Association (NIA-AA) research framework³⁴ recommends a biomarker-based classification system for studying the preclinical, MCI, and dementia stages of AD based on “β-amyloid deposition, pathologic tau, and neurodegeneration [AT(N)].” Our study shows that dynamic PiB-PET can be used to gain insight into two of these three categories, with DVR yielding a measure of amyloid burden and R₁ yielding a measure of functional neurodegeneration.

Certain kinetic properties, such as a high first-pass influx rate, of PiB and other amyloid radiotracers allow for the use of their early frames and radiotracer delivery parameters as correlates of rCBF.¹² These kinetic properties stem from these radiotracers’ moderate lipophilicity, which yields optimal passive transfer across the blood-brain barrier.³⁵ Our findings may be generalizable to other radiotracers with similar properties.

There is a tight coupling between rCBF and glucose metabolism, with rCBF increasing in response to increased glucose demand by functionally activated normal cerebral tissue,³⁶ which constitutes the biological basis for the observed correlations between ¹⁵O-H₂O and FDG-PET imaging.^37,38 Consistent with the association between CBF and glucose metabolism, our results agree with previous cross-sectional reports demonstrating that early amyloid time frames and PiB-R₁ can serve as surrogates of glucose metabolism.^4–6,8,9,28

Studies have shown that the early amyloid time frames are strongly correlated with the relative flow parameter R₁,⁸ and to a lesser extent, with the absolute flow parameter K₁^5,14 computed from amyloid PET. Direct comparisons of early amyloid frames and kinetic parameters reflecting radiotracer delivery with regard to their reflection of glucose metabolism or CBF have been limited. Ottoy et al. demonstrated good cross-sectional agreement among relative flow measures computed from 1-min ¹⁵O-H₂O, ¹⁸F-FDG, and ¹⁸F-FBP-PET, with the 0–6-min early amyloid frames yielding stronger correlations with ¹⁵O-H₂O compared to FBP-R₁ in most regions.²⁸ We found minimal differences between the EA and R₁ estimates with regard to their correlations with rCBF. However, the intraclass correlation coefficients accounting for overall longitudinal decline that we reported for PiB-R₁ were higher than the ICCs of rCBF computed from ¹⁵O-H₂O-PET and of PiB-EA. This suggests that PiB-R₁ is a more reliable measure, which is an important consideration, especially for longitudinal studies. The low-to-moderate ICCs of rCBF measures may in part explain the low correlations we observed for longitudinal rates of change between rCBF- and PiB-based measures.

We found similar associations with age and time for each of the three PET measures in the AD metabolic signature meta-ROI. This finding confirmed that each PET measure can capture age-related cross-sectional associations as well as longitudinal declines. Gietl et al. found negative associations between age and early PiB frames in the hippocampus, as well as negative age associations in frontal, temporal, and parietal GM regions among amyloid+ cognitively normal and mild cognitively impaired individuals.³⁹

Given that the correlations between rCBF and PiB-based measures were not perfect, we would expect to find some differences in statistical inference using different PET measurements. This was indeed the case for the analyses we conducted using the meta-ROI measures in relation to their associations with amyloid status and sex. We found that rCBF and R₁ were lower among amyloid+ individuals, but this effect was not statistically significant for EA. Given the weak positive correlation between PiB-EA and PiB-DVR among amyloid+ individuals, it is possible that amyloid-binding signal is present to an extent in the EA measurements, which could explain the lack of a statistically significant amyloid status association for EA. Similar to our study, the study by Gietl et al. did not find associations between early PiB frames and amyloid positivity.³⁹ Another study including individuals with dementia due to AD, as well as cognitively normal and mild cognitively impaired individuals, reported lower early PiB frame measures only in posterior cingulate and parietal cortices among amyloid+ individuals, whereas FDG-PET demonstrated wide-spread differences in cortical and subcortical GM.⁵ The sex effect we observed for the PiB-based measures, but not for rCBF, is harder to explain. This suggests possible differences between cerebral delivery or clearance of PiB and ¹⁵O-H₂O.

Our study has several limitations. Our previous investigations of partial volume correction on ¹⁵O-H₂O-PET scans yielded rCBF values with low signal-to-noise ratios due to increased noise, and therefore, we did not use partial volume corrected data. However, given that both our ¹⁵O-H₂O and PiB-PET scans were acquired on the same scanner, the same scanner point spread function applies, and thus, differences due to partial volume effects should be minimal. We could not evaluate differences in PET measures by clinical diagnosis given our small sample size. We did not investigate different methods for PiB kinetic parameter estimation, which could influence agreement of R₁ values with rCBF. We do not obtain arterial input functions given the longitudinal nature of our PET study, and as a result, had to rely on reference region approaches for estimating PET parameters and could not compute absolute delivery parameters K₁. A strength of our study was the ability to compare rates of change in rCBF to change in PiB-based surrogate measures, due to the extensive longitudinal PET data available in the BLSA.

Our conclusion that the PiB-R₁ kinetic parameter can constitute a surrogate measure of rCBF will be useful for other studies involving dynamic PiB-PET imaging. The longitudinal agreement we observed among these measures indicates that changes in neuronal activity over time can be assessed with longitudinal PiB-PET, without the acquisition of an additional PET scan. Our study demonstrates the value of a dynamic 0–70-min rather than a delayed PiB-PET acquisition. A better understanding of longitudinal changes in neuronal activity will help assess disease risk factors, identify protective mechanisms, and enable better predictors of cognitive progression.

In conclusion, the relative delivery parameter R₁ computed from dynamic PiB-PET can serve as a reasonable surrogate of relative regional cerebral blood flow computed from ¹⁵O-H₂O-PET. We found a moderate level of cross-sectional agreement between PiB-R₁ and rCBF, and a low-to-moderate level of agreement in longitudinal rates of change. Both PET measures exhibited similar age and amyloid status associations cross-sectionally and showed marked longitudinal declines. The higher longitudinal stability of PiB-R₁ compared to rCBF is an important property for longitudinal studies. Our results demonstrate that dynamic PiB-PET scans can provide a measure of neuronal activity in addition to amyloid burden. This dual-imaging property may obviate the need to acquire a ¹⁵O-H₂O or ¹⁸F-FDG-PET scan, thereby reducing patient burden, radioactivity exposure, and cost.

Supplemental Material

Supplemental material for Longitudinal evaluation of surrogates of regional cerebral blood flow computed from dynamic amyloid PET imaging

Supplemental Material for Longitudinal evaluation of surrogates of regional cerebral blood flow computed from dynamic amyloid PET imaging by Murat Bilgel, Lori Beason-Held, Yang An, Yun Zhou, Dean F Wong and Susan M Resnick in Journal of Cerebral Blood Flow & Metabolism

Footnotes

Acknowledgments

We thank the BLSA participants, BLSA staff, and the JHU PET facility faculty and staff for their assistance and dedication to these studies.

Authors’ contributions

MB conceived the study concept and design, analysis and interpretation of data, and drafting of manuscript; LB-H and YZ were involved in the revision of manuscript for intellectual content; YA took part in analysis and interpretation of data, and revision of manuscript for intellectual content; DFW was involved in the study design and data acquisition supervision, revision of manuscript for intellectual content; SMR took part in the study concept and design, oversight of data acquisition, interpretation of data, revision of manuscript for intellectual content.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The Federal Ministry of Education and Research (BMBF), Center for Stroke-research Berlin (CSB) and by grants from the Deutsche Forschungsgemeinschaft (DFG: VA 244/10-1) FOR 2325 NVI.

Declaration of conflicting interests

The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: This study was supported by the Intramural Research Program of the NIA, NIH. Dr. Wong receives contract funding via JHU from AVID/Lilly, Roche Neuroscience, Lundbeck and Five Eleven Pharma, and NIH.

Supplemental material

Supplemental material for this paper can be found at the journal website:

References

Villemagne

Burnham

Bourgeat

, et al. Amyloid β deposition, neurodegeneration, and cognitive decline in sporadic Alzheimer’s disease: a prospective cohort study. Lancet Neurol 2013; 12: 357–367.

Jansen

Ossenkoppele

Knol

, et al. Prevalence of cerebral amyloid pathology in persons without dementia. JAMA 2015; 313: 1924–1938.

Klunk

Engler

Nordberg

, et al. Imaging brain amyloid in Alzheimer’s disease with Pittsburgh Compound-B. Ann Neurol 2004; 55: 306–319.

Rostomian

Madison

Rabinovici

, et al. Early 11C-PIB frames and 18F-FDG PET measures are comparable: a study validated in a cohort of AD and FTLD patients. J Nucl Med 2011; 52: 173–179.

Forsberg

Engler

Blomquist

, et al. The use of PIB-PET as a dual pathological and functional biomarker in AD. Biochim Biophys Acta – Mol Basis Dis 2012; 1822: 380–385.

Liu

Zhang

, et al. Comparison of dual-biomarker PIB-PET and dual-tracer PET in AD diagnosis. Eur Radiol 2014; 24: 2800–2809.

Tiepolt

Hesse

Patt

, et al. Early [18F]florbetaben and [11C]PiB PET images are a surrogate biomarker of neuronal injury in Alzheimer’s disease. Eur J Nucl Med Mol Imaging 2016; 43: 1700–1709.

Rodriguez-Vieitez

Carter

Chiotis

, et al. Comparison of early-phase 11C-Deuterium-L-Deprenyl and 11C-Pittsburgh Compound B PET for assessing brain perfusion in Alzheimer disease. J Nucl Med 2016; 57: 1071–1077.

Meyer

Hellwig

Amtage

, et al. Dual-biomarker imaging of regional cerebral amyloid load and neuronal activity in dementia with PET and 11C-labeled Pittsburgh compound B. J Nucl Med 2011; 52: 393–400.

10.

Hsiao

Huang

Hsieh

, et al. Correlation of early-phase 18F-florbetapir (AV-45/Amyvid) PET images to FDG images: preliminary studies. Eur J Nucl Med Mol Imag 2012; 39: 613–620.

11.

Daerr

Brendel

Zach

, et al. Evaluation of early-phase [18F]-florbetaben PET acquisition in clinical routine cases. NeuroImage Clin 2017; 14: 77–86.

12.

Blomquist

Engler

Nordberg

, et al. Unidirectional influx and net accumulation of PIB. Open Neuroimag J 2008; 2: 114–125.

13.

Gjedde

Aanerud

Braendgaard

, et al. Blood-brain transfer of Pittsburgh compound B in humans. Front Aging Neurosci 2013; 5: 1–9.

14.

Chen

Rosario

Mowrey

, et al. Relative 11C-PiB delivery as a proxy of relative CBF: quantitative evaluation using single-session 15O-water and 11C-PiB PET. J Nucl Med 2015; 56: 1199–1205.

15.

Rodell

Aanerud

Braendgaard

, et al. Washout allometric reference method (WARM) for parametric analysis of [(11)C]PIB in human brains. Front Aging Neurosci 2013; 5: 45.

16.

Magistretti

Pellerin

. Astrocytes couple synaptic activity to glucose utilization in the brain. News Physiol Sci 1999; 14: 177–182.

17.

Attwell

Laughlin

. An energy budget for signaling in the grey matter of the brain. J Cereb Blood Flow Metab 2001; 21: 1133–1145.

18.

Silverman

Small

Chang

, et al. Positron emission tomography in evaluation of dementia: Regional brain metabolism and long-term outcome. JAMA 2001; 286: 2120–2127.

19.

Mormino

Betensky

Hedden

, et al. Synergistic effect of β-amyloid and neurodegeneration on cognitive decline in clinically normal individuals. JAMA Neurol 2014; 71: 1379–1385.

20.

Soldan

Pettigrew

Cai

, et al. Hypothetical preclinical Alzheimer disease groups and longitudinal cognitive change. JAMA Neurol 2016; 73: 698–705.

21.

Zhao

Tudorascu

Lopez

, et al. Amyloid β deposition and suspected non-Alzheimer pathophysiology and cognitive decline patterns for 12 years in oldest old participants without dementia. JAMA Neurol 2018; 75: 88–96.

22.

Burnham

Bourgeat

Doré

, et al. Clinical and cognitive trajectories in cognitively healthy elderly individuals with suspected non-Alzheimer’s disease pathophysiology (SNAP) or Alzheimer’s disease pathology: a longitudinal study. Lancet Neurol 2016; 15: 1044–1053.

23.

Wirth

Mormino

, et al. The effect of amyloid β on cognitive decline is modulated by neural integrity in cognitively normal elderly. Alzheimer’s Dement 2013; 9: 687–698.

24.

Doshi

Erus

, et al. MUSE: MUlti-atlas region Segmentation utilizing Ensembles of registration algorithms and parameters, and locally optimal atlas selection. Neuroimage 2016; 127: 186–195.

25.

Erus

Doshi

, et al. Longitudinally and inter-site consistent multi-atlas based parcellation of brain anatomy using harmonized atlases. Neuroimage 2018; 116: 71–78.

26.

Landau

Harvey

Madison

, et al. Associations between cognitive, functional, and FDG-PET measures of decline in AD and MCI. Neurobiol Aging 2011; 32: 1207–1218.

27.

Jenkinson

Bannister

Brady

, et al. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 2002; 17: 825–841.

28.

Ottoy

Miedema

De Puydt

, et al. Early frame 18F-AV45 and 18F-FDG-PET as proxies of CBF: Comparison to 15O-H2O PET data. Alzheimer’s Dement 2017; 13: P763–P764.

29.

Zhou

Resnick

, et al. Using a reference tissue model with spatial constraint to quantify [11C]Pittsburgh compound B PET for early diagnosis of Alzheimer’s disease. Neuroimage 2007; 36: 298–312.

30.

R Core Team. R: A language and environment for statistical computing 2017, https://www.R-project.org.

31.

Pinheiro JC, Bates DM, DebRoy S, et al. nlme: Linear and nonlinear mixed effects models. R package version 3.1-137, 2018, https://CRAN.R-project.org/package= nlme.

32.

Wickham

. ggplot2: Elegant graphics for data analysis, New York, NY: Springer-Verlag, 2016.

33.

Donohue MC and Edland SD. longpower: Power and sample size calculators for longitudinal data. R package version 1.0-16.1, 2016, https://cran.r-project.org/package=longpower.

34.

Jack

Bennett

Blennow

, et al. NIA-AA research framework: toward a biological definition of Alzheimer’s disease. Alzheimer’s Dement 2018; 14: 535–562.

35.

Pike

. PET radiotracers: crossing the blood-brain barrier and surviving metabolism. Trends Pharmacol Sci 2009; 30: 431–440.

36.

Paulson

Hasselbalch

Rostrup

, et al. Cerebral blood flow response to functional activation. J Cereb Blood Flow Metab 2010; 30: 2–14.

37.

Wong

CYO

Thie

Gaskill

, et al. A statistical investigation of normal regional intra-subject heterogeneity of brain metabolism and perfusion by F-18 FDG and O-15 H2O PET imaging. BMC Nucl Med 2006; 6: 2–7.

38.

Bentourkia

Bol

Ivanoiu

, et al. Comparison of regional cerebral blood flow and glucose metabolism in the normal brain: effect of aging. J Neurol Sci 2000; 181: 19–28.

39.

Gietl

Warnock

Riese

, et al. Regional cerebral blood flow estimated by early PiB uptake is reduced in mild cognitive impairment and associated with age in an amyloid-dependent manner. Neurobiol Aging 2015; 36: 1619–1628.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

1.04 MB

Longitudinal evaluation of surrogates of regional cerebral blood flow computed from dynamic amyloid PET imaging

Abstract

Keywords

Introduction

Material and methods

Participants

Image acquisition and processing

Structural imaging

Transmission scan

15O-H2O-PET imaging

PiB-PET imaging

Statistical analyses

Intraclass correlation coefficients

Within-ROI correlations

Cross-sectional and longitudinal associations

Power calculations

Software

Results

Intraclass correlation coefficients

Within-ROI correlations

Cross-sectional and longitudinal associations

Power calculations

Discussion

Supplemental Material

Supplemental material for Longitudinal evaluation of surrogates of regional cerebral blood flow computed from dynamic amyloid PET imaging

Footnotes

Acknowledgments

Authors’ contributions

Funding

Declaration of conflicting interests

Supplemental material

References

Supplementary Material

¹⁵O-H₂O-PET imaging