Abstract
Background:
Inflammation has been described as a key pathogenic event in Alzheimer’s disease (AD), downstream of amyloid and tau pathology. Preclinical and clinical data suggest that the cholinergic basal forebrain may moderate inflammatory response to different pathologies.
Objective:
To study the association of cholinergic basal forebrain volume and functional connectivity with measures of neuroinflammation in people from the AD spectrum.
Methods:
We studied 261 cases from the DELCODE cohort, including people with subjective cognitive decline, mild cognitive impairment, AD dementia, first degree relatives, and healthy controls. Using Bayesian ANCOVA, we tested associations of MRI indices of cholinergic basal forebrain volume and functional connectivity with cerebrospinal fluid (CSF) levels of sTREM2 as a marker of microglia activation, and serum levels of complement C3. Using Bayesian elastic net regression, we determined associations between basal forebrain measures and a large inflammation marker panel from CSF and serum.
Results:
We found anecdotal to moderate evidence in favor of the absence of an effect of basal forebrain volume and functional connectivity on CSF sTREM2 and serum C3 levels both in Aβ42/ptau-positive and negative cases. Bayesian elastic net regression identified several CSF and serum markers of inflammation that were associated with basal forebrain volume and functional connectivity. The effect sizes were moderate to small.
Conclusion:
Our data-driven analyses generate the hypothesis that cholinergic basal forebrain may be involved in the neuroinflammation response to Aβ42 and phospho-tau pathology in people from the AD spectrum. This hypothesis needs to be tested in independent samples.
INTRODUCTION
Alzheimer’s disease (AD) is characterized by a complex interplay of upstream pathologies, such as cerebral amyloid accumulation, and downstream effects on neuronal structural and functional integrity, partly mediated by neuroinflammatory response [1–3]. Interestingly, these effects seem to work in both directions. Thus, changes in neuronal functional activity and inflammation result from abnormal amyloid accumulation [4, 5]. At the same time, microglia derived neuroinflammatory factors can facilitate or even induce amyloid-β (Aβ) production [6, 7].
The cholinergic system is affected early in AD [8], and it is believed to moderate the inflammatory response to different pathologies [9]. Selective stimulation or suppression of basal forebrain cholinergic neuron activity in mice affected the peripheral expression of TNFα levels in response to endotoxemia [10]. This finding indicated that central action of acetylcholine on muscarinic acetylcholine receptors in the brain may suppress circulating TNFα and other proinflammatory cytokines. Similarly, the centrally acting acetylcholinesterase inhibitor galantamine suppressed an acid aspiration-induced acute respiratory syndrome in rabbits [11]. Consistent with the suspected association of cholinergic system integrity with inflammatory response, higher baseline expression of complement C3 in blood and a larger longitudinal increase of soluble TREM2 (sTREM2) levels in the cerebrospinal fluid (CSF) were associated with higher rates of atrophy of basal forebrain volume in cognitively healthy people with an increased CSF tau to amyloid ratio [12]. These effects were interpreted to indicate microglia involvement into clearance of amyloid and tau and resulting inflammatory responses. Taken together with the preclinical data on a modulatory role of cholinergic activity [10, 11], these data suggest that loss of structural integrity of the cholinergic basal forebrain may unleash the microglial inflammatory response to amyloid and tau accumu-lation.
Here, we leveraged structural and resting state functional MRI, CSF, and blood data from the DELCODE cohort [13] to extend these previous findings to the AD spectrum. First, in a hypothesis-driven approach we determined the effect of basal forebrain volume [14, 15] on levels of sTREM2 in the CSF and complement C3 in serum, controlling for Aβ and ptau status as well as clinical diagnosis across the AD spectrum. Secondly, extending previous evidence we determined associations between sTREM2 and complement factor C3 with resting state functional connectivity of the basal forebrain [16, 17]. Thirdly, we determined the association of a broad neuroinflammatory marker panel in CSF and serum with basal forebrain volume and functional connectivity. For this data-driven analysis we chose a Bayesian approach with an elastic net penalty [18] to directly estimate the plausibility of the presence or absence of an effect and its strength, and to avoid overfitting in the presence of large number of predictors. The data of these analyses can serve a better understanding of the possible role of cholinergic basal forebrain in the central and peripheral inflammatory response to Aβ and tau pathology in the AD spectrum.
MATERIAL AND METHODS
Participants
We used the baseline data, including 1,079 participants of the multicenter DELCODE study, conducted by the German Center for Neurodegenerative Diseases (DZNE) [13]. We excluded cases without available CSF inflammatory biomarkers, leaving 309 cases. After removing cases without an available MRI scan, we included 261 cases for further analysis. The sample group consisted of older healthy controls, first-degree relatives of a person with a documented diagnosis of AD dementia, and participants with AD, mild cognitive impairment (MCI), or subjective cognitive decline (SCD) (Table 1). DELCODE excluded participants with a current major depressive episode, past or present major psychiatric disorders, neurological diseases other than AD, or unstable medical condition [13]. SCD was defined as a persistent self-perceived cognitive decline in the absence of objective cognitive impairment as measured by the CERAD test battery, lasting at least for 6 months and being unrelated to an acute event [19]. The MCI patients met the core clinical criteria for MCI according to National Institute on Aging-Alzheimer’s Association (NIA-AA) workgroup guidelines [20]. The AD patients had a clinical diagnosis of probable AD dementia according to the NIA-AA workgroups guidelines [21]. The control participants had no objective cognitive impairment in cognitive tests, no history of neurological or psychiatric disease and did not report self-perceived cognitive decline. All participants or their representatives provided written informed consent. The study protocol was approved by the local institutional review boards and ethical committees of the participating centers. It was conducted in accord with the Helsinki Declaration of 1975.
Patient demographics
HC, healthy controls; SCD, subjective cognitive decline; MCI, mild cognitive impairment; AD, Alzheimer’s disease; MMSE, Mini-Mental State Examination. 1Bayes factor in favor of a dependence between diagnosis and sex (BF10) = 11.4; i.e., the dependency between sex and diagnosis is 11.4 times more likely than the independence of both factors. 2Bayes factor in favor of a group effect (BF10) = 7.2 * 106. 3Bayes factor in favor of a group effect (BF10) = 0.14; i.e., the absence of a group effect is 7 times more likely than the presence of such effect. 4Bayes factor in favor of a group effect (BF10) = 3.7*1039.
Biomaterial sampling and fluid biomarker acquisition
Biomaterial sampling procedures and biomarker measurements for DELCODE have been described elsewhere [13]. In brief, material was sampled by trained study personal following standard operating procedures (SOP) for collection and storage of samples. Both serum and CSF samples were aliquoted after centrifugation and stored at –80°C until analysis.
Image acquisition
The data were acquired from nine Siemens 3.0 Tesla MRI scanners (4 Verio, 1 Skyra, 3 TimTrio, and 1 Prisma system) using identical acquisition parameters and harmonized instructions. To ensure high image quality throughout the acquisition phase, all scans had to pass a semi-automated quality check during the study conduction, so that protocol deviations could be reported to the study sites, and the acquisition at the respective site could be adjusted. Functional MRI was based on a T2*-weighted echo-planar imaging (EPI) sequence using a 64×64 image matrix with 47 axial slices (thickness 3.5 mm, no gap) and interleaved acquisition. The field of view was 224×224×165 mm, isotropic voxel size of 3.5 mm, echo time 30 ms, repetition time 2,580 ms, flip angle 80°, and parallel imaging acceleration factor 2. The sequence took 7 min 54 s. High-resolution T1-weighted anatomical images were obtained using a sagittal magnetization-prepared rapid gradient echo (MPRAGE) sequence (field of view 256×256 mm, matrix size 256×256, isotropic voxel size 1 mm, echo time 4.37 ms, flip angle 7°, repetition time 2500 ms, number of slices 192, parallel imaging acceleration factor 2). The duration of the sequence was 5 min 8 s.
CSF AD biomarker assessment
Biomarker data were provided by the DELCODE data management. AD markers had been determined using commercially available kits according to vendor specifications: V-PLEX Aβ Peptide Panel 1 (6E10) Kit (K15200E) and V-PLEX Human Total Tau Kit (K151LAE) (Mesoscale Diagnostics LLC, Rockville, USA), and Innotest Phospho-Tau(181P) (81581; Fujirebio Germany GmbH, Hannover, Germany). For amyloid and tau isoforms, cut-off values for amyloid/tau isoforms were based on Gaussian mixture modeling of the DELCODE data independent of any group assignments using the R package flexmix (version 2.3–15) [22]: Amyloid ratio (A) Aβ42/Aβ40 0.08; tau pathology (T) by p-tau-181 73.65 pg/ml; neurodegeneration (N) by t-tau 510.9 pg/ml. We calculated the Aβ42/ptau ratio for each participant to classify cases into AD pathology positive and negative cases, respectively, using a a cutoff for the Aβ42/ptau ratio of 9.68.
CSF and serum inflammatory marker panel assessment
The inflammatory panel contained proteins with robust detection in both CSF and serum, related to inflammation and other mechanisms (ferritin, ApoE, sTrem2, sAXL, sTyro3, IL-6, IL-18, MCP-1, IP-10, macrophage migration inhibitory factor (MIF), YKL-40, CRP, complement factors C1q, C3, C3b, C4, factor B, factor H) as well as neurogranin and FABP-3 as non-tau related markers of neurodegeneration. Data of experimental biomarkers were provided by the DELCODE data management. In brief, all fluid biomarkers had been determined by use of commercially available enzyme-linked immunosorbent assays (ELISA) utilizing different assays formats and detection technologies optimized to the respective target’s dilution factors, detection range and multiplexing options, as described in [23]. Assay performance was controlled by use of internal control samples, with data determined in duplicates with a maximum coefficient of variance (CV) of 20%.
Imaging data processing
The T1-weighted anatomical images were initially coregistered to the mean functional images and subsequently preprocessed using the Computational Anatomy Toolbox (CAT12, v9.6/r7487 [24] for Statistical Parametric Mapping 12 (SPM12, v12.6/r1450, Wellcome Centre for Human Neuroimaging, London, UK). The images were segmented into grey matter (GM), white matter (WM), and CSF, followed by spatial normalization to the default CAT12 brain template in Montreal Neurological Institute (MNI) reference space using the Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra (DARTEL) algorithm. During this step, the images were resliced to an isotropic voxel size of 1.5 mm and modulated to adjust for expansion and shrinkage of the grey matter tissue. We then applied a mask containing the cholinergic nuclei of the basal forebrain [14] to derive the raw basal forebrain volumes. The raw volumes were divided by the total intracranial volume to adjust for head size. The resulting normalized basal forebrain volumes were entered in the statistical models.
Functional MRI data processing was carried out using Data Processing Assistant for Resting-State fMRI Advanced (DPARSFA 4.3) [25]. After the removal of the first ten images, we applied a series of steps including slice timing correction and realignment to eliminate the head motion. Scans with excessive head motion, i.e., a framewise displacement of greater than 0.5 mm for more than 30% of the scan duration, were excluded. We regressed out linear detrend, 24 head motion parameters, WM and CSF signal as nuisance regressors, to reduce the influence of noise. The functional images were band pass filtered to 0.1–0.01 Hz and spatially normalized to MNI space by applying the deformation fields estimated for the coregistered T1-weighted images beforehand. During warping, the functional images were resliced to an isotropic voxel resolution of 3 mm and finally smoothed with a 6 mm full-width-at-half-maximum (FWHM) Gaussian kernel. Based on the previous study [16], we defined the two functional subdivisions of the basal forebrain, anterior and posterior basal forebrain, as seeds for the functional connectivity analyses. Global functional connectivity was determined based on Pearson correlation coefficient of the seed region’s and all other GM voxel’s time series [26]. The resulting voxel-wise functional connectivity maps were Fisher z-transformed, thresholded at z≥0.3 (p < 0.001) to define the respective functional connectivity networks and averaged to produce the global functional connectivity values. Note that we included only positive correlation coefficients for computing the global connectivity, because positive and negative correlations otherwise may cancel each other out when averaging the correlation coefficients. Global connectivity values, one for each individual and functional division of the basal forebrain, were entered into the subsequent statistical analysis.
Statistical analysis
Statistical analyses were conducted in a Bayesian framework to allow estimation of model plausibility and determining effect sizes with credibility intervals. The analyses were conducted in two steps:
The first step was a hypothesis driven analysis of the effect of basal forebrain volumes on CSF TREM2 and serum complement C3, following previous evidence [12]. Consequently, we determined plausibility of regression models using CSF TREM2 or serum complement C3 as dependent variables, respectively, and basal forebrain volume as predictor with age, sex, education years, and dummy-coded diagnosis as covariates, separately for Aβ42/ptau-positive and Aβ42/ptau-negative cases. These index models were compared with the null models only containing the covariates. We used Bayes factor (BF) hypothesis testing to compare the index models containing the marker of interest against the null model. This approach allowed us to accept the best possible hypothesis for the data [27, 28]. We used Jeffreys’ Amazing Statistics Program (JASP Version 0.11), available at jasp-stats.org, to calculate the models. We report the Bayes Factor (BF10) quantifying evidence against the null hypotheses. To address potential issues with non-normally distributed residuals in the multiple regressions we applied Markov-Monte Carlo chain sampling to each analysis 1,000 times. We used the JASP default JZS prior. We decided a priori to follow up the size of the marker effect only for index models with at least moderate plausibility, i.e., BF10 > 3.
The second step was a data-driven candidate marker selection based on the 20 CSF and 20 serum inflammatory markers. Previous to the subsequent analyses all variables were scaled to zero mean and unit standard deviation.
Due to the large number of predictors and their expected collinearity we used regularized regression. Regularization aims to avoid overfitting of a model with a large number of predictors [29]. The penalty term used for regularization shrinks parameters towards zero and more so for small than large parameters [30], see Supplementary Figure 4 for an example. The Bayesian framework offers a wide variety of shrinkage priors, such as lasso, elastic net, or ridge [31]. We chose an elastic net penalty term that provides a reasonable compromise between both ridge and lasso penalties [18, 32] and allows selecting groups of correlated features from a high number of candidate markers [18]. The full Bayesian approach, as implemented in the library “bayesreg” in R (accessed through R Studio, version 1.1.463), allows for simultaneous estimation of both penalty parameters and model parameters [30, 31]. We used one half of the data as training data to estimate the parameter weights and the other half as the test data to determine the para-meter fit.
From the regularized regression models, we selected the most relevant parameters based on their credibility intervals. Thus, a parameter was kept in the model if its credibility interval did not contain zero and was excluded otherwise. The credibility interval provides the probability that the population parameters lie between the particular upper and lower bounds [33]. The credibility interval thresholds leading to most accurate feature selection with the Bayes elastic net have extensively been studied using simulations [30]. A credibility interval of 95% was found to be too wide, whereas a credibility interval of 50% was found to be most efficient in providing a set of variables with the highest fit of the true data [30]. Consequently, a credibility interval of 50% is the default in penalized regression in “bayesreg” or other libraries like “rstanarm”, and was used here as well.
RESULTS
Demographic and basal forebrain between group differences
As shown in Table 1, participants were highly likely to differ in age, sex, and, as expected, in MMSE scores, but not in years of education across diagnostic groups.
We found extreme evidence for an effect of diagnosis on basal forebrain volume after controlling for age and sex, with a Bayes factor of 1.3*1032 in favor of an effect of diagnosis. Evidence was in favor of no effect of diagnosis for anterior and posterior basal forebrain division functional connectivity (BF10 = 0.065 and 0.132, respectively). Further details on the post hoc comparisons across diagnostic groups can be found in the Supplementary Material, including Supplementary Table 1 and Supplementary Figures 1–3.
Basal forebrain volume and functional connectivity versus CSF TREM2 and serum C3
Using Bayesian ANCOVA in JASP, we compared the null model, including age, sex, education years, and diagnosis, with the index models, additionally including the volume or functional connectivity measures as predictors.
In the Aβ42/ptau-positive cases, the model predicting CSF-TREM2 by basal forebrain volume yielded a Bayes factor of 0.438, indicating that this model was 2.28 times less likely than the null model only containing the covariates. Similarly, in the Aβ42/ptau-negative cases, the model with basal forebrain volume was 4.0 times less likely than the null model (BF10 = 0.250).
In the Aβ42/ptau-positive cases, the model predicting serum complement C3 by basal forebrain volume was 2.0 times less likely than the null model (BF10 = 0.498). In the Aβ42/ptau-negative cases, the serum complement C3 model was 3.0 times less likely than the null model (BF10 = 0.331).
In 17 amyloid positive SCD cases and controls, similarly to the full sample analysis, we found anecdotal level of evidence in favor of the absence of an effect of basal forebrain volume on CSF sTREM2 levels (BF01 = 1.9) and serum complement C3 (BF01 = 1.7).
When using anterior and posterior basal forebrain functional connectivity as predictor variables, respectively, none of the models showed more than anecdotal evidence for CSF TREM2 or serum complement C3 as dependent variables; in fact, most models showed anecdotal to moderate evidence in favor of the null model (data not shown).
Since we found no evidence in favor of the index models, we conducted no further follow-up effect size estimates.
A complementary Bayesian implementation of robust regression, an approach, which is less sensitive to outliers than classical ANCOVA models, yielded no evidence for relevant effects of basal forebrain volume on CSF sTREM2 and serum complement C3 levels as well (see the Supplementary Material).
Basal forebrain volume and connectivity and central and peripheral inflammatory markers
Figure 1 provides an overview of the predictors that had relevant associations with basal forebrain volume or connectivity, separately for Aβ42/ptau positive and negative cases and central and peripheral inflammatory markers, yielding twelve different models. Relevant associations indicated that the 50% credibility interval of the respective parameter’s posterior distribution excluded zero. The posterior distributions of the most salient parameters are shown in Figs. 2 through 4, an overview over all parameters is shown in Supplementary Figure 5.

Predictors according to amyloid status and basal forebrain measures. For each of the twelve penalized regression models (Aβ42/ptau positive (_p) versus negative (_n), CSF versus serum, basal forebrain volume versus posterior basal forebrain connectivity versus anterior basal forebrain connectivity) we depicted those markers that showed a relevant effect, i.e., the 50% credibility interval of the marker’s posterior distribution excluded zero. Red markers had a negative, blue markers a positive association.

Posterior distribution of most salient CSF and serum parameters for basal forebrain volume. Posterior distributions of the parameters whose 50% CI excluded zero. The shaded area within the frequency distribution curve from the posteriors of each marker indicate the 50% CI, the blue vertical line indicates the mean. Figures show the CSF predictors (A) and the serum predictors (B) for Aβ42/ptau positive cases as well as the CSF predictors (C) and the serum predictors (D) for Aβ42/ptau negative cases.
For basal forebrain volume, in the Aβ42/ptau-positive cases the elastic net identified a negative effect of age, a diagnosis of MCI or a diagnosis of AD and a positive effect of a diagnosis of SCD (compared to a diagnosis of healthy control). CSF levels of YKL40, and Il18 were negatively, and CSF levels of AXL and MCP1 were positively associated with basal forebrain volume. The posterior distributions of the most salient parameters of the CSF model are shown in Fig. 2A. For serum levels relevant positive effects were found for CRP, C3, Factor-H, and MIF, and relevant negative effects were found for C4, Il6, and IL18. The corresponding posterior distributions of the most salient parameters can be found in Fig. 2B.
In the Aβ42/ptau-negative cases the elastic net identified ten CSF measures and eight serum measures that were positively or negatively associated with basal forebrain volume, among them CSF sTREM2, CSF FABP3, CSF YKL40, and CSF IL-6, as well as serum complement C3 and serum IL-6. Details of the posterior distributions of all salient parameters can be found in Fig. 2C and 2D, in Supplementary Figure 5C and 5D for all parameters.
Results for anterior basal forebrain connectivity in the Aβ42/ptau-positive cases were sparser. Only four CSF markers and only one serum marker had relevant associations with anterior basal forebrain connectivity, among them CSF ApoE protein and serum Il-18 (Fig. 3A, B), but also in the Aβ42/ptau-negative cases relevant associations with anterior basal forebrain connectivity were only few with four CSF and three serum markers (Fig. 3C, D).

Posterior distribution of CSF or serum parameters for anterior basal forebrain connectivity. The posterior distributions of the parameters whose 50% CI excluded zero predicting anterior basal forebrain connectivity. The shaded area within the frequency distribution curve from the posteriors of each marker indicate the 50% CI, the blue vertical line indicates the mean. Figures show the CSF predictors (A) and the serum predictors (B) for Aβ42/ptau positive cases as well as the CSF predictors (C) and the serum predictors (D) for Aβ42/ptau negative cases.
For posterior basal forebrain connectivity, in the Aβ42/ptau-positive cases we found relevant associations with six CSF and eight serum markers (Fig. 4A, B). In the Aβ42/ptau-negative cases there were no relevant associations of posterior basal forebrain connectivity with CSF markers, but only with nine serum markers (Fig. 4C, D).

Posterior distribution of CSF or serum parameters for posterior basal forebrain connectivity. The posterior distributions of the parameters whose 50% CI excluded zero predicting posterior basal forebrain connectivity. The shaded area within the frequency distribution curve from the posteriors of each marker indicate the 50% CI, the blue vertical line indicates the mean. Figures show the CSF predictors (A) and the serum predictors (B) for Aβ42/ptau positive cases as well as the CSF predictors (C) and the serum predictors (D) for Aβ42/ptau negative cases.
Estimates of the strengths of the effects
To estimate the added value of the CSF/serum makers which were identified in the previous data-driven analyses, we determined how much the strength of correlations between predictors and volume/functional connectivity measures increased from only the demographic predictors (including age, sex, education years and diagnosis) to the demographic predictors plus the most salient CSF/serum markers. The markers contributing are those reported in Fig. 1.
An overview of the results can be found in Fig. 5. Substantial increases of demographic correlations by the additional markers were only found for Aβ42/ptau-positive cases in three constellations: CSF markers and serum markers with basal forebrain volume as outcome, and serum markers with basal forebrain posterior connectivity as outcome (at least 75% probability that effects were higher with than without CSF markers, see Supplementary Figure 6 for an example). In addition, CSF markers substantially contributed to prediction of posterior basal forebrain connectivity; however, the correlation with demographics alone was only r < 0.20, rendering the overall effect of little relevance. In all other constellations, CSF and serum markers contributed only modest effects.

Correlations between volumes and demographic versus demographics plus salient CSF/serum parameters and probability of superiority. For each of the twelve penalized regression models (Aβ42/ptau positive versus negative, CSF versus serum, basal forebrain volume versus posterior basal forebrain connectivity versus anterior basal forebrain connectivity) we depicted the correlation for only the demographic predictors (red bars), the correlations for the demographic predictors plus the CSF/serum parameters (green bars) as well as the likelihood (in percent) that the correlation coefficient was higher with the CSF/serum parameters than with the demographic predictors alone (blue bars).
DISCUSSION
We investigated associations of inflammatory CSF and serum markers with volumes and functional connectivity of cholinergic basal forebrain in a group of older adults spanning the AD spectrum. Cases were stratified according to the CSF Aβ42/ptau ratio, indicating absence or presence of AD pathology.
Contrary to our a priori hypothesis which we had based on [12], we found evidence in favor of the absence of an effect of basal forebrain volumes and functional connectivity on CSF levels of sTREM2 and serum levels of complement C3, both in Aβ42/ptau ratio positive and negative cases. sTREM2 had been found involved in seeding of amyloid plaques [34], and increased with in vivo measures of microglia activation [35] in transgenic mouse models, suggesting a role in microglia-related phagocytosis of amyloid. Complement C3 as key factor in the complement cascade was found upregulated and associated with amyloid plaques containing dystrophic neurites in AD brains [36]. APP/PS1 mutation carrying mice with C3 knock-out demonstrated partly rescued cognitive phenotype and protection from neuron and synapse loss compared to APP/PS1 transgenic animals [37]. Peripheral and central cholinergic activity was found to be an endogenous factor to check innate immune response to different pathologies [9–11]. Our findings partly disagree with a previous study in preclinical AD cases [12]. This study had found that serum complement C3 at baseline was associated with a subsequently steeper decline of basal forebrain volume [12]. Based on sTREM2 and complement C3 alone we could not contribute convincing evidence from human data for the hypothesis that the cholinergic system modulated the microglia and complement-mediated response to Aβ and tau pathology. Of note, this conclusion is based on cross-sectional data only so that this does not preclude effects of baseline brain volume on longitudinal change in inflammatory markers. The higher power of longitudinal compared with cross-sectional analyses has been shown before empirically using genetic and neuroimaging markers as trait and state markers of disease, respectively [38]. In addition, the group comparisons of basal forebrain volume followed a pattern that would indicate that higher volume reflects more preserved structural integrity within the AD spectrum. However, for functional connectivity evidence was in favor of no group effect, suggesting compensatory effects in early stages of AD and questioning its usefulness as a direct marker of basal forebrain functional integrity in the AD spectrum.
Beyond sTREM2 and C3, a previous study had found no direct association of CSF levels of YKL40 und MCP-1 as markers of microglial activation with voxel-based measures of WM microstructural integrity from diffusion tensor imaging in cognitively healthy middle-aged to older people from the WRAP cohort [39]. This study reported, however, an interaction effect of MCP-1 with CSF Aβ42 on regionally restricted white matter regions in frontal and inferior temporal lobes [39]. A panel of anti- and pro-inflammatory serum markers was not associated with regional cortical and subcortical grey matter volumes in cognitively healthy older adults from the SNAC-K cohort [40]. Both studies did not include the basal forebrain as target region, and the study using a larger panel of markers [40] did not employ methods of variable selection.
In conclusion, preclinical data suggest an effect of cholinergic system integrity on the inflammatory response to danger-associated molecular patterns, but evidence from human studies, mostly focusing on a small number of selected inflammatory markers, is very limited so far. We would argue that a single inflammatory marker from CSF or serum or an unsystematic selection of markers from a larger panel may not be sufficient to capture the complex central and peripheral immune response to Aβ and ptau pathology. Therefore, we went beyond the previous studies and used a data driven approach to explore sets of inflammatory markers from CSF and serum that may be associated with cholinergic system changes in Aβ/ptau positive and negative cases. We used penalized regression to avoid overfitting. We chose elastic net regularization as this penalty can handle groups of correlated variables within a larger number of variables ([32], page 56), a feature which also characterizes the current sets of CSF and serum markers (as shown in Fig. 1 of [23]).
Subsequently, we will discuss the models where the inflammatory markers substantially contributed to the predictive value of the demographic characteristics (Fig. 5). In the other models the size of the effect appeared small so that discussing the contributing markers may be premature lacking further confirmation.
For the basal forebrain volume in the Aβ/ptau positive cases, likelihood was at least 75% that the serum and CSF markers contributed to the model fit compared with the demographic predictors alone. As a side note, this analysis illustrates the beauty of the Bayesian approach to enable a direct estimate for the relative likelihood of the data under the null and the alternative hypothesis. CSF and serum levels of IL-18 were negatively associated with basal forebrain volume. IL-18 serum levels had been found increased in prodromal and clinical stages of AD in previous studies [41, 42], while an earlier study with very small sample size showed no effects [43]. Functionally, IL-18 is released as proinflammatory cytokine downstream of NLRP3 inflammasome activation [44] and part of a cascade inducing amyloid plaque formation [45, 46]. The negative association of basal forebrain volume and levels of IL-18 in CSF and serum would be consistent with an inhibiting effect of cholinergic stimulation on microglial response to Aβ and ptau pathology [47]. We want to point out again, however, that given the exploratory nature of this analysis, this is a pure post-hoc explanation. In concert with IL-18, CSF levels of YKL40 were negatively associated with basal forebrain volume in the Aβ/ptau positive cases. A previous study used voxel-based analysis of grey matter volume changes over two years in cognitively healthy people. They found a significant positive association of YKL40 baseline CSF levels with a small inferior parietal cluster of change in grey matter volume over time [48]. YKL40 is expressed by reactive astrocytes [49] with elevated CSF levels in prodromal and early dementia stages of AD [50]. Our finding suggests lower levels of CSF YKL40 in presence of preserved cholinergic basal forebrain volume. In addition, in the Aβ/ptau positive cases, CSF levels of soluble Axl receptor tyrosine kinase (sAXL) and MCP-1 were positively associated with basal forebrain volume. MCP-1 is a chemokine attracting microglia to brain regions affected by pathology [51]. Co-expression of MCP-1 increased Aβ pathology in APP transgenic mice [52]. MCP-1 levels were increased in the CSF of AD patients in some [53], but not all [54, 55] previous studies. In MCI, MCP-1 levels in CSF were associated with faster cognitive decline [55]. Increase was mainly found at the dementia, but not the preclinical [53] or prodromal stages of AD [56]. Culture of microglia from the brains of AD patients compared with non-demented older people showed release of MCP-1 in response to Aβ peptide exposure which was, however, not different between AD and controls [57]. Together, these data suggest a potential disease accelerating role of MCP-1 starting at the prodromal and early clinical, but not the preclinical stages of AD. Therefore, the positive association of MCP-1 with basal forebrain volume in our Aβ/ptau positive cases is unexpected, since more preserved basal forebrain volume would indicate less affection by pathology.
Soluble AXL levels were associated with longitudinal decrease of CSF Aβ42 levels in cognitively healthy controls with normal Aβ42 levels at baseline, but not in healthy controls with abnormal Aβ42 at baseline and in AD patients [58]. These data may indicate that sAXL levels are increased in very early stages of AD pathogenesis, even preceding the onset of decreased Aβ42 levels in CSF and hence its parenchymal deposition. Therefore, our data may indicate sAXL increase with more intact cholinergic basal forebrain volume in the presence of abnormal Aβ42/ptau levels in CSF.
For posterior basal forebrain connectivity in the Aβ42/ptau positive cases, likelihood was 96% that the serum markers contributed to the model fit compared with the demographic predictors alone. These markers included CRP, Ferritin, sAXL, complement C4, IP10, MIF, IL-18, and ApoE protein in serum. Posterior basal forebrain connectivity was found decreased in a regionally very restricted manner with higher global amyloid levels from amyloid sensitive PET measures in cognitively normal people with subjective memory decline [17]. In the current DELCODE cohort we had shown that pattern of basal forebrain connectivity was spatially consistent between diagnostic groups with small reductions in Aβ42 positive cases with subjective cognitive decline, MCI and AD dementia compared with controls [59]. Here, we found a positive association of posterior basal forebrain connectivity with proinflammatory markers, such as sAXL, complement C4 or IL-18, and inflammation markers such as CRP or Ferritin in serum. IP-10 is a downstream marker of Th1 cell response that was found unaltered in CSF and serum of MCI and AD dementia patients compared with controls [60–62]. CSF levels of MIF were found only weakly associated with Aβ42, but more strongly with tau pathology [23], It is also considered an important opponent of glucocorticoid activity [63]. Its role in AD is unresolved with evidence for both pathogenetic as well as protective potential [64]. The basal forebrain posterior connectivity was positively associated with higher levels of proinflammatory and negatively with levels of potentially protective inflammatory markers in blood. This effect is consistent across markers. Functional connectivity may reflect a compensatory response to higher levels of inflammation. Of note, this is a post hoc interpretation that needs confirmation in an independent study.
Strengths of our study are the use of a comprehensive panel of CSF and serum markers of inflammatory response together with a systematic approach of feature selection using elastic net penalized regression. Bayesian analysis allowed us to directly assess levels of plausibility for the presence or the absence of an effect.
Limitations of our approach are the relatively small number of cases after stratification in Aβ42/ptau negative and positive cases that precluded us from analyzing effects within diagnostic subgroups. Therefore, we used clinical diagnosis as a co-variate in our analyses. Cross-sectional data can only provide limited insight into the complex interplay of brain atrophy and peripheral and central inflammatory marker up- and downregulation that may even follow a non-linear trajectory over time. Once available, we will use the longitudinal data of the DELCODE cohort to assess if basal forebrain volume is a predictor for the change in sTREM2 and complement C3 levels as well as of other neuroinflammatory markers over time across the AD spectrum. Of note, effect sizes for each single marker considered alone were small to moderate with direction of effects sometimes going in an unexpected direction. This may reflect the fact that each single marker is only a proxy for an underlying but unobservable mechanism of neuroinflammation. An alternative approach to make such presumed unobservable mechanisms explicit would be the use of structural equation modelling where neuroinflammation would be treated as a latent variable. This factor cannot directly be measured but each marker level in CSF or serum would be considered a realization of a measurement related to this latent factor. This may help to generate more consistency of findings across different cohorts, however, at the cost of allowing only indirect inference on possibly underlying specific pathways of neuroinflammation. Finally, different approaches are available to determine functional connectivity of brain regions from fMRI data. Here, we used a method that we had established in two independent cohorts before [16, 17], but we did not aim at exploring differences of results between different approaches of functional connectivity measurement.
In summary, with our cross-sectional data we could not confirm the assumption that sTREM2 and complement C3 release in response to Aβ42/ptau pathology would be moderated by volume or functional connectivity of the cholinergic basal forebrain. Data driven analysis across a large panel of CSF and serum markers of inflammatory response revealed several proinflammatory markers that were associated with basal forebrain volume and posterior connectivity in Aβ42/ptau positive cases. Due to the exploratory nature of this analyses, these data require independent confirmation. Together our findings provide preliminary evidence that cholinergic basal forebrain may be involved in the neuroinflammatory response to Aβ42 and ptau pathology in AD. The effect sizes, however, were small, suggesting that in vivo measures of neuroinflammatory markers in human studies may carry too much noise or may be of too little specificity to establish firm evidence for the role of cholinergic basal forebrain in the inflammatory response in AD.
Footnotes
ACKNOWLEDGMENTS
The DELCODE study (Study-ID: BN012) was supported and conducted by the German Center for Neurodegenerative Diseases (DZNE). The data samples were provided by the DELCODE study group. Details and participating sites can be found at
. The DELCODE study was supported in respect to the MR imaging by Max Delbrück Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin; Center for Cognitive Neuroscience Berlin (CCNB) at Freie Universität Berlin; Bernstein Center for Computational Neuroscience (BCCN), Berlin; Core Facility MR-Research in Neurosciences, University Medical Center Goettingen; Institute for Clinical Radiology, Ludwig Maximilian University, Munich; Institute of Diagnostic and Interventional Radiology, Pediatric Radiology and Neuroradiology, Rostock University Medical Center; and Magnetic Resonance research center, University Hospital Tuebingen.
