Next Article in Journal
Spike-Dependent Dynamic Partitioning of the Locus Coeruleus Network through Noradrenergic Volume Release in a Simulation of the Nucleus Core
Next Article in Special Issue
Characteristics of Families with Members Going Out for Work and Associated Factors among Persons with Schizophrenia in Rural China
Previous Article in Journal
Integrating Distribution-Based and Anchor-Based Techniques to Identify Minimal Important Change for the Tinnitus Functional Index (TFI) Questionnaire
Previous Article in Special Issue
The Impact of Electroconvulsive Therapy on Negative Symptoms in Schizophrenia and Their Association with Clinical Outcomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Frequency-Specific Analysis of the Dynamic Reconfiguration of the Brain in Patients with Schizophrenia

College of Information and Computer, Taiyuan University of Technology, No. 209, Daxue Street, Yuci District, Jinzhong 030024, China
*
Author to whom correspondence should be addressed.
Brain Sci. 2022, 12(6), 727; https://doi.org/10.3390/brainsci12060727
Submission received: 24 March 2022 / Revised: 1 May 2022 / Accepted: 28 May 2022 / Published: 1 June 2022

Abstract

:
The analysis of resting-state fMRI signals usually focuses on the low-frequency range/band (0.01–0.1 Hz), which does not cover all aspects of brain activity. Studies have shown that distinct frequency bands can capture unique fluctuations in brain activity, with high-frequency signals (>0.1 Hz) providing valuable information for the diagnosis of schizophrenia. We hypothesized that it is meaningful to study the dynamic reconfiguration of schizophrenia through different frequencies. Therefore, this study used resting-state functional magnetic resonance (RS-fMRI) data from 42 schizophrenia and 40 normal controls to investigate dynamic network reconfiguration in multiple frequency bands (0.01–0.25 Hz, 0.01–0.027 Hz, 0.027–0.073 Hz, 0.073–0.198 Hz, 0.198–0.25 Hz). Based on the time-varying dynamic network constructed for each frequency band, we compared the dynamic reconfiguration of schizophrenia and normal controls by calculating the recruitment and integration. The experimental results showed that the differences between schizophrenia and normal controls are observed in the full frequency, which is more significant in slow3. In addition, as visual network, attention network, and default mode network differ a lot from each other, they can show a high degree of connectivity, which indicates that the functional network of schizophrenia is affected by the abnormal brain state in these areas. These shreds of evidence provide a new perspective and promote the current understanding of the characteristics of dynamic brain networks in schizophrenia.

1. Introduction

As a complex psychotic disorder, schizophrenia (SZ) usually occurs in adolescence or adulthood and is related to abnormal integration between distal brain regions [1]. It is considered a chronic and diverse genetic disease with abnormal perception, emotion, and brain connections [2]. Clinically, it is emotional and cognitive dysfunction, accompanied by symptoms such as hallucinations and delusions [1,3,4]. Therefore, SZ is increasingly considered a disease caused by brain network dysfunction.
In functional magnetic resonance imaging (fMRI) studies, traditional brain network construction methods are generally during the resting-state. Current dynamic network analyses have confirmed that fluctuations in functional connections exist, which has attracted increasing attention in the academic world [5,6]. Studies have shown that dynamic network analysis can better detect fluctuations in the brain’s functional connections. For example, Gifford and colleagues used a modeling method called the dynamic modular organization to investigate better inter-group differences in dynamic community structure in SZ [7]. Cui and colleagues also believe that dynamic functional connectivity can reflect the time-varying characteristics of brain networks and capture the changes of network topology and cognitive behavior over time [8]. The most common method of dynamic network research is to use sliding window technology to divide blood oxygen level-dependent (BOLD) signals into shorter time intervals or windows and obtain functional brain networks from each time interval [9,10]. These intervals provide more reliable results for the dynamic analysis of SZ.
Dynamic reconfiguration of multilayer brain networks can be effective to analyze SZ caused by cognitive impairment and is increasingly used to find abnormal neural activities. Previous studies have revealed that low frequency is generally associated with abnormal dynamic functional connections of the brain network. To illustrate, Debo Dong and colleagues [11] comprehensively studied the dynamic reconfiguration of the dynamic reorganization of resting-state functional connectivity (rsFC) in SZ when the brain frequency was less than 0.1 Hz and found patients who showed consistently increased rsFC variability in sensory and perceptual systems. Xiao Wang and colleagues believed there are significant differences between slow4 and slow5’s SZ functional connectivity densities [12]. In conclusion, studies on the dynamic reconfiguration of multilayer brain networks will provide evidence for further exploration of the mechanism of schizophrenia.
Previous studies on SZ’s multilayer dynamic brain network generally focused on low frequency [12,13] because this frequency band is usually associated with the abnormal dynamic functional connection of the brain network. Still, it cannot completely cover the complex neural activities in the resting-state brain. Some research studies have shown that brain regions in different frequency bands have different resting-state functional connections. In other words, physiological signals in the same brain network may compete or cooperate with each other in various frequency bands [1]. Recently, researchers have found that the high-frequency band can also provide valuable information for the diagnosis of SZ. For example, Gifford and colleagues found that the flexibility of a satisfactory score is higher at 0.08–0.25 Hz, which supports the view that SZ involves dynamic brain network changes [7]. Yu et al. and colleagues showed that considering different high-frequency bands helps measure brain activity in SZ [14]. That is to say, frequency is closely related to the study of SZ and the abnormal structure of large-scale brain networks.
Based on the above research, we assumed that the brain network connection of SZ is frequency dependent. Therefore, we used the following five different frequency bands to explore the changes between brain regions (full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz)), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz) and slow2 (0.198–0.25 Hz). A multilayer community detection algorithm is applied to identify the temporal community and calculate recruitment, integration, and module allegiance [15]. These two indicators reflect the dynamic interactions between and within different brain regions and quantify community changes over time [8] (Figure 1). In conclusion, our purpose was to study the dynamic reconfigurability of the brain network structure of SZ in different frequency bands and provide new ideas for future research.

2. Materials and Methods

2.1. Participants

The data were selected from the UCLA Consortium for Neuropsychiatric Phenomics LA5c Study project, and the UCLA Institutional Review Board approved the study. The data were obtained via a public database, open fMRI (https://openfmri.org/dataset/ds000030/, accesssed on 1 September 2020). Our study comprised 82 subjects, including 42 schizophrenia patients (SZ) and 40 normal controls (NC). There was no significant difference in age and gender between the NC and SZ (T-test). It is worth noting that 42 SZ in 82 subjects participated in the evaluation of SAPS. Specific demographic characteristics are shown in Table 1.

2.2. Imaging Acquisition and Preprocessing

All subjects underwent structural and functional MRI scans on a 3-T scanner with a 32 channel head coil at the China University of Electronic Science and Technology. Soft foam and earplugs were used to fix the subjects’ head and reduce scanning noise during the scanning process. They were asked to stay relaxed, open their eyes, and avoid intentional thinking activities. Resting-state fMRI data were collected using a sequence of T2*-weighted echo-planar imaging (EPI). The parameters are as follows: repetition time (TR) = 2 s, echo time (TE) = 30 ms, slice thickness = 4 mm, slices = 34, flip angle = 90°, field of view (FOV) = 192 mm, and matrix = 64 × 64. The resting-state fMRI scan lasted a total of 304 s. For the structural scan, a T1-weighted high-resolution anatomical scan was collected with the following parameters: slice thickness = 1 mm, slices = 176, repetition time (TR) = 1.9 s, echo time(TE) = 2.26 ms, matrix = 256 × 256, and field of view (FOV) = 250 mm.
Data preprocessing was carried out by the DPABI toolbox (http://rfmri.org/dpabi, accessed on 1 September), the first ten volumes of the signal were discarded, and the data of the other 142 time points were preprocessed as follows: (1) slice timing correction; (2) realignment; (3) normalization: the image space was standardized to the Montreal Neurological Institute (MNI) head anatomy template and resampled with 3 × 3 × 3 mm3 voxels [16]. (4) Filtering: considering the frequency-specific, five frequency bands are divided here: full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz) [17]; (5) Smoothing: the images were spatially smoothed using a Gaussian filter with a full-width at half-maximum (FWHM) of 6 mm [18]; (6) the covariates were removed, and the brain was divided into 90 regions using the automatic anatomical marker template (AAL) [19] and residual time-series were extracted for each voxel.
According to the time series extracted by the AAL template, we divided brain regions into five functional networks [20]: somatosensory/motor and auditory network (SMN), visual network (VN), attention network (AN), default mode network (DMN) and limbic/paralimbic and subcortical network (LSN).

2.3. Formatting of Mathematical Components

After the above processing, each subject was divided into five frequency bands. The sliding time window is divided for each frequency band, and the time series is divided into smaller time intervals. According to Leonardi and Van De Ville [21], the minimum window size that can be used in dynamic network studies is 1 / f m i n ( f m i n denotes the minimum frequency of the data included). Here, we select the minimum value of the whole frequency band (0.01 Hz) as the frequency for dividing the time window. Hence, the window length we used was 100 s (50 TRs). Each window was shifted 2 s (1 TR), resulting in 93 overlapping windows.
To better consider the changes in network organization over time, we introduced a multilayer network model [22]. We connected nodes in one time window to themselves in adjacent time windows to represent time dependence. The supra-adjacency matrix of an f-layered multilayer network can be expressed as:
A = [ A 1 H 1 f H f 1 A f ]
where A α is the intra-layer network adjacency matrix of layer α , 1 α f . H k l is a diagonal matrix and a square matrix corresponding to the inter-layer connection matrix between layers k and l .

2.4. Multilayer Community Detection

The community detection algorithm provides a method called “module”, which can decompose the network into dense node groups. Here, we use a multilayer community detection algorithm called Genlouvain to identify the brain community of each participant [23,24]. It divides the communities in the multilayer network by optimizing the multilayer modular quality function Q and getting the information about the nodes. It can be defined as:
Q = 1 2 μ i j l r { ( A i j l γ l B i j l ) δ l r + δ i j ω j l r } δ ( g i l , g j r )
where μ is the total edge weight of the network and δ ( g i l , g j r ) is 1 if nodes belong to the same community and 0 otherwise [25]. A i j l represents the adjacency matrix between node i and node j in layer l , γ l is the structural resolution parameter of layer l and element B i j l represents the component of layer l matrix corresponding to the optimized zero model. The element ω j l r gives the connection strength from node j of layer r to node j of layer l , which is called “inter-layer coupling parameter”. Here, we set the structure resolution parameter and the inter-layer coupling parameter to the default value of 1 [8,26,27].
We used the Genlouvain MATLAB toolbox [23,24] to calculate the community allocations. This community assignment represents the evolution of assumed functional modules in the brain network over time. In addition, the algorithm is a generalized greedy-like algorithm [28], whose output may change at each run due to the random nature of the optimized partition function. Therefore, we repeat the detection algorithm 50 times, and then take the average of these measures as the final estimate value.

2.5. Dynamic Network Statistics

2.5.1. Module Allegiance

In the module allegiance matrix P , each element P i j represents the probability that node i and node j are assigned to the same community during the whole scanning process [29]. The matrix P consists of N × N and is the number of brain regions (In this study, N = 90). P i j can be written as follows:
P i j = 1 O T o = 1 O t = 1 T a i , j k , o
where O is the number of iterations of the multilayer community detection algorithm, T is the number of layers. For each optimization O and layers T , if they are in the same community network, the value of module loyalty is 1 (the values on the main diagonal of the matrix are all 1); otherwise, it is 0.

2.5.2. Recruitment and Integration

Recruitment: We use the module allegiance matrix to evaluate the dynamic role of cognitive systems in task execution [15]. According to the above module allegiance matrix, recruitment and integration are calculated. Recruitment is defined as the probability that its region and the region from the same system appear in the community together. For node i in the community network S , the recruitment is specifically expressed as follows:
R i S = 1 n s j S P i j
where n S is the number of nodes in network S . P i j represents the number of times that nodes i and j are assigned to the same module.
Integration: Integration is defined as the probability that its region will co-appear in the community with regions from other systems. For node i in the community network S , the integration can be written as follows:
I i S = 1 N n S j S P i j
where N is the total number of brain regions.

2.6. Statistical Analysis

An independent sample T-test was used for statistical tests to quantify the differences in integration and recruitment between NC and SZ. For each metric, we averaged the metrics obtained from the network. Spearman’s correlation was used to assess the correlation between SZ and SAPS scale scores. The Benjamini and Hochberg error discovery rate (BH_FDR) method was used to calibrate all the results, and the threshold value for the significant difference was set to <0.05.

3. Results

3.1. Group Comparisons of the Whole-Brain Level at Different Frequencies

In this part, we average the recruitment and integration of the whole brain region of each subject and discuss the dynamic difference between brain networks in full frequency and slow5, slow4, slow3 and slow2. See Table 2 for the specific number of modules As shown in Figure 2A, there was no significant difference in the integration between the two groups, but the recruitment decreased significantly in full frequency (T(80) = −2.600, P = 0.011) and slow3 (T(80) = −3.090, P = 0.003) (Figure 2B).

3.2. Group Comparisons of RSN Level at Different Frequencies

According to the AAL template, we divide 90 brains into five functional networks in the preprocessing. Here, we will study the differences in functional networks of SZ in five frequency bands. As shown in Figure 3A, by comparing the integration between NC and SZ, we found no significant difference in the 0.01–0.25 Hz and slow2, slow4, slow5. However, in slow3, the VN results showed a significant increase (T(80) = 2.263, P = 0.026). As shown in Figure 3B, the recruitment results showed a significant difference in the 0.01–0.25 Hz and slow3.In 0.01–0.25 Hz, the recruitment of VN (T(80) = −2.840, P = 0.006) and AN (T(80) = −2.392, P = 0.019) decreased. In slow3, except for SMN, the recruitment of the other four RSNs are significant. Among them, the recruitment of VN (T(80) = −4.101, P = 0.000) and AN (T(80) = −3.577, P = 0.001) decreased significantly. The recruitment of DMN (T(80) = −2.178, P= 0.032) decreased, while that of LSN (T(80) = 2.456, P = 0.016) increased. The specific brain network information is shown in Table 3.

3.3. Group Comparisons of RSN to RSN Integration at Different Frequencies

After proving the abnormality of SZ at the RSN level, we are still interested in the details between networks. To further investigate the group differences between SZ and NC in five functional networks, we decomposed these values into the RSN × RSN matrix of mean module allegiance. The results showed that the SZ recruitment is lower than that of NC in the full frequency. In slow3, only the recruitment of SZ with LSM is higher than that of NC, and the other modules are lower than that of NC. In addition, we also observed that both SZ and NC prefer to arrange nodes in their functional network rather than through the whole system; that is, they have the characteristics of high recruitment and low integration (Figure 4).
After knowing that the integration of VN is obvious, we still study which pair of networks has a significant difference in integration. In slow3, it is found that the integration between VN and AN (t (80) = 3.429, P = 0.001) (Figure 5A), VN and DMN (t (80) = 3.000, P = 0.004) (Figure 5B) is more obvious in the five brain networks of RSN to RSN after FDR correction. Moreover, the scores of SZ are higher than NC. The specific RSN to RSN integration information is shown in Table 4. (All information values can be viewed in Table 5).

3.4. Group Comparisons of Node Level at Different Frequencies

In addition to the results on the whole brain and RSN levels, we also compared the differences in integration and recruitment at the node level between the NC and SZ. When controlling for multiple comparisons using FDR correction, for integration, in 0.01–0.25 Hz, no significant node differences were observed (Figure 6A); in slow3, they were significantly different in four RSNs modules with a total of sixteen brain regions (Figure 6A). Specifically, the SMN has eight regions. The VN has three regions. The DMN has two regions, and LSN has three regions (Table 6). For recruitment, in 0.01–0.25 Hz, there were found two RSNs modules with a total of six brain regions (Figure 6B). Specifically, the VN has four regions, and LSN has two regions; in slow3, there were thirty-three brain regions in five RSNs modules that were significantly different (Figure 6B). Specifically, the SMN module has one region. The VN has four regions, AN has thirteen regions, DMN has five regions, and LSN has ten regions. The specific node information is shown in Table 7.

3.5. Correlation between Network Measures and SAPS Scores

In order to further understand whether integration and recruitment are related to the severity of SZ, we investigated the association between Scale for Assessment of Positive Symptoms (SAPS) and indicators. In the slow3, based on the Spearman correlations, we found there was a negative correlation between the average recruitment score and the SAPS score in the whole brain level (r = −0.305; p = 0.05) (Figure 7A). In RSN level, there was a negative correlation between the AN and SAPS score (r = −0.305; p = 0.05) (Figure 7B). In node level, there were significant negative correlation between four nodes: inferior frontal gyrus, triangular part (IFGtriang. L) (r = −0.378; p = 0.014), inferior frontal gyrus, triangular part (IFGtriang. R) (r = −0.338; p = 0.029)) and SAPS score (Figure 7C), middle front gyrus (MFG. L) (r = −0.365; p = 0.017), middle frontal gyrus (MFG. R) (r = −0.319; p = 0.04) (Figure 7D).

4. Discussion

To study the frequency-specific dynamic reconfiguration in SZ, we use a sliding time window to construct multilayer brain networks for each band, respectively. Recruitment and integration are introduced to study the differences between NC and SZ. The results showed that compared with the NC, the dynamic reconfiguration of the whole brain in the full frequency of SZ decreased. A more significant decrease is shown in slow3, which means the dynamic reconfiguration of the multilayer brain network of schizophrenia has frequency specificity. These findings may provide a new perspective for explaining the underlying pathological mechanism of SZ.

4.1. Reduced Recruitment in SZ Patients

In terms of whole-brain, the recruitment of SZ decreased significantly in the frequency bands of 0.01–0.25 Hz and slow3. Our results showed that the recruitment ability of SZ is reduced, which affects the cognitive ability of patients. Moreover, this decline can be explained by flexibility. Studies have shown that the dynamic brain network reconfiguration of SZ shows higher flexibility [7,8,9]. In a survey of the genetic risk of dynamic brain network reconfiguration in SZ, this higher flexibility was taken as a manifestation of the lack of organization and stability in the patient’s network module organization [9]. These will lead to instability in SZ.
In addition, at the whole-brain level of the slow3, recruitment is negatively correlated with the SAPS score. This imbalance indicates that the frequency band may regulate the dynamic abnormality of SZ, and the network is disintegrated in these patients, which reduces the alliance preference within the system and the communication between the systems. Besides, SAPS is mainly used to evaluate the positive symptoms of schizophrenia, including hallucinations, delusions, and other symptoms. The higher the SAPA score, the more serious the SZ illness and the worse the SZ cognitive ability.

4.2. Abnormal Brain Networks/Regions of Dynamic Reconfiguration in SZ Patients at slow3

Our results also compared several groups of dynamic network reconfiguration at the RSN level. Firstly, for visual network (VN), we found that the recruitment decreased, and the integration increased at the slow3. The evidence of continuous fusion shows that the collection and transmission of high-frequency visual information may contribute to the bottom-up information integration of the brain [1]. According to previous studies, the reduced activation of the region in the visual magnetic plasma path is associated with deficits in motor processing in SZ [30]. Furthermore, visual stimulation has been proved to increase the correlation between functionally related regions and decrease the correlation between unrelated regions simultaneously [31]. In addition, we found that LING and FFG nodes were remarkably correlated in the slow3, indicating that these nodes were more sensitive to the severity of SZ. Lying between the calcarine sulcus, the LING, alone with the FFG [32] next to it, is directly connected to limbic/paralimbic and subcortical network (LSN) [33,34]. The study of SZ found that the LING is a brain region supporting visual memory [35], and the FFG is related to the identification of facial information or other objects [36,37]. JY Jung and colleagues [34] found abnormal filtering of irrelevant information in visual cortices and altered functional connectivity between the frontoparietal network and visual cortices in SZ. Damage to this area can lead to visual memory dysfunction and visual edge disconnection syndrome. Based on these findings, we believe that the abnormal state of the visual cortex will affect the cognitive function of SZ and lead to a mental defect.
Secondly, for attention network (AN), we find that the recruitment falls off at the slow3, but the integration between attention networks and visual networks is more pronounced. Some studies have shown that the reorganization of brain network modules might contribute to attention processing [38], and visual processing and attention deficit are the main causes of disability in SZ [39]. As a part of the attention network, MFG plays a vital role in the dynamic network changes involved in attention processing. Some studies have shown that changing the critical node MFG of SZ affecting the time-varying brain network will reconfiguration brain network modules [40]. Penghui and colleagues found that the MFG interconnected network in the left hemisphere of SZ plays a leading role in visual and attention network topology [38]. In other words, MFG will promote the network connection between the attention network and other regions to improve the ability of continuous attention. We also found notable differences in recruitment between node MFG and IFG, and there was a negative correlation with SAPA scores. The regions of MFG and IFG are usually related to decision making, action inhibition, and conflict monitoring [41]. In SZ, the IFG plays a crucial role in executive functions such as cognitive inhibition and semantic and linguistic functions [42,43] and is associated with language production and linguistic working memory [44]. We believe that the impairment of brain function in these regions will affect the impairment of the functional cortical network. The resulting brain dysfunction may lead to abnormal clinical and cognitive measures [45].
Thirdly, for the default mode network (DMN), we find that the recruitment decreased at the slow3. In addition, our results also indicate the importance of visual network and default mode network in the relationship between dynamic brain networks. DMN, as a distributed network of brain regions [46], and has a high degree of functional connectivity. In SZ, DMN is usually over-activated and over-connected, leading to cognitive impairment, hallucination, and delusions [47]. Susan and colleagues [47] suggest that the uniqueness of neuropsychiatric disorders may reflect the interaction between DMN and other brain networks. Researchers have also found that the DMN regularly disintegrates into many components in the resting-state, which can act synchronously with the sensorimotor and attention networks [5]. Therefore, we believe that the DMN is easily connected with other brain regions in dynamic reconfiguration.
Finally, for limbic/paralimbic and subcortical network (LSN), we find that the recruitment of node DCG and AMYG in slow3 are significantly correlated. DCG is a significant component of the limbic system and is associated with memory and spatial orientation [48]. This indicates that the decrease in connectivity strength between DCG brain regions in psychiatric patients may affect brain memory. AMYG is nestled in the medial temporal lobe, located anterior to the hippocampus [49]. As a part of the limbic system, AMYG is an all-cortical region involved in emotional processes, learning, and memory. The AMYG plays an essential role in social information processing, particularly in associating emotional salience to sensory stimuli [50] through connections with multiple brain regions such as the prefrontal cortex, motor, and sensory regions [51,52]. We believe that the abnormal association between the limbic system and other systems may lead to the pathophysiology and symptomatology of SZ. Ho and colleagues demonstrated that amygdala orbitofrontal functional connectivity decreased in schizophrenic patients in the study of resting functional connectivity [49].

4.3. Frequency-Specificity of Multilayer Brain Networks in SZ Patients

Our results showed that the dynamic reconfiguration of multilayer brain networks of SZ is different in full frequency. Moreover, the difference in slow3 is more evident than that in full frequency, and some information which is not observable in the full frequency band can even be observed. Therefore, we speculate that the brain network reconstruction of schizophrenia is related to frequency. Previous studies have indicated that higher low-frequency signals may improve the stability of supervisory regulation during the processing of information and high-frequency signals, and frequency band information should be considered in future studies on SZ [1]. Additionally, previous studies on SZ focused on low-frequency (0.01–0.1 Hz) [12,14], which may ignore the information on high-frequency dynamic weighted functional connectivity. Recent neuroimaging studies have also shown that high-frequency (>0.1 Hz) fMRI signals also alter spontaneous neural activity [53].
We suggest that functional integration between brain regions may occur in a specific frequency domain based on these results. Research Indicators (recruitment and integration) are also more observable in specific frequency domains. Suril R. Gohel’s and colleagues’ results show that the functional integration between brain regions measured by BOLD signal correlation occurs over a wider frequency band than was examined in previous studies. This functional integration between brain regions within the same network is specific to one frequency band [54]. This is consistent with our results and may provide new ideas for future SZ research.
By studying the characteristics of each frequency band, we found that the slow3 was more pronounced for SZ. At the same time, Wang and colleagues obtained higher accuracy in slow3 by using the hierarchical sparse learning method to diagnose schizophrenia (SZ) in resting-state functional magnetic resonance imaging (rs-fMRI) [3]. Besides, Zuo and colleagues analyzed the classification weight distribution and concluded that slow3 had the best classification performance while slow2 had the worst [4]. It is considered that the signals in slow3 and slow2 mainly reflect white matter signals and high-frequency physiological noise, while the signs in slow4 and slow5 mainly reflect gray matter signals [55]. Therefore, we believe that slow3 has a specific research value. This proved our conjecture.

5. Conclusions

This study analyzed the dynamic reconfiguration of multilayer brain networks in SZ and normal controls at different frequencies. We found that SZ showed differences in full frequency, which was more significant in slow3. The functional integration between brain regions in SZ was more likely to occur in specific frequency domains. We also found that the recruitment of SZ was significantly abnormal. The decrease in the visual network (VN) and dynamic functional connectivity strength of the attention network (AN) in SZ can decrease the internal links in patients. And this may explain why the RSN level mainly occurred in VN and AN. However, the integration abnormalities in SZ at the RSN level primarily focus on VN. This may be because the collection and transmission of high-frequency visual information contributes to the bottom-up information integration of the brain. In the future, our results may provide potential implications for exploring the neuropathological mechanisms of SZ.

Author Contributions

Conceptualization, Y.Y. and X.C. (Xiaohong Cui); methodology, X.C. (Xiaohong Cui); software, J.X.; validation, Y.Z., X.C. (Xueting Cheng) and T.L.; formal analysis, B.W.; investigation, D.L.; resources, X.C. (Xiaohong Cui) and Y.Z.; data curation, X.C. (Xiaohong Cui) and Y.Z.; writing—original draft preparation, X.C. (Xiaohong Cui) and Y.Z.; writing—review and editing, X.C. and Y.Z.; visualization, X.C. (Xiaohong Cui); supervision, X.C. (Xiaohong Cui); project administration, X.C. (Xiaohong Cui); funding acquisition, J.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by research grants from the National Natural Science Foundation of China (62176177, 61873178, 61906130), Shanxi Province Applied Basic Research Plan (201901D211075, 20210302123112, 20210302123099, 20210302124550, 20210302124166), Scientific research funding project for Returned Overseas students in Shanxi Province (2021-039). Science and Technology Innovation Plan Project of Shanxi Province (2021L045, 2021L046). The sponsors had no role in the design or execution of the study, the collection, management, analysis, and interpretation of the data, or preparation, review, and approval of the manuscript.

Informed Consent Statement

Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

https://openfmri.org/dataset/ds000030/ (accessed on 23 March 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Luo, Y.; He, H.; Duan, M.; Huang, H.; Luo, C.J. Dynamic Functional Connectivity Strength Within Different Frequency-Band in Schizophrenia. Front. Psychiatry 2020, 10, 995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Luo, Y.; He, H.; Duan, M.; Huang, H.; Luo, C.J. Resting-state network connectivity and metastability predict clinical symptoms in schizophrenia. Schizophr. Res. 2018, 201, 208–216. [Google Scholar]
  3. Wang, M.; Hao, X.; Huang, J.; Wang, K.; Shen, L.; Xu, X.; Liu, M.J.N. Hierarchical Structured Sparse Learning for Schizophrenia Identification. Neuroinformatics 2019, 18, 43–57. [Google Scholar] [CrossRef]
  4. Zou, H.; Yang, J.J. Multi-frequency Dynamic Weighted Functional Connectivity Networks for Schizophrenia Diagnosis. Appl. Magn. Reson. 2019, 50, 847–859. [Google Scholar] [CrossRef]
  5. Allen, E.A.; Damaraju, E.; Plis, S.M.; Erhardt, E.B.; Eichele, T.; Calhoun, V.D.J.C.C. Tracking Whole-Brain Connectivity Dynamics in the Resting State. Cereb. Cortex 2014, 24, 663–676. [Google Scholar] [CrossRef]
  6. Farinha, M.; Amado, C.; Morgado, P.; Cabral, J.J.F. Increased excursions to functional networks in schizophrenia in the absence of task. Front. Neurosci. 2022, 16, 821179. [Google Scholar] [CrossRef]
  7. Gifford, G.; Crossley, N.; Kempton, M.J.; Morgan, S.; Mcguire, P.J.N.C. Resting State fMRI Based Multilayer Network Configuration in Patients with Schizophrenia. NeuroImage Clin. 2020, 25, 102169. [Google Scholar] [CrossRef]
  8. Cui, X.; Ding, C.; Wei, J.; Xue, J.; Xiang, J.J. Analysis of Dynamic Network Reconfiguration in Adults with Attention-Deficit/Hyperactivity Disorder Based Multilayer Network. Cereb. Cortex 2021, 31, 4945–4957. [Google Scholar] [CrossRef]
  9. Braun, U.; SchaFer, A.; Bassett, D.S.; Rausch, F.; Schweiger, J.I.; Bilek, E.; Erk, S.; Romanczuk-Seiferth, N.; Grimm, O.; Geiger, L.S. Dynamic brain network reconfiguration as a potential schizophrenia genetic risk mechanism modulated by NMDA receptor function. Proc. Natl. Acad. Sci. USA 2016, 113, 12568–12573. [Google Scholar] [CrossRef] [Green Version]
  10. Xiaosong, H.; Bassett, D.S.; Ganne, C.; Sperling, M.R.; Lauren, K.; Tracy, J.I.N. Disrupted dynamic network reconfiguration of the language system in temporal lobe epilepsy. Brain 2018, 141, 1375–1389. [Google Scholar]
  11. Dong, D.; Duan, M.; Wang, Y.; Zhang, X.; Jia, X.; Li, Y.; Xin, F.; Yao, D.; Luo, C. Reconfiguration of Dynamic Functional Connectivity in Sensory and Perceptual System in Schizophrenia. Brain 2018, 29, 3577–3589. [Google Scholar] [CrossRef] [PubMed]
  12. Xiao, W.; Yan, Z.; Long, Z.; Zheng, J.; Zhang, Y.; Han, S.; Zhao, J.J. Frequency-specific alteration of functional connectivity density in antipsychotic-naive adolescents with early-onset schizophrenia. J. Psychiatr. Res. 2017, 95, 68–75. [Google Scholar]
  13. Han, S.; Zong, X.; Hu, M.; Yu, Y.; Wang, X.; Long, Z.; Chen, H.J. Frequency-selective alteration in the resting-state corticostriatal-thalamo-cortical circuit correlates with symptoms severity in first-episode drug-naive patients with schizophrenia. Schizophr. Res. 2017, 189, 175–180. [Google Scholar] [CrossRef] [PubMed]
  14. Yu, R.; Chien, Y.; Wang, H.S.; Liu, C.; Liu, C.; Hwang, T.; Tseng, W.Y. Frequency-specific alternations in the amplitude of low-frequency fluctuations in schizophrenia. Brain Mapp. 2012, 35, 627–637. [Google Scholar] [CrossRef]
  15. Mattar, M.G.; Cole, M.W.; Thompson-Schill, S.L.; Bassett, D.S. A Functional Cartography of Cognitive Systems. PLoS Comput. Biol. 2015, 11, e1004533. [Google Scholar] [CrossRef] [Green Version]
  16. Cordes, D.; Haughton, V.M.; Arfanakis, K.; Carew, J.D.; Turski, P.A.; Moritz, C.H.; Wendt, G.J. Mapping functionally related regions of brain with functional connectivity in the cerebral cortex in “resting-state” data. Am. J. Neuroradiol. 2001, 21, 1636–1644. [Google Scholar]
  17. Zuo, X.N.; Martino, A.D.; Kelly, C.; Shehzad, Z.E.; Gee, D.G.; Klein, D.F.; Milham, M.P.J.N. The oscillating brain: Complex and reliable. NeuroImage 2010, 49, 1432–1445. [Google Scholar] [CrossRef] [Green Version]
  18. Stam, C.J. Modern network science of neurological disorders. Nat. Rev. Neurosci. 2014, 15, 683. [Google Scholar] [CrossRef]
  19. Tzourio-Mazoyer, N.; Landeau, B.; Papathanassiou, D.; Crivello, F.; Etard, O.; Delcroix, N.; Joliot, M.J.N. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. NeuroImage 2002, 15, 273–289. [Google Scholar] [CrossRef]
  20. Yong, H.; Wang, J.; Liang, W.; Chen, Z.J.; Yan, C.; Yang, H.; Zang, Y.J. Uncovering Intrinsic Modular Organization of Spontaneous Brain Activity in Humans. PLoS ONE 2009, 4, e5226. [Google Scholar]
  21. Leonardi, N.; Ville, D.J.N. On spurious and real fluctuations of dynamic functional connectivity during rest. Neuroimage 2015, 104, 430–436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Kivela, M.; Arenas, A.; Barthelemy, M.; Gl Ee Son, J.P.; Moreno, Y.; Porter, M.A. Multilayer Networks. J. Complex Netw. 2013, 2, 203–271. [Google Scholar] [CrossRef] [Green Version]
  23. Mucha, P.J.; Richardson, T.; Macon, K.; Porter, M.A.; Onnela, J.P.J.S. Community Structure in Time-Dependent, Multiscale, and Multiplex Networks. Science 2010, 328, 876–878. [Google Scholar] [CrossRef]
  24. Jutla, I.S.; Jeub, L.G.; Mucha, P.J. A Generalized Louvain Method for Community Detection Implemented in MATLAB. 2011. Available online: http://netwiki.amath.unc.edu/GenLouvain/GenLouvain (accessed on 23 March 2022).
  25. Bassett, D.S.; Wymbs, N.F.; Rombach, M.P.; Porter, M.A.; Mucha, P.J.; Grafton, S.T. Task-Based Core-Periphery Organization of Human Brain Dynamics. PLoS Comput. Biol. 2013, 9, e1003171. [Google Scholar] [CrossRef] [PubMed]
  26. Hl, A.; Ke, H.; Ypd, E.; Xt, D.; Meng, W.; Bo, M.A.; Ald, G. Dynamic Reconfiguration of Human Brain Networks Across Altered States of Consciousness. Behav. Brain Res. 2021, 419, 113685. [Google Scholar]
  27. Ding, C.; Xiang, J.; Cui, X.; Wang, X.; Li, D.; Cheng, C.; Wang, B.J. Abnormal Dynamic Community Structure of Patients with Attention-Deficit/Hyperactivity Disorder in the Resting State. J. Atten. Disord. 2022, 26, 34–47. [Google Scholar] [CrossRef]
  28. Blondel, V.D.; Guillaume, J.L.; Lambiotte, R.; Lefebvre, E. Fast unfolding of communities in large networks. J. Stat. Mech. Theory Exp. 2008, 2008, P10008. [Google Scholar] [CrossRef] [Green Version]
  29. Chai, L.R.; Mattar, M.G.; Asher, B.I.; Evelina, F.; Bassett, D.S. Functional Network Dynamics of the Language System. Cereb. Cortex 2016, 26, 4148–4159. [Google Scholar] [CrossRef] [Green Version]
  30. Kim, D.; Wylie, G.; Pasternak, R.; Butler, P.D.; Javitt, D.C. Magnocellular contributions to impaired motion processing in schizophrenia. Schizophr. Res. 2006, 82, 1–8. [Google Scholar] [CrossRef] [Green Version]
  31. Nir, Y.; Hasson, U.; Levy, I.; Yeshurun, Y.; Malach, R. Widespread functional connectivity and fMRI fluctuations in human visual cortex in the absence of visual stimulation. NeuroImage 2006, 30, 1313–1324. [Google Scholar] [CrossRef]
  32. Mendoza, J.E.; Foundas, A.L. The Somatosensory Systems. In Clinical Neuroanatomy: A Neurobehavioral Approach; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  33. Conrad, C.D.; Stumpf, W.E. Direct visual input to the limbic system: Crossed retinal projections to the nucleus anterodorsalis thalami in the tree shrew. Exp. Brain Res. 1975, 23, 141–149. [Google Scholar] [CrossRef] [PubMed]
  34. Jung, J.Y.; Kang, J.; Won, E.; Nam, K.; Ham, B.J. Impact of lingual gyrus volume on antidepressant response and neurocognitive functions in Major Depressive Disorder: A voxel-based morphometry study. J. Affect. Disord. 2014, 169, 179–187. [Google Scholar] [CrossRef] [PubMed]
  35. Zhang, L.; Qiao, L.; Chen, Q.; Yang, W.; Xu, M.; Yao, X.; Dong, Y.J. Gray Matter Volume of the Lingual Gyrus Mediates the Relationship between Inhibition Function and Divergent Thinking. Front. Psychol. 2016, 7, 1532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Liu, B.; Feng, Y.; Yang, M.; Chen, J.Y.; Li, J.; Huang, Z.C.; Zhang, L.L. Functional Connectivity in Patients With Sensorineural Hearing Loss Using Resting-State MRI. Am. J. Audiol. 2015, 24, 145. [Google Scholar] [CrossRef] [PubMed]
  37. Takahashi, T.; Zhou, S.Y.; Nakamura, K.; Tanino, R.; Suzuki, M.J. Psychiatry, B. A follow-up MRI study of the fusiform gyrus and middle and inferior temporal gyri in schizophrenia spectrum. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 2011, 35, 1957–1964. [Google Scholar] [CrossRef] [PubMed]
  38. Song, P.; Lin, H.; Liu, C.; Jiang, Y.; Lin, Y.; Xue, Q.; Xu, P.; Wang, Y. Transcranial Magnetic Stimulation to the Middle Frontal Gyrus During Attention Modes Induced Dynamic Module Reconfiguration in Brain Networks. Front. Neuroinform. 2019, 13, 22. [Google Scholar] [CrossRef] [PubMed]
  39. Arkin, S.C.; Ruiz-Betancourt, D.; Jamerson, E.C.; Smith, R.T.; Patel, G.H. Deficits and compensation: Attentional control cortical networks in schizophrenia. NeuroImage Clin. 2020, 27, 102348. [Google Scholar] [CrossRef] [PubMed]
  40. Juha, G.; Rasmus, Z.; Mikko, N.; Antti, P.; Synnve, C.J.C.C. Neural Substrate for Metacognitive Accuracy of Tactile Working Memory. Cereb. Cortex 2017, 11, 5343–5352. [Google Scholar]
  41. Maxim, K.; Alexander, K.; Natalia, M.; Ruslan, M.; Svyatoslav, M.J.F. Deceptive but Not Honest Manipulative Actions Are Associated with Increased Interaction between Middle and Inferior Frontal gyri. Front. Neurosci. 2017, 11, 482. [Google Scholar]
  42. Swick, D.; Ashley, V.; Turken, A.U.J.B.N. Left inferior frontal gyrus is critical for response inhibition. BMC Neurosci. 2008, 9, 102. [Google Scholar] [CrossRef] [Green Version]
  43. Kuperberg, G.R.; Broome, M.; Mcguire, P.K.; David, A.S.; Eddy, M.; Ozawa, F.; Kouwe, A.J.A. Regionally Localized Thinning of the Cerebral Cortex in Schizophrenia. Arch. Gen. Psychiatry 2003, 60, 878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Paulesu, E.; Frith, C.D.; Frackowiak, R.J.N. The neural correlates of the verbal component of working memory. Nature 1993, 362, 342–345. [Google Scholar] [CrossRef] [PubMed]
  45. Kim, S.; Kim, Y.W.; Shim, M.; Jin, M.J.; Lee, S.H.J.F. Altered Cortical Functional Networks in Patients With Schizophrenia and Bipolar Disorder: A Resting-State Electroencephalographic Study. Front. Psychiatry 2020, 11, 661. [Google Scholar] [CrossRef] [PubMed]
  46. Fox, M.D.; Snyder, A.Z.; Vincent, J.L.; Corbetta, M.; Van Essen, D.C.; Raichle, M.E.J.P. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc. Natl. Acad. Sci. USA 2005, 102, 9673–9678. [Google Scholar] [CrossRef] [Green Version]
  47. Whitfield-Gabrieli, S.; Ford, J.M. Default Mode Network Activity and Connectivity in Psychopathology, in Annual Review of Clinical Psychology, Vol 8, S. NolenHoeksema, Editor. Annu. Rev. Clin. Psychol. 2012, 8, 49–76. [Google Scholar] [CrossRef]
  48. Song, K.; Li, J.; Zhu, Y.; Ren, F.; Huang, Z. Altered small-world functional network topology in patients with optic neuritis: A resting-state fMRI study. Dis. Markers 2020, 2021, 9948751. [Google Scholar] [CrossRef]
  49. Ho, N.F.; Chong, P.; Lee, D.R.; Chew, Q.H.; Chen, G.; Sim, K.J.H.R. The Amygdala in Schizophrenia and Bipolar Disorder: A Synthesis of Structural MRI, Diffusion Tensor Imaging, and Resting-State Functional Connectivity Findings. Harv. Rev. Psychiatry 2019, 27, 150–164. [Google Scholar] [CrossRef]
  50. Katharina, S.; Stephan, B.; Tim, V.; Andrea, F.; Roland, W.; Müri, R.; Sebastian, W.J.S.B. Limbic Interference During Social Action Planning in Schizophrenia. Schizophr. Bull. 2018, 44, 359–368. [Google Scholar]
  51. Zhang, M.; Yang, F.; Fan, F.; Wang, Z.; Hong, L.E.J.N.C. Abnormal amygdala subregional-sensorimotor connectivity correlates with positive symptom in schizophrenia. NeuroImage Clin. 2020, 26, 102218. [Google Scholar] [CrossRef]
  52. Berman, R.A.; Gotts, S.J.; McAdams, H.M.; Greenstein, D.; Lalonde, F.; Clasen, L.; Clasen, L.; Watsky, R.E.; Shora, L.; Ordonez, A.E.; et al. Disrupted sensorimotor and social-cognitive networks underlie symptoms in childhood-onset schizophrenia. Brain 2016, 139, 276–291. [Google Scholar] [CrossRef] [Green Version]
  53. Yu, X.; Yuan, B.; Cao, Q.; An, L.; Wang, P.; Vance, A.; Sun, L.J. Frequency-specific abnormalities in regional homogeneity among children with attention deficit hyperactivity disorder: A resting-state f MRI study. Sci. Bull. 2016, 61, 682–692. [Google Scholar] [CrossRef] [Green Version]
  54. Gohel, S.R.; Biswal, B.B. Functional integration between brain regions at rest occurs in multiple-frequency bands. Brain Connect. 2015, 5, 23–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Xiao, H.; Ying-Hang, W.U.; Ping, N.I.; Li-Yuan, F.U.; Hui, L.I.; Chen, Z.Q. Resting-State Functional MRI Study on the Low-Frequency Fluctuation within Different Band of Amplitude of Alzheimer’s Disease. China Med. Devices 2014, 11, 5–10. [Google Scholar]
Figure 1. Schematic overview of analysis strategy. The resting-state fMRI data preprocessing results were first divided into 90 brain regions using the existing standard brain profile in the AAL template. Each brain region represented a node in the network. The brain’s physiological signals are decomposed into five frequency bands. The sliding window technique divides the time series into shorter time intervals. The functional connections within each layer are estimated using Pearson correlation. They connect the same nodes in adjacent periods and build a multilayer network for each participant. The dynamic community structure is detected by maximizing the multilayer modular quality function. We calculated the module allegiance matrix, recruitment, and integration to analyze the difference between NC and SZ.
Figure 1. Schematic overview of analysis strategy. The resting-state fMRI data preprocessing results were first divided into 90 brain regions using the existing standard brain profile in the AAL template. Each brain region represented a node in the network. The brain’s physiological signals are decomposed into five frequency bands. The sliding window technique divides the time series into shorter time intervals. The functional connections within each layer are estimated using Pearson correlation. They connect the same nodes in adjacent periods and build a multilayer network for each participant. The dynamic community structure is detected by maximizing the multilayer modular quality function. We calculated the module allegiance matrix, recruitment, and integration to analyze the difference between NC and SZ.
Brainsci 12 00727 g001
Figure 2. The difference of recruitment and integration from full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz) in the whole brain level between NC and SZ. (A) Integration. (B) Recruitment. Asterisk indicates pairwise group differences; * represents p < 0.05; ** represents p < 0.01.
Figure 2. The difference of recruitment and integration from full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz) in the whole brain level between NC and SZ. (A) Integration. (B) Recruitment. Asterisk indicates pairwise group differences; * represents p < 0.05; ** represents p < 0.01.
Brainsci 12 00727 g002
Figure 3. The difference of integration and recruitment from full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz) in the RSN level between NC and SZ. (A) Integration. (B) Recruitment. Asterisk indicates pairwise group differences; *** represents p < 0.001, ** represents p < 0.01, * represents p < 0.05.
Figure 3. The difference of integration and recruitment from full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz) in the RSN level between NC and SZ. (A) Integration. (B) Recruitment. Asterisk indicates pairwise group differences; *** represents p < 0.001, ** represents p < 0.01, * represents p < 0.05.
Brainsci 12 00727 g003
Figure 4. Matrix of module allegiance in NC and SZ. (A) full frequency (0.01–0.25 Hz). (B) slow3 (0.073–0.198 Hz). P represents the mean value of module allegiance. The main diagonal represents the recruitment coefficient, and the upper/lower triangle represents the integration coefficient.
Figure 4. Matrix of module allegiance in NC and SZ. (A) full frequency (0.01–0.25 Hz). (B) slow3 (0.073–0.198 Hz). P represents the mean value of module allegiance. The main diagonal represents the recruitment coefficient, and the upper/lower triangle represents the integration coefficient.
Brainsci 12 00727 g004
Figure 5. Group differences of integration for each pair of RSN to RSN in slow3 (0.073–0.198 Hz). (A) the integration between Visual network(VN) and Attention network(AN). (B) the integration between Visual network(VN) and Default mode network(DMN). Asterisk indicates pairwise group differences; *** denotes p < 0.001, ** denotes p < 0.01.
Figure 5. Group differences of integration for each pair of RSN to RSN in slow3 (0.073–0.198 Hz). (A) the integration between Visual network(VN) and Attention network(AN). (B) the integration between Visual network(VN) and Default mode network(DMN). Asterisk indicates pairwise group differences; *** denotes p < 0.001, ** denotes p < 0.01.
Brainsci 12 00727 g005
Figure 6. Difference between the NC and SZ in node vulnerability at full frequency (0.01–0.25 Hz) and slow3 (0.073–0.198 Hz). (A) Integration. (B) Recruitment.
Figure 6. Difference between the NC and SZ in node vulnerability at full frequency (0.01–0.25 Hz) and slow3 (0.073–0.198 Hz). (A) Integration. (B) Recruitment.
Brainsci 12 00727 g006
Figure 7. Spearman correlations between dynamic properties and ASRS scores in slow3 (0.073–0.198 Hz). (A) Recruitment (B) Recruitment in the attention network (AN). (C) Recruitment in the inferior frontal gyrus, triangular part (IFGtriang) and (D) middle front gyrus (MFG).
Figure 7. Spearman correlations between dynamic properties and ASRS scores in slow3 (0.073–0.198 Hz). (A) Recruitment (B) Recruitment in the attention network (AN). (C) Recruitment in the inferior frontal gyrus, triangular part (IFGtriang) and (D) middle front gyrus (MFG).
Brainsci 12 00727 g007
Table 1. Demographic and clinical characteristics.
Table 1. Demographic and clinical characteristics.
CharacteristicSZNCStatistical Test
Number of subjects4240--
Age (years)35.19 ± 8.3732.25 ± 8.81P = 0.125
Sex (male/female)30/1225/15P = 0.396
SAPS30.67 ± 22.26----
Note: The values are denoted as mean ± standard deviation.
Table 2. The number of modules in full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz).
Table 2. The number of modules in full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz).
NCSZ
0.01–0.25 HzSlow5Slow4Slow3Slow20.01–0.25 HzSlow5Slow4Slow3Slow2
5.26.467.44.65.874.85.44.6
7.65.476.67.26.666.87.65.4
6.66.26.25.86.25.27.66.26.44.8
68.66.25.45.26.4775.66
77.86.87.66.26.26.266.44.6
6.654.66.85.26666.84.4
77.87.47.86.8665.27.85.8
5.84.86.48.26.85.27.67.454.8
5.87.4675.85.65.45.46.27
7.47.67.46.46.47.276.47.86.8
6.26.877.25.46.64.26.46.24.8
675.46.64.85.26.43.66.85
7.467.28.25.665.85.25.86
7.47.475.87.86.45.26.27.65.6
6.86.28.26.25.87.46.67.29.47
7.267.24.86.46.46.477.25
7.47.26.87.27.685.6794.8
7.88.67.65.476.665.65.47
6.46.686.26.26.46.86.87.65.2
6.26.86.27.85.27.28.66.874.2
65.85.66.46.87.65.27.26.46.4
75.85.27.46.25.65.246.45.4
75.254.85.26.26.27.48.26.2
5.87.47.257.67.86.46.666.2
6.275.27.67.256.65.45.86.6
5.45.85.47.46.26.67.66.884.8
7.87.27.47.44.47.87.26.86.24.6
6.86.86.47.66.865.85.67.46.8
7.66.675.85.46.25.25.666.8
5.66.68.26.64.86.8776.46.6
6.66.47.27.47.45.66.44.69.24.6
5.67.65.46.66.26.2786.44.6
76.85.46.25.277.45.25.45.2
6.27.66.67.26.46.27.26.47.26.8
77.87.47.67.26.27.47.65.45.4
6.46.24.85.256.67.47.27.66
6.69.67.88.467.47.46.86.84.6
7.28.6875.66.86.877.84
6.87.28.26.666.86.67.47.84.6
5.67.865.876.88.26.474.4
\\\\\6.26.258.46.6
\\\\\86.888.25.6
Table 3. Integration and Recruitment of brain networks in full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz).
Table 3. Integration and Recruitment of brain networks in full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz).
Characteristic0.01–0.25 HzSlow5Slow4Slow3Slow2
IntegrationSMNT = −1.277
P = 0.205
T = 0.478
P = 0.634
T = 0.171
P = 0.865
T = −0.246
P = 0.807
T = 0.073
P = 0.942
VNT = 1.012
P = 0.314
T = 0.279
P = 0.781
T = 2.263
P = 0.026 *
T = 0.522
P = 0.603
T = 0.522
P = 0.603
ANT = −0.781
P = 0.437
T = 0.639
P = 0.525
T = 0.378
P = 0.707
T = 0.109
P = 0.914
T = 1.124
P = 0.264
DMNT = −1.010
P = 0.316
T = 0.455
P = 0.651
T = 0.048
P = 0.962
T = −0.564
P = 0.574
T = 1.394
P = 0.167
LSNT = −1.610
P = 0.111
T = −0.041
P = 0.967
T = 1.457
P = 0.149
T = −0.576
P = 0.566
T = 1.031
P = 0.305
RecruitmentSMNT = −0.783
P = 0.436
T = 1.761
P = 0.082
T = −0.677
P = 0.501
T = −0.196
P = 0.845
T = 0.889
P = 0.377
VNT = −2.840
P = 0.006 **
T = −0.701
P = 0.485
T = −4.101
P = 0.000 ***
T = −0.830
P = 0.409
T = −0.272
P = 0.787
ANT = −2.392
P = 0.019 *
T = 0.07
P = 0.945
T = −3.557
P = 0.001 ***
T = −1.266
P = 0.209
T = 0.195
P = 0.846
DMNT = −1.123
P = 0.265
T = 1.352
P = 0.180
T = −2.178
P = 0.032 *
T = 0.390
P = 0.698
T = −0.803
P = 0.424
LSNT = −0.713
P = 0.478
T = 1.489
P = 0.140
T = 2.456
P = 0.016 *
T = −0.431
P = 0.668
T = 0.163
P = 0.871
Asterisk indicates pairwise group differences; *** represents p < 0.001, ** represents p < 0.01, * represents p < 0.05.
Table 4. RSN-to-RSN integration scores with significant differences in slow3 (0.073–0.198 Hz).
Table 4. RSN-to-RSN integration scores with significant differences in slow3 (0.073–0.198 Hz).
RSN1RSN2NC (SD)SZ (SD)P (FDR)
Visual network (VN)Attention network (AN)0.194 (0.055)0.235 (0.055)0.000
Visual network (VN)Default mode network (DMN)0.212 (0.049)0.240 (0.036)0.001
Table 5. RSN-to-RSN integration scores in full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz).
Table 5. RSN-to-RSN integration scores in full frequency (0.01–0.25 Hz), slow5 (0.01–0.027 Hz), slow4 (0.027–0.073 Hz), slow3 (0.073–0.198 Hz), and slow2 (0.198–0.25 Hz).
FrequencyRSN1RSN2NC (SD)SZ (SD)P (FDR)
0.01–0.25 Hzvisual network (VN)somatosenery/motor and auditory network (SMN)0.258 (0.115)0.267 (0.132)0.731
visual network (VN)attention network (AN)0.168 (0.082)0.183 (0.062)0.367
visual network (VN)default mode network (DMN)0.201 (0.071)0.213 (0.055)0.774
visual network (VN)limbic/paralimbic and subcortical network (LSN)0.165 (0.069)0.180 (0.056)0.312
Slow5visual network (VN)somatosenery/motor and auditory network (SMN)0.348 (0.110)0.344 (0.126)0.858
visual network (VN)attention network (AN)0.292 (0.088)0.295 (0.069)0.839
visual network (VN)default mode network (DMN)0.321 (0.078)0.329 (0.068)0.639
visual network (VN)limbic/paralimbic and subcortical network (LSN)0.294 (0.067)0.307 (0.089)0.469
Slow4visual network (VN)somatosenery/motor and auditory network (SMN)0.311 (0.119)0.299 (0.108)0.613
visual network (VN)attention network (AN)0.210 (0.075)0.215 (0.071)0.706
visual network (VN)default mode network (DMN)0.254 (0.069)0.260 (0.069)0.666
visual network (VN)limbic/paralimbic and subcortical network (LSN)0.232 (0068)0.236 (0.071)0.774
Slow3visual network (VN)somatosenery/motor and auditory network (SMN)0.271 (0.074)0.271 (0.071)0.552
visual network (VN)attention network (AN)0.194 (0.055)0.235 (0.055)0.000
visual network (VN)default mode network (DMN)0.212 (0.049)0.240 (0.036)0.001
visual network (VN)limbic/paralimbic and subcortical network (LSN)0.220 (0.052)0.215 (0.043)0.797
Slow2visual network (VN)somatosenery/motor and auditory network (SMN)0.313 (0.055)0.311 (0.061)0.856
visual network (VN)attention network (AN)0.294 (0.035)0.297 (0.054)0.759
visual network (VN)default mode network (DMN)0.284 (0.041)0.290 (0.044)0.458
visual network (VN)limbic/paralimbic and subcortical network (LSN)0.269 (0.042)0.280 (0.051)0.271
Table 6. Brain map of regions with significant differences on integration in slow3 (0.073–0.198 Hz).
Table 6. Brain map of regions with significant differences on integration in slow3 (0.073–0.198 Hz).
FrequencyNetworkNameAbbROINC (SD)SZ (SD)P
(FDR)
Slow3SMNRolandic_OperROL.L170.289 (0.053)0.245 (0.053)0.001
ROL.R180.279 (0.041)0.243 (0.053)0.011
InsulaINS.L290.300 (0.055)0.258 (0.062)0.015
INS.R300.300 (0.055)0.260 (0.060)0.015
HeschlHES.L790.283 (0.049)0.244 (0.052)0.011
HES.R800.280 (0.044)0.244 (0.048)0.011
Temporal_SupSTG.L810.286 (0.051)0.253 (0.055)0.036
STG.R820.284 (0.049)0.243 (0.047)0.000
VNLingualLING.L470.212 (0.053)0.248 (0.050)0.015
LING.R480.208 (0.052)0.239 (0.054)0.045
Fusiform_RFFG.R560.241 (0.058)0.288 (0.047)0.000
DMNCingulum_AntACG.L310.297 (0.052)0.260 (0.052)0.015
ACG.R320.297 (0.054)0.263 (0.052)0.032
LSNAmygdala_RAMYG.R420.272 (0.056)0.233 (0.046)0.011
Temporal_Pole_MidTPOmid.L870.280 (0.052)0.248 (0.046)0.028
TPOmid.R880.294 (0.058)0.247 (0.041)0.000
Table 7. Brain map of regions with significant differences on recruitment in full frequency (0.01–0.25 Hz) and slow3 (0.073–0.198 Hz).
Table 7. Brain map of regions with significant differences on recruitment in full frequency (0.01–0.25 Hz) and slow3 (0.073–0.198 Hz).
FrequencyNetworkNameAbbROINC (SD)SZ (SD)P
(FDR)
0.01–0.25 HzVNLingual_LLING.L470.743 (0.143)0.634 (0.169)0.037
Occipital_SupSOG.L490.769 (0.114)0.676 (0.145)0.037
SOG.R500.769 (0.118)0.661 (0.158)0.037
Fusiform_RFFG.R560.534 (0.247)0.371 (0.219)0.037
LSNCingulum_MidDCG.L330.287 (0.131)0.200 (0.097)0.037
DCG.R340.281 (0.122)0.199 (0.106)0.037
Slow3SMNInsula_LINS.L290.456 (0.116)0.395 (0.110)0.048
VNLingual_LLING.L470.610 (0.129)0.477 (0.143)0.001
Occipital_Inf_LIOG.L530.494 (0.157)0.405 (1.173)0.048
FusiformFFG.L550.312 (0.166)0.232 (0.112)0.040
FFG.L560.349 (0.176)0.233 (0.115)0.004
ANFrontal_MidMFG.L70.421 (0.960)0.355 (0.093)0.010
MFG.R80.429 (0.088)0.357 (0.075)0.002
Frontal_Inf_OperIFGoperc.L110.384 (0.112)0.328 (0.093)0.048
IFGoperc.R120.375 (0.107)0.317 (0.086)0.031
Frontal_Inf_TriIFGtriang.L130.434 (0.091)0.357 (0.098)0.003
IFGtriang.R140.423 (0.095)0.345 (0.086)0.002
Frontal_Inf_OrbORBinf.L150.394 (0.085)0.346 (0.088)0.048
ORBinf.R160.399 (0.101)0.334 (0.082)0.011
Parietal_InfIPL.L610.365 (0.101)0.297 (0.068)0.004
IPL.R620.369 (0.106)0.317 (0.067)0.034
AngularANG.L650.365 (0.101)0.313 (0.074)0.034
ANG.R660.372 (0.099)0.316 (0.082)0.023
Temporal_Inf_LITG.L890.309 (0.093)0.260 (0.089)0.048
DMNFrontal_Sup_MedialSFGmed.L230.442 (0.096)0.385 (0.078)0.017
SFGmed.R240.446 (0.091)0.387 (0.081)0.011
Cingulum_Post_LPCG.L350.340 (0.111)0.287 (0.081)0.048
Temporal_MidMTG.L850.341 (0.108)0.268 (0.076)0.004
MTG.R860.324 (0.112)0.258 (0.072)0.010
LSNCingulum_MidDCG.L330.270 (0.107)0.184 (0.074)0.001
DCG.R340.267 (0.101)0.189 (0.073)0.002
HippocampusHIP.L370.396 (0.075)0.479 (0.099)0.001
HIP.R380.389 (0.080)0.483 (0.105)0.001
ParaHippocampalPHG.L390.395 (0.071)0.468 (0.103)0.003
PHG.R400.389 (0.074)0.470 (0.106)0.002
AmygdalaAMYG.L410.392 (0.082)0.466 (0.108)0.005
AMYG.R420.386 (0.089)0.469 (0.108)0.003
Temporal_Pole_MidTPOmid.L870.324 (0.085)0.423 (0.131)0.002
TPOmid.R880.336 (0.099)0.409 (0.123)0.017
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, Y.; Zhang, Y.; Xiang, J.; Wang, B.; Li, D.; Cheng, X.; Liu, T.; Cui, X. Frequency-Specific Analysis of the Dynamic Reconfiguration of the Brain in Patients with Schizophrenia. Brain Sci. 2022, 12, 727. https://doi.org/10.3390/brainsci12060727

AMA Style

Yang Y, Zhang Y, Xiang J, Wang B, Li D, Cheng X, Liu T, Cui X. Frequency-Specific Analysis of the Dynamic Reconfiguration of the Brain in Patients with Schizophrenia. Brain Sciences. 2022; 12(6):727. https://doi.org/10.3390/brainsci12060727

Chicago/Turabian Style

Yang, Yanli, Yang Zhang, Jie Xiang, Bin Wang, Dandan Li, Xueting Cheng, Tao Liu, and Xiaohong Cui. 2022. "Frequency-Specific Analysis of the Dynamic Reconfiguration of the Brain in Patients with Schizophrenia" Brain Sciences 12, no. 6: 727. https://doi.org/10.3390/brainsci12060727

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop