Hostname: page-component-788cddb947-r7bls Total loading time: 0 Render date: 2024-10-14T14:48:26.710Z Has data issue: false hasContentIssue false

Identification of common neural substrates with connectomic abnormalities in four major psychiatric disorders: A connectome-wide association study

Published online by Cambridge University Press:  03 December 2020

Pei-Chi Tu
Affiliation:
Department of Medical Research, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Institute of Philosophy of Mind and Cognition, National Yang-Ming University, Taipei, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan
Mu-Hong Chen
Affiliation:
Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan Institute of Brain Science, National Yang-Ming University, Taipei, Taiwan
Wan-Chen Chang
Affiliation:
Department of Medical Research, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Biomedical Engineering, National Yang-Ming University, Taipei, Taiwan
Zih-Kai Kao
Affiliation:
Department of Medical Research, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan
Ju-Wei Hsu
Affiliation:
Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan
Wei-Chen Lin
Affiliation:
Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan Institute of Brain Science, National Yang-Ming University, Taipei, Taiwan
Cheng-Ta Li
Affiliation:
Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan Institute of Brain Science, National Yang-Ming University, Taipei, Taiwan
Tung-Ping Su*
Affiliation:
Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan Institute of Brain Science, National Yang-Ming University, Taipei, Taiwan Department of Psychiatry, Cheng Hsin General Hospital, Taipei, Taiwan
Ya-Mei Bai*
Affiliation:
Department of Psychiatry, Taipei Veterans General Hospital, Taipei112, Taiwan Department of Psychiatry, Faculty of Medicine, National Yang-Ming University, Taipei, Taiwan Institute of Brain Science, National Yang-Ming University, Taipei, Taiwan
*
Tung-Ping Su, E-mail: tomsu0402@gmail.com
Ya Mei Bai, E-mail: ymbi@mail2000.com.tw

Abstract

Background

Recent imaging studies of large datasets suggested that psychiatric disorders have common biological substrates. This study aimed to identify all the common neural substrates with connectomic abnormalities across four major psychiatric disorders by using the data-driven connectome-wide association method of multivariate distance matrix regression (MDMR).

Methods

This study analyzed a resting functional magnetic resonance imaging dataset of 100 patients with schizophrenia, 100 patients with bipolar I disorder, 100 patients with bipolar II disorder, 100 patients with major depressive disorder, and 100 healthy controls (HCs). We calculated a voxel-wise 4,330 × 4,330 matrix of whole-brain functional connectivity (FC) with 8-mm isotropic resolution for each participant and then performed MDMR to identify structures where the overall multivariate pattern of FC was significantly different between each patient group and the HC group. A conjunction analysis was performed to identify common neural regions with FC abnormalities across these four psychiatric disorders.

Results

The conjunction of the MDMR maps revealed that the four groups of patients shared connectomic abnormalities in distributed cortical and subcortical structures, which included bilateral thalamus, cerebellum, frontal pole, supramarginal gyrus, postcentral gyrus, lingual gyrus, lateral occipital cortex, and parahippocampus. The follow-up analysis based on pair-wise FC of these regions demonstrated that these psychiatric disorders also shared similar patterns of FC abnormalities characterized by sensory/subcortical hyperconnectivity, association/subcortical hypoconnectivity, and sensory/association hyperconnectivity.

Conclusions

These findings suggest that major psychiatric disorders share common connectomic abnormalities in distributed cortical and subcortical regions and provide crucial support for the common network hypothesis of major psychiatric disorders.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - SA
This is an Open Access article, distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike licence (http://creativecommons.org/licenses/by-nc-sa/4.0/), which permits non-commercial reuse, distribution, and reproduction in any medium, provided the same Creative Commons licence is included and the original work is properly cited. The written permission of Cambridge University Press must be obtained for commercial re-use.
Copyright
© The Author(s), 2021. Published by Cambridge University Press on behalf of the European Psychiatry Association

Introduction

Making a diagnosis based on a cluster of clinical symptoms according to the diagnostic classification system and prescribing appropriate psychotropics play a central role in modern psychiatric practice. However, the overlap of clinical symptoms in patients with different disorders suggests that psychiatric disorders may share common neural substrates. For example, affective symptoms are common in patients with schizophrenia; depressive and mania symptoms were reported in 80 and 20% of patients with schizophrenia, respectively [Reference an der Heiden, Könnecke, Maurer, Ropeter and Häfner1]. A cross-sectional study also revealed a high prevalence of current hallucinations in hospitalized psychiatric patients with mixed bipolar depression (22.9%), bipolar manic depression (11.2%), bipolar depression (10.5%), and unipolar depression (5.9%) [Reference Baethge, Baldessarini, Freudenthal, Streeruwitz, Bauer and Bschor2]. Nevertheless, most previous structural or functional imaging studies have focused on patients with one psychiatric disorder and have rarely tested the aforementioned hypothesis directly in the same study. Meta-analysis of imaging studies across different diagnostic groups is an important approach for testing the hypothesis of the common neural substrates. A meta-analysis of 193 voxel-based morphometry studies comprising 15,892 individuals across 6 diverse diagnostic groups (schizophrenic disorder [SZ], bipolar disorder [BD], major depressive disorder [MDD], substance addiction, obsessive–compulsive disorder, and anxiety disorder) with matched controls revealed that gray matter loss converged across diagnoses in three regions: dorsal anterior cingulate cortex and bilateral insula [Reference Goodkind, Eickhoff, Oathes, Jiang, Chang and Jones-Hagata3]. Functional connectivity (FC) analysis indicated that these three regions form a tightly interconnected network called the salience network. Another meta-analysis of 283 whole-brain functional magnetic resonance imaging (fMRI) studies during cognitive control tasks [Reference McTeague, Huemer, Carreon, Jiang, Eickhoff and Etkin4] in 5,493 patients with various disorders (SZ, BD, MDD, anxiety disorders, and substance abuse disorders) found transdiagnostic abnormal activations in several frontal and parietal regions in the multidemand network and salience network. Another meta-analysis of 226 task-related functional imaging studies on mood disorders, posttraumatic stress disorder, and anxiety disorders also revealed that the most consistent transdiagnostic abnormalities in task-related brain activity converge in regions that are primarily associated with inhibitory control and salience processing [Reference Janiri, Moser, Doucet, Luber, Rasgon and Lee5]. These results provided preliminary evidences for the hypothesis of common neural substrates in major psychiatric disorders.

Recently, increasing numbers of studies adopted a transdiagnostic approach, which recruited patients with various psychiatric disorders and controls in the same study, and showed the potential to identify the common neural substrates across various psychiatric disorders directly. Several studies have successfully identified common deficits of cognitive controls [Reference McTeague, Goodkind and Etkin6], dopamine synthesis [Reference Jauhar, Nour, Veronese, Rogdaki, Bonoldi and Azis7], gray matter volume [Reference Gong, Scarpazza, Dai, He, Xu and Shi8], and network-level FC [Reference Gong, Hu, Pettersson-Yeo, Xu, Lui and Crossley9] among various psychiatric disorders. Our previous FC study focused on the thalamus and identified the common patterns of thalamocortical dysconnectivity in four major psychiatric disorders (SZ, BD-I, BD-II, and MDD) [Reference Tu, Bai, Li, Chen, Lin and Chang10]. However, it is unclear whether common FC abnormalities also exist in structures other than the thalamus. To determine all the common structures with FC abnormalities in major psychiatric disorders by using FC analyses, a brain-wise map of FC abnormalities in each disorder is required. In this study, we proposed that the recently developed connectome-wide association method of multivariate distance matrix regression (MDMR) [Reference Shehzad, Kelly, Reiss, Cameron Craddock, Emerson and McMahon11] was suitable for the purpose and able to identify all common neural substrates with FC abnormalities across major psychiatric disorders. MDMR analysis was developed for high-dimensional data sets related to genetics or imaging, and it enables researchers to relate p variables to additional m factors collected from N individuals, where p > N [Reference Zapala and Schork12]. Shehzad et al. [Reference Shehzad, Kelly, Reiss, Cameron Craddock, Emerson and McMahon11] applied MDMR to a connectome-wide association study and provided a detailed evaluation of its statistical properties. Following the initial introduction of the method, MDMR has been used to identify connectivity deficits in youths with psychosis-spectrum symptoms [Reference Satterthwaite, Vandekar, Wolf, Bassett, Ruparel and Shehzad13], individuals at a clinically high risk of psychosis [Reference Colibazzi, Yang, Horga, Chao-Gan, Corcoran and Klahr14], and those with a transdiagnostic risk of mental illness [Reference Elliott, Romer, Knodt and Hariri15]. This method was also used to identify the FC features in different groups according to data-driven phenotypic categories [Reference Van Dam, O'Connor, Marcelle, Ho, Cameron Craddock and Tobe16]. In these studies, a pair-wise FC matrix of each voxel with variable resolutions was first calculated, and MDMR was then performed to identify voxels whose FC patterns with all other voxels were significantly associated with a clinical phenotype of interest. Subsequently, a follow-up seed-based FC analysis of the structures identified by using MDMR was usually adopted to understand the nature of FC abnormalities and to compare them with those in the literature. Importantly, as the analysis is based on the statistical inference of voxels and typically provides a complete brain-wise map of connectomic abnormalities for the effect of interest, neural substrates with FC abnormalities across different psychiatric disorders can be identified.

Here, this study aimed to demonstrate a brain-wise map of common neural structures with connectomic abnormalities in four major psychiatric disorders, namely SZ, BD-I, BD-II, and MDD, through a data-driven connectome-wide association study. We focused on these four psychiatric disorders because FC abnormalities in these disorders were the major research interests of our laboratory and we had collected a single-site resting fMRI dataset of these disorders with enough sample sizes for MDMR analysis. First, we performed MDMR to identify regions where the overall multivariate pattern of FC was significantly different between the patients of these four groups and controls. A conjunction analysis was performed to identify common neural structures with connectomic abnormalities in the four groups of patients. Then, we focused on the pair-wise FCs of these identified regions and conducted a follow-up analysis to understand the major pattern of FC abnormalities. The results of this study may provide a brain-wise map of common connectomic abnormalities in the four major psychiatric disorders and further evidence for the common network hypothesis of major psychiatric illness.

Methods

Participants

This study included 100 patients with SZ, 100 patients with BD-I, 100 patients with BD-II, 100 patients with MDD, and 100 healthy controls (HCs) (Table 1). The diagnoses were made according to the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition, Text Revision. These participants were selected from a single-site resting fMRI dataset (see Supplementary Table 1 for demographic properties of the dataset) to match the age and sex of the five groups as much as possible. The mean ages of the patients with SZ, BD-1, BD-II, MDD, and HC were 36.0, 38.3, 37.5, 37.7, 36.4, and 36.4 years, respectively, and no significant difference was noted among these five groups. However, the BD-II group in the dataset was the smallest, and hence, our BD-II sample could not be matched to the other groups by sex. The mean illness durations of the patients with SZ, BD-1, BD-II, and MDD were 12.5, 12.4, 12.1, and 7.9 years, respectively. The SZ, BD-I, and BD-II groups were matched by the illness duration. The MDD group has shorter illness duration than the other three groups because of the late onset. All participants were recruited from outpatient and inpatient units of Taipei Veterans General Hospital in Taiwan. The participants with the following conditions were excluded: (a) substance abuse or dependence over the preceding 6 months; (b) a history of head injuries that resulted in sustained loss of consciousness, cognitive sequelae, or both; and (c) neurological illnesses or any other disorder that affects cerebral metabolism. HCs were recruited through advertisements. An experienced psychiatrist used the MINI for screening to exclude candidates with major psychiatric illnesses. Furthermore, healthy candidates whose first-degree relatives had axis-I disorders, including SZ, MDD, and BD, were excluded. The clinical status of the patients with SZ was characterized using the Positive and Negative Syndrome Scale (PANSS). The clinical assessments used for the BD-I, BD-II, and MDD patients were the Young Mania Rating Scale (YMRS) and Montgomery Åsberg Depression Rating Scale (MADRS) (Table 1). However, some of patients with BD-1 and BD-II only used self-rated visual analog mood scales to rate their mood states, and their MDRS and YMRS ratings were not available. These patients were receiving treatment with various atypical antipsychotics, antidepressants, and mood stabilizers before participating in the experiment. All procedures were approved by the Institutional Review Board of Taipei Veterans General Hospital in accordance with the Declaration of Helsinki, and all participants provided written informed consent after they were informed of the experimental procedures.

Table 1. The enrolled participants’ demographic data.

Abbreviations: BDI, Beck Depression Inventory; HAMD-17, 17-item Hamilton Depression Scales; MADRS, Montgomery Åsberg Depression Rating Scale; PANSS, Positive and Negative Syndrome Scale for Schizophrenia; WM, working memory; YMRS, Young Mania Rating Scale.

Magnetic resonance image acquisition

Magnetic resonance images were acquired using a 3.0-Tesla GE Discovery 750 whole-body high-speed imaging device with an 8-channel high-resolution brain coil. Head stabilization was achieved with cushioning, and all the participants wore earplugs (29-dB rating) to attenuate noise. Automated shimming procedures were performed, and scout images were obtained. The resting-state functional images were collected using a gradient echo T2* weighted sequence (repetition time [TR]/echo time [TE]/Flip = 2,500 ms/30 ms/90°). Forty-seven contiguous horizontal slices parallel to the intercommissural plane (voxel size: 3.5 × 3.5 × 3.5 mm) were acquired and interleaved. These slices covered the cerebellum of each participant. During the functional scans, the participants were instructed to keep their eyes open (each scan lasted for 8 min and 24 s across 200 time points). In addition, a high-resolution structural image was acquired in the sagittal plane using a high-resolution sequence (TR = 2,530 ms, echo spacing = 7.25 ms, TE = 3 ms, and flip angle = 7°) and an isotropic 1 mm voxel (FOV 256 × 256).

Quality check of functional images

Regarding head motion during image acquisition, we used the method of scrubbing within regression (spike regression) suggested in a previous study [Reference Yan, Cheung, Kelly, Colcombe, Craddock and Di Martino17] to minimize the effect of head motion on the FC measurement. This method identifies “bad” time points using a threshold of framewise displacement > 0.2 mm as well as 1 back and 2 forward neighbors (as reported by [Reference Power, Barnes, Snyder, Schlaggar and Petersen18] and then models each “bad” time point as a separate regressor in the regression models [Reference Lemieux, Salek-Haddadi, Lund, Laufs and Carmichael19, Reference Satterthwaite, Elliott, Gerraty, Ruparel, Loughead and Calkins20]. We also calculated the number of contaminated volumes in each group, and no significant difference was observed between the groups (Supplementary Table 2).

FC preprocessing

All preprocessing was performed using the Data Processing Assistant for Resting-State fMRI (http://www.restfmri.net), which is based on Statistical Parametric Mapping (http://www.fil.ion.ucl.ac.uk/spm) and the Resting-State fMRI Data Analysis Toolkit (http://www.restfmri.net). The functional scans received slice-timing correction and motion correction and were normalized to a standard anatomical space (Montreal Neurological Institute). To prepare the data for FC analysis, the following additional preprocessing steps were used: (a) spatial smoothing by using a Gaussian kernel (6-mm full width at half-maximum), (b) temporal filtering (0.009 Hz < f < 0.08 Hz), and (c) removal of spurious or nonspecific sources of variance through the regression of the following variables: (a) Six head motion parameters and autoregressive models of motion: Six head motion parameters, Six head motion parameters one time point before, and the 12 corresponding squared items [Reference Friston, Williams, Howard, Frackowiak and Turner21] (Friston 24-parameter model); (b) the mean whole-brain signal; (c) the mean signal within the lateral ventricles; and (d) the mean signal within a white matter mask. Furthermore, the regressors used in the method of scrubbing within regression were included to minimize the effect of head motion on the FC measurement. The regression of each of these signals was computed simultaneously, and the residual time course was then retained for the correlation analysis.

Connectome-wide association analysis by using MDMR

We created a brain mask including all cortical and subcortical structures of the Montreal Neurological Institute (MNI152) atlas, except regions in the ventricles. The brain mask was subsampled twice into 4,330 isotropic voxels with 8-mm resolutions to decrease the computation load. The connectome-wide association analysis for each patient group was conducted in three steps [Reference Shehzad, Kelly, Reiss, Cameron Craddock, Emerson and McMahon11]. First, we calculated a voxel-wise whole-brain 4,330 × 4,330 FC matrix based on Pearson’s correlation coefficient of low-frequency fMRI fluctuations with Fisher’s r-to-z transformation [Reference Vincent, Snyder, Fox, Shannon, Andrews and Raichle22]. In the second step, for each voxel in the connectome, we calculated the similarity between the connectivity maps for all possible pairings of the participants by using a spatial correlation and converted them to a distance metric$ \sqrt{2\ast \left(1-r\right)} $, yielding a 200 × 200 matrix (100 patients and 100 controls). Finally, MDMR was used to test how well the diagnosis of interest explained the relation between each participant’s pattern of connectivity, with age, sex, and education level added as the covariates of no interest. For each voxel’s connectivity pattern, MDMR yielded a pseudo-F statistic whose significance was assessed using 100,000 iterations of a permutation test. The threshold was set at p < 0.01 (corrected for false discovery rate). Any voxel exceeding the threshold was considered significant. An MDMR analysis was separately performed for each diagnostic group of patients. Then, we adopted a conjunction analysis to identify common areas that showed significant FC abnormalities across the four major psychiatric disorders. The clusters in conjunction analysis were selected as regions of interest (ROIs) for follow-up analysis.

Follow-up analysis

Although MDMR identified clusters where the overall multivariate pattern of connectivity is dimensionally related to SZ diagnosis, it did not provide the specific pattern of connectivity driving the significant result. We conducted post hoc ROI-based analyses from the clusters through conjunction analysis. The pair-wise FC values were calculated and compared between each group of patients and HCs by using independent sample t-tests, with age, sex, and education level as the covariates of no interest. The FC abnormalities with a significance of p < 0.05 (uncorrected) were demonstrated to visualize the major patterns of FC abnormalities. Importantly, the follow-up analysis was performed to show the major patterns of FC abnormalities between these ROIs and did not focus on the statistical significance of the specific FCs between ROI pairs.

Result

Demographic results

The demographic data and clinical ratings of the participants are shown in Table 1. Except for the BD-II group, all the other groups were age and sex matched. The BD-II group had more women. In all the four groups, the patients had a low-education level.

MDMR analysis

Compared with HC, the four groups of patients showed FC abnormalities in distributed cortical and subcortical regions (Figure 1). The conjunction analysis identified the common neural substrates showing connectomic abnormalities in all the four patient groups (Figure 1, Table 2). The structures with the largest cluster size were the left and right postcentral gyrus. Other sensorimotor-related cortices included precentral gyrus, lingual gyrus, lateral occipital cortex, and fusiform cortex. Furthermore, large clusters in the subcortical structures of the cerebellum and thalamus were also found. At last, the common structures also included several smaller clusters in fronto-parietal association cortices, including bilateral frontal pole, supramarginal gyrus, precuneus, and posterior cingulate. We also quantified the numbers of voxels showing FC abnormalities in each diagnostic group and observed severity in the order of SZ > BD-I > BD-II > MDD (Figure 1B).

Figure 1. (A) The neural substrates with connectomic abnormalities in each diagnosis were identified by MDMR. Then, the common structures with connectomic abnormalities among the four diagnostic groups were identified through conjunction analysis and included the thalamus, cerebellum, postcentral gyrus, and association cortices in frontal and parietal regions. (B) The total numbers of voxels with connectomic abnormalities in each diagnostic group were calculated and the results showed a graded severity of schizophrenia (SZ) > bipolar I disorder (BD-1) > bipolar II disorder (BD-2) > MDD.

Table 2. Conjunction analysis of functional connectivity patterns for four groups.

Abbreviations: L, left; R, right.

Follow-up analysis

The results of follow-up analyses based on pair-wise FCs of these structures were provided in Figure 2. We grouped ROIs according to the functional classification of subcortical structures, sensory-related cortices, and association cortices. The 4 groups of patients showed a similar pattern of FC abnormalities characterized by sensory/subcortical hyperconnectivity, association/subcortical hypoconnectivity, and sensory/association hyperconnectivity. In addition to the three major dysconnectivity patterns, hypoconnectivities among the structures belonging to the same category were also observed. The intracalcarine gyrus was an exception and showed hyperconnectivities with other sensory-related cortices despite its role in visual processing.

Figure 2. The result of the follow-up analysis based on pair-wise FCs of the common structures identified using MDMR. The order of structures in the figure was arranged according to the classification of subcortical structures, sensory-related cortices, and association cortices. The four groups of patients shared similar patterns of FC abnormalities, characterized by sensory/subcortical hyperconnectivity, association/subcortical hypoconnectivity, and sensory/association hyperconnectivity. Furthermore, a decrease in FC was observed within each structure category were noted.

Control analysis

Because shorter duration of illness in patients with MDD was a major confounder in this study, we performed a control analysis to evaluate the effect of illness duration on our findings. We selected 100 patients with MDD (age: 41.8 ± 12.5 years; duration of illness: 11.4 ± 8.5 years) with a similar duration of illness as the other three diagnostic groups from our resting fMRI dataset and performed the same MDMR and conjunction analysis. The control analysis revealed slightly more voxels with connectomic abnormalities and additional clusters in the bilateral superior temporal gyrus (see Supplementary Table 3 for detailed structures and cluster size); however, the pattern of graded severities (Figure 1B) and the characteristic patterns of dysconnectivities (Figure 2) did not change significantly.

Discussion

What are the core neural substrates with FC abnormalities in major psychiatric disorders? And what is the major pattern of FC abnormalities? This study attempted to answer the questions through a large-scale connectome-wide association study. We used a fully data-driven survey of the functional connectome to identify common neural substrates with connectomic abnormalities across the four major psychiatric disorders. These structures included the thalamus, cerebellum, sensorimotor-related cortices, and association cortices in the frontal and parietal regions. The follow-up analysis demonstrated that the four major psychiatric disorders also shared a similar pattern of FC abnormalities among these structures; the abnormalities were characterized by sensory/subcortical hyperconnectivity, association/subcortical hypoconnectivity, and sensory/association hyperconnectivity. The hypoconnectivity among structures within the same category was also observed. The results extended our previous findings of common FC abnormalities in thalamus in four major psychiatric disorders and successfully delineated all the common neural substrates with connectomic abnormalities across these major psychiatric disorders.

Cognitive dysmetria theory is the first to theoretically emphasize the role of cortico-thalamic-cerebellar-cortical (CCTC) circuits in the pathogenesis of patients with SZ [Reference Andreasen, Paradiso and O'Leary23]. Our finding that the thalamus and cerebellum are the common subcortical structures involved across psychiatric disorders is in accordance with the theory and suggests that the CCTC circuit deficits are shared across the four major psychiatric disorders we investigated. Moreover, intrinsic abnormalities in the cerebello-thalamo-cortical circuitry were also observed in individuals at a clinically high risk of psychosis in a previous study based on the North American Prodrome Longitudinal Study [Reference Bernard, Orr and Mittal24]. In our follow-up analysis, these two subcortical structures showed a similar pattern of bidirectional FC changes, marked by hyperconnectivity with sensory-related cortices and hypoconnectivity with the association cortex. The bidirectional changes were first found in FC studies of thalamocortical connectivity in SZ [Reference Woodward, Karbasforoushan and Heckers25, Reference Anticevic, Cole, Repovs, Murray, Brumbaugh and Winkler26], but also applied to other subcortical structures of basal ganglion and cerebellum in a recent study [Reference Ji, Diehl, Schleifer, Tamminga, Keshavan and Sweeney27]. Although the subcortical structures of the putamen, hippocampus, and midbrain were not identified in our conjunction analysis, MDMR revealed that they exhibited connectomic abnormalities in those with SZ or BD; these three structures also shared the same pattern of bidirectional FC changes. These findings suggested that that the patterns of sensory/subcortical hyperconnectivity and association/subcortical hypoconnectivity are fundamental to FC abnormalities of subcortical structures in major psychiatric disorders.

The largest cortical cluster of the common neural substrate is the postcentral gyrus, which extends to the precentral gyrus and involves the whole primary sensorimotor cortex. In addition to primary sensorimotor cortices, the MDMR analysis identified common neural substrates in other sensory-related areas, including structures in the occipital and temporal cortex for primary and higher visual processing. The follow-up analysis demonstrated that the two major patterns of FC abnormalities were sensory/subcortical and sensory/association hyperconnectivity. Our findings that four major psychiatric disorders shared connectomic abnormalities in primary sensorimotor cortex and other sensory-related cortices suggested the important role of sensorimotor processing in pathogenesis of major psychiatric disorders. It is in accordance with a recent transdiagnostic study that associated three latent components of psychopathology with connectivity alterations within the somatosensory-motor network and its connectivity with subcortical structures and association cortices in several cortical executive networks [Reference Kebets, Holmes, Orban, Tang, Li and Sun28]. Functional deficits in visual and primary sensorimotor processing received more clinical attentions recently. The deficits at different levels of visual processing were found in patients with SZ [Reference Kogata and Iidaka29] and BD-I [Reference Jahshan, Wynn, McCleery, Glahn, Altshuler and Green30, Reference Fernandes, Silverstein, Almeida and Santos31]. The excitability deficit in primary sensorimotor deficit measured using transcranial magnetic stimulation has been well studied, and a review suggested a general alteration in motor cortical inhibition in mental illness, rather than disease-specific changes [Reference Bunse, Wobrock, Strube, Padberg, Palm and Falkai32]. Although the auditory-related superior temporal gyrus reached significance only in SZ and BD-I in the MDMR analysis and therefore was not identified in the conjunction analysis, they also shared a similar pattern of FC abnormalities as the primary sensorimotor cortex.

Despite smaller cluster size, our analysis also identified several fronto-parietal structures, majorly in lateral frontal and parietal regions. The FC deficits in the association cortices of frontal and parietal regions have received special interest in the past 10 years in FC studies of SZ and other psychiatric disorders because of their important roles in cognitive controls. One crucial functional classification of the association cortices is the distinction between task-negative and task-positive networks. The task-negative network was less active during various cognitive-demanding tasks and consisted of the default mode network [Reference Greicius, Srivastava, Reiss and Menon33], whereas task-positive networks were more active during various cognitive-demanding tasks and majorly included the dorsal attention network [Reference Fox, Snyder, Vincent, Corbetta, Van Essen and Raichle34], salience network, and executive control network [Reference Dosenbach, Fair, Miezin, Cohen, Wenger and Dosenbach35, Reference Seeley, Menon, Schatzberg, Keller, Glover and Kenna36]. Among them, the FCs of the task-negative default mode network initially received great research interest and has been widely studied in major psychiatric disorders; however, the results have not been completely consistent. Hyperconnectivity [Reference Whitfield-Gabrieli, Thermenos, Milanovic, Tsuang, Faraone and McCarley37] and hypoconnectivity [Reference Bluhm, Miller, Lanius, Osuch, Boksman and Neufeld38, Reference Camchong, MacDonald, Bell, Mueller and Lim39] have both been reported in FC studies of schizophrenia. By contrast, numerous studies have reported that task-positive control networks have FC impairments [Reference White, Joseph, Francis and Liddle40Reference Manoliu, Meng, Brandl, Doll, Tahmasian and Scherr43]. In this study, the structures in association cortices identified in our MDMR analysis included major clusters in the lateral frontal and parietal cortices and a small cluster in the precuneus. As the lateral frontal parietal cortices were functionally connected to the task-positive salience or executive control networks, our findings suggest that the common neural substrates in major psychiatric disorders may involve the task-positive networks. It is consistent the two meta-analyses that suggested structural deficits in the salience network and functional activation deficits in task-positive multi-demand network in major psychiatric disorders [Reference Goodkind, Eickhoff, Oathes, Jiang, Chang and Jones-Hagata3, Reference McTeague, Huemer, Carreon, Jiang, Eickhoff and Etkin4].

Our findings have several methodological limitations and alternative conceptualizations. First, the illness durations were not matched in MDD patients. Due to a late illness onset in the MDD group, controlling the age and illness duration simultaneously was impossible. We performed a control analysis by selecting the MDD groups with a similar duration of illness and the main findings were not significantly changed. Second, our patient groups were exposed to various psychotropic medications. Since previous studies have demonstrated the effects of antipsychotics [Reference Kraguljac, White, Hadley, Hadley, Ver Hoef and Davis44, Reference Wang, Tang, Fan, Zhang, Geng and Jiang45], anti-depressants [Reference Wang, Xia, Li, Zeng, Su and Dai46Reference Cullen, Klimes-Dougan, Vu, Westlund Schreiner, Mueller and Eberly48], and mood stabilizers [Reference Pavuluri, Ellis, Wegbreit, Passarotti and Stevens49] on brain FCs, the use of various psychotropics in participants with psychiatric disorders may be confounding our findings in this study. Third, not all of the patients with BD-1 and BD-II received full evaluations of clinical severity by using the MADRS and YMRS. Although these patients were all in euthymic states clinically, the lack of complete evaluations may have limited our interpretations in this study. Lastly, our follow-up analysis focused only on the structures identified using the MDMR analysis. Shehzed et al. [Reference Shehzad, Kelly, Reiss, Cameron Craddock, Emerson and McMahon11] compared the method of MDMR and traditional seed-based FC analysis and demonstrated an overlap, but the findings of these two methods were not in total agreement. Traditional univariate FC analysis may reveal other connectomic abnormalities that were not identified through the MDMR analysis in this study.

Conclusion

This data-driven connectome-wide analysis successfully identifies the common neural substrates with connectomic abnormalites, majorly in thalamus, cerebellum, sensorimotor-related cortices, and fronto-parietal association cortices in four major psychiatric disorders. These disorders also share a similar pattern of dysconnectiviies among these identified neural structures. These findings provide considerable support for the common networks hypothesis of major psychiatric disorders.

Acknowledgments

We gratefully thank all the participants who took part in this research and all the research assistants and staff who facilitated their involvement.

Financial Support

The study was supported by grants from Taipei Veterans General Hospital (V99C1-040, V101C1-159, V104C-039, V105C-119, V106C-091, and V107C-100) and Taiwan Ministry of Science and Technology (NSC99-2628-B-010-021-MY2, MOST105-2314-B-075-056-MY2, MOST103-2314-B-075-065-MY2, and MOST108-2314-B-075-038). The funders of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

Conflict of Interest

The authors declare no competing interests.

Data Availability Statement

The data that support the findings of this study are available upon request, Prof. Pei-Chi Tu ().

Supplementary Materials

To view supplementary material for this article, please visit http://dx.doi.org/10.1192/j.eurpsy.2020.106.

References

an der Heiden, W, Könnecke, R, Maurer, K, Ropeter, D, Häfner, H. Depression in the long-term course of schizophrenia. Eur Arch Psychiatry Clin Neurosci. 2005;255(3):174184.CrossRefGoogle ScholarPubMed
Baethge, C, Baldessarini, RJ, Freudenthal, K, Streeruwitz, A, Bauer, M, Bschor, T. Hallucinations in bipolar disorder: characteristics and comparison to unipolar depression and schizophrenia. Bipolar Disord. 2005;7(2):136145.CrossRefGoogle Scholar
Goodkind, M, Eickhoff, SB, Oathes, DJ, Jiang, Y, Chang, A, Jones-Hagata, LB, et al. Identification of a common neurobiological substrate for mental illness. JAMA Psychiatry. 2015;72(4):305315.CrossRefGoogle ScholarPubMed
McTeague, LM, Huemer, J, Carreon, DM, Jiang, Y, Eickhoff, SB, Etkin, A. Identification of common neural circuit disruptions in cognitive control across psychiatric disorders. Am J Psychiatry. 2017;174(7):676685.Google ScholarPubMed
Janiri, D, Moser, DA, Doucet, GE, Luber, MJ, Rasgon, A, Lee, WH, et al. Shared neural phenotypes for mood and anxiety disorders: a meta-analysis of 226 task-related functional imaging studies. JAMA Psychiatry. 2020;77(2):172179.Google ScholarPubMed
McTeague, LM, Goodkind, MS, Etkin, A. Transdiagnostic impairment of cognitive control in mental illness. J Psychiatr Res. 2016;83:3746.CrossRefGoogle ScholarPubMed
Jauhar, S, Nour, MM, Veronese, M, Rogdaki, M, Bonoldi, I, Azis, M, et al. A test of the transdiagnostic dopamine hypothesis of psychosis using positron emission tomographic imaging in bipolar affective disorder and schizophrenia. JAMA Psychiatry. 2017;74(12):12061213.Google Scholar
Gong, Q, Scarpazza, C, Dai, J, He, M, Xu, X, Shi, Y, et al. A transdiagnostic neuroanatomical signature of psychiatric illness. Neuropsychopharmacology. 2019;44(5):869875.CrossRefGoogle ScholarPubMed
Gong, Q, Hu, X, Pettersson-Yeo, W, Xu, X, Lui, S, Crossley, N, et al. Network-level dysconnectivity in drug-naïve first-episode psychosis: dissociating transdiagnostic and diagnosis-specific alterations. Neuropsychopharmacol: Off Publ Am College Neuropsychopharmacol. 2017;42(4):933940.CrossRefGoogle ScholarPubMed
Tu, PC, Bai, YM, Li, CT, Chen, MH, Lin, WC, Chang, WC, et al. Identification of common thalamocortical dysconnectivity in four major psychiatric disorders. Schizophr Bull. 2019;45(5):11431151.CrossRefGoogle ScholarPubMed
Shehzad, Z, Kelly, C, Reiss, PT, Cameron Craddock, R, Emerson, JW, McMahon, K, et al. A multivariate distance-based analytic framework for connectome-wide association studies. Neuroimage. 2014;93(Pt 1):7494.CrossRefGoogle ScholarPubMed
Zapala, MA, Schork, NJ. Statistical properties of multivariate distance matrix regression for high-dimensional data analysis. Front Genet. 2012;3:190.CrossRefGoogle ScholarPubMed
Satterthwaite, TD, Vandekar, SN, Wolf, DH, Bassett, DS, Ruparel, K, Shehzad, Z, et al. Connectome-wide network analysis of youth with psychosis-spectrum symptoms. Mol Psychiatry. 2015;20(12):15081515.CrossRefGoogle ScholarPubMed
Colibazzi, T, Yang, Z, Horga, G, Chao-Gan, Y, Corcoran, CM, Klahr, K, et al. Aberrant temporal connectivity in persons at clinical high risk for psychosis. Biol Psychiatry Cogn Neurosci Neuroimaging. 2017;2(8):696705.Google ScholarPubMed
Elliott, ML, Romer, A, Knodt, AR, Hariri, AR. A connectome-wide functional signature of transdiagnostic risk for mental illness. Biol Psychiatry. 2018;84(6):452459.CrossRefGoogle ScholarPubMed
Van Dam, NT, O'Connor, D, Marcelle, ET, Ho, EJ, Cameron Craddock, R, Tobe, RH, et al. Data-driven phenotypic categorization for neurobiological analyses: beyond DSM-5 labels. Biol Psychiatry. 2017;81(6):484494.CrossRefGoogle ScholarPubMed
Yan, C-G, Cheung, B, Kelly, C, Colcombe, S, Craddock, RC, Di Martino, A, et al. A comprehensive assessment of regional variation in the impact of head micromovements on functional connectomics. NeuroImage. 2013;76:183201.CrossRefGoogle ScholarPubMed
Power, JD, Barnes, KA, Snyder, AZ, Schlaggar, BL, Petersen, SE. Steps toward optimizing motion artifact removal in functional connectivity MRI; a reply to Carp. NeuroImage. 2013;76:439441. doi: 10.1016/j.neuroimage.2012.03.017.CrossRefGoogle ScholarPubMed
Lemieux, L, Salek-Haddadi, A, Lund, TE, Laufs, H, Carmichael, D. Modelling large motion events in fMRI studies of patients with epilepsy. Magn Reson Imaging. 2007;25(6):894901.Google ScholarPubMed
Satterthwaite, TD, Elliott, MA, Gerraty, RT, Ruparel, K, Loughead, J, Calkins, ME, et al. An improved framework for confound regression and filtering for control of motion artifact in the preprocessing of resting-state functional connectivity data. NeuroImage. 2013;64:240256. doi: 10.1016/j.neuroimage.2012.08.052.CrossRefGoogle ScholarPubMed
Friston, KJ, Williams, S, Howard, R, Frackowiak, RS, Turner, R. Movement-related effects in fMRI time-series. Magn Reson Med. 1996;35(3):346355.CrossRefGoogle ScholarPubMed
Vincent, JL, Snyder, AZ, Fox, MD, Shannon, BJ, Andrews, JR, Raichle, ME, et al. Coherent spontaneous activity identifies a hippocampal-parietal memory network. J Neurophysiol. 2006;96(6):35173531.CrossRefGoogle ScholarPubMed
Andreasen, NC, Paradiso, S, O'Leary, DS. “Cognitive dysmetria” as an integrative theory of schizophrenia: a dysfunction in cortical-subcortical-cerebellar circuitry? Schizophr Bull. 1998;24(2):203218.Google ScholarPubMed
Bernard, JA, Orr, JM, Mittal, VA. Cerebello-thalamo-cortical networks predict positive symptom progression in individuals at ultra-high risk for psychosis. NeuroImage Clin. 2017;14:622628.CrossRefGoogle ScholarPubMed
Woodward, ND, Karbasforoushan, H, Heckers, S. Thalamocortical dysconnectivity in schizophrenia. Am J Psychiatry. 2012;169(10):10921099.CrossRefGoogle Scholar
Anticevic, A, Cole, MW, Repovs, G, Murray, JD, Brumbaugh, MS, Winkler, AM, et al. Characterizing thalamo-cortical disturbances in schizophrenia and bipolar illness. Cereb Cortex. 2014;24(12):31163130.CrossRefGoogle ScholarPubMed
Ji, JL, Diehl, C, Schleifer, C, Tamminga, CA, Keshavan, MS, Sweeney, JA, et al. Schizophrenia exhibits bi-directional brain-wide alterations in cortico-striato-cerebellar circuits. Cereb Cortex. 2019;29(11):44634487.CrossRefGoogle ScholarPubMed
Kebets, V, Holmes, AJ, Orban, C, Tang, S, Li, J, Sun, N, et al. Somatosensory-motor dysconnectivity spans multiple transdiagnostic dimensions of psychopathology. Biol Psychiatry. 2019;86(10):779791.Google ScholarPubMed
Kogata, T, Iidaka, T. A review of impaired visual processing and the daily visual world in patients with schizophrenia. Nagoya J Med Sci. 2018;80(3):317328.Google ScholarPubMed
Jahshan, C, Wynn, JK, McCleery, A, Glahn, DC, Altshuler, LL, Green, MF. Cross-diagnostic comparison of visual processing in bipolar disorder and schizophrenia. J Psychiatr Res. 2014;51:4248.CrossRefGoogle Scholar
Fernandes, TP, Silverstein, SM, Almeida, NL, Santos, NA. Visual impairments in type 1 bipolar disorder. World J Biol Psychiatry. 2019;19.Google ScholarPubMed
Bunse, T, Wobrock, T, Strube, W, Padberg, F, Palm, U, Falkai, P, et al. Motor cortical excitability assessed by transcranial magnetic stimulation in psychiatric disorders: a systematic review. Brain Stimul. 2014;7(2):158169.CrossRefGoogle ScholarPubMed
Greicius, MD, Srivastava, G, Reiss, AL, Menon, V. Default-mode network activity distinguishes Alzheimer’s disease from healthy aging: evidence from functional MRI. Proc Natl Acad Sci USA. 2004;101(13):46374642.CrossRefGoogle ScholarPubMed
Fox, MD, Snyder, AZ, Vincent, JL, Corbetta, M, Van Essen, DC, Raichle, ME. The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci USA. 2005;102(27):96739678.CrossRefGoogle ScholarPubMed
Dosenbach, NU, Fair, DA, Miezin, FM, Cohen, AL, Wenger, KK, Dosenbach, RA, et al. Distinct brain networks for adaptive and stable task control in humans. Proc Natl Acad Sci USA. 2007;104(26):1107311078.CrossRefGoogle ScholarPubMed
Seeley, WW, Menon, V, Schatzberg, AF, Keller, J, Glover, GH, Kenna, H, et al. Dissociable intrinsic connectivity networks for salience processing and executive control. J Neurosci: Off J Soc Neurosci. 2007;27(9):23492356.CrossRefGoogle ScholarPubMed
Whitfield-Gabrieli, S, Thermenos, HW, Milanovic, S, Tsuang, MT, Faraone, SV, McCarley, RW, et al. Hyperactivity and hyperconnectivity of the default network in schizophrenia and in first-degree relatives of persons with schizophrenia. Proc Natl Acad Sci USA. 2009;106(4):12791284.CrossRefGoogle ScholarPubMed
Bluhm, RL, Miller, J, Lanius, RA, Osuch, EA, Boksman, K, Neufeld, RW, et al. Spontaneous low-frequency fluctuations in the BOLD signal in schizophrenic patients: anomalies in the default network. Schizophr Bull. 2007;33(4):10041012.CrossRefGoogle ScholarPubMed
Camchong, J, MacDonald, AW III, Bell, C, Mueller, BA, Lim, KO. Altered functional and anatomical connectivity in schizophrenia. Schizophr Bull. 2011;37(3):640650.CrossRefGoogle Scholar
White, TP, Joseph, V, Francis, ST, Liddle, PF. Aberrant salience network (bilateral insula and anterior cingulate cortex) connectivity during information processing in schizophrenia. Schizophr Res. 2010;123(2):105115.CrossRefGoogle Scholar
Pu, W, Li, L, Zhang, H, Ouyang, X, Liu, H, Zhao, J, et al. Morphological and functional abnormalities of salience network in the early-stage of paranoid schizophrenia. Schizophr Res. 2012;141(1):1521.CrossRefGoogle ScholarPubMed
Baker, JT, Holmes, AJ, Masters, GA, Yeo, BTT, Krienen, F, Buckner, RL, et al. Disruption of cortical association networks in schizophrenia and psychotic bipolar disorder. JAMA Psychiatry. 2014;71(2):109118.Google ScholarPubMed
Manoliu, A, Meng, C, Brandl, F, Doll, A, Tahmasian, M, Scherr, M, et al. Insular dysfunction within the salience network is associated with severity of symptoms and aberrant inter-network connectivity in major depressive disorder. Front Human Neurosci. 2014;7:930.CrossRefGoogle ScholarPubMed
Kraguljac, NV, White, DM, Hadley, N, Hadley, JA, Ver Hoef, L, Davis, E, et al. Aberrant hippocampal connectivity in unmedicated patients with schizophrenia and effects of antipsychotic medication: a longitudinal resting state functional MRI study. Schizophr Bull. 2016;42(4):10461055.CrossRefGoogle ScholarPubMed
Wang, Y, Tang, W, Fan, X, Zhang, J, Geng, D, Jiang, K, et al. Resting-state functional connectivity changes within the default mode network and the salience network after antipsychotic treatment in early-phase schizophrenia. Neuropsychiatr Dis Treat. 2017;13:397406.CrossRefGoogle ScholarPubMed
Wang, L, Xia, M, Li, K, Zeng, Y, Su, Y, Dai, W, et al. The effects of antidepressant treatment on resting-state functional brain networks in patients with major depressive disorder. Hum Brain Mapp. 2015;36(2):768778.CrossRefGoogle ScholarPubMed
Yang, R, Gao, C, Wu, X, Yang, J, Li, S, Cheng, H. Decreased functional connectivity to posterior cingulate cortex in major depressive disorder. Psychiatry Res. 2016;255:1523.CrossRefGoogle ScholarPubMed
Cullen, KR, Klimes-Dougan, B, Vu, DP, Westlund Schreiner, M, Mueller, BA, Eberly, LE, et al. Neural correlates of antidepressant treatment response in adolescents with major depressive disorder. J Child Adolesc Psychopharmacol. 2016;26(8):705712.CrossRefGoogle ScholarPubMed
Pavuluri, MN, Ellis, JA, Wegbreit, E, Passarotti, AM, Stevens, MC. Pharmacotherapy impacts functional connectivity among affective circuits during response inhibition in pediatric mania. Behav Brain Res. 2012;226(2):493503.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. The enrolled participants’ demographic data.

Figure 1

Figure 1. (A) The neural substrates with connectomic abnormalities in each diagnosis were identified by MDMR. Then, the common structures with connectomic abnormalities among the four diagnostic groups were identified through conjunction analysis and included the thalamus, cerebellum, postcentral gyrus, and association cortices in frontal and parietal regions. (B) The total numbers of voxels with connectomic abnormalities in each diagnostic group were calculated and the results showed a graded severity of schizophrenia (SZ) > bipolar I disorder (BD-1) > bipolar II disorder (BD-2) > MDD.

Figure 2

Table 2. Conjunction analysis of functional connectivity patterns for four groups.

Figure 3

Figure 2. The result of the follow-up analysis based on pair-wise FCs of the common structures identified using MDMR. The order of structures in the figure was arranged according to the classification of subcortical structures, sensory-related cortices, and association cortices. The four groups of patients shared similar patterns of FC abnormalities, characterized by sensory/subcortical hyperconnectivity, association/subcortical hypoconnectivity, and sensory/association hyperconnectivity. Furthermore, a decrease in FC was observed within each structure category were noted.

Supplementary material: File

Tu et al. supplementary material

Tables S1-S3

Download Tu et al. supplementary material(File)
File 22.8 KB
Submit a response

Comments

No Comments have been published for this article.