Age-Related Decline in the Variation of Dynamic Functional Connectivity: A Resting State Analysis
Yuanyuan Chen1,2 
Weiwei Wang2,3 Xin Zhao2,3* Miao Sha2,3 Ya’nan Liu2,3 Xiong Zhang2,3 Jianguo Ma1 Hongyan Ni4 
Dong Ming2,3*- 1College of Microelectronics, Tianjin University, Tianjin, China
- 2Tianjin International Joint Research Center for Neural Engineering, Academy of Medical Engineering and Translational Medicine, Tianjin University, Tianjin, China
- 3Department of Biomedical Engineering, College of Precision Instruments and Optoelectronics Engineering, Tianjin University, Tianjin, China
- 4Department of Radiology, Tianjin First Center Hospital, Tianjin, China
Normal aging is typically characterized by abnormal resting-state functional connectivity (FC), including decreasing connectivity within networks and increasing connectivity between networks, under the assumption that the FC over the scan time was stationary. In fact, the resting-state FC has been shown in recent years to vary over time even within minutes, thus showing the great potential of intrinsic interactions and organization of the brain. In this article, we assumed that the dynamic FC consisted of an intrinsic dynamic balance in the resting brain and was altered with increasing age. Two groups of individuals (N = 36, ages 20–25 for the young group; N = 32, ages 60–85 for the senior group) were recruited from the public data of the Nathan Kline Institute. Phase randomization was first used to examine the reliability of the dynamic FC. Next, the variation in the dynamic FC and the energy ratio of the dynamic FC fluctuations within a higher frequency band were calculated and further checked for differences between groups by non-parametric permutation tests. The results robustly showed modularization of the dynamic FC variation, which declined with aging; moreover, the FC variation of the inter-network connections, which mainly consisted of the frontal-parietal network-associated and occipital-associated connections, decreased. In addition, a higher energy ratio in the higher FC fluctuation frequency band was observed in the senior group, which indicated the frequency interactions in the FC fluctuations. These results highly supported the basis of abnormality and compensation in the aging brain and might provide new insights into both aging and relevant compensatory mechanisms.
Introduction
Normal aging in the human brain refers to degradation phenomena that occur in brain structures, brain function and brain morphology with increasing age, indicating that a certain degree of senior brain dysfunction will occur (Hedden and Gabrieli, 2004; Fjell et al., 2014). Considering the increasing size of the aging population, the incidences of diseases that are highly associated with age, such as Alzheimer’s (Kern and Behl, 2009; Mosher and Wyss-Coray, 2014) and Parkinson’s (Xu et al., 2012; Reeve et al., 2014), are also increasing. Until now, the mechanism of aging has remained unclear, and further investigation of brain aging could greatly help in managing problems associated with both aging and disease.
Functional connectivity (FC) based measures of the resting state functional magnetic resonance imaging (rs-fMRI), which reflect the coherence between temporal fluctuations across brain regions, are organized into distinct systems or networks (Damoiseaux et al., 2008; Zuo et al., 2010). Most studies of FC have focused on the decline of specific functional systems, such as the default mode network (DMN; Raichle et al., 2001; Damoiseaux et al., 2008; Wang et al., 2010; Ferreira and Busatto, 2013), or have focussed on other specific brain networks or regions, such as the language system (Zou et al., 2012), subcortical regions (Yi et al., 2015) or the motor system (Coynel et al., 2010; De Vico Fallani et al., 2013). Increasing evidence has shown that the decline in cognitive function associated with aging is related to changes in communication between different brain regions and subsystems (Andrews-Hanna et al., 2007; Sambataro et al., 2010), even in the resting state (Shehzad et al., 2009; Meindl et al., 2010; Guo et al., 2012; Zuo et al., 2013). Despite this progress, how brain systems cooperate to handle aging-associated declines remains unclear, especially considering the averaging of complex spatiotemporal phenomena during a period of time (Hutchison et al., 2013a).
Traditionally, functional connectivities derived from fMRI data are computed using signals across the entire scan time; it is assumed that the functional connectivities among the brain regions are static during the duration of the resting time (Handwerker et al., 2012; Zuo et al., 2013). However, recent work has shown that FC is temporally dynamic (Chang and Glover, 2010; Calhoun et al., 2014) even at rest. This dynamic FC, which varies over a timeframe of seconds, could be highly related to unconstrained mental activity during the resting state (Hutchison et al., 2013a; Allen et al., 2014; Zalesky and Breakspear, 2015) and even under anesthesia (Hutchison et al., 2013b). A widely applied method for analyzing temporal dynamics is the sliding window correlation method (Sakolu et al., 2010; Hutchison et al., 2013a; Di and Biswal, 2015). A series of FC matrices was obtained using this method, which showed the time-varying connectivity network. Thus, researchers have started to perform dynamic FC investigations of mild cognition impairment (Wee et al., 2016), epilepsy (Liu et al., 2016), schizophrenia (Rashid et al., 2014; Du et al., 2016), major depressive disorder (Demirta et al., 2016) and normal development (Sakolu et al., 2010; Rashid et al., 2014; Qin et al., 2015).
The convergent results of previous studies have suggested that the dynamic resting state FC is highly intrinsic and physiologically relevant. Several studies have reported that the FC states revealed by changes in connectivity over the course of the scan can be sensitive to changes related to neurological disorders (Sakolu et al., 2010; Li et al., 2014; Leonardi and Van De Ville, 2015; Ou et al., 2015; Shakil et al., 2016). Increasing efforts have been directed toward using functional microstates and their transmissions to depict the working mechanisms of the brain (Allen et al., 2014; Shakil et al., 2016). Microstate transmissions could be the bases of integration and segregation between different brain networks or cognitive resources (Hansen et al., 2015; Yu et al., 2015; Shakil et al., 2016). Previous work has shown that aging impacts not only within-network connectivity but also the integration and segregation of different brain networks (Ferreira and Busatto, 2013). Advancing age induces increased reorganization to establish compensatory mechanisms or plasticity that counteract the aging process (Meunier et al., 2014; Sala-Llonch et al., 2015; Sugiura, 2016). Segregation and integration are the bases of reorganization of the brain connectivity network. The investigation of dynamic FC in the resting state may provide new insights into communication between various cognitive resource pools in the aging brain. The derived patterns of temporal variation in FC thus reflect the interactions of the brain functional networks and are therefore expected to facilitate our understanding of the mechanisms that underlie mental diseases.
We expected that the fluctuations of resting FC comprise a dynamic balance that maintains the intrinsic connectivity patterns in the brain. The dynamic balance of FC allows us to capture the interactions between all of the subsystems and the basic states of brain connectivity. Since the FC levels and patterns are age-related, this dynamic balance and the serial connectivity networks must also change with increasing age. A previous study (Leonardi and Van De Ville, 2015) suggested that the spontaneous fluctuations in the FC have frequency dependence and result from the interactions of various frequency components associated with neural activities. Limited resources and decreased processing speed can indicate performance during aging; thus, we hypothesized that the dynamic FC can provide clues for the capacity and efficiency of the connectivity states that transfer and present aging features. This article focuses on revealing the effects of aging on the time-varying FC of the brain in the resting state. Using the sliding window correlation method, the resting state fMRI data from two groups of young and senior healthy individuals were processed to construct dynamic FC matrices. We expected that the variation and the frequency spectrum of the FC fluctuations were the important bases of the dynamic balance and were highly related to aging.
Materials and Methods
Participants and fMRI Data Acquisition
All resting-state fMRI data used in this study were obtained from the NKI-Rockland Sample (NKI-RS1), which is provided by the Nathan Kline Institute (NKI, Orangeburg, NY, USA) and is available online in a public database. To study the changes in the dynamic characteristics of FC that resulted from normal brain aging, we collected fMRI data from 68 healthy subjects who were organized into two groups: 36 young subjects were assigned to one group (mean age, 28.1 years; range, 20–35 years; 24 male), and 32 senior subjects were assigned to the other group (mean age, 70.6 years; range, 60–85 years; 15 male). According to the demographic information provided by the NKI-RS data set, there was a remarkable difference in the age of the participants, but no significant differences in either gender or hand dominance between the two groups.
Resting-state fMRI data were collected using an echo-planar imaging (EPI) sequence on a 3.0 T SIMENS Trio scanner. The scanning parameter settings were as follows: TR/TE = 2500/30 ms, flip angle (FA) = 80°, field of view (FOV) = 216 × 216 mm2, voxel size = 3 × 3 × 3 mm3, number of slices = 38, scan time = 650 s time points = 260. During the data acquisition, the subjects were instructed to keep their eyes closed and to stay awake. High-resolution T1-weighted images were also acquired using the magnetization-prepared rapid gradient echo (MPRAGE) sequence. The acquisition parameter settings were as follows: TR/TE = 2500/3.5 ms, FA = 8°, FOV = 256 × 256 mm2, voxel size = 1 × 1 × 1 mm3, slice = 192.
Data Preprocessing
Functional images were preprocessed using the Connectome Computation System (CCS2). The CCS designed by Zuo et al. (2013) provides a computational platform for multimodal neuroimaging brain connectomics computations by integrating the functionalities of AFNI, FSL and FreeSurfer (Zuo et al., 2013; Betzel et al., 2014; Cao et al., 2014). The functional preprocessing included the following: the first ten functional volumes were discarded to allow for signal equilibration; slice timing was corrected using the middle slice as the reference frame; 3D geometrical displacement was used to correct for head motion; and 4D global mean-based intensity correction was performed. In addition, the Friston-24 model was used to remove micro-level motion artifacts (Friston et al., 1996) and nuisance regressors; for instance, the individual white matter and the cerebrospinal fluid (CSF) mean signals were regressed out. The functional data were also temporal band-pass filtered (0.01–0.1 Hz) and detrended (both linear and quadratic trends). Finally, spatial smoothing was performed with a Gaussian filter kernel (FWHM = 6 mm). The structural processing steps were as follows: the image noise was removed using a spatially adaptive non-local means filter and brain surface reconstruction; the individual functional space was spatially normalized to the MNI152 standard brain space; a customized group T1 template in the standard space was generated to reduce the error term that resulted from the image registration and bias in the template selection; and the fMRI images in the native space of each subject was registered to the standard space with a final resolution of 3 mm.
Dynamic Functional Connectivity Network Construction
Because time-varying FC is complicated and differs from static FC, the recruitment of more regions in the associated networks could help to provide more precise information. A total of 142 regions that covered the cingulo-opercular network (CON), DMN, fronto-parietal network (FPN), occipital network (OCC) and sensorimotor network (SMN) as defined by Dosenbach et al. (2010) were selected. The cerebellum network was neglected because we sought to examine only the effects of aging on the higher-order brain network interactions and the dynamics of brain cognition and perception. Among these networks, OCC and SMN are involved in perception and primary visual and motion processing, respectively; the other three networks are important to higher-order cognitive functions. In each of these brain regions, time courses were extracted and averaged over a spherical region of interest (ROI) with a diameter of 6 mm. Then, a dynamic FC network was estimated using the sliding window Pearson correlation method, which yielded a series of 142 × 142 correlation matrices. We used a fixed-length rectangle window (width = 24 × TRs = 60 s), and the window was shifted by 1 TR. The obtained correlation series were then Fisher-Z transformed and low-pass filtered with a cut-off frequency of 1/w Hz. All of these network matrices were vectorized to simplify the analysis.
Phase Randomization
As suggested previously (Hutchison et al., 2013a; Hindriks et al., 2016), phase randomization analysis was used to explore the dependability of the dynamic FC fluctuation. The processed rs-fMRI time courses from the senior and young groups were phase randomized into new time courses in which the frequency spectra of the bold signals were invariable. We called the phase processed data the null group, which was then compared with the senior and young groups. Phase randomization was conducted for all parameters except amplitude (Friston et al., 1994; Handwerker et al., 2012), which could preserve the temporal correlation properties. These steps were taken to allow for assessing the dependability of the dynamic fluctuations and to verify whether the FC fluctuations over the rs-fMRI involved specific neural activities.
Functional Connectivity Variation (FCV)
The dynamic functional connectivity variation (FCV) was calculated as the standard variation of the dynamic FC series. In this approach, the stability of the FC fluctuation over time is quantitatively measured and compared between brain region pairs. Previous studies of resting-state fMRI have demonstrated that some intrinsic neural activities are related to the variations in FC. These findings likely suggest the internal mechanism of the resting-state fMRI; thus, the FCV matrix was calculated for each subject. The original FCV matrices and phase randomized FCV matrices of both groups were statistically analyzed and compared with the averages of each group using a one-sample t-test. Thus, we could easily examine network modularization and the FCV of each network or connection.
Frequency Spectrum Analysis
A sliding rectangle window was used with a low-pass filtering effect on the functional fluctuations. The cut-off frequency was 1/w (1/60 Hz = 0.018 Hz). We assumed that the frequency spectrum of the dynamic FC fluctuation would change with aging; thus, we sought to specify the age-related changes in frequency or energy as age increased. The frequency band of the FC fluctuation was divided equally into two frequency bands, 0–1/2w Hz and 1/2w–1/w Hz, within which the fluctuation energies were calculated with a Fourier transform. Then, the energy ratio of the two frequency bands was calculated as the energy of the lower frequency band dividing the energy of the higher frequency band.
Permutation Tests
To obtain robust results on aging-related variations within and between groups, a non-parametric permutation test with 5000 randomizations of the group labels was utilized for all measures described above. We defined the t-statistics between the two groups as the difference measurement, yielding a distribution of t-statistics after 5000 randomizations. Then, p = 0.001 was set as the threshold of significance. To better understand the results, both connectivity-based and network-averaged indices were examined. With the non-parametric approach, the size of the Type I error was guaranteed to be set at the prescribed significance level. The permutation test demonstrated an excellent ability to differentiate between different profiles, even when those profiles appeared to be highly similar.
Sliding Window Length Analysis
Previous works (Hutchison et al., 2013a; Hindriks et al., 2016) have much discussed the influence of sliding window parameter settings. However, there have been no definite conclusions about the optimal window length. In this article, we also carefully assessed the influence of the window length on the variation and frequency spectrum of the dynamic FC. A sequence of window lengths from 2 to 256 time points was selected to examine the FC variation between the mean dynamic FC and the static FC. Another sequence of window lengths from 10 to 70 time points was selected to obtain the frequency characteristics of the dynamic function connectivity time series.
Results
Phase Randomized and Within-Group Analyses
The within group and between group comparisons were conducted after passing the normality tests. Both groups lost the network organization pattern after phase randomization, and all of the within- and between-network connections showed similar FC variations (Figures 1A,D). From the original data obtained from both groups, clear modularization could be identified from the lower within-network variation and the higher between-network variation (Figures 1B,E). In the young group, higher variation of the dynamic FC was found in the DMN-related inter-network connections and the connections between the FPN and OCC. By contrast, in the senior group, only the DMN-related inter-network connections showed high variation. In both groups, however, the inter-network variations of the connections between the CON and SMN were clearly lower than the corresponding group averaged values; a similar result was obtained for the CON and OCC. Both the within- and between-network variations based on the original data were significantly higher than those of the phase randomized data for almost all of the connections (Figures 1C,F).
Figure 1. Illustration of the dynamic functional connectivity (FC) variation patterns within groups. Pictures (A,D) show the one-sample t-test results of senior and young null groups with phase randomization processed data (FDR p-value < 0.05; actual p-value < 0.0010); Pictures (B,E) show the one-sample t-test results of senior and young groups with the original data (FDR p-value < 0.05; actual p-value < 0.0011); Pictures (C,F) show two-sample t-test results between the original and phase randomized data within the senior and young groups (FDR p-value < 0.05; actual p-value < 0.0008). Red indicates a higher-than-average level in the one-sample t-test and a higher-than-null group in the two-sample t-test.
Age-Related Changes in Static FC
Compared with the young group, all connections indicated decreased within-network FC in the senior group, and fewer within-network functional connections showed increases within the CON and OCC (Figure 2A). Most of the between-network connections, especially between the CON, DMN, FPN and OCC, showed increased FC in the senior group compared with those of the young group. From the 3D view in Figure 2B, more of the connections crossing the cerebral hemispheres were changed compared with the connections within one side. Most of the connections that both crossed hemispheres and occurred within a hemisphere were located between the posterior-anterior brain; in these connections, several distinct, intensively connected nodes were found. In the networks that were averaged and assessed, all or most of the between-network connectivities increased and the inner network functional connectivities decreased in the senior group compared with those of the young group (Figure 2C). However, only the DMN and SMN showed significantly decreased within-network FC.
Figure 2. Illustration of the two-sample t-test results of the static FC compared between the senior and young groups. Pictures (A,C) depict the connections level and the network averaging level, respectively; Picture (B) is a 3D view of the results shown in picture (A) (constructed using the BrainNet Viewer). Red indicates an increase in the senior group; blue indicates a decrease in the senior group. The significant level is a p-value < 0.001.
Age-Related Changes in FC Variation
Most of the significantly changed connections shown in Figure 3 were between-network connections. The few connections with increased connectivity were located between the hemispheres, and the changed connections shared a similar location with the FC connection, which crossed both hemispheres, showed a posterior-anterior distribution and was also changed (Figure 3B). Several intensively connected nodes were also obvious in the prefrontal and occipitotemporal regions. The averaging analysis of the networks showed that these decreases occurred only in the inter-networks, including all of the FPN-associated inter-networks, the OCC-DMN inter-network and the OCC-CON inter-network. The solid black lines in Figures 2–4 mark the subcortical regions.
Figure 3. Illustration of the two-sample t-test results of the dynamic FC variation (FCV) between the senior and young groups. Pictures (A,C) depict the connections level and network averaging level, respectively; Picture (B) is a 3D view of the results shown in picture (A) (constructed using the BrainNet Viewer). Red indicates an increase in the senior group; blue indicates a decrease in the senior group. The significance level is a p-value < 0.001.
Figure 4. Illustration of the two-sample t-test results of the energy percentage in the low fluctuation frequency band compared between the senior and young groups. Pictures (A,C) depict the connections level and the network averaging level, respectively; Picture (B) is a 3D view of the results shown in picture (A) (constructed using the BrainNet Viewer). Red indicates an increase in the senior group; blue indicates a decrease in the senior group. The significance level is a p-value < 0.001.
Age-Related Changes in Frequency Spectra
The connections that covered all networks showed increased energy ratios between the energies of the higher and lower frequency bands (Figure 4A). On average, the within-network connections, which included the DMN, OCC and SMN, and the inter-network connections CON-OCC and OCC-SMN indicated significantly increased ratios in the senior group. The 3D view revealed that these changed connections had both cross-hemisphere and anterior-posterior distributions (Figure 4B).
Influence of the Sliding Window Length
The absolute difference between the mean dynamic FC and the static FC over the scan time was calculated and illustrated (Figure 5). Five network differences between the dynamic FC and the static FC are shown in the chart in Figure 5, in which the red central lines are the mean values and the gray lines indicate the values for different subjects. All five networks followed similar trends of variance, and when the window length was increased, the dynamic deviation of the static FC also varied. At approximately 50–60 s this difference reached a minimum, and at approximately 300 s the difference was at its maximum.
Figure 5. Illustration of the difference (absolute value of the difference) between the mean dynamic FC and the static FC. The results of five networks are individually averaged for all connections, which are plotted as gray curves; the red curves represent the means of all subjects from both groups.
Accordingly, the frequency spectrums of the five networks from one typical subject also varied when the sliding window length increased from 10 to 70 time points (25–175 s Figure 6). The white line in Figure 6 is at 1/60 Hz and indicates the cut-off frequency of the low-pass filtering on the dynamic FC time series. With the increasing length of the sliding window, the energy of the higher frequency attenuated faster than the energy of the lower frequency. At frequencies above 1/60 Hz, the energy was very small, and at window lengths below approximately 60 s the energy difference was sufficiently stable for examination.
Figure 6. Illustration of the frequency spectrum, which varied with the sliding window length, for all five networks of one typical subject. The results of five networks were calculated from the individually averaged FC sequence of the within-network connections.
Discussion
In this study, we showed the pattern of age-related changes in the dynamic connectivity profile between and within the whole-brain resting state networks. The originality of this study consisted of characterizing the age-related variation and frequency transition of whole-brain dynamic FC with the sliding window correlation approach. Several interesting findings were as follows: (1) in both groups, the FC variation showed distinct organization and age-related modularization, which were missing after phase randomization; (2) the FC variation indicated a significant decrease between networks, which increased with age and was dramatic in the FPN-associated and OCC-associated inter-networks; and (3) at a higher frequency, the dynamic FC showed an increased energy ratio in the senior group. These results not only shed light on the mechanisms of dynamic FC but also add to our understanding of normal brain aging. Here, we carefully discuss the results described above as well as the methodology of the sliding window correlation.
Similar to the functional specialization in FC, temporal variation of the connectivity also revealed a similar pattern or modularization, as illustrated in Figure 1. The whole-brain variations in the dynamic FC were not uniform across all systems and fell into two sides of the mean level within each group. The within-network connectivity showed a lower variation, and most of the between-network connectivity showed higher variation compared with the mean except for the CON-SMN and CON-OCC connections. Many previous studies have reported that the states of the dynamic FC matrices showed high variation in the between-network connectivity (Allen et al., 2014; Hansen et al., 2015). Allen et al. (2014) reported that seven reproducible states could be differentiated by connectivity between the DMN regions, indicating great variation. Those findings are consistent with the current results (Figures 1B,E), in which the inter-networks of the DMN showed high variation in both groups. No type of pattern among the networks was indicated in the phase randomized analysis (Figures 1A,D). This modularization of the FC variation was most likely associated with the intrinsic neural activities and interactions between the systems.
The pattern and modularization of the dynamic FC variation showed a clear decline in the senior group. The connections with higher and lower levels of variation both tended toward the mean, which indicated that less diversity in the variation between the subsystems occurred with increased age. A previous study (Du et al., 2016) reported that the FC of some connections showed less variation in schizophrenic patients. The changed modularizations of the FC variation in the two groups provided clear evidence that the modularization of the connectivity variation also declined with age, which also reflected functional integration and segmentation in aging brains (Hagmann et al., 2008; La Corte et al., 2016). Dedifferentiation of cognitive functions occurs in the aging brain, and many regions are reconfigurable to compensate for declines in other regions that occur with increased age (Sleimen-Malkoun et al., 2014). The dynamic variation in connectivity reflected that reconfiguration occurs all the time and follows some patterns, and when age increases, this pattern slowly changes. Previous studies have also reported that the brain regions dynamically participate or reconfigure into different modules during the scan time at rest (Bassett et al., 2011; Schaefer et al., 2014). The declining modularization of the FC variation in the senior group was likely related to aging and compensatory mechanisms.
The presence of FC variation, even within a single brain state (including the resting state), has been increasingly recognized (Chang and Glover, 2010; Hutchison et al., 2013a; Allen et al., 2014) and has been established as clinically relevant (Damaraju et al., 2014; Kucyi and Davis, 2014; Elton and Gao, 2015). A recent direct comparison of the awake, resting state with the anesthetized state has revealed a dramatic reduction in the connectivity variation during unconsciousness, which suggests that the connectivity variation is at least partly related to conscious operations (Barttfeld et al., 2015). One potential source of connectivity variation during consciousness is “mind wandering”, in which the brain consciously engages in different mental operations that produce fluctuations in the FC. However, we expected that the FC variation is a capacity of elasticity or operations to maintain states or connectivity transitions. This concept is similar to the cognition capacity resource, in which cognitive resources are limited and reduced in the aging brain. High variation of FC leads to the increased possibility of reconfiguration to address the loss of some resources. Qin et al. (2015) used the amplitude of the low frequency fluctuations (ALFF) of the dynamic FC to predict brain maturation between 7 years and 30 years of age. The findings of that study suggested that the increased variation was highly related to maturation. The ALLF was essentially the same with respect to the FC variation, both of which could be associated with the plasticity of the brain.
All decreased dynamic FC variations in the senior group were located in the between-network connections, especially in the FPN-linked and OCC-linked inter-networks. The inter-network between FPN and OCC was also indicated to be missing in the senior group (Figures 1B,E). These findings might indicate a posterior-anterior shift in aging (Davis et al., 2008; Vinette and Bray, 2015), which has been interpreted as compensatory in that higher-order cognitive processes are recruited to offset deficits in sensory processing. In all of the 3D views of changed connections, many posterior-anterior connections were also obvious. Cognitive and perceptual changes could be linked because they are susceptible to the same age-related factors, and a perceptual system decline could have an impact on the cognition outcome (Allen and Roberts, 2016). In young adults, visual learning engages an extended network of occipito-temporal, parietal and frontal regions, which is known to be involved in perceptual decisions (Kim and Shadlen, 1999; Shadlen and Newsome, 2001; Heekeren et al., 2004; Mayhew and Kourtzi, 2013). The subcortical regions, including the basal ganglia and thalamus, have been shown to be associated in cognition processing (Koziol and Budding, 2009). The results for the subcortical regions showed mainly increased connectivity in the DMN and decreased variation in the FPN and OCC. Considering the central or hub roles of both functional and structural networks, the subcortical regions might regulate passing signals or communication between the perception and cognitive systems, such as the OCC and FPN, the OCN and OCC, or the CON and SMN (Marchand et al., 2011). The increased fluctuation energy ratio of the high-frequency connected networks of the subcortical region might further reflect an age-related dynamic regulation, which requires further research in terms of behavioral and task fMRI experiments. With limited cognitive resources of the cognitive and perceptual systems in the aging brain, the interaction between these two systems declines as brain processing capacity is reduced.
Except for the declining capacity of cognitive resources, the changed FC fluctuation energy ratio of the higher frequencies could allow for a deeper inspection of the dynamic communication between regions. Shakil et al. (2015) found that the dynamic FC has a frequency dependance. Fluctuations within the frequency of 0–1/w were also suggested, which implied a real, physiologically dynamic connectivity (Leonardi and Van De Ville, 2015). The connectivity fluctuation is a type of intrinsic interaction of different frequency components involved in different neural activities. The speed of processing in the aging brain (Park and McDonough, 2013) is also an important property that reflects cognitive resources. Aging-related decreases in the amplitudes of low-frequency BOLD signal fluctuations have been observed, suggesting that the low-frequency fluctuations of neuroactivities are more vulnerable to aging-associated declines (Hu et al., 2014). We thought that the increased energy ratio of the fluctuating high frequencies was most likely caused by increased damage to the fluctuations of the low frequencies. The decreased variation and increased speed of fluctuations in the dynamic FC would result in disorder of the dynamic communications between different brain regions in senior individuals.
Logically, in the sliding window correlation method, the window length should have a substantial effect on the captured connectivity fluctuation. This factor was the most important consideration in terms of the overall accuracy of the technique (Leonardi and Van De Ville, 2015; Hindriks et al., 2016; Shakil et al., 2016). However, there is still no clearly determined standard for window length selection. Convergent results suggested that a window length between 50 s and 60 s is optimal. As we expected, the dynamic FC during rest was found to be a dynamic balance that maintained intrinsic connectivity patterns and even vigilance for cognitive tasks. Moreover, the variation in the FC was predicted to be centered around the static connectivity level. This expectation was supported by the appearance of a minimum on the difference curve between the mean dynamic and the static FC, as shown in Figure 5. The maximum was located at approximately 300 s (5 min), which was consistent with the suggestion of a longer scan time for reliably detecting resting state connectivity (Heekeren et al., 2004; Zuo et al., 2013). The individual frequency spectrum changed with the window length, as shown in Figure 6; this result indicated that a window shorter than 60 s retained potentially important energy for all of the networks, and the lower frequency energy contributed most of the connectivity fluctuation. A previous study (Leonardi and Van De Ville, 2015) suggested a meaningful frequency of under 1/w, which was even lower than the commonly considered frequencies for BOLD fluctuations. A fluctuation frequency of 0–1/w was suggested previously (Leonardi and Van De Ville, 2015) and employed in a later study (Qin et al., 2015; Liu et al., 2016). In the present article, the frequency was even lower than the commonly considered 0.01–0.10 Hz of the BOLD signal fluctuation, which was in accordance with the neural activity. This type of mismatch could be understood given that unlike BOLD signal fluctuations, which indicate direct neural metabolic activity, the fluctuation in the functional connectivities reflected interactions between the different regions. The interactions were also revealed in the current results, which showed that more of the affected connections were located in inter-network and inter-hemispheric connections. Furthermore, we believe that the FC variation must depended more on the fluctuation frequency, which is influenced by various factors, including anatomy, cognition, physiology and disease. The more that we understand the physiology of dynamic FC and the methodologies used to study these phenomena, the more insights we will have into the mechanisms of aging and disease.
Conclusion
In conclusion, this study presented a resting-stated dynamic FC analysis of normal brain aging. All of the results converged to expound compensation and reorganization of the networks during aging. We examined the modularization of the FC variations in the brain and decreased modularization in the aging brain. Additionally, decreased variation and increased damage to the low frequency fluctuations of the dynamic FC with aging were detected; these changes were interpreted to be associated with declining cognitive resources and limited processing speeds in the senior brain. In this article, we provided and applied new insights into FC analysis for use in aging research. The results indicated that the dynamic features of the resting-state FC were actually the intrinsic interactions between regions and cognitive resources. When some cognitive resources were reduced in aging, this type of dynamic mechanism acts to reconfigure or even train a new cognitive resource. Our conclusions in this article were fully supported by dependable results; we believe that the dynamic FC can potentially capture the intrinsic rules of compensatory processes in the aging brain, and that the present results will promote insightful understanding of spontaneous fluctuations in FC as well as aging mechanisms.
Author Contributions
MS, YL and XiongZ processed all the image data and conducted some analysis work; YC and WW were in charge of the analysis work and wrote the manuscript; JM and HN provide some useful guidance and ideas; YC, XinZ and DM designed and provided the original idea; XinZ and DM sponsored the whole research.
Funding
This research was supported by National Natural Science Foundation of China (No. 81571762, 81630051, 91520205), National Key Technology R and D Program of the Ministry of Science and Technology of China (No. 2012BAI34B02) and Tianjin Key Technology R and D Program (No. 15ZCZDSY00930).
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.
Footnotes
References
Allen, E. A., Damaraju, E., Plis, S. M., Erhardt, E. B., Eichele, T., and Calhoun, V. D. (2014). Tracking whole-brain connectivity dynamics in the resting state. Cereb. Cortex 24, 663–676. doi: 10.1093/cercor/bhs352
Allen, H. A., and Roberts, K. L. (2016). Editorial: perception and cognition: interactions in the aging brain. Front. Aging Neurosci. 8:130. doi: 10.3389/fnagi.2016.00130
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
Barttfeld, P., Uhrig, L., Sitt, J. D., Sigman, M., Jarraya, B., and Dehaene, S. (2015). Signature of consciousness in the dynamics of resting-state brain activity. Proc. Natl. Acad. Sci. U S A 112, 887–892. doi: 10.1073/pnas.1418031112
Bassett, D. S., Wymbs, N. F., Porter, M. A., Mucha, P. J., Carlson, J. M., and Grafton, S. T. (2011). Dynamic reconfiguration of human brain networks during learning. Proc. Natl. Acad. Sci. U S A 108, 7641–7646. doi: 10.1073/pnas.1018985108
Betzel, R. F., Byrge, L., He, Y., Goñi, J., Zuo, X. N., and Sporns, O. (2014). Changes in structural and functional connectivity among resting-state networks across the human lifespan. Neuroimage 102, 345–357. doi: 10.1016/j.neuroimage.2014.07.067
Calhoun, V. D., Miller, R., Pearlson, G., and Adali, T. (2014). The chronnectome: time-varying connectivity networks as the next frontier in fMRI data discovery. Neuron 84, 262–274. doi: 10.1016/j.neuron.2014.10.015
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
Chang, C., and Glover, G. H. (2010). Time-frequency dynamics of resting-state brain connectivity measured with fMRI. Neuroimage 50, 81–98. doi: 10.1016/j.neuroimage.2009.12.011
Coynel, D., Marrelec, G., Perlbarg, V., Pélégrini-Issac, M., Van de Moortele, P. F., Ugurbil, K., et al. (2010). Dynamics of motor-related functional integration during motor sequence learning. Neuroimage 49, 759–766. doi: 10.1016/j.neuroimage.2009.08.048
Damaraju, E., Allen, E. A., Belger, A., Ford, J. M., McEwen, S., Mathalon, D. H., et al. (2014). Dynamic functional connectivity analysis reveals transient states of dysconnectivity in schizophrenia. Neuroimage Clin. 5, 298–308. doi: 10.1016/j.nicl.2014.07.003
Damoiseaux, J. S., Beckmann, C. F., Arigita, E. J. S., Barkhof, F., Scheltens, P., Stam, C. J., et al. (2008). Reduced resting-state brain activity in the “default network” in normal aging. Cereb. Cortex 18, 1856–1864. doi: 10.1093/cercor/bhm207
Davis, S. W., Dennis, N. A., Daselaar, S. M., Fleck, M. S., and Cabeza, R. (2008). Que PASA? the posterior-anterior shift in aging. Cereb. Cortex 18, 1201–1209. doi: 10.1093/cercor/bhm155
De Vico Fallani, F., Pichiorri, F., Morone, G., Molinari, M., Babiloni, F., Cincotti, F., et al. (2013). Multiscale topological properties of functional brain networks during motor imagery after stroke. Neuroimage 83, 438–449. doi: 10.1016/j.neuroimage.2013.06.039
Demirta, M., Tornador, C., Falcón, C., López-Sàla, M., Hernández-Ribas, R., Pujol, J., et al. (2016). Dynamic functional connectivity reveals altered variability in functional connectivity among patients with major depressive disorder. Hum. Brain Mapp. 37, 2918–2930. doi: 10.1002/hbm.23215
Di, X., and Biswal, B. B. (2015). Characterizations of resting-state modulatory interactions in human brain. J. Neurophysiol. 114, 2785–2796. doi: 10.1152/jn.00893.2014
Dosenbach, N. U. F., Nardos, B., Cohen, A. L., Fair, D. A., Power, J. D., Church, J. A., et al. (2010). Prediction of individual brain maturity using fMRI. Science 329, 1358–1361. doi: 10.1126/science.1194144
Du, Y., Pearlson, G. D., Yu, Q., He, H., Lin, D., Sui, J., et al. (2016). Interaction among subsystems within default mode network diminished in schizophrenia patients: a dynamic connectivity approach. Schizophr. Res. 170, 55–65. doi: 10.1016/j.schres.2015.11.021
Elton, A., and Gao, W. (2015). Task-related modulation of functional connectivity variability and its behavioral correlations. Hum. Brain Mapp. 36, 3260–3272. doi: 10.1002/hbm.22847
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., McEvoy, L., Holland, D., Dale, A. M., Walhovd, K. B., and Alzheimer’s Disease Neuroimaging Initiative. (2014). What is normal in normal aging? effects of aging, amyloid and Alzheimer’s disease on the cerebral cortex and the hippocampus. Prog. Neurobiol. 117, 20–40. doi: 10.1016/j.pneurobio.2014.02.004
Friston, K. J., Holmes, A., Poline, J. B., Price, C. J., and Frith, C. D. (1996). Detecting activations in PET and fMRI: levels of inference and power. Neuroimage 4, 223–235. doi: 10.1006/nimg.1996.0074
Friston, K. J., Tononi, G., Reeke, G. N. Jr., Sporns, O., and Edelman, G. M. (1994). Value-dependent selection in the brain: simulation in a synthetic neural model. Neuroscience 59, 229–243. doi: 10.1016/0306-4522(94)90592-4
Guo, C. C., Kurth, F., Zhou, J., Mayer, E. A., Eickhoff, S. B., Kramer, J. H., et al. (2012). One-year test-retest reliability of intrinsic connectivity network fMRI in older adults. Neuroimage 61, 1471–1483. doi: 10.1016/j.neuroimage.2012.03.027
Hagmann, P., Cammoun, L., Gigandet, X., Meuli, R., Honey, C. J., Wedeen, V. J., et al. (2008). Mapping the structural core of human cerebral cortex. PLoS Biol. 6:e159. doi: 10.1371/journal.pbio.0060159
Handwerker, D. A., Roopchansingh, V., Gonzalez-Castillo, J., and Bandettini, P. A. (2012). Periodic changes in fMRI connectivity. Neuroimage 63, 1712–1719. doi: 10.1016/j.neuroimage.2012.06.078
Hansen, E. C. A., Battaglia, D., Spiegler, A., Deco, G., and Jirsa, V. K. (2015). Functional connectivity dynamics: modeling the switching behavior of the resting state. Neuroimage 105, 525–535. doi: 10.1016/j.neuroimage.2014.11.001
Hedden, T., and Gabrieli, J. D. E. (2004). Insights into the aging mind: a view from cognitive neuroscience. Nat. Rev. Neurosci. 5, 87–96. doi: 10.1038/nrn1323
Heekeren, H. R., Marrett, S., Bandettini, P. A., and Ungerleider, L. G. (2004). A general mechanism for perceptual decision-making in the human brain. Nature 431, 859–862. doi: 10.1038/nature02966
Hindriks, R., Adhikari, M. H., Murayama, Y., Ganzetti, M., Mantini, D., Logothetis, N. K., et al. (2016). Can sliding-window correlations reveal dynamic functional connectivity in resting-state fMRI? Neuroimage 127, 242–256. doi: 10.1016/j.neuroimage.2015.11.055
Hu, S., Chao, H. H. A., Zhang, S., Ide, J. S., and Li, C. S. R. (2014). Changes in cerebral morphometry and amplitude of low-frequency fluctuations of BOLD signals during healthy aging: correlation with inhibitory control. Brain Struct. Funct. 219, 983–994. doi: 10.1007/s00429-013-0548-0
Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., et al. (2013a). Dynamic functional connectivity: promise, issues and interpretations. Neuroimage 80, 360–378. doi: 10.1016/j.neuroimage.2013.05.079
Hutchison, R. M., Womelsdorf, T., Gati, J. S., Everling, S., and Menon, R. S. (2013b). Resting-state networks show dynamic functional connectivity in awake humans and anesthetized macaques. Hum. Brain Mapp. 34, 2154–2177. doi: 10.1002/hbm.22058
Kern, A., and Behl, C. (2009). The unsolved relationship of brain aging and late-onset Alzheimer disease. Biochim. Biophys. Acta 1790, 1124–1132. doi: 10.1016/j.bbagen.2009.07.016
Kim, J.-N., and Shadlen, M. N. (1999). Neural correlates of a decision in the dorsolateral prefrontal cortex of the macaque. Nat. Neurosci. 2, 176–185. doi: 10.1038/5739
Koziol, L. F., and Budding, D. E. (2009). Subcortical Structures and Cognition: Implications for Neuropsychological Assessment. New York, NY: Springer.
Kucyi, A., and Davis, K. D. (2014). Dynamic functional connectivity of the default mode network tracks daydreaming. Neuroimage 100, 471–480. doi: 10.1016/j.neuroimage.2014.06.044
La Corte, V., Sperduti, M., Malherbe, C., Vialatte, F., Lion, S., Gallarda, T., et al. (2016). Cognitive decline and reorganization of functional connectivity in healthy aging: the pivotal role of the salience network in the prediction of age and cognitive performances. Front. Aging Neurosci. 8:204. doi: 10.3389/fnagi.2016.00204
Leonardi, N., and Van De Ville, D. (2015). On spurious and real fluctuations of dynamic functional connectivity during rest. Neuroimage 104, 430–436. doi: 10.1016/j.neuroimage.2014.09.007
Li, X., Zhu, D., Jiang, X., Jin, C., Zhang, X., Guo, L., et al. (2014). Dynamic functional connectomics signatures for characterization and differentiation of PTSD patients. Hum. Brain Mapp. 35, 1761–1778. doi: 10.1002/hbm.22290
Liu, F., Wang, Y., Li, M., Wang, W., Li, R., Zhang, Z., et al. (2016). Dynamic functional network connectivity in idiopathic generalized epilepsy with generalized tonic-clonic seizure. Hum. Brain Mapp. 38, 957–973. doi: 10.1002/hbm.23430
Marchand, W. R., Lee, J. N., Suchy, Y., Garn, C., Johnson, S., Wood, N., et al. (2011). Age-related changes of the functional architecture of the cortico-basal ganglia circuitry during motor task execution. Neuroimage 55, 194–203. doi: 10.1016/j.neuroimage.2010.12.030
Mayhew, S. D., and Kourtzi, Z. (2013). Dissociable circuits for visual shape learning in the young and aging human brain. Front. Hum. Neurosci. 7:75. doi: 10.3389/fnhum.2013.00075
Meindl, T., Teipel, S., Elmouden, R., Mueller, S., Koch, W., Dietrich, O., et al. (2010). Test-retest reproducibility of the default-mode network in healthy individuals. Hum. Brain Mapp. 31, 237–246. doi: 10.1002/hbm.20860
Meunier, D., Stamatakis, E. A., and Tyler, L. K. (2014). Age-related functional reorganization, structural changes and preserved cognition. Neurobiol. Aging 35, 42–54. doi: 10.1016/j.neurobiolaging.2013.07.003
Mosher, K. I., and Wyss-Coray, T. (2014). Microglial dysfunction in brain aging and Alzheimer’s disease. Biochem. Pharmacol. 88, 594–604. doi: 10.1016/j.bcp.2014.01.008
Ou, J., Xie, L., Jin, C., Li, X., Zhu, D., Jiang, R., et al. (2015). Characterizing and differentiating brain state dynamics via hidden markov models. Brain Topogr. 28, 666–679. doi: 10.1007/s10548-014-0406-2
Park, D. C., and McDonough, I. M. (2013). The dynamic aging mind: revelations from functional neuroimaging research. Perspect. Psychol. Sci. 8, 62–67. doi: 10.1177/1745691612469034
Qin, J., Chen, S.-G., Hu, D., Zeng, L.-L., Fan, Y.-M., Chen, X.-P., et al. (2015). Predicting individual brain maturity using dynamic functional connectivity. Front. Hum. Neurosci. 9:418. doi: 10.3389/fnhum.2015.00418
Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U S A 98, 676–682. doi: 10.1073/pnas.98.2.676
Rashid, B., Damaraju, E., Pearlson, G. D., and Calhoun, V. D. (2014). Dynamic connectivity states estimated from resting fMRI identify differences among schizophrenia, bipolar disorder and healthy control subjects. Front. Hum. Neurosci. 8:897. doi: 10.3389/fnhum.2014.00897
Reeve, A., Simcox, E., and Turnbull, D. (2014). Aging and Parkinson’s disease: why is advancing age the biggest risk factor? Aging Res. Rev. 14, 19–30. doi: 10.1016/j.arr.2014.01.004
Sakolu, U., Pearlson, G. D., Kiehl, K. A., Wang, Y. M., Michael, A. M., and Calhoun, V. D. (2010). A method for evaluating dynamic functional network connectivity and task-modulation: application to schizophrenia. MAGMA 23, 351–366. doi: 10.1007/s10334-010-0197-8
Sala-Llonch, R., Bartrés-Faz, D., and Junqué, C. (2015). Reorganization of brain networks in aging: a review of functional connectivity studies. Front. Psychol. 6:663. doi: 10.3389/fpsyg.2015.00663
Sambataro, F., Murty, V. P., Callicott, J. H., Tan, H. Y., Das, S., Weinberger, D. R., et al. (2010). Age-related alterations in default mode network: impact on working memory performance. Neurobiol. Aging 31, 839–852. doi: 10.1016/j.neurobiolaging.2008.05.022
Schaefer, A., Margulies, D. S., Lohmann, G., Gorgolewski, K. J., Smallwood, J., Kiebel, S. J., et al. (2014). Dynamic network participation of functional connectivity hubs assessed by resting-state fMRI. Front. Hum. Neurosci. 8:195. doi: 10.3389/fnhum.2014.00195
Shadlen, M. N., and Newsome, W. T. (2001). Neural basis of a perceptual decision in the parietal cortex (area LIP) of the rhesus monkey. J. Neurophysiol. 86, 1916–1936.
Shakil, S., Keilholz, S. D., and Lee, C. H. (2015). “On frequency dependencies of sliding window correlation”, in Proceeding IEEE International Conference on Bioinformatics and Biomedicine (BIBM), (Washington, DC), 363–368. doi: 10.1109/BIBM.2015.7359708
Shakil, S., Lee, C. H., and Keilholz, S. D. (2016). Evaluation of sliding window correlation performance for characterizing dynamic functional connectivity and brain states. Neuroimage 133, 111–128. doi: 10.1016/j.neuroimage.2016.02.074
Shehzad, Z., Kelly, A. M. C., Reiss, P. T., Gee, D. G., Gotimer, K., Uddin, L. Q., et al. (2009). The resting brain: unconstrained yet reliable. Cereb. Cortex 19, 2209–2229. doi: 10.1093/cercor/bhn256
Sleimen-Malkoun, R., Temprado, J.-J., and Hong, S. L. (2014). Aging induced loss of complexity and dedifferentiation: consequences for coordination dynamics within and between brain, muscular and behavioral levels. Front. Aging Neurosci. 6:140. doi: 10.3389/fnagi.2014.00140
Sugiura, M. (2016). Functional neuroimaging of normal aging: declining brain, adapting brain. Aging Res. Rev. 30, 61–72. doi: 10.1016/j.arr.2016.02.006
Vinette, S. A., and Bray, S. (2015). Variation in functional connectivity along anterior-to-posterior intraparietal sulcus and relationship with age across late childhood and adolescence. Dev. Cogn. Neurosci. 13, 32–42. doi: 10.1016/j.dcn.2015.04.004
Wang, L., Laviolette, P., O’Keefe, K., Putcha, D., Bakkour, A., Van Dijk, K. R. A., et al. (2010). Intrinsic connectivity between the hippocampus and posteromedial cortex predicts memory performance in cognitively intact older individuals. Neuroimage 51, 910–917. doi: 10.1016/j.neuroimage.2010.02.046
Wee, C.-Y., Yang, S., Yap, P.-T., Shen, D., and Alzheimer’s Disease Neuroimaging Initiative (2016). Sparse temporally dynamic resting-state functional connectivity networks for early MCI identification. Brain Imaging Behav. 10, 342–356. doi: 10.1007/s11682-015-9408-2
Xu, Y., Yan, J., Zhou, P., Li, J., Gao, H., Xia, Y., et al. (2012). Neurotransmitter receptors and cognitive dysfunction in Alzheimer’s disease and Parkinson’s disease. Prog. Neurobiol. 97, 1–13. doi: 10.1016/j.pneurobio.2012.02.002
Yi, L. Y., Liang, X., Liu, D. M., Sun, B., Ying, S., Yang, D. B., et al. (2015). Disrupted topological organization of resting-state functional brain network in subcortical vascular mild cognitive impairment. CNS Neurosci. Ther. 21, 846–854. doi: 10.1111/cns.12424
Yu, Q., Erhardt, E. B., Sui, J., Du, Y., He, H., Hjelm, D., et al. (2015). Assessing dynamic brain graphs of time-varying connectivity in fMRI data: application to healthy controls and patients with schizophrenia. Neuroimage 107, 345–355. doi: 10.1016/j.neuroimage.2014.12.020
Zalesky, A., and Breakspear, M. (2015). Towards a statistical test for functional connectivity dynamics. Neuroimage 114, 466–470. doi: 10.1016/j.neuroimage.2015.03.047
Zou, L., Abutalebi, J., Zinszer, B., Yan, X., Shu, H., Peng, D., et al. (2012). Second language experience modulates functional brain network for the native language production in bimodal bilinguals. Neuroimage 62, 1367–1375. doi: 10.1016/j.neuroimage.2012.05.062
Zuo, X. N., Di Martino, A., Kelly, C., Shehzad, Z. E., Gee, D. G., Klein, D. F., et al. (2010). The oscillating brain: complex and reliable. Neuroimage 49, 1432–1445. doi: 10.1016/j.neuroimage.2009.09.037
Keywords: aging, dynamic functional connectivity, functional connectivity variation, functional connectivity fluctuation frequency, resting-stated fMRI
Citation: Chen Y, Wang W, Zhao X, Sha M, Liu Y, Zhang X, Ma J, Ni H and Ming D (2017) Age-Related Decline in the Variation of Dynamic Functional Connectivity: A Resting State Analysis. Front. Aging Neurosci. 9:203. doi: 10.3389/fnagi.2017.00203
Received: 01 December 2016; Accepted: 06 June 2017;
Published: 30 June 2017.
Edited by:
Ana B. Vivas, CITY College, International Faculty of the University of Sheffield, GreeceReviewed by:
Ana Maria Buga, University of Medicine and Pharmacy of Craiova, RomaniaCharis Styliadis, Aristotle University of Thessaloniki, Greece
Copyright © 2017 Chen, Wang, Zhao, Sha, Liu, Zhang, Ma, Ni and Ming. 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) or licensor 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: Xin Zhao, zhaoxin@tju.edu.cn
Dong Ming, richardming@tju.edu.cn