Mapping Functional Connectome Changes of the Human Brain across the Life Span
M. Cao,Zhengjia Dai,Fengmei Fan,Lili Jiang,X. Cao,Mingrui Xia,N. Shu,XiNian Zuo,Yong He
2011-01-01
Abstract:Introduction: Uncovering the lifespan changes of the human brain is definitely thrilling and a fundamental goal in discovery neuroscience. The lifespan trajectory coding both normal developing and aging stages of the brain has reflected dramatic changes both its structure [1-3] and function [4-6]. It has been commonly accepted that human brain is structurally and functionally organized into a complex network allowing the segregation and integration of information processing [7,8]. In particular, the brain network, describing the interaction between different brain regions, has demonstrated consistent non-trivial properties and reliable changes throughout normal development, aging and in various pathological conditions. Beyond these advances, how the brain network is reorganized through the whole lifespan or its lifespan trajectory has been rarely studied. Here, armed with the graph theory, we aim to chart the lifespan trajectory of human whole-brain functional networks based on resting-state functional MRI (R-fMRI) datasets from 150 subjects collected from the same scanner. Materials and Methods: Participants and Data acquisition: A public R-fMRI dataset including total 150 right-handed healthy people (males, 87; females, 63; age range, 7-85 years) from the NKI/Rockland sample (http://fcon_1000.projects.nitrc.org/indi/pro/nki.html) was employed. All subjects have no history of neurological and psychiatric disorders. The resting state fMRI scans were performed on the Siemens Trio 3.0 Tesla MRI scanner using an echo-planar imaging (EPI) sequence (TR/TE=2500/30 ms, FA=80°, FOV=216 mm, matrix=256×256, slices=38, thickness=3.0 mm, 260 volumes). High-resolutionT1-weighted images were also collected by magnetization-prepared rapid gradient echo (MPRAGE) sequence (TR/TE=2500/3.5 ms, FA=8°, thickness=1.0 mm, slices=192, matrix=256×256, FOV=256mm). Functional data preprocessing: Data preprocessing was similar to as implemented in [6], producing individual residual timeseries in the MNI152 standard space. Network Construction: A group mask in standard space was made which includes all voxels with non-zeros standard deviation of bold time series and gray matter tissue probability greater than 20%. A high-resolution parcellation template was generated by randomly partitioning the group mask into 1024 contiguous and uniform regions using an algorithm described in [9]. Resting state functional connectivity (RSFC) for each pair of 1024 regions of interest was obtained by computing Pearson correlation coefficients of their representative (i.e., mean across all voxels) time series. Each subject’s weighted functional connectivity network was then constructed based on a Bonforroni correction procedure. Network properties analysis: We studied the age-related variation trends of the following regional characteristic and connection property: 1) the regional correlation strength which measures the average correlation extent by which a given region, i, is connected to the rest of the network. Regions with correlation strength greater than the mean and two standard deviations of the global network correlation strength were defined as hubs in every person. 2) The anatomical distance of the connection between region i and region j is defined as the Euclidean distance between centroids of the two regions. The percentages of long-distance (>75mm) and short-distance (≤75mm) connection of the functional connectome were computed for each subject. Results: The probability map of hubs in population was gotten and several specific regions including default-mode network regions, attention network regions, motor-related regions and visual cortex were detected to be functional hubs consistently (Figure 1.A). Besides, we observed significantly linear decreases of RSFC strength with age primarily located in several default-mode (e.g., medial frontal and parietal cortex, lateral temporal and parietal cortex) and attention network regions (e.g., insula, dorsal anterior cingulate cortex, lateral frontal cortex and temporal-parietal conjunction) while few regions including SMA and temporal pole showed significantly linear increases (p<0.05, FDR corrected) (Figure 1.B). Meanwhile, significant quadratic trajectories of RSFC strength on age were also detected within the left hippocampus and bilateral superior temporal gyrus (p<0.01, uncorrected) (Figure 1.C). The network sparsity (i.e., number of graph edges or network connections showing significant RSFC) between brain regions demonstrated a linear trend of decrease with age (not significant) through the lifespan (p=0.10). The percentage of long-distance connections was found increase up to age 42 and then decrease while the short-distance connections showed an inverse trajectory (p=0.01) (Figure 2).