White Matter Structural Connectivity Is Not Correlated to Cortical Resting-State Functional Connectivity over the Healthy Adult Lifespan

Introduction

It is widely accepted that the normal human aging process involves changes in the brain's structural and functional connections. Understanding these changes will greatly improve our ability to diagnose and treat age-related neurodegenerative diseases, such as Alzheimer's Disease (AD), amyotrophic lateral sclerosis (ALS), and Parkinson's Disease (PD) (Pievani et al., 2014; Iturria-Medina and Evans, 2015; Gao and Wu, 2016). Non-invasive neuroimaging techniques including diffusion tensor imaging (DTI) and resting-state functional MRI (rs-fMRI) permit the investigation of white and gray matter connectivity in the brain. Metrics derived from DTI are used to quantify the white matter (WM) microstructure [termed structural connectivity (SC)], and correlations of the blood oxygen level dependent (BOLD) time signals computed from rs-fMRI are used to quantify the strength of resting state functional connections between distinct gray matter (GM) regions [termed functional connectivity (FC)].

Several studies have independently used DTI and rs-fMRI to demonstrate changes in SC and FC over the healthy human lifespan. In general, DTI studies have observed a non-linear inverted U-shaped trajectory association between age and fractional anisotropy (FA), and an U-shaped trajectory (opposite to FA) for axial, mean, and radial diffusivity (AD, MD, RD, respectively; Westlye et al., 2010; Lebel et al., 2012; Chen et al., 2013). Furthermore, previous studies have also used DTI to demonstrate that the degree of age-related cognitive decline correlates with WM microstructural alterations (Madden et al., 2012; Hawkins et al., 2015). On the other hand, rs-fMRI studies have reported both negative and positive (as well as both linear and non-linear) associations between age and FC, which were dependent on the brain region under investigation (Wang et al., 2012; Cao et al., 2014; Fjell et al., 2015a). In addition, rs-fMRI studies have also demonstrated that cognitive decline is related to decreased FC in the salience network (Onoda et al., 2012). As functionally linked cortical regions are connected anatomically via the underlying WM architecture (van den Heuvel et al., 2009), investigating SC and FC simultaneously to determine their interrelationship has the potential to provide a better, more comprehensive, understanding of the brain changes associated with aging.

The relationship between SC and FC, however, is not straightforward. For example, in one case following complete commissurotomy, FC was preserved across hemispheres between regions associated with the default mode network (Uddin et al., 2008). Similarly, another study showed no statistical differences in inter-hemispheric FC between subjects with complete agenesis of the corpus callosum and normal controls of comparable age, gender and IQ (Tyszka et al., 2011). Hence, multiple underlying (and possibly indirect) structural architectures must exist to support functional networks. Indeed, studies of healthy subjects have demonstrated strong FC between cortical regions with direct structural (i.e., WM) connections, as well as between regions in the absence of a direct WM pathway (Koch et al., 2002; Honey et al., 2009). Furthermore, Honey et al. (2009) also demonstrated that FC between indirectly connected regions was mediated by WM tract distance. Hence, these studies provide evidence that functionally connected cortical regions of a resting-state network (RSN) are either connected anatomically via a direct WM pathway or indirectly via WM tracts through one or more intermediate cortical or subcortical regions. However, whether the change in SC and FC measures with age are correlated within multiple large-scale RSNs across the healthy lifespan remains unknown or poorly understood.

The first study to adopt a multi-modal analysis of SC and FC demonstrated that FC between medial prefrontal cortex and the posterior cingulate/retrosplenial cortex (regions associated with the default mode network) was positively correlated with mean FA of the superior longitudinal fasciculus and cingulum WM tracts in elderly subjects (Andrews-Hanna et al., 2007). Subsequent multi-modal studies that compared data between two groups of subjects (i.e., young vs. elderly) demonstrated that FA and MD were significantly correlated with FC (Fjell et al., 2015b; Marstaller et al., 2015). To our knowledge, only three recent studies have examined SC and FC data from healthy subjects across the adult lifespan (Betzel et al., 2014; Lee et al., 2015; Fjell et al., 2016). Using a graph theory analysis approach, Betzel et al. showed that on average, FC remained relatively constant over the adult lifespan for regions with direct structural connections, but the change in FC with age was progressively greater as the structural connection distance between regions increased. Lee et al. demonstrated significant increases in SC and FC with age between prefrontal cortex and posterior regions of the parietal and temporal lobes, suggesting the brain adapts to neural challenges during normal aging. The study by Fjell et al. demonstrated a weak relationship between SC and FC measures for certain major WM tracts and their associated regions of the default mode network. While these studies provided important insights into SC and FC changes associated with normal aging, it remains unclear if WM microstructural changes over the adult lifespan are correlated with the corresponding FC changes within multiple large-scale RSNs. As there are only few studies in literature that investigated the relationship between SC and FC over the adult lifespan, and the results reported are related to certain specific GM regions and WM tracts, a study that investigates more broadly across multiple commonly described RSNs and the associated WM tracts is warranted in normal subjects over a wide age span.

In this study, we hypothesized that WM SC (using measurements of MD, FA, mean tract length, and number of streamlines) derived from DTI data across the adult lifespan are correlated with corresponding FC measures within seven commonly described large-scale RSNs. The aims of this work were to investigate (1) the relationship in SC and FC measures with age, (2) sex differences of these measures with age, and (3) the relationship between SC and FC measures, within seven commonly described RSNs in healthy participants across the adult lifespan.

Materials and Methods

Participants

As part of an on-going normative study (the Calgary Normative Study) that was approved by the University of Calgary Research Ethics Board, healthy community-dwelling participants were recruited, initially screened over the phone, and only those who indicated no known neurological diseases and no contraindications to MR imaging were enrolled in the study. Informed written consent and basic medical history were obtained from each eligible participant prior to imaging. Participants were excluded from the analysis if there were medically significant incidental findings found on their MR images. In addition to MR imaging, the Montreal Cognitive Assessment (MoCA) was administered to each participant as a brief screening tool for mild cognitive impairment or dementia. Two hundred and twenty-one participants provided data for this study. Five subjects were excluded from our analysis due to incidentally discovered, potentially medically significant findings (1 subject), incomplete scan or missing data (2 subjects), or poor quality DTI data (2 subjects). Furthermore, 39 (male/female = 17/22) participants who obtained scores of <26 (out of 30) on the MoCA were excluded, as this falls outside the normal range. Subsequently, 177 subjects (aged 18–87 years; male/female = 69/108; Table 1) were included in the study.

TABLE 1

Table 1. Summary of study subject characteristics.

Image Acquisition

MR imaging was performed on a 3.0 T clinical scanner (Discovery MR750; GE Healthcare, Waukesha, WI) using a 12-channel phased-array head coil. The image acquisition protocol included DTI, rs-fMRI, T2-weighted FLAIR, and T1-weighted imaging sequences. DTI acquisition employed a single-shot spin-echo echo-planar imaging (EPI) sequence [echo time (TE) = 80 ms; repetition time (TR) = 9,000 or 10,000 ms; 48–52 contiguous 3-mm thick slices; field of view (FOV) = 240 × 240 mm; acquired matrix = 80 × 80 interpolated to 256 × 256; reconstructed in-plane resolution = 0.94 × 0.94 mm] with diffusion sensitizing gradients applied in 31 non-collinear directions (b = 1000 s/mm²) and 4 b = 0 s/mm² volumes. rs-fMRI acquisition consisted of a single-shot gradient-echo EPI sequence (TE = 30 ms; TR = 2,000 ms; 37 contiguous 3.8-mm thick slices; acquired matrix = 64 × 64; FOV = 240 × 240 mm; reconstructed voxel size = 3.8 mm isotropic) and acquired 150 whole brain volumes over a 5-min interval. T2-weighted FLAIR images were acquired using an inversion recovery prepared fast spin echo sequence [flip angle = 111°; inversion time (TI) = 2,250 ms; TE = 141.4 ms; TR = 9,000 ms; 48 contiguous 3-mm thick slices, FOV = 240 × 240 mm; reconstructed voxel size = 0.94 × 0.94 mm]. T1-weighted anatomical images were acquired using a 3D inversion recovery prepared spoiled gradient-echo sequence [flip angle = 8°; inversion time (TI) = 650 ms; TE = 2.5 ms; TR = 6.3 ms; acquired matrix size = 256 × 256 × 166; phase FOV = 85%; reconstructed voxel size = 1 mm isotropic].

Image Processing

Identification of Cortical Regions

Images from each subject were processed using an in-house automated pipeline developed from freely available software packages and a semi-automated tool (Cerebra-WML; Gobbi et al., 2012) for WM hyper-intensity mask identification (see Figure 1). A cortical parcellation atlas (Yeo et al., 2011) was used to define cortical regions related to seven resting-state networks (RSNs). The atlas was constructed by processing rs-fMRI data acquired from 1,000 healthy participants and used a clustering algorithm to parcellate the cortex into multiple RSNs. There are other whole-brain atlases available with cortical parcellation including (Auzias et al., 2016; Fan et al., 2016) that can also be used in such multi-modal analysis, but the template by Yeo et al. was chosen as it is widely adopted by many previous studies. The coarse-resolution seven RSNs parcellation was chosen from the selected template atlas over the fine-resolution 17 networks parcellation simply to reduce computing resources and processing times. Nevertheless, the coarse parcellation accurately reflects seven distinct and commonly identified RSNs. Regions within each network were extracted based on four pre-specified anatomical lobes (i.e., frontal, parietal, temporal, occipital) and resulted in the identification of 40 cortical regions across both hemispheres. Specifically, cortical regions in each hemisphere for the seven networks were

(1) visual network: occipital, intraparietal, inferior temporal regions;

(2) somato-motor network: posterior frontal, anterior parietal, superior, and anterior temporal regions;

(3) dorsal attention network: superior frontal, superior parietal, posterior temporal regions;

(4) ventral attention network: inferior medial frontal, inferior lateral and superior medial parietal, superior temporal regions;

(5) limbic network: inferior prefrontal, inferior temporal regions;

(6) frontal-parietal network: lateral frontal, superior posterior parietal, inferior temporal regions; and

(7) default mode network: medial lateral frontal, inferior parietal, lateral temporal regions.

FIGURE 1

Figure 1. In-house developed pipeline used to process MR imaging data of each subject. Data processing and analysis through this pipeline was performed on an iMac (2.9GHz quad-core Intel Core i5; 32GB 1600MHz DDR3 memory) and high performance computing clusters (https://www.westgrid.ca; used to generate cortical surface labels—explained above). SC and FC measures were derived for each resting-state network (RSN). Subsequently, the values of each connectivity measure were averaged in each network.

Cortical surface labels using T1-weighted images were obtained from each subject (FreeSurfer; http://surfer.nmr.mgh.harvard.edu) and were used to transform each region from the atlas space to the subject native space. Subsequently, the 40 cortical regions were transformed to the subject DTI and rs-fMRI spaces for analysis.

Measurement of Structural Connectivity (SC)

DTI data were first corrected for motion and eddy current distortion using FSL (FMRIB Software Library, version 5.0.8; http://www.fmrib.ox.ac.uk/fsl; Jenkinson et al., 2012). Maps of MD and FA were computed from the DTI data (Diffusion Toolkit; http://trackvis.org/dtk/), as well as statistics of mean tract length and number of streamlines were extracted from the tractography algorithm. Whole brain WM tracts were delineated by deterministic tractography using the second order Runge-Kutta algorithm (Basser et al., 2000) with the FA threshold set to 0.20 to exclude gray matter voxels and the angle threshold set to 35° to exclude tracks with sharp curvature. The 40 cortical regions from the template were first dilated and used as seeding/target regions to delineate WM tracts that either originate or terminate at each cortical GM region of the seven RSNs (TrackVis; http://trackvis.org). Example of the cortical seeding regions used for the visual network and the associated WM tracts delineated is shown in Figure 2. WM hyper-intensity voxels were excluded from FA and MD maps using the masks defined from FLAIR images.

FIGURE 2

Figure 2. Delineation of WM tracts using TrackVis in the visual resting-state network (RSN). The red, yellow, and green regions are the occipital, parietal, and temporal regions of the visual RSN, and the voxels in blue represent the WM voxels of all the tracts that are part of the RSN.

Measurement of Functional Connectivity (FC)

The rs-fMRI data were first processed using the FSL package and included skull stripping (Brain Extraction Tool, BET; Smith, 2002), interleaved slice timing correction and motion correction (MCFLIRT algorithm; Jenkinson et al., 2002), spatial smoothing (6-mm full width at half maximum), and temporal high-pass filtering (>0.01 Hz) to eliminate low frequency artifacts. Noise components in the pre-processed data were removed using independent component analysis-based methods (FSL Xnoiseifier, FIX; Salimi-Khorshidi et al., 2014). In addition, the time-points of large motion perturbations in the original four-dimensional resting-state time series dataset were identified (FSL Motion Outliers) using a threshold of 0.2 mm (a stringent threshold for scrubbing; Power et al., 2014) applied for frame-wise displacement. A confound matrix was created for the large motion time-points and was included as additional event variable in the analysis to remove nuisance variables from the resting-state dataset (FMRI Expert Analysis Tool, FEAT). Cerebrospinal fluid (CSF) and WM masks were manually drawn on T1-weighted images from each subject and then transformed into the rs-fMRI image space. Six motion parameters and the average time series from CSF and WM masks were regressed out as nuisance variables (FEAT) from the pre-processed noise reduced (FIX) four-dimensional rs-fMRI dataset. The time-point volumes with motion greater than the threshold were removed from the “cleaned” rs-fMRI dataset. Average rs-fMRI-obtained BOLD signal of all voxels within individual regions from this processed dataset was computed for all time points for subsequent analysis (MATLAB R2015b; MathWorks, Natick, MA). Pearson correlation coefficients (r) of the averaged BOLD signal time series between pairs of regions in each network were converted to z-scores using Fisher's r-to-z transformation {z = 0.5 × ln [(1 + r)/(1 − r)]}.

Statistical Analyses

Average SC and FC measures were computed for each network, and these data were used for all the statistical analyses (SPSS version 22.0; IBM Corp, Armonk, NY). In all analyses, the critical value was chosen as α = 0.05 and multiple comparison corrections were applied across the seven networks using the Bonferroni method (Holm, 1979). Therefore, p < α/7 ≈ 0.007 were considered to be significant.

The following statistical tests (T1 to T4) were performed to address the three objectives of this study:

T1: Pearson correlations of each SC and FC measure with age in each network.

T2: Pearson correlations of each SC and FC measure with age in each network for male and female subjects separately. Following this, the slope of the linear regression line for each connectivity measure was tested for sex differences using t-tests provided in the Real Statistics Resource Pack software (Release 4.3; Zaiontz, 2013–2015).

T3: In addition to testing the age relationship with SC and FC measures above, we tested whether sex is a significant predictor of SC and FC using multiple linear regression. Both age and sex were added into a model (SC or FC = β₀ + β₁ * age + β₂ * sex) using stepwise selection input method.

T4: Pearson correlations between each SC measure and FC in each network. In addition, partial correlations were performed to control for the effect of age that may affect both SC and FC.

Results

SC and FC Changes with Age (T1)

A representative example of the relationships between SC and FC measure changes with age for the ventral attention network is shown in Figure 3. All four SC measures (i.e., FA, MD, mean tract length, number of streamlines) were significantly correlated with age in all networks (Table 2). FA, mean tract length, and number of streamlines were decreased with age, while MD was increased with age. On the other hand, there was a general trend of FC decrease with age in all networks. FC was significantly negatively correlated with age in five out of seven RSNs (somato-motor, dorsal attention, ventral attention, limbic, and frontal-parietal).

FIGURE 3

Figure 3. Representative plots of SC and FC measures for all subjects (blue asterisk, male; red circle, female) for the ventral attention network. The blue and red dash lines in each plot represent the linear trajectory to model the change of SC and FC measures with age for all male and female subjects, respectively. The Pearson correlation coefficient (r), and the corresponding uncorrected p-value, for each connectivity measure with age by sex is also shown.

TABLE 2

Table 2. Pearson correlation coefficients (r) and slopes of the linear regression lines of SC (i.e., MD, FA, mean tract length, number of streamlines) and FC (i.e., Fisher's r-to-z transformed Pearson correlation coefficient) measures with age in each resting-state network (RSN).

Sex Differences in SC and FC Changes with Age (T2)

The slope of the linear regression lines for all SC and FC measures with age were not statistically different between male and female subjects in any network (Table 3).

TABLE 3

Table 3. Pearson correlation coefficients (r) and slopes of the linear regression lines of SC (i.e., MD, FA, mean tract length, number of streamlines) and FC (i.e., Fisher's r-to-z transformed Pearson correlation coefficient) measures with age for male and female subjects in each resting-state network (RSN).

Other Predictors of SC and FC (T3)

Multiple linear regression analysis found that both the age and sex terms were significant for FA, mean tract length, and number of streamlines in all networks except in the visual network (only the age term was significant for FA). For MD and FC, only the age term was significant in all cases except in the visual network (FC remained unchanged with age; Table 4).

TABLE 4

Table 4. Multiple linear regression of SC (i.e., MD, FA, mean tract length, number of streamlines) and FC (i.e., Fisher's r-to-z transformed Pearson correlation coefficient) measures in each resting-state network (RSN).

Relationships between SC and FC (T4)

An example of Pearson correlations between SC and FC in the somato-motor network is shown in Figure 4. The four SC measures were not correlated with FC in all networks, except that the mean tract length in the somato-motor and ventral attention networks, and number of streamlines in the frontal-parietal network, were significantly related to FC (Table 5). Partial correlations did not reveal any significant relationship between SC and FC after controlling for the effect of age.

FIGURE 4

Figure 4. Plots showing the linear relationship between SC (i.e., MD, FA, mean tract length, number of streamlines) and FC measures for the somato-motor network. The Pearson correlation coefficient (r), and the corresponding uncorrected p-value, for each SC measure with FC is also shown. In this network, only mean tract length was significantly correlated with FC.

TABLE 5

Table 5. Pearson and partial (removing the age effect) correlation coefficients (r) of SC (i.e., MD, FA, tract length, num streamlines) with FC (i.e., Fisher's r-to-z transformed Pearson correlation coefficient) measures in each resting-state network (RSN).

Discussion

In this study, the four SC measures demonstrated significant age-related changes in all seven RSNs across the healthy adult lifespan, while FC demonstrated significant age-related changes in four of the seven networks. In general, SC measures were not related to FC suggesting that WM microstructure as inferred from the SC measures derived from DTI data do not correlate with the corresponding cortical FC changes within a RSN.

SC and FC Changes with Age

The observed correlations of SC and FC measures with age among all seven networks are consistent with published independent DTI (Westlye et al., 2010; Lebel et al., 2012; Chen et al., 2013) and rs-fMRI (Mevel et al., 2013; Cao et al., 2014; Fjell et al., 2015a) human brain aging studies. In these studies, other trajectories, such as quadratic or Poisson trajectories, were used to model age-related changes in SC (FA, MD, AD, RD) or FC across wider age ranges that included children and adolescents. Visual inspection of these published results for SC and FC changes in adulthood (≥18 years, the age range of participants in this study) shows the trend to more closely resemble a linear trajectory, which serves as a good approximation for SC and FC measures in this work. However, longitudinal and/or larger studies that include child and adolescent participants may better be able to elucidate the exact trajectories. Results for the other two SC measures in this study (i.e., mean tract length and number of streamlines) are also consistent with literature. A recent study that analyzed DTI data from 121 subjects between age 4 and 40 years demonstrated significant decrease in the number of streamlines, and the loss of streamlines occurred earlier in females than in males (Lim et al., 2015). The mean tract length decrease observed in this study was also consistent with two earlier studies that demonstrated WM fiber bundle length decrease in healthy adults over 50 years of age (Baker et al., 2014; Behrman-Lay et al., 2015).

Sex Differences in SC and FC

A significant effect of sex in the regression model was observed for FA, tract length, and number of streamlines across all networks except for FA in the visual network (Table 4). These three SC measures demonstrated significantly higher mean values in males. There was no significant difference between male and female subjects observed for MD. Our results agree with sex differences in FA (higher in male over female) of certain major WM tracts reported in previous studies (Hsu et al., 2008; Lebel et al., 2012). Higher mean FC was observed for males, though the sex term in the regression model was not significant for FC in any network. Sex differences in FC in healthy adults have been reported in previous studies using graph theory (Cao et al., 2014; Scheinost et al., 2015), however, not in previous rs-fMRI aging studies in normal adults that used the same FC metric as this work (i.e., z-transformed correlation of the averaged BOLD time signal; Wang et al., 2012; Fjell et al., 2015a).

Relationship between SC and FC

In general, SC measures for WM tracts were not correlated with FC. However, some measures had correlations that were significant at a trend level, before multiple comparison correction (i.e., p < 0.05), suggesting that a weak relationship may exist between SC and FC in this sample (Table 5). Spatial averaging of metrics within a network could mask possible relationships between SC and FC in more spatially localized areas. Further analysis was performed to address this concern. All combinations of pairs of regions in each network were used as seed regions in the tractography algorithm, but only those with delineated WM tracts connecting them were considered in this analysis. The more spatially specific correlations of SC (i.e., FA, MD, tract length, number of streamlines) and FC for relevant pairs of regions in each network are shown in Table 6. Both network-averaged and more spatially specific results showed a lack of significant relationship between SC and FC, which is contrary to our hypothesis, and suggest that the changes in WM microstructure do not play a significant role in the corresponding changes in FC within large-scale RSNs.

TABLE 6

Table 6. Pearson correlation coefficients (r) of SC with FC measures for individual pairs of regions in each resting stating network (RSN).

The relationship between SC and FC has been explored previously within younger and elderly healthy participants in different studies (Andrews-Hanna et al., 2007; Fjell et al., 2015b; Marstaller et al., 2015; Hirsiger et al., 2016), however, results have been inconsistent. The first study by Andrews-Hanna et al. demonstrated a significant positive linear relationship between FC and FA (FC measured for the prefrontal cortex and retrosplenial/posterior cingulate cortex in the default mode network, and FA measured in a large WM region that included tracks connecting anterior to posterior regions) in elderly subjects. Similarly, subsequent studies by Fjell et al. and Marstaller et al. demonstrated that in younger subjects, FA of the uncinate was negatively correlated with FC between hippocampus and cortical regions (Fjell et al.), and global FA was negatively correlated with FC in the prefrontal regions of frontal-parietal and salience networks (FPN and SN; Marstaller et al.). In elderly subjects, MD of the cingulate bundle was positively correlated with FC between caudate and cortical regions (Fjell et al.), and global MD was positively correlated with FC in the prefrontal regions of FPN and SN (Marstaller et al.). While these studies demonstrated significant relationships between FC and FA or MD, a recent study by Hirsiger et al. did not find any significant relationships between FC (measured between posterior cingulate cortex and medial prefrontal cortex) and either AD, MD, RD, or FA (measured from the cingulum bundle) in healthy elderly subjects. During the preparation of this manuscript, another study similar to our present study was published (Fjell et al., 2016) and demonstrated modest relationship between SC and FC measures only within certain regions of the default mode network. It should be pointed out that the results reported by Fjell et al. were essentially derived from data obtained between two age groups of participants between 20 and 40 years and above 60 years, with only one participant around 50 years of age.

The results from the present cross-sectional study across the adult lifespan with relatively even number of participants in each decade of adulthood, rather than comparing between groups of elderly and younger subjects, show that the change in WM microstructure is not significantly related to the corresponding change in FC within the seven RSNs tested, contrary to our expectation.

Study Limitations

This study has a number of limitations. First there is a gender imbalance, which is most pronounced in the middle aged (40–49 years) and elderly (>60 years) groups. Sex differences were observed in SC in most networks, but there were no significant age-sex interactions. Therefore, we believe that results are generalizable to both women and men, though future studies should investigate sex differences further. Second, the question related to whether WM microstructural changes in normal aging precede FC changes cannot be addressed in this cross-sectional study. Future longitudinal studies will need to test age-related changes within individuals to better elucidate the relationships of SC and FC. Third, we used template-defined RSNs that do not allow for potential changes in the topography of these networks with age. It would be an interesting avenue for future study to incorporate spatial and temporal lifespan changes into networks defined by age group or at an individual level. Furthermore, we have not comprehensively investigated the relationship of SC and FC in all RSNs, for example, the salience network that was not part of the template was not included in the analysis. Future studies are needed to explore other networks that have been omitted in this work. Finally, despite the importance of DTI based tractography algorithms to provide quantitative measures to characterize WM microstructure integrity and architecture, these algorithms are limited to provide accurate delineation of the anatomical structural connections (Thomas et al., 2014). This latter observation may lead to spurious or missed WM tracts that belong to RSNs. Furthermore, one difference across studies of SC and FC relates to the choice of tractography algorithm used. It was shown recently that both deterministic and probabilistic tractography algorithms yielded similar relationships between SC with FC in regions associated with the default mode network (Khalsa et al., 2014). However, another study showed that both deterministic and probabilistic tractography algorithms underestimated the corticospinal tract connections to the sensorimotor cortex but a more complex algorithm based on constrained spherical deconvolution (CSD) reliably delineate the tracts that closely resembled to known anatomy of that brain region (Farquharson et al., 2013). Future studies should compare different tractography approaches and evaluate the impact on the relationships between SC and FC in normal aging.

Conclusions

A multi-modal analysis approach using DTI and rs-fMRI data was used to investigate SC and FC within seven commonly described RSNs. SC measures demonstrated significant age-related changes in all networks, while FC demonstrated significant age-related changes in four of the seven networks. Despite significant age correlations in both SC and FC parameters, however, these were in general not significantly related to each other, suggesting that the change in WM microstructure measures with age is too weak to reflect the corresponding cortical FC change in resting-state networks. These results help further understand healthy brain aging, and lay the foundation for future studies to investigate age-related changes in connectivity in adults with neurodegenerative diseases.

Ethics Statement

This study was carried out in accordance with the recommendations of University of Calgary Conjoint Health Research Ethics Board (CHREB) with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the CHREB.

Author Contributions

AT performed analyses and wrote the manuscript. CL, SB, BG, and RS provided guidance in DTI tractography and resting state fMRI analyses. MH assisted in data processing. CM and RF designed study and reviewed manuscript.

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 an operating grant from the Canadian Institutes of Health Research (CIHR; MOP-333931 and MOP-106571) (RF). AT received post-doctoral fellowship awards from the NSERC CREATE International Industrial Imaging (I3T) program and the Hotchkiss Brain Institute (HBI). RF is the Hopewell Professor of Brain Imaging. We acknowledge Linda B Andersen, Ph.D. for subject recruitment, the technical support provided by the Calgary Image Processing and Analysis Centre (CIPAC) and the high-performance computing resources provided by Compute Canada/WestGrid.

References

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Auzias, G., Coulon, O., and Brovelli, A. (2016). MarsAtlas: a cortical parcellation atlas for functional mapping. Hum. Brain Mapp. 37, 1573–1592. doi: 10.1002/hbm.23121

CrossRef Full Text | Google Scholar

Baker, L. M., Laidlaw, D. H., Conturo, T. E., Hogan, J., Zhao, Y., Luo, X., et al. (2014). White matter changes with age utilizing quantitative diffusion MRI. Neurology 83, 247–252. doi: 10.1212/WNL.0000000000000597

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Behrman-Lay, A. M., Usher, C., Conturo, T. E., Correia, S., Laidlaw, D. H., Lane, E. M., et al. (2015). Fiber bundle length and cognition: a length-based tractography MRI study. Brain Imaging Behav. 9, 765–775. doi: 10.1007/s11682-014-9334-8

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Errangi, B., Li, L., Glasser, M. F., Westlye, L. T., Fjell, A. M., et al. (2013). Brain aging in humans, chimpanzees (Pan troglodytes), and rhesus macaques (Macaca mulatta): magnetic resonance imaging studies of macro- and microstructural changes. Neurobiol. Aging 34, 2248–2260. doi: 10.1016/j.neurobiolaging.2013.03.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, L., Li, H., Zhuo, J., Zhang, Y., Wang, J., Chen, L., et al. (2016). The human brainnetome atlas: a new brain atlas based on connectional architecture. Cereb. Cortex 26, 3508–3526. doi: 10.1093/cercor/bhw157

PubMed Abstract | CrossRef Full Text | Google Scholar

Farquharson, S., Tournier, J. D., Calamante, F., Fabinyi, G., Schneider-Kolsky, M., Jackson, G. D., et al. (2013). White matter fiber tractography: why we need to move beyond DTI. J. Neurosurg. 118, 1367–1377. doi: 10.3171/2013.2.JNS121294

PubMed Abstract | CrossRef Full Text | Google Scholar

Fjell, A. M., Sneve, M. H., Grydeland, H., Storsve, A. B., Amlien, I. K., Yendiki, A., et al. (2016). Relationship between structural and functional connectivity change across the adult lifespan: a longitudinal investigation. Hum. Brain Mapp. 38, 561–573. doi: 10.1002/hbm.23403

PubMed Abstract | CrossRef Full Text | Google Scholar

Fjell, A. M., Sneve, M. H., Grydeland, H., Storsve, A. B., de Lange, A. G., Amlien, I. K., et al. (2015a). Functional connectivity change across multiple cortical networks relates to episodic memory changes in aging. Neurobiol. Aging 36, 3255–3268. doi: 10.1016/j.neurobiolaging.2015.08.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Fjell, A. M., Sneve, M. H., Storsve, A. B., Grydeland, H., Yendiki, A., and Walhovd, K. B. (2015b). 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, L. L., and Wu, T. (2016). The study of brain functional connectivity in Parkinson's disease. Transl. Neurodegener. 5:18. doi: 10.1186/s40035-016-0066-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gobbi, D., Lu, Q., Frayne, R., and Salluzzi, M. (2012). “Cerebra-WML: a rapid workflow for quantification of white matter hyperintensities,” in Canadian Stroke Congress (Calgary, AB).

Hawkins, K. M., Goyal, A. I., and Sergio, L. E. (2015). Diffusion tensor imaging correlates of cognitive-motor decline in normal aging and increased Alzheimer's disease risk. J. Alzheimers. Dis. 44, 867–878. doi: 10.3233/JAD-142079

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Holm, S. (1979). A simple sequentially rejective multiple test procedure. Scand. J. Stat. 6, 65–70.

Google Scholar

Honey, C. J., Sporns, O., Cammoun, L., Gigandet, X., Thiran, J. P., Meuli, R., et al. (2009). Predicting human resting-state functional connectivity from structural connectivity. Proc. Natl. Acad. Sci. U.S.A. 106, 2035–2040. doi: 10.1073/pnas.0811168106

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, J. L., Leemans, A., Bai, C. H., Lee, C. H., Tsai, Y. F., Chiu, H. C., et al. (2008). Gender differences and age-related white matter changes of the human brain: a diffusion tensor imaging study. Neuroimage 39, 566–577. doi: 10.1016/j.neuroimage.2007.09.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Iturria-Medina, Y., and Evans, A. C. (2015). On the central role of brain connectivity in neurodegenerative disease progression. Front. Aging Neurosci. 7:90. doi: 10.3389/fnagi.2015.00090

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkinson, M., Beckmann, C. F., Behrens, T. E., Woolrich, M. W., and Smith, S. M. (2012). FSL. Neuroimage 62, 782–790. doi: 10.1016/j.neuroimage.2011.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Khalsa, S., Mayhew, S. D., Chechlacz, M., Bagary, M., and Bagshaw, A. P. (2014). The structural and functional connectivity of the posterior cingulate cortex: comparison between deterministic and probabilistic tractography for the investigation of structure-function relationships. Neuroimage 102 (Pt 1), 118–127. doi: 10.1016/j.neuroimage.2013.12.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Koch, M. A., Norris, D. G., and Hund-Georgiadis, M. (2002). An investigation of functional and anatomical connectivity using magnetic resonance imaging. Neuroimage 16, 241–250. doi: 10.1006/nimg.2001.1052

CrossRef Full Text | Google Scholar

Lebel, C., Gee, M., Camicioli, R., Wieler, M., Martin, W., and Beaulieu, C. (2012). Diffusion tensor imaging of white matter tract evolution over the lifespan. Neuroimage 60, 340–352. doi: 10.1016/j.neuroimage.2011.11.094

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, A., Ratnarajah, N., Tuan, T. A., Chen, S. H., and Qiu, A. (2015). Adaptation of brain functional and structural networks in aging. PLoS ONE 10:e0123462. doi: 10.1371/journal.pone.0123462

PubMed Abstract | CrossRef Full Text | Google Scholar

Lim, S., Han, C. E., Uhlhaas, P. J., and Kaiser, M. (2015). Preferential detachment during human brain development: age- and sex-specific structural connectivity in diffusion tensor imaging (DTI) data. Cereb. Cortex 25, 1477–1489. doi: 10.1093/cercor/bht333

PubMed Abstract | CrossRef Full Text | Google Scholar

Madden, D. J., Bennett, I. J., Burzynska, A., Potter, G. G., Chen, N. K., and Song, A. W. (2012). Diffusion tensor imaging of cerebral white matter integrity in cognitive aging. Biochim. Biophys. Acta 1822, 386–400. doi: 10.1016/j.bbadis.2011.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Marstaller, L., Williams, M., Rich, A., Savage, G., and Burianova, H. (2015). Aging and large-scale functional networks: white matter integrity, gray matter volume, and functional connectivity in the resting state. Neuroscience 290, 369–378. doi: 10.1016/j.neuroscience.2015.01.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Mevel, K., Landeau, B., Fouquet, M., La Joie, R., Villain, N., Mezenge, F., et al. (2013). Age effect on the default mode network, inner thoughts, and cognitive abilities. Neurobiol. Aging 34, 1292–1301. doi: 10.1016/j.neurobiolaging.2012.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Onoda, K., Ishihara, M., and Yamaguchi, S. (2012). Decreased functional connectivity by aging is associated with cognitive decline. J. Cogn. Neurosci. 24, 2186–2198. doi: 10.1162/jocn_a_00269

PubMed Abstract | CrossRef Full Text | Google Scholar

Pievani, M., Filippini, N., van den Heuvel, M. P., Cappa, S. F., and Frisoni, G. B. (2014). Brain connectivity in neurodegenerative diseases–from phenotype to proteinopathy. Nat. Rev. Neurol. 10, 620–633. doi: 10.1038/nrneurol.2014.178

PubMed Abstract | CrossRef Full Text | Google Scholar

Power, J. D., Mitra, A., Laumann, T. O., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2014). Methods to detect, characterize, and remove motion artifact in resting state fMRI. Neuroimage 84, 320–341. doi: 10.1016/j.neuroimage.2013.08.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Salimi-Khorshidi, G., Douaud, G., Beckmann, C. F., Glasser, M. F., Griffanti, L., and Smith, S. M. (2014). Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers. Neuroimage 90, 449–468. doi: 10.1016/j.neuroimage.2013.11.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Scheinost, D., Finn, E. S., Tokoglu, F., Shen, X., Papademetris, X., Hampson, M., et al. (2015). Sex differences in normal age trajectories of functional brain networks. Hum. Brain Mapp. 36, 1524–1535. doi: 10.1002/hbm.22720

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, C., Ye, F. Q., Irfanoglu, M. O., Modi, P., Saleem, K. S., Leopold, D. A., et al. (2014). Anatomical accuracy of brain connections derived from diffusion MRI tractography is inherently limited. Proc. Natl. Acad. Sci. U.S.A. 111, 16574–16579. doi: 10.1073/pnas.1405672111

PubMed Abstract | CrossRef Full Text | Google Scholar

Tyszka, J. M., Kennedy, D. P., Adolphs, R., and Paul, L. K. (2011). Intact bilateral resting-state networks in the absence of the corpus callosum. J. Neurosci. 31, 15154–15162. doi: 10.1523/JNEUROSCI.1453-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Uddin, L. Q., Mooshagian, E., Zaidel, E., Scheres, A., Margulies, D. S., Kelly, A. M., et al. (2008). Residual functional connectivity in the split-brain revealed with resting-state functional MRI. Neuroreport 19, 703–709. doi: 10.1097/WNR.0b013e3282fb8203

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Heuvel, M. P., Mandl, R. C., Kahn, R. S., and Hulshoff Pol, H. E. (2009). Functionally linked resting-state networks reflect the underlying structural connectivity architecture of the human brain. Hum. Brain Mapp. 30, 3127–3141. doi: 10.1002/hbm.20737

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Su, L., Shen, H., and Hu, D. (2012). Decoding lifespan changes of the human brain using resting-state functional connectivity MRI. PLoS ONE 7:e44530. doi: 10.1371/journal.pone.0044530

PubMed Abstract | CrossRef Full Text | Google Scholar

Westlye, L. T., Walhovd, K. B., Dale, A. M., Bjornerud, A., Due-Tonnessen, P., Engvig, A., et al. (2010). Life-span changes of the human brain white matter: diffusion tensor imaging (DTI) and volumetry. Cereb. Cortex 20, 2055–2068. doi: 10.1093/cercor/bhp280

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeo, B. T., Krienen, F. M., Sepulcre, J., Sabuncu, M. R., Lashkari, D., Hollinshead, M., et al. (2011). The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. 106, 1125–1165. doi: 10.1152/jn.00338.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaiontz, C. (2013–2015). Real Statistics Using Excel. Available online at: http://www.real-statistics.com/ (Accessed Aug 10, 2016).