Gray matter in the brain, which is reduced as we age, is responsible for:

Introduction

Aging is associated with widespread cerebral grey-matter atrophy (Sigurdsson et al. 2012). Age-related loss of grey matter is more pronounced in men (up to − 0.70% per year) than in women (up to − 0.55% per year) (Alexander et al. 2005; Enzinger et al. 2005; Sigurdsson et al. 2012) and differs among brain regions (Fjell and Walhovd 2010). Hence, global measures of brain atrophy do not depict areas of increased vulnerability due to aging and disregard the notion of the brain being organized in networks maturing and dying together (Alexander-Bloch et al. 2013). Age-related atrophy is most pronounced in the hippocampus, caudate nucleus, association cortex, cerebellum, and the medial temporal lobe, while little volume loss is seen in other cortical regions such as the entorhinal cortex and the primary visual cortex (Salat et al. 2004; Alexander et al. 2005; Raz et al. 2005; Zielinski et al. 2010; Jiang et al. 2014). It has been recognized using cortical thickness networks that intra-individual differences in the structure of a brain region covary with the structure of other brain regions (structural covariance, SCN) (Lerch et al. 2006).

DuPre and Spreng (2017) showed that age-related disintegration follows these structural covariance networks. Indeed, studies examining the structural covariance have shown that disintegration in several canonical networks, such as the default mode network, the dorsal attention network, the fronto-parietal network, the somatomotor network, the ventral attention network, or a language-related semantic network are strongly associated with age, while others, such as a temporal network, an auditory network or cerebellar networks remained relatively preserved during the aging process (Montembeault et al. 2012; Hafkemeijer et al. 2014; Foster-Dingley et al. 2016; DuPre and Spreng 2017). A network involving areas for transmodal processing, including the lateral prefrontal cortex, the frontal eye field, the intraparietal cortex, the superior temporal sulcus, the posterior cingulate cortex, and the medial temporal lobe, has shown an inverted u-shaped age trajectory and vulnerability in different pathologies (i.e., Alzheimer’s disease and schizophrenia) and an association to memory performance (Douaud et al. 2014). Hence, a grey-matter network-dependent vulnerability creates a biological basis for cognitive deterioration during aging (Seeley et al. 2009).

Considering the association with cognition, the previous studies indeed found impaired age-related structural covariance to be associated with cognitive dysfunction (Brickman et al. 2007; Steffener et al. 2013; Spreng and Turner 2013; Douaud et al. 2014) and fine motor skills (Hoogendam et al. 2014). We here extend the previous work by examining network disintegration and its association with cognitive functioning and fine motor skills over a wide age range in a large community-dwelling cohort.

We identified those grey-matter SCNs that are affected by aging, determined their domain-specific association with cognitive function, and assessed if associations between age and cognitive function and fine motor skills are mediated by grey-matter network disintegration.

Methods

Compliance with ethical standards

The study was approved by the ethics committee of the Medical University of Graz, Austria. Informed consent was obtained from all individual participants included in the study. Disclosures of authors are quoted at the end of the manuscript.

Subjects and assessments

Overall, 257 non-demented subjects were included in the study. Two-hundred sixteen of them were elderly participants of the Austrian Stroke Prevention Study (ASPS), with a mean age of 68.2 ± 10.7 years, and 41 were younger individuals who have been invited to serve as healthy participants in a different study with a mean age of 26.9 ± 4.7 years. In total, there were 113 men and 144 women. All subjects underwent an identical imaging protocol. The 41 young subjects were included to broaden the age range of the ASPS from 39 to 87 years for the age stratum of 20–38 years (for age distribution, see Fig. 1). The ASPS is a single-center, prospective, follow-up study on the cerebral effects of vascular and genetic risk factors in the normal elderly population of Graz, Austria (Schmidt et al. 2003, 2005). Individuals were excluded from the ASPS study if they had a history of neuropsychiatric disease, including cerebrovascular attacks and dementia, or an abnormal neurologic examination determined on the basis of a structured clinical interview and physical and neurologic examinations. The same selection criteria were applied to the 41 subjects that were recruited to extend the range of the study down to 20 years of age.

Fig. 1

Gray matter in the brain, which is reduced as we age, is responsible for:

Age distribution of the total study sample (n = 257)

Full size image

All ASPS participants, but not the 41 younger subjects, underwent neuropsychological testing including tests on memory, executive functions, and fine motor skills. The tests employed have been widely used in the German speaking area and were always applied in the same order and under unchanged laboratory conditions. Intermediate memory recall and learning ability was assessed by the “Bäumler’s Lern- und Gedächtnistest” (LGT-3) (Bäumler 1974), a highly demanding paper–pencil procedure consisting of six subtests. Three subtests (word and digit association tasks, and story recall) screen for verbal memory, and two subtests (trail and design recall) screen for figural memory. The sum of weighted scores from these subtests and of an image recognition paradigm results in the total learning and memory performance score. The subtests were weighted based on the reliability scores provided in the test manual. The stimulus sets of the word association task (German–Turkish word pairs), the story (facts about construction of a library), and design recall (core symbol and frame), and the recognition paradigm (objects) consist of 20 items each. A trail in an abstracted city map serves as the trail recall test. These sets of stimuli were presented to the person being tested for 1 min. Two minutes were given for learning the 13 items of the digit association task (three-digit telephone numbers and names of extension holders). During a learning phase, the six sets of stimuli are subsequently presented to the person being tested. The recall phase starts immediately thereafter and follows the same order. The delay between presentation and recall for a given subtest ranges between 7 and 11 min. Executive functions were tested by the part B of the trail making test (Department 1944) and the digit span forward and backwards, which is part of the Wechsler adult intelligence scale, revised (Tewes 1991). Fine motor skills were evaluated by the Purdue pegboard test (Tiffin and Asher 1948). All scores were z-transformed for harmonization. The Mini Mental State Examination (MMSE) was used for the evaluation of the presence of Mild Cognitive Impairment (MCI). A subject was classified as mild cognitively impaired when the subject’s test score was below the 25th percentile of age and education specific norm data (Crum et al. 1993).

Image acquisition

Magnetic resonance imaging was performed on a 3T whole body scanner (TimTrio; Siemens Healthcare, Erlangen, Germany) with a 12-channel head coil. For each participant, a high-resolution T1-weighted three-dimensional anatomical image with magnetization preparation (MPRAGE) and whole brain coverage (TR = 1900 ms, TE = 2.19 ms, inversion time = 900 ms, flip angle = 9°, isotropic resolution of 1 mm) was acquired. For the identification of white matter hyperintensities, an axial T2-weighted FLAIR sequence (TR = 10,000 ms, TE = 69 ms, inversion time = 2500 ms, number of slices = 40, slice thickness = 3 mm, in-plane resolution = 0.86 mm × 0.86 mm) was used.

Image analysis

To ensure sufficient data quality, all T1-weighted images were visually checked to exclude potentially artefact afflicted scans. Imaging data analysis was performed using the FMRIB’s Software Library (FSL, 5.0.9, Oxford, UK) (Smith et al. 2004). For a detailed description of processing procedures, see (Hafkemeijer et al. 2014). In short, the following pre-processing steps were conducted: brain extraction from T1-weighted images using a semi-automated tool as implemented in FSL (Smith 2002), and tissue-type segmentation into grey matter, white matter, and cerebrospinal fluid, using a voxel-based morphometric analysis (Ashburner and Friston 2000). Again, visual inspection of all grey-matter segmented images was performed to ensure data quality. The individual grey-matter images were aligned to the grey-matter MNI152 standard space (Montreal Neurological Institute, Montreal, QC, Canada) (Jenkinson et al. 2002), followed by a non-linear registration (Andersson et al. 2007). The resulting images were averaged to create a study-specific grey-matter template. Then, all native grey-matter images were non-linearly re-registered to this study-specific template and “modulated” to correct for local expansion (or contraction) due to the non-linear component of the spatial transformation. The modulated grey-matter images were then smoothed with an isotropic Gaussian kernel with a sigma of 3 mm.

The modulated images of all 257 subjects were used as a four-dimensional data set on which an independent component analysis was applied using multivariate exploratory linear optimized decomposition into independent components (Beckmann et al. 2005). This procedure decomposes the signals into spatial component maps of maximal statistical independence (Beckmann and Smith 2004). When applied on grey-matter images, this method defines fully automatically spatial components based on the covariation of grey-matter intensities among subjects (i.e., SCNs), without a predefined region of interest. As there exists no consensus on the best number of components, we restricted the number to 20. A mixture model was used to assign significance to individual voxels within a spatial map, using a standard threshold level of 0.5 (Beckmann and Smith 2004). To obtain the integrity score of the SCNs, the four-dimensional data set of grey-matter images was spatially regressed against the 20 SCN probability maps using a general linear model approach integrated in FSL (Filippini et al. 2009). This method calculates beta values (positive and negative scores) for each network. The beta scores, representing measures for the network integrity, were used in all statistical analyses. Anatomical locations were determined using the Harvard–Oxford cortical, subcortical and the cerebellar atlases implemented in FSL. An equivalent approach was used in (Hafkemeijer et al. 2014, 2016; Foster-Dingley et al. 2016).

Global brain volume, normalized for subject head size, was estimated with SIENAX (Smith et al. 2001), part of FSL (Smith et al. 2004). SIENAX starts by extracting brain and skull images from the single whole-head input data (Smith 2002). The brain image is then affine-registered to MNI152 space (Jenkinson and Smith 2001) (using the skull image to determine the registration scaling); this is primarily to obtain the volumetric scaling factor, to be used as a normalization for head size. Next, tissue-type segmentation with partial volume estimation is carried out (Zhang et al. 2001) to calculate total volume of brain tissue (including separate estimates of volumes of grey matter, white matter, peripheral grey matter, and ventricular CSF).

Confounding variables

To adjust for potential confounders on cognition, the following variables were tested in simple regression analyses for their association with cognitive functions, and considered in the statistical models if p was < 0.1. Besides the variables age, sex, and education, also the variables hypertension [present yes/no, 22.6% (n = 58) of the subjects], lacunes [present yes/no, 10.9% (n = 28) of the subjects], white matter hyperintensity score (WMHs), and normalized global brain volume showed associations with memory, executive function, and fine motor skills. Visual rating of WMHs was rated according to our scheme (Fazekas et al. 1987, 1993) into absent [grade 0, 17.9% (n = 46) of the subjects], punctuate [grade 1, 36.1% (n = 93) of the subjects], early confluent [grade 2, 17.1% (n = 44) of the subjects], and confluent [grade 3, 10.1% (n = 26) of the subjects]. No association with cognition was found for diabetes [present yes/no, 4.3% (n = 11) of the subjects], cholesterol and smoking (yes/no, M = 203.7 mg/dl, SD = 41.5).

Statistical analysis

For statistical analyses, IBM SPSS Statistics Version 22, IBM Corp., USA and R (version 3.2.4) (R-Core-Team 2016) was used. We examined the association between age and grey-matter network disintegration adjusted for sex and education in all 257 subjects (SCN (Y) ~ [age + sex + education]). We aimed to assess the predictive value of all structural covariance networks, risk factors, age-related brain abnormalities (WMH score, lacunes) as well as age, sex, and education for memory, executive function, and fine motor skills. We did not perform multiple linear regression to avoid the risks of overfitting, overadjustment, and ultimately biased estimation (Ranucci et al. 2010; Yoo et al. 2014) due to multicollinearity (intercorrelation of variables). Instead, we used random forest regression, which assesses the explanatory power of variables while accounting for all other variables. One major independent variables, even in the presence of complex interactions and multicollinearity when applying conditional inference trees (Strobl et al. 2009). We calculated 1001 conditional inference trees with unbiased variable selection using the standard parameters (5 randomly preselected variables for each split, unbiased resampling scheme) using the R package ‘party’ (version 1.0-25) (Strobl et al. 2007). From these trees, we next calculated a conditional permutation importance (following the permutation principle of the ‘mean decrease in accuracy’ importance measure) (Strobl et al. 2008) for each variable together with a 95% confidence interval from 100 repetitions. To ensure that the results were not driven by subjects with MCI, the random forest analysis was repeated excluding subjects classified as MCI.

To test if age effects on cognition or on fine motor skills are mediated by grey-matter disintegration, we applied simple mediation models for estimating indirect effect size (Hayes 2013). Mediation was evaluated separately in those grey-matter SCNs that disintegrated with aging and also showed a significant association with cognitive impairment or fine motor skills in the random forest model. In total, eight networks fulfilled these requirements. As some networks were associated with more than one cognitive domain, 25 mediation analyses were performed. Mediator effect size and 95% confidence intervals were estimated using a bootstrap-based method developed by Preacher and Hayes (Preacher and Hayes 2008). If the 95% confidence interval of the indirect effect does not contain 0, a significant mediation effect is probable, whereas no mediation is present if 0 is included in the 95% confidence interval. Mediation analyses was performed with the macro PROCESS [http://www.afhayes.com (Hayes 2013)] implemented in SPSS.

Results

Structural covariance networks and aging

Independent component analysis revealed 14 supratentorial covariance networks (Fig. 2) and six infratentorial networks (Fig. 3) to be present in community-dwelling subjects free of stroke and dementia. The areas included in each of the supratentorial networks are described in Table 1. Those of the infratentorial networks are described in Table 2.

Fig. 2

Gray matter in the brain, which is reduced as we age, is responsible for:

Supratentorial structural covariance networks of non-demented subjects. For areas included in each of these networks, refer to Table 1. Box–whisker plots indicate grey-matter integrity per decade, and p value indicates significance of age as predictor of network disintegration (adjusted for sex and education, and corrected for multiple comparisons; FDR, q < 0.05, significant networks are marked with an asterisk). A linear decline of grey-matter integrity with aging was found in networks (ah), in one network (i) the association was u-shaped, and in another (j), it was inversely u-shaped. The number of subjects within the age ranges is indicated in parenthesis

Full size image

Fig. 3

Gray matter in the brain, which is reduced as we age, is responsible for:

Infratentorial structural covariance networks of non-demented subjects. For areas included in each of these networks, refer to Table 2. Box–whisker plots indicate grey-matter integrity per decade, and p value indicates significance of age as predictor of network disintegration (adjusted for sex, and education, and corrected for multiple comparisons; FDR, q = 0.05, significant networks are marked with an asterisk). Networks ux show a reduction of grey-matter network integrity with increasing age. The number of subjects within the age ranges is indicated in parenthesis

Full size image

Table 1 Areas included in supratentorial structural covariance networks (cluster size k > 1000 voxels)

Full size table

Table 2 Areas included in infratentorial structural covariance networks (cluster size k > 1000 voxels)

Full size table

Ten supratentorial (Fig. 2a–j) and four infratentorial networks (Fig. 3u–x) lost integrity with advancing age. The association with aging was linear in all but two networks (Fig. 2i, j).

Structural covariance network disintegration and cognition

As can be seen from Fig. 4, in our study, population advancing age was associated with performance decrease in all cognitive domains.

Fig. 4

Gray matter in the brain, which is reduced as we age, is responsible for:

Cognitive performance changes with aging per decade in memory, executive functions, and fine motor skills. All R2’s are based on a hierarchical linear regression (step 1: sex, education; step 2: age)

Full size image

Correlation matrices of the structural covariance networks revealed a high intercorrelation among the networks (data not shown). Therefore, to assess the contribution of each network to memory, executive function, and fine motor skills, while accounting for intercorrelations (multicollinearity), we applied random forest regression and calculated the conditional variable importance (VI). Education and age had the highest variable importance (Fig. 5; Table 3). SCNs which contributed independently to memory were the secondary somatosensory network (Fig. 2h), the temporal network (Fig. 2e), the limbic network (Fig. 2d), the fronto-parietal network (Fig. 2g), the fronto-occipital network (Fig. 2i), and the cuneal network (Figs. 2a, 5; Table 3). Those contributing to executive function were the secondary somatosensory network (Fig. 2h), the temporal network (Fig. 2e), the sensorimotor network (Fig. 2j), the limbic network (Fig. 2d), the fronto-parietal network (Fig. 2g), and a cerebellar network (Figs. 3v, 5; Table 3). The secondary somatosensory network (Fig. 2h), the temporal network (Fig. 2e), the sensorimotor network (Fig. 2j), the limbic network (Fig. 2d), the fronto-parietal network (Fig. 2g), a cerebellar network (Fig. 3v), and the fronto-occipital network contributed independently to fine motor skills of study participants (Figs. 2i, 5; Table 3). The WMH score had predictive value only for fine motor skills. Normalized global brain volume, hypertension, and lacunes were not associated with any outcome (Fig. 5).

Fig. 5

Gray matter in the brain, which is reduced as we age, is responsible for:

Descriptive ranking of variable importance (VI) for the prediction of cognitive and motor skills, as determined by random forest regression. Only networks changing integrity with age were included in the model. In addition, demographical variables (education, age, and sex), age-related brain abnormalities (WMH score, lacunes), normalized global brain volumes (nGBV), and hypertension were included. Spheres and vertical lines indicate VI mean and 95% confidence interval. Only if the confidence interval did not overlap with zero, a variable was considered to have a significant contribution in predicting the performance (filled dot). Education and age showed the highest variable importance in all domains. Additional association with the domains were found for the secondary somatosensory network (h), the temporal network (e), the sensorimotor network (j), the limbic network (d), the fronto-parietal network (g), a cerebellar network (v), the fronto-occipital network (i), and the cuneal network (a)

Full size image

Table 3 Conditional variable importance mean (point estimator) and 95% confidence interval of structural covariance networks, age-related brain abnormalities, risk factors and demographics for cognitive outcomes and fine motor skills related to Fig. 5

Full size table

By means of MMSE norm scores, 34 subjects were classified as being mild cognitively impaired. Repetition of the random forest analyses without these subjects revealed almost unchanged results. Equivalently to the above results, age, sex, and education had the utmost association with cognitive and motor domains. Similar, memory was associated with the limbic network (Fig. 2d), the temporal network (Fig. 2e), the fronto-parietal network (Fig. 2g), the secondary somatosensory network (Fig. 2h), and additionally two cerebellar networks (Fig. 3v, w). SCNs contributing to executive functions were the temporal network (Fig. 2e), the sensorimotor network (Fig. 2j), the limbic network (Fig. 2d), the fronto-parietal network (Fig. 2g), and a cerebellar network (Fig. 3v). The normalized global brain volume and all but the fronto-occipital network of the total sample were associated with fine motor skills.

Mediation analyses

Mediation analyses were performed for each SCN disintegrating with age and showing an association with cognition or fine motor skills. Of 25 mediation analyses, in 13 cases, the effect of age on cognition or fine motor skills was mediated by an SCN. Seven of the tested SCNs were among these mediators (Table 4). A negative association of age on cognition and fine motor skills, ranging between − 0.1612 and − 0.5196 (total effect), was found. The total effect gets reduced by the indirect effect [− 0.0030, 0.0206] of the mediator, i.e., the SCN, reducing the effect of age on cognition and fine motor skills [direct effect, − 0.0275, − 0.4813]. The ratio of the indirect to the total effect ranges between 7.38 and 23.48%.

Table 4 Mediation models assessing the effect of grey-matter disintegration as measured by the SCNs on the relationship between age and cognitive performance and fine motor skills

Full size table

Discussion

In this study, we identified structural covariance networks that disintegrate with increasing age and related them to cognitive function independently of risk factors and age-related brain abnormalities including vascular lesions and brain atrophy. Of 14 networks that lost integrity with aging, eight were related to either cognitive or motor function, of which seven mediated the effect between age and cognition and motor function. With the exception of two networks, the loss of network integrity was linear with advancing age. Besides age and education, disintegration of the temporal, limbic, fronto-parietal, fronto-occipital, sensorimotor, secondary somatosensory, cuneal and a cerebellar network showed strongest association with cognitive or motor function. Importantly, these results remained largely unchanged when excluding subjects with Mild Cognitive Impairment (MCI) potentially showing incipient neuropathological alterations. The grey-matter networks identified in our investigation showed an overlap with covariance networks described by prior literature (Hafkemeijer et al. 2014, 2016), and also partly overlapped with intrinsic functional networks identified with BOLD-contrast imaging (Smith et al. 2009). This is not implausible, since direct anatomical connections between areas are related to functional connectivity but not mandatory, i.e., functional connectivity can arise in the absence of structural connectivity (Alexander-Bloch et al. 2013). Moreover, the results of the mediation analyses confirmed the partial mediating role of grey-matter disintegration in the relation between age and cognition and fine motor skills, suggesting a neuronal basis accounting for at least part of the association between age and cognition and fine motor skills. The mediator could account for up to roughly 23% of the total effect (Table 4, ratio). The results of the current study are in line with the view that growth and degeneration of the cortex occurs at the level of networks, rather than in a region specific manner (Alexander-Bloch et al. 2013; Hafkemeijer et al. 2014). According to our data network disintegration is seen already above the age of 30, a finding in keeping with previous literature reporting grey-matter loss soon after adolescence (Courchesne et al. 2000; Giedd 2004; Alexander-Bloch et al. 2013).

There exists little information on age-related cortical network degeneration so far. Previous investigations used subprofile scaling models (Brickman et al. 2007; Steffener et al. 2013) or partial least squares (DuPre and Spreng 2017) and reported that networks that best differentiate between younger and older subjects were neurocognitive networks linked to attention (Brickman et al. 2007; Steffener et al. 2013), language (Brickman et al. 2007), memory (Brickman et al. 2007; Steffener et al. 2013), executive functions (Brickman et al. 2007), and fluid abilities (Steffener et al. 2013). In our study, on top of age, sex, and education, the grey-matter covariance networks presented up to 10.08% of importance for memory, 3.11% for executive function, and 15.08% for fine motor skills. Normalized global brain volume, lacunes, and hypertension did show an association with cognitive and motor abilities, and WMH score was associated with fine motor skills only. Of the eight networks that were identified to be important for memory, executive function, and fine motor skill in the random forest analyses, seven networks revealed a partly mediating effect between age and cognitive and motor outcome. Hence, the disintegration of grey-matter covariance networks constitutes an additional potential factor for cognitive decline in the aging brain independent of normalized global brain volume or other confounders. Determining the causes for this selective vulnerability of disintegration of certain networks with aging leads to the developmental and maturing processes underlying these covariance networks. Coordinated neurodevelopment has been suggested to constitute the basis of morphological covariance which may be induced in different ways. Provisional theories suggest that the correlation of phenotypic traits could be evoked by shared genetic influences, common environmental factors, inductive signaling from one developing tissue to another or simultaneous exposure to signals from third party, timing of development or the sharing of a developmental precursor (Riska 1986; Alexander-Bloch et al. 2013). Hence, coordinated neurodevelopment might form the scaffold of structural covariance. Preliminary imaging results underpin this hypothesis. A strong covariance has been identified for growth and disintegration in volume of networks (DuPre and Spreng 2017). In a longitudinal study, it was shown that areas showing structural covariance were also correlated in their rate of change considering cortical thickness (Alexander-Bloch et al. 2013).

Increasing age differentially affected the integrity of networks determined in our study. While most networks showed a linear disintegration with age, some remained stable over the entire age range between 20 and 87 years, and two had an u-shaped trajectory. Stable networks were found to connect the superior temporal gyri, the cerebellum and the hippocampus (Fig. 2k), frontal areas with the occipital pole (Fig. 2l), the inferior temporal gyrus and the cerebellum (Fig. 2m) and the supplementary motor area with the insular, the precuneus, the temporal occipital fusiform cortex, and the cerebellum (Fig. 2n). Networks with decreasing connectivity included a cuneal network, mainly showing covariance between the cuneus and the paracingulate gyrus, an anterior cingulate network, with associations between the anterior cingulate cortex, the occipital fusiform gyrus and the insular cortex, a precuneal network, showing associations between the precuneus and the cerebellum, a limbic network, connecting subcortical and cortical areas, a temporal network, showing major associations between the inferior temporal gyrus, the temporal pole, and the paracingulate gyrus, a subcortical network, mainly comprising the amygdala, a fronto-parietal network, connecting the frontal poles with the lateral occipital cortices, a secondary sensorimotor network, including both parietal opercular cortices and the posterior cingulate gyrus (Fig. 2a–h), and finally four cerebellar networks (Fig. 3u–x). Notably, two networks comprising the superior frontal gyrus and the thalamus (Fig. 2i) and another showing connection between the pre- and postcentral gyri, the cerebellum, and the hippocampi (Fig. 2j) showed (inverted) u-shaped trajectories. While the inverted u-shape in network j, which is primarily a sensorimotor network, indicates that increasing network integration in this important functional network may occur up to the age of 60, it is difficult to explain the u-shaped trajectory of network i. We cannot exclude that this was a chance finding, although all statistics in our investigation were corrected for multiple comparisons. It needs to be emphasized that the results of current cross-sectional investigation need to be replicated by longitudinal studies which ideally follow the same individuals over their life span.

Our results revealed an association between network disintegration and higher order cognitive functions and fine motor skills, respectively. Memory has been associated with the grey-matter integrity of the temporal, secondary somatosensory, limbic, cuneal and the fronto-occipital network. The strongest association with executive function was seen for the fronto-parietal, temporal, sensorimotor, limbic, secondary somatosensory, and the cerebellar network. The sensorimotor, secondary somatosensory, temporal, limbic, fronto-occipital, fronto-parietal, and cerebellar network were identified as being the most important determinants for fine motor skills. The temporal, limbic, fronto-occipital, and the secondary somatosensory network have been associated with more than one cognitive domain or fine motor skill, respectively.

Four cerebellar networks showed an association with age, but only one with fine motor skills. This was somewhat surprising, since the cerebellum has repeatedly been shown to be involved in cognitive and motor function (Stoodley 2012). In addition, this is in contrast to a prior study on subjects > 75 years reporting no association between cerebellar grey-matter networks and age, but an association with psychomotor speed (Hafkemeijer et al. 2014; Foster-Dingley et al. 2016). An explanation for this discrepancy might be found in our networks which encompass widespread cerebellar regions instead of rather small circumscribed areas. A prior study using functional imaging revealed a specialization of cerebellar subfields depending on the task performed (Stoodley et al. 2011) highlighting narrow circumscribed cerebellar areas. It is, therefore, thinkable that our networks over-represent a specific domain.

We here examined grey-matter network changes in a cross-sectional design. This approach does not allow inferences on intra-individual trajectories of grey-matter network disintegration over time limiting the interpretability of individual association between risk factors, network disintegration, and the effect on cognition. Longitudinal designs are considered optimal to examine intra-individual changes, but imaging studies spanning long time periods also face problems such as advancements in sequence technology or even hardware changes hampering comparability.

The inclusion of individuals between 20 and 38 years of age allowed to broaden the age range of our study and to assess covariance of networks from post-adolescence up to the ninth decade of life. However, these subjects have not been cognitively tested and led to a three-modal distribution of age with an underrepresentation of subjects in the age range between 35 and 45 years. This group of subjects was added from a different study with the same imaging protocol, but without neuropsychological assessment. We are aware of this drawback, but inclusion of this group of individuals allowed to broaden the age range of our study and to assess covariance of networks from post-adolescence up to the ninth decade of life. Third, the number of independent components was chosen arbitrarily. In functional imaging, most frequently, a number between 10 and 20 are chosen (Barkhof et al. 2014). To obtain a better sub-network segmentation, we chose a higher number. Finally, the mediation analyses were not corrected for multiple comparisons as with using confidence intervals for significant result interpretation, and no p values are calculated by the software PROCESS by Andrew Hayes. Hence, we assumed that the null is always true, and therefore, the probability of not making any type I error would be 0.9513 = 0.51, which means that 1 of the 13 analyses has a probability of 0.49 of being a false positive. Thus, one of the significant 13 results of the mediation analyses has a 49% chance of being incorrect.

Our study cannot determine whether the observed association between network disintegration and cognitive impairment with advancing age is at least partly due to evolving neurodegenerative disease. Future longitudinal studies will have to determine if any of the identified covariance networks have the prognostic potential to identify individuals with a high risk for developing mild cognitive impairment or conversion to dementia beyond what can be expected for the measurement of regional brain atrophy alone. They will also have to determine as to how the prognostic value of structural covariance networks compares to other dementia biomarkers such as CSF amyloid and tau or amyloid PET.

References

  • Alexander G, Chen K, Merkley T et al (2005) Regional network of magnetic resonance imaging gray matter volume in healthy aging. Neuroreport 17:951–956

    Article  Google Scholar 

  • Alexander-Bloch A, Giedd JN, Bullmore E (2013) Imaging structural co-variance between human brain regions. Nat Rev Neurosci 14:322–336

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • Andersson J, Jenkinson M, Smith S (2007) Non-linear registration aka Spatial normalization. FMRIB Tech Rep TR07JA2

  • Ashburner J, Friston KJ (2000) Voxel-based morphometry—the Methods. Neuroimage 11:805–821

    Article  PubMed  CAS  Google Scholar 

  • Barkhof F, Haller S, Rombouts SARB. (2014) Resting-state functional MR imaging: a new window to the brain. Radiology 272:29–49

    Article  PubMed  Google Scholar 

  • Bäumler G (1974) Lern und Gedächtnistest (LGT). Dr CJ Hogrefe Verlag für Psychologie, Göttingern

    Google Scholar 

  • Beckmann C, Smith S (2004) Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans Med Imaging 23:137–152

    Article  PubMed  Google Scholar 

  • Beckmann CF, DeLuca M, Devlin JT, Smith SM (2005) Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond B Biol Sci 360:1001–1013

    Article  PubMed  PubMed Central  Google Scholar 

  • Brickman AM, Habeck C, Zarahn E et al (2007) Structural MRI covariance patterns associated with normal aging and neuropsychological functioning. Neurobiol Aging 28:284–295

    Article  PubMed  Google Scholar 

  • Courchesne E, Chisum HJ, Townsend J et al (2000) Normal brain development and aging: quantitative analysis at in vivo MR imaging in healthy volunteers. Radiology 216:672–682

    Article  PubMed  CAS  Google Scholar 

  • Crum RM, Anthony JC, Bassett SS, Folstein MF (1993) Population-based norms for the mini-mental state examination by age and educational level. JAMA 269:2386–2391

    Article  PubMed  CAS  Google Scholar 

  • Department USW (1944) Army individual test battery manual of directions and scoring. Adjutant General’s Office, Washington, DC

    Google Scholar 

  • Douaud G, Groves AR, Tamnes CK et al (2014) A common brain network links development, aging, and vulnerability to disease. Proc Natl Acad Sci 111:17648–17653

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • DuPre E, Spreng RN (2017) Structural covariance networks across the life span, from 6 to 94 years of age. Netw Neurosci 1:302–323

    Article  PubMed  PubMed Central  Google Scholar 

  • Enzinger C, Fazekas F, Matthews PM et al (2005) Risk factors for progression of brain atrophy in aging: 6-year follow-up of normal subjects. Neurology 64:1704–1711

    Article  PubMed  CAS  Google Scholar 

  • Fazekas F, Chawluk J, Alavi A et al (1987) MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging. AJR 149:351–356

    Article  PubMed  CAS  Google Scholar 

  • Fazekas F, Kleinert R, Offenbacher H et al (1993) Pathologic correlates of incidental MRI white matter signal hyperintensities. Neurology 43:1683–1689

    Article  PubMed  CAS  Google Scholar 

  • Filippini N, MacIntosh BJ, Hough MG et al (2009) Distinct patterns of brain activity in young carriers of the APOE-epsilon4 allele. Proc Natl Acad Sci USA 106:7209–7214

    Article  PubMed  PubMed Central  Google Scholar 

  • Fjell AM, Walhovd KB (2010) Structural brain changes in aging: courses, causes and cognitive consequences. Rev Neurosci 21:187–221

    Article  PubMed  Google Scholar 

  • Foster-Dingley JC, Hafkemeijer A, van den Berg-Huysmans AA et al (2016) Structural covariance networks and their association with age, features of cerebral small-vessel disease, and cognitive functioning in older persons. Brain Connect 6:681–690

    Article  PubMed  Google Scholar 

  • Giedd JN (2004) Structural Magnetic Resonance Imaging of the Adolescent Brain. Ann N Y Acad Sci 1021:77–85

    Article  PubMed  Google Scholar 

  • Hafkemeijer A, Altmann-Schneider I, de Craen AJM et al (2014) Associations between age and gray matter volume in anatomical brain networks in middle-aged to older adults. Aging Cell 13:1068–1074

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • Hafkemeijer A, Möller C, Dopper E et al (2016) Differences in structural covariance brain networks between behavioral variant frontotemporal dementia and Alzheimer’s disease. Hum Brain Mapp 37:978–988

    Article  PubMed  Google Scholar 

  • Hayes AF (2013) Introduction to mediation, moderation, and conditional process analysis: a regression-based approach. The Guilford Press, New York

    Google Scholar 

  • Hoogendam YY, van der Lijn F, Vernooij MW et al (2014) Older age relates to worsening of fine motor skills: a population-based study of middle-aged and elderly persons. Front Aging Neurosci 6:259

    Article  PubMed  PubMed Central  Google Scholar 

  • Jenkinson M, Smith S (2001) A global optimisation method for robust affine registration of brain images. Med Image Anal 5:143–156

    Article  PubMed  CAS  Google Scholar 

  • Jenkinson M, Bannister P, Brady M, Smith S (2002) Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17:825–841

    Article  PubMed  Google Scholar 

  • Jiang J, Sachdev P, Lipnicki D et al (2014) A longitudinal study of brain atrophy over two years in community-dwelling older individuals. Neuroimage 86:203–211

    Article  PubMed  CAS  Google Scholar 

  • Lerch J, Worsley K, Shaw W et al (2006) Mapping anatomical correlations across cerebral cortex (MACACC) using cortical thickness from MRI. Neuroimage 31:993

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  • Preacher KJ, Hayes AF (2008) Asymptotic and resampling strategies for assessing and comparing indirect effects in multiple mediator models. Behav Res Methods 40:879–891

    Article  PubMed  Google Scholar 

  • Ranucci M, Castelvecchio S, Menicanti L et al (2010) Accuracy, calibration and clinical performance of the EuroSCORE: can we reduce the number of variables? Eur J Cardio-Thoracic Surg 37:724–729

    Article  Google Scholar 

  • Raz N, Lindenberger U, Rodrigue KM et al (2005) Regional brain changes in aging healthy adults: general trends, individual differences and modifiers. Cereb Cortex 15:1676–1689

    Article  PubMed  Google Scholar 

  • R-Core-Team (2016) A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna

    Google Scholar 

  • Riska B (1986) Some models for development, growth, and morphometric correlation. Evolution 40:1303–1311

    Article  PubMed  Google Scholar 

  • Salat D, Buckner R, Snyder A et al (2004) Thinning of the cerebral cortex in aging. Cereb Cortex 14:721–730

    Article  PubMed  Google Scholar 

  • Schmidt R, Enzinger C, Ropele S et al (2003) Progression of cerebral white matter lesions: 6-year results of the Austrian stroke prevention study. Lancet 361:2046–2048

    Article  PubMed  Google Scholar 

  • Schmidt R, Ropele S, Enzinger C et al (2005) White matter lesion progression, brain atrophy, and cognitive decline: the Austrian stroke prevention study. Ann Neurol 58:610–616

    Article  PubMed  Google Scholar 

  • Seeley WW, Crawford RK, Zhou J et al (2009) Neurodegenerative diseases target large-scale human brain networks. Neuron 62:42–52

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • Sigurdsson S, Aspelund T, Forsberg L et al (2012) Brain tissue volumes in the general population of the elderly: the AGES-Reykjavik study. Neuroimage 59:3862–3870

    Article  PubMed  Google Scholar 

  • Smith SM (2002) Fast robust automated brain extraction. Hum Brain Mapp 17:143–155

    Article  PubMed  PubMed Central  Google Scholar 

  • Smith SM, De Stefano N, Jenkinson M, Matthews PM (2001) Normalized accurate measurement of longitudinal brain change. J Comput Assist Tomogr 25:466–475

    Article  PubMed  CAS  Google Scholar 

  • Smith SM, Jenkinson M, Woolrich MW et al (2004) Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23(Suppl 1):S208–S219

    Article  PubMed  Google Scholar 

  • Smith SM, Fox PT, Miller KL et al (2009) Correspondence of the brain’s functional architecture during activation and rest. Proc Natl Acad Sci USA 106:13040–13045

    Article  PubMed  PubMed Central  Google Scholar 

  • Spreng RN, Turner GR (2013) Structural covariance of the default network in healthy and pathological aging. J Neurosci 33:15226–15234

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  • Steffener J, Brickman A, Habeck C et al (2013) Cerebral blood flow and gray matter volume covariance patterns of cognition in aging. Hum Brain Mapp 34:3267–3279

    Article  PubMed  Google Scholar 

  • Stoodley CJ (2012) The cerebellum and cognition: evidence from functional imaging studies. Cerebellum 11:352–365

    Article  PubMed  Google Scholar 

  • Stoodley CJ, Valera EM, Schmahmann JD (2011) Functional topography of the cerebellum for motor and cognitive tasks: an fMRI study. Neuroimage 59:1560–1570

    Article  PubMed  Google Scholar 

  • Strobl C, Boulesteix A-L, Zeileis A, Hothorn T (2007) Bias in random forest variable importance measures: illustrations, sources and a solution. BMC Bioinform 8:25

    Article  CAS  Google Scholar 

  • Strobl C, Boulesteix A-L, Kneib T et al (2008) Conditional variable importance for random forests. BMC Bioinform 9:307

    Article  CAS  Google Scholar 

  • Strobl C, Malley J, Tutz G (2009) An introduction to recursive partitioning: rationale, application, and characteristics of classification and regression trees, bagging, and random forests. Psychol Methods 14:323–348

    Article  PubMed  PubMed Central  Google Scholar 

  • Tewes U (1991) Hamburg–Wechsler Intelligenztest für Erwachsene, Revised edn. Verlag Hans Huber, Bern

    Google Scholar 

  • Tiffin J, Asher EJ (1948) The Purdue Pegboard: norms and studies of reliability and validity. J Appl Psychol 32(3):234

    Article  PubMed  CAS  Google Scholar 

  • Yoo W, Mayberry R, Bae S et al (2014) A study of effects of multicollinearity in the multivariable analysis. Int J Appl Sci Technol 4:9–19

    PubMed  PubMed Central  Google Scholar 

  • Zhang Y, Brady M, Smith S (2001) Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging 20:45–57

    Article  PubMed  CAS  Google Scholar 

  • Zielinski BA, Gennatas ED, Zhou J, Seeley WW (2010) Network-level structural covariance in the developing brain. Proc Natl Acad Sci USA 107:18191–18196

    Article  PubMed  PubMed Central  Google Scholar 

Download references

What happens if grey matter is reduced?

Since the neuronal cells of the grey matter are constantly working, they require a high supply of oxygen to function efficiently. Therefore when the grey matter does not have access to oxygen, the cells will begin to die, leading to possible irreversible brain damage and loss of function.

What happens to grey matter as we age?

Gray matter volume decreases with age in networks containing subcortical structures, sensorimotor structures, posterior, and anterior cingulate cortices. Gray matter volume in temporal, auditory, and cerebellar networks remains relatively unaffected with advancing age.

Why does grey matter decrease with age?

Gray matter volume decline after middle age (about 40 years), however, is thought to primarily reflect degenerative changes in the gray matter (Sowell et al., 2003), due to shrinkage or loss of neurons (Terry et al., 1987) and loss of dendritic arborization (Jacobs et al., 1997).

What does grey matter in the brain do?

The grey matter serves to process information in the brain. Structures within the grey matter process signals generated in the sensory organs or other areas of the grey matter. This tissue directs sensory (motor) stimuli to nerve cells in the central nervous system where synapses induce a response to the stimuli.