Abstract
Introduction:
Convergent evidence has demonstrated a shared rich-club reorganization across multiple major psychiatric conditions. However, previous studies assessing altered functional couplings between rich-club regions have typically focused on the mean time series from entire functional magnetic resonance imaging (fMRI) scanning session, neglecting their time-varying properties.
Methods:
In this study, we aim to explore the common and/or unique alterations in the temporal variability of rich-club organization among schizophrenia (SZ), bipolar disorder (BD), and attention deficit/hyperactivity disorder (ADHD). We employed a temporal rich-club (TRC) approach to quantitatively assess the propensity of well-connected nodes to form simultaneous and stable structures in a temporal network derived from resting-state fMRI data of 156 patients with major psychiatric disorders (SZ/BD/ADHD = 71/45/40) and 172 healthy controls. We executed the TRC workflow at both whole-brain and subnetwork scales across varying network sparsity, sliding window strategies, lengths and steps of sliding windows, and durations of TRC coefficients.
Results:
The SZ and BD groups displayed significantly decreased TRC coefficients compared to corresponding HC groups at the whole-brain scale and in most subnetworks. In contrast, the ADHD group exhibited reduced TRC coefficients in longer durations, as opposed to shorter durations, which markedly differs from the SZ and BD groups. These findings reveal both transdiagnostic and illness-specific patterns in temporal variability of rich-club organization across SZ, BD, and ADHD.
Discussion:
TRC may serve as an effective metric for detecting brain network disruptions in particular states, offering novel insights and potential biomarkers into the neurobiological basis underpinning the behavioral and cognitive deficits observed in these disorders.
1. Introduction
Mental disorders, characterized by behavioral or mental patterns causing significant distress or impairment of personal functioning, affect 1 in 8 people and 1 in 4 family worldwide (1, 2). Due to overlapping features, accurate diagnoses of mental disorders, such as schizophrenia (SZ) and bipolar disorder (BD), present challenges (3, 4). Similarly, the differential diagnosis of attention deficit hyperactivity disorder (ADHD) and BD remains difficult due to shared symptoms and high comorbidity rates (5, 6). Understanding the homogeneity and heterogeneity of various mental disorders could elucidate their neurobiological underpinnings and inform the development of targeted diagnostic and treatment strategies.
Recent advances in functional magnetic resonance imaging (fMRI) have highlighted the importance of functional connectivity (FC) in the pathophysiology of mental disorders. Abnormalities in FC have been observed in patients with mental disorders at regional, subnetwork, and whole-brain scales (7–9). Employing graph theory-based approaches, researchers have identified disturbances in brain network topology that serve as valuable classification features for distinguishing patients with mental disorders from healthy controls (10, 11). Rich-club organization, a core feature of the brain networks, has garnered increasing attention in network neuroscience research (12, 13) and has provided novel insights into mental disorders (14–16). This organization refers to a set of brain hub regions with disproportionately high number of edges and interconnections, facilitating efficient communication between brain regions. Altered rich-club organization may indicate disrupted brain function and information transmission in specific disease states (17). Moreover, rich-club organization can reflect the trade-off between costs and benefits in human brain function, indicating the efficiency of brain operation (13). However, the topological properties of rich-club organization in mental disorders remain underexplored, particularly in terms of temporality and simultaneity.
Many previous studies have constructed static functional brain networks, aggregating FC data over entire scanning sessions and failing to capture the temporal dynamics of rich-club organization (15, 16). Static networks may inaccurately represent rich-club organization, as they could include edges active at unrelated times rather than simultaneously (18, 19). Consequently, contemporary research has shifted focus to dynamic FC, capturing time-varying patterns in brain connectivity (20–22). Despite the identification of time-varying features associated with mental disorders, the temporality, simultaneity, and time span of rich-club organization interactions remain largely unexplored.
This study aims to investigate the transdiagnostic and/or illness-specific disruptions of time-varying rich-club organization properties across SZ, BD, and ADHD. We introduce a novel temporal rich-club (TRC) analysis workflow to assess the tendency of well-connected nodes to form simultaneous and stable structures in temporal networks. Our analysis encompasses whole-brain and subnetwork scales and considers various network sparsity levels, sliding window strategies, sliding window lengths, sliding window steps, and duration of TRC coefficients. We hypothesize that: (1) the novel TRC workflow will reveal altered rich-club coefficients in major psychiatric conditions, providing a more nuanced and dynamic perspective across the entire time series and (2) different major psychiatric conditions may be distinguishable based on global or local attributes of TRC coefficients.
2. Materials and methods
2.1. Participants
All participants included in this study were obtained from two open-access datasets. Specifically, the patients with SZ and corresponding heathy controls (HCs) were collected from the Center for Biomedical Research Excellence (COBRE) database.1 The COBRE database contains raw anatomical and functional MR data from 74 patients with SZ and 74 HCs (ages ranging from 18 to 65 in each group). In total, 71 patients with SZ and 74 HCs from the dataset were used for our subsequent experiments, as the class labels of the other three SZ participants were not provided. Similar inclusion criteria can be found in some previous studies (9, 23, 24). The BD, ADHD patients and the corresponding HCs were collected from the Consortium for Neuropsychiatric Phenomics (CNP) database2 (25). After quality control, age- and sex- matching, 45 BD patients, 40 ADHD patients and 98 HCs were included in subsequent analyses. Detailed demographic information can be found in Tables 1, 2.
Table 1
Detailed demographic information for COBRE database.
Two-sample t-test.
Chi-square two-tailed test.
Table 2
Detailed demographic information for CNP database.
Two-sample t-test.
Chi-square two-tailed test.
The p-values-1 are calculated based on ADHD group vs. HC group, while p-values-2 are calculated based on BD group vs. HC group.
2.2. Imaging data preprocessing
All fMRI images were preprocessed using Data Processing & Analysis for Brain Imaging (DPABI) (26). The preprocessing procedure included the removal of the first 10 volumes of functional runs due to fMRI signal instability. Slice time correction, head-motion correction, and co-registration of T1-weighted MRI images and fMRI images were performed for the remaining volumes. Subsequently, all fMRI images were normalized to Montreal Neurological Institute (MNI) space and re-sampled to mm voxels. Smooth ( mm FWHM) and band-pass filter (0.01–0.1 Hz) were applied to the images transformed to the MNI space.
To construct FC matrices, all brain images were parcellated into 160 regions by registering images to Dosenbach’s 160 atlas after data preprocessing (27). For ease of comparison, a canonical division of brain regions into cognitive systems was used (28), dividing the 160 brain regions of each participant into six functional subnetworks: visual network (VSN), sensory-motor network (SMN), dorsal attention network (DAN), ventral attention network (VAN), frontoparietal network (FPN), default mode network (DMN), and subcortical network (SBN). According to the division scheme in Yan et al. (28), there were 14, 36, 28, 32, 21, and 22 functional brain regions for DAN, DMN, FPN, SMN, VAN, and VSN, respectively. We mainly analyzed and discussed the TRC coefficients in DAN, DMN, FPN, SMN, VAN, and VSN in later sections, while the remaining seven brain regions in SBN are not included.
2.3. Construction of FC matrix
This study constructed two types of FC matrices: whole-brain-scale FC matrix and subnetwork-scale FC matrix. For all matrices, the level of FC between each pair of nodes (brain regions) was computed as the Pearson correlation coefficient (PCC) between their averaged regional time series. This classical method can be found in many recent studies (7, 16, 24). The sizes of whole-brain-scale FC matrix and subnetwork-scale FC matrix are represented as and , respectively, where denotes the total number of brain regions of the whole brain, and denotes the total number of brain regions in a specific subnetwork. By setting an absolute or proportional threshold (9, 10), if the connection strength (e.g., based on the correlation values of bold signals) between two brain regions is higher than the threshold (e.g., for PCC), then it is considered that there is an edge connection between the two brain regions (each edge corresponds to a degree).
2.4. Definition of temporal rich-club effect
Before introducing the temporal rich-club effect, let us review the classical rich-club effect. The rich-club effect is defined as the density of edges in the subset of the nodes with degree larger than k, that is, , where is the number of edges among . A larger indicates that nodes have a disproportionately high number of edges and many edges between each other, termed the “rich-club” effect (16, 19, 29).
Many previous brain network studies of fMRI analyze the rich-club effects in an average manner (i.e., using average values across the whole time series or a duration of the time series by dividing sliding windows). However, static networks are often aggregated representations of the resulting temporal networks (30). Thus, the rich-club structures found in such static networks could be formed by edges that were active at unrelated times (18). Since the human brain is a dynamic system, analyzing the dynamic and temporal properties of rich-club organizations may better reflect the complex neuroscience mechanisms (e.g., the hub nodes in the rich-club organizations may change continuously) during a duration of time or the whole time series.
The temporal rich-club (TRC) phenomenon has been observed recently by Pedreschi et al. (19) using different levels of network datasets such as the US air transportation network dataset. Given a temporal network, the TRC quantifies whether nodes that interact with increasing numbers of other nodes tend to interact with each other simultaneously and in a stable way (19). To consider temporality, Pedreschi et al. (19) propose to define at each time the Δ −cohesion . It represents the number of ties (between the nodes of ) that remain stable over , normalized by its maximal possible value . Note that is the instantaneous density between the nodes of . The in this paper refers to duration equals to one fMRI timepoint. Therefore, the TRC coefficient can be easily defined as the maximal density of temporal edges observed in a stable way for a duration Δ among nodes of aggregated degree larger than :where the M(k, Δ) quantifies whether the static rich-club patterns correspond to a structure that actually existed at some instant and is formed by links that appeared at unrelated times and not in a simultaneous way. A M(k, Δ) increasing with k denotes the TRC effects in a simultaneous way for a duration of at least Δ. Given the consideration of complex brain dynamic patterns, performing the TRC analyses at the regional scale may shed light on the understanding of functional temporality of major psychiatric groups.
2.5. Calculation of temporal rich-club coefficients
In this study, we performed both whole-brain-scale TRC analyses and subnetwork-scale TRC analyses. The whole-brain-scale TRC analyses of fMRIs were calculated based on the correlation matrix constructed by PCC (i.e., in this study). To investigate the dynamic TRC relationships, we divided the BOLD signals of the whole time series into different windows (i.e., a duration of bold signals). In later sections, we set k as half of the number of brain regions in the whole brain or subnetworks for whole-brain-scale and subnetwork-scale TRC analyses, respectively. As there is no universally appropriate window size and step of window shifting for fMRI data, we extensively explored adopting both the non-overlapping sliding window and overlapped sliding window methods for the TRC analyses with different window sizes and sliding steps. For non-overlapping sliding window methods, we evaluated the following parameter settings: a. window size = 5 TR, sliding step = 5 TR, sparsity threshold ; b. window size = 5 TR, sliding step = 5 TR, sparsity threshold ; c. window size = 5 TR, sliding step = 5 TR, sparsity threshold; d. window size = 5 TR, sliding step = 5 TR, sparsity threshold ; e. window size = 10 TR, sliding step = 10 TR, sparsity threshold ; f. window size = 10 TR, sliding step = 10 TR, sparsity threshold ; g. window size = 10 TR, sliding step = 10 TR, sparsity threshold ; h. window size = 10 TR, sliding step = 10 TR, sparsity threshold .
For overlapped sliding windows, we evaluated the following parameter settings (all adopting sparsity threshold ): a. window size = 10 TR, sliding step = 2 TR; b. window size = 10 TR, sliding step = 4 TR; c. window size = 20 TR, sliding step = 2 TR; d. window size = 20 TR, sliding step = 4 TR; e. window size = 40 TR, sliding step = 2 TR; f. window size = 40 TR, sliding step = 4 TR; g. window size = 80 TR, sliding step = 2 TR; h. window size = 80 TR, sliding step = 4 TR. In our initial experiments, we have already validated that the impact of window size, sliding step and threshold on TRC (also see Figures 1, 2) are significantly smaller relative to the Delta (Δ) itself. Therefore, to detect the subtle heterogeneity of these major psychiatric groups in stably maintaining the connections between the nodes (brain regions) of , we also performed whole-brain-scale TRC analyses with the altered duration Delta (Δ) from Delta = 1 to Delta = 8, based on the non-overlapping sliding window method with window size = 10 TR, sliding step = 10 TR, and sparsity threshold . A larger duration Δ denotes the connections are maintained for more fMRI timepoints.
Figure 1
Figure 2
Since the results in non-overlapping sliding window strategies and overlapped sliding window strategies in whole-brain-scale TRC analyses display strong patterns in the consistency of results (see Figures 1, 2), the subnetwork-scale TRC analyses were directly performed on non-overlapping sliding window strategies (window size = 10 TR, sliding step = 10 TR) adopting the sparsity threshold from to . Our purpose here is mainly to avoid excessive display of redundant results. The basic flow of this study is depicted in Figure 3. Please note that in scenarios involving multiple comparisons, this study uniformly adopts as the multiple comparison correction threshold.
Figure 3
3. Results
3.1. Whole-brain TRC analysis of major psychiatric groups vs. HC group using non-overlapping sliding window methods
Utilizing a window size of 10 TR, a sliding step of 10 TR, and sparsity threshold of , we constructed group-scale average FC matrices for non-overlapping sliding windows for the COBRE (i.e., SZ and HC groups) and CNP (ADHD, BD and HC groups) databases (see in the Supplementary Figure S1). Based on these sliding windows, we computed whole-brain TRC coefficients for major psychiatric groups relative to their corresponding HC groups (Figure 1). Notably, differences in whole-brain TRC coefficients between major psychiatric groups and HC groups remain consistent for R values ranging from 0.2 to 0.8. Specifically, SZ and BD groups show significantly reduced average whole-brain TRC coefficients compared to their respective HC groups in most windows, indicating impaired TRC organizations. In contrast, TRC relationships between ADHD and HC groups are more complex, with whole-brain TRC coefficients significantly differing only in specific sliding windows.
To verify whether the TRC would be a powerful biomarker to identify psychiatric disorders, we performed subject-scale classification experiments for each major psychiatric group with their corresponding HC group. Table 3 displays all of the quantitative classification metrics based on an AdaBoost ensemble classifier (31) with 3-fold cross validation protocol (2 folds for training and 1 fold for test in each experiment). This protocol has been repeatedly performed for 10 times to obtain the mean results. In order to ensure the adequacy of the experiments, we conducted comparison experiments for the sparsity threshold for R ranging from 0.2 to 0.8. According to Table 3, we can easily note that the ability to recognize SZ (best ACC = 0.868 ± 0.07) using TRC is significantly stronger compared to ADHD (best ACC = 0.709 ± 0.10) or BD (best ). Additionally, the classification performance is better when R is approaching 0.2 rather than 0.8, which indicates that excessive network sparsity should be avoided when conducting TRC-based classification. These findings are consistent with the between-group difference results presented in Figure 1. However, we should also note that compared to the COBRE database, the CNP database has a more biased sample distribution, which may be one of the reasons for the worse classification results for ADHD and BD.
Table 3
| Groups | ACC | SEN | SPE | PPV | NPV |
|---|---|---|---|---|---|
| SZ vs. HC (R = 0.2) | |||||
| SZ vs. HC (R = 0.4) | 0.07 | ||||
| SZ vs. HC (R = 0.6) | |||||
| SZ vs. HC (R = 0.8) | |||||
| ADHD vs. HC (R = 0.2) | |||||
| ADHD vs. HC (R = 0.4) | 0.22 | ||||
| ADHD vs. HC (R = 0.6) | |||||
| ADHD vs. HC (R = 0.8) | |||||
| BD vs. HC (R = 0.2) | |||||
| BD vs. HC (R = 0.4) | |||||
| BD vs. HC (R = 0.6) | |||||
| BD vs. HC (R = 0.8) |
Subject-scale classification performance of the three major psychiatric groups vs. their corresponding HC groups using TRC with non-overlapping sliding window methods.
ACC, accuracy; SEN, sensitivity; SPE, specificity; PPV, positive predictive value; NPV, negative predictive value. All of the values are denoted by Mean with SD.
3.2. Whole-brain TRC analysis of major psychiatric groups vs. HC group using overlapped sliding window methods
We further assessed TRC coefficients for major psychiatric groups relative to HC groups using overlapping sliding window methods, exploring various window sizes (e.g., W = 10, 20, 40, and 80 s) and sliding steps (e.g., or ). Correlation matrices were constructed based on a sparsity threshold of , yielding comparable results with and . As shown in Figure 2, average TRC coefficients for SZ and BD groups were consistently lower than those for corresponding HC groups across all sliding windows (p < 0.05 in all paired t-tests). This pattern persisted across different window sizes (e.g., from to ) and sliding steps (e.g., 2 s or 4 s). Similar to the non-overlapping sliding window methods, TRC coefficients in the between-group analysis of ADHD group vs. HC groups were complex and varied across different sliding windows.
3.3. Whole-brain TRC analysis of major psychiatric groups vs. HC group using different duration
To examine the homogeneity and heterogeneity of the three psychiatric disorders, we calculated TRC coefficients based on different duration (we denote it as Δ in figures). Figure 4 displays TRC coefficients for participants with psychiatric disorders relative to their corresponding HCs. Larger Δ values represent greater stability in the density of connected edges over time. Both SZ and BD groups consistently demonstrated smaller whole-brain average TRC coefficients than their corresponding HC groups (value of p < 0.05), with these patterns persisting from Δ = 1 to Δ = 8. In contrast, significantly smaller whole-brain TRC coefficients for the ADHD group compared to HC group were only observed for Δ = 8 and not for other smaller Δ values. These results suggest the underlying mechanisms of ADHD differ from those of SZ and BD, with reduced TRC coefficients only evident in longer durations.
Figure 4
3.4. Subnetwork-scale TRC analysis of major psychiatric groups vs. HC groups using non-overlapping sliding window methods
In addition to whole-brain scale TRC analyses, we conducted subnetwork-scale TRC analyses to identify unique TRC alteration patterns in specific major psychiatric groups at a finer scale. For clarity, we illustrated the subnetwork-scale correlation matrices for the six subnetworks at non-overlapping window 1 and 14 (based on a window size = 10 TR and a sliding step of 10 TR) for the COBRE database in Figure 5. The sizes of FC matrices for DAN, DMN, FPN, SMN, VAN, and VSN were , , , , , and , respectively.
Figure 5
The dynamic change patterns of all subnetwork-scale TRC coefficients for COBRE database and CNP database were also displayed in Figure 6. To ensure the stability of the subnetwork-scale TRC results, we conducted experiments using subnetwork-scale correlation matrices with three sparsity thresholds ranging from to . Specifically, the SZ group exhibited significantly decreased TRC coefficients compared to the HC group for most sliding windows in DAN, DMN, SMN, VAN, and VSN, except for FPN; the ADHD group displayed fluctuating in TRC coefficients for single sliding windows compared to HC group in all subnetworks. However, the BD group showed significantly decreased TRC coefficients compared to the HC group for most sliding windows in DAN, DMN, FPN, VAN, and VSN, with a few exceptions in the SMN. These results indicate that the pathological patterns of the SZ group and the BD group are more similar to each other than the HC group in terms of subnetwork-scale TRC differences.
Figure 6
We further illustrated the subnetwork-scale mean TRC alterations for SZ, ADHD and BD groups in Figure 5, which help to demonstrate the overall trends of subnetwork-scale TRC coefficients for major psychiatric groups relative to corresponding HC groups. The top three subnetworks with the highest TRC coefficients were VSN (), SMN (), and DAN () for the SZ group; VSN (), SMN (), and VAN () for the ADHD group; and VSN (), SMN (), and VAN () for the BD group. Results indicate that all major psychiatric groups had significantly higher TRC coefficients in VSN and SMN than their corresponding HC groups ( for t-tests). Quantitative comparisons can be found in Figures 5D–F, which reveal that both the SZ and BD groups exhibited significantly decreased average TRC coefficients in all six subnetworks compared to their corresponding HC groups (value of p < 0.01). However, for the ADHD group, only the TRC coefficients of VSN were significantly higher relative to HC group (value of p < 0.05).
4. Discussion
In this study, we present a novel TRC analysis approach to uncover the altered brain network patterns in ADHD, BD and SZ at both whole-brain and subnetwork scales. We considered various parameters such as network sparsity, sliding window strategies, window length, window step, and TRC coefficient duration. Our findings suggest that TRC coefficients may serve as effective biomarker to differentiate major psychiatric groups with similar phenotypes, and can complement the dynamic properties absent in existing graph theory metrics, particularly in the context of mental disorders.
At the whole-brain scale, we observed that SZ and BD groups displayed significantly decreased average whole-brain TRC coefficients in comparison to their corresponding HC groups across most windows, reflecting impaired TRC organizations. These findings were consistent across different network sparsity levels, sliding window strategies, window lengths, window steps, and TRC coefficient durations. In contrast, the ADHD group exhibited reduced TRC coefficients primarily in longer durations, indicating a fundamentally different TRC mechanism compared to SZ and BD groups (32). The results of subject-scale classification experiments indicate that TRC can serve as an important biomarker to distinguish psychiatric patients from the corresponding HCs. However, when the sample distribution is biased (e.g., ADHD and BD subjects in CNP database), the classification performance may significantly degrade. Collectively, the above results suggest that the three major psychiatric disorders are characterized by a diminished capacity for dynamic information processing, reflecting distinct cognitive processes.
At the subnetwork scale, we found that the SZ group exhibited significantly decreased TRC coefficients compared to the HC group across most sliding windows in DAN, DMN, SMN, VAN, and VSN, with the exception of FPN. This suggests that the SZ group’s brain network is less interconnected in a simultaneous and stable manner, indicating instability and fragility. The BD group, on the other hand, displayed significantly decreased TRC coefficients compared to the HC group across most sliding windows in DAN, DMN, FPN, VAN, and VSN, with a few exceptions in SMN. The presence of highly connected TRCs implies efficient information integration across various brain regions; therefore, significantly decreased TRC coefficients in BD individuals may critically impair the exchange of information between subnetworks and the coordination of cognitive processes (33, 34). Our results highlight the potential of subnetwork-scale TRC coefficients in FPN and SMN as useful biomarkers to differentiate SZ from the BD groups. Moreover, subnetwork-scale TRC coefficients for the ADHD group were also distinguishable, displaying significantly increased TRC coefficients in VSN compared to the HC group, a marked contrast to the decreased subnetwork-scale TRC coefficients observed in most subnetworks for SZ and BD groups.
Our findings align well with previous studies investigating FC in these disorders. For instance, prior SZ studies reported decreased FCs in brain regions primarily located in the DAN (35), FPN (36, 37), and SMN (38); BD studies found aberrant FCs in the DMN (39), FPN (40), and VAN (41); and ADHD studies identified hyperactivity in the VSN (32, 42). If the previous research results are compared to “pictures,” the results of this study are more like the “videos” corresponding to these pictures, which may own more information. For example, in the “picture” direction, both the SZ and BD groups displayed significantly decreased average TRC coefficients in all of the six subnetworks than their corresponding HC groups (see Figures 5D–F), but the dynamic change patterns (in the “video” direction) of TRC coefficients are different for SZ and BD groups: SZ group showed persistently decreased TRC coefficients in SMN and VSN, while the BD group showed persistently decreased TRC coefficients only in VSN (see Figures 5A,B). Also, in the “picture” direction, ADHD group showed significantly increased average TRC coefficients in VSN (see Figure 5E), but as for the dynamic change patterns (in the “video” direction), the ADHD group was usually indistinguishable from the HC group (see Figure 5B). All these evidences have shown the complementarity of TRC results to previous studies. Compared to these previous results, our study offers a more dynamic and informative perspective on brain network alterations in mental disorders. The complementarity of TRC results to previous research is evident in both static and dynamic patterns.
Altered subnetwork-scale TRC coefficients may accurately reflect disorder-specific impairments in neural circuit mechanisms. For example, previous studies has shown that decreased connections in VSN are common features in both SZ and BD groups, and that individuals with BD may compensate for disrupted VSN connectivity while those with SZ do not (43). Our results show that mean TRC gaps between BD and HC groups are , which are significantly smaller than mean TRC gaps between SZ and HC groups (). Obviously, smaller TRC gaps are easier to be compensated. Similar patterns can also be found in the SMN of SZ group and BD group. Results in Figures 5A–C show that all SZ, ADHD and BD groups have the highest TRCs in VSN, which points out the neural circuit bias of visual processing ability in mental disorders (32, 43). Our findings support previous observations and reveal complementary information on the dynamic and static aspects of subnetwork-scale TRCs, underscoring the importance of including both in analyses of mental disorders.
Despite its contributions, this study has certain limitations that should be considered when interpreting the results. First, although the TRC analysis workflow has been performed at both whole-brain and subnetwork scales across different parameter settings, there are many properties of TRC still under explore. For example, in future study, the TRC analysis could be further applied to regional scale to display the hub distribution patterns (e.g., hub nodes and feeder nodes). Second, due to the limited data, we fail to correlate TRC with some clinical or cognitive scales, thus impairing the clinical interpretability. Third, the TRC analysis workflow employed in this study was mainly based on the linear correlation that one brain region may influence another, ignoring the inherently nonlinear properties of fMRI signals, future studies may benefit from examining the differences in non-linear relationships among rich-club regions.
5. Conclusion
We present the TRC analysis workflow as a comprehensive method to uncover altered rich-club patterns in major psychiatric conditions at both whole-brain and subnetwork levels. Our results reveal that SZ and BD groups exhibit significantly decreased TRC coefficients compared to their corresponding HC groups at the whole-brain scale and in most subnetworks. Conversely, the ADHD group exhibits reduced TRC coefficients in longer durations, rather than shorter duration, which significantly differs from the SZ and BD groups. Our findings highlight the potential of TRC coefficients as effective biomarkers to distinguish between psychiatric groups with similar phenotypes, and emphasize the importance of considering both static and dynamic aspects of subnetwork-scale TRCs in the study of mental disorders.
Funding
This study was supported by the Zhejiang Xinmiao Talents Program (2023R401197), the Innovation Fund of the Department of Education of Gansu Province (2022B-023), and the Fund of The First Hospital of Lanzhou University. China (Ldyyyn2021-73).
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Ethics statement
The imaging data and phenotypic information was collected and shared by the Mind Research Network and the University of New Mexico funded by a National Institute of Health Center of Biomedical Research Excellence grant 1P20RR021938-01A2 and the Consortium for Neuropsychiatric Phenomics (NIH Roadmap for Medical Research grants UL1-DE019580, RL1MH083268, RL1MH083269, RL1DA024853, RL1MH083270, RL1LM009833, PL1MH083271, and PL1NS062410). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
MN and HG contributed to conception and statistical analysis of the study. YF wrote the first draft of the manuscript, performed validation, and conducted the whole project. ZZ revised the first draft of the manuscript and performed validation. All authors contributed to the article and approved the submitted version.
Conflict of interest
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2023.1226143/full#supplementary-material
References
1.
HymanS.ChisholmD.KesslerR.PatelV.WhitefordH., (2006). In: Disease control priorities related to mental, neurological, developmental and substance abuse disorders. Oxford University Press. 1–20.
2.
Organization, W.H. World mental health report: transforming mental health for all. Regional Office for the Eastern Mediterranean: World Health Organization. (2019).
3.
YamadaYMatsumotoMIijimaKSumiyoshiT. Specificity and continuity of schizophrenia and bipolar disorder: relation to biomarkers. Curr Pharm Des. (2020) 26:191–200. doi: 10.2174/1381612825666191216153508
4.
WolfersTRokickiJAlnæsDBerthetPAgartzIKiaSMet al. Replicating extensive brain structural heterogeneity in individuals with schizophrenia and bipolar disorder. Hum Brain Mapp. (2021) 42:2546–55. doi: 10.1002/hbm.25386
5.
O’ConnellKSShadrinABahramiSSmelandOBBettellaFFreiOet al. Identification of genetic overlap and novel risk loci for attention-deficit/hyperactivity disorder and bipolar disorder. Mol Psychiatry. (2021) 26:4055–65. doi: 10.1038/s41380-019-0613-z
6.
SalviVRibuoliEServasiMOrsoliniLVolpeU. ADHD and bipolar disorder in adulthood: clinical and treatment implications. Medicina. (2021) 57:466. doi: 10.3390/medicina57050466
7.
YaoZShiJZhangZZhengWHuTLiYet al. Altered dynamic functional connectivity in weakly-connected state in major depressive disorder. Clin Neurophysiol. (2019a) 130:2096–104. doi: 10.1016/j.clinph.2019.08.009
8.
DuMZhangLLiLJiEHanXHuangGet al. Abnormal transitions of dynamic functional connectivity states in bipolar disorder: a whole-brain resting-state fMRI study. J Affect Disord. (2021) 289:7–15. doi: 10.1016/j.jad.2021.04.005
9.
FuYNiuMGaoYDongSHuangYZhangZet al. Altered nonlinear granger causality interactions in the large-scale brain networks of patients with schizophrenia. J Neural Eng. (2022a) 19:066044. doi: 10.1088/1741-2552/acabe7
10.
YaoZZouYZhengWZhangZLiYYuYet al. Structural alterations of the brain preceded functional alterations in major depressive disorder patients: evidence from multimodal connectivity. J Affect Disord. (2019b) 253:107–17. doi: 10.1016/j.jad.2019.04.064
11.
IngalhalikarMShindeSKarmarkarARajanARangaprakashDDeshpandeG. Functional connectivity-based prediction of autism on site harmonized ABIDE dataset. IEEE Trans Biomed Eng. (2021) 68:3628–37. doi: 10.1109/TBME.2021.3080259
12.
Van Den HeuvelMPSpornsO. Rich-club organization of the human connectome. J Neurosci. (2011) 31:15775–86. doi: 10.1523/JNEUROSCI.3539-11.2011
13.
CollinGSpornsOMandlRCVan Den HeuvelMP. Structural and functional aspects relating to cost and benefit of rich club organization in the human cerebral cortex. Cereb Cortex. (2014) 24:2258–67. doi: 10.1093/cercor/bht064
14.
RaySMillerMKaralunasSRobertsonCGraysonDSCaryRPet al. Structural and functional connectivity of the human brain in autism spectrum disorders and attention-deficit/hyperactivity disorder: a rich club-organization study. Hum Brain Mapp. (2014) 35:6032–48. doi: 10.1002/hbm.22603
15.
LiuXHeCFanDZhuYZangFWangQet al. Disrupted rich-club network organization and individualized identification of patients with major depressive disorder. Prog Neuro-Psychopharmacol Biol Psychiatry. (2021) 108:110074. doi: 10.1016/j.pnpbp.2020.110074
16.
FuY.HuangY.NiuM.XueL.DongS.GuoS.et al. (2022b). Active index: an integrated index to reveal disrupted brain network organizations of major depressive disorder patients. In: 2022 IEEE 19th international symposium on biomedical imaging (ISBI). Kolkata: IEEE, 1–5.
17.
WangYDengFJiaYWangJZhongSHuangHet al. Disrupted rich club organization and structural brain connectome in unmedicated bipolar disorder. Psychol Med. (2019) 49:510–8. doi: 10.1017/S0033291718001150
18.
BrahaDBar-YamY. From centrality to temporary fame: dynamic centrality in complex networks. Complexity. (2006) 12:59–63. doi: 10.1002/cplx.20156
19.
PedreschiNBattagliaDBarratA. The temporal rich club phenomenon. Nat Phys. (2022) 18:931–8. doi: 10.1038/s41567-022-01634-8
20.
HutchisonRMWomelsdorfTAllenEABandettiniPACalhounVDCorbettaMet al. Dynamic functional connectivity: promise, issues, and interpretations. NeuroImage. (2013) 80:360–78. doi: 10.1016/j.neuroimage.2013.05.079
21.
HindriksRAdhikariMHMurayamaYGanzettiMMantiniDLogothetisNKet al. Can sliding-window correlations reveal dynamic functional connectivity in resting-state fMRI?NeuroImage. (2016) 127:242–56. doi: 10.1016/j.neuroimage.2015.11.055
22.
ChengJCAnzolinABerryMHonariHPaschaliMLazaridouAet al. Dynamic functional brain connectivity underlying temporal summation of pain in fibromyalgia. Arthritis Rheumatol. (2022) 74:700–10. doi: 10.1002/art.42013
23.
De DomenicoMSasaiSArenasA. Mapping multiplex hubs in human functional brain networks. Front Neurosci. (2016) 10:326. doi: 10.3389/fnins.2016.00326
24.
XiangYWangJTanGWuF-XLiuJ. Schizophrenia identification using multi-view graph measures of functional brain networks. Front Bioeng Biotechnol. (2020) 7:479. doi: 10.3389/fbioe.2019.00479
25.
PoldrackRACongdonETriplettWGorgolewskiKKarlsgodtKMumfordJet al. A phenome-wide examination of neural and cognitive function. Sci Data. (2016) 3:1–12. doi: 10.1038/sdata.2016.110
26.
YanC-GWangX-DZuoX-NZangY-F. DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics. (2016) 14:339–51. doi: 10.1007/s12021-016-9299-4
27.
DosenbachNUNardosBCohenALFairDAPowerJDChurchJAet al. Prediction of individual brain maturity using fMRI. Science. (2010) 329:1358–61. doi: 10.1126/science.1194144
28.
YanC-GChenXLiLCastellanosFXBaiT-JBoQ-Jet al. Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. Proc Natl Acad Sci. (2019) 116:9078–83. doi: 10.1073/pnas.1900390116
29.
BertoleroMAYeoBTD’EspositoM. The diverse club. Nat Commun. (2017) 8:1277. doi: 10.1038/s41467-017-01189-w
30.
HolmePSaramäkiJ. Temporal network theory.Nature Switzerland AG: Springer (2019).
31.
FuYZhangJLiYShiJZouYGuoHet al. A novel pipeline leveraging surface-based features of small subcortical structures to classify individuals with autism spectrum disorder. Prog Neuropsychopharmacol Biol Psychiatry. (2021) 104:109989. doi: 10.1016/j.pnpbp.2020.109989
32.
GuanSWanDZhaoRCanarioEMengCBiswalBB. The complexity of spontaneous brain activity changes in schizophrenia, bipolar disorder, and ADHD was examined using different variations of entropy. Hum Brain Mapp. (2023) 44:94–118. doi: 10.1002/hbm.26129
33.
XuNHuggonBSaundersKE. Cognitive impairment in patients with bipolar disorder: impact of pharmacological treatment. CNS Drugs. (2020) 34:29–46. doi: 10.1007/s40263-019-00688-2
34.
Zarp PetersenJVaroCSkovsenCFOttCVKjærstadHLVietaEet al. Neuronal underpinnings of cognitive impairment in bipolar disorder: a large data-driven functional magnetic resonance imaging study. Bipolar Disord. (2022) 24:69–81. doi: 10.1111/bdi.13100
35.
HummerTAYungMGGoñiJConroySKFrancisMMMehdiyounNFet al. Functional network connectivity in early-stage schizophrenia. Schizophr Res. (2020) 218:107–15. doi: 10.1016/j.schres.2020.01.023
36.
JiaWZhuHNiYSuJXuRJiaHet al. Disruptions of frontoparietal control network and default mode network linking the metacognitive deficits with clinical symptoms in schizophrenia. Hum Brain Mapp. (2020) 41:1445–58. doi: 10.1002/hbm.24887
37.
MaQTangYWangFLiaoXJiangXWeiSet al. Transdiagnostic dysfunctions in brain modules across patients with schizophrenia, bipolar disorder, and major depressive disorder: a connectome-based study. Schizophr Bull. (2020) 46:699–712. doi: 10.1093/schbul/sbz111
38.
WatskyREGottsSJBermanRAMcAdamsHMZhouXGreensteinDet al. Attenuated resting-state functional connectivity in patients with childhood-and adult-onset schizophrenia. Schizophr Res. (2018) 197:219–25. doi: 10.1016/j.schres.2018.01.003
39.
BradyROJrTandonNMastersGAMargolisACohenBMKeshavanMet al. Differential brain network activity across mood states in bipolar disorder. J Affect Disord. (2017) 207:367–76. doi: 10.1016/j.jad.2016.09.041
40.
WangYGaoYTangSLuLZhangLBuXet al. Large-scale network dysfunction in the acute state compared to the remitted state of bipolar disorder: a meta-analysis of resting-state functional connectivity. EBioMedicine. (2020) 54:102742. doi: 10.1016/j.ebiom.2020.102742
41.
ZhangHZhouZDingLWuCQiuMHuangYet al. Divergent and convergent imaging markers between bipolar and unipolar depression based on machine learning. IEEE J Biomed Health Inform. (2022) 26:4100–10. doi: 10.1109/JBHI.2022.3166826
42.
CorteseS. Understanding the network bases of ADHD: an overview of the fMRI evidence In: Brain Network dysfunction in neuropsychiatric illness: methods, applications, and implications. Eds. DiwadkarV. A.EickhoffS. B.. Cham: Springer International Publishing. (2021). 333–43.
43.
JimenezAMRiedelPLeeJReavisEAGreenMF. Linking resting-state networks and social cognition in schizophrenia and bipolar disorder. Hum Brain Mapp. (2019) 40:4703–15. doi: 10.1002/hbm.24731
Summary
Keywords
major psychiatric disorders, rich-club organization, temporal properties, functional connectivity, brain networks
Citation
Niu M, Guo H, Zhang Z and Fu Y (2023) Abnormal temporal variability of rich-club organization in three major psychiatric conditions. Front. Psychiatry 14:1226143. doi: 10.3389/fpsyt.2023.1226143
Received
20 May 2023
Accepted
07 August 2023
Published
31 August 2023
Volume
14 - 2023
Edited by
Wentian Li, Wuhan Mental Health Center, China
Reviewed by
Dandan Li, Taiyuan University of Technology, China; Xinyuan Yan, University of Minnesota Twin Cities, United States
Updates
Copyright
© 2023 Niu, Guo, Zhang and Fu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhe Zhang, zhangz@hznu.edu.cnYu Fu, yufu1994@zju.edu.cn
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.