Abstract
Hippocampal subfields have different vulnerability to the degenerative processes related to aging, amnestic mild cognitive impairment (MCI) and Alzheimer’s disease (AD), but the temporal evolution in Parkinson’s disease (PD) is unknown. The purposes of the current work are to describe regional hippocampal changes over time in a sample of PD patients classified according to their baseline cognitive status and to relate these changes to verbal memory loss. T1-weighted images and verbal memory assessment were obtained at two separate time points (3.8 ± 0.4 years apart) from 28 PD with normal cognition (PD-NC), 16 PD with MCI (PD-MCI) and 21 healthy controls (HCs). FreeSurfer 6.0 automated pipeline was used to segment the hippocampus into 12 bilateral subregions. Memory functions were measured with Rey’s Auditory Verbal learning test (RAVLT). We found significant reductions in cornu ammonis 1 (CA1) over time in controls as well as in PD subgroups. Right whole-hippocampal volumes showed time effects in both PD groups but not in controls. PD-NC patients also displayed time effects in the left hippocampal tail and right parasubiculum. Regression analyses showed that specific hippocampal subfield volumes at time 1 predicted almost 60% of the variability in RAVLT delayed-recall score decline. Changes in several hippocampal subregions also showed predictive value for memory loss. In conclusion, CA1 changes in PD were similar to those that occur in normal aging, but PD patients also had more decline in both anterior and posterior hippocampal segments with a more pronounced atrophy of the right hemisphere. Hippocampal segments are better predictors of changes in memory performance than whole-hippocampal volumes.
Introduction
Hippocampal atrophy is a key finding in neurodegenerative diseases (Camicioli et al., ; Small et al., 2011; Bartsch and Wulff, ; Yang and Yu, 2017), although it is also present in healthy aging (Fjell et al., ). In neuroimaging studies, the hippocampus has traditionally been assessed as a single component, but more advanced techniques have allowed studying the hippocampus as a complex structure with specific regional vulnerability to aging and subtypes of dementia (Small et al., 2011).
Extensive previous literature consistently reports region 1 of the cornu ammonis (area CA1), the subiculum (Mueller et al., ) and area CA3 (Pereira et al., ; Wisse et al., 2014) as the regions that are most vulnerable to degeneration in normal aging and Alzheimer’s disease (AD; de Flores et al., ). In Parkinson’s disease (PD), hippocampal atrophy has been associated with dementia (Junqué et al., ; Summerfield et al., 2005; Ibarretxe-Bilbao et al., ), although volume reductions can also be detectable in non-demented PD (Junqué et al., ; Pereira et al., ) and even in unmedicated patients (Noh et al., ).
The detection of regional hippocampal atrophy and its association with memory decline is of high interest in PD since memory impairment has been described as a risk factor for dementia (Levy et al., ). Total hippocampal volumes correlated with learning tasks (Pereira et al., ); recognition memory, on the other hand, has been associated with left hippocampal atrophy (Camicioli et al., ). More recently, volume reductions in some subregions such as areas CA2–3 and CA4 and the dentate gyrus (DG) have been linked to verbal learning impairment in PD (Engvig et al., ; Pereira et al., ). Moreover, CA2–3 atrophy has been found to discriminate healthy controls (HCs) from amnestic mild cognitive impairment (MCI) patients better than global hippocampal volumes (Hanseeuw et al., ).
In the last 3 years, thanks to the development of automated segmentations tools, it has become possible to divide the hippocampus into 12 bilateral segments based on a statistical atlas built upon ultra-high resolution ex-vivo MRI data (Iglesias et al., ). To our knowledge, only one published study investigated differences in percentage change over a 1.5-year follow-up between PD patients with normal cognition (PD-NC) and with MCI (PD-MCI) in hippocampal subfields also using this automated segmentation pipeline (Foo et al., ). However, because it did not include a HC group, this study could not distinguish hippocampal atrophy due to normal aging from that due to PD degeneration.
The aims of the present study were: (1) to investigate longitudinal changes in hippocampal segments in a sample of PD patients classified according to their baseline cognitive status over a 4-year follow-up; (2) to examine the predictive utility of specific hippocampal subfield volumes as well as total volumes at time 1 to determine changes in memory test scores over time in the PD subgroups; and (3) to investigate the relationship between hippocampal changes over time and memory performance decline.
Based on the previous literature on aging, we would expect that CA1 would be one of the segments atrophied over time, but we would also expect to observe changes in other subfields more specific of PD such as CA2–3. We also hypothesized that the changes in total hippocampal volumes as well as specific segments would explain progressive memory decline.
Materials and Methods
Participants
Forty-four PD patients (PD-NC = 28; PD-MCI = 16) from the PD and Movement Disorders Unit, Hospital Clinic (Barcelona, Spain) and 21 HC from the Aging Institute in Barcelona were assessed twice at an interval of 3.8 ± 0.4 years (range: 3.1–5.3).
At time 1, 90 PD patients and 32 HC were recruited between October 2010 and March 2012. Detailed information of the sample can be found in our previous work (Segura et al., 2014). In the present study, only subjects who underwent comprehensive neuropsychological and MRI acquisition at both times were included.
At time 2, two patients underwent deep brain stimulation, five patients and one HC died, 12 PD patients and two controls refused to participate or had moved at follow-up, three PD patients and three controls had developed neurological/psychiatric comorbidities, 15 PD patients had functional impairment and reduced mobility that prevented going to the hospital for MRI scanning, six patients and three HC had MRI motion artifacts or could not finish the scanning protocol and three patients and two HC were excluded due to problems in longitudinal image preprocessing.
Inclusion criteria for patients at time 1 were: (i) fulfilling the UK PD Society Brain Bank diagnostic criteria for PD (Hughes et al., ); and (ii) no surgical treatment with deep-brain stimulation. Exclusion criteria for PD patients and HC were: (i) dementia according to the Movement Disorders Society (MDS) criteria (Emre et al., ) and to clinical assessment performed by a clinical neurologist (MM, FV, YC); (ii) red flags for atypical parkinsonisms; (iii) Hoehn and Yahr (H&Y) scale (Hoehn and Yahr, ) score >3; (iv) young-onset PD; (v) age below 50 years; (vi) presence of severe psychiatric or neurological comorbidity; (vii) low global intellectual quotient estimated by the Vocabulary subtest of the Wechsler Adult Intelligence Scale (scalar score ≤ 7); (viii) Mini Mental State Examination (MMSE) score (Folstein et al., ) below 25; (ix) claustrophobia; (x) pathological MRI findings other than mild white matter hyperintensities in the FLAIR sequence; and (xi) MRI artifacts. At time 2, a diagnosis of dementia, H&Y score >3 and MMSE scores below 25 were not considered as exclusion criteria.
Motor symptoms were assessed with the Unified PD Rating Scale motor section (UPDRS-III, Fahn and Elton, ). All PD patients were taking antiparkinsonian drugs, consisting of different combinations of L-DOPA, catechol-O-methyltransferase inhibitors, monoamine oxidase inhibitors, dopamine agonists and amantadine. In order to standardize doses, the L-DOPA equivalent daily dose (LEDD) was calculated (Tomlinson et al., 2010).
Written informed consent was obtained from all study participants after full explanation of the procedures. The study was approved by the institutional Ethics Committee from the University of Barcelona (IRB00003099).
Neuropsychological and Clinical Assessment
The diagnosis of PD-MCI was established in line with MDS task force recommendations (Litvan et al., ) as previously described in Segura et al. (2014). The memory domain was assessed with Rey’s Auditory Verbal learning test (RAVLT; Lezak et al., ) using total learning (RAVLT total), delayed recall (RAVLT recall) and recognition (RAVLT recognition) scores. Initially, z-scores for each test and for each subject were calculated based on the control group’s means and standard deviations (SDs) from time 1. Expected z-scores adjusted for age, sex and education for each test and each subject were calculated based on a multiple regression analysis performed in the HC group (Aarsland et al., ).
Neuropsychiatric symptoms were evaluated with the Beck Depression Inventory-II (Beck et al., ), Starkstein’s Apathy Scale (Starkstein et al., 1992) and Cumming’s Neuropsychiatric Inventory (Cummings et al., ).
Preprocessing of Longitudinal Imaging
MRI data were acquired with a 3T scanner (MAGNETOM Trio, Siemens, Germany) at both times. The scanning protocol included high-resolution 3-dimensional T1-weighted images acquired in the sagittal plane (TR = 2,300 ms, TE = 2.98 ms, TI = 900 ms, 240 slices, FOV = 256 mm; 1 mm isotropic voxel) and an axial FLAIR sequence (TR = 9,000 ms, TE = 96 ms).
Cross sectional preprocessing of both times was estimated using the automated FreeSurfer stream (version 5.11). Detailed description of FreeSurfer procedures is reported in Segura et al. (2014). In addition, to extract reliable volume and thickness estimates, images were automatically processed with FreeSurfer’s longitudinal stream (Reuter et al., ). Specifically, an unbiased within-subject template space and image is created using robust, inverse consistent registration (Reuter et al., ). Several processing steps, such as skull stripping, Talairach transforms, atlas registration as well as spherical surface maps and parcellations are then initialized with common information from the within-subject template, significantly increasing reliability and statistical power (Reuter et al., ).
After longitudinal preprocessing, FreeSurfer version 6.0 was used to segment the hippocampal subfields2. For a visual representation of the hippocampal segments, see Figure 1.
Figure 1
Ratios were calculated for all hippocampal segment volumes to global hippocampal volumes ((lh or rh segments/lh or rh hippocampus)*100). Global hippocampal to estimated total intracranial volume ratios (eTIV, (lh or rh hippocampus/eTIV)*100) were also calculated.
Statistical Analysis
Cross-Sectional Analyses
Group differences in demographic variables and disease outcomes were analyzed with Kruskal-Wallis tests followed by Mann-Whitney-Wilcox’s pairwise comparisons and Bonferroni correction for quantitative measures. Pearson’s chi-squared test was used where appropriate for categorical measures. These analyses were conducted using RStudio Version 1.1.419 (RStudio Team, ); information on the libraries and functions can be found in Supplementary Methods S1.
A general linear model and Monte Carlo permutation testing with 10,000 iterations were applied to perform group comparisons of hippocampal volumes ratios at time 1 using Matlab R2017a (The MathWorks Inc., Natick, MA, USA). To control type-I errors, a Bonferroni correction was applied. Age and years of education were included as covariates of no interest.
Repeated Measures Analyses
Repeated measures analyses were also conducted with Matlab as described above. Main effects of time and group-by-time interaction were tested on clinical variables such as UPDRS, LEDD and neuropsychiatric symptoms, on memory performance scores and hippocampal subfield volumes between PD groups and HC. Age at time 1 and years of education were used as covariates of no interest in longitudinal hippocampal subfield analyses. For repeated memory score analyses, we used the z-scores adjusted for age, education and sex as described above. Bonferroni correction was applied to all analyses.
Multiple Regression Analyses in the PD Patient Sample
Two different multiple linear regression analyses were performed using two models. As a response variable, both models included the difference between time 2 minus time 1 RAVLT raw scores in total learning, recall and recognition. The first model included age at time 1, years of education and hippocampal segments as predictors. The second model included age at time 1, years of education and whole hippocampal volumes as predictors.
First, we assessed the predictive utility of hippocampal ratios at time 1 to explain the variability in memory performance changes. Second, we included the change in hippocampal ratios (time 2 − time 1) as explanatory variables of memory change.
A stepwise model selection by Akaike information criterion (AIC) was applied on the multiple linear regression models described above. This method picks the best-fitted model that most adequately describes an unknown, high dimensional reality (Zhang, 2016). Resulting hippocampal structures that best described prediction of changes in memory performance can be found in Supplementary Methods S2.
Finally, only multiple regression models with statistical significance are reported. Within the RAVLT recall models, an ANOVA was used to test if there were significant differences between the segments model and the global volumes model.
Results
Demographic Characteristics and Clinical Evolution
There were no significant differences in scan interval between groups (H = 0.013; P = 0.994). Thus, involution of the hippocampus and memory can be directly compared. Moreover, the groups had similar disease duration and H&Y staging scores. Although not significant at p < 0.05, subjects in the HC group were older than those in the PD subgroups; for this reason, age was included as a covariate in group analyses and as a variable of interest in multiple regression models as described in the “Materials and Methods” section (Table 1).
Table 1
| PD-NC n = 28 | PD-MCI n = 16 | PD whole sample n = 44 | Controls n = 21 | Test stats | P-value | |
|---|---|---|---|---|---|---|
| Age, years, median (IQ range) | ||||||
| Time 1 | 59.0 (11.8) | 61.5 (17.3) | 60.0 (10.5) | 67.0 (13.0) | H = 6.2261 U = 593.02 | 0.0441 0.0672 |
| Time 2 | 63.5 (12.5) | 66.0 (15.8) | 64.5 (12.0) | 70.0 (12.0) | H = 5.7951 U = 590.02 | 0.0551 0.0732 |
| Education, years, median (IQ range) | 13.0 (8.3) | 10.0 (7.5) | 12.0 (8.3) | 10.0 (7.0) | H = 1.5831 U = 405.02 | 0.4531 0.4262 |
| Sex, female, n (%) | 8 (28.6) | 6 (37.5) | 14 (31.8) | 10 (47.6) | X2 = 1.8721 X2 = 0.9212 | 0.3921 0.3372 |
| Disease duration, years, median (IQ range) | ||||||
| Time 1 | 6.0 (5.0) | 6.0 (9.5) | 6.0 (7.5) | NA | U = 197.5 | 0.759 |
| Time 2 | 9.0 (4.0) | 8.0 (8.3) | 9.0 (6.3) | NA | U = 248.0 | 0.565 |
| Age of onset, years, median (IQ range) | 53.5 (15.8) | 55.3 (15.8) | 54.0 (15.8) | NA | U = 180.0 | 0.294 |
| Hoehn & Yahr stage, n 1/1.5/2/2.5/3/4 | ||||||
| Time 1 | 11/1/12/2/2/0 | 6/0/8/1/1/0 | 17/1/20/3/3/0 | NA | X2 = 0.718 | 0.949 |
| Time 2 | 3/0/14/0/11/0 | 3/0/5/0/7/1 | 6/0/19/0/18/1 | NA | X2 = 3.111 | 0.375 |
Demographic and clinical characteristics of the sample at both times.
IQ range, interquartile range; NA, not applicable; PD-MCI, Parkinson’s disease with mild cognitive impairment; PD-NC, Parkinson’s disease with normal cognition. P-values are from Kruskal-Wallis test followed by Mann-Whitney pairwise test and Bonferroni correction for continuous variables and chi-squared test for categorical variables. 1Test stats and P-values of PD subgroups and HC comparisons. 2Test stats and P-values of the PD collapsed sample and HC comparisons.
Table 2 summarizes the time effects observed for the clinical measures. The collapsed PD sample had significant decline over time in global cognition scores, increased motor severity as measured by the UPDRS-III and increased neuropsychiatric symptoms. All the PD groups had increased neuropsychiatric symptoms and PD-NC also showed increased depression scores. No significant changes were seen in apathy scale scores.
Table 2
| PD-NC n = 28 | PD-MCI n = 16 | PD all sample n = 44 | Controls n = 21 | Time effects | ||||
|---|---|---|---|---|---|---|---|---|
| PD-NC | PD-MCI | PD all sample | Controls | |||||
| Mini Mental State Examination, mean (SD) | ||||||||
| Time 1 | 29.5 (0.7) | 28.7 (1.5) | 29.2 (1.1) | 29.8 (0.4) | t = 1.366; P = 0.055 | t = 2.957; P = 0.050 | t = 2.843; P = 0.007* | t = 1.291; P = 0.052 |
| Time 2 | 29.1 (1.0) | 27.6 (3.4) | 28.6 (2.3) | 29.3 (0.9) | ||||
| UPDRS part III, mean (SD) | ||||||||
| Time 1 | 13.9 (9.2) | 11.8 (11.0) | 13.1 (9.8) | NA | t = 1.533; P = 0.052 | NS | t = 2.073; P = 0.022* | NA |
| Time 2 | 17.8 (9.0) | 17.4 (12.7) | 17.7 (10.3) | NA | ||||
| LEDD, mg, mean (SD) | ||||||||
| Time 1 | 700.8 (470.6) | 675.6 (535.2) | 691.6 (489.0) | NA | NS | NS | NS | NA |
| Time 2 | 720.5 (388.2) | 693.3 (481.1) | 710.6 (419.0) | NA | ||||
| Beck Depression Inventory II, mean (SD) | ||||||||
| Time 1 | 7.3 (4.9) | 11.1 (6.3) | 8.7 (5.7) | 7.0 (5.4) | t = 1.508; P = 0.043* | NS | NS | NS |
| Time 2 | 8.5 (7.4) | 10.0 (6.3) | 9.1 (7.0) | 5.3 (4.7) | ||||
| Starkstein’s Apathy Scale, mean (SD) | ||||||||
| Time 1 | 10.6 (5.7) | 14.3 (8.7) | 11.9 (7.0) | 9.1 (5.6) | NS | NS | NS | NS |
| Time 2 | 11.4 (6.6) | 14.1 (9.4) | 12.4 (7.7) | 9.3 (5.6) | ||||
| Cummings’ Neuropsychiatric Inventory, mean (SD) | ||||||||
| Time 1 | 4.5 (7.2) | 9.1 (11.8) | 6.2 (9.3) | 1.9 (3.6) | t = 1.993; P = 0.041* | t = 2.343; P = 0.041* | t = 3.016; P = 0.006* | NS |
| Time 2 | 7.1 (8.8) | 13.1 (11.6) | 9.3 (10.2) | 2.3 (2.6) | ||||
Repeated clinical measures.
LEDD, L-DOPA equivalent daily dose; NA, not applicable; PD-MCI, Parkinson’s disease with mild cognitive impairment; PD-NC, Parkinson’s disease with normal cognition; SD, standard deviation; UPDRS part III, Unified Parkinson’s Disease Rating Scale motor section. NS, non-significant. *P-values < 0.05 after Bonferroni correction.
Longitudinal Changes in Hippocampal Segments
Longitudinally, both PD-NC and PD-MCI as well as the PD collapsed sample showed a significant time effect in the right whole hippocampus. Regarding time effects in hippocampal segments, the right CA1 displayed a significant effect of time in all groups of PD patients and HC. Moreover, the left hippocampal tail (HP_tail) and right parasubiculum had significant decreases in the PD collapsed sample and PD-NC patients. Significant group-by-time interaction was seen in the right parasubiculum in the contrast HC vs. PD-NC and PD-NC vs. PD-MCI. Means and SDs of hippocampal segments can be found in Supplementary Table S1; test stats and uncorrected P-values can be found in Table 3. After Bonferroni correction, P-values were not significant.
Table 3
| Time effects | Interaction time × group | |||||||
|---|---|---|---|---|---|---|---|---|
| PD-NC | PD-MCI | PD all sample | Controls | HC > PD-NC | HC > PD-MCI | PD-NC > PD-MCI | HC > PD all sample | |
| Left hippocampal tail | t = 1.923; P = 0.019 | NS | t = 1.764; P = 0.028 | NS | ||||
| Right CA1 | t = 1.682; P = 0.030 | t = 2.082; P = 0.019 | t = 1.864 P = 0.027 | t = 2.037 P = 0.023 | ||||
| Right parasubiculum | t = 1.985 P = 0.035 | NS | t = 1.932 P = 0.026 | NS | t = 1.498; P = 0.031 | t = −2.574 P = 0.020 | ||
| Right whole hippocampus | t = 1.614; P = 0.031 | t = 1.539; P = 0.040 | t = 1.626; P = 0.030 | NS | ||||
Time effects and group-by-time interaction for significant hippocampal segments.
CA, Cornu Ammonis; HC, healthy controls; NS, non-significant; PD-MCI, Parkinson’s disease mild cognitive impairment; PD-NC, Parkinson’s disease normal cognition.
Predictive Utility of Hippocampal Volumes in Memory Decline
Memory Decline
Regarding memory performance, all z-scores were lower at follow-up. The collapsed PD sample showed significant decline in all variables mainly due to progressive impairment in the PD-NC group. PD-NC had a significant decrease in RAVLT total learning and recognition (Table 4). For RAVLT total learning, there was a significant group-by-time interaction between PD-NC and HC (t = 2.301; P = 0.013; P-corrected < 0.05). For RAVLT recognition, the interaction was significant for PD-NC and HC (t = 2.969; P < 0.001; P-corrected < 0.05) and for all PD sample vs. controls (t = 2.713; P < 0.001; P-corrected = 0.05).
Table 4
| PD-NC n = 28 mean (SD) | PD-MCI n = 16 mean (SD) | PD all sample n = 44 mean (SD) | Controls n = 21 mean (SD) | Time effects | ||||
|---|---|---|---|---|---|---|---|---|
| PD-NC | PD-MCI | PD all sample | Controls | |||||
| RAVLT total learning | ||||||||
| Time 1 | 0.14 (1.12) | −1.52 (0.80) | −0.47 (1.29) | 0.00 (0.81) | t = 3.836; P = 0.001* | NS | t = 3.315; P = 0.001* | NS |
| Time 2 | −0.58 (1.30) | −1.66 (1.25) | −0.97 (1.37) | −0.06 (0.73) | ||||
| RAVLT delayed recall | ||||||||
| Time 1 | −0.15 (1.12) | −1.30 (0.94) | −0.57 (1.19) | 0.06 (0.90) | NS | NS | t = 1.956; P = 0.034 | NS |
| Time 2 | −0.46 (1.38) | −1.70 (1.10) | −0.91 (1.41) | −0.04 (0.74) | ||||
| RAVLT recognition | ||||||||
| Time 1 | 0.20 (1.04) | −0.95 (2.23) | −0.22 (1.65) | −0.17 (0.76) | t = 4.059; P < 0.001* | NS | t = 4.178; P < 0.001* | NS |
| Time 2 | −1.38 (2.85) | −1.76 (1.86) | −1.51 (2.52) | 0.01 (0.94) | ||||
Repeated measures analysis of memory performance.
NS, not significant; PD-MCI, Parkinson’s disease with mild cognitive impairment; PD-NC, Parkinson’s disease with normal cognition; RAVLT, Rey’s auditory verbal learning test. Means are adjusted z-scores by age, years of education and sex. *Indicates statistical significance (P-values < 0.05) after Bonferroni correction.
Hippocampal Volume Ratios at Time 1 as Predictors of Memory Change Over Time
The first multiple regression approach investigated whether hippocampal volume ratios at time 1 can be good predictors of memory performance change.
For RAVLT total changes over time, using whole hippocampal volume ratios as predictors, the right whole hippocampus was a significant predictive variable (R2 = 0.16; adjusted R2 = 0.14; F = 8.074; P = 0.007). However, the segments model for change over time was not significant (R2 = 0.26; adjusted R2 = 0.07; F = 1.332; P = 0.257).
For RAVLT recall change, the significant predictive variables were age, left CA3, right CA4, left parasubiculum, right subiculum, right fimbria, right HP_tail, right fissure and left molecular layer (R2 = 0.69; adjusted R2 = 0.56; F = 5.634; P < 0.001). Considering that the bilateral whole hippocampus was a significant predictor (R2 = 0.23; adjusted R2 = 0.19; F = 5.960; P = 0.005), there were significant differences between the segments and the global volumes models (F = 4.539; P = 0.001).
Regarding RAVLT recognition change, the bilateral CA1, right CA4, bilateral subiculum and right Hippocampal Amygdala Transition Area (HATA) were significant predictive variables (R2 = 0.44; adjusted R2 = 0.27; F = 2.557; P = 0.021).
Detailed information of the multiple regression models for each test can be found in Supplementary Table S2.
Relationship Between Hippocampal Volume Ratio Change and Memory Decline
The second multiple regression approach aimed to investigate whether changes in hippocampal volume ratios can explain changes in memory performance.
Changes in right fimbria, right HP_tail and left fissure were significant explanatory variables of changes in RAVLT recall scores (R2 = 0.67; adjusted R2 = 0.45; F = 3.093; P = 0.005). In the global model, the left hippocampus was the only significant variable (R2 = 0.20; adjusted R2 = 0.14; F = 3.240; P = 0.032). There were significant differences between the two models (F = 2.659; P = 0.015).
Finally, for RAVLT recognition, when considering the whole structure (R2 = 0.14; adjusted R2 = 0.12; F = 6.589; P = 0.014), the left hippocampus was a significant predictor.
Information regarding the multiple regression models for each test can be found in Supplementary Table S3.
Discussion
The main findings of the present study were: (1) the right CA1 was sensitive to time effects in normal aging and in PD with NC and with MCI; (2) volume decrements in right whole hippocampus volume as well as specific regional volumes were only found in PD; and (3) hippocampal subfields were better predictors of delayed verbal memory recall decline than global hippocampal volumes.
The right CA1 showed a significant time effect for all PD groups and for HC. No significant group-by-time interactions were found. Therefore, the changes observed seem to be due to aging effects rather than specific of PD. CA1 has been reported as one of the regions with the earliest and strongest involvement over time in AD (Small et al., 2011), being useful to discriminate healthy subjects from those with MCI (Mueller et al., ) and to predict conversion from MCI to AD (Apostolova et al., ). Indeed, early neuropathological studies have described a high susceptibility to the accumulation of amyloid-β in CA1 both in a mouse model and humans (Furcila et al., ).
Cross-sectional studies have reported that the head of the hippocampus is the most vulnerable region in normal aging (Ta et al., 2012), in non-demented PD (Ibarretxe-Bilbao et al., ) and in demented PD patients (Bouchard et al., ; Ibarretxe-Bilbao et al., ) but demented patients also revealed posterior hippocampal atrophy (Ibarretxe-Bilbao et al., ). In the present longitudinal study, in addition to the decrements described above, we also found significant volume decrements in the left HP_tail for the PD-NC group, suggesting that specific posterior hippocampal atrophy takes place at earlier stages of the disease.
When considering global hippocampal volumes, the right hippocampus had specific time effects in all PD subgroups. This pattern is different from what occurs in amnestic MCI and AD. In a meta-analysis of 14 studies it has been reported that although in both MCI and AD there are progressive bilateral reductions, the effect size is greater for the left hemisphere when compared with the right (Shi et al., 2009). A recent work by Yue et al. (2018) reported hippocampal asymmetry in MCI patients and individuals with subjective cognitive decline compared with HC where the left hemisphere was more atrophic than the right.
The right parasubiculum was also sensitive to time effects in the PD-NC group. Volume decrements in this region were significantly higher than those observed in the HC and the PD-MCI group as demonstrated by the significant group-by-time interaction. The parasubiculum is a small hippocampal structure that is usually studied together with the subiculum and parasubiculum. Therefore, there is not much previous literature using MRI techniques describing the implication of this structure in aging or neurodegenerative processes. However, we could speculate that the volume decrements found in the right parasubiculum might be related to parietal atrophy through transneuronal degeneration. The parasubiculum has direct projections to medial parieto-temporal regions involved in visuospatial processes (Dalton and Maguire, ). In PD, there is a structural temporo-parietal atrophy suggested as a marker of cognitive decline (Segura et al., 2014). More specifically, medial parietal atrophy has been linked to visuospatial impairment in PD patients (Garcia-Diaz et al., ).
Regarding memory performance, hippocampal segment volumes as well as whole hippocampal volumes at baseline have been reported as significant predictors of verbal memory delayed recall changes (Beyer et al., ). In our study, the segments model explained 56% of the variability, whereas the whole hippocampal volumes model only explained 19%. In line with this, the right fimbria, right HP_tail and left fissure volume changes over time were linked to changes in RAVLT recall changes in a model that explained almost 50% of the variance. By contrast, the whole hippocampal model was more useful to predict changes in RAVLT total learning scores, although these models did not explain much variability. We could speculate that declines in verbal memory learning would be related to hippocampal-neocortical connectivity more than based on structural changes in the hippocampal formation per se (Fjell et al., ). Finally, models including RAVLT recognition changes also explained less than 50% of the variance, with the segments model explaining more variance than the global volumes model.
The strengths of the present study are: the inclusion of a control group allowed us to compare hippocampal atrophy in PD with atrophy that occurs in elderly healthy subjects as part of the aging process. Also, the use of novel neuroimaging automated pipelines to accurately segment the hippocampus to investigate regional vulnerability. FreeSurfer’s pipeline is based on ex vivo 7T images to manually segment the hippocampus in order to create the statistical atlas; and it has been recently proved to have a good test-retest reliability over time (Worker et al., 2018). However, the authors (Iglesias et al., ) recommend caution on the interpretation of results involving the internal subfields such as the CA4, molecular layer or the Granule cells in the molecular layer of the DG. To the best of our knowledge, there is no study that compares the subfields overlap with other manual or automated segmentations methods.
The most important limitation would be the small sample size due to a high attrition rate in the Parkinson’s cohort. To overcome this problem, which affects generalization of the results, larger multicentric studies in PD should help clarify the progressive pattern of degeneration in the hippocampus. There are few longitudinal studies in PD cohorts performing MRI assessments over more than 1.5 years of follow-up. The work of Ulla et al. (2013) followed a cohort of PD patients over 3 years and they also reported an attrition rate of 50%. Multicentric longitudinal initiatives are particularly common in AD, such as the ADNI database3 or the AIBL initiative. In the longitudinal AIBL cohort, almost 60% of the participants returned to MRI and PET follow-up scans (Doré et al., ). This relatively higher percentage could be explained because elder controls and AD patients have less motor impairment than PD patients. Moreover, we would like to highlight as a frequent limitation of longitudinal studies that participants with worse disease prognosis, with more depressive/apathetic symptoms or with greater functional impairment in their daily living are more likely to be lost to follow-up.
It could be also mentioned that, due to the exploratory nature of the hippocampal subfield study, we report uncorrected P-values. These results should thus be interpreted with caution. This limitation is common to all studies using these new hippocampal segmentations; however, exploratory analyses are necessary to progress in neuroimaging research.
In conclusion, besides regional vulnerability in hippocampus degeneration dependent in part of aging, we found specific hippocampal regions that were more sensitive to time effects in PD. The right global hippocampus also seems to be more vulnerable than the left. Finally, specific hippocampal segment volumes were found to be good markers of verbal delayed recall performance decline over time.
Statements
Author contributions
CJ contributed in the design of the study. CU, AA and AC contributed to the analysis of the data and CU, BS, HB, AA, AC, MM, FV, YC and CJ contributed to the interpretation of the data. CU and CJ contributed to the draft of the article. CU, BS, HB, AA, AC, MM, FV, YC, NB and CJ revised the manuscript critically for important intellectual content and approved the final version of the manuscript.
Funding
CU was supported by a fellowship from 2014, Spanish Ministry of Economy and Competitiveness (Ministerio de Economía y Competitividad; BES-2014-068173) and cofinanced by the European Social Fund (ESF), AC was supported by APIF predoctoral fellowship from the University of Barcelona (2017–2018) and AA was supported by a 2016 fellowship from the Departament d’Empresa i Coneixement de la Generalitat de Catalunya, Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR; 2016FI_B 00360). This study was sponsored by the Spanish Ministry of Economy and Competitiveness (PSI2013-41393-P; PSI2017-86930-P cofinanced by Agencia Estatal de Investigación (AEI) and the European Regional Development Fund), by Generalitat de Catalunya (2017SGR 98) and by Fundació La Marató de TV3 in Spain (20142310).
Acknowledgments
We thank the cooperation of the patients, their families and control subjects. We are also indebted to the Magnetic Resonance Imaging core facility of the IDIBAPS for the technical support, especially to C. Garrido, G. Lasso, V. Sanchez and A. Albaladejo; and we would also like to acknowledge the CERCA Programme/Generalitat de Catalunya.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2018.00325/full#supplementary-material
Footnotes
1.^https://surfer.nmr.mgh.harvard.edu
2.^https://surfer.nmr.mgh.harvard.edu/fswiki/LongitudinalHippocampalSubfields
References
1
AarslandD.BrønnickK.LarsenJ. P.TysnesO. B.AlvesG. (2009). Cognitive impairment in incident, untreated parkinson disease: the norwegian parkwest study. Neurology72, 1121–1126. 10.1212/01.wnl.0000338632.00552.cb
2
ApostolovaL. G.ThompsonP. M.GreenA. E.HwangK. S.ZoumalanC.JackC. R.et al. (2010). 3D comparison of low, intermediate, and advanced hippocampal atrophy in MCI. Hum. Brain Mapp.31, 786–797. 10.1002/hbm.20905
3
BartschT.WulffP. (2015). The hippocampus in aging and disease: from plasticity to vulnerability. Neuroscience309, 1–16. 10.1016/j.neuroscience.2015.07.084
4
BeckA. T.SteerR. A.BrownG. K. (1996). Manual for the Beck Depression Inventory.San Antonio, TX: Psychological Corporation.
5
BeyerM. K.BronnickK. S.HwangK. S.BergslandN.TysnesO. B.LarsenJ. P.et al. (2012). Verbal memory is associated with structural hippocampal changes in newly diagnosed Parkinson’s disease. J. Neurol. Neurosurg. Psychiatry84, 23–28. 10.1136/jnnp-2012-303054
6
BouchardT. P.MalykhinN.MartinW. R. W.HanstockC. C.EmeryD. J.FisherN. J.et al. (2008). Age and dementia-associated atrophy predominates in the hippocampal head and amygdala in Parkinson’s disease. Neurobiol. Aging29, 1027–1039. 10.1016/j.neurobiolaging.2007.02.002
7
CamicioliR.MooreM. M.KinneyA.CorbridgeE.GlassbergK.KayeJ. A. (2003). Parkinson’s disease is associated with hippocampal atrophy. Mov. Disord.18, 784–790. 10.1002/mds.10444
8
CummingsJ. L.MegaM.GrayK.Rosenberg-ThompsonS.CarusiD. A.GornbeinJ. (1994). The neuropsychiatric inventory: comprehensive assessment of psychopathology in dementia. Neurology44, 2308–2308. 10.1212/wnl.44.12.2308
9
DaltonM. A.MaguireE. A. (2017). The pre/parasubiculum: a hippocampal hub for scene-based cognition?Curr. Opin. Behav. Sci.17, 34–40. 10.1016/j.cobeha.2017.06.001
10
de FloresR.La JoieR.ChételatG. (2015). Structural imaging of hippocampal subfields in healthy aging and Alzheimer’s disease. Neuroscience309, 29–50. 10.1016/j.neuroscience.2015.08.033
11
DoréV.VillemagneV. L.BourgeatP.FrippJ.AcostaO.ChetélatG.et al. (2013). Cross-sectional and longitudinal analysis of the relationship between aβ deposition, cortical thickness, and memory in cognitively unimpaired individuals and in Alzheimer disease. JAMA Neurol.70:903. 10.1001/jamaneurol.2013.1062
12
EmreM.AarslandD.BrownR.BurnD. J.DuyckaertsC.MizunoY.et al. (2007). Clinical diagnostic criteria for dementia associated with Parkinson’s disease. Mov. Disord.22, 1689–1707. 10.1002/mds.21507
13
EngvigA.FjellA. M.WestlyeL. T.SkaaneN. V.SundsethØ.WalhovdK. B. (2012). Hippocampal subfield volumes correlate with memory training benefit in subjective memory impairment. Neuroimage61, 188–194. 10.1016/j.neuroimage.2012.02.072
14
FahnS.EltonR. (1987). “UPDRS program member unified Parkinson’s disease rating scale,” in Recent Developments in Parkinson’s Disease, (Vol. 2) eds FahnD.MarsdenS.GoldsteinC. D.CalneM. (Florham Park, NJ: Macmillan Healthcare Information), 153–163.
15
FjellA. M.SneveM. H.StorsveA. B.GrydelandH.YendikiA.WalhovdK. B. (2016). Brain events underlying episodic memory changes in aging: a longitudinal investigation of structural and functional connectivity. Cereb. Cortex26, 1272–1286. 10.1093/cercor/bhv102
16
FjellA. M.WestlyeL. T.GrydelandH.AmlienI.EspesethT.ReinvangI.et al. (2014). Accelerating cortical thinning: unique to dementia or universal in aging?Phys. Rev. E24, 919–934. 10.1093/cercor/bhs379
17
FolsteinM. F.FolsteinS. E.McHughP. R. (1975). “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J. Psychiatr. Res.12, 189–198. 10.1016/0022-3956(75)90026-6
18
FooH.MakE.ChanderR. J.NgA.AuW. L.SitohY. Y.et al. (2016). Associations of hippocampal subfields in the progression of cognitive decline related to Parkinson’s disease. Neuroimage Clin.14, 37–42. 10.1016/j.nicl.2016.12.008
19
FurcilaD.DeFelipeJ.Alonso-nanclaresL. (2018). A study of amyloid-B and phosphotau in plaques and neurons in the hippocampus of Alzheimer’s disease patients. J. Alzheimers Dis.64, 417–435. 10.3233/jad-180173
20
Garcia-DiazA. I.SeguraB.BaggioH. C.MartiM. J.ValldeoriolaF.ComptaY.et al. (2014). Structural MRI correlates of the MMSE and pentagon copying test in Parkinson’s disease. Park. Relat. Disord.20, 1405–1410. 10.1016/j.parkreldis.2014.10.014
21
HanseeuwB. J.Van LeemputK.KavecM.GrandinC.SeronX.IvanoiuA. (2011). Mild cognitive impairment: differential atrophy in the hippocampal subfields. Am. J. Neuroradiol.32, 1658–1661. 10.3174/ajnr.A2589
22
HoehnM. M.YahrM. D. (1967). Parkinsonism: onset, progression and mortality. Neurology17, 427–442. 10.1212/wnl.17.5.427
23
HughesA. J.DanielS. E.KilfordL.LeesA. J. (1992). Accuracy of clinical diagnosis of idiopathic Parkinson’s disease: a clinico-pathological study of 100 cases. J. Neurol. Neurosurg. Psychiatry55, 181–184. 10.1136/jnnp.55.3.181
24
Ibarretxe-BilbaoN.Ramírez-RuizB.TolosaE.MartíM. J.ValldeoriolaF.BargallóN.et al. (2008). Hippocampal head atrophy predominance in Parkinson’s disease with hallucinations and with dementia. J. Neurol.255, 1324–1331. 10.1007/s00415-008-0885-8
25
IglesiasJ. E.AugustinackJ. C.NguyenK.PlayerC. M.PlayerA.WrightM.et al. (2015). A computational atlas of the hippocampal formation using ex vivo, ultra-high resolution MRI: application to adaptive segmentation of in vivo MRI. Neuroimage115, 117–137. 10.1016/j.neuroimage.2015.04.042
26
JunquéC.Ramírez-RuizB.TolosaE.SummerfieldC.MartíM. J.PastorP.et al. (2005). Amygdalar and hippocampal MRI volumetric reductions in Parkinson’s disease with dementia. Mov. Disord.20, 540–544. 10.1002/mds.20371
27
LevyG.JacobsD. M.TangM. X.CôtéL. J.LouisE. D.AlfaroB.et al. (2002). Memory and executive function impairment predict dementia in Parkinson’s disease. Mov. Disord.17, 1221–1226. 10.1002/mds.10280
28
LezakM. D.HowiesonD. B.LoringD. W.HannayH. J.FischerJ. S.BiglerE. D.et al. (2012). Neuropsychological Assessment.New York, NY: Oxford University Press. Available online at: https://global.oup.com/academic/product/neuropsychological-assessment-9780195395525?cc=it&lang=en& [Accessed on July 11, 2018].
29
LitvanI.GoldmanJ. G.TrösterA. I.SchmandB. A.WeintraubD.PetersenR. C.et al. (2012). Diagnostic criteria for mild cognitive impairment in Parkinson’s disease: movement disorder society task force guidelines. Mov. Disord.27, 349–356. 10.1002/mds.24893
30
MuellerS. G.SchuffN.YaffeK.MadisonC.MillerB.WeinerM. W. (2010). Hippocampal atrophy patterns in mild cognitive impairment and Alzheimer’s disease. Hum. Brain Mapp.31, 1339–1347. 10.1002/hbm.20934
31
NohS. W.HanY. H.MunC. W.ChungE. J.KimE. G.JiK. H.et al. (2014). Analysis among cognitive profiles and gray matter volume in newly diagnosed Parkinson’s disease with mild cognitive impairment. J. Neurol. Sci.347, 210–213. 10.1016/j.jns.2014.09.049
32
PereiraJ. B.JunquéC.Bartrés-FazD.Ramírez-RuizB.MartiM. J.TolosaE. (2013). Regional vulnerability of hippocampal subfields and memory deficits in Parkinson’s disease. Hippocampus23, 720–728. 10.1002/hipo.22131
33
PereiraJ. B.Valls-PedretC.RosE.PalaciosE.FalcónC.BargallóN.et al. (2014). Regional vulnerability of hippocampal subfields to aging measured by structural and diffusion MRI. Hippocampus24, 403–414. 10.1002/hipo.22234
34
ReuterM.RosasH. D.FischlB. (2010). Highly accurate inverse consistent registration: a robust approach. Neuroimage53, 1181–1196. 10.1016/j.neuroimage.2010.07.020
35
ReuterM.SchmanskyN. J.RosasH. D.FischlB. (2012). Within-subject template estimation for unbiased longitudinal image analysis. Neuroimage61, 1402–1418. 10.1016/j.neuroimage.2012.02.084
36
RStudio Team. (2015). RStudio: Integrated Development. Available online at: http://www.rstudio.com/
37
SeguraB.BaggioH. C.MartiM. J.ValldeoriolaF.ComptaY.Garcia-DiazA. I.et al. (2014). Cortical thinning associated with mild cognitive impairment in Parkinson’s disease. Mov. Disord.29, 1495–1503. 10.1002/mds.25982
38
ShiF.LiuB.ZhouY.YuC.JiangT. (2009). Hippocampal volume and asymmetry in mild cognitive impairment and Alzheimer’s disease: meta-analyses of MRI studies. Hippocampus19, 1055–1064. 10.1002/hipo.20573
39
SmallS. A.SchobelS. A.BuxtonR. B.WitterM. P.BarnesC. A. (2011). A pathophysiological framework of hippocampal dysfunction in ageing and disease. Nat. Rev. Neurosci.12, 585–601. 10.1038/nrn3085
40
StarksteinS. E.MaybergH. S.PreziosiT. J.AndrezejewskiP.LeiguardaR.RobinsonR. G. (1992). Reliability, validity, and clinical correlates of apathy in Parkinson’s disease. J. Neuropsychiatry Clin. Neurosci.4, 134–139. 10.1176/jnp.4.2.134
41
SummerfieldC.JunquéC.TolosaE.Salgado-PinedaP.Gómez-AnsónB.MartíM. J.et al. (2005). Structural brain changes in parkinson disease with dementia: a voxel-based morphometry study. Arch. Neurol.62, 281–285. 10.1001/archneur.62.2.281
42
TaA. T.HuangS.-E.ChiuM.-J.HuaM.-S.TsengW.-Y. I.ChenS.-H. A.et al. (2012). Age-related vulnerabilities along the hippocampal longitudinal axis. Hum. Brain Mapp.33, 2415–2427. 10.1002/hbm.21364
43
TomlinsonC. L.StoweR.PatelS.RickC.GrayR.ClarkeC. E. (2010). Systematic review of levodopa dose equivalency reporting in Parkinson’s disease. Mov. Disord.25, 2649–2653. 10.1002/mds.23429
44
UllaM.BonnyJ. M.OuchchaneL.RieuI.ClaiseB.DurifF. (2013). Is R2* a new MRI biomarker for the progression of Parkinson’s disease? A longitudinal follow-up. PLoS One8:e57904. 10.1371/journal.pone.0057904
45
WisseL. E. M.BiesselsG. J.HeringaS. M.KuijfH. J.KoekD. L.LuijtenP. R.et al. (2014). Hippocampal subfield volumes at 7T in early Alzheimer’s disease and normal aging. Neurobiol. Aging35, 2039–2045. 10.1016/j.neurobiolaging.2014.02.021
46
WorkerA.DimaD.CombesA.CrumW. R.StrefferJ.EinsteinS.et al. (2018). Test-retest reliability and longitudinal analysis of automated hippocampal subregion volumes in healthy ageing and Alzheimer’s disease populations. Hum. Brain Mapp.39, 1743–1754. 10.1002/hbm.23948
47
YangW.YuS. (2017). Synucleinopathies: common features and hippocampal manifestations. Cell. Mol. Life Sci.74, 1485–1501. 10.1007/s00018-016-2411-y
48
YueL.WangT.WangJ.LiG.WangJ.LiX.et al. (2018). Asymmetry of hippocampus and amygdala defect in subjective cognitive decline among the community dwelling chinese. Front. Psychiatry9:226. 10.3389/fpsyt.2018.00226
49
ZhangZ. (2016). Variable selection with stepwise and best subset approaches. Ann. Transl. Med.4:136. 10.21037/atm.2016.03.35
Summary
Keywords
hippocampal subfields, MRI, Parkinson’s disease, longitudinal assessment, memory impairment, aging
Citation
Uribe C, Segura B, Baggio HC, Campabadal A, Abos A, Compta Y, Marti MJ, Valldeoriola F, Bargallo N and Junque C (2018) Differential Progression of Regional Hippocampal Atrophy in Aging and Parkinson’s Disease. Front. Aging Neurosci. 10:325. doi: 10.3389/fnagi.2018.00325
Received
19 July 2018
Accepted
26 September 2018
Published
11 October 2018
Volume
10 - 2018
Edited by
Daniel Ortuño-Sahagún, Universidad de Guadalajara, Mexico
Reviewed by
Richard Camicioli, University of Alberta, Canada; Liuqing Yang, Johns Hopkins Medicine, United States
Updates
Copyright
© 2018 Uribe, Segura, Baggio, Campabadal, Abos, Compta, Marti, Valldeoriola, Bargallo and Junque.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Carme Junque cjunque@ub.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.