Chandan Shah1†
Jia Liu1†Peilin Lv2Huaiqiang Sun1Yuan Xiao1Jieke Liu1
Youjin Zhao1Wenjing Zhang1Li Yao1Qiyong Gong1
Su Lui1*- 1Huaxi MR Research Center, Department of Radiology, West China Hospital of Sichuan University, Chengdu, China
- 2Department of Anesthesiology, West China Hospital of Sichuan University, Chengdu, China
Introduction: There are still uncertainties about the true nature of age related changes in topological properties of the brain functional network and its structural connectivity during various developmental stages. In this cross- sectional study, we investigated the effects of age and its relationship with regional nodal properties of the functional brain network and white matter integrity.
Method: DTI and fMRI data were acquired from 458 healthy Chinese participants ranging from age 8 to 81 years. Tractography was conducted on the DTI data using FSL. Graph Theory analyses were conducted on the functional data yielding topological properties of the functional network using SPM and GRETNA toolbox. Two multiple regressions were performed to investigate the effects of age on nodal topological properties of the functional brain network and white matter integrity.
Result: For the functional studies, we observed that regional nodal characteristics such as node betweenness were decreased while node degree and node efficiency was increased in relation to increasing age. Perversely, we observed that the relationship between nodal topological properties and fasciculus structures were primarily positive for nodal betweenness but negative for nodal degree and nodal efficiency. Decrease in functional nodal betweenness was primarily located in superior frontal lobe, right occipital lobe and the global hubs. These brain regions also had both direct and indirect anatomical relationships with the 14 fiber bundles. A linear age related decreases in the Fractional anisotropy (FA) value was found in the callosum forceps minor.
Conclusion: These results suggests that age related differences were more pronounced in the functional than in structural measure indicating these measures do not have direct one-to-one mapping. Our study also indicates that the fiber bundles with longer fibers exhibited a more pronounced effect on the properties of functional network.
Introduction
Aging is considered as one of the greatest risk factors for neurodegenerative diseases where there is progressive loss of either structure or function of the neurons (Morrison and Hof, 1997). It is well-stipulated that marked transformations in white matter structure occur as early as third trimester (Ball et al., 2014) which persists throughout childhood and adolescence (Lebel et al., 2008; Hagmann et al., 2010; Tamnes et al., 2010; Dubois et al., 2014; van den Heuvel et al., 2015; Wierenga et al., 2016). These transformations reflects an increase in axonal conduction speed (Baumann and Pham-Dinh, 2001), thereby resulting in improvement of information transfer (van der Knaap et al., 1991). Similarly, functional studies have also investigated the effects of age on functional brain connectivity showing older adults to have lower functional connectivity between regions of the default mode network compared to younger adults (Ferreira and Busatto, 2013; Dennis and Thompson, 2014). However, little is known about these transformations beyond the young adults until the sixth or seventh decades of life, which is more vulnerable to the effects neurodegenerative diseases that show its effects abruptly after a possible long preclinical period (Grady et al., 1988; Berg et al., 1992; Morris et al., 1993). Although recent advances in network analysis have provided new insights into both structural and functional connectivity pattern, which is considered as the core of the brain activity (Sporns et al., 2005; Bullmore and Bassett, 2011; Sporns, 2012), it is still not fully understood about the relationships between them across lifespan, the evaluation of which might help us elucidate the possible mechanisms behind age related neurodegenerative diseases.
Among the broad range of network analysis approaches, graph theory approach is considered one of the favorable as it applies to both structural and functional connectivity, which is a natural framework for the exact mathematical representation of complex networks. In graph theory, a network is defined as a set of nodes with edges between them. Some frequently used network measures include: global network properties such as small-world properties [clustering coefficient (Cp) and characteristics path length (Lp)]; efficiency metrics (global efficiency and local efficiency); and regional nodal properties (nodal betweenness, nodal degree, and nodal efficiency) (Bullmore and Sporns, 2009). Similarly, patterns of in vivo structural brain connectivity can be investigated by using diffusion tensor imaging (DTI), which is a MRI modality based on principle of water diffusion measures (Beaulieu, 2002). One of the most frequently examined aspects of structural connectivity is fractional anisotropy (FA), where higher FA is interpreted as greater structural connectivity between regions.
In recent past, increasing number of studies have used graph theory to describe large scale topological organization of various structural and functional brain networks such as small world property, network efficiency, modular structure, and rich club architecture (Achard et al., 2006; He et al., 2007; Bassett and Bullmore, 2009; Stam, 2010; van den Heuvel and Sporns, 2011, 2013). Using a graph theoretical approach Wierenga and colleagues measured whole brain connectivity in childhood and adolescence and observed a sequential maturational model where they observed that connections between unimodal regions strengthen in childhood, followed by connection of these unimodal regions to association regions, adolescence was characterized by the strengthening of connections between association regions within the frontal and parietal cortex (Wierenga et al., 2016). In fact, their studies showed that white matter strengthening was not homogenous throughout the brain during childhood and adolescence showing strengthening of short association fibers rather than long association fibers. Interestingly, Cao and colleagues in their study (Cao et al., 2014) found an increase in network efficiency during early adulthood which showed a decreasing trend with increasing age. They also observed the proportions of short-distance fibers to be higher than those of long-distance fibers in older adults did. This was also supported by the findings of Sala-Lincon and colleagues, who observed higher average clustering coefficient as well as higher shortest path length with increasing age (Sala-Llonch et al., 2014). Similarly, age related differences in structural connectivity have also observed wide–spread decreases in fractional anisotropy in older compared to younger adults (Damoiseaux and Greicius, 2009; Burzynska et al., 2010) which is also supported by other studies (Gong et al., 2009b; Otte et al., 2015; Zhao et al., 2015). It can therefore be asserted that, although strong functional connectivity between brain regions can exist in absence of strong structural connectivity (Damoiseaux and Greicius, 2009; Zimmermann et al., 2016), similar patterns of age related differences in structural and functional connectivity metrics greatly suggests probable association between these two measures (Andrews-Hanna et al., 2007; Betzel et al., 2014; Fjell et al., 2016; Zimmermann et al., 2016).
Despite of many advances in neuroscience, there is still much to learn about the relationship between white matter integrity and nodal topological properties of structural and functional network. So far, Cao et al. (2014) has employed graph theory across life span ranging 7–85 years. However, they have only used resting state functional fMRI data to examine the topological age related effects. Other study by Hirsiger et al. (2016) (n = 165 age range 64–85) have examined the association of structural and functional connectivity of the cingulum bundle to probe the cognitive and motor performance, where they have suggested that only structural connectivity but not resting state functional connectivity was significantly associated with age. Since, the characterization on structural and functional connectivity is essential to determine the relationship between them. We, therefore have made an effort to ascertain the nature of functional connections by examining the strength of structural connections and vice-versa in a large sample of healthy subjects (n = 458) with an age range of 8–81 years. We hypothesize that changes in both white matter integrity and functional connectivity would be observed in certain networks, and these changes in nodal topological properties of the functional network would be associated with the changes in FA.
Materials and Methods
Participants
DTI and fMRI data were acquired on 3T MRI systems from 458 healthy Chinese Han subjects (right handed; 229 males, 229 females; age range, 8–81 years) at West China Hospital of Sichuan University. The subjects were recruited via local advertising. Among these, 443 participants had fMRI data, 346 participants had DTI data, and 331 participants had both fMRI and DTI data (Table 1). Participants had no history of brain injury, neurological, or psychiatric diseases. Informed written consent was obtained from all participants including parent/guardian consent for younger participants. Research protocol was approved by the local ethics committee.
Table 1. Participant characteristics.
Data Acquisition
All subjects were scanned in a single session without changing their position and were instructed to remain motionless as possible and relax their minds with their eyes open. All participants confirmed they remained awake and alert throughout the scanning session.
Functional MRI
The subjects underwent a resting-state fMRI scan in one of two 3T MRI systems [General Electric (EXCITE, Millwaukee, USA) or Siemens (Trio a Tim, Erlangen, Germany)]. Parameters of General Electric machine were; repetition time = 2,000 ms, echo time = 30 ms, field of view (FOV) = 24 × 24 cm2, flip angle = 90°, slice thickness = 5.0 mm (no gap), voxel size, 3.75 × 3.75 × 5 mm3; matrix, 64 × 64. Parameters of Siemens machine were: repetition time = 2,000 ms, echo time = 30 ms, field of view (FOV) = 24 × 24 cm2, flip angle = 90°, slice thickness = 5.0 mm (no gap), voxel size 3.75 × 3.75 × 5 mm3, matrix 64 × 64.
Diffusion Tensor Imaging
DTI data were also acquired from one of two 3T MRI systems. Parameters of General Electric machine were: TR = 10,000 ms, TE = 70.8 ms, field of view (FOV) = 24 × 24 cm2, resolution matrix 128 × 128, flip angle = 90°, slice thickness = 3.0 mm, A total of 672 slices were acquired for b values of b = 0 and b = 1,000 mm2/s, which were obtained by applying gradients along 15 non-collinear directions. Parameters of Siemens machine were; TR = 6,800 ms, TE = 93 ms, field of view (FOV) = 24 × 24 cm2, resolution matrix 128 × 128, flip angle = 90°, slice thickness = 3.0 mm. 42 slices were acquired for b = 0 and b = 1,000 mm2/s; these were obtained by applying gradients along 30 non-collinear directions.
Data Preprocessing
Preprocessing contained both functional and structural preprocessing steps.
Functional preprocessing was carried out using Statistical Parametric Mapping (SPM 8) (SPM8, http://www.fil.ion.ucl.ac.uk/spm) and GRETNA (He et al., 2008). Briefly, preprocessing was done by (i) discarding first 10 functional volumes for signal equilibration, (ii) slice timing correction for timing offsets, (iii) head motion correction by 3D geometrical displacement, and (iv) normalization to Montreal Neurological Institute (MNI) space. All data used in this study satisfied the criteria of spatial movement in any direction <1.5 mm or degree. Subjects demonstrated no significant group differences in head-motion parameters. Furthermore, linear detrend and band-pass filtering (0.01–0.08 Hz) was performed to reduce the effects of low-frequency drift and high-frequency noise. Subsequently, several nuisance signals including head motion, global mean, and signals from the cerebrospinal fluid and white matter were regressed from the data. For structural preprocessing, raw DTI images were preprocessed using the FSL (FMRIB Software Library, FMRIB, Oxford, UK) (Smith et al., 2004) software package. For each DTI dataset, all diffusion weighted images were affinely coregistered to the b0 image using FLIRT (FMRIB's Linear Image Registration Tool) (Jenkinson and Smith, 2001) with 12 degrees of freedom to correct for eddy current-induced distortion and subtle head motion. Brain mask was created from the b0 image using the BET (Brain extraction Tool) (Smith, 2002) with a fractional intensity threshold of 0.2; FDT (FMRIB's Diffusion Toolbox) (Behrens et al., 2003) was used to fit the tensor model.
For fiber tract identification, a MATLAB-based open source software termed “automatic fiber quantification (AFQ),” which implements both algorithms proposed by Hua and Zhang (Hua et al., 2008; Zhang et al., 2008) was used. Identification procedure included three primary steps; firstly, whole-brain fiber trajectory was performed on the preprocessed tensor images. Secondly, fiber tract segmentation was performed based on the waypoint region of interest (ROI) procedure as described by Wakana et al. (2007). The waypoint ROIs set developed in Mori's lab was warped into individual space from the MNI template space via non-linear transformation. Each fiber was defined as a candidate to a particular fiber group if it crossed through two-way point ROIs that were used to define a specific fiber tract. Thirdly, fiber refinement was accomplished by comparing each candidate fiber to the fiber tract probability map proposed by Hua et al. (2008). Fiber tract probability maps were also transformed into an individual's native space; candidate fibers for a particular fiber groups were assigned scores according to the probability values of the voxels that they passed through. Candidate fibers with low scores were discarded. Additionally, an iterative procedure was used to filter fibers that were aberrantly longer than the mean fiber length or that were distant from the core of the fiber tract. Fourteen fiber bundles were identified according to the predefined ROIs and probability maps. These bundles were the bilateral thalamic radiation, corticospinal tract, inferior fronto-occipital fasciculus, inferior longitudinal fasciculus, superior longitudinal fasciculus as well as the uncinate, arcuate, genu and splenium of the corpus callosum (Figure 1).
Figure 1. The 14 main fiber bundles.
Construction of Functional Brain Network
Topological properties of corresponding brain networks were examined using the GRETNA toolbox. Small world behavior was assessed by coefficient σ (Humphries et al., 2006) which uses a ratio of network clustering and path length to contrast with the same metrics from an equivalent random network. For σ > 1, a network is considered to have small world behavior (Yang et al., 2015). Brain areas (except pons and cerebellum) were parceled into 90 brain regions according to the Automated Anatomical Labeling atlas (Tzourio-Mazoyer et al., 2002) and were used to define the nodes of the functional network. Then, global network properties including small-world properties, efficiency metrics, and regional nodal properties were simultaneously calculated. Briefly, small-world network is an anatomical network containing specialized (segregated) modules and a robust number of intermodular (integrating) links, that has clustering coefficient (Cp), and characteristics path length (Lp) as its properties. Some of the other metrics frequently used are efficiency metrics such as global efficiency (which represents the inverse of the average shortest path length between nodes in the entire network) and local efficiency (which represents the inverse of the average shortest path length between all nearest neighbors of a node). Regional nodal properties such as nodal betweenness (an indicator of a nodes centrality in a network, which is equal to the number of shortest paths, form all vertices to all others that pass through a node). Nodal degree (the degree of an individual node, which is equal to the number of links connected to that node), and nodal efficiency (the average smallest path weight between a given node and all other nodes in the network) (Bullmore and Sporns, 2009). Equations used to calculate these metrics could be found elsewhere (Rubinov and Sporns, 2010).
Statistical Analysis
All statistical analyses were performed using SPSS (version 19). Two multiple regressions were performed to investigate the effects of age on nodal topological properties of the functional brain network and white matter integrity with nodal topological properties and mean FA values as dependent variables while; age, gender, and education kept as independent variables. Additionally, multiple regression model was also used to explore the association between nodal topological properties and white matter integrity by using the value of each nodal topological property as an independent variable and mean FA value of each fiber bundle, maintaining age, gender, and education as dependent variables. Finally, two subgroup analyses were performed to investigate the interactive effect between gender and nodal topological properties or white matter integrity.
Linear Regression of Topological Properties in the Functional Brain Network According to Age
We performed linear regression to analyze the effect of age on the topological properties (global network properties and regional network properties as dependent variables with statistical threshold set as P < 0.05).
Linear Regression of FA Value and Age
At first, 14 fiber bundles from the DTI data and 100 nodal FA values of each fasciculus (Figure 1), then the mean FA values for each of the 14 fiber bundles were calculated. Finally, linear regression was performed to explore the effect of age on white matter integrity with mean FA of each fiber bundle as the dependent variable.
Linear Regression of Nodal Topological Properties According to Mean FA Value
Linear regression model was also used to explore the association between nodal topological properties and white matter integrity. While, the mean FA value of each fiber bundle was used as the independent variable, the value of each nodal topological property was used as the dependent variable.
Results
Linear Regression of the Topological Properties of the Functional Brain Network According to Age
Linear age-related changes in the node betweenness were revealed in 21 brain regions. Positive correlation between age and node degree and node efficiency was observed in bilateral middle and inferior frontal gyrus, bilateral anterior and median cingulated cortex, bilateral postcentral gyrus, superior parietal gyrus, bilateral superior temporal gyrus, middle temporal gyrus, bilateral caudate nucleus and left lenticular nucleus (Table 2, Figures 2B,C, red regions). In addition, negative correlation between age and node degree and node efficiency was observed in bilateral inferior temporal gyrus and medial superior frontal gyrus (Figures 2B,C, blue regions). Likewise, age showed positive correlation with node betweenness in olfactory cortex, fusiform gyrus, superior and middle temporal gyrus (Figure 2A, red regions). However, negative correlation was observed in bilateral inferior and frontal gyrus; left superior frontal gyrus, bilateral insula, right superior and middle occipital gyrus (Figure 2A, blue regions). Figure 3 shows the linear regression results of the effects of age on the small-world properties and efficiency metrics wherein, linear positive age-related changes (p < 0.05) were calculated for global efficiency (β = 0.009, p = 0.004), local efficiency (β = 0.01, p = 0.002), and the clustering coefficient (β = 0.011, p = 0.002). A linear age-related reduction in path length was also observed (β = −0.074, p = 0.004).
Table 2. Regression models assessing the effect of age on nodal topological properties.
Figure 2. Effect of age on regional nodal properties. Significant positive and negative correlations between age and regional nodal parameters are shown in red and blue, respectively. The node sizes indicate the values of the regional nodal parameters. The distributions of nodes showing altered regional properties were visualized with the Brain Net Viewer (http://www.nitrc.org/projects/bnv/). (A) Node betweenness, (B) Node degree, and (C) Node efficiency.
Figure 3. Effect of age on global network properties. Age is positively correlated with normalized global efficiency (Eglob) (β = 0.009, p = 0.004), local efficiency (Eloc) (β = 0.01, p = 0.002), and clustering coefficient (Cp) (β = 0.011, p = 0.002). Age is negatively correlated with path length (Lp) (β = −0.074, p = 0.004).
Linear Regression of FA Value and Age
Linear age-related decreases in the FA value were observed in callosum forceps minor (β = −0.001, p < 0.001) (Figure 4). However, the mean FA values of other fiber bundles did not demonstrate linear relationship with age.
Figure 4. Change with age in the callosum forceps minor fasciculus fractional anisotropy (FA) values. Age is negatively correlated with the mean FA values for the callosum forceps minor (β = −0.001, p < 0.001).
Linear Regression of Nodal Topological Properties and Fasciculus Structure Connectivity
As observed in Supplementary Table 1, all of the mean FA values had either positive or negative association with node betweenness of 49 brain regions. Among these, only 14 showed direct positive anatomical relationship with fiber bundles, which were mainly located in the frontal and occipital lobes. Likewise, the mean FA values of 13 of the fiber bundles had either positive or negative linear relationship with node degree of 52 brain regions, while 14 of the fiber bundles had either positive or negative linear relationship with node efficiency in 60 brain regions. Among these regions, the node degree of 14 and the node efficiency of 16 had a direct anatomical relationship with the fiber bundles (Supplementary Table 2); these relationships were primarily negative, and located in the frontal and occipital lobes. We also observed that the long fiber bundles had significant relationships with more brain regions than the short fiber bundles.
Discussion
In this study, we investigated age related changes in topological properties of brain functional network and structural connectivity. We observed that the nodes and topological properties that changed with increasing age were limited to specific regions, and FA values of most fiber bundles did not alter with increasing age, suggesting that age related changes are more pronounced in functional rather than structural connectivity measures indicating these measures do not have a direct one-to-one mapping. It was also noted that fiber bundles with longer fibers exhibited more pronounced effect on the properties of functional network rather than structural network.
For regional characteristics of functional properties, we observed that age related decreases in node betweenness were primarily located in the superior frontal lobe, right occipital lobe and global hubs. Among 21 brain regions that had significant relationship with age, 11 global hubs identified crucial to efficient communication (Achard et al., 2006; Chen et al., 2008; Iturria-Medina et al., 2008; Gong et al., 2009a; He et al., 2009; Wu et al., 2012) were mostly association cortices (13 out of 44) with no identifiable primary regions. This result supports the hypothesis that age-related changes are characteristics of association cortices as opposed to primary cortices (Albert, 1994). Decreases in local betweenness were primarily located in the superior frontal lobe and right occipital lobe, which is consistent with previous studies (Wu et al., 2012, 2013). Association regions contribute to the integrity of multiple functional systems, such, as memory and attention systems, and are mainly involved intelligent processing and maintenance of superior spiritual activity (Mesulam, 1998). Vulnerability of frontal regions with advancing age might explain the reason behind the cognitive function decline in many elderly populations (Jernigan et al., 2001). Moreover, identification of association cortex supports the hypothesis that age-related changes are the characteristics of association cortex but primary cortices (Park et al., 2008), which also in consistent with “last –in-first-out” hypothesis indicating that late-maturing regions (such as heteromodal association cortices) have damaging effects of aging (Kalpouzos et al., 2009; Terribilli et al., 2011). Age related increases in node degree and node efficiency were predominantly located in the posterior frontal lobe and parietal lobe. The alteration in nodal degree and nodal efficiency were very similar to each other. Nearly half of the 90 brain regions were observed to have changes with regard to these two nodal properties. Among these regions, alteration ratio of sub cortical (8 of 10), limbic (2 of 4), primary (4 of 8), and paralimbic (10 of 24) were similar, but alteration ratio of the association regions (16 of 44) in nodal degree was slightly larger than that of nodal efficiency (15 of 44), suggesting a tight relationship between these two regional nodal properties.
On the contrary to the functional studies, we found that the relationship between nodal topological properties and fasciculus structures were primarily positive for nodal betweenness, and negative for nodal degree and nodal efficiency. These brain regions also had both direct and indirect anatomical relationships with the 14 fiber bundles. Since, the normalized betweenness measures the ability of a node relative to information flow between other nodes within the network, a positive relationship between the fasciculus FA and node betweenness in most regions suggests that increased structural connectivity of the fasciculus may improve communication between a node and other nodes in the network. Likewise, both direct and indirect relationships with the 14 fiber bundles might be due to the co-activation of the regions even when there is no direct structural connection between them (Rubinov and Sporns, 2010), which is also shown by previous studies indicating strong relationship between structural integrity and functional connectivity in resting state networks direct one-to-one relationship between structural and functional connectivity (Bullmore and Sporns, 2009; Damoiseaux and Greicius, 2009).
Other finding of our study was that the FA values of most of the fiber bundles did not alter with increasing age except that of callosum forceps minor. This age-dependent decrease in FA in the callosum forceps minor was in accordance with findings from several recent studies on the age-related alterations in the white matter microstructure (Park et al., 2004; Salat et al., 2005; Ardekani et al., 2007). This alteration in the organization of the corpus callosum during the aging process may explain the reason behind aging to be vulnerable to neurodegenerative disorders such as in Alzheimer's disease, as there might be an interruption of information from the sensory neocortex to the prefrontal neocortex. However, a study by Voineskos et al. (2012), healthy individuals (n = 48) age ranging from 18 to 85 indicated an anterioposterior gradient of age related decline in corpus callosum fibers, where a potential role of regional white matter damage (i.e., posterior fibers of the corpus callosum) in influencing different cognitive performances in healthy subjects was noted. One of the reasons for different findings might be due to difference in sample size, where a larger sample size might have permitted retention of more paths in their study. In addition, larger number of sample subjects between the age group 20 and 40 in our study might also play its role in different findings between these two studies. Other reason includes the methodological differences in these studies where Voineskos and colleagues did not study all white matter tracts (focusing primarily on cortico-cortical white matter tracts). The current study also indicates that the FA value of the most of the fiber bundles did not substantially change with age. However, a direct inference on anatomical connectivity differences cannot be based on this observation alone. The unchanged FA value might also be due to a balance created among alterations in fiber size, density, and myelination or fiber coherence. However, these differences in findings also suggest that white matter alterations are variable throughout the brain.
Apart from the findings above, functional brain network exhibited economical small-world properties in all healthy individuals. In this study, we demonstrated that there was increase in global efficiency (β = 0.009, p = 0.004), local efficiency (β = 0.01, p = 0.002), and clustering coefficient (β = 0.011, p = 0.002) (Figure 2) in functional brain network with the advancing age. An economical small-world offers a topological substrate for specialized or modular processing in local neighborhoods and distributed or integrated processing over the entire network with the combination of both high clustering and low characteristic path length (Sporns and Zwi, 2004; Stam, 2004; Achard et al., 2006; Achard and Bullmore, 2007). Small-world properties have been demonstrated by the studies using fMRI in the human brain functional networks (Eguiluz et al., 2005; Salvador et al., 2005; Achard et al., 2006; He et al., 2007)and have also demonstrated that these properties are the characteristics of large-scale anatomical networks of the human cerebral cortex. Thus, our findings also support the notion that there is the presence of an efficient network structure across the development process.
Limitations
Several issues must be addressed while considering the results of our study. First, the population distribution of our study was not uniform. Proportions of subjects between 20 and 40 year was larger than that of other age groups. Second, we used a DTI-based streamline tractography approach (Mori et al., 1999; Basser et al., 2000) to define the edges of the structural network. Although, this is the most widely applied tractography method primarily due to its simplicity, robustness, and speed (Cheng et al., 2012; Griffa et al., 2013), such tractography method, should be used cautiously in order to resolve crossing fiber bundles (Tournier et al., 2011; Jones et al., 2013). Third, due to limited number of subjects with neuropsychological assessment in our study, we were restricted in determining the relationship between white matter tracts and cognitive ratings in our findings. Further studies will benefit from the use of assessment scales to quantify the association between white matter integrity and cognitive performance. Finally, we adopted Automated Anatomical Labeling (AAL) template as a parcellation scheme. The AAL template is based on sulcal patterns from only one subject. Primary advantage of using AAL template for nodal parcellation is that, it can support a direct comparison of results to previous connectome studies using the same AAL template in healthy adults (Gong et al., 2009a) and patient populations (Zalesky et al., 2011). It is important to consider that although the atlas that we used was carefully checked for registration errors, a probabilistic atlas of the human brain might be better for obtaining regional parcellation or defining individual brain regions through a combination of DTI with fMRI due to the inter-individual variability of anatomical structures (Sporns et al., 2005).
Conclusion
Current study indicates that, age related changes are more pronounced in the functional than in structural measure indicating these measures do not have a direct one-to-one mapping. Frontal regions are more vulnerable with advancing age and fiber bundles with longer fibers exhibited pronounced effect on the functional network properties. Future longitudinal studies would be useful to investigate the changes in the functional and anatomical neural networks that occur with normal aging.
Author Contributions
Study concept and design: SL and QG. Acquisition, analysis, or interpretation of data: HS and JiaL. Drafting of the manuscript: CS and JiaL. Critical revision of the manuscript for important intellectual content: SL and PL. Statistical analysis: HS, JieL, and WZ. Obtained funding: SL and QG. Adminstrative, technical, or material support: CS, YX, YZ, and LY. Study supervision: SL.
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 study was supported by National Natural Science Foundation of China (grant numbers 81371527, 81671664, 81621003, 81220108013, and 81030027) and Program for Changjiang Scholars and Innovative Research Team in University of China (PCSIRT, grant number IRT16R52). QG would like to acknowledge the support from Changjiang Scholar Professorship Awards of China (award number T2014190). SL would like to acknowledge the support from Chang Jiang Scholars of China (award number Q2015154) and the National Program for Support of Top-notch Young Professionals (award number W02070140). The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2018.00318/full#supplementary-material
References
Achard, S., and Bullmore, E. (2007). Efficiency and cost of economical brain functional networks. PLoS Comput. Biol. 3:e17. doi: 10.1371/journal.pcbi.0030017
Achard, S., Salvador, R., Whitcher, B., Suckling, J., and Bullmore, E. (2006). A resilient, low-frequency, small-world human brain functional network with highly connected association cortical hubs. J. Neurosci. 26, 63–72. doi: 10.1523/JNEUROSCI.3874-05.2006
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–935. doi: 10.1016/j.neuron.2007.10.038
Ardekani, S., Kumar, A., Bartzokis, G., and Sinha, U. (2007). Exploratory voxel-based analysis of diffusion indices and hemispheric asymmetry in normal aging. Magn. Reson. Imaging 25, 154–167. doi: 10.1016/j.mri.2006.09.045
Ball, G., Aljabar, P., Zebari, S., Tusor, N., Arichi, T., Merchant, N., et al. (2014). Rich-club organization of the newborn human brain. Proc. Natl. Acad. Sci. U.S.A. 111, 7456–7461. doi: 10.1073/pnas.1324118111
Basser, P. J., Pajevic, S., Pierpaoli, C., Duda, J., and Aldroubi, A. (2000). In vivo fiber tractography using DT-MRI data. Magn. Reson. Med. 44, 625–632. doi: 10.1002/1522-2594(200010)44:4<625::AID-MRM17>3.0.CO;2-O
Bassett, D. S., and Bullmore, E. T. (2009). Human brain networks in health and disease. Curr. Opin. Neurol. 22, 340–347. doi: 10.1097/WCO.0b013e32832d93dd
Baumann, N., and Pham-Dinh, D. (2001). Biology of oligodendrocyte and myelin in the mammalian central nervous system. Physiol. Rev. 81, 871–927. doi: 10.1152/physrev.2001.81.2.871
Beaulieu, C. (2002). The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed. 15, 435–455. doi: 10.1002/nbm.782
Behrens, T. E. J., Woolrich, M. W., Jenkinson, M., Johansen-Berg, H., Nunes, R. G., Clare, S., et al. (2003). Characterization and propagation of uncertainty in diffusion-weighted MR imaging. Magn. Reson. Med. 50, 1077–1088. doi: 10.1002/mrm.10609
Berg, L., Miller, J. P., Baty, J., Rubin, E. H., Morris, J. C., and Figiel, G. (1992). Mild senile dementia of the Alzheimer type. 4. Evaluation of intervention. Ann. Neurol. 31, 242–249. doi: 10.1002/ana.410310303
Betzel, R. F., Byrge, L., He, Y., Goni, J., Zuo, X. N., and Sporns, O. (2014). Changes in structural and functional connectivity among resting-state networks across the human lifespan. Neuroimage 102 (Pt 2), 345–357. doi: 10.1016/j.neuroimage.2014.07.067
Bullmore, E. T., and Bassett, D. S. (2011). Brain graphs: graphical models of the human brain connectome. Annu. Rev. Clin. Psychol. 7, 113–140. doi: 10.1146/annurev-clinpsy-040510-143934
Bullmore, E., and Sporns, O. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci. 10, 186–198. doi: 10.1038/nrn2575
Burzynska, A. Z., Preuschhof, C., Backman, L., Nyberg, L., Li, S. C., Lindenberger, U., et al. (2010). Age-related differences in white matter microstructure: region-specific patterns of diffusivity. Neuroimage 49, 2104–2112. doi: 10.1016/j.neuroimage.2009.09.041
Cao, M., Wang, J. H., Dai, Z. J., Cao, X. Y., Jiang, L. L., Fan, F. M., et al. (2014). Topological organization of the human brain functional connectome across the lifespan. Dev. Cogn. Neurosci. 7, 76–93. doi: 10.1016/j.dcn.2013.11.004
Chen, Z. J., He, Y., Rosa-Neto, P., Germann, J., and Evans, A. C. (2008). Revealing modular architecture of human brain structural networks by using cortical thickness from MRI. Cereb. Cortex18, 2374–2381. doi: 10.1093/cercor/bhn003
Cheng, H., Wang, Y., Sheng, J., Kronenberger, W. G., Mathews, V. P., Hummer, T. A., et al. (2012). Characteristics and variability of structural networks derived from diffusion tensor imaging. Neuroimage 61, 1153–1164. doi: 10.1016/j.neuroimage.2012.03.036
Damoiseaux, J. S., and Greicius, M. D. (2009). Greater than the sum of its parts: a review of studies combining structural connectivity and resting-state functional connectivity. Brain Struct. Funct. 213:525–533. doi: 10.1007/s00429-009-0208-6
Dennis, E. L., and Thompson, P. M. (2014). Functional brain connectivity using fMRI in aging and Alzheimer's disease. Neuropsychol. Rev. 24, 49–62. doi: 10.1007/s11065-014-9249-6
Dubois, J., Dehaene-Lambertz, G., Kulikova, S., Poupon, C., Huppi, P. S., and Hertz-Pannier, L. (2014). The early development of brain white matter: a review of imaging studies in fetuses, newborns and infants. Neuroscience 276, 48–71. doi: 10.1016/j.neuroscience.2013.12.044
Eguiluz, V. M., Chialvo, D. R., Cecchi, G. A., Baliki, M., and Apkarian, A. V. (2005). Scale-free brain functional networks. Phys. Rev. Lett. 94:018102. doi: 10.1103/PhysRevLett.94.018102
Ferreira, L. K., and Busatto, G. F. (2013). Resting-state functional connectivity in normal brain aging. Neurosci. Biobehav. Rev. 37, 384–400. doi: 10.1016/j.neubiorev.2013.01.017
Fjell, A. M., Sneve, M. H., Storsve, A. B., Grydeland, H., Yendiki, A., and Walhovd, K. B. (2016). Brain events underlying episodic memory changes in aging: a longitudinal investigation of structural and functional connectivity. Cereb. Cortex 26, 1272–1286. doi: 10.1093/cercor/bhv102
Gong, G., He, Y., Concha, L., Lebel, C., Gross, D. W., Evans, A. C., et al. (2009a). Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diffusion tensor imaging tractography. Cereb. Cortex 19, 524–536. doi: 10.1093/cercor/bhn102
Gong, G., Rosa-Neto, P., Carbonell, F., Chen, Z. J., He, Y., and Evans, A. C. (2009b). Age- and gender-related differences in the cortical anatomical network. J. Neurosci. 29, 15684–15693. doi: 10.1523/JNEUROSCI.2308-09.2009
Grady, C. L., Haxby, J. V., Horwitz, B., Sundaram, M., Berg, G., Schapiro, M., et al. (1988). Longitudinal study of the early neuropsychological and cerebral metabolic changes in dementia of the Alzheimer type. J. Clin. Exp. Neuropsychol. 10, 576–596. doi: 10.1080/01688638808402796
Griffa, A., Baumann, P. S., Thiran, J. P., and Hagmann, P. (2013). Structural connectomics in brain diseases. Neuroimage 80, 515–526. doi: 10.1016/j.neuroimage.2013.04.056
Hagmann, P., Sporns, O., Madan, N., Cammoun, L., Pienaar, R., Wedeen, V. J., et al. (2010). White matter maturation reshapes structural connectivity in the late developing human brain. Proc. Natl. Acad. Sci. U.S.A. 107, 19067–19072. doi: 10.1073/pnas.1009073107
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. doi: 10.1093/cercor/bhl149
He, Y., Chen, Z., and Evans, A. (2008). Structural insights into aberrant topological patterns of large-scale cortical networks in Alzheimer's disease. J. Neurosci. 28, 4756–4766. doi: 10.1523/JNEUROSCI.0141-08.2008
He, Y., Wang, J., Wang, L., Chen, Z. J., Yan, C., Yang, H., et al. (2009). Uncovering intrinsic modular organization of spontaneous brain activity in humans. PLoS ONE 4:e5226. doi: 10.1371/journal.pone.0005226
Hirsiger, S., Koppelmans, V., Merillat, S., Liem, F., Erdeniz, B., Seidler, R. D., et al. (2016). Structural and functional connectivity in healthy aging: associations for cognition and motor behavior. Hum. Brain Mapp. 37, 855–867. doi: 10.1002/hbm.23067
Hua, K., Zhang, J., Wakana, S., Jiang, H., Li, X., Reich, D. S., et al. (2008). Tract probability maps in stereotaxic spaces: analyses of white matter anatomy and tract-specific quantification. Neuroimage 39, 336–347. doi: 10.1016/j.neuroimage.2007.07.053
Humphries, M. D., Gurney, K., and Prescott, T. J. (2006). The brainstem reticular formation is a small-world, not scale-free, network. Proc. Biol Sci. 273, 503–511. doi: 10.1098/rspb.2005.3354
Iturria-Medina, Y., Sotero, R. C., Canales-Rodriguez, E. J., Aleman-Gomez, Y., and Melie-Garcia, L. (2008). Studying the human brain anatomical network via diffusion-weighted MRI and Graph Theory. Neuroimage 40, 1064–1076. doi: 10.1016/j.neuroimage.2007.10.060
Jenkinson, M., and Smith, S. (2001). A global optimisation method for robust affine registration of brain images. Med. Image Anal. 5, 143–156. doi: 10.1016/S1361-8415(01)00036-6
Jernigan, T. L., Archibald, S. L., Fennema-Notestine, C., Gamst, A. C., Stout, J. C., Bonner, J., et al. (2001). Effects of age on tissues and regions of the cerebrum and cerebellum. Neurobiol. Aging 22, 581–594. doi: 10.1016/S0197-4580(01)00217-2
Jones, D. K., Knosche, T. R., and Turner, R. (2013). White matter integrity, fiber count, and other fallacies: the do's and don'ts of diffusion MRI. Neuroimage 73, 239–254. doi: 10.1016/j.neuroimage.2012.06.081
Kalpouzos, G., Chetelat, G., Baron, J. C., Landeau, B., Mevel, K., Godeau, C., et al. (2009). Voxel-based mapping of brain gray matter volume and glucose metabolism profiles in normal aging. Neurobiol. Aging 30, 112–124. doi: 10.1016/j.neurobiolaging.2007.05.019
Lebel, C., Walker, L., Leemans, A., Phillips, L., and Beaulieu, C. (2008). Microstructural maturation of the human brain from childhood to adulthood. Neuroimage 40, 1044–1055. doi: 10.1016/j.neuroimage.2007.12.053
Mesulam, M. M. (1998). From sensation to cognition. Brain 121 (Pt 6), 1013–1052. doi: 10.1093/brain/121.6.1013
Mori, S., Crain, B. J., Chacko, V. P., and van Zijl, P. C. (1999). Three-dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Ann. Neurol. 45, 265–269. doi: 10.1002/1531-8249(199902)45:2<265::AID-ANA21>3.0.CO;2-3
Morris, J. C., Edland, S., Clark, C., Galasko, D., Koss, E., Mohs, R., et al. (1993). The consortium to establish a registry for Alzheimer's disease (CERAD). Part IV. Rates of cognitive change in the longitudinal assessment of probable Alzheimer's disease. Neurology 43, 2457–2465. doi: 10.1212/WNL.43.12.2457
Morrison, J. H., and Hof, P. R. (1997). Life and death of neurons in the aging brain. Science 278, 412–419.
Otte, W. M., van Diessen, E., Paul, S., Ramaswamy, R., Subramanyam Rallabandi, V. P., Stam, C. J., et al. (2015). Aging alterations in whole-brain networks during adulthood mapped with the minimum spanning tree indices: the interplay of density, connectivity cost and life-time trajectory. Neuroimage 109, 171–189. doi: 10.1016/j.neuroimage.2015.01.011
Park, H. J., Kim, J. J., Lee, S. K., Seok, J. H., Chun, J., Kim, D. I., et al. (2008). Corpus callosal connection mapping using cortical gray matter parcellation and DT-MRI. Hum. Brain Mapp. 29, 503–516. doi: 10.1002/hbm.20314
Park, H. J., Westin, C. F., Kubicki, M., Maier, S. E., Niznikiewicz, M., Baer, A., et al. (2004). White matter hemisphere asymmetries in healthy subjects and in schizophrenia: a diffusion tensor MRI study. Neuroimage 23, 213–223. doi: 10.1016/j.neuroimage.2004.04.036
Rubinov, M., and Sporns, O. (2010). Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52, 1059–1069. doi: 10.1016/j.neuroimage.2009.10.003
Sala-Llonch, R., Junque, C., Arenaza-Urquijo, E. M., Vidal-Pineiro, D., Valls-Pedret, C., Palacios, E. M., et al. (2014). Changes in whole-brain functional networks and memory performance in aging. Neurobiol. Aging 35, 2193–2202. doi: 10.1016/j.neurobiolaging.2014.04.007
Salat, D. H., Tuch, D. S., Greve, D. N., van der Kouwe, A. J. W., 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. doi: 10.1016/j.neurobiolaging.2004.09.017
Salvador, R., Suckling, J., Coleman, M. R., Pickard, J. D., Menon, D., and Bullmore, E. (2005). Neurophysiological architecture of functional magnetic resonance images of human brain. Cereb. Cortex 15, 1332–1342. doi: 10.1093/cercor/bhi016
Smith, S. M. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E. J., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23, S208–S219. doi: 10.1016/j.neuroimage.2004.07.051
Sporns, O. (2012). From simple graphs to the connectome: networks in neuroimaging. Neuroimage 62, 881–886. doi: 10.1016/j.neuroimage.2011.08.085
Sporns, O., and Zwi, J. D. (2004). The small world of the cerebral cortex. Neuroinformatics 2, 145–162. doi: 10.1385/NI:2:2:145
Sporns, O., Tononi, G., and Kotter, R. (2005). The human connectome: a structural description of the human brain. PLoS Comput. Biol. 1:e42. doi: 10.1371/journal.pcbi.0010042
Stam, C. J. (2004). Functional connectivity patterns of human magnetoencephalographic recordings: a ‘small-world’ network? Neurosci. Lett. 355, 25–28. doi: 10.1016/j.neulet.2003.10.063
Stam, C. J. (2010). Characterization of anatomical and functional connectivity in the brain: a complex networks perspective. Int. J. Psychophysiol. 77, 186–194. doi: 10.1016/j.ijpsycho.2010.06.024
Tamnes, C. K., Ostby, Y., Fjell, A. M., Westlye, L. T., Due-Tonnessen, P., and Walhovd, K. B. (2010). Brain maturation in adolescence and young adulthood: regional age-related changes in cortical thickness and white matter volume and microstructure. Cereb. Cortex 20, 534–548. doi: 10.1093/cercor/bhp118
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. doi: 10.1016/j.neurobiolaging.2009.02.008
Tournier, J. D., Mori, S., and Leemans, A. (2011). Diffusion tensor imaging and beyond. Magn. Reson. Med. 65, 1532–1556. doi: 10.1002/mrm.22924
Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage 15, 273–289. doi: 10.1006/nimg.2001.0978
van den Heuvel, M. P., and Sporns, O. (2011). Rich-club organization of the human connectome. J. Neurosci. 31, 15775–15786. doi: 10.1523/JNEUROSCI.3539-11.2011
van den Heuvel, M. P., and Sporns, O. (2013). An anatomical substrate for integration among functional networks in human cortex. J. Neurosci. 33, 14489–14500. doi: 10.1523/JNEUROSCI.2128-13.2013
van den Heuvel, M. P., Kersbergen, K. J., de Reus, M. A., Keunen, K., Kahn, R. S., Groenendaal, F., et al. (2015). The neonatal connectome during preterm brain development. Cereb. Cortex 25, 3000–3013. doi: 10.1093/cercor/bhu095
van der Knaap, M. S., Valk, J., Bakker, C. J., Schooneveld, M., Faber, J. A., Willemse, J., et al. (1991). Myelination as an expression of the functional maturity of the brain. Dev. Med. Child Neurol. 33, 849–857. doi: 10.1111/j.1469-8749.1991.tb14793.x
Voineskos, A. N., Rajji, T. K., Lobaugh, N. J., Miranda, D., Shenton, M. E., Kennedy, J. L., et al. (2012). Age-related decline in white matter tract integrity and cognitive performance: a DTI tractography and structural equation modeling study. Neurobiol. Aging 33, 21–34. doi: 10.1016/j.neurobiolaging.2010.02.009
Wakana, S., Caprihan, A., Panzenboeck, M. M., Fallon, J. H., Perry, M., Gollub, R. L., et al. (2007). Reproducibility of quantitative tractography methods applied to cerebral white matter. Neuroimage 36, 630–644. doi: 10.1016/j.neuroimage.2007.02.049
Wierenga, L. M., van den Heuvel, M. P., van Dijk, S., Rijks, Y., de Reus, M. A., and Durston, S. (2016). The development of brain network architecture. Hum. Brain Mapp. 37, 717–729. doi: 10.1002/hbm.23062
Wu, K., Taki, Y., Sato, K., Hashizume, H., Sassa, Y., Takeuchi, H., et al. (2013). Topological organization of functional brain networks in healthy children: differences in relation to age, sex, and intelligence. PLoS ONE 8:e55347. doi: 10.1371/journal.pone.0055347
Wu, K., Taki, Y., Sato, K., Kinomura, S., Goto, R., Okada, K., et al. (2012). Age-related changes in topological organization of structural brain networks in healthy individuals. Hum. Brain Mapp. 33, 552–568. doi: 10.1002/hbm.21232
Yang, S., Hua, P., Shang, X., Cui, Z., Zhong, S., Gong, G., et al. (2015). Deficiency of brain structural sub-network underlying post-ischaemic stroke apathy. Eur. J. Neurol. 22, 341–347. doi: 10.1111/ene.12575
Zalesky, A., Fornito, A., Seal, M. L., Cocchi, L., Westin, C. F., Bullmore, E. T., et al. (2011). Disrupted axonal fiber connectivity in schizophrenia. Biol. Psychiatry 69, 80–89. doi: 10.1016/j.biopsych.2010.08.022
Zhang, W., Olivi, A., Hertig, S. J., van Zijl, P., and Mori, S. (2008). Automated fiber tracking of human brain white matter using diffusion tensor imaging. Neuroimage 42, 771–777. doi: 10.1016/j.neuroimage.2008.04.241
Zhao, T., Cao, M., Niu, H., and Zuo, X. N. (2015). Age-related changes in the topological organization of the white matter structural connectome across the human lifespan. Hum. Brain Mapp. 36, 3777–3792. doi: 10.1002/hbm.22877
Keywords: healthy participants, magnetic resonance imaging, topological properties, fractional anisotropy, white matter
Citation: Shah C, Liu J, Lv P, Sun H, Xiao Y, Liu J, Zhao Y, Zhang W, Yao L, Gong Q and Lui S (2018) Age Related Changes in Topological Properties of Brain Functional Network and Structural Connectivity. Front. Neurosci. 12:318. doi: 10.3389/fnins.2018.00318
Received: 28 November 2017; Accepted: 24 April 2018;
Published: 15 May 2018.
Edited by:
Giuseppe Pignataro, University of Naples Federico II, ItalyReviewed by:
Raffaella Gozzelino, Centro de Estudos de Doenças Crónicas (CEDOC), PortugalFrancesca Trojsi, Università degli Studi della Campania ‘Luigi Vanvitelli’ Naples, Italy
Copyright © 2018 Shah, Liu, Lv, Sun, Xiao, Liu, Zhao, Zhang, Yao, Gong and Lui. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Su Lui, lusuwcums@tom.com
†These authors have contributed equally to this work.