Abstract
Background:
Plasma biomarkers of Alzheimer’s disease (AD) constitute a non-invasive tool for diagnosing and classifying subjects. They change even in preclinical stages, but it is necessary to understand their properties so they can be helpful in a clinical context.
Objective:
With this work we want to study the evolution of p-tau231 plasma levels in the preclinical stages of AD and its relationship with both cognitive and imaging parameters.
Methods:
We evaluated plasma phosphorylated (p)-tau231 levels in 146 cognitively unimpaired subjects in sequential visits. We performed a Linear Mixed-effects Model to analyze their rate of change. We also correlated their baseline levels with cognitive tests and structural and functional image values. ATN status was defined based on cerebrospinal fluid biomarkers.
Results:
Plasma p-tau231 showed a significant rate of change over time. It correlated negatively with memory tests only in amyloid-positive subjects. No significant correlations were found with any imaging measures.
Conclusions:
Increases in plasma p-tau231 can be detected at one-year intervals in cognitively healthy subjects. It could constitute a sensitive marker for detecting early signs of neuronal network impairment by amyloid.
INTRODUCTION
Alzheimer’s disease (AD) is a neurodegenerative process characterized by the accumulation of extracellular plaques of amyloid-β (Aβ) and intracellular neurofibrillary tangles of phosphorylated tau (p-tau) in the brain. These pathological changes happen decades before the first symptoms and are reflected in the cerebrospinal fluid (CSF) [1]. Thus, along the AD continuum, the concentration of CSF total tau (t-tau) and p-tau levels increase, and both Aβ42 and Aβ42/Aβ40 ratio levels decrease [2]. The preclinical stages, when we can already detect these changes, but the subject does not have symptoms yet, constitute a valuable therapeutic window to administer disease-modifying drugs before irreversible neurodegenerative damage develops [2, 3]. However, cost-effective, non-invasive, and population-level useful biomarkers are needed to identify suitable patients.
In the last few years, the field of AD has undergone a revolution with the development of high-sensitivity assays that identify plasma markers such as Aβ42/Aβ40, phosphorylated tau at Thr181 (p-tau181), Thr217 (p-tau217), Thr231 (p-tau231), neurofilament light, and glial fibrillary acid protein [4–9], opening up the possibility of avoiding other invasive, more expensive, and less available diagnostic methods such as lumbar puncture or positron emission tomography (PET) scans. Plasma markers have been shown to be highly accurate in distinguishing between AD pathology and other dementias [6, 10], and several studies suggest that the different p-tau species in plasma are sensitive markers of cerebral amyloid deposition [7].
With disease-modifying treatments on the verge of approval, it becomes necessary to understand these markers better, as they can help select suitable treatment candidates and monitor their response [11]. Moreover, they might provide us with accurate and reliable information in the preclinical stages of AD. In this regard, plasma markers offer valuable opportunities. Their accessibility makes it possible for sequential assessment, which may increase their diagnostic potential [12]. Thus, we can monitor changes in biomarker levels within a single subject rather than comparing them with the population. Thus, by understanding the evolution of plasma markers in the AD continuum, we could choose the earliest ones as a population screening tool in cognitively unimpaired (CU) subjects and also select the most appropriate markers to monitor the response to treatment.
However, these markers also present challenges, and many questions remain to be answered. It is unknown which marker is most useful at which stage of the disease, what time intervals are necessary to detect representative changes, and the relationship of these changes over time with neuropsychological or structural and functional imaging measures.
P-tau231 has emerged as one of the most promising plasma biomarkers. Previous studies in CSF and, more recently in plasma, suggest that p-tau231 alters very early, along with the first Aβ changes [7, 13]. Although its performance has been evaluated in patients with mild cognitive impairment and AD in cross-sectional studies, sequential changes of p-tau231 levels in the preclinical stages of AD need to be assessed. In this longitudinal study, we aim to provide information about the short-term evolution of p-tau231 in the preclinical stages of AD and its relationship with other parameters, such as neuropsychological or imaging markers.
MATERIALS AND METHODS
Participants
We performed this research with volunteers of the ‘Valdecilla Cohort for the study of memory and brain aging’ from the Memory Unit of the Marqués de Valdecilla University Hospital [14, 15]. The Valdecilla cohort is designed to help us better understand the preclinical phases of AD. It comprises CU volunteers who responded to an open call in the local media. Inclusion criteria are 1) age≥55 years; 2) signed consent for the collection and storage of biological samples and the performance of invasive techniques. Exclusion criteria are 1) cognitive impairment (determined by Clinical Dementia Rating (CDR)>0 [16]; 2) major psychiatric pathology; 3) significant systemic disease or sensory deprivation impairing the performance of the cognitive tests; 4) any contraindication to perform the complementary tests such as anticoagulation or claustrophobia.
At baseline, all participants were assessed with a comprehensive neuropsychological (NPS) battery conducted by neuropsychologists specialized in dementia, a brain magnetic resonance imaging (MRI), and an 18-fluorodeoxyglucose positron emission tomography (FDG-PET). DNA samples were extracted for genetic analyses, and CSF AD markers were determined (Aβ42, Aβ40, p-tau181, and t-tau). Annual follow-ups include plasma collection and a NPS assessment. One hundred and forty-six Caucasian participants were studied for baseline plasma p-tau231, one hundred and twenty-three had data from the first follow-up visit and sixteen for a second follow-up visit. All participants remained CU during the study.
This research has been conducted in accordance with the Declaration of Helsinki and has been approved by the ethics committee of the Hospital Universitario Marqués de Valdecilla. Title: Valdecilla Cohort for the study of memory and brain aging. Internal code: 2018.111. All subjects have given their signed consent to participate.
Plasma and CSF collection
Baseline plasma and CSF extractions were performed with the subjects fasting, on the same day, between 9 and 10 AM and a difference of less than 30 min between them. Our institution is part of the Alzheimer’s Association quality control program, so we follow international recommendations for plasma and CSF collection and storage [17, 18]. Aβ, p-tau, and t-tau values were determined using Fujirebio’s automated immunoassay analyzer Lumipulse G600 II [19] with the kits Lumipulse G β-Amyloid 1-40 (lot 4YX3085), Lumipulse G β-Amyloid 1-42 (lot 7ZX3084), Lumipulse G p-tau181 (lot 5DX3055) and Lumipulse G t-tau (lot 6BX3064).
Plasma samples were obtained following the standardized operating procedure described elsewhere [20]. In short, the blood is stored in 10 ml EDTA tubes and kept cold until processing within the next three hours. Then, it is centrifuged at 1800 g for 10 min. The supernatant is stored in polypropylene tubes in volumes of 500μl and frozen at -80°C until analyzed at the Clinical Neurochemistry Laboratory at the University of Gothenburg, Sweden. Plasma P-tau231 concentration was measured using an in-house ultrasensitive Single molecule array (Simoa) assay on an HD-X Analyzer (Quanterix, Billerica, MA, USA). Lower limit of quantification (2 pg/ml) and the rest of the analytical parameters can be consulted in previous studies [13].
CSF biomarkers
Participants were categorized according to ATN classification [21] based on their CSF biomarkers at baseline. As our cohort consists of CU subjects, instead of using the cut-off points established by the manufacturer in the kits mentioned above, they were established using an unbiased Gaussian mixture modelling [22]. We have dichotomized these variables and considered Aβ-positive (A+) when CSF Aβ42/40 ratio < 0.076, tau-positive (T+) when CSF p-tau181 > 73.2 pg/ml, and neurodegeneration-positive (N+) if CSF t-tau>543 pg/ml.
Magnetic resonance and PET imaging
All MRI scans were performed on the same 3T Philips Medical Systems MRI scanner with an 8-channel head coil. To determine the volume of the different structures, a sagittal T1-weighted MPRAGE sequence was used (170 slices, 1.2 mm voxel size, 9° of flip angle and shortest echo and repetition times) with subsequent processing as described in previous publications [14]. To segment the hippocampus, we used the automated FreeSurfer protocol, version 6.0 (http://surfer.nmr.mgh.harvard.edu). This protocol includes skull stripping, labeling of the volumes of each segmentation, and normalization of voxel intensities. Cortical and subcortical volume measurements were derived using the surface and subcortical segmentation pipelines [23]. Subcortical measurements were automatically derived from the subcortical processing flow. The quality check was performed using the ENIGMA Consortium quality control protocol (http://enigma.ini.usc.edu/).
The synthesis of 18F-FDG and the acquisition of PET images were carried out in our hospital’s Department of Nuclear Medicine. 18F was obtained using the 18O(p,n)18F nuclear reaction in a cyclotron PET Trace (General Electric Healthcare, Wisconsin, USA) and 18F-FDG scans were acquired in a Siemens Biograph LSO Pico 3D equipment (Siemens Healthcare Molecular Imaging, IL, USA). Subjects were administered intravenously with 3–4 MBq/kg 18F-FDG. Image acquisition consisted of a static image acquired 30-45 min after tracer injection. Images were reconstructed on a 128 x 128 matrix using the ordered subsets expectation maximization iterative method (six iterations and 16 subsets, zoom 2.5, filter FWHM 2 mm and sinogram trim factor 1.5). The axial slices were reoriented parallel to the frontal-occipital axis. The semiquantitative analysis of each study was performed using the Scenium ® program integrated into the Siemens Syngo.via Neurology software (Siemens CTI Molecular Imaging, Knoxville, TN, USA). The average standardized uptake value (SUV) ratio values in regions of interest (ROIs) from the Automated Anatomical Labelling Atlas [24] were used. We have taken the cerebellum as a reference.
APOE testing
Apolipoprotein E (APOE) genotype was determined by TaqMan SNP genotyping assays (Applied Biosystems, Foster City, CA, United States). Subjects carrying≥1 copy of the ɛ4 allele were considered ɛ4 + and the rest were considered ɛ4-.
Neuropsychological tests
Participants in our cohort performed a comprehensive neuropsychological battery designed to detect early signs of AD. It includes a global assessment through the CDR and the Mini-Mental State Examination (MMSE) [25], and several tests to evaluate the different cognitive domains. For memory assessment, we used the validated Spanish version of the Free and Cued Selective Reminding Test (FCSRT) [26], a test designed to measure verbal learning and memory using a list learning. We have considered different items such as the total free recall (TFR), total recall (TR), delayed free recall (DFR), and delayed total recall (DTR). For narrative memory, we used the logical memory (LM) subtest of the Weschler Memory Scale-III [27], which is based on the telling of two stories that are recalled in an immediate and delayed form (DR). In this test, we considered the total delayed units (LMDU). We also performed the Spanish version of the Face Name Associative Memory Exam (FNAME), a highly demanding face-name associative memory exam [28]. It consists of asking subjects to recall sixteen faces associated with a name and an occupation each. The first section includes initial recall of face-occupation and face-name pairs, followed by immediate recall of the name and occupation associated with each face. After 30 min, the subject is asked to recall as much of the information associated with each face as possible. The total value (TFNAME) is the number of all names and occupations recalled throughout the test and it is the measure we have considered for this study. The executive functions were assessed by means of the two parts of the Trail Making Test (TMT-A and TMT-B) [29] and the Digit-Symbol Substitution Test [26]. The Preclinical Alzheimer’s cognitive composite (PACC) [30] was also constructed with the results of the NPS assessment (we have used the PACC5 version, which considers the MMSE, the DR of the Logical Memory subtest, the Digit-Symbol test, the DFR and DTR of the FCSRT, and categorical fluency).
Statistical analysis
The assessment of normality was conducted using Kolmogorov-Smirnov test, and Shapiro-Wilk test. We used Student’s t-test for paired samples to analyze the differences between p-tau231 concentrations at baseline and follow-up and provided 95% confidence intervals (CI). We also used Student’s test-test for paired samples after stratifying subjects according to their ATN group followed by the Bonferroni correction to adjust p-values (p-values presented in the text are adjusted p-values). Additionally, we performed the analysis of variance (ANOVA) followed by Tukey’s test to examine the differences between the ATN groups. Given the small sample size of the second follow-up visit (N = 16), for this analysis we have only considered the baseline visit and the first follow-up. The effect size of the mean analysis was analyzed using Cohen’s d.
We performed a Linear Mixed-effects Model (LMM) to analyze the rate of change of the p-tau231 levels. LMM are useful when we have data with more than one source of random variability. The repeated measures taken over time allow to account for both within-person and across-person variability. The models performed considered the plasma p-tau231 concentrations as the dependent variable and number of days between plasma collections as the time variable. A random intercept and slope model was performed, both conditions allow everyone to have a different level of p-tau231 at baseline and different slope over time. Also, the ATN group was considered as an explanatory variable and age, sex, and APOE ɛ4 status as covariates. Furthermore, p-tau231 levels were skewed and transformed (tfp-tau231) by operating the Box-Cox transformation to normalize the distributions. We have used this technique to improve the accuracy of the predictions made by LMM since it stabilizes the variance and makes it easier to interpret the data.
As p-tau231 levels did not follow a normal distribution, we used Spearman’s Rho to study the correlation between the plasma biomarker level and CSF, neuropsychological and imaging markers. Thus, we correlated baseline p-tau231 values with baseline CSF biomarkers (Aβ40, Aβ42, Aβ42/Aβ40 ratio, p-tau181, and t-tau), MRI measures (right and left hippocampal volume), PET-FDG (regional metabolism in right and left hippocampus, right and left precuneus, and globally) and with the main outputs of FCSRT, FNAME, LM, and TMT test. We also correlated baseline p-tau231 with cognitive tests stratified by the ATN group. Afterward, we selected the significant results and adjusted for covariates using a univariate generalized linear model with cognitive measures as the dependent variable. Missing data have been handled by case removal.
Statistical analyses were performed using SPSS Statistics V.20.0 (IBM, NY, USA) and R Statistical Software (version 4.3.0; R Foundation for Statistical Computing, Vienna, Austria).
RESULTS
Sample description
The sample includes 146 cognitively unimpaired subjects (average MMSE score of 28.8±1.5 and average DTR of 14.7±1.95), with a mean age of 64.7±8.33 years, and female predominance (71.9%). Ninety-one participants were A-T-N- (74.7% females), 40 A + T-N- (70% females), and 15 A + T+N+(60% females). The mean age was significantly higher in the A + T-N- group than the A-T-N- group (66.8 and 63.7 years, respectively; p = 0.007) and in the A + T+N+group than the A-T-N- (70.1 and 63.7 years respectively; p = 0.0004). However, it was not significantly different between the A + T-N- and A + T+N+subjects (66.8 and 70.1 years; p = 0.06). The proportion of APOE ɛ4 carriers was 29.4%, and 37.7% of subjects had at least one positive AD marker in the CSF (Table 1). The average time between the collection of the baseline plasma samples and the first follow-up visit was 267.3±86.5 days, and between the first and second follow-ups was 343.7±62.1 days.
Baseline characteristics of the subjects
n, number of subjects; SD, standard deviation; MMSE, Mini-Mental State Examination; FCSRT, Free and Cued Selective Reminding Test; CSF, cerebrospinal fluid; Aβ, amyloid-beta; A, amyloid; T, tau; N, neurodegeneration.
Longitudinal analysis
Plasma p-tau231 was analyzed in 146 subjects at baseline, 123 at first follow-up, and 16 individuals at second determination. Correlations between baseline plasma p-tau231 and CSF biomarkers can be found in Supplementary Table 1. The changing of the p-tau231 concentrations in our sample was examined in two ways. First, we tested whether the average of p-tau231 was higher as the visits progressed (only at the baseline and first follow-up, due to the small sample size of the second follow-up); and second, we analyzed the rate of change of the biomarkers over time (days elapsed). The former analysis was performed to test the hypotheses that the sample experienced a higher average of p-tau231 in each collection, while the latter allows us to estimate the rate of change of p-tau231 concentrations per day elapsed.
Plasma p-tau231 differences between timepoints
The baseline p-tau231 values deviated from normality (Shapiro Wilk’s test p = 0.0003) and so did the levels at first follow-up (Shapiro Wilk’s test p = 0.006) and the change in p-tau231 from baseline to first follow-up (Shapiro Wilk’s test p = 0.0003). The mean baseline p-tau231 was 11.4 pg/ml±5.98 and 13.3 pg/ml±6.04 at the first follow-up. The difference between baseline and first follow-up levels was statistically significant (1.99 pg/ml; 95% CI 1.22 to 2.77; p < 0.0001; Cohen’s d = 0.36) (Fig. 1).

Longitudinal changes in p-tau231. The Y axis represents plasma p-tau231 concentration expressed in pg/ml and the X axis, the sequential visits. In the boxplot, the boxes show the interquartile range (the upper boundary is the Q3, and the lower boundary is the Q1). The line inside the box corresponds to the median of the sample and the whiskers represent the maximum (upper) and minimum (lower) values. 146 subjects were studied at baseline and 123 at first follow-up. The mean baseline was 11.40 pg/ml±5.98, and at first follow-up it was 13.28 pg/ml±6.04. The mean difference between baseline and first visit was significant (p = 0.0000013).
At baseline, the different ATN groups showed significantly different values of p-tau231 (p < 0.001). The difference between the A-T-N- and A + T-N- groups was not significant (1.88 pg/ml; p = 0.19); but it was significant between the A-T-N- and A + T+N+groups (6.27 pg/ml; p = 0.0003; Cohen’s d = 1.11) and between the A + T-N- and A + T+N+groups (4.38 pg/ml; p = 0.03; Cohen’s d = 0.73). The same happened at the first follow-up visit, with significant differences between the different ATN groups (p = 0.0001). Between the A-T-N- and A + T-N- groups there were no significant differences (1.74 pg/ml; p = 0.26) but there were among the A-T-N- and A + T+N+subjects (7.32 pg/ml; p < 0.0001; Cohen’s d = 1.47) and between A + T-N- and A + T+N+(5.58 pg/ml; p = 0.007; Cohen’s d = 0.77).
When we have assessed the differences between time points and stratified by ATN group, in the A-T-N-, the mean baseline p-tau231 was 10.1 pg/ml (±5.00), and at first follow-up visit it was 11.9 pg/ml (±4.51). In this group, the mean difference was 1.79 pg/ml (95% CI = 2.82 to 0.76; p = 0.03; Cohen’s d = 0.35). In the A + T-N- group, the mean baseline p-tau231 value was 11.8 pg/ml (±5.24) and 14.2 pg/ml (±6.4) at follow-up. The difference was also significant in this group, with a value of 2.36 pg/ml (95% CI = 3.56 to 1.16; p = 0.001; Cohen’s d = 0.37). Finally, in the A + T+N+group, the mean value of p-tau231 at baseline was 17.4 pg/ml (±8.0) and at follow-up 19.6 pg/ml (±8.88). The average difference in this group was not significant, with a value of 2.24 pg/ml (95% CI = 5.83 to 1.34; p = 0.58; Cohen’s d = 0.46) (Fig. 2). The p-tau231 values per time point and ATN group and the differences between them can be found in Supplementary Tables 2 and 3, respectively.

Changes in p-tau231 stratified by ATN group. The Y axis represents plasma p-tau231 concentration expressed in pg/ml and the X axis, the baseline and first visit. The boxes show the interquartile range (the upper boundary is the Q3, and the lower boundary is the Q1). The line inside the box corresponds to the median of the sample and the whiskers represent the maximum (upper) and minimum (lower) values. The dots outside the boxes indicate outliers. White boxes correspond to A-T-N- group, grey boxes to A + T-N- and the black ones to the A + T+N+group. An increase in p-tau231 values is observed between visits. The mean differences in the A-T-N- and A + T+N+groups were significant in both baseline and first visit (p = 0.0003 and p < 0.0001 respectively; asterisks with lines). The same happened between the A + T-N- and A + T+N+groups (p = 0.03 at baseline and p = 0.007 in the first visit; asterisks).
Longitudinal trajectory of Plasma p-tau231 concentrations over time
We analyzed the longitudinal change of tfp-tau231 concentrations considering the influence of variables such as age, sex, presence of APOE ɛ4 and ATN group. We estimated a tfp-tau231 increase of 0.0025 pg/ml/day (p < 0.0001). Sex did not significantly influence the tfp-tau231 basal levels, nor did the age (Estimate=-0.46 pg/ml, p = 0.22 and Estimate = 0.05 pg/ml, p = 0.08, respectively). APOE ɛ4 carriers showed a predicted increase in baseline tfp-tau231 levels of 0.8 pg/ml over non-carriers (p = 0.03), and the A + T + N + group presented an estimated baseline 1.62 pg/ml higher than the A-T-N- group (p = 0.01). However, the ATN and ApoE groups showed no interactions with the growth rate of tfp-tau231 over time (Fig. 3). Standardized estimates can be found in Supplementary Table 3.

Longitudinal trajectory of plasma p-tau231 concentrations over time (days) by ATN group, segmented by APOE, and adjusted by sex and age at baseline. In this linear mixed model plot the dots represent each individual value. The shaded area around the regression lines (continuous lines) represents the 95% confidence interval. The Y axis represents the Box-cox transformed plasma p-tau231 values. Dashed lines represent the individual trajectories of each subject. A-T-N- subjects are shown in red, A + T-N- in green, and A + T+N+in blue. In the X axis we show number of days between samples. We divided the plot between non-carriers of the ApoE4 allele (left side) and carriers of at least one copy (right side). There are no differences between slopes, but there are differences in the general level. The A + T+N+group showed to have an estimated baseline 1.62 units higher than the A-T-N- group (p = 0.014). The ATN groups showed no interactions with the different follow-up visits in the growth rate of p-tau231.
Correlations with neuropsychological and imaging measurements
We correlated p-tau231 plasma levels and their rate of change between baseline and the first follow-up visit with a wide range of neuropsychological tests and imaging measurements (Table 2).
Correlation between plasma p-tau231 and both neuropsychological and imaging measurements
NPS, neuropsychological; MRI, magnetic resonance imaging; FDG PET, fluorodeoxyglucose positron emission tomography; N, number of subjects; p, statistical significance; Rho, correlation coefficient; LMDU, total delayed units in Logical Memory test; DFR, delayed free recall in Free and Cued Selective Reminding test; DTR, Delayed Total Recall; TMT-A, Trail Making Test part A; TMT-B, Trail Making Test part B; TFNAME, total value of Face Name Associative Memory Exam; PACC5, preclinical Alzheimer cognitive composite.
P-tau231 levels at baseline showed no correlation with the total right (r=-0.01; p = 0.85) or left (r=-0.12; p = 0.23) hippocampal volume. Nor did the change in p-tau231 over time correlate significantly with left (r = 0.01; p = 0.94) or right (r = 0.01; p = 0.88) hippocampal volume. Considering the effect of age, the uptake of FDG-PET did not correlate significantly with baseline p-tau231 values in different regions of interest, such as right precuneus (r=-0.03; p = 0.68), left precuneus (r=-0.03; p = 0.76), right hippocampus (r=-0.09; p = 0.31), left hippocampus (r=-0.12; p = 0.17) and neither with global average uptake (r=-0.05; p = 0.59).
However, baseline p-tau231 levels showed a correlation with different neuropsychological tests, such as the DTR of the FCSRT (r=-0.25; p = 0.002) and with the TFNAME (r=-0.25; p = 0.002). These tests did not correlate with the change of p-tau231 over time (r=-0.026; p = 0.77 and r = 0.03; p = 0.76 respectively). Neither PACC5 was significantly correlated with baseline p-tau231 levels (r=-0.08; p = 0.35) nor with difference between samples (r=-0.03; p = 0.76).
To better understand our results, we stratified the baseline correlations of p-tau231 with neuropsychological tests by ATN groups. P-tau231 levels did not correlate significantly with TFNAME in the A-T-N- group (r=-0.03; p = 0.75), but it did in the amyloid-positive groups, in which the magnitude of the correlation increased progressively: A + T-N- (r=-0.37; p = 0.02) and A + T+N+(r=-0.67; p = 0.01) (Table 3). We further examined the association between baseline p-tau231 and TFNAME by ATN group accounting for age, sex, and APOE4 status by a generalized linear model. Using p-tau231 as a dependent variant, we found that in the A + T-N- group, the relationship with TFNAME remained significant (Estimate=-4.2; p = 0.04) after adjusting for age, sex, and APOE. This did not happen in the A + T+N+group (Estimate=-9.1; p = 0.16), though the number of individuals included in this analysis was limited (N = 14).
Correlation between baseline plasma p-tau231 and cognitive tests stratified by ATN group
p, statistical significance; Rho, correlation coefficient; LMDU, total delayed units in Logical Memory test; DFR, delayed free recall in Free and Cued Selective Reminding test; DTR, Delayed Total Recall; TMT-A, Trail Making Test part A; TMT-B, Trail Making Test part B; TFNAME, total value of Face Name Associative Memory Exam; PACC5, preclinical Alzheimer cognitive composite; A, amyloid; T, tau; N, neurodegeneration.
PACC5 did not show significant correlations with baseline p-tau231 in A-T-N- subjects (r=-0.1; p = 0.3), A + T-N- (r=-0.2; p = 0.18) nor in the A + T+N+group (r=-0.002; p = 0.9).
DISCUSSION
Novel plasma biomarkers for AD are a promising tool for diagnosing and classifying patients in hospital and primary care settings. Among them, p-tau231 rises in plasma very early and allows the detection of subjects at risk for AD [7, 13]. Previous studies in cognitively unimpaired individuals modelling the trajectory of p-tau231 showed that it significantly increased after amyloid burden passed the threshold of 26.4 centiloids on Aβ PET [7]. In the same study, p-tau217 levels also raised very early but required a higher threshold of amyloid pathology. On the other hand, p-tau181 is a very specific marker of AD and predicts cognitive impairment and hippocampal atrophy within one year, but it only increases significantly in plasma with much higher amyloid loads [7, 31]. The results of our study show that plasma p-tau231 levels progressively increase in CU subjects with a mean age of 64.8 years, in samples obtained with one year of difference on average. This change is appreciable even in A-T-N- subjects who would be considered as not in the AD continuum.
Although it could constitute an ultra-early sign of amyloid dysmetabolism, its significance is still unclear and does not correlate with any cognitive, functional, or structural trait. To clarify this aspect of our findings, longer-term studies determining p-tau231 in A-T-N- subjects are necessary. Evaluating the evolution of p-tau231 in subjects who convert from A- to A + during follow-up would be a relevant issue. We have been unable to assess this point because we do not perform serial CSF measurements in our cohort. A similar trend of increase of p-tau231 from baseline to the follow-up visit in A + T-N- and A + T+N+was found, though this was only significant in A + T-N-, which could be explained by the low numbers of the A + T+N+group (N = 12). This is a limitation of our study; therefore, we might be underpowered for some analysis.
The linear mixed model estimated that the growth per day was 0.0025 units of tfp-tau231, and we did not find an interaction with the ATN groups; therefore, the slopes of the increase we observed in the three ATN groups were not significantly different. This would point toward a linear increase in our population of subjects around 65 years of age on average independent of their ATN status. However, as mentioned above, some of our subgroups (especially the A + T+N+) had a very low number of participants, and we might be underpowered to detect interactions. In this sense, although subjects showed significantly older ages in the progressive ATN groups, age did not show an effect on baseline p-tau231 levels in our LMM, but the A + T+N+group did show a baseline increase of 1.62 units compared to the A-T-N- (and a difference of 6.27 pg/ml in the study of means), pointing out that plasma p-tau231 levels might be more related to AD pathology than to age. However, the standard deviation in the A + T+N+group was 8.01. This fact makes it difficult to interpret the changes in p-tau231 along the AD continuum, thus requiring longer-term studies and larger sample sizes to confirm which periods are necessary to identify changes in the marker greater than its intrinsic variability.
A relevant aspect of our cohort is that it encompasses the age spectrum in which AD is known to begin to increase in incidence [32]. Even though age was not significantly associated with p-tau231, we speculate that our population might be suitable to detect early p-tau231 levels increase, reflecting an early stage in the neurodegeneration process.
Longitudinal studies with post-mortem anatomopathological examination show that p-tau231 levels also increase in the late stages of AD and are higher than those of subjects with mild amyloid pathology [33]. This gradual correlation with the main pathological signatures of AD in clinicopathological studies suggests that the p-tau231 increase might be related to early AD changes. In this sense, p-tau231 has shown to be able to discriminate between Braak 0 and Braak I–II stages [34]. However, recent evidence indicates that, although p-tau217 increases progressively along the AD continuum, p-tau231 reaches a plateau and may not be useful for monitoring disease progression or long-term treatments efficacy [32]. P-tau181 appears to act similarly to p-tau217, as it increases progressively along the AD continuum, including the preclinical stages, and its levels correlate well with amyloid pathology [35, 36].
To characterize better the early preclinical phase of AD, we correlated p-tau231 levels with neuropsychological tests and measures of brain volume and metabolism. We found that the change from baseline to the follow-up did not correlate with any of these measures. However, baseline levels were associated with cognitive measurements. In contrast to what appears to happen with p-tau217 and p-tau181, which are associated with cortical atrophy in CU subjects [37, 38], we did not detect any significant correlation with structural measures such as hippocampal volume or glucose metabolism of different regions of interest for AD.
However, baseline levels of p-tau231 correlated with poorer cognitive outcomes, measured by sensitive neuropsychological memory tests such as the FNAME and FCSRT. Interestingly, when we stratified by ATN, this happened only from the A + T-N- group onwards and progressively increased along the AD continuum, and this correlation remained significant after adjusting for covariates. These findings are in line with another previous study that has shown a correlation between baseline p-tau231 levels and longitudinal changes in cognition, determined by the PACC at three years, only in Aβ+subjects [7]. Our interpretation of this finding is that we can only detect an association between tau pathology (plasma levels of p-tau231) and cognitive impairment in individuals with significant amyloid deposition (positive Aβ42/40 ratio in CSF).
For its part, in more recent studies, baseline levels of plasma p-tau217 have been shown to correlate with a worsening of memory (word list delayed recall) in a group of CU presenilin-1 E280A carriers [39] and increases in this biomarker also correlate with a worsening in modified PACC and MMSE over as long as six years in Aβ+CU subjects [35]. P-tau181 is not altered as early as p-tau217 and p-tau231 but has shown to be a good predictor of progression to AD and is longitudinally related to cognitive decline through parameters such as the MMSE, CDR and PACC [40].
Some limitations of our study should be pointed out. First, we only have CSF from the baseline visit, so we cannot confirm whether any of the subjects changed their ATN group during the follow-up time. Furthermore, although the changes in p-tau231 between timepoints are significant, the effect size is small, and this may difficult its application in clinical practice. On the other hand, even though our main findings are statistically sound, as we have mentioned above, the small sample size is a limitation for some of the stratified analyses so our results should be taken with caution and replications, ideally in larger and more ethnically heterogeneous samples, would be needed to confirm our conclusions.
Conclusions
Our results suggest that time intervals as short as one year could be meaningful in detecting significant changes in plasma p-tau231 values at a group level in CU participants, even in A-T-N- subjects. Sequential testing of plasma biomarkers could play an important role when selecting people for clinical trials and, in the near future, for initiating disease-modifying treatments in an efficient and non-invasive way. Still, long-term longitudinal studies are needed to give us a better understanding of its temporal evolution.
AUTHOR CONTRIBUTIONS
Francisco Martínez Dubarbie (Conceptualization; Formal analysis; Investigation; Writing – original draft); Sara López-García (Investigation; Writing – review & editing); Carmen Lage (Investigation; Writing – review & editing); Guglielmo Di Molfetta (Investigation; Writing – review & editing); Marta Fernández-Matarrubia (Investigation; Writing – review & editing); Ana Pozueta-Cantudo (Investigation; Writing – review & editing); María García-Martínez (Investigation; Writing – review & editing); Andrea Corrales-Pardo (Investigation; Writing – review & editing); María Bravo (Investigation; Writing – review & editing); Julio Jiménez-Bonilla (Investigation; Writing – review & editing); Remedios Quirce (Investigation; Writing – review & editing); Enrique Marco de Lucas (Investigation; Writing – review & editing); Marta Drake-Pérez (Investigation; Writing – review & editing); Diana Tordesillas (Investigation; Writing – review & editing); Marcos López-Hoyos (Investigation; Writing – review & editing); Juan Irure-Ventura (Investigation; Writing – review & editing); Elizabeth Valeriano-Lorenzo (Formal analysis; Software; Writing – review & editing); Kaj Blennow (Investigation; Resources; Writing – review & editing); Nicholas J. Ashton (Investigation; Resources; Writing – review & editing); Henrik Zetterberg (Investigation; Resources; Writing – review & editing); Eloy Rodríguez-Rodríguez (Conceptualization; Investigation; Methodology; Supervision; Writing – original draft); Pascual Sánchez-Juan (Conceptualization; Formal analysis; Investigation; Methodology; Supervision; Writing – original draft).
Footnotes
ACKNOWLEDGMENTS
We would like to thank all the Valdecilla Cohort volunteers for their participation.
FUNDING
KB is supported by the Swedish Research Council (#2017-00915 and #2022-00732), the Swedish Alzheimer Foundation (#AF-930351, #AF-939721 and #AF-968270), Hjärnfonden, Sweden (#FO2017-0243 and #ALZ2022-0006), the Swedish state under the agreement between the Swedish government and the County Councils, the ALF-agreement (#ALFGBG-715986 and #ALFGBG-965240), the European Union Joint Program for Neurodegenerative Disorders (JPND2019-466-236), the Alzheimer’s Association 2021 Zenith Award (ZEN-21-848495), and the Alzheimer’s Association 2022-2025 Grant (SG-23-1038904 QC).
CONFLICT OF INTEREST
KB has served as a consultant and at advisory boards for Acumen, ALZPath, BioArctic, Biogen, Eisai, Julius Clinical, Lilly, Novartis, Ono Pharma, Prothena, Roche Diagnostics, and Siemens Healthineers; has served at data monitoring committees for Julius Clinical and Novartis; has given lectures, produced educational materials and participated in educational programs for Biogen, Eisai and Roche Diagnostics; and is a co-founder of Brain Biomarker Solutions in Gothenburg AB (BBS), which is a part of the GU Ventures Incubator Program, outside the work presented in this paper.
Henrik Zetterberg is an Editorial Board Member of this journal but was not involved in the peer-review process of this article nor had access to any information regarding its peer-review.
All other authors have no conflict of interest to report.
DATA AVAILABILITY
The data supporting the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
