2413

Optimal Time-Dependent Window-Size Reveals a More Accurate Picture of Dynamic Functional Connectivity
Xiaowei Zhuang1, Zhengshi Yang1, Brent Bluett1, Sarah Banks1, and Dietmar Cordes1,2

1Cleveland Clinic Lou Ruvo Center for Brain Health, Las Vegas, NV, United States, 2University of Colorado, Boulder, CO, United States

Synopsis

We have introduced a new method to determine the optimal time-dependent window-size for calculating sliding-window correlations between two non-stationary time series. The time-dependent window-size is calculated from the local information of intrinsic mode functions of each time series computed using empirical mode decomposition. Results from simulation demonstrate that the running-correlation computed with a time-dependent window-size is able to capture local transients without creating unstable fluctuations. By incorporating the optimal window-size in a whole-brain dynamic functional connectivity analysis, we are able to view differences in whole-brain temporal dynamics between normal control subjects and PD subjects more precisely.

Introduction

Temporal dynamics of brain’s intrinsic networks have been widely studied using a sliding-window method1, 2. A technical challenge of this method is the proper choice of the window-size. Specifically, a large window may miss existing temporal transients whereas a small window may produce unstable results3. In this study, we describe a new method to determine a time-dependent window-size in a sliding-window analysis. This window-size is calculated from the instantaneous period and energy of each intrinsic mode function (IMF) obtained from empirical mode decomposition (EMD) of each time series4, 5. The IMFs track the local periodic changes of non-stationary time series and can be used to compute an optimum window-size at each time-point. We further incorporate the optimum window-size to explore whole-brain dynamic functional connectivity2,6 changes in normal controls (NCs) and subjects with Parkinson’s disease (PD).

Methods

Time-dependent window-size: Given two time series y1(t) and y2(t) where t denotes each time point, a time-dependent window-size is determined as follows: Each time course (yp(t),p=1,2) is first decomposed into different IMFs (cpi(t)) using EMD4,5, i.e.yp(t)=Σicpi(t)+rp(t), where rp(t) are the residuals. The Hilbert transform4 is then performed on each IMF to compute the instantaneous period (Tpi(t)) and energy density (Epi(t)). Next, the time-dependent period (Tp(t)) at every t is determined as an average of Tpi(t) weighted by the energy density Epi(t) , i.e.Tp(t)=(1/ΣiEpi(t))*Σi(Tpi(t)*Epi(t)). The final time-dependent window-size of y1(t) and y2(t) to obtain an optimal sliding-window correlation is chosen to be Td(t)=max(T1(t),T2(t)). Simulation: Two non-stationary time series, y1 and y2, were simulated at TR = 1s (Fig. 1(A)), with a static correlation of -0.02. The dynamic correlations between y1 and y2 were calculated and compared using the sliding-window method with two fixed window-sizes: 10TR and 50TR, as well as the time-dependent window-size Td(t). Real fMRI data: We studied data from 18 NCs (14 Males; 64.25±9.50 years) and 20 de-novo PD patients (11 Males; 58.03±11.54 years; UPDRS-III: 15.05±7.43) obtained from Parkinson’s Progression Markers Initiative (PPMI) database7. Resting-state fMRI were performed on 3T Siemens scanners (TR/TE/Flip Angle/Resolution=2.4s/25ms/80deg/3.3mm3, 210 time-frames) and went through standard preprocessing steps. Whole-brain dynamic functional connectivity analysis with the time-dependent window-size: A flow chart of our analysis is presented in Fig. 2. First, a group ICA with 100 components was carried out using fMRI time-series from both PD and NC subjects. 72 components were identified as network-related components and the corresponding subject-specific ICA maps and time courses were calculated using dual-regression. The whole-brain dynamic connectivity matrices were calculated using sliding-window correlations between pairwise subject-specific ICA time series with the time-dependent window-size Td(t). The connectivity matrices in each sliding-window were then concatenated in time for each subject, and further stacked for both PD and NC subjects. K-means clustering was performed on the concatenated connectivity matrix from all subjects to estimate dynamic functional states for both groups. The optimum cluster number k was determined by the leave-one-out cross-validation. Finally, the time spent in each state and the frequency of state alternations were calculated for every subject separately and used to compare the temporal dynamics between the PD and NC groups.

Results

The dynamic correlations between the two simulated time series (Fig.1(A)) were calculated using different window-sizes (Fig.1(B)). Compared to the fixed window-size (dashed lines, Fig.1(C)), the time-dependent window (solid line) captures local transients and avoids unstable fluctuations in calculating correlations. Fig. 2 shows the flow chart of our whole-brain dynamic connectivity analysis with time-dependent window-size. Static correlation within and between each network (Fig.3(A)) are revealed by Fisher’s z statistics and shown in Fig. 3(B). Three dynamic functional states are determined from the leave-one-out cross-validation in k-means clustering for both PD and NC subjects (Fig. 4(A)). Fig. 4(B) shows that NC subjects spend significantly more time in State II, which has stronger connections, both between and within networks, whereas PD patients tend to stay longer in weaker connected functional states I and III. Furthermore, significant reduced frequency of state alternations (p=0.0046) is found in PD group (Fig. 4(D)), which is not observed when repeat the same analysis with a fixed window-size of 13TR (~30secs) or 20TR (~50secs), as used in previous studies2,3.

Discussion and Conclusion

Our method determined a time-dependent window-size in a sliding-window method based on local information of IMFs from each ICA time series. The running-correlation computed with a time-dependent window-size was able to capture local transients without creating unstable fluctuations (Fig. 1). By incorporating the optimum window-size in whole-brain dynamic functional connectivity analysis (Fig. 2), we found altered whole-brain temporal dynamics in PD subjects (Fig. 4). This finding corroborates with previously reported electrophysiology data in PD8.

Acknowledgements

The study is supported by the National Institutes of Health (grant number 1R01EB014284 and P20GM109025). PPMI is sponsored and partially funded by The Michael J. Fox Foundation for Parkinson’s Research (MJFF). Other funding partners include a consortium of industry players, non-profit organizations and private individuals (for a full list see http://www.ppmi-info.org/about-ppmi/who-we-are/study-sponsors/).

References

[1]. Chang and Glover. 2010. Time–frequency dynamics of resting-state brain connectivity measured with fMRI. NeuroImage.

[2]. Allen et al. 2014. Tracking whole-brain connectivity dynamics in the resting state. Cerebral.

[3]. Hutchison et al. 2013. Dynamic functional connectivity: promise, issues, and interpretations. NeuroImage.

[4]. Huang et al. 1998. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. In Proceedings of the Royal Society of London A: mathematical, physical and engineering sciences. The Royal Society.

[5]. Chen et al. 2010. The time-dependent intrinsic correlation based on the empirical mode decomposition. Advances in Adaptive Data Analysis.

[6]. Yu et al. 2015. Assessing dynamic brain graphs of time-varying connectivity in fMRI data: application to healthy controls and patients with schizophrenia. NeuroImage.

[7]. Parkinson Progression Marker Initiative. 2011. The Parkinson Progression Marker Initiative (PPMI). Prog Neurobiol.

[8]. Yanagisawa et al. 2012. Regulation of motor representation by phase–amplitude coupling in the sensorimotor cortex. The Journal of Neuroscience.

Figures

Figure 1. Simulation. (A) Two non-stationary simulated time series: y1 (blue) and y2 (red). (B). Instantaneous period of y1 (dashed blue) and y2 (dashed red), and the time-dependent window-size (green) at each time point. (C). The dynamic correlation between two synthetic time series was calculated using sliding-window method with different window-sizes: 10TR (dashed blue), 50TR (dashed red) and time-dependent window-size (solid green). The static correlation between these two time series are indicated by the dashed purple line.

Figure 2. Whole-brain dynamic functional connectivity analysis flow chart.

Figure 3. (A). Spatial components from 72 network-related ICA components. (B). Whole brain static functional connectivity matrix.

Figure 4. (A). Three dynamic states obtained from whole-brain dynamic functional connectivity analysis using a time-dependent window-size. (B). Average window fractions of time spent in each state (window in each state / total number of windows) for PD (orange) and NC (blue) groups. (C). Statistical comparisons of window fractions spent in each state between PD and NC groups. * denotes p<0.05. (D). Statistical comparison of frequency of state alternations between PD and NC groups.

Proc. Intl. Soc. Mag. Reson. Med. 26 (2018)
2413