An Effective Method to Identify Adolescent Generalized Anxiety Disorder by Temporal Features of Dynamic Functional Connectivity

Introduction

Adolescent generalized anxiety disorder (GAD) is one of the fairly common anxiety disorders among youth. Typical symptoms of adolescent GAD include excessive, uncontrolled and lasting anxiety over common things and events in daily life. Although high risk of adult GAD, social phobia and major depressive disorder (MDD) exists in adolescents with GAD (Pine et al., 1998), adolescent GAD has been rarely investigated in neurophysiological aspect, and by now the diagnosis of adolescent GAD has mainly depended on clinical symptoms and signs. Researches aiming at revealing the pathophysiology and finding relatively objective biomarkers of adolescent GAD are therefore of great importance at current stage.

Recently, resting-state functional magnetic resonance imaging (fMRI) has been a widely-used and valid technique to reveal the functional abnormalities in brain caused by various mental diseases. Functional connectivity (FC), defined as the temporal dependence of neuronal activity patterns of anatomically separated brain regions (Aertsen et al., 1989), is able to describe the functional communication between brain regions during task-involved or resting state. Many studies have investigated FC to explore the abnormal functional organization of brain in patients with specific psychiatric illnesses, such as major depression disorder, social anxiety disorder (SAD), schizophrenia and Alzheimer’s disease (AD; Greicius et al., 2007; Stam et al., 2007; Hahn et al., 2011). By present, task-involved studies have reported the disease-related FC in GAD patients, and most of them were associated with amygdala, prefrontal cortex (PFC), anterior cingulate cortex (ACC) and concentrated in the structures of “fear circuit” which deals with fear and other negative emotions (Monk et al., 2008; Etkin et al., 2010). Resting state FC was mainly calculated on the basis of predefined regions of interest (ROI). Amygdala was the ROI widely investigated in GAD studies, and abnormal FC between amygdala and PFC has been frequently reported in these researches (Etkin et al., 2009; Liu W. J. et al., 2015). Disrupted amygdalar subregions FC, increased FC between amygdala and cerebellum, insula, superior temporal gyrus and putamen in GAD patients have also been found (Etkin et al., 2009; Liu W. J. et al., 2015). Additionally, by the whole brain FC analysis, increased FC between hippocampus and fusiform in GAD patients has been revealed by Cui et al. (2016). These studies tried to estimate the FC that was able to describe the functional relevance between brain regions during scanning, which was based on the assumption that the FC remained stationary (static FC) during rest or task-processed state. However, the varying levels of attention, mind-wandering, even mood-swings may take place during scanning, leaving the observed blood oxygenation level-dependent (BOLD) signals non-stationary and deviating from the assumption. Moreover, the static FC abandons dynamic properties of FC which might be thought to be related to specific diseases (Sakoglu et al., 2010; Jones et al., 2012) and potential to serve as biomarkers.

Consequently, increasing number of articles has been published to explore the temporal properties of FC. Different from the model of static FC, dynamic FC is evaluated as time-varying covariance of neural signals between brain regions, making it able to describe the collaboration of brain regions in a precise way. Temporal characteristics of dynamic FC are therefore promising to explore the temporal alterations related to psychiatric illnesses. Several previous studies have investigated the dynamic FC in mental diseases (Damaraju et al., 2014; Rashid et al., 2014; Yu et al., 2015). Damaraju et al. (2014) found that schizophrenia patients failed to maintain states typified by strong, large-scale connectivity, and the abnormal connections in schizophrenia patients during the states deviated more from normal levels. By investigating dynamic FC, Rashid et al. (2014) observed fewer transitions to specific states in schizophrenia and bipolar disorder patients. The findings mentioned above revealed the transient abnormal connectivity patterns which might be absent in the researches of static FC. Additionally, variance in the dynamic graph metrics of brain were introduced by Yu et al. (2015), and the result showed the less changeable states transition in schizophrenia patients. In summary, the dynamic analysis of FC appears to be useful for acquiring additional measurements to investigate the alterations caused by mental diseases. To the best of our knowledge, the temporal features of dynamic FC in GAD patients have not been explored, and the alterations caused by GAD have remained unclear. In the present study, Hurst exponent (H) and variance of dynamic FC were calculated to explore the sensitivities to GAD and possibilities of biomarkers. Additionally, previous classification studies for GAD have focused on clinical scales and questionnaires, meaning that our work would be the pilot work on identifying adolescent GAD by neuroimaging data.

Both healthy adolescents and those suffering from GAD were included in our study. A data-driven method independent component analysis (ICA) was utilized to parcellate cortex into intrinsic connectivity networks (ICNs). To classify adolescents with GAD from healthy controls, heterogeneous features including H, variance of dynamic FC and static FC were calculated. In order to find out the impact of adolescent GAD on the connections, the activation pattern analysis proposed in the previous study (Haufe et al., 2014) was applied to the classification model, and the regions associated with the connections with lager impact of GAD were considered GAD-related. We hypothesized that the temporal features of dynamic FC could be helpful to classification and potential to serve as biomarkers of GAD.

Materials and Methods

Subjects

In our study, 31 adolescents diagnosed with GAD and 28 demographically similar normal controls (NC) were recruited from local high schools in Hunan province by advertisement or school notice. We had fully explained the study to each adolescent and his or her legal guardians before the written informed consent was obtained from them. The adolescents were diagnosed by the same trained clinician using the fourth version of Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) criteria and the Schedule for Affective Disorders and Schizophrenia for School Age Children-Present and Lifetime (K-SADS-PL) version. All of subjects were right-handed, non-medicated, and voluntary to the study with the permission of their legal guardians. After excluding the subjects with excessive head motion (>2 mm translation or >2 degree rotation), the remaining 27 patients and 25 NC kept gender-, age-, IQ-matched and did not differ in head motion according to framewise displacement (FD) measurement (Jenkinson et al., 2002; see Table 1). Every recruited patient was diagnosed with current first episode GAD without comorbidity disorder by clinical psychiatrists. Exclusion criteria for all subjects were the same as previous studies (Zhang et al., 2013; Liao et al., 2014), including seizures history, neurological abnormalities, head trauma or unconsciousness, physical disease and use of psychoactive substances. This study was carried out in accordance with the recommendations of ICH-GCP, “China-GCP”, related regulation and law of China and the Ethics Committee at the Second Xiangya Hospital of Central South University 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 Ethics Committee at the Second Xiangya Hospital of Central South University.

TABLE 1

Table 1. Demographical characteristics of the generalized anxiety disorder (GAD) patients and normal controls (NCs).

Image Acquisition

Subjects were instructed to keep their eyes closed, relaxed, awake and mindless during scanning. Foam paddings and earplugs were used to limit head motion. After scanning, subjects were asked whether they had fallen asleep or opened their eyes during scanning, and those responded yes were excluded. All of the participants’ high-resolution T1-weighted anatomic images and resting-state fMRI were performed on a 3.0-Tesla Philips MRI scanner. The anatomic images were acquired by magnetization-prepared rapid gradient-echo sequence (repetition time (TR) = 7.5 ms, echo time (TE) = 3.7 ms, slice thickness = 1 mm, filed of view = 256 × 256 mm², flip angle = 8°, number of slices = 180). Resting-state functional images were obtained using an echo-planar imaging (EPI) sequence and the parameters included TR = 3000 ms, TE = 30 ms, filed of view = 240 × 240 mm², flip angle = 90°, number of slices = 36, slice thickness = 4 mm, in-plane matrix = 64 × 64 and total volume = 180.

Image Preprocessing

All participants’ images were preprocessed with the Statistical Parametric Mapping (SPM8) based toolkit Data Processing Assistant Resting-State fMRI (DPARSFA; Chao-Gan and Yu-Feng, 2010)¹. The first 10 volumes for each subject were discarded to allow for equilibration of the magnetic field. Slices signal acquisition time correction, head motion correction, and realignment were applied to the remaining volumes. To reduce the potential influence of head motion, subjects with head motion exceeding 2 mm or 2 degree were excluded, which leaved 25 healthy adolescents and 27 GAD patients for the further analysis. All images were then spatially normalized to Montreal Neurological Institute (MNI) space with 3 mm × 3 mm × 3 mm resolution. Resulting images were spatially smoothed with a Gaussian kernel of 8 mm full width at half maximum (FWHM).

Group ICA and Post Processing

Group ICA

Instead of using a priori template, we utilized the data-driven technique group ICA to avoid the deviation which might be induced by pre-defined templates. All preprocessed fMRI data were entered to perform a group ICA implemented in GIFT toolbox². To allow investigation of FC between subsystems of ICNs (Allen et al., 2014), a relatively high-order ICA with 100 independent components (ICs) was executed. In the data dimension reduction step, principle components analysis (PCA) was performed to obtain 150 principal components for subject-specific data and 100 principal components for group data maintained by expectation-maximization (EM) algorithm (Roweis, 1998). And then Infomax group ICA algorithm (Calhoun et al., 2001) was used to retain 100 ICs. Spatial maps and corresponding time serials for every subject were acquired after the back reconstruction step. Consistent with previous studies (Menon, 2011; Rashid et al., 2014), eight ICNs consisting of 44 ICs were selected from the 100 ICs. Shown in Figure 1, the eight selected ICNs corresponded to auditory (AUD), basal ganglia (BG), sensorimotor (SM), posterior insula (PINS), visual (VIS), cognitive control (CC), cerebellar (CB), and default mode (DM) networks.

FIGURE 1

Figure 1. Spatial distributions of non-artificial intrinsic connectivity networks (ICNs). Spatial distribution of each ICN was in a white frame, and spatial distributions of the same components were in the row within the black box.

Post Processing

Additional steps were carried out to further eliminate artificial noise of the time serials in ICNs. The linear, quadratic and cubic trends of time serials were removed, and the six rigid realignment parameters and their temporal derivatives were regressed out. Outliers of time serials were detected and removed by 3DDESPIKE. A low-pass filter with the cutoff frequency of 0.15 Hz was used to remove the high frequency components of signals (Allen et al., 2014).

Dynamic FC and Temporal Features Estimation

Dynamic FC Calculation

Dynamic connection between any pair of ICs was estimated by instantaneous phase difference using Hilbert transform. The previous study (Glerean et al., 2012) showed that this method could be reliable and comparable with correlation-based methods but with higher time resolution. Recently, this method has been used to detect community structures in brain, and the resulting synchronization communities were similar to those acquired by well-established approach such as ICA (Ponce-Alvarez et al., 2015). The computational model used in the recent study (Demirtas et al., 2016) was adopted in our work. Every post processed signal of selected IC was Hilbert transformed to get the corresponding analytic signal with instantaneous phases. Then the instantaneous phase difference between each pair of ICs was calculated and normalized between 0 and 1 indicating perfect synchronization and perfect anti-synchronization, respectively (Ponce-Alvarez et al., 2015).

Temporal Features Estimation

Vaillancourt and Newell (2002) argued that abnormal complexity of a physiological or behavioral system might be the result of aging and disease. Temporal complexity, defined as the difficulties arising when describing or predicting a signal (Lu et al., 2008), has been applied to fMRI researches. Hurst exponent widely used to reflect the “long-term memory” of time series, is directly related to the fractal dimension of signals and enable to measure the temporal complexity of physiological signal. Previous study (Maxim et al., 2005) has revealed the disease-related temporal complexity in several regions of cortex by comparing the H of BOLD signals between healthy controls and AD patients. One purpose of the present study is to figure out the extent to which GAD affects the H of dynamic FC. R/S method (Annis and Lloyd, 1976) was adopted to calculate H in our study. The calculation procedure is described below:

1. Dynamic FC is divided into d subseries of length n;

2. For each subseries X_m, m = 1…d, X_m is normalized to $Z_m=X_m-\overline{X_m},$ where $\overline{X_m}$ is the mean value of subseries X_m;

3. Calculate cumulative time series $Y_{j,m}={\displaystyle {\sum }_{k=1}^j{Z}_{k,m}},$ j = 1…n for every normalized subseries, where Z_k,m is the k-th value of Z_m;

4. Calculate range R_m = max(Y_1,m, Y_2,m, … Y_n,m) − min(Y_1,m, Y_2,m, … Y_n,m);

5. Rescale the range R_m/S_m, where S_m is the standard deviation of subseries X_m;

6. Calculate the mean value over all rescaled ranges of length n: ${(R/S)}_n=\frac{1}{d}{\displaystyle {\sum }_{m=1}^d{R}_m/S_m}$.

Presented in the previous study (Mandelbrot and Wallis, 1969), the relation between length n and corresponding rescaled range (R/S)_n is: (R/S)_n~cn^H. Therefore, the estimation of H can be acquired by linear regression over a sample of increasing time horizons: log(R/S)_n = logc + Hlogn. Due to the relatively small number of time points of dynamic FC, The (R/S)_n estimated in our study is corrected with the method mentioned in the previous study (Weron, 2002).

The other temporal feature of dynamic FC, variance independent to H, is investigated in the present study as well. The variance of dynamic FC reflects the temporal variability in FC across time. Larger variances indicate more changeable functional connections between brain regions, and are related to the frequently changing states of information processing in brain (Yu et al., 2015; Marusak et al., 2017). The most used unbiased estimator of variance was used in our study. In addition, we calculated static FC of all pairs of ICs as Pearson correlation coefficients of the corresponding signals. Fisher’s Z transformed static FC was then used as features to classify GAD patients from healthy controls to test the validity of being a biomarker.

Feature Selection

To pick out the disease-related connections and acquire better classification performance, nested feature selection procedure was conducted to gather the features with high discriminative power before training a classifier. H, variance of dynamic FC and static FC were regarded as three types of features to be selected for classifying healthy controls and GAD patients. First, to compare the effectiveness of every type of features, each type of features was selected and used alone to identify adolescent GAD. Subsequently, the union set of the temporal features including H and variance and that of all the three types were treated as features respectively to find out the highest accuracy. Of note, H and variance of dynamic FC were normalized by Z score method, and static FC was Fisher’s Z transformed before feature selection.

Not being limited by the uncertain distribution of data, measuring the relevance to class labels is a robust method to evaluate the discriminative power of features. We utilized the Kendall (1948) tau rank correlation coefficient to quantify the relevance to classification of features, with larger absolute values indicating greater discriminative powers. In the present study, the features were sorted by the absolute value of Kendall tau rank correlation coefficient in descending order, and the top features with coefficients over a threshold were selected to be the input of the classifier. In order to find out the best number of features for classification, the number of features to be selected in every type increased from 1 to 400 with the step length of 1, and the feature set with the highest accuracy was adopted.

Since the feature set is consisted of multiple types of features, we selected top features from every type of features respectively, and the union set of all selected features from each type were input into classifier.

Support Vector Machine Learning

The support vector machine (SVM) was applied to classification in our study. SVM is a supervised machine learning algorithm introduced by Vapnik (2013). This algorithm is pretty suitable for the pattern cognition problem with small sample number, thus in theory being able to handle the classification with relatively small number of samples in the present study. To reduce the risk of overfitting, the linear kernel SVM was chosen in our work. We adopted the LIBLINEAR toolbox (Fan et al., 2008) on MATLAB to construct classifier, and the default parameters were used in the current study.

Leave-One-Out Cross-Validation and Activation Pattern Analysis

Since relatively few samples were included in our study, we used the leave-one-out cross-validation (LOOCV) strategy to estimate the generalization ability of the classifiers. In every LOOCV iteration, one sample was chosen as testing set and the rest samples were used for training set. The feature selection and classifier training steps were only carried out with the data in training set, without using the information of the testing set. The selected features of testing sample were then entered to the trained classifier to predict whether the testing subject was GAD patient or healthy control. This procedure was repeated 52 iterations to ensure every subject was predicted.

In order to find out the altered features caused by adolescent GAD, activation pattern analysis was performed in our study. Shown in previous study (Haufe et al., 2014), the weights of multivariate linear classifiers could not directly interpret the relationship between features and the sample labels, because significant nonzero weights might also be obtained at features which are statistically independent of the class labels, and the transformation converting weights vector to activation patterns could achieve the desirable interpretation. Proven in Haufe et al. (2014), the activation pattern of a linear classifier could be obtained by: A = cov(X) × W × cov(S)⁻¹, where X is the matrix consisted of selected features in training set, S is the vector of the sample labels in training set, W is the weights of selected features, and cov() represents covariance of a given variable. The values in A indicate the effect directions and strengths of the class label in the selected features, which enable the desirable interpretation. In the present study, the activation pattern analysis was conducted when H and variance were combined. In each LOOCV iteration, the activation value of every selected feature was calculated, and we set the values of the features not selected to zero. Thus, by averaging the absolute values of activation patterns across all iterations, we obtained every feature’s standard activation value which was able to quantify the GAD’s influence on the feature. We defined the activation value of every IC by summing the standard activation values of features associated with the IC. The significantly important ICs were defined as regions whose activation values were not less than the sum of two standard deviations and mean value of all ICs’ activation values. Additionally, the features selected in all iterations of LOOCV were regarded stably helpful to classification, and were defined as consensus connections.

Results

Classification Accuracy

Table 2 showed the accuracies, sensitivities, specificities and the area under the curve (AUC) obtained from receiver operating characteristic (ROC) analysis of all cases of classifications with different feature sets. The highest accuracy of 88.46% and the largest AUC of 0.8889 were acquired when both H and variance of dynamic FC were used. In this case, the top 211 ranked variance and the top 138 ranked H were selected in each LOOCV iteration. It should be noted that when H, variance of dynamic FC and static FC were all used, the highest accuracy could be 88.46% as well, but the number of selected features in static FC could hardly affect this accuracy once there were around 211 and 138 features selected in variance and H respectively. Therefore, we considered the best accuracy was achieved when only H and variance used.

TABLE 2

Table 2. Classification performances of different types of features.

Both H and variance of dynamic FC performed relatively better than static FC when used as features alone. When only H of dynamic FC was used, the highest accuracy was 84.62% with the corresponding AUC of 0.8726, and the top 57 ranked H were selected in each LOOCV iteration. When it came to variance of dynamic FC, the highest accuracy of 82.69% with the corresponding AUC of 0.8800 were achieved, and the top 62 ranked variance were selected in each LOOCV iteration. Relatively low accuracy of 61.54% was acquired when only static FC was used to classification. This difference in accuracy might indicate that temporal characteristics of dynamic FC were affected by GAD and suitable to serve as biomarkers of GAD.

Figure 2 showed the ROC analysis of the four classifiers with different types of features used. The trend of AUC was a little different as that of accuracy, which meant the largest AUC of 0.8889 was obtained when temporal features were combined, the second was 0.8800 when only variance was used, the third was 0.8726 when only H was used, and the smallest was 0.5511 when only static FC was used.

FIGURE 2

Figure 2. Receiver operating characteristic (ROC) curves of different types of features. Different colors were used to represent the ROC curves of all four types of features.

Activation Pattern Analysis

In the case of best accuracy where the top 211 ranked variance and the top 138 ranked H were selected in each LOOCV iteration, the activation pattern analysis was conducted. Since the implication of H differs from that of variance, the important ICs were identified respectively from the standard activation values of H and those of variance. The distribution of consensus features and their standard activation values were shown in Figure 3.

FIGURE 3

Figure 3. Distribution of consensus connections. Nodes were colored by ICN, and their sizes were weighted by activation values. The red lines indicate the decreased trend of value in generalized anxiety disorder (GAD) adolescents, and the blue lines indicate the increased trend of value in GAD adolescents. The widths of lines are weighted by their standard activation values. (A) Region activation values and consensus connections of H features. (B) Region activation values and consensus connections of variance features.

The important IC identified from standard activation values of H features was DMN2. CB1 and DMN2 were identified as important ICs associated with variance features. DMN2 was shared by the two kinds of important ICs, which suggested that the neural activity in this network might be greatly influenced by GAD. The spatial distributions of the two kinds of important ICs were shown in Table 3.

TABLE 3

Table 3. Spatial distributions of important independent components (ICs).

Discussions

In the present work, we managed to classify adolescents with GAD from healthy controls using information from resting FC. To our knowledge, studies exploring the GAD-related FC alterations have been based on static FC, and there has been no study investigating the temporal features of dynamic FC in GAD patients. Results in our study showed that H and variance of dynamic FC could be treated as features to identify adolescents with GAD even when used alone (see Table 2). Furthermore, we acquired a relatively high accuracy of 88.46% when H and variance of dynamic FC were combined, which suggested that they were potential to be used as diagnostic criteria for GAD. Moreover, by activation pattern analysis, we identified specific brain regions affected by GAD: DMN2 for H features; CB1 and DMN2 for variance features. DMN2 was the region important to both temporal features.

Classification Performance

When used alone as features, both H and variance of dynamic FC achieved a higher accuracy (larger than 61.54%) compared with static FC, and the performance of H was slightly better than variance. It’s not surprising that the combination of H and variance improved the performance of classification, which implies that these two temporal characteristic of dynamic FC are complementary to each other, and it is necessary to take them into account to identify GAD patients. Of note, when static FC was the only used type of feature to discriminate subjects, accuracy decreased to 61.54% that is much lower than that achieved by temporal features of dynamic FC. The further analysis showed that adding features from static FC failed to promote performance of classifier. The reason might be that static FC abandoned the dynamic information of FC, and was not precise enough to describe the communication between brain regions, which might veil the disease-related details of brain activity.

ICNs with Greater Magnitude of Activation Value

In the present study, by activation pattern analysis, DMN was found greatly affected by GAD. DMN has been considered as a critical role in monitoring the external environment and supporting internal mentation (Mitchell et al., 2006; Gilbert et al., 2007; Zhong et al., 2014). Recently, there have been several studies revealing the abnormal activities of DMN in patients with anxiety disorder, including less deactivation in medial prefrontal cortex (MPFC) and greater deactivation in posterior cingulate cortex (PCC; Zhao et al., 2007), decreased functioning (Sylvester et al., 2012), abnormal connectivity between posterior hippocampal and DMN (Chen and Etkin, 2013). In the current study, we identified the GAD-affected region, DMN2 by both the activation pattern analyses of H and variance features. DMN2 is mainly located in MPFC which is the core region of DMN. MPFC is known to be associated with self-referential processing (Kelley et al., 2002; Northoff et al., 2006) and has the ability to understand the mental state of oneself and others (Allen et al., 2008). The abnormal FC dynamics associated with MPFC revealed the disrupted information exchange related to this region, and might result in the failure to sense how “self and others think, feel, perceive, imagine, react, attribute, infer and so on.” (Sharp et al., 2011), which could make the patients with GAD unable to control their worry and anxiety as sensitively as healthy controls. Our results might provide additional evidence for the notion that DMN plays a key role in the pathophysiology of GAD.

Besides the IC in DMN network, cerebellum was also characterized as GAD-related regions by activation pattern analysis. It’s not surprising that cerebellum is included in the GAD-affected regions for there are evidence proving the role cerebellum plays in certain non-motor functions, including emotion and cognitive processing regulation (Dolan, 1998; Schmahmann and Sherman, 1998; Schmahmann and Caplan, 2006; Hu et al., 2008). Previous study showed the impairment to cerebellum existed in patients with GADs (Abadie et al., 1999). Also, several studies (Critchley et al., 2000; Bonne et al., 2003; Sakai et al., 2005) on anxiety symptoms reported a hyperactivity of the cerebellum. Since the hyperactivity of cerebellum was found in patients with MDD as well, a recent study (Phillips et al., 2015) suggested that the abnormality in cerebellum might be the reason of the attention impairments of the two disorders. The abnormal dynamics associated with cerebellum found in our study is consistent with the previous findings and reflect the impaired information exchange between cerebellum and other regions, which might cause the extreme autonomic reactions coupled with anxiety (American Psychiatric Association, 2013) and the impaired cognitive control functions (Bögels and Zigterman, 2000; Stefanopoulou et al., 2014).

Comparison with Other Classification Analyses

To the best of our knowledge, no previous study has explored the framework of classification of GAD. Therefore, we compared the proposed method in our study with three other methods for SAD classification. The sample size of all the four studies are comparable, and the leave-one-out cross validation was adopted in all of these studies. Pantazatos et al. (2014) used the static FC during face processing task as feature to discriminate individuals with SAD from healthy controls. Liu F. et al. (2015) calculated the static FC between each pair of regions based on automated anatomical labeling atlas, and resulting FC were used as feature to identify patients with SAD. Zhang et al. (2015) adopted the regional homogeneity of brain voxels as feature for SAD classification. As we can see in Table 4, the performance of the method proposed in the present study is among the best, which proves the efficacy of our proposed method.

TABLE 4

Table 4. Comparison on classification performance of social anxiety disorder (SAD) methods.

Limitations

In our study, temporal features of dynamic FC were used successfully to identify adolescent GAD at a relatively high accuracy. Although results showed that both H and variance of dynamic FC were able to serve as features to classify GAD adolescents from healthy controls, there were limitations in the present study. The first was the relatively smaller sample (52 subjects in total) of subjects included in our work. With a small number of samples included in the present study, leave-one-out cross validation was adopted to prevent the training set from deviating too much from the overall population. This process is helpful to discover the critical features derived from dynamic FC, but it also increases the risk of overfitting. Second, the selected ICs were not complete and accurate enough to parcellate the whole brain, resulting in the loss of information and details of brain activity. Given these limitations, larger number of subjects and more precise atlas of brain are needed in the future study.

Conclusion

To sum up, relatively higher accuracy of differentiating adolescents with GAD from healthy controls was obtained when temporal properties of dynamic FC were combined, proving the potential of dynamic FC as a biomarker for adolescent GAD. Moreover, by activation pattern analysis, we discovered the GAD-related regions, including the subregions of CB and DMN. Furthermore, our results suggested that MPFC was the prominent region with more sensitivity in identification of adolescent GAD.

Author Contributions

ZY, ML, BH and LL conceived and designed the experiments. ZY, TH and YZ organized and analyzed the raw data. ZY, TH, ML and LL participated in the statistical analysis and interpretation of data. ZY and TH wrote the article. ZZ, BH, FZ, JG and DM revised the manuscript.

Funding

This work was supported by the National Basic Research Program of China (973 Program) (No. 2014CB744600), the National Natural Science Foundation of China (Grant No. 61632014, No. 61210010), Program of Beijing Municipal Science & Technology Commission (No. Z171100000117005), and the Program of International S&T Cooperation of MOST (No. 2013DFA11140).

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.

The reviewer BD and handling Editor declared their shared affiliation.

Footnotes

1. ^ http://www.restfmri.net
2. ^ http://icatb.sourceforge.net/

References

Abadie, P., Boulenger, J. P., Benali, K., Barre, L., Zarifian, E., and Baron, J. C. (1999). Relationships between trait and state anxiety and the central benzodiazepine receptor: a PET study. Eur. J. Neurosci. 11, 1470–1478. doi: 10.1046/j.1460-9568.1999.00556.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Aertsen, A. M., Gerstein, G. L., Habib, M. K., and Palm, G. (1989). Dynamics of neuronal firing correlation: modulation of “effective connectivity”. J. Neurophysiol. 61, 900–917.

PubMed Abstract | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen, J. G., Fonagy, P., and Bateman, A. W. (2008). Mentalizing in Clinical Practice. Washington, DC: American Psychiatric Publishing.

Google Scholar

Annis, A., and Lloyd, E. (1976). The expected value of the adjusted rescaled Hurst range of independent normal summands. Biometrika 63, 111–116. doi: 10.1093/biomet/63.1.111

CrossRef Full Text | Google Scholar

American Psychiatric Association. (2013). Diagnostic and Statistical Manual of Mental Disorders (DSM-5^®). 5th Edn. Arlington, VA: American Psychiatric Publishing.

Google Scholar

Bögels, S. M., and Zigterman, D. (2000). Dysfunctional cognitions in children with social phobia, separation anxiety disorder, and generalized anxiety disorder. J. Abnorm. Child Psychol. 28, 205–211. doi: 10.1023/A:1005179032470

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonne, O., Gilboa, A., Louzoun, Y., Brandes, D., Yona, I., Lester, H., et al. (2003). Resting regional cerebral perfusion in recent posttraumatic stress disorder. Biol. Psychiatry 54, 1077–1086. doi: 10.1016/s0006-3223(03)00525-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Calhoun, V. D., Adali, T., Pearlson, G. D., and Pekar, J. (2001). A method for making group inferences from functional MRI data using independent component analysis. Hum. Brain Mapp. 14, 140–151. doi: 10.1002/hbm.1048

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao-Gan, Y., and Yu-Feng, Z. (2010). DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, A. C., and Etkin, A. (2013). Hippocampal network connectivity and activation differentiates post-traumatic stress disorder from generalized anxiety disorder. Neuropsychopharmacology 38, 1889–1898. doi: 10.1038/npp.2013.122

PubMed Abstract | CrossRef Full Text | Google Scholar

Critchley, H. D., Corfield, D. R., Chandler, M. P., Mathias, C. J., and Dolan, R. J. (2000). Cerebral correlates of autonomic cardiovascular arousal: a functional neuroimaging investigation in humans. J. Physiol. 523, 259–270. doi: 10.1111/j.1469-7793.2000.t01-1-00259.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, H., Zhang, J., Liu, Y., Li, Q., Li, H., Zhang, L., et al. (2016). Differential alterations of resting-state functional connectivity in generalized anxiety disorder and panic disorder. Hum. Brain Mapp. 37, 1459–1473. doi: 10.1002/hbm.23113

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Demirtas, M., Tornador, C., Falcón, C., López-Solà, 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Dolan, R. J. (1998). A cognitive affective role for the cerebellum. Brain 121, 545–546. doi: 10.1093/brain/121.4.545

PubMed Abstract | CrossRef Full Text | Google Scholar

Etkin, A., Prater, K. E., Hoeft, F., Menon, V., and Schatzberg, A. F. (2010). Failure of anterior cingulate activation and connectivity with the amygdala during implicit regulation of emotional processing in generalized anxiety disorder. Am. J. Psychiatry 167, 545–554. doi: 10.1176/appi.ajp.2009.09070931

PubMed Abstract | CrossRef Full Text | Google Scholar

Etkin, A., Prater, K. E., Schatzberg, A. F., Menon, V., and Greicius, M. D. (2009). Disrupted amygdalar subregion functional connectivity and evidence of a compensatory network in generalized anxiety disorder. Arch. Gen. Psychiatry 66, 1361–1372. doi: 10.1001/archgenpsychiatry.2009.104

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, R.-E., Chang, K.-W., Hsieh, C.-J., Wang, X.-R., and Lin, C.-J. (2008). LIBLINEAR: a library for large linear classification. J. Mach. Learn. Res. 9, 1871–1874. doi: 10.1145/1390681.1442794

CrossRef Full Text | Google Scholar

Gilbert, S. J., Dumontheil, I., Simons, J. S., Frith, C. D., and Burgess, P. W. (2007). Comment on “Wandering minds: the default network and stimulus-independent thought”. Science 317:43, author reply 43. doi: 10.1126/science.1140801

PubMed Abstract | CrossRef Full Text | Google Scholar

Glerean, E., Salmi, J., Lahnakoski, J. M., Jaaskelainen, I. P., and Sams, M. (2012). Functional magnetic resonance imaging phase synchronization as a measure of dynamic functional connectivity. Brain Connect. 2, 91–101. doi: 10.1089/brain.2011.0068

PubMed Abstract | CrossRef Full Text | Google Scholar

Greicius, M. D., Flores, B. H., Menon, V., Glover, G. H., Solvason, H. B., Kenna, H., et al. (2007). Resting-state functional connectivity in major depression: abnormally increased contributions from subgenual cingulate cortex and thalamus. Biol. Psychiatry 62, 429–437. doi: 10.1016/j.biopsych.2006.09.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Hahn, A., Stein, P., Windischberger, C., Weissenbacher, A., Spindelegger, C., Moser, E., et al. (2011). Reduced resting-state functional connectivity between amygdala and orbitofrontal cortex in social anxiety disorder. Neuroimage 56, 881–889. doi: 10.1016/j.neuroimage.2011.02.064

PubMed Abstract | CrossRef Full Text | Google Scholar

Haufe, S., Meinecke, F., Görgen, K., Dahne, S., Haynes, J. D., Blankertz, B., et al. (2014). On the interpretation of weight vectors of linear models in multivariate neuroimaging. Neuroimage 87, 96–110. doi: 10.1016/j.neuroimage.2013.10.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, D., Shen, H., and Zhou, Z. (2008). Functional asymmetry in the cerebellum: a brief review. Cerebellum 7, 304–313. doi: 10.1007/s12311-008-0031-2

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

Jones, D. T., Vemuri, P., Murphy, M. C., Gunter, J. L., Senjem, M. L., Machulda, M. M., et al. (2012). Non-stationarity in the “resting brain’s” modular architecture. PLoS One 7:e39731. doi: 10.1371/journal.pone.0039731

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelley, W. M., Macrae, C. N., Wyland, C. L., Caglar, S., Inati, S., and Heatherton, T. F. (2002). Finding the self? An event-related fMRI study. J. Cogn. Neurosci. 14, 785–794. doi: 10.1162/08989290260138672

PubMed Abstract | CrossRef Full Text | Google Scholar

Kendall, M. G. (1948). Rank Correlation Methods. London: Griffin.

Google Scholar

Liao, M., Yang, F., Zhang, Y., He, Z., Su, L., and Li, L. (2014). White matter abnormalities in adolescents with generalized anxiety disorder: a diffusion tensor imaging study. BMC Psychiatry 14:41. doi: 10.1186/1471-244X-14-41

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, F., Guo, W., Fouche, J. P., Wang, Y., Wang, W., Ding, J., et al. (2015). Multivariate classification of social anxiety disorder using whole brain functional connectivity. Brain Struct. Funct. 220, 101–115. doi: 10.1007/s00429-013-0641-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W. J., Yin, D. Z., Cheng, W. H., Fan, M. X., You, M. N., Men, W. W., et al. (2015). Abnormal functional connectivity of the amygdala-based network in resting-state FMRI in adolescents with generalized anxiety disorder. Med. Sci. Monit. 21, 459–467. doi: 10.12659/MSM.893373

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, S., Chen, X., Kanters, J. K., Solomon, I. C., and Chon, K. H. (2008). Automatic selection of the threshold value R for approximate entropy. IEEE Trans. Biomed. Eng. 55, 1966–1972. doi: 10.1109/TBME.2008.919870

PubMed Abstract | CrossRef Full Text | Google Scholar

Mandelbrot, B. B., and Wallis, J. R. (1969). Robustness of the rescaled range R/S in the measurement of noncyclic long run statistical dependence. Water Resour. Res. 5, 967–988. doi: 10.1029/wr005i005p00967

CrossRef Full Text | Google Scholar

Marusak, H. A., Calhoun, V. D., Brown, S., Crespo, L. M., Sala-Hamrick, K., Gotlib, I. H., et al. (2017). Dynamic functional connectivity of neurocognitive networks in children. Hum. Brain Mapp. 38, 97–108. doi: 10.1002/hbm.23346

PubMed Abstract | CrossRef Full Text | Google Scholar

Maxim, V., Sendur, L., Fadili, J., Suckling, J., Gould, R., Howard, R., et al. (2005). Fractional Gaussian noise, functional MRI and Alzheimer’s disease. Neuroimage 25, 141–158. doi: 10.1016/j.neuroimage.2004.10.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Menon, V. (2011). Large-scale brain networks and psychopathology: a unifying triple network model. Trends Cogn. Sci. 15, 483–506. doi: 10.1016/j.tics.2011.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, J. P., Macrae, C. N., and Banaji, M. R. (2006). Dissociable medial prefrontal contributions to judgments of similar and dissimilar others. Neuron 50, 655–663. doi: 10.1016/j.neuron.2006.03.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Monk, C. S., Telzer, E. H., Mogg, K., Bradley, B. P., Mai, X., Louro, H. M., et al. (2008). Amygdala and ventrolateral prefrontal cortex activation to masked angry faces in children and adolescents with generalized anxiety disorder. Arch. Gen. Psychiatry 65, 568–576. doi: 10.1001/archpsyc.65.5.568

PubMed Abstract | CrossRef Full Text | Google Scholar

Northoff, G., Heinzel, A., de Greck, M., Bermpohl, F., Dobrowolny, H., and Panksepp, J. (2006). Self-referential processing in our brain—a meta-analysis of imaging studies on the self. Neuroimage 31, 440–457. doi: 10.1016/j.neuroimage.2005.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pantazatos, S. P., Talati, A., Schneier, F. R., and Hirsch, J. (2014). Reduced anterior temporal and hippocampal functional connectivity during face processing discriminates individuals with social anxiety disorder from healthy controls and panic disorder and increases following treatment. Neuropsychopharmacology 39, 425–434. doi: 10.1038/npp.2013.211

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, J. R., Hewedi, D. H., Eissa, A. M., and Moustafa, A. A. (2015). The cerebellum and psychiatric disorders. Front. Public Health 3:66. doi: 10.3389/fpubh.2015.00066

PubMed Abstract | CrossRef Full Text | Google Scholar

Pine, D. S., Cohen, P., Gurley, D., Brook, J., and Ma, Y. (1998). The risk for early-adulthood anxiety and depressive disorders in adolescents with anxiety and depressive disorders. Arch. Gen. Psychiatry 55, 56–64. doi: 10.1001/archpsyc.55.1.56

PubMed Abstract | CrossRef Full Text | Google Scholar

Ponce-Alvarez, A., Deco, G., Hagmann, P., Romani, G. L., Mantini, D., and Corbetta, M. (2015). Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneity. PLoS Comput. Biol. 11:e1004100. doi: 10.1371/journal.pcbi.1004100

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Roweis, S. T. (1998). “EM algorithms for PCA and SPCA,” in Advances in Neural Information Processing Systems, eds M. I. Jordan, M. J. Kearns and S. A. Solla (Cambridge, MA: MIT Press), 626–632.

Google Scholar

Sakai, Y., Kumano, H., Nishikawa, M., Sakano, Y., Kaiya, H., Imabayashi, E., et al. (2005). Cerebral glucose metabolism associated with a fear network in panic disorder. Neuroreport 16, 927–931. doi: 10.1097/00001756-200506210-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakoglu, 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

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmahmann, J. D., and Caplan, D. (2006). Cognition, emotion and the cerebellum. Brain 129, 290–292. doi: 10.1093/brain/awh729

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmahmann, J. D., and Sherman, J. C. (1998). The cerebellar cognitive affective syndrome. Brain 121, 561–579. doi: 10.1093/brain/121.4.561

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharp, C., Pane, H., Ha, C., Venta, A., Patel, A. B., Sturek, J., et al. (2011). Theory of mind and emotion regulation difficulties in adolescents with borderline traits. J. Am. Acad. Child Adolesc. Psychiatry 50, 563.e1–573.e1. doi: 10.1016/j.jaac.2011.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Stam, C. J., Jones, B. F., Nolte, G., Breakspear, M., and Scheltens, P. (2007). Small-world networks and functional connectivity in Alzheimer’s disease. Cereb. Cortex 17, 92–99. doi: 10.1093/cercor/bhj127

PubMed Abstract | CrossRef Full Text | Google Scholar

Stefanopoulou, E., Hirsch, C. R., Hayes, S., Adlam, A., and Coker, S. (2014). Are attentional control resources reduced by worry in generalized anxiety disorder? J. Abnorm. Psychol. 123, 330–335. doi: 10.1037/a0036343

PubMed Abstract | CrossRef Full Text | Google Scholar

Sylvester, C. M., Corbetta, M., Raichle, M. E., Rodebaugh, T. L., Schlaggar, B. L., Sheline, Y. I., et al. (2012). Functional network dysfunction in anxiety and anxiety disorders. Trends Neurosci. 35, 527–535. doi: 10.1016/j.tins.2012.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Vaillancourt, D. E., and Newell, K. M. (2002). Changing complexity in human behavior and physiology through aging and disease. Neurobiol. Aging 23, 1–11. doi: 10.1016/s0197-4580(01)00247-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Vapnik, V. (2013). The Nature of Statistical Learning Theory. New York, NY: Springer Science and Business Media.

Google Scholar

Weron, R. (2002). Estimating long-range dependence: finite sample properties and confidence intervals. Phys. A Stat. Mechan. Appl. 312, 285–299. doi: 10.1016/s0378-4371(02)00961-5

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Li, L., Yu, R., Liu, J., Tang, J., Tan, L., et al. (2013). White matter integrity alterations in first episode, treatment-naive generalized anxiety disorder. J. Affect. Disord. 148, 196–201. doi: 10.1016/j.jad.2012.11.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, W., Yang, X., Lui, S., Meng, Y., Yao, L., Xiao, Y., et al. (2015). Diagnostic prediction for social anxiety disorder via multivariate pattern analysis of the regional homogeneity. Biomed. Res. Int. 2015:763965. doi: 10.1155/2015/763965

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X. H., Wang, P. J., Li, C. B., Hu, Z. H., Xi, Q., Wu, W. Y., et al. (2007). Altered default mode network activity in patient with anxiety disorders: an fMRI study. Eur. J. Radiol. 63, 373–378. doi: 10.1016/j.ejrad.2007.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, Y., Huang, L., Cai, S., Zhang, Y., von Deneen, K. M., Ren, A., et al. (2014). Altered effective connectivity patterns of the default mode network in Alzheimer’s disease: an fMRI study. Neurosci. Lett. 578, 171–175. doi: 10.1016/j.neulet.2014.06.043

PubMed Abstract | CrossRef Full Text | Google Scholar