Study participants
The current study was based on data from the ADNI database. ADNI is a multicentre study launched in 2003 as a public-private partnership, led by Principal Investigator Michael W. Weiner, MD. ADNI’s primary goal is to test whether the combination of neuroimaging and biochemical biomarkers and clinical and neuropsychological assessments can be used for early detection and monitoring of AD dementia[15]. The ADNI study was approved by local Institutional Review Boards of all of the participating institutions, and informed written consent was provided by enrolled participants at each site. Full information regarding the ADNI inclusion/exclusion criteria is described elsewhere[16]. ADNI is a prospective cohort study that continues to recruit participants; this study was based on participants with available plasma p-tau181 data (data downloaded in June 2020).
The study population was classified into two diagnostic groups: cognitively unimpaired (CU) and cognitively impaired (CI) individuals. The CU classification was based on a CDR of 0; participants who had no cognitive dysfunction but reported subjective cognitive decline were analyzed together with CU, as per the National Institute of Aging-Alzheimer’s Association’s biological AD research framework[2]. The CI group consisted of individuals that were clinically defined as having MCI or AD dementia. MCI and AD dementia classification followed the criteria described elsewhere.[15,17] CSF p-tau181 and [18F]FDG PET data were matched with plasma p-tau181 data collected on the same ADNI study visit. Cross-sectional analysis was based on CSF and plasma measures of p-tau181 and PET measures of glucose metabolism ([18F]FDG) in a subset of n=823 participants (CU, n=262; CI, n=561 [MCI, n=426; AD, n=135]). Longitudinal analysis was based on a subset of participants with a baseline and longitudinal (up to 24 months) plasma p-tau181 and [18F]FDG PET assessment, which consisted of n=389 participants (CU, n=138; CI, n=251 [MCI, n=213; AD, n=38]). A description of the cross-sectional and longitudinal sample selections can be found in Additional file 1. The first available plasma p-tau181 measurement was used as the baseline time point for longitudinal analyses, as well as for age and diagnostic classification for cross-sectional and longitudinal analyses.
Plasma p-tau181 measurement
Blood samples were collected, shipped, and stored as described by the ADNI Biomarker Core Laboratory[18]. Plasma p-tau181 was analyzed with the Single Molecule Array (Simoa) technique, using a clinically validated in-house assay described previously[9]. Plasma p-tau181 was measured on Simoa HD-X instruments (Quanterix, Billerica, MA, USA) in April 2020 at the Clinical Neurochemistry Laboratory, University of Gothenburg, Mölndal, Sweden. Plasma p-tau181 data was collected over 47 analytical runs. Assay precision was assessed by measuring two different quality control samples at the start and end of each run, resulting in within-run and between-run coefficients of variation of 3.3%-11.6% and 6.4%-12.7%, respectively. Out of 3762 ADNI samples, four were removed due to inadequate volumes. The remaining 3758 all measured above the assay’s lower limit of detection (0.25 pg/ml), with only six below the lower limit of quantification (1.0 pg/ml). Plasma p-tau181 measurements were downloaded from the ADNI database (accessed 2020-06-20).
CSF p-tau181 measurement
CSF samples were collected by lumbar puncture, shipped, and stored as described by the ADNI Biomarker Core Laboratory[18]. CSF concentrations of p-tau181 were quantified using fully automated Elecsys immunoassays (Roche Diagnostics) at the ADNI Biomarker Laboratory at the University of Pennsylvania. The lower and upper technical limits for CSF p-tau181 were 8 and 120 pg/mL. Procedures have been described in detail previously[19,20].
MRI acquisition and processing
Pre-processed 3T MRI T1-weighted magnetization-prepared rapid acquisition gradient echo images were downloaded from the ADNI database; full information regarding ADNI acquisition and pre-processing protocols of MRI data can be found elsewhere [21,22]. Images underwent linear and non-linear registration to the ADNI template space, and all images were visually inspected to ensure proper alignment to the ADNI template.
PET acquisition and processing
Pre-processed [18F]FDG and [18F]Florbetapir PET images were downloaded from the ADNI database; full information regarding ADNI acquisition and pre-processing protocols of PET data can be found elsewhere[23]. Images underwent spatial normalization to the ADNI standardized space using the automatic registration of PET images to their corresponding T1-weighted image space as well as the linear and non-linear transformations from the T1-weighted image space to the ADNI template space. PET images were spatially smoothed to achieve a final resolution of 8 mm full width at half maximum (FWHM) and were visually inspected to ensure proper alignment to the ADNI template.
[18F]FDG and [18F]Florbetapir standardized uptake value ratio (SUVR) maps were generated using the pons and the full cerebellum as the reference region, respectively. For each participant, a global [18F]FDG SUVR value was estimated by averaging the SUVR from the angular gyrus, posterior cingulate, and inferior temporal cortices[24]. A global [18F]Florbetapir SUVR value was similarly estimated using the precuneus, prefrontal, orbitofrontal, parietal, temporal, anterior, and posterior cingulate cortices[24]. Amyloid-β (Aβ) positivity was determined for each participant by a global [18F]Florbetapir SUVR exceeding 1.11[25].
Statistical analyses
All nonimaging statistical analyses were performed using R v4.0.0. Voxelwise imaging statistical analyses were executed using the VoxelStats toolbox [26] in MATLAB version 9.4. Subjects were considered outliers if their baseline plasma p-tau181 value was three standard deviations above the population mean, and their data were excluded. Comparing demographic and clinical characteristics between diagnostic groups was done using χ2 test with continuity correction for categorical variables, Mann-Whitney U test for non-normal continuous variables, and one-way ANOVA for normal continuous variables. Correlations between plasma p-tau181 levels and demographic and clinical characteristics used Pearson’s correlation coefficient (r). All p values were two-tailed and p values <0.05 were considered significant.
Cross-sectional data were evaluated with correlations between CSF and plasma p-tau181 concentrations using Pearson’s correlation coefficient, with subjects stratified by diagnostic group and Aβ status. Voxelwise linear regression models tested the cross-sectional associations between [18F]FDG PET uptake and both CSF and plasma p-tau181 concentrations, adjusting for age and sex, in diagnostic groups (with and without Aβ status stratification).
Longitudinal analyses investigated the associations between baseline plasma p-tau181 levels and longitudinal metabolic decline. Annual rates of change were calculated both for global [18F]FDG SUVR and voxelwise for [18F]FDG images by subtracting the baseline value from the follow-up value and normalizing by time difference between time points, in years. Correlations and voxelwise linear regression models then tested the associations between annual rate of change in metabolic decline (using [18F]FDG SUVR and images, respectively) and baseline concentration of plasma p-tau181 and, adjusting for age and sex. Log-transformation of CSF and plasma p-tau181 measurements in pg/mL was used in all voxelwise analyses in order to reduce the skew of the distribution. Random field theory with a cluster threshold of p < 0.001 was used to correct voxelwise analyses for multiple comparisons[27].
Characteristic
|
Cross-sectional dataset (n = 823)
|
Longitudinal dataset (n = 389)
|
Diagnostic group
|
CU
|
CI
|
CU
|
CI
|
n
|
262
|
561
|
138
|
251
|
Age (median [IQR])
|
73.00 [68.52, 78.56]
|
72.81 [67.09, 77.60]
|
75.04 [69.93, 80.38]
|
71.47 [65.98, 77.30]
|
Males (n, %)†
|
120 (45.8)
|
317 (56.5)
|
75 (54.3)
|
138 (55.0)
|
Education (median [IQR])†
|
16.00 [15.00, 18.00]
|
16.00 [14.00, 18.00]
|
16.00 [16.00, 19.00]
|
16.00 [14.00, 18.00]
|
APOE ε4 carriers (n, %)†
|
75 (28.6)
|
296 (52.8)
|
35 (25.4)
|
127 (50.6)
|
MMSE (median [IQR])†
|
29.00 [29.00, 30.00]
|
28.00 [25.00, 29.00]
|
29.00 [28.25, 30.00]
|
28.00 [26.00, 29.00]
|
CDRSB (median [IQR])†
|
0.00 [0.00, 0.00]
|
1.50 [1.00, 3.00]
|
0.00 [0.00, 0.00]
|
1.50 [1.00, 2.50]
|
Plasma p-tau181 (median [IQR])§†
|
13.54 [9.32, 18.15]
|
18.08 [11.91, 24.39]
|
13.93 [9.72, 19.08]
|
15.51 [10.94, 24.80]
|
CSF p-tau181 (median [IQR])§†
|
19.67 [15.49, 26.61]
|
25.27 [18.04, 36.12]
|
-
|
-
|
Aβ+ (n, %)†
|
52 (19.8)
|
282 (50.3)
|
28 (20.3)
|
114 (45.4)
|
[18F]Florbetapir SUVR (median [IQR])†
|
0.96 [0.89, 1.07]
|
1.11 [0.94, 1.28]
|
0.96 [0.90, 1.06]
|
1.07 [0.94, 1.24]
|
[18F]FDG SUVR (mean (SD))†
|
1.80 (0.20)
|
1.67 (0.26)
|
1.78 (0.19)
|
1.72 (0.24)
|
Table 1: Demographic and clinical characteristics of the samples
† Statistically significant difference between groups (p < 0.05)
CU = cognitively unimpaired; CI = cognitively impaired; MMSE = Mini-Mental State Examination; CDR = Clinical Dementia Rating; Aβ = amyloid-β; SUVR = standardized uptake value ration; SD = standard deviation; IQR = interquartile range.
Table 1: The demographic and clinical characteristics for participants in the cross-sectional and longitudinal datasets are presented, stratified by cognitive status. Normal variables were summarized using mean and standard deviation, while non-normal variables were summarized using median and interquartile range. Statistical differences between the cognitively unimpaired and impaired groups were tested for both datasets, using χ2 test with continuity correction for categorical variables, Mann-Whitney U test for non-normal continuous variables, and one-way ANOVA for normal continuous variables.