Abstract
Episodic memory declines with advancing age. Neuroimaging studies have associated such decline to age-related changes in general cognitive-control networks as well as to changes in process-specific encoding or retrieval networks. To assess the specific influence of aging on encoding and retrieval processes and associated brain systems, it is vital to dissociate encoding and retrieval from each other and from shared cognitive-control processes. We used multivariate partial-least-squares to analyze functional magnetic resonance imaging data from a large population-based sample (n = 292, 25–80 years). The participants performed a face-name paired-associates task and an active baseline task. The analysis revealed two significant network patterns. The first reflected a process-general encoding-retrieval network that included frontoparietal cortices and posterior hippocampus. The second pattern dissociated encoding and retrieval networks. The anterior hippocampus was differentially engaged during encoding. Brain scores, representing whole-brain integrated measures of how strongly an individual recruited a brain network, were correlated with cognitive performance and chronological age. The scores from the general cognitive-control network correlated negatively with episodic memory performance and positively with age. The encoding brain scores, which strongly reflected hippocampal functioning, correlated positively with episodic memory performance and negatively with age. Univariate analyses confirmed that bilateral hippocampus showed the most pronounced activity reduction in older age, and brain structure analyses found that the activity reduction partly related to hippocampus atrophy. Collectively, these findings suggest that age-related structural brain changes underlie age-related reductions in the efficient recruitment of a process-specific encoding network, which cascades into upregulated recruitment of a general cognitive-control network.
Introduction
Episodic memory (Tulving, 1983) is the most age-sensitive long-term memory system (Nyberg et al., 1996b; Rönnlund et al., 2005). Imaging studies have linked episodic memory impairment in aging to several key regions of a large-scale episodic memory network, including medial temporal lobe (MTL; Grady et al., 1995), parietal (Nyberg et al., 2003), and frontal (Cabeza et al., 1997) regions. However, the process-impure nature of cognitive tasks (Foster and Jelicic, 1999) makes it difficult to say whether observations in imaging studies of age-related changes in functional networks is indicative of a core episodic memory deficit or whether it reflects more general working-memory and/or executive deficits (Gazzaley et al., 2008).
Moreover, it remains unclear whether aging selectively impairs encoding or retrieval processes, or whether it has more process-general mnemonic influences. Some studies have found that certain MTL regions are jointly activated during encoding and retrieval (Nyberg et al., 1996c; Greicius et al., 2003), but there is also evidence for functional specialization (Gabrieli, 1997; Rombouts et al., 2001; Giovanello et al., 2009). Similarly, there is evidence that various episodic encoding and retrieval conditions engage overlapping frontal regions (Nyberg et al., 1996c), but also that specific frontal regions are differentially engaged during encoding versus retrieval (Cabeza and Nyberg, 2000).
To specifically assess the influence of aging on episodic memory encoding and retrieval processes and associated brain systems, it is vital to dissociate encoding and retrieval from each other as well as from shared cognitive-control processes. In the present study, we used a multivariate partial-least-squares (PLS) analysis method (McIntosh and Lobaugh, 2004) to approach this issue. The PLS method has previously been used to identify common as well as process-specific patterns of brain regions involved in episodic encoding and retrieval (Nyberg et al., 1996c). PLS accommodates indices of both similarities and differences by considering all experimental variables, simultaneously. Thus, if encoding and retrieval engage many of the same brain regions, PLS should reveal a common network. Alternatively, or in addition, if some brain regions are differentially engaged during encoding and retrieval, PLS can reveal process-specific encoding and retrieval networks.
PLS was applied to data from a large-scale sample [N = 292, 25–80 years of age (Salami et al., 2012)] from a population-based study [Betula (Nilsson et al., 1997, 2004; Nyberg et al., 2010)]. The participants were scanned with functional magnetic resonance imaging (fMRI) while they engaged in a face-name paired-associates (FN-PA) encoding and retrieval task (Sperling et al., 2001; Zeineh et al., 2003) as well as an active baseline task (Stark and Squire, 2001). They also took part in a detailed cognitive test session that included tests of episodic memory, semantic memory, and visuospatial ability. We correlated individual “brain scores” from the PLS analysis, representing an integrated measure of how strongly an individual recruited a certain brain network, with memory performance and chronological age (Grady et al., 2006). We expected process-specific episodic networks to be positively related to performance on episodic tasks (but not other cognitive tasks) and negatively related to chronological age.
Materials and Methods
Participants.
All participants in the present analyses (N = 292, 152 females) were part of the Betula Prospective Cohort Study on memory, health, and aging. Participants' ages ranged from 25 to 80 years (mean = 59.5 ± 13.1 years). From the initial sample of 376 participants, 84 were excluded. Reasons for exclusion included low task performance (lower than chance level; 62 participants, 73.8%), not performing the task at all (9 participants, 10.7%), technical error in E-prime (1 participant, 1.1%), misunderstanding the task (1 participant, 1.1%), pathology (1 participant, 1.1%), and corrupted T1-weighted structural image (10 participants, 11.9%). All remaining participants were Swedish native speakers, had normal or corrected to normal vision, and had no history of severe neurological illness or events that might cause dementia. The older adults were examined by a physician and were generally healthy with no evidence of severe health problems (e.g., heart or kidney diseases). Informed consent was signed by all participants, in accordance with guidelines of the Swedish Council for Research in the Humanities and Social Sciences.
fMRI activation task.
The FN-PA task was used, which consists of 20 blocks (in total) of encoding, retrieval, and control tasks, alternating during the scanning session. Within each encoding and retrieval block (six blocks each), four stimuli were presented for 4 s each, with a randomized interstimulus interval of 1.5, 2.5, 3, or 4.5 s. Block-specific instructions were provided to participants before each block for 4 s. During the encoding blocks, the face-name stimuli were presented and participants were instructed to remember the name associated with each face. They were also instructed to press a button, using the index finger, to indicate that they saw a face-name pair. The faces (an equal number of males and females) were digital color photographs portraying persons, unknown to the participants. The names were fictional first names, printed in white on a black background, that were paired with the faces by the investigators.
During the retrieval blocks, each face was presented along with three letters of which one corresponded to the first letter of the previously encoded face-name pair. The task was to indicate the letter corresponding to the name that was previously encoded with the face. The top, middle, and bottom letters corresponded to the index, the ring, and the middle fingers, respectively. Participants were instructed to respond by guessing if they could not remember the face-name pair.
An active baseline task was used to dissociate different aspects of episodic memory from sensory motor components. Also, it has been shown that activity in the MTL is substantially higher during rest than during active baseline conditions such as the one adopted in the current study (Stark and Squire, 2001). Given that we expected to capture reliable hippocampal activation relevant to episodic memory functions, it was critical to accommodate a suitable baseline task. During each control block (a total of eight blocks), a fixation cross appeared in the center of the visual field for 1–2.5 s and was briefly replaced with a circle that appeared for 0.5 s. The circle was followed by another fixation cross, which was presented for 2–3.5 s. Participants were instructed to indicate, as quickly as possible, when the fixation cross changed to a circle.
Image acquisition.
fMRI data were collected on a 3 T Discovery MR750 (General Electric) scanner with a 32-channel head coil. Scanner parameters for the gradient-echo EPI sequence were as follows: TR = 2.0 s, 37 slices (3.9 mm thick), 96 × 96 matrix, FOV = 250 mm, flip angle = 80°, TE = 30 ms. Ten dummy scans were collected to allow for equilibration of the fMRI signal. Stimuli were presented on a computer screen that was seen by participants through a tilted mirror attached to the head coil. E-prime (Psychology Software Tools; www.pstnet.com/eprime) was used for presentation of the stimuli and the responses were made on an MR-compatible response pad (Current Designs). Structural high-resolution T1-weighted images (180 slices, 1 mm thickness, FOV = 25 × 25 cm, voxel size = 0.5 × 0.5 × 1 mm3) were collected following the functional images.
Data analysis.
Functional and structural images were preprocessed using statistical parametric mapping software (SPM8; Wellcome Department of Imaging Science, Functional Imaging Laboratory). All functional images were first corrected for acquisition time differences between slices within each volume. The slice-timing corrected images were then rigidly aligned to the first image volume to correct for head movement. A within-subject rigid registration was conducted to align functional and structural images together. T1-weighted images were segmented into gray matter (GM) and white matter (WM) using a new segmentation algorithm implemented in SPM8 [modified version of unified segmentation by (Ashburner and Friston, 2005)]. Then, a group-specific template was created using Diffeomorphic Anatomical Registration Using Exponentiated Lie Algebra (DARTEL) (Ashburner, 2007). This was done by first importing tissue class images (e.g., GM, WM) into the DARTEL space using the normalization parameter yielded during the segmentation step followed by resampling to isotropic voxels (1.5 × 1.5 × 1.5 mm). Then, the imported images went through an iterative procedure that began by producing an initial template as a mean of GM/WM across all participants (n = 292). Deformation from the initial template to each of the subject-specific GM/WM images was computed and the inverse of the deformation was applied to each of the subject-specific GM/WM images. A second template was then created as the mean of the deformed subject-specific GM/WM images across all participants and this procedure was repeated until a sixth template was created. Finally, the coregistered fMRI images and segmented GM/WM images were nonlinearly normalized, subject by subject, to the sample-specific template (using a subject-specific flow field), affine aligned into a Montreal Neurological Institute template, and finally smoothed using a 8.0 mm full-width at half-maximum Gaussian filter. Thus, both preprocessed structural MRI and fMRI images were in the same space and had the same voxel size (1 × 1 × 1 mm).
The preprocessed fMRI image data were analyzed with spatiotemporal PLS (McIntosh and Lobaugh, 2004; McIntosh et al., 2004), which determines time-varying distributed patterns of brain activity associated with either some aspects of the experimental design or a behavioral measure. This brain activity encompasses those voxels whose signal change covaries with the experimental conditions/behavioral measure in the same way. PLS assesses activation changes across all regions of the brain simultaneously rather than tessellations of regions, thereby eliminating the need for regional multiple-comparison correction. Additionally, such an analysis accommodates all task conditions into the analyses simultaneously, which facilitates the identification of both common and unique brain activity across conditions. A detailed description of the spatiotemporal analyses for (event-related) fMRI data has been given previously (Salami et al., 2010). In short, a data matrix was created by including each individual's fMRI block onsets and durations for each condition. We also used SPM's hemodynamic response function (HRF) model to create the data matrix. A cross-block covariance matrix between changes in brain activity and experimental conditions was then subjected to mean centering by subtracting the mean of each column from each value of that column. Singular value decomposition (SVD) was conducted on the mean-centered data matrix to reveal orthogonal latent variables (LVs), which reflect cohesive patterns of brain activity related to the experimental conditions. The first and successive LVs account for the greatest and progressively less amount of the cross-covariance matrix, respectively. Within each LV, voxel and design saliencies represent the relation of each voxel and experimental design to the LVs, correspondingly. In addition, brain scores and design scores were computed for each LV. The former was obtained as a dot product of each subject's image volume and voxel saliencies of each LV to derive an estimate of how strong each subject contributes to the pattern of each LV. The latter was computed as a dot product of each subject's image volume and design saliencies of each LV to derive an estimate of how much each particular task is related to weighted voxels of each LV.
In correspondence with the number of tasks, three LVs were extracted from the SVD. To assess the statistical significance of each LV, 500 permutation tests were conducted to determine whether the overall pattern in each LV was significantly different from randomness. The statistical cutoff was set to p < 0.001. In addition, 100 bootstrap tests were performed to derive an estimate of how reliably each voxel salience contributed to the pattern of LVs. The reliability of each voxel was assessed by the bootstrap ratio (BSR), which is the ratio of the voxel salience to the SE of that salience. Voxels with BSR > 4, which approximate a p value of 0.0001, were considered reliable. In addition, bootstrap resampling was used to generate 95% confidence intervals (CI) around brain scores, which provide a reliable measure to facilitate the interpretation. For example, significant difference between brain scores in different conditions is indicated by nonoverlapping CIs. Similarly, brain scores in a condition were deemed to be unreliable if the CI for that condition crossed zero (McIntosh and Lobaugh, 2004). All the reported clusters consisted of at least 20 contiguous voxels with a minimum distance of 10 mm between each cluster.
To investigate age-related functional differences during the episodic tasks, a voxelwise general linear model (GLM) was set up for each participant to generate a subject-specific contrast (encoding or retrieval vs baseline). To this end, each condition (encoding, retrieval, and baseline) was modeled as a box-car function that was convolved with the HRF. The six realignment parameters extracted from the realignment algorithm were included in the model as covariates of no interest to account for motion artifacts. Subject-specific contrasts were taken into a second level random-effects model using multiple regressions with age and the square of the orthogonalized age as covariates. Local maxima with p < 0.05 (familywise error corrected), and an extent threshold of 20 contiguous voxels (k > 20) were considered significant for both main and age effects. The age effect was investigated within a mask created from the main effect, thereby reducing the number of comparisons to correct for.
To investigate age-related structural differences, a multiple regression analysis was performed with smoothed normalized GM volumes as a dependent variable and age and the square of orthogonalized age as independent variables. To adjust for intergender variation in brain volume, a subject-specific total intracranial volume was computed by summing the GM and WM segments and was set as a covariate of no interest in the multiple regression analysis.
To explore whether age-related functional differences were driven by local GM atrophy, the Biological Parametric Mapping (BPM) toolbox was used (Casanova et al., 2007). BPM enables multimodal image analyses that model changes in one imaging modality as a function of another modality using the GLM framework. In this context (using GLM), BPM is connatural to SPM, but in addition can use image regressors. Thus, the design matrix in the BPM analysis is voxel specific in contrast to the identical design matrix for all voxels in the traditional SPM analysis. Local maxima with p < 0.05 (false discovery rate corrected) were considered to be significant.
Results
Behavioral results
During the retrieval task, 65.64% (15.75 ± 3.13 items) of the faces were correctly retrieved whereas 21.66% (5.19 ± 2.82 items) and 12.70% (3.04 ± 2.59 items) were associated to incorrect responses and nonresponses, respectively. The number of correctly retrieved names during scanning declined as a function of increasing age (r = −0.47; p < 0.0001). Also, mean reaction time for correctly retrieved items increased as a function of increasing age (r = 0.49; p < 0.0001).
fMRI results
Two significant LVs were identified in the task PLS analysis. LV1 (p < 0.0001; Fig. 1A) accounted for 93.73% of the cross-block covariance and identified brain regions differentiating both task conditions from baseline, with the largest differentiation between the retrieval and baseline conditions (Fig. 1B). Peak voxels from regions contributing to the pattern in LV1 are given in Table 1. The signal changes for selected regional maxima are shown in Figure 1, C and D. Regions showing greater activity during both tasks compared with baseline included the inferior parietal lobule, the middle and inferior frontal gyrus, the fusiform gyrus, the cerebellum, the posterior hippocampus, the occipital gyrus, the lingual gyrus, and the insula, bilaterally. Brain regions with relatively more activity during baseline included bilateral superior and middle temporal gyrus, bilateral putamen, bilateral superior and medial frontal gyrus, bilateral fusiform gyrus, and right thalamus.
LV2 (p < 0.0001; Fig. 1E) accounted for 6.37% of the cross-block covariance. This LV identified brain regions whose activity differentiated between the encoding and retrieval conditions, whereas the baseline condition did not reliably contribute to this pattern (CI crossing zero; Fig. 1F). Table 2 summarizes the significant regional saliencies in LV2. The positive saliencies correspond to greater activity during encoding than during retrieval and included bilateral parietal cortex, bilateral middle frontal gyrus (more extended in the left hemisphere), bilateral anterior hippocampus, bilateral superior occipital gyrus, left middle occipital gyrus, and right fusiform gyrus. The negative saliencies reflected greater activity during retrieval than during encoding and included bilateral cerebellum, right lingual gyrus, left middle occipital gyrus, right precuneus, the right insula, and the right inferior frontal gyrus. The signal change for specific regional maxima is shown in Figure 1, G and H, and confirms previous reports of asymmetric encoding-retrieval activations in the prefrontal cortex (PFC; Nyberg et al., 1996a; Habib et al., 2003).
The PLS analysis identified distinct encoding and retrieval networks (LV2), and separated these from a process-general network (LV1). Figure 2, A and C, illustrates the process-general recruitment of the posterior MTL and the encoding-specific recruitment of anterior hippocampus, respectively (Lepage et al., 1998; Schacter and Wagner, 1999; Persson et al., 2011). It should, though, be stressed that the pattern as a whole provides a better relation to the experimental conditions than any individual region. To further characterize the encoding and retrieval networks, individual brain scores were correlated with age, memory performance from the retrieval phase of the FN-PA, and performance on the different measures in the off-line cognitive battery (see Table 3).
LV1 identified a process-general effect, common to both encoding and retrieval, and the brain scores from encoding and retrieval were therefore averaged for the correlation analyses. The brain scores showed a significant positive correlation with age (r = 0.16; p < 0.01) as well as significant negative correlations with episodic memory performance as measured during scanning (r = −0.21; p < 0.01) and in the off-line battery (r = −0.16; p < 0.01). Thus, the brain network identified by LV1 was recruited to a greater extent by older adults, and stronger engagement of this network correlated with lower episodic memory performance.
Next, brain scores for encoding and retrieval, identified by LV2, were related to age and cognitive performance. Encoding brain scores were negatively correlated with age (r = −0.13; p < 0.01), showing that older adults tended to recruit this network to a lesser degree. Also, the encoding brain scores showed a selective positive correlation with episodic memory performance (scanner performance: r = 0.16, p < 0.01; off-line battery: r = 0.18, p < 0.01), but not with the performance on nonepisodic tasks (Table 3). Thus, stronger engagement of this network during encoding was related to higher levels of episodic memory performance. For retrieval, the correlations were nonsignificant.
The opposing relation between age and recruitment of the process-general (LV1) and process-specific (LV2) networks indicated that an age effect on the hippocampus-based encoding network cascades into upregulated recruitment of the general network. By this view, elderly with more pronounced decline in the episodic encoding network should engage the general network to a particularly high degree. We tested this prediction by identifying a subgroup of elderly (>60 years old, n = 90; mean age: 64.40 ± 7.36 years; mean hit rate: 15.21 ± 3.05 items) with lowest engagement of the process-specific encoding network according to a median-split analysis. We then examined whether this group showed higher recruitment of the general network as compared with age- and performance-matched elderly counterparts (n = 91; mean age: 64.12 ± 6.69 years; mean hit rate:15.04 ± 2.85 items) who engaged the process-specific encoding network to a greater degree. A significant difference was observed, confirming that lower process-specific network recruitment predicted higher engagement of the process-general network (t(179) = 2.52, p < 0.01; Fig. 3).
The bilateral hippocampus was a salient component of the encoding network. We supplemented the analyses of correlations between overall brain scores and age/memory performance, with correlations of β values (Table 3) for bilateral hippocampus. In keeping with the correlation patterns for the overall brain scores, the values for bilateral hippocampus were significantly correlated with age, scanner memory performance, and off-line episodic memory performance. In addition, significant correlations were found for letter digit task performance.
Given the significant correlation between encoding brain scores and age, as well as the correlations of bilateral hippocampal scores with age, a multiple regression analysis across the whole brain was performed to assess age-related functional decline during episodic encoding. The results confirmed a significant association of age and bilateral hippocampus recruitment (left: xyz = −15 −13 −19, t(290) = 4.59, k = 989; right: xyz = 18 −10 −18, t(290) = 4.56, k = 910) as well as the left lingual gyrus (xyz = −13 −93 −13; t(290) = 4.16; k = 1043). Figure 4 illustrates the pattern of age-related functional decline, showing a nonlinear trend, with accelerated hippocampal functional decline in the oldest age groups [left: r2 (linear) = 0.67; r2 (nonlinear) = 0.73; Δr2 = 6%; right: r2 (linear) = 0.59; r2 (nonlinear) = 0.69; Δr2 = 10%]. Control analyses showed that the age-related diminished hippocampal activation persisted after controlling for age-related changes in memory performance during scanning (left: t(289) = 3.89, p < 0.0001; right: t(289) = 3.33, p < 0.0005).
The plot of signal changes in the bilateral hippocampus identified some elderly with very low and nonsignificant hippocampal activation (Fig. 4, insert). These “low activators” were screened for ApoE status and Mini Mental State Examination (MMSE) performance (Folstein et al., 1975). Only two were ApoE e4-carriers and all scored >26 on the MMSE. Moreover, their performance on the tasks in the Betula test battery was within normal levels. Collectively, these results suggest that the low activators were unlikely to be in a preclinical phase of dementia (and all passed clinical screening for dementia before scanning).
Structure–function integration
The multiple regression analysis of GM volume revealed marked age-related volume reductions in several regions, including the thalamus, amygdala, putamen, hippocampus, the inferior frontal gyrus, and the middle temporal gyrus, bilaterally (Table 4).
In a final set of analyses, functional and structural data were integrated using BPM. Of chief interest was to investigate whether the observation of diminished hippocampal functional responses in aging could be accounted for by age-related local changes in brain GM. The multiple regression approach showed that age-related diminished hippocampal activation persisted after controlling for GM volume, but with smaller extent and a reduced peak value (left: xyz = −21 −16 −18, t(290) = 3.91, k = 815; right: xyz = 20 −10 −18, t(290) = 3.89, k = 585) (Fig. 5).
Discussion
This study provides novel evidence for opposing effects of age on memory encoding and cognitive control at the systems level. In contrast, most previous studies investigated regional effects of age on memory encoding, retrieval, and cognitive control (i.e., limited evidences at the network level). Using multivariate PLS, we dissociated a process-general (encoding/retrieval) cognitive network from process-specific encoding and retrieval networks. The general network involved large portions of frontoparietal cortex and also cerebellum and posterior hippocampus. Middle and inferior frontal gyrus have been associated with several general cognitive-control processes, including top-down attentional control (Banich et al., 2000; Hopfinger et al., 2000; Naghavi and Nyberg, 2005) and attention to targets (Corbetta et al., 2002), as well as response inhibition and selection (Liddle et al., 2001; Corbetta and Shulman, 2002; Aron et al., 2004). Moreover, the inferior and middle frontal gyrus as well as the inferior parietal lobule have been implicated in spatial attention (Naghavi and Nyberg, 2005). The posterior hippocampal activation was common to both conditions, suggesting that this part of the hippocampus is involved in processes that are taxed during both encoding and retrieval. The present face-name task should engage associative operations at both processing stages (Nyberg et al., 2000), and previous imaging findings relate associative memory operations to posterior segments of the hippocampus (Henke et al., 1997, 1999; Kirwan and Stark, 2004).
The second LV revealed distinct networks for encoding and retrieval. Differential patterns of prefrontal recruitment contributed to this effect. Consistent with the hemispheric encoding/retrieval asymmetry (HERA) model (Tulving et al., 1994; Nyberg et al., 1996a; Habib et al., 2003), a left PFC region was involved during episodic encoding (Fig. 1H) whereas a right PFC region was more involved during episodic retrieval (Fig. 1G). Notably, in keeping with previous findings [(Schacter and Wagner, 1999; Persson et al., 2011], bilateral anterior hippocampus was differentially more engaged during episodic encoding. This effect might reflect binding processes that are of particular relevance for the formation of new associations in memory (Chua et al., 2007).
Brain scores from the general network and the process-specific networks were related to age and cognitive performance and strikingly different patterns were observed. The process-general brain network was recruited to a greater extent by older adults (positive correlation between brain scores and chronological age). Stronger engagement of this network predicted lower performance on episodic memory tasks. In contrast, the encoding network was recruited to a lesser degree by older adults, and greater engagement of this network during encoding was related to higher performance levels for episodic memory task. Moreover, elderly with lesser engagement of the process-specific encoding network recruited the process-general network more than elderly with greater engagement of the process-specific encoding network. Together, these findings suggest that age-related episodic memory impairment may be caused by inefficient use of the encoding network (Naveh-Benjamin et al., 2003), which in turn cascades into a greater reliance on general cognitive-control processes (Gazzaley et al., 2008).
The analyses converged to show that age-related functional changes in the hippocampus contributed to inefficient encoding. Both simple correlational and multiple regression analyses (across the entire brain) confirmed that the strongest age-related functional decline was seen in the anterior hippocampus. This effect remained after controlling for age-related performance differences. It has been debated whether hippocampal structure and function is affected in normal, nonpathological aging (for review, see Buckner, 2004; Hedden and Gabrieli, 2004). Some previous studies reported age-related hippocampus hypoactivation (Gutchess et al., 2005; Daselaar et al., 2006), whereas others found preservation of hippocampus responses across the examined age span (Persson et al., 2011). The present results indicate that age-related reductions in hippocampal functional responses are relatively minor before age 70, suggesting that differences in age ranges across studies may account for some between-studies discrepancies. Elderly who displayed weak hippocampus engagement were scrutinized for ApoE status, MMSE score, and cognitive status. No evidence was obtained to support the view that hippocampal under-recruitment was related to pathological processes such as dementia.
The structural analyses provided evidence for age-related atrophy in several brain regions, including the hippocampus. The observed pattern was broadly consistent with that seen in previous studies (Raz et al., 2005; Fjell et al., 2009). We have previously observed that structural changes in frontal cortex were related to age-related reduction of prefrontal activation during cognitive tasks (Nyberg et al., 2010). Here, when functional and structural data were inter-related, overlap was observed in left anterior hippocampus, indicating that structural changes in the hippocampus, at least in part, contributed to the age-related under-recruitment of the functional encoding network. However, the age-related structural change did not completely mediate the age-related changes in hippocampus activation (i.e., the effect persisted after controlling for GM volume; Fig. 5), and several other factors might be of relevance, including age-related molecular changes (Chételat et al., 2011) and changes in the pattern of the hypothalamic-pituitary adrenal axis (Miller and O'Callaghan, 2003).
The PLS analysis of the data from the present large-scale dataset was instrumental for the identification of a general network along with specific episodic memory networks. Importantly, PLS opened up for the use of subject-specific brain scores as indicators of a person's large-scale brain network status. Such integrated scores could be valuable in imaging genetics studies, in particular in efforts to combine imaging and genetic data from genome-wide association studies. A potential limitation is that the findings were derived from a blocked statistical model, as PLS is not a suitable tool for fast event-related designs (McIntosh et al., 1996; McIntosh and Lobaugh, 2004). Also, in the current study, only a subsample of participants fulfilled the criterion of having a sufficient number of items (remembered vs forgotten). However, the event-related control analysis using a multiple regression approach in SPM revealed similar encoding and retrieval networks, as well as a comparable age-related functional decline in the anterior hippocampus. Another, perhaps more serious, limitation is that the present inferences on age-related changes were based on cross-sectional rather than longitudinal data. Previous studies show that cross-sectional estimates of true age changes may be misleading (Rönnlund et al., 2005; Nyberg et al., 2010). However, cross-sectional and longitudinal results tend to converge for older cohorts (Rönnlund et al., 2005), indicating that the present findings of hippocampal functional reduction that primarily was expressed in the oldest cohorts may be valid.
In conclusion, the results of this study support a model where age-related structural brain changes, in the hippocampus and elsewhere, translates into age-related reductions in the efficient recruitment of a process-specific encoding network along with upregulated recruitment of a general cognitive control network. As such, the findings suggest that age-related decline in memory can be linked to encoding processes as well as to general cognitive-control operations.
Footnotes
This study was supported by the Göran Gustafsson Award in Medicine (L.N.), a grant from the Swedish Science Council (L.N.), and a Wallenberg Scholar Grant from the Knut and Alice Wallenberg Foundation (L.N.). The BETULA Project is supported by a grant from the Swedish Science Council (L.-G. Nilsson and L.N.). We thank the staff of the BETULA Project, R. Adolfsson and L.-G. Nilsson, and the staff at the Umeå Center for Functional Brain Imaging (especially Micael Andersson).
- Correspondence should be addressed to Alireza Salami, Department of Integrative Medical Biology, Umeå Center for Functional Brain Imaging, Umeå University, S-901 87 Umeå, Sweden. alireza.salami{at}physiol.umu.se