Age-related changes in brain structural covariance networks

Introduction

The majority of neuroimaging studies of aging have reported a consistent pattern of gray matter (GM) volumetric reductions in the human cortex, involving mainly prefrontal regions, parietal, and temporal association cortices, and the insula and cingulum (Resnick et al., 2003; Sowell et al., 2003; Raz et al., 2005; Du et al., 2006; Terribilli et al., 2011). However, more and more studies have demonstrated that cerebral changes with normal aging are not simply characterized by regional alterations but rather by the reorganization of cortical connectivity patterns (O'Sullivan et al., 2001; Koch et al., 2010; Wu et al., 2011b; Zhu et al., 2012). Using diffusion tensor imaging (DTI), several studies have consistently reported a diffuse loss of axonal integrity in senior populations (Salat et al., 2005; Pagani et al., 2008; Madden et al., 2009), which allowed for inferences regarding changes of structural connectivity in older people compared with younger adults. Building on complex network analysis methods, Gong et al. (2009) reported a reduction in overall cortical connectivity, decreased local efficiency, and a shift in regional efficiency from parietal and occipital to frontal and temporal neocortex in older brains.

The investigation of structural covariance networks (SCNs) using structural magnetic resonance imaging (sMRI) is another useful method to explore structural brain networks. This approach mainly characterizes the pattern of structural covariance in GM morphology (e.g., volume, thickness and surface area) between brain regions using a general linear model (GLM) framework (Mechelli et al., 2005; Lerch et al., 2006; Nosarti et al., 2010; Zielinski et al., 2010; Montembeault et al., 2012; Soriano-Mas et al., 2012). Many studies have demonstrated underlying relationships among brain areas using structural correlation by sMRI, anatomical connectivity by DTI, and functional correlation by resting-state functional MRI. For example, He and colleagues found that structural networks based on cortical thickness measurements were compatible with known functional networks (He et al., 2007). Greicius et al. measured functional connectivity using independent component analysis and anatomical connectivity using DTI and found there existed white matter tract structural connections between functionally connected regions (Greicius et al., 2009). A recent study reported agreement in the correlations in GM thickness and underlying fiber connections across brain areas, but more information was included for the thickness network than the fiber network (Gong et al., 2012). In addition, Seeley and colleagues demonstrated that SCNs using voxel-based morphometry (VBM) were able to recapitulate the functional connectivity network topologies (Seeley et al., 2009). The covariance of different cortical regions in their GM volumes was considered to be the result of mutual trophic influences (Ferrer et al., 1995) or common experience-related plasticity (Draganski et al., 2004; Mechelli et al., 2004). The consistency among these three networks provides substantial support for SCNs serving a measure of network integrity in cross-sectional studies.

Of note, the analysis of SCNs using VBM has been successfully applied to map the eight SCNs and explore how neural systems build large-scale structural covariance during development (Zielinski et al., 2010). One previous aging study using structural covariance approaches compared two populations (younger vs. older subjects) and reported reduced structural associations in order adults, specifically in high-order cognitive networks (Montembeault et al., 2012). However, because subjects in only two age categories were involved in that experiment, it could not determine the aging trajectories of these sensorimotor and high-order cognitive networks. Zielinski et al. (2010) found that there were non-linear trajectories of primary visual, auditory, and sensorimotor networks during development. Moreover, these networks were provisionally established by early childhood, but underwent significant expansion in early adolescence before contraction or pruning in late adolescence. Many studies have shown that age-related atrophy of some neural regions follows variable, non-linear patterns (Allen et al., 2005; Kennedy et al., 2009; Terribilli et al., 2011). However, it remains unknown whether and how age-related changes of primary sensorimotor and high-order cognitive SCNs exhibit non-linear trajectories during normal aging.

Here, we used an SCN mapping approach to investigate eight large-scale networks in 240 healthy participants aged 18–89. These participants were subdivided into young (18–23 years), middle aged (30–58 years), and older (61–89 years) subjects. Given that SCNs subserving language, social–emotional, and executive control functions have shown gradual deterioration through normal aging, we expected to observe age-related trajectories in the deterioration of these large-scale networks. To test our hypothesis, we first investigated the aging trajectories of three sensorimotor and five high-order cognitive SCNs in three groups composed of an equal number of subjects at different ages, and then compared these SCNs differences between the groups.

Materials and Methods

Participants

Three hundred sixteen right-handed, healthy subjects were selected from the Open Access Series of Imaging Studies (OASIS) cross-sectional database (http://www.oasis-brains.org) (Marcus et al., 2007). Data from seven subjects were excluded from further analysis because of poor image quality or image preprocessing. Because of relatively few subjects aged 30–60 years, we first selected 80 subjects from this age range as the middle-aged group (30 males and 50 females). Then, we selected 80 separate subjects around 75 years of age as the old group. Finally, we selected the youngest 80 subjects matched for gender as the young group. The names and characteristics of the groups are shown in Table 1. All subjects were evaluated using the Mini-Mental State Examination (MMSE) (Folstein et al., 1975) and Clinical Dementia Rating (CDR) scales (Morris, 1993; Morris et al., 2001). MMSE scores were higher than 29, and CDR scores were all zero. For demographic data on all subjects, see Marcus et al. (2007). This dataset has been used in several previous studies (Bakkour et al., 2009; Fjell et al., 2009; Salat et al., 2009; Li et al., 2011).

TABLE 1

Table 1. Characteristics of the subjects in this study.

Image Acquisition

For each subject, three to four individual T1-weighted magnetization-prepared rapid gradient echo (MP-RAGE) images were acquired on a 1.5T Vision scanner (Siemens, Erlangen, Germany) within a single session. Head movement was minimized with cushioning and a thermoplastic facemask. Images were motion corrected and averaged to create a single image with a high contrast-to-noise ratio. MP-RAGE parameters were empirically optimized for gray/white contrast: TR = 9.7 ms; TE = 4 ms; flip angle = 10°; slice number = 128; resolution = 256 × 256 (1 × 1 mm); thickness = 1.25 mm.

Image Preprocessing

Structural MR images were processed using a technical computing software program (MATLAB 2010; The MathWorks Inc., Natick, Mass) and Statistical Parametric Mapping software (SPM 8; The Wellcome Department of Imaging Neuroscience, London, UK). Following the inspection of image artifacts, image preprocessing was performed with the VBM8 toolbox (http://dbm.neuro.uni-jena.de/vbm/). Briefly, all native-space MRIs were segmented to extract GM based on an adaptive maximum a posteriori technique (Rajapakse et al., 1997) and partial volume estimation method (Tohka et al., 2004) without the need for a priori tissue probability information. In addition, a spatially adaptive non-local denoising filter (Manjon et al., 2010) and a hidden Markov random field model (Rajapakse et al., 1997) were applied to minimize the level of noise in the resulting GM segments. Subsequently, the high-dimensional DARTEL (diffeomorphic anatomical registration using exponentiated lie algebra) approach provided non-linear deformation to normalize the images to the DARTEL template in Montreal Neurological Institute (MNI) space, which was derived from 550 healthy control subjects (Ashburner, 2007). Non-brain tissue was also removed. Additionally, the Jacobian determinants derived from the spatial normalization were used to modulate the GM value for each voxel to preserve the total amount of GM from the original images (Good et al., 2001). We used non-linear components only, which allowed us to analyze relative (i.e., corrected for individual brain size) differences in regional GM volume. Finally, the resulting modulated and normalized images were smoothed with a 12 mm full width at half maximum isotropic Gaussian kernel.

Extraction of Seed Volumes

To assess the structural covariance pattern of each large-scale network, we extracted individual GM volumes from eight seed regions of interest (ROIs). These ROIs included the primary visual, auditory, and motor cortex, as well as language-related speech and semantic areas, areas related to salience and executive control, and the default-mode network (DMN). We based our ROIs on previous studies (Zielinski et al., 2010; Montembeault et al., 2012). For each region, seeds were defined with the MarsBar ROI toolbox (http://marsbar.sourceforge.net/) as 4-mm radial spheres centered at the following MNI coordinates: right calcarine sulcus (9, −81, 7), right Heschl's gyrus (46, −18, 10), right precentral gyrus (28, −16, 66), left inferior frontal gyrus, pars opercularis (IFGo) (−50, 18, 7), left temporal pole (−38 10, −28), right frontoinsular cortex (38, 26, −10), right dorsolateral prefrontal cortex (DLPFC) (44, 36, 20), and right angular gyrus (46, −59, 23). In addition, we selected contralateral seeds by changing the sign on each seed's x coordinate to perform similar SCNs analyses (Table A1, and Figure A1).

Statistical Analysis

A voxel-wise statistical analysis was performed on the GM images using the GLM as applied in SPM8. Eight multiple regression models were used to test the strength of the structural covariance between each seeds and all other regions across whole-brain GM for each age group. In each regression model, the extracted mean GM volume from each ROI was entered as a covariate of interest, and gender as a confounding covariate. Because of the unequal sample size across the genders, we removed the effects of gender on the SCN patterns in the correlation analysis. These statistical analyses identified voxels that had a positive covariance with each a priori selected ROI in each group. The resulting correlation maps were using height and extent thresholds at P < 0.05 with family-wise error (FWE) correction. They were displayed on the MNI template to allow qualitative comparisons between the age groups using BrainNet Viewer software (http://www.nitrc.org/projects/bnv/). To quantify differences in normal aging trajectories across networks, we calculated the total number of significant positive ipsilateral, contralateral, and whole-brain voxels and plotted these across the age groups.

To further examine the effects of age on specific regional covariance, we performed between-group difference analyses of SCN patterns for any two groups according to the scientific literature (Lerch et al., 2006). For any pair of voxels in two groups, their structural correlation may have different slopes, and the difference in slope may represent the difference in their structural association. The difference in slope was tested using a classic interaction linear model:

$V_i={\beta }_0+{\beta }_1{V}_j+{\beta }_2\text{Group}+{\beta }_3(V_j\times \text{Group})+\epsilon$

where V_i and V_j represented the GM volumes of a pair of voxels in two groups. The Group component was modeled using treatment contrasts, and significance tested using the Student's t statistic. Specific t-value contrasts were established to map the voxels that expressed a significantly different structural association between any two groups. The threshold for the resulting statistical parametric maps was given at P < 0.05 for the height and extent thresholds, with FWE multiple comparisons correction.

Results

For each network, the age-related SCN trajectory was identified by the spatial extent of each SCN. These results were demonstrated by functional domain, as described in the following sections. We also explored regions that significantly differed in their structural associations between any two groups.

Primary Sensory and Motor Networks

Seeds within primary visual cortex (right calcarine sulcus) produced SCNs with a relatively preserved pattern (Figure 1). However, there were small changes in the distributions of the covariance maps. In the young and middle-aged groups, the covariance regions mainly included the bilateral calcarine sulcus, lingual gyrus, cuneus and right lateral occipital gyrus, whereas only included bilateral calcarine sulcus, lingual gyrus, and cuneus in the old group. Primary auditory cortex (i.e., the right Heschl's gyrus) covaried with the bilateral insula and precuneus, the right posterior cingulate cortex (PCC), parahippocampus, and inferior frontal gyrus, and the left orbital-frontal cortex in the young participants, and underwent significant contraction in the middle-aged and old groups to include only the bilateral insula regions. There was a flat transition in the covariance maps during aging. Primary motor cortex (right precentral gyrus) correlated with the ipsilateral precentral regions in the young participants, which progressed at middle age to include the contralateral precentral and supplementary motor areas. In old participants, these structural associations were similar to the young participants.

FIGURE 1

Figure 1. Age-related changes in primary sensory and motor structural covariance networks. (A) Statistical maps of regions significantly correlated with the seed region in each group. The results are presented as correlation coefficient values (P < 0.05, FWE corrected). (B) The plots of voxel counts by group indicate small covariance changes throughout the life-span. The y-axis represents the voxel number (×10⁴). Abbreviations: Calcs, calcarine sulcus; HG, Heschl's gyrus; PreCG, precentral gyrus; R, right; L, left; Y, young group; M, middle-aged group; O, old group; CC, correlation coefficient.

In addition, within the primary sensory and motor networks, there were no significant differences observed in the structural association between any two groups when we compared the slopes of the structural associations on a voxel-by-voxel basis.

Language-Related Speech and Semantic Networks

The two language-related speech and semantic networks followed similar variation patterns (Figure 2). From the young participants to the middle group, the covariance regions showed abrupt decreases. They then showed smaller changes at the following ages. For the language-related speech network, the covariance map involved the bilateral anterior insula, medial frontal, cingulate cortex, and temporal regions in the young participants, and then contracted to include seed autocorrelation in the older participants. The semantic regions correlated with the bilateral temporal cortices, cingulate gyrus, insula and frontal regions in the young participants, and shrank sharply to include only the anterior temporal cortices during normal aging.

FIGURE 2

Figure 2. Age-related changes in speech and semantic structural covariance networks. (A) Statistical maps of regions significantly correlated with the seed region in each group. The results are presented as correlation coefficient values (P < 0.05, FWE corrected). (B) The plots of voxel counts by group indicate abrupt contraction in the middle-aged group and mild changes in the old group. The y-axis represents the voxel number (x× 10⁴). Abbreviations: IFGo, inferior frontal gyrus, pars opercularis; TPole, temporal pole; L, left; R, right; Y, young group; M, middle-aged group; O, old group; CC, correlation coefficient.

Comparing the other SCNs, language-related speech and semantic networks showed more age-related changes. Specifically, the left supplementary motor and superior temporal areas showed significant positive associations with the speech seed region in young adults, whereas this covariance disappeared in the middle-aged group (Table 2 and Figure 3). Decreased positive associations between the left superior temporal region and left temporal pole were found in the middle-aged group compared with the young group (Table 2 and Figure 4). There were significant positive associations between the left anterior and the PCC and left IFG in the young participants, and this covariance became a negative correlation in the old group (Table 2 and Figure 4). When compared with the old group, there was a significant positive association between the right PCC and the semantic seed region in the young participants (Table 2 and Figure 5).

TABLE 2

Table 2. Contrast analysis of structural covariance network trajectories.

FIGURE 3

Figure 3. Age-related changes in salience, executive control, and default-mode structural covariance networks. (A) Statistical maps of regions significantly correlated with the seed region in each group. The results are presented as correlation coefficient values (P < 0.05, FWE corrected). (B) The plots of voxel counts by group indicate abrupt contraction in the middle-aged group and mild changes in the old group. The y-axis represents the voxel number (× 10⁴). Abbreviations: DMN, default-mode network; DLPFC, dorsolateral prefrontal cortex; FI, frontoinsular cortex; ANG, angular gyrus; R, right; L, left; Y, young group; M, middle-aged group; O, old group; CC, correlation coefficient.

FIGURE 4

Figure 4. Group differences for the young vs. the middle-aged subjects. Significant between-group differences within the DMN (A), speech (B), and semantic (C) structural covariance networks were found. For each network, the region of interest (upper) and region showing the most significant structural association (lower) are presented on the right, and a plot of slop differences between the seed region and a 4-mm radius sphere centered on the peak voxel are presented in the left. Voxels with P_FWE < 0.05 are displayed. Abbreviations: DMN, default-mode network; ANG, angular gyrus; IFGo, inferior frontal gyrus, pars opercularis; TPole, temporal pole; PFC, prefrontal cortex; SMA, supplementary motor area; STC, superior temporal cortex; R, right; L, left; Y, young group; M, middle-aged group.

FIGURE 5

Figure 5. Group differences for the young vs. the old subjects. Significant between-group differences within the DMN (A), speech (B), semantic (C), and executive (D) structural covariance networks were found. For each network, the region of interest (upper) and region showing the most significant structural association (lower) are presented on the right, and a plot of slop differences between the seed region and a 4-mm radius sphere centered on the peak voxel are presented in the left. Voxels with P_FWE < 0.05 are displayed. Abbreviations: DMN, default-mode network; ANG, angular gyrus; IFGo, inferior frontal gyrus, pars opercularis; TPole, temporal pole; DLPFC, dorsolateral prefrontal cortex; PCC, posterior cingulate cortex; ACC, anterior cingulate cortex; R, right; L, left; Y, young group; O, old group.

Salience, Executive Control, and Default-Mode Networks

The three networks associated with social–emotional and cognitive function showed similar patterns to the language-related speech and semantic networks, with more distributed structural covariance in the young participants and limited covariance patterns throughout the later stages (Figure 5). Specifically, the right frontoinsular cortex anchored covariance maps included extensive areas of the bilateral lateral and medial frontal cortex, temporal cortices, and cingulate regions. This SCN underwent significant shrinkage in the middle-aged group to include only bilateral insular regions. In the older subjects, this network somewhat extended to include the bilateral medial prefrontal regions. The right DLPFC seed covaried with the bilateral frontal, temporal, and cingulate cortices in the young participants, but shrank into a more focal distribution in the middle-aged and older subjects. A seed in right angular gyrus produced an aging SCN, including the bilateral angular gyrus, middle temporal, cingulum, and prefrontal regions in the young participants, and contracting to include the bilateral PCC and angular gyrus in the middle-aged group, and only the bilateral angular gyrus in the old group.

There were no significant differences observed in the structural associations between any two groups within the salience network. The left PCC showed a significant positive association with the executive control seed region and right angular region (DMN seed) in the young adults, whereas this covariance disappeared in the old group (Table 2 and Figure 5). Furthermore, a less robust positive association between the right prefrontal region and DMN seed region was found in the middle-aged group compared with the young participants (Table 2 and Figure 4).

Discussion

In this study, we investigated the age-related trajectories of eight large-scale networks using an SCN mapping approach in three distinct age groups. All networks exhibited non-linear patterns across normal aging in terms of the spatial extent of the network. Except for the primary motor network, these networks showed a more distributed topology in young participants, which shrank sharply to a more localized topology in the middle-aged group, and maintained this localized topology in the old group. Primary sensory and motor networks showed fewer age-related changes compared with high-order cognitive networks. Moreover, there were significant between-group differences in language-related speech and semantic networks, the executive control network, and the DMN. Taken together, our results provide evidence of variations in the topological organization of SCNs during normal aging.

Comparison of Structural Covariance Networks and Other Networks

The cerebral cortex is organized into networks of functionally complementary areas. New advances in modern neuroimaging techniques and quantitative analysis of complex networks have made the investigation of brain network topological organization possible. DTI is a useful tool for non-invasively mapping cortico-cortical anatomical connections by examining axonal integrity. Insight into the structural covariance of GM morphology (e.g., volume, thickness, and surface area) between brain regions is provided by sMRI. Resting-state fMRI has allowed for assessments of the strength of functional connections within a network by quantifying correlated activity [i.e., spontaneous, low-frequency fluctuations in the blood oxygen level-dependent (BOLD) signal] between brain regions at rest.

The current studies have demonstrated the underlying relationships among structural correlations, anatomical connectivity, and functional correlations. For example, structural networks based on cortical thickness measurements were consistent with known functional networks (He et al., 2007). A recent study reported agreement in the correlations in GM thickness and underlying fiber connections across brain areas, but more information was included for the thickness correlation network than the fiber connection network (Gong et al., 2012). In addition, SCNs determine using VBM recapitulated the canonical intrinsic connectivity networks topologies (Seeley et al., 2009). The consistency among the SCNs and other networks provides substantial support for the use of SCNs as a measure of network integrity for cross-sectional studies.

Age-Related Changing Trajectory of SCNs

In this study, we employed three age groups to map the trajectory of SCN changes over age. Except the primary motor network, these SCNs appeared the similar non-linear pattern that had a distributed topology in young participants, a sharply localized topology in middle aged participants, and were relatively stable in older participants. This trajectory was similar to age-related changes of integrated local efficiency of brain network (Wu et al., 2011b). They employed the graph theory analysis method and reported that the local efficiency in the young group was significantly larger than those of the middle and old groups, whereas no significant difference was found between the middle and old groups. The shrinkage of SCNs as well as the reduction of local efficiency might be explained by the regionally distributed pattern of GM atrophy (Bergfield et al., 2010). As for primary motor network, we found inverted V-curve tendency among three age groups in the SCN of primary motor cortex. The previous study reported increased inter-regional correlations between bilateral primary motor cortex in the aging brain (Chen et al., 2011). In our study, we found the structural covariance of primary motor cortex between two hemispheres in the middle group. Alternatively, we found that the volume of right precentral gyrus reduced with age increasing (Figure A2), which might reflect a compensation mechanism that a network may need to work harder by becoming overactive due to regional volumetric atrophy with the network (Reuter-Lorenz and Cappell, 2008).

Intragroup Patterns and between-Group Differences of SCNs

Primary sensory and motor networks

In this study, we found that seeds within the primary visual, auditory, and motor cortices produced SCNs with smaller spatial extents (i.e., total voxel count exhibiting a significant correlation at the corrected threshold) compared with other networks. This is consistent with previous work (Lerch et al., 2006) with mapped anatomical correlations across cerebral cortex (MACACC) using cortical thickness. Similarly, they used the number of cortical region showing significant correlations with seed region to describe the strength of MACACC. They found the primary motor, sensorimotor, and visual areas had the lowest strengths of correlation, whereas the association cortices had the highest strength. This could be explained by the functions of association cortices, because these cortices receive and integrate inputs from multiple cortical and non-cortical sources, and distribute information to multiple areas.

In addition, primary sensory and motor networks showed fewer age-related changes compared with high-order cognitive networks. There were no significant differences observed in the structural association between any two groups when the slopes of the structural associations were contrasted on a voxel-by-voxel basis within the primary sensory and motor networks. This finding is in line with a recent study (Montembeault et al., 2012) in which between-group differences were only observed in high-order cognitive SCNs between young and old subjects, rather than primary sensory and motor networks.

Language-related speech and semantic networks

The language-related speech and semantic networks showed significant age-related changes, with a more localized topology with increasing age. These results could account for the decline of language abilities across in normal aging. Federmeier and colleagues observed age-related changes in the timing with which message-level information impacted semantic analysis using event-related potentials, which indicated that aging affects high-order language processes (Federmeier et al., 2003). Wierenga and colleagues also reported that healthy older adults had increased difficulty in word retrieval with unchanged semantic knowledge (Wierenga et al., 2008).

Within the language-related speech SCN, a reduced structural association was observed the between inferior frontal gyrus and the supplementary motor area in the middle-aged group compared with the young participants. The speech SCN linked the language and motor systems that enable speech fluency (Seeley et al., 2009). This reduced structural covariance may explain the decline of motor speech skills with the transference of language to a speech code in older adults. Similarly, there was a significant decrease in the covariance between the seed region and the anterior cingulate cortex (ACC) in the old group compared with the young participants. This could be related to the functional degradation of these areas. Takahashi et al. reported that significant age-related reductions in regional cerebral blood flow in the left IFG and the bilateral medial frontal gyri and ACC using single-photon emission tomography (Takahashi et al., 2005). In addition, as the ACC plays a critical role in the control of speech responses (Paus et al., 1993), this reduced covariance may explain some speech difficulties in senior populations.

In semantic SCN analysis, a reduced structural association between the left temporal pole and left superior temporal cortex was found in the middle-aged group compared with the young participants. Many studies have shown that the superior temporal cortex is associated with understanding spoken words (Demonet et al., 1992; Chao et al., 2002; Okada and Hickok, 2006). Thus, this reduced structural association may be related to a decline in auditory word comprehension with increasing age. In the comparisons of the young and old participants, the main differences we found were in the structural associations between the left temporal pole and PCC. The results of functional neuroimaging studies in young, healthy adults provide compelling evidence for the involvement of the PCC in memory retrieval. For example, PCC activation was elicited during recognition of thematic narrative information learned during training sessions (Maguire et al., 1999). The reduced structural association between these two regions may explain the decline of semantic memory in senior populations.

Salience, executive control, and default-mode networks

The three high-order cognitive SCNs showed a similar pattern to the language SCNs, showing a distributed topology in young participants that shrank sharply to a localized topology in the middle-aged group, and maintained a localized topology in the old group. Cepeda and colleagues found age-related decline of executive control processes by examining task-switching performance (Cepeda et al., 2001). Kelly and colleagues also demonstrated the stability of executive control decreased in older subjects (West et al., 2002). Several studies have shown that the DMN is altered in old subjects (Andrews-Hanna et al., 2007; Koch et al., 2010; Sambataro et al., 2010), especially in anterior regions (Damoiseaux et al., 2008). Similarly, Chen et al. found the decreased inter-regional correlations with the DMN module (Chen et al., 2011). Cognitive decline with age, such as cognitive control (Persson et al., 2007) and working memory (Sambataro et al., 2010), is associated with decreased DMN connectivity.

In the executive-control SCN analysis, the main differences we found were in the structural association between the right DLPFC and PCC, which exhibited a positive correlation in the young participants and no association in the old group. The DLPFC shows increased activity when experimental stimuli are presented, and is thought to support on-task processing in attention tasks. The PCC, in contrast, shows decreased activity during stimulus presentation, and is thought to support off-task processing in attention tasks. Some studies have reported that activity in the DLPFC shows a negative correlation with activity in the PCC during the resting state (Greicius et al., 2003; Fransson, 2005), and the negative coupling between these regions is lower in older adults (Sambataro et al., 2010). Our findings are inconsistent with functional connectivity between these areas, which should be further investigated in future studies.

Finally, for the DMN, we found a reduction in the structural association between the right angular gyrus and right prefrontal cortex from the young to the middle-aged participants. In the comparison of the young and old group, we found age-related differences in the structural association between the right angular gyrus and PCC. Similarly, Wu and colleagues reported that the left angular gyrus has a significantly reduced correlation with the PCC in older compared with younger participants in a resting-state fMRI study (Wu et al., 2011a). As the PCC is a prominent region within the DMN, the reduced structural association between the seed regions and the PCC indicates that the DMN SCN shrinks with age.

Laterality of the SCNs

We found similar SCNs across both hemispheres (Figure A1), consistent with previous SCN analyses of the developing brain (Zielinski et al., 2010). However, the auditory and speech networks showed a notable exception. The SCN derived from the left seed of Heschl's gyrus was distinctly smaller than the one derived from the right seed. In contrast, the SCN derived from the left seed of the IFGo was distinctly larger than the one derived from the right seed, which may explain why speech networks are left-dominated (Powell et al., 2006; Xiang et al., 2010).

Further Considerations

To build upon this study, several issues need to be addressed. First, we assigned subjects in their 30s to the middle-aged group because there were fewer subjects in their 40s and 50s in our sample. This resulted in wide age range in the middle-aged group. In future, this age range should be refined by including more subjects. Second, a previous study has reported the brain functional connectivity pattern could be affected by the sex factor (Biswal et al., 2010). In future, it is interesting to explore the sexual differences when the SCNs change with age. Third, it should be highlighted that this study used ROIs based on previous studies (Zielinski et al., 2010; Montembeault et al., 2012). Future studies using ROIs extracted from ICNs of the same subjects may obtain more accurate definitions of seed regions.

Conclusions

In this study, we used an SCN mapping approach to investigate eight large-scale networks in a large cohort of healthy participants. Our results show age-related changes in the topological organization of SCNs and provide insights into the normal aging process of the human brain.

Conflict of Interest Statement

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.

Acknowledgments

This work was supported by the National Science Foundation of China (Nos. 81171403, 81071212, and 30700182) and the State Key Laboratory of Software Development Environment (No. SKLSDE-2011ZX-10).

References

Allen, J. S., Bruss, J., Brown, C. K., and Damasio, H. (2005). Normal neuroanatomical variation due to age: the major lobes and a parcellation of the temporal region. Neurobiol. Aging 26, 1245–1260.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Andrews-Hanna, J. R., Snyder, A. Z., Vincent, J. L., Lustig, C., Head, D., Raichle, M. E., et al. (2007). Disruption of large-scale brain systems in advanced aging. Neuron 56, 924.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage 38, 95–113.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bakkour, A., Morris, J. C., and Dickerson, B. C. (2009). The cortical signature of prodromal AD: regional thinning predicts mild AD dementia. Neurology 72, 1048–1055.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bergfield, K. L., Hanson, K. D., Chen, K., Teipel, S. J., Hampel, H., Rapoport, S. I., et al. (2010). Age-related networks of regional covariance in MRI gray matter: reproducible multivariate patterns in healthy aging. Neuroimage 49, 1750–1759.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Biswal, B. B., Mennes, M., Zuo, X. N., Gohel, S., Kelly, C., Smith, S. M., et al. (2010). Toward discovery science of human brain function. Proc. Natl. Acad. Sci. U.S.A. 107, 4734–4739.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Cepeda, N. J., Kramer, A. F., and Gonzalez De Sather, J. (2001). Changes in executive control across the life span: examination of task-switching performance. Dev. Psychol. 37, 715.

Pubmed Abstract | Pubmed Full Text

Chao, L. L., Weisberg, J., and Martin, A. (2002). Experience-dependent modulation of category-related cortical activity. Cereb. Cortex 12, 545–551.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Chen, Z. J., He, Y., Rosa-Neto, P., Gong, G. L., and Evans, A. C. (2011). Age-related alterations in the modular organization of structural cortical network by using cortical thickness from MRI. Neuroimage 56, 235–245.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Damoiseaux, J., Beckmann, C., Arigita, E. J. S., Barkhof, F., Scheltens, P., Stam, C., et al. (2008). Reduced resting-state brain activity in the “default network” in normal aging. Cereb. Cortex 18, 1856–1864.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Demonet, J. F., Chollet, F., Ramsay, S., Cardebat, D., Nespoulous, J. L., Wise, R., et al. (1992). The anatomy of phonological and semantic processing in normal subjects. Brain 115, 1753–1768.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Draganski, B., Gaser, C., Busch, V., Schuierer, G., Bogdahn, U., and May, A. (2004). Neuroplasticity: changes in grey matter induced by training. Nature 427, 311–312.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Du, A. T., Schuff, N., Chao, L. L., Kornak, J., Jagust, W. J., Kramer, J. H., et al. (2006). Age effects on atrophy rates of entorhinal cortex and hippocampus. Neurobiol. Aging 27, 733–740.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Federmeier, K. D., Van Petten, C., Schwartz, T. J., and Kutas, M. (2003). Sounds, words, sentences: age-related changes across levels of language processing. Psychol. Aging 18, 858.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Ferrer, I., Blanco, R., Carulla, M., Condom, M., Alcantara, S., Olive, M., et al. (1995). Transforming growth factor-alpha immunoreactivity in the developing and adult brain. Neuroscience 66, 189–199.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fjell, A. M., Westlye, L. T., Amlien, I., Espeseth, T., Reinvang, I., Raz, N., et al. (2009). High consistency of regional cortical thinning in aging across multiple samples. Cereb. Cortex 19, 2001–2012.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Folstein, M. F., Folstein, S. E., and McHugh, P. R. (1975). “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J. Psychiatr. Res. 12, 189–198.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Fransson, P. (2005). Spontaneous low-frequency BOLD signal fluctuations: an fMRI investigation of the resting-state default mode of brain function hypothesis. Hum. Brain Mapp. 26, 15–29.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gong, G., He, Y., Chen, Z. J., and Evans, A. C. (2012). Convergence and divergence of thickness correlations with diffusion connections across the human cerebral cortex. Neuroimage 59, 1239–1248.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Gong, G., Rosa-Neto, P., Carbonell, F., Chen, Z. J., He, Y., and Evans, A. C. (2009). Age- and gender-related differences in the cortical anatomical network. J. Neurosci. 29, 15684–15693.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Good, C. D., Johnsrude, I. S., Ashburner, J., Henson, R. N., Friston, K. J., and Frackowiak, R. S. (2001). A voxel-based morphometric study of ageing in 465 normal adult human brains. Neuroimage 14, 21–36.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc. Natl. Acad. Sci. U.S.A. 100, 253.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Greicius, M. D., Supekar, K., Menon, V., and Dougherty, R. F. (2009). Resting-state functional connectivity reflects structural connectivity in the default mode network. Cereb. Cortex 19, 72–78.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

He, Y., Chen, Z. J., and Evans, A. C. (2007). Small-world anatomical networks in the human brain revealed by cortical thickness from MRI. Cereb. Cortex 17, 2407–2419.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kennedy, K. M., Erickson, K. I., Rodrigue, K. M., Voss, M. W., Colcombe, S. J., Kramer, A. F., et al. (2009). Age-related differences in regional brain volumes: a comparison of optimized voxel-based morphometry to manual volumetry. Neurobiol. Aging 30, 1657–1676.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Koch, W., Teipel, S., Mueller, S., Buerger, K., Bokde, A., Hampel, H., et al. (2010). Effects of aging on default mode network activity in resting state fMRI: does the method of analysis matter? Neuroimage 51, 280–287.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lerch, J. P., Worsley, K., Shaw, W. P., Greenstein, D. K., Lenroot, R. K., Giedd, J., et al. (2006). Mapping anatomical correlations across cerebral cortex (MACACC) using cortical thickness from MRI. Neuroimage 31, 993–1003.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Li, S., Xia, M., Pu, F., Li, D., Fan, Y., Niu, H., et al. (2011). Age-related changes in the surface morphology of the central sulcus. Neuroimage 58, 381–390.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Madden, D. J., Bennett, I. J., and Song, A. W. (2009). Cerebral white matter integrity and cognitive aging: contributions from diffusion tensor imaging. Neuropsychol. Rev. 19, 415–435.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Maguire, E. A., Frith, C. D., and Morris, R. (1999). The functional neuroanatomy of comprehension and memory: the importance of prior knowledge. Brain 122, 1839–1850.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Manjon, J. V., Coupe, P., Marti-Bonmati, L., Collins, D. L., and Robles, M. (2010). Adaptive non-local means denoising of MR images with spatially varying noise levels. J. Magn. Reson. Imaging 31, 192–203.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Marcus, D. S., Wang, T. H., Parker, J., Csernansky, J. G., Morris, J. C., and Buckner, R. L. (2007). Open Access Series of Imaging Studies (OASIS): cross-sectional MRI data in young, middle aged, nondemented, and demented older adults. J. Cogn. Neurosci. 19, 1498–1507.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mechelli, A., Crinion, J. T., Noppeney, U., O'Doherty, J., Ashburner, J., Frackowiak, R. S., et al. (2004). Neurolinguistics: structural plasticity in the bilingual brain. Nature 431:757. doi: 10.1038/431757a

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Mechelli, A., Friston, K. J., Frackowiak, R. S., and Price, C. J. (2005). Structural covariance in the human cortex. J. Neurosci. 25, 8303–8310.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Montembeault, M., Joubert, S., Doyon, J., Carrier, J., Gagnon, J. F., Monchi, O., et al. (2012). The impact of aging on gray matter structural covariance networks. Neuroimage 63, 754–759.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Morris, J. C. (1993). The Clinical Dementia Rating (CDR): current version and scoring rules. Neurology 43, 2412–2414.

Pubmed Abstract | Pubmed Full Text

Morris, J. C., Storandt, M., Miller, J. P., McKeel, D. W., Price, J. L., Rubin, E. H., et al. (2001). Mild cognitive impairment represents early-stage Alzheimer disease. Arch. Neurol. 58, 397–405.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Nosarti, C., Mechelli, A., Herrera, A., Walshe, M., Shergill, S. S., Murray, R. M., et al. (2010). Structural covariance in the cortex of very preterm adolescents: a voxel-based morphometry study. Hum. Brain Mapp. 32, 1615–1625.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Okada, K., and Hickok, G. (2006). Identification of lexical-phonological networks in the superior temporal sulcus using functional magnetic resonance imaging. Neuroreport 17, 1293–1296.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

O'Sullivan, M., Jones, D. K., Summers, P. E., Morris, R. G., Williams, S. C., and Markus, H. S. (2001). Evidence for cortical “disconnection” as a mechanism of age-related cognitive decline. Neurology 57, 632–638.

Pubmed Abstract | Pubmed Full Text

Pagani, E., Agosta, F., Rocca, M. A., Caputo, D., and Filippi, M. (2008). Voxel-based analysis derived from fractional anisotropy images of white matter volume changes with aging. Neuroimage 41, 657–667.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Paus, T., Petrides, M., Evans, A. C., and Meyer, E. (1993). Role of the human anterior cingulate cortex in the control of oculomotor, manual, and speech responses: a positron emission tomography study. J. Neurophysiol. 70, 453–469.

Pubmed Abstract | Pubmed Full Text

Persson, J., Lustig, C., Nelson, J. K., and Reuter-Lorenz, P. A. (2007). Age differences in deactivation: a link to cognitive control? J. Cogn. Neurosci. 19, 1021–1032.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Powell, H., Parker, G. J. M., Alexander, D. C., Symms, M. R., Boulby, P. A., Wheeler-Kingshott, C. A. M., et al. (2006). Hemispheric asymmetries in language-related pathways: a combined functional MRI and tractography study. Neuroimage 32, 388–399.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rajapakse, J. C., Giedd, J. N., and Rapoport, J. L. (1997). Statistical approach to segmentation of single-channel cerebral MR images. IEEE Trans. Med. Imaging 16, 176–186.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Raz, N., Lindenberger, U., Rodrigue, K. M., Kennedy, K. M., Head, D., Williamson, A., et al. (2005). Regional brain changes in aging healthy adults: general trends, individual differences and modifiers. Cereb. Cortex 15, 1676–1689.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Resnick, S. M., Pham, D. L., Kraut, M. A., Zonderman, A. B., and Davatzikos, C. (2003). Longitudinal magnetic resonance imaging studies of older adults: a shrinking brain. J. Neurosci. 23, 3295–3301.

Pubmed Abstract | Pubmed Full Text

Reuter-Lorenz, P. A., and Cappell, K. A. (2008). Neurocognitive aging and the compensation hypothesis. Curr. Dir. Psychol. Sci. 17, 177–182.

Salat, D. H., Greve, D. N., Pacheco, J. L., Quinn, B. T., Helmer, K. G., Buckner, R. L., et al. (2009). Regional white matter volume differences in nondemented aging and Alzheimer's disease. Neuroimage 44, 1247–1258.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Salat, D. H., Tuch, D. S., Greve, D. N., van der Kouwe, A. J., Hevelone, N. D., Zaleta, A. K., et al. (2005). Age-related alterations in white matter microstructure measured by diffusion tensor imaging. Neurobiol. Aging 26, 1215–1227.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sambataro, F., Murty, V. P., Callicott, J. H., Tan, H. Y., Das, S., Weinberger, D. R., et al. (2010). Age-related alterations in default mode network: impact on working memory performance. Neurobiol. Aging 31, 839.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Seeley, W. W., Crawford, R. K., Zhou, J., Miller, B. L., and Greicius, M. D. (2009). Neurodegenerative diseases target large-scale human brain networks. Neuron 62, 42–52.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Soriano-Mas, C., Harrison, B. J., Pujol, J., Lopez-Sola, M., Hernandez-Ribas, R., Alonso, P., et al. (2012). Structural covariance of the neostriatum with regional gray matter volumes. Brain Struct. Funct. doi: 10.1007/s00429-012-0422-5. [Epub ahead of print].

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sowell, E. R., Peterson, B. S., Thompson, P. M., Welcome, S. E., Henkenius, A. L., and Toga, A. W. (2003). Mapping cortical change across the human life span. Nat. Neurosci. 6, 309–315.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Takahashi, K., Yamaguchi, S., Kobayashi, S., and Yamamoto, Y. (2005). Effects of aging on regional cerebral blood flow assessed by using technetium Tc 99 m hexamethylpropyleneamine oxime single-photon emission tomography with 3D stereotactic surface projection analysis. AJNR Am. J. Neuroradiol. 26, 2005–2009.

Pubmed Abstract | Pubmed Full Text

Terribilli, D., Schaufelberger, M. S., Duran, F. L., Zanetti, M. V., Curiati, P. K., Menezes, P. R., et al. (2011). Age-related gray matter volume changes in the brain during non-elderly adulthood. Neurobiol. Aging 32, 354–368.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tohka, J., Zijdenbos, A., and Evans, A. (2004). Fast and robust parameter estimation for statistical partial volume models in brain MRI. Neuroimage 23, 84–97.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

West, R., Murphy, K. J., Armilio, M. L., Craik, F. I. M., and Stuss, D. T. (2002). Lapses of intention and performance variability reveal age-related increases in fluctuations of executive control. Brain Cogn. 49, 402–419.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wierenga, C. E., Benjamin, M., Gopinath, K., Perlstein, W. M., Leonard, C. M., Rothi, L. J. G., et al. (2008). Age-related changes in word retrieval: role of bilateral frontal and subcortical networks. Neurobiol. Aging 29, 436–451.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wu, J. T., Wu, H. Z., Yan, C. G., Chen, W. X., Zhang, H. Y., He, Y., et al. (2011a). Aging-related changes in the default mode network and its anti-correlated networks: a resting-state fMRI study. Neurosci. Lett. 504, 62–67.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wu, K., Taki, Y., Sato, K., Kinomura, S., Goto, R., Okada, K., et al. (2011b). Age-related changes in topological organization of structural brain networks in healthy individuals. Hum. Brain Mapp. 33, 552–568.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Xiang, H. D., Fonteijn, H. M., Norris, D. G., and Hagoort, P. (2010). Topographical functional connectivity pattern in the perisylvian language networks. Cereb. Cortex 20, 549.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zhu, W., Wen, W., He, Y., Xia, A., Anstey, K. J., and Sachdev, P. (2012). Changing topological patterns in normal aging using large-scale structural networks. Neurobiol. Aging 33, 899–913.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Zielinski, B. A., Gennatas, E. D., Zhou, J., and Seeley, W. W. (2010). Network-level structural covariance in the developing brain. Proc. Natl. Acad. Sci. U.S.A. 107, 18191–18196.

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Appendix

TABLE A1

Table A1. Contrast analysis of structural covariance network trajectories derived from contralateral seeds.

FIGURE A1

Figure A1. Trajectories of the structural covariance networks derived from contralateral seeds. The plots of voxel counts by group (P < 0.05, FWE) indicate a similar pattern to the ipsilateral seed, except for the auditory and speech networks, which show differences in their spatial extents. The y-axis represents the voxel number (×10⁴). Abbreviations: Calcs, calcarine sulcus; HG, Heschl's gyrus; PreCG, precentral gyrus; IFGo, inferior frontal gyrus, pars opercularis; TPole, temporal pole; ANG, angular gyrus; DLPFC, dorsolateral prefrontal cortex; FI, frontoinsular cortex; R, right; L, left; Y, young group; M, middle-aged group; O, old group.

FIGURE A2

Figure A2. Trajectories of the gray matter volume for each seed. The y-axis represents the mean gray matter volume for each age group. Abbreviations: Calcs, calcarine sulcus; HG, Heschl's gyrus; PreCG, precentral gyrus; IFGo, inferior frontal gyrus, pars opercularis; TPole, temporal pole; ANG, angular gyrus; DLPFC, dorsolateral prefrontal cortex; FI, frontoinsular cortex; R, right; L, left; Y, young group; M, middle-aged group; O, old group.