Abstract
INTRODUCTION
Alzheimer’s disease (AD) is a neurodegenerative disease that leads to progressive cognitive impairment and dementia. The pathology of AD, which starts decades before the appearance of clinical symptoms [1], is characterized by extracellular amyloid plaques containing amyloid-β (Aβ) and intracellular neurofibrillary tangles (NFTs) comprised of phosphorylated tau (P-tau) [2, 3]. The most widely accepted hypothesis regarding AD pathogenesis proposes that the deposition of Aβ leads to formation of NFTs, neuronal dysfunction, and dementia [4]. Although abundant evidence places Aβ and tau pathology at the center of AD pathogenesis, the mechanisms linking the two together and to eventual neuronal dysfunction and death are still unclear [4, 5]. As Aβ and tau pathology levels are relatively constant in the symptomatic stages of the disease, they are not considered suitable as progression markers [6–9]. Biomarkers reflecting other aspects of AD pathology are therefore needed. Genome-wide association studies have identified about 30 risk genes with a high proportion related to lipid metabolism, immune response, or both [10, 11]. These findings suggest that dysregulation of lipids and inflammatory proteins play an essential part in the pathogenesis of AD.
Lipids play various roles in the human body, both as structural components of cell membranes and in diverse biochemical processes, including membrane trafficking and cell signaling [12]. A dozen of major lipid classes are found within eukaryotic organisms, each comprising hundreds of individual molecular species [13]. Major classes, including glycerophospholipids, sphingolipids, fatty acids, and cholesterol [14], are abundantly expressed in the brain. These lipids are utilized in different compartments of glial cells and neurons [15]. In recent years, the role of lipid metabolism defects in AD pathogenesis has gained increased attention. Several molecular mechanisms have been identified, which connect membrane lipids to the generation and aggregation of Aβ. Pathological forms of Aβ proteins are formed by proteolytic cleavage of the transmembrane protein AβPP by β- and γ-secretases [16]. Lipid membrane structure and organization can affect the activity of these transmembrane enzymes, and thus AβPP processing and Aβ production [17]. Furthermore, secretases, AβPP, and its derivatives also appear to affect the activity of lipid metabolic enzymes and subcellular trafficking, thereby changing the membrane lipid composition. Lipids might, therefore, play a role in the initiation and progression of AD pathogenesis [16].
Studies have established neuroinflammation as a contributing factor in the pathogenesis and progression of AD and other neurodegenerative diseases [18, 19]. Aβ plaques induce an immune response by activation of microglia and astrocytes [20–22], which in turn is thought to play a role in the formation of NFTs, contributing to neuronal dysfunction and loss [23]. The glial proteins chitinase-3-like-1 protein (YKL-40), calcium-binding protein S100B, and glial fibrillary acidic protein (GFAP) have been associated with AD pathology [24]. All proteins are expressed primarily (YKL-40 and S100B) [25, 26], or exclusively (GFAP) [27] in astrocytes within the central nervous system (CNS). YKL-40, a chitin-binding glycoprotein [28], has been reported to be a promising candidate biomarker of glial activation in AD. Previous studies have detected positive relationships between YKL-40 and the neurodegeneration markers tau [29–32] and neurofilament light (NFL) [33] in cerebrospinal fluid (CSF), demonstrating an association between glial activation and neurodegeneration [34]. NFL is mainly located in myelinated axons, with recent studies indicating a potential for this protein as both a diagnostic and a progression marker in AD and other neurodegenerative diseases [35, 36]. S100B, a calcium-binding protein, exerts both intracellular and extracellular functions and has been found to be upregulated in AD tissues [37, 38]. GFAP is an intermediate filament protein and a marker for astrocyte activation, which has both been associated with amyloid plaque load and the number of NFTs [39–41].
In recent years, a paradigm shift has occurred from clinical to biological definition of AD based on
MATERIALS AND METHODS
Subjects
Subjects from The Icelandic MCI study cohort (

Flow diagram of sample selection and lipid profiling techniques applied to the analysis of CSF samples.
The study has been approved by the National Research Ethics Committee of Iceland (VSN-14-028) and all subjects signed an informed consent. The study was conducted in accordance with the Helsinki Declaration latest revision of 2013.
CSF collection and analysis
Collection of CSF was done via lumbar puncture with a 22-gauge spinal needle at the L3/4 or L4/5 interspace. Samples, uncentrifuged, were frozen in 2 ml polypropylene tubes and stored at –80°C. Levels of all proteins were determined using commercially available sandwich enzyme-linked immuno-sorbent assays (ELISAs) and performed according to manufacturer‘s instructions. Levels of T-tau (IBL International, Hamburg, Germany), P-tau181 (INNOTEST, Gent, Belgium), and Aβ42 (IBL International, Hamburg, Germany), were measured in the ISO 15189 accredited medical laboratory MVZ Labor P.D. Dr. Volkmann und Kollegen GbR (Karls-ruhe, Germany). Levels of NFL (Uman Diagnostics, Umeå, Sweden), YKL-40 (Quantikine ELISA Human Chitinase-3–like 1; R&D systems, MN, USA), S100B (BioVendor GmbH, Heidelberg, Germany), and GFAP (BioVendor GmbH, Heidelberg, Germany) were measured in a laboratory at the University of Iceland. All assays had mean Intra-assay CV < 10% and Inter-assay CV < 15%.
Subject grouping based on CSF measures
Each subject was classified based on CSF T-tau and Aβ42 values, independently of clinical diagnosis. A cut-off of 0.52 for T-tau/Aβ42 ratio was selected based on results from a large memory clinic cohort study [49]. T-tau/Aβ42 ratio > 0.52 was defined as a signature CSF AD profile. The CSF AD profile group had a total of 34 subjects (20 with a clinical diagnosis of AD dementia, 10 with MCI, three with SCI and one with Lewy body dementia) while the non-AD profile had a total of 26 subjects (13 with MCI, 10 with SCI, two with LBD, and one with Parkinson’s disease). The same ratio cut-off point was also used as a part of the clinical diagnosis of AD, explaining full concordance with the CSF AD profile.
Neuropsychological tests
A detailed neuropsychological assessment, for the evaluation of different cognitive domains, was performed by licensed psychologists under the supervision of a clinical neuropsychologist. A significant impairment in episodic memory is commonly the earliest clinical symptom of AD [50], and therefore of specific interest here. Two tests were used for the evaluation of verbal episodic memory, The Rey Auditory Verbal Learning Test (RAVLT) [51], and a Story test based on the Logical Memory test of the Wechsler Memory Scale-Revised [52]. RAVLT consists of 15 nouns presented across five consecutive trials, with each trial followed by a free-recall test (immediate recall). A score for RAVLT immediate recall was calculated by summing up the number of words recalled from trials 1 through 5 (0 to 75 points). After a 30 min delay, subjects were required to recall the words without being reread the list (delayed recall). A point was given for each correct word (0 to 15 points). The second test was composed of an orally presented story, which included 25 ideas. Right after the presentation, the subject was asked to repeat what they remembered without being given any clues. After a 30 min delay, there was another recall without the story being reread to the subject. For both immediate and delayed recall, a point was given for each idea (0 to 25 points).
Sample preparation and scanning UPLC-MS analysis
CSF sample extraction was based on the method used in Bird et al. [53]. Briefly, C12 ceramide and C17 sphingomyelin (SM d18:1/17:0) were purchased from Avanti Polar Lipids (Alabaster, AL, USA) and added to 30μl of CSF as internal standards prior to lipid extraction. Dried lipid extracts were resuspended in 300μL of ACN/IPA/H2O (65:30:5 v/v/v) and stored at –80°C prior to analyses. C18 SM (d18:1/18:0) and C18 ceramide (d18:1/18:0) were run alongside samples for reference as both SMs and ceramides have consistently been associated with AD [54]. Ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS) lipidomic analysis was carried out as described in Kotronoulas et al. [55]. The analytical instrumentation used was an ACQUITY UPLC system (UPLC ACQUITY, Waters Corporation, Milford, MA) coupled to a traveling-wave ion mobility (IMS) qTOF mass spectrometer (Synapt G2 HDMS, Waters Corporation, Manchester, UK). The chromatographic gradient separation was performed on an ACQUITY CSH C18 column (2.1 mm× 150 mm, 1.7μm particle size, Waters Corporation) at 60°C (Supplementary Figure 1). Mobile phase A was H2O/2-isopropanol (80:20) and mobile phase B 2-isopropanol/ACN/H2O (90:9.1), both with 0.05% of formic acid and 5 mM ammonium formate. Injection volume was 10μL, flow rate was 0.4 mL/min, and the run time was 17 min. The following gradient pattern (solvent B) was used: 0 min, 40% B; 1 min, 40% B; 3 min, 60% B; 10 min, 100% B; 13.5 min, 100% B; 14 min, 40% B; 17 min, 40 % B. Both positive (+) and negative (-) electrospray ionization (ESI) modes were acquired. The capillary and cone voltage were 2.5 kV and 30 V, respectively. The source and desolvation temperature were 120 and 500°C and the desolvation gas flow was 800 L/h. During High Definition MSE (HDMSE) experiments, the collision energy in the trap cell was off, and in the transfer cell, it ranged from 20 to 30 eV for the positive mode and from 25 to 40 eV for the negative mode. The resulting MS data were analyzed further within the MS-Dial program.
Spectral processing
Waters RAW files were converted to ABF format using the Reifycs Abf Converter tool. All patient ABF files and pooled samples from either positive or negative controls were loaded into the MS-Dial application [56]. MS-Dial was used for sample alignment and peak detection. Selected mass-to-charge ratio (m/z) peaks, hereafter referred to as m/z features, were manually curated to minimize the effect of sample drift. All m/z features were normalized to the internal standards. M/z feature annotation was performed with MS-Dial and compared to MS/MS spectra in the LipidBlast Library [57]. A total of 1013 m/z features were detected based on their peak mass. M/z features with more than 5% missing values were excluded from the analysis (
Statistical analysis
All CSF measures (detected m/z features and proteins) were log2-transformed to fit a Gaussian distribution. After autoscaling, each CSF measure had an average and a standard deviation of 0 and 1, respectively. Mann-Whitney U non-parametric tests were performed to compare levels of different variables between CSF profile groups. The sample (
RESULTS
Table 1 presents the demographic, pathophysiological, and cognitive characteristics of subjects, both within the whole sample and divided by CSF profile (AD and non-AD). No statistical differences (
Subject demographics, CSF protein levels and cognitive scores by CSF profile
AD, Alzheimer’s disease; CSF, Cerebrospinal fluid; LBD, Lewy body dementia; MCI, mild cognitive impairment; MMSE, Mini-Mental State Examination; N/A, Not applicable; PD, Parkinson‘s disease; RAVLT, Rey Auditory Verbal Learning Test; SCI, subjective cognitive impairment. Values are shown as median (range) or as numbers per group. aMann-Whitney U non-parametric test used for continuous variables and Chi-Square test for the categorical variable (gender). bRaw values (neither log2-transformed nor autoscaled). cAnalysis based on 59 subjects, one missing value. d
In order to relate changes in lipid species levels to typical AD pathology and related measures, we performed a lipidomic analysis of CSF from patients described in Table 1. The selection of the m/z features best distinguishing between CSF profile groups was done using Mann-Whitney U non-parametric tests. Of the 1008 features detected, eight significantly differed in levels between groups within both the discovery and validation sets in more than 200 of the 1000 (20%) bootstrap replicates (Fig. 2) and were selected for further identification. Eight m/z features reached significance in 10–20% of the replicates, 23 in 2–10% of the replicates and 969 in less than 2% of the replicates. The distribution of all m/z features after 1000 boostrap replications is depicted in Supplementary Figure 2.

Frequencies of m/z features significantly distinguishing between CSF profiles (AD and non-AD) in both discovery and validation sets (
Table 2 presents the lipid species annotations that were automatically assigned to the eight selected m/z features using the LipidBlast spectral library. It also includes comparisons in levels of each feature between CSF profiles. Annotations of compounds searching the library are based on similarity matching of elements, including m/z peaks and retention time (RT). The lipid species selected belonged to four different lipid categories; glycerophospholipids (PE 16:23/16:4, PMeOH 20:3/22:6), sphingolipids (ceramide d18:1/18:0), branched fatty acid esters of hydroxy fatty acids (FAHFA 18:5/22:2, 18:4/20:3 and 20:3/22:4), and monogalactosyldiacylglycerol (MGDG 10:0/10:0). Of the eight lipid species annotated via LipidBlast library search, a reference standard was only run for C18 ceramide (d18:1/18:0). A comparison of the MS/MS fragmentation spectra of the C18 ceramide standard confirmed the annotation of C18 ceramide in the CSF samples (Supplementary Figure 3). The other seven annotated lipid species remain unconfirmed, as a comparison of each MS/MS fragmentation spectra to the LipidBlast reference spectra did not clearly confirm the correct chemical structure. All measured m/z features are listed in Supplementary Table 1A (positive ionization mode) and 1B (negative ionization mode).
LipidBlast annotations of the eight selected m/z features. Levels for each feature compared between CSF AD and non-AD profile groups
m/z, mass-to-charge ratio; RT, retention time; Cer, ceramide; FAHFA, fatty acid ester of hydroxyl fatty acid; MGDG, monogalactosyldiacyl-glycerol; PE, phosphatidylethanolamines; PMeOH, phosphatidyl methanol. aBased on normalized peak area arbitrary units (A.U.), log2-transformed and autoscaled before analysis. bMann-Whitney U non-parametric test used.
As the structure for C18 ceramide was confir-med with high confidence, it was selected for fur-ther analysis. Figure 3 presents the levels of CSF C18 ceramide by different CSF profiles. As can be observed, levels of C18 ceramide were elevated in the CSF AD profile group compared to the non-AD group (

Comparison in levels of CSF C18 ceramide by CSF profile (non-AD and AD). *
Linear regression (unadjusted and adjusted for gender, age, and education) was performed to estimate the relationships of established AD (Aβ42, T-tau, P-tau), inflammatory (YKL-40, S100B, GFAP) and neuronal degeneration (NFL) markers with CSF C18 ceramide (Table 3). Levels of Aβ42 were negatively associated with C18 ceramide when adjusted for age, gender, and education (st. β= –0.36,
Linear regression estimates (unadjusted and adjusted) for the association between various measures and levels of CSF C18 ceramide
A.U., arbitrary units; RAVLT, Rey Auditory Verbal Learning Test. Numbers present standardized beta coefficients (st.β). *Adjusted for age, gender, and years of education. aValues log2-transformed and autoscaled before analysis. bAnalysis based on 59 subjects, with one missing.
The statistically significant relationships between the levels of proteins and C18 ceramide from (

Pearson’s correlations between CSF levels of a) Aβ42, b) T-tau, c) T-tau/ Aβ42, d) S100B and CSF C18 ceramide. Normalized peak area arbitrary units (A.U.). All CSF measures were log2-transformed and autoscaled before analysis. C18 ceramide levels significantly correlated with levels of the core AD markers Aβ42, T-tau, and inflammatory marker S100B.
Pearson‘s correlations between selected CSF markers and levels of CSF C18 ceramide within CSF profiles
A.U., arbitrary units. aValues log2-transformed and autoscaled before analysis.
DISCUSSION
We used an untargeted lipidomic approach to identify CSF lipid species associated with a signature CSF profile reflecting AD pathology. A total of 1008 m/z features were detected, with eight selected as candidate markers. Out of these, one was fully confirmed as corresponding to the lipid species C18 ceramide. Our results showed that C18 ceramide levels were higher among subjects with a CSF AD profile. Relationships were also detected between established AD markers (Aβ42, T-tau) and C18 ceramide. Higher levels of C18 ceramide associated with lower levels of Aβ42 and higher levels of T-tau. In addition, levels of the inflammatory marker S100B positively related to C18 ceramide. Overall, our results indicate that CSF C18 ceramide levels could increase during pathological changes in the early stages of AD.
Ceramides, the core constituents of sphingolipid metabolism, are composed of a sphingosine backbone linked to a fatty acid chain of varying carbon atom length (C14-C26) [59]. Alterations in sphingolipid metabolism have been observed in healthy aging and in neurodegenerative diseases, including AD [60]. They play essential roles in the structural stability of membranes and as signaling molecules affect differentiation, proliferation, inflammation, and apoptosis [61, 62]. Ceramides are synthesized via two main pathways in eukaryotic cells [63]. In the salvage pathway, sphingomyelin is hydrolyzed through sphingomyelinase (SMase) to produce ceramide, which can be further metabolized to sphingosine by ceramidase. Ceramides can also be generated through the de novo pathway via anabolism of serine and palmitate. Results from cellular and animal studies suggest that both direct and indirect mechanisms by which ceramides can contribute to an increase in Aβ levels and AD pathogenesis [64]. Ceramides stimulate Aβ generation by stabilizing β-secretase enzyme BACE1 and increasing its half-life [65, 66]. Furthermore, soluble and fibrillar forms of Aβ can induce degradation of SM to ceramide by SMases [67, 68] through oxidative stress-mediated mechanisms [69]. This positive loop of ceramide production possibly contributes to immune activation and neuronal loss in AD. Long-chain ceramides, specifically C18 ceramide, have also been linked to tau phosphorylation through modulation of PP2A activity [70–74]. Several postmortem studies have found increased levels of ceramide in brain tissues of AD patients compared to healthy controls [60, 75–77]. Two of the studies [60, 76] examined different ceramide species, with long-chain ceramide levels (C18 and C24) being significantly elevated in the AD group. Han et al. [75] observed the highest elevation of total ceramide in the brains of patients with mild AD, compared to those with severe AD and to controls, indicating early changes in the pathological processes of AD. To further confirm the implications of ceramide metabolism in AD, upregulation [78] and increase in activity of enzymes [77] controlling ceramide synthesis have been observed in brain areas (temporal and frontal cortices) affected in the early stages of the disease.
While postmortem studies are essential,
There are several limitations to our study. First, our sample size was relatively small, with a vast number of features compared to the number of subjects, resulting in little power and a considerable risk of type-II errors. Corrections for multiple comparisons were, therefore, not performed when comparing levels of features between groups in the discovery set. A more significant validation set would also have allowed for evaluation of accuracy in distinguishing between CSF profile groups using multiple lipid species simultaneously as predictors. Second, our study did not include a healthy control group or patients with moderate to severe AD. It also included very few participants with other dementias. It would be of interest to examine the relationships between ceramides and AD-related markers in CSF in a more diverse cohort, for a better evaluation of C18 ceramide as a marker of AD progression and severity, preferably in a longitudinal study. Third, information about the
In summary, our results indicate that CSF C18 ceramide levels associate with established markers of AD pathology (Aβ42 and T-tau) and inflammation (S100B) at the symptomatic pre- and early stages of dementia. These findings suggest that ceramide metabolism could influence the pathophysiological processes during the early stages of AD. Furthermore, ceramides could potentially serve as therapeutic targets, with strategies aiming at reducing ceramide levels to slow down the progression of the disease. Longitudinal studies are, however, needed to validate the pathological implications of these results.
Footnotes
ACKNOWLEDGMENTS
This study was supported by the St. Josef's Hospital Fund, Reykjavik, Iceland, the Landspitali University Hospital Research Fund and the Icelandic Research Fund of the Icelandic Centre for Research (163172-051). The authors thank all the subjects of The Icelandic MCI study for their participation. They also wish to thank Kristin H. Hannesdottir for managing participant administration and the staff of the LUH Memory Clinic.
