1. Introduction {#s0005}
===============

Alzheimer\'s disease (AD) is the main cause of dementia in the elderly, responsible for about half of the nearly 47 million dementia cases worldwide ([@bb0220]). Although the disease is defined by the accumulation of amyloid beta plaques and neurofibrillary tau tangles in the brain ([@bb0040]), other pathological processes can be identified such as oxidative stress, vascular dysfunction, and inflammation ([@bb0005]; [@bb0255]; [@bb0135]). Additionally, defects in essential metabolic processes for energy supply and membrane function have been implicated in AD ([@bb0125]; [@bb0165]).

The human brain is highly vulnerable to disturbances in energy metabolism, due to its relatively large energy consumption. Previous AD studies demonstrated alterations in global and cellular energy metabolism. FDG-PET studies have shown glucose hypometabolism in the retrosplenial cortex (RSC) and medial temporal lobe in people with mild cognitive impairment (MCI) and AD patients, as well as in cognitively normal carriers of the APOE ε4 allele ([@bb0160]; [@bb0190]; [@bb0210]; [@bb0215]). Furthermore, the enzyme creatine kinase (CK), obtained from post-mortem AD tissue, shows reduced activity compared with samples that are free from neurological disease ([@bb0010]; [@bb0065]). The CK reaction, key to balance brain energy metabolism, can quickly replenish adenosine triphosphate (ATP) from the energy buffer phosphocreatine (PCr), when local energy demands suddenly increase ([@bb0120]). This reaction also enhances the efficiency of mitochondrial oxidative phosphorylation (OXPHOS) by keeping adenosine diphosphate (ADP) sufficiently available ([@bb0230]) and prevents acidification by maintaining pH (for a review see [@bb0275]).

One of the earliest pathological changes in AD is synaptic dysfunction, which correlates well with cognitive dysfunction and disease severity ([@bb0245]). Alterations in neuronal membrane composition, vital for synaptic transmission, have been linked to synaptic dysfunction in AD ([@bb0200]). For instance, post mortem studies found reduced levels of the major membrane components phosphatidylcholine, phosphatidylethanolamine and phosphatidylinositol ([@bb0165]; [@bb0185]), as well as altered activity of catabolic and anabolic enzymes, suggesting compensatory metabolic changes to reduce the rate of neuronal membrane loss ([@bb0225]).

High energy phosphates, such as ATP and PCr, and metabolites of phospholipid membrane metabolism can be assessed in vivo by phosphorus magnetic resonance spectroscopy (^31^P--MRS). This technique has been applied in AD (e.g. [@bb0080]; [@bb0130]; [@bb0180]; [@bb0235]), but previous studies, essentially focusing on single brain regions, were severely hampered by small sample sizes and low spectral resolution, resulting in a wide disparity in findings. Furthermore, differences in the control groups that were used, in the brain regions that were studied and in the disease stages and medication statuses of the patients, may have contributed to apparent inconsistencies in the literature.

The aim of the current study was to investigate whether abnormal phospholipid and energy metabolism can be detected non-invasively in early stage AD patients by ^31^P MR spectroscopic imaging (^31^P--MRSI). This investigation was designed to overcome the limitations of previous studies. First, we performed the measurements at 3 T field strength with proton (^1^H) decoupling to increase spectral resolution and sensitivity. Secondly, by applying ^31^P--MRS in a 3-dimensional imaging mode we could investigate multiple brain regions simultaneously. We specifically addressed regions that are of interest in AD, namely, the hippocampus, the RSC, and the anterior cingulate cortex (ACC). Finally, we selected a well-defined drug-naïve patient group and an age- and gender-matched control group.

2. Materials and methods {#s0010}
========================

2.1. Subjects and design {#s0015}
------------------------

All visits of patients and controls to the Radboud university medical center (Nijmegen, the Netherlands) took place between 2012 and 2015. Patients aged ≥50 years with a diagnosis of possible or probable AD according to the revised National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer\'s Disease and Related Disorders Association 2011 criteria ([@bb0140]) with evidence of the pathophysiological process (i.e. from structural MRI or cerebrospinal fluid biomarker assays) and with a minimum Mini Mental State Examination (MMSE) score of 20, were recruited from the hospital\'s memory clinic or by referral from regional hospitals. All AD patients participated in a randomized controlled trial on the effect of a medical food on brain phospholipid metabolism ([@bb9000]; Dutch Trial Register 3346). Only baseline data were used in the current study. Healthy control subjects, age- and gender-matched to the AD group, were screened by telephone before being invited to the clinic. All subjects were drug-naïve for AD medication (cholinesterase inhibitors and NMDA-antagonists) and were free from neurological or psychiatric disorders (other than dementia, for the AD patients). Patients nor healthy controls used nutritional supplements containing docosahexaenoic acid, eicosapentaenoic acid or \>200% of the recommended daily allowance of vitamins B, C, or E. All subjects were screened for MRI contra-indications before inclusion in the study.

Written informed consent was obtained from all subjects and from the informal caregivers of the AD patients. The local ethics committee reviewed and approved the protocol and the study was conducted in accordance with the Declaration of Helsinki.

Medical history, medication use and MMSE were recorded from all subjects, as well as date of birth, sex, ethnicity, smoking habits, alcohol consumption, and family history of AD. Level of education was classified according to Verhage ([@bb0270]) (comparable with the International Standard Classification of Education) and categorized in three groups: primary education or lower (low), junior vocational training (middle), and senior vocational or academic training (high).

2.2. MR acquisition {#s0020}
-------------------

MRI and MRS were performed on a Magnetom Tim Trio 3T MR system (Siemens, Erlangen, Germany) with a dual-tuned ^1^H/^31^P volume head coil (Rapid Biomedical, Wuerzburg, Germany). High resolution MR images were acquired with a T1-weighted magnetisation-prepared rapid gradient-echo sequence (TR = 2300 ms, TE = 3.16 ms, TI = 1100 ms, 15° flip-angle, 176 sagittal slices, slice-matrix size = 256 × 256 mm, slice thickness = 1 mm, voxel-size = 1 × 1 × 1 mm, TA = 6:25 min).

^31^P-MR spectra were acquired with whole brain 3D MRSI (a pulse-acquire sequence with 40° flip-angle, acquisition delay of 0.7 ms including 0.3-ms gradient phase encoding time, 512 ms acquisition duration with WALTZ4 ^1^H decoupling during the first 256 ms with γB1 = 250 Hz ensuring (γB1) ^2^ \> (J~PH~)^2^, TR = 2000 ms, FOV 260 × 260 × 260 mm; matrix size = 10 × 10 × 10, acquisition time 13:08 min). Prior to acquisition, an automatic procedure was applied to optimize the magnetic field homogeneity of a rectangular shim volume maximizing brain coverage while bound by the skull. Subsequently, if necessary, the field homogeneity was further adjusted manually until a value of \<25 Hz was obtained for the full width at half maximum of the absolute-value ^1^H MR water signal. K-space was sampled with a weighted elliptical phase-encoding scheme with four averages. The volume of interest was centered on the midline and parallel to the line from the anterior commissure to the posterior commissure. Spatial post-processing consisted of zerofilling to a matrix size of 16 × 16 × 16. The nominal volume of the measured voxels is about 17.5 cm^3^, the effective voxel at 64% of the spatial response function has a spherical shape with a volume of about 40 cm^3^ ([@bb0195]).

2.3. MRS and MRI data analyses {#s0025}
------------------------------

Four regions of interest (ROI) were selected for analysis from the 3D ^31^P--MRSI data: ACC, RSC, and left and right hippocampus (HL and HR). Guided by the T1-weighted images, the MRSI grid was shifted in the x, y and z dimension to position voxels in the ROI\'s (by AR; [Fig. 1](#f0005){ref-type="fig"}). For the ACC, a voxel was selected anterior to the genu of the corpus callosum, with the inferior border aligned with the line from the anterior commissure to the posterior commissure (AC-PC line), and on the midsagittal plane. For the RSC, a voxel on the midsagittal plane was selected posterior to the splenium of the corpus callosum, with the inferior border aligned with the AC-PC line. For the HL and HR, a voxel was selected at the anterior end of each hippocampus, with the inferior and lateral border of the voxel aligned with the inferior and lateral side of the hippocampus.Fig. 1Voxel selection of ^31^P MR spectra displayed on sagittal (left) and coronal (right) anatomical images: anterior cingulate cortex (blue), retrosplenial cortex (yellow), left and right hippocampus (purple). The nominal voxel size is indicated by squares and an approximation of the effective spherical voxel size is indicated by circles. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)Fig. 1

The software package Metabolite Report (Siemens Healthcare, Erlangen, Germany) was used for post-processing (i.e. 100 ms exponential filter, zero-filling, phase correction, baseline correction) and for automatic fitting of the spectra in the time-domain using prior knowledge appropriate for ^31^P MR spectra (own spectral data and from [@bb0095]). Eleven well-resolved resonance peaks were fitted: the phospholipid metabolites phosphoethanolamine (PEth), phosphocholine (PCh), glycerophosphoethanolamine (GPEth) and glycerophosphocholine (GPCh), the high-energy phosphorus molecules PCr, ATP (α-ATP, β-ATP and γ-ATP), nicotinamide adenine dinucleotide (NAD(H)) and inorganic phosphate (Pi), and membrane-bound phospholipid (details of this procedure are described in the Appendix).

Both a quantitative evaluation of the fitting results and a visual quality control were performed. Quantitatively, only metabolite fits with a Cramer--Rao lower bound of ≤30% were considered reliable. Qualitatively, two MR spectroscopists (AR and MvdG) independently checked the spectra by visual inspection of the original spectra and the fitting results. If a metabolite peak was visually present, and its fit was assigned to the correct resonance---giving a minimal residue in the subtraction spectrum---the fitting result was accepted.

 The integral of each metabolite resonance was expressed as a percent area of the total phosphorous signal in the corresponding spectrum. Hereby, the data were normalized, and corrected for cerebrospinal fluid (CSF) content and thus for differences in local atrophy. Phosphomonoesters (PME) were calculated from the sum of PEth and PCh, phosphodiesters (PDE) from GPEth and GPCh and total ATP (tATP) from the sum of α-ATP, β-ATP and γ-ATP. In addition, the ratios PME/PDE, PE/GPE, PC/GPC and PCr/Pi were computed. Intracellular pH was determined from the chemical shift difference between the PCr and Pi resonance peaks ([@bb0175]).

The T~1~-weighted images were segmented into grey matter (GM), white matter (WM) and CSF using automatic segmentation software (by AR; SPM8, Welcome Trust Centre for Neuroimaging, London, UK; VBM8, Structural Brain Mapping Group, Jena, USA). Binary tissue maps (matrices) of GM, WM, CSF were convolved with the, digitally sampled, mathematically-modelled three dimensional point-spread function of the MRSI, using Matlab 2012b (The Mathworks, Inc., Matrick, MA), to obtain theoretically correct values of the contribution of tissue type to each ROI (each consisting of one MRS voxel). GM-to-WM ratio (GM/WM) was calculated for each ROI for each subject.

2.4. Statistical analyses {#s0030}
-------------------------

Two-factorial mixed ANOVAs were used to test differences in tissue content (GM, WM, and CSF) and GM/WM in each ROI between mild AD patients and controls. Group and brain region (ROI) differences were analyzed for each spectroscopy outcome parameter using a mixed model with group (AD or control), sex, brain region and its interaction with group as fixed factors, considering brain region as a within subject factor, and adjusting for age. An unstructured variance-covariance matrix for brain region was selected. The group-by-brain region interaction was only kept in the model if *p* \< 0.10. In the covariate analysis, GM/WM in each ROI was added to the final primary model of each parameter, to investigate the effect of tissue composition on the effects of group and brain region. Although the inclusion of GM/WM in the final model led to a small change (10--20%) in many parameter estimates, this was mainly without a change in main and/or interaction effects. Therefore we report the primary model in tables and figures and refer to the effect of GM/WM only where it led to a change in the main and/or interaction effects, and with \>20% change in the parameter estimates of group or brain region. Goodness of fit was assessed by comparing Bayesian Information Criteria (BIC) and homoscedasticity and normality of residuals were assessed by visual inspection. Statistical analyses were performed using SAS 9.2 Software for Windows (SAS Institute Inc., Cary, NC, USAf). Statistical significance was set at *p* \< 0.05 without correction for multiple testing.

3. Results {#s0035}
==========

3.1. Subjects {#s0040}
-------------

A total of 33 mild AD patients and 35 healthy control subjects were enrolled in the study and subjected to MRI and MRS measurements. Two AD patients had MR spectra of low quality (due to severe motion and a dental implant) and four control subjects were found to meet an exclusion criterion during or after the visit (two subjects took nutritional supplements, two subjects presented with structural neurological damage from ischemic infarcts and trauma). Therefore, 31 AD patients and 31 controls were included in the final analysis. The groups were similar with respect to age and gender (see [Table 1](#t0005){ref-type="table"} for subject characteristics).Table 1Subject characteristics.Table 1AD\
*n* = 31Controls\
n = 31Age73.4 years (6.8)73.5 years (6.3)Gender (%)13 males (42%)15 males (48%)MMSE23.2 (2.0)28.1 (1.4)Time since AD diagnosis1.6 months \[0.2--15.3\]--BMI25.4 kg/m^2^ (3.8)26.3 kg/m^2^ (4.5)Educational level (low/middle/high)9/16/64/9/18Smoking (never/past/current)13/17/16/22/3≥1 relatives with AD (n)197[^2]

3.2. ROI tissue contribution {#s0045}
----------------------------

The interaction between group and brain region for GM and CSF contributions to the total volume of each ROI was not significant, but GM and CSF differed significantly between groups, and between brain regions (all *p* \< 0.001). GM was lower and CSF was higher in AD patients than in healthy control subjects (mean % ± SD, GM:36.3 ± 4.7 vs. 40.5 ± 4.0; CSF: 22.5 ± 4.7 vs. 18.1 ± 3.5). overall, GM was highest in the right hippocampus, followed by left hippocampus, and then followed by equal amounts in RSC and ACC (data not shown). CSF in RSC and ACC was higher than in both hippocampi. There was a significant interaction between group and brain region for WM (*p* = 0.004) and GM/WM (*p* = 0.021). WM was higher in AD patients compared with controls in the right hippocampus (41.5 ± 3.0 vs. 40.0 ± 2.6), but not in the other brain regions. GM/WM was lower in AD patients than in controls in ACC, and both hippocampi, but not in RSC (data not shown).

3.3. MRS results {#s0050}
----------------

High quality MR spectra were obtained with well-resolved PME, PDE, Pi, PCr and ATP signals for the brain regions evaluated in this study ([Fig. 2](#f0010){ref-type="fig"}). The linewidth of PCr varied between 8 and 13 Hz in the different brain regions (ACC: mean ± SD 13.5 ± 3.3 Hz, RSC: 8.1 ± 2.6 Hz, HL: 11.2 ± 2.6 Hz, HR: 12.5 ± 3.0 Hz). We did not observe a difference in the linewidth of PCr between AD and controls (mean ± SD 11.5 ± 3.9 Hz and 11.1 ± 3.0 Hz, not tested). ROI tissue contributions showed typical differences with lower GM and higher CSF in AD patients compared with controls.Fig. 2Representative spectra from the anterior cingulate cortex (ACC), retrosplenial cortex (RSC), left hippocampus (HL) and right hippocampus (HR) of an Alzheimer\'s disease patient (73 years old). Zero filling to 4096 data points and an 8 Hz Gaussian filter were applied for display purposes. PEth, phosphoethanolamine; PCh, phosphocholine; Pi, inorganic phosphate; GPEth, glycerophosphoethanolamine; GPCh, glycerophosphocholine; MP, membrane phospholipids; PCr, Phosphocreatine; ATP, adenosine triphosphate; NAD(H), nicotinamide adenine dinucleotide; ppm, parts per million.Fig. 2

In the statistical analysis, the interaction between group and brain region was kept in the model for PCr, γATP and βATP, because it survived our criterion of *p* \< 0.10. For all other outcome parameters, the interaction term had a *p*-value \>0.10 and was removed from the model.

3.4. Differences in energy metabolites and pH between groups {#s0055}
------------------------------------------------------------

A significant interaction between group and brain region was found for PCr, showing that normalized PCr signals were higher in mild AD patients compared with healthy control subjects in the RSC (Least Square (LS) mean difference = 0.86, SEM = 0.29, *p* = 0.004), HL (LS mean difference = 1.52, SEM = 0.52, *p* = 0.005) and HR (LS mean difference = 1.47, SEM = 0.38, *p* \< 0.001), but not in the ACC (*p* = 0.947) ([Fig. 3](#f0015){ref-type="fig"}A). The interaction between group and brain region was not significant for γATP and βATP. A main effect of group was found for PCr/Pi, showing that this ratio was higher in mild AD patients compared with healthy controls (LS mean difference = 0.22, SEM = 0.11, *p* = 0.046) ([Fig. 3](#f0015){ref-type="fig"}B). Also pH was increased in AD patients compared with controls (LS mean difference = 0.008, SEM = 0.004, *p* = 0.032) ([Fig. 3](#f0015){ref-type="fig"}C). No group differences were found for γATP, βATP, αATP, total ATP, Pi, and NAD(H) (all *p* \> 0.05). For a summary of group effects see [Table 2](#t0010){ref-type="table"}.Fig. 3Normalized phosphocreatine (PCr; A), ratio of PCr to inorganic phosphate (PCr/Pi; B) and pH (C) in patients with Alzheimer\'s disease (AD) and healthy control (HC) subjects. PCr was normalized by expressing the integral as a percent area of the total phosphorous signal in the corresponding spectrum. Data represent Least Square means±SEM, \**p* \< 0.05 AD vs HC. ACC, anterior cingulate cortex; RSC, retrosplenial cortex; HL, left hippocampus; HR, right hippocampus.Fig. 3Table 2Main effect of group (healthy controls vs Alzheimer\'s disease) on phosphorus metabolite signals and tissue pH for all investigated regions.Table 2ADControlp-ValuePCr\*19.04 ± 0.18\
*n* *=* *124*18.08 ± 0.18\
*n* *=* *124***\<0.001**Pi5.30 ± 0.11\
*n* *=* *121*5.47 ± 0.11\
*n* *=* *123*0.240PCr/Pi3.71 ± 0.09\
*n* *=* *121*3.49 ± 0.09\
*n* *=* *123***0.046**total ATP42.88 ± 0.27\
*n* *=* *124*43.17 ± 0.27\
*n* *=* *124*0.445NAD(H)3.96 ± 0.11\
*n* *=* *103*4.09 ± 0.11\
*n* *=* *96*0.388PEth8.19 ± 0.11\
*n* *=* *119*8.33 ± 0.11\
*n* *=* *122*0.375PCh3.90 ± 0.07\
*n* *=* *112*3.93 ± 0.07\
*n* *=* *114*0.757GPEth4.32 ± 0.11\
*n* *=* *109*4.31 ± 0.11\
*n* *=* *117*0.942GPCh5.96 ± 0.15\
*n* *=* *108*6.08 ± 0.15\
*n* *=* *109*0.526pH7.037 ± 0.003\
*n* *=* *121*7.028 ± 0.003\
*n* *=* *123***0.032**[^3]

3.5. Differences in phospholipid metabolites between groups {#s0060}
-----------------------------------------------------------

No interaction or group difference between mild AD patients and healthy control subjects were found for any of the phospholipid metabolites (PEth, PCh, GPEtn, GPCh, PME, and PDE; [Table 2](#t0010){ref-type="table"}) or their ratios (PEth/GPEth, PCh/GPCh, and PME/PDE) (all *p* \> 0.05).

3.6. Differences in phosphorus metabolites and pH between brain regions {#s0065}
-----------------------------------------------------------------------

Brain region significantly affected all outcome parameters (all *p* \< 0.05), except PCh and PME (*p* \> 0.05). The normalized signals of total ATP, αATP, βATP, γATP, NAD(H), PEth, GPEtn, GPCh, and PDE were lower in the RSC than in the ACC and both hippocampi (all p \< 0.05), while the pattern was reversed for the ratios PCh/GPCh and PME/PDE (*p* \< 0.001). We found no difference between the left and right hippocampus for any of the energy and phospholipid metabolites or pH (p \> 0.05), except for αATP (LS mean difference = 0.65, SEM = 0.26, *p* = 0.015). See [Table 3](#t0015){ref-type="table"} for an overview of the differences between brain regions.Table 3Main effect of brain region on phosphorus signal integrals and tissue pH.Table 3ACCRSCHLHRp-ValuePCr\*17.75 ± 0.18^†,‡,§^\
*n* *=* *62*19.16 ± 0.14^¶^\
*n* *=* *62*18.47 ± 0.26\
*n* *=* *62*18.87 ± 0.19\
*n* *=* *62***\<0.001**Pi5.34 ± 0.12\
*n* *=* *61*5.65 ± 0.10^¶^\
*n* *=* *62*5.21 ± 0.16\
*n* *=* *61*5.34 ± 0.18\
*n* *=* *60***0.027**PCr/Pi3.42 ± 0.09^‡,§^\
*n* *=* *61*3.46 ± 0.07^¶,\#^\
*n* *=* *62*3.76 ± 0.14\
*n* *=* *61*3.76 ± 0.12\
*n* *=* *60***0.009**total ATP42.58 ± 0.38^†,‡,§^\
*n* *=* *62*40.06 ± 0.33^¶,\#^\
*n* *=* *62*45.24 ± 0.38\
*n* *=* *62*44.21 ± 0.41\
*n* *=* *62***\<0.001**NAD(H)4.39 ± 0.12^†^\
*n* *=* *49*3.47 ± 0.11^¶,\#^\
*n* *=* *61*4.17 ± 0.16\
*n* *=* *45*4.08 ± 0.14\
*n* *=* *44***\<0.001**PEth8.45 ± 0.18^†^\
*n* *=* *61*7.81 ± 0.12^¶,\#^\
*n* *=* *62*8.31 ± 0.15\
*n* *=* *60*8.46 ± 0.18\
*n* *=* *58***\<0.001**PCh3.99 ± 0.10\
*n* *=* *55*3.96 ± 0.10\
*n* *=* *61*3.94 ± 0.09\
*n* *=* *56*3.76 ± 0.11\
*n* *=* *54*0.429GPEth4.76 ± 0.15^†,,§^\
*n* *=* *51*3.72 ± 0.09^¶,\#^\
*n* *=* *61*4.42 ± 0.13\
*n* *=* *57*4.36 ± 0.11\
*n* *=* *57***\<0.001**GPCh6.69 ± 0.26^†^\
*n* *=* *46*5.01 ± 0.12^¶,\#^\
*n* *=* *61*6.11 ± 0.15\
*n* *=* *56*6.28 ± 0.16\
*n* *=* *54***\<0.001**pH7.023 ± 0.004^‡,§^\
*n* *=* *61*7.017 ± 0.002^¶,\#^\
*n* *=* *62*7.039 ± 0.004\
*n* *=* *61*7.050 ± 0.005\
*n* *=* *60***\<0.001**[^4]

3.7. Influence of tissue content {#s0070}
--------------------------------

In the covariate analysis GM/WM was added to the final model. This did not change the effects of group, brain region or the interaction between group and brain region on PCr, Pi, total ATP, αATP, βATP, γATP, NAD(H), pH or any of the phospholipid metabolite signals. However, for PCr/Pi the main effect of brain region was reduced to a trend (*p* = 0.062) and the parameter estimates changed substantially, when GM/WM was included in the model. There was no change in the effect of group on PCr/Pi.

4. Discussion {#s0075}
=============

In this study we acquired high-quality ^31^P-MRSI data of the brain in a large group of mild AD patients and in a well-matched healthy elderly population. High spectral resolution enabled us to gather information on high-energy phosphates, as well as on specific phospholipid metabolites. MR spectra of the four brain regions that are among the most studied in AD were investigated and for the first time, these areas were measured simultaneously by ^31^P-MRSI. We detected increased phosphocreatine and pH, but no changes in phospholipid metabolite signals in mild AD patients. Interestingly, the increase in PCr was observed in the RSC and in the hippocampi, regions that both show early pathological changes in AD, but not in the ACC, a region known to be involved in the disease at a later stage. This suggests that the PCr abnormality expands across brain regions in a similar fashion as other AD pathologies, such as tau ([@bb0040]). Furthermore, in both groups pH, high energy phosphates, and most phospholipid metabolites differed across brain regions.

As in vivo human ^31^P MRS brain studies are nearly always performed under signal saturating conditions, changes in T1 relaxation time may affect metabolite signal intensities ([@bb0115]). Although we cannot exclude that such an effect has contributed to the increased PCr signal, it would require a rather drastic decrease in PCr\'s T1 value from about 3 to 2 s. Moreover, no increases in PME and PDE signals were detected, while those would be affected more severely by a T1 decrease ([@bb0110]). A potential cause for relaxation enhancement is a higher brain iron accumulation in AD patients ([@bb0240]), but this would increase the signal linewidth, which was not observed. Thus, the most likely explanation for the increased PCr signal is a higher tissue concentration of this compound.

Elevated PCr signals have been reported before for AD patients in comparison with young controls ([@bb0130]), but as PCr also increases with age in healthy persons ([@bb0055]; [@bb0075]; [@bb0115]) it is unclear whether this finding was due to an aging effect. The current results showed increased PCr, and PCr/Pi, in AD patients compared with elderly controls of the same age. In contrast, two earlier studies showed a decrease in prefrontal PCr in mild, but not severe, AD ([@bb0180]) and decreased PCr/Pi in the frontal lobe ([@bb0235]). However, these studies were performed at 1.5 T with only surface coil localization and involved small sample sizes.

As the concentration of high energy phosphates depends on the balance between energy production and energy utilization, increases in steady state PCr levels and PCr/Pi reflect a redistribution in the content of metabolites involved in the connected CK and ATPase pseudo-equilibria ([@bb0070]). This may be the consequence of reduced utilization of ATP ([@bb0075]), caused by synaptic dysfunction or loss of other energy requiring functionalities ([@bb0050]; [@bb0205]). Although these changes are expected to lower ATP production and thus PCr, the final balance between reduced energy production and consumption may still result in a higher PCr level. To maintain equilibrium, the increased PCr should be accompanied by changes in the levels of other substrates of the CK reaction, such as a higher ATP concentration. Previous studies reported an increase in (γ)ATP ([@bb0130]; [@bb0145]; [@bb0180]), although no alterations in ATP levels were found in the current nor in several other studies (e.g. [@bb0080]; [@bb0090]; [@bb0235]). Alternatively, a higher PCr may result from an increased total amount of cellular creatine (Cr). However, in vivo measurements of total Cr content in MCI and AD patients provide no clear evidence that Cr is increased in the disease ([@bb0100]; [@bb0260]). Finally, decreases in ADP and/or H^+^ concentration may accompany the PCr increase. Indeed, a slight increase in pH was detected in mild AD patients, but it does not fully compensate for the increase in PCr when equilibrium is assumed. Changes in PCr levels have also been associated with decreased CK activity in the AD brain ([@bb0075]; [@bb0130]). However, it is unlikely to be a direct effect of this decrease as even in cytosolic CK knockout mice the remaining CK activity (\<20%) is sufficient to preserve equilibrium ([@bb0285]). Clearly, a strongly affected reaction rate will affect processes that require rapid (local) energy supply, and a higher PCr level may represent a flux adaptation to this condition.

There is some evidence for involvement of PCr in processes beyond the PCr-CK system. Since glutamate uptake by synaptic vesicles is stimulated by PCr, independent of ATP or CK ([@bb0280]), increased PCr may prevent glutamate excitotoxicity ([@bb0025]). Furthermore, PCr may be able to bind to and stabilize phospholipid membranes ([@bb0250]) and may function as an osmolyte ([@bb0030]) and a neurotransmitter ([@bb0015]). In addition, differences in metabolite levels may be explained by variations in tissue and cell type ([@bb0105]; [@bb0265]), although our results remained unchanged when adjusting for grey and white matter. A prominent neuropathological feature in AD is astrogliosis ([@bb0170]). As this alters the neuron-to-glia ratio within grey and white matter ([@bb0020]), it is possible that morphological differences between AD patients and healthy aged controls play a role in the alterations in PCr and pH. Furthermore, atrophy may introduce slight differences in the brain area that the measured volumes encompass, which effect may not be entirely captured by adjusting for grey and white matter. However, ROIs were placed such that volume of the actual area of interest was maximized and the measured signals are dominated by the reported brain regions of interest.

Contrary to our expectation, we did not observe differences between mild AD patients and controls in any of the phospholipid metabolites, although preclinical and post-mortem studies show alterations in membrane function and lipid composition in AD ([@bb0165]; [@bb0225]). No previous in vivo studies have reported on PEth, PCh, GPEth or GPCh, but several studies reported changes of their total amounts, i.e. PME, PDE or their ratio, in AD (e.g. [@bb0085]; [@bb0090]; [@bb0145]), although others reported no differences in these metabolites ([@bb0035]; [@bb0045]; [@bb0155]; [@bb0235]). The most consistent finding seems to be an increase in total PME in the prefrontal cortex in AD compared with elderly controls ([@bb0060]; [@bb0085]; [@bb0180]). In the present study, part of the frontal cortex was covered in the ACC measurement, but we did not find changes in PEth, PCh, nor total PME. However, changes in the complex phospholipid cycle may be too small at this disease stage to be detected with the current methodology. Although we had sufficient spectral resolution to resolve the phosphomonoester signal of PEth and PCh, and the phosphodiesters signal of GPEth and GPCh, these may still reflect metabolites connected to different phospholipid pools. Moreover, differential expression of phosphatidylethanolamine and phosphatidylcholine species across hippocampal subfields was recently shown in AD ([@bb0150]).

MRSI revealed that high energy phosphates, pH and most phospholipid metabolites vary extensively across brain regions that are of interest in AD. These differences existed in both mild AD patients and controls, and were not explained by regional variation in grey and white matter content, with the ratio PCr/Pi as the sole exception. This indicates that the investigated regions are metabolically different from each other. It also underlines that regional differences should be taken into account in any cross-sectional or longitudinal metabolic study. Moreover, as shown here for PCr, changes in metabolite levels in disease may not be the same for all brain regions.

The current study only reports on relative metabolite levels, as obtaining absolute quantifications requires several assumptions, which adds complexity and reduces reliability. Moreover, when the variation in atrophy between and within groups is large, relative measurements are preferred over absolute. A second limitation may be that our results have not been corrected for multiple testing, because of the exploratory nature of this study. Replication of these results in an independent population would strengthen our findings. In future research it would be valuable to include preclinical AD patients or an at-risk population to clarify the disease stage in which abnormalities in PCr and pH appear, as well as perform longitudinal (intervention) studies to follow these abnormalities. In addition, magnetization transfer experiments to determine the CK and ATPase reaction rates with measurements of absolute concentrations of PCr and ATP to determine fluxes ([@bb0070]), would elucidate whether the dynamics of the PCr-CK and ATPases systems are affected in AD.

4.1. Conclusion {#s0080}
---------------

To the best of our knowledge, this is the largest study to date investigating brain phospholipid and energy metabolism in AD by ^31^P MRS. Increased PCr levels were found in regions that show early degeneration in AD, but not in the ACC, a region known to be involved in this disease at a later stage. Together with an increased pH, this indicates that energy metabolism is altered in mild AD. Although most phospholipid metabolites, like the high energy phosphates, differed between the investigated brain regions in both groups, no alterations in phospholipid metabolism were found specific to mild AD.

Appendix A. Supplementary data {#s0105}
==============================

Supplementary materialImage 1

Acknowledgements {#s0090}
================

We thank all participants and their caregivers for participating in this study. The authors are indebted to Jack van Asten, Miriam Lagemaat and Bart Philips (department of Radiology and Nuclear Medicine, Radboud university medical center) for their help in developing the acquisition and post-processing protocols and Rianne de Heus (department of Geriatric Medicine, Radboud university medical center) for her assistance during data acquisition and analysis. We thank the Institute of Nuclear Sciences Applied to Health and the Visual Neurosciences Lab, Institute for Biomedical Imaging and Life Sciences (University of Coimbra, Coimbra, Portugal) for providing us with spectroscopic co-registry software.

Funding {#s0095}
=======

This work was supported by Nutricia Research. They had no role in the analysis and interpretation of the data, or in writing this report.

Declaration of interest {#s0100}
=======================

Conflicts of interest: none

Supplementary data to this article can be found online at <https://doi.org/10.1016/j.nicl.2018.01.031>.

[^1]: Shared last author.

[^2]: MMSE, mini mental state examination; AD, Alzheimer\'s disease; BMI, body mass index. Data are presented as mean (standard deviation) or median \[range\] unless otherwise indicated.

[^3]: LS means ± SEM of metabolite levels, as percentage of total phosphorus signal in the corresponding spectrum, in patients with Alzheimer\'s disease and healthy control subjects. Significant group effects (p \< 0.05, AD versus Control) indicated in bold. \*PCr also has a significant interaction between group and brain region. LS, least squares; SEM, standard error of the mean; n, number of observations; AD, Alzheimer\'s disease; HC, healthy control; PCr, Phosphocreatine; Pi, inorganic phosphate; ATP, adenosine triphosphate; NAD(H), nicotinamide adenine dinucleotide; PEth, phosphoethanolamine; PCh, phosphocholine; GPEth, glycerophosphoethanolamine; GPCh, glycerophosphocholine.

[^4]: LS means ± SEM of metabolite levels, as percentage of total phosphorus signal in the corresponding spectrum, in anterior cingulate cortex, retrosplenial cortex and hippocampi. Significant (p \< 0.05) main effect of brain region indicated in bold. \*PCr also has a significant interaction between group and brain region. *p* \< 0.05, ACC vs †RSC, ^‡^HL, or ^§^HR; RSC vs ^¶^HL, or ^\#^HR. LS, least squares; SEM, standard error of the mean; n, number of observations; ACC, anterior cingulate cortex; RSC, retrosplenial cortex; HL, left hippocampus; HR, right hippocampus; PCr, Phosphocreatine; Pi, inorganic phosphate; ATP, adenosine triphosphate; NAD(H), nicotinamide adenine dinucleotide; PEth, phosphoethanolamine; PCh, phosphocholine; GPEth, glycerophosphoethanolamine; GPCh, glycerophosphocholine.
