Striatum-related spontaneous coactivation patterns predict treatment response on positive symptoms of drug-naive first-episode schizophrenia with risperidone monotherapy

Background: Evidence from functional magnetic resonance imaging (fMRI) studies of schizophrenia suggests that interindividual variation in the stationary striatal functional circuit may be correlated with antipsychotic treatment response. However, little is known about the role of the dynamic striatum-related network in predicting patients’ clinical improvement. The spontaneous coactivation pattern (CAP) technique has recently been found to be important for elucidating the non-stationary nature of functional brain networks.

Methods: Forty-two drug-naive first-episode schizophrenia patients underwent fMRI and T1W imaging before and after 8 weeks of risperidone monotherapy. The striatum was divided into 3 subregions, including the putamen, pallidum, and caudate. Spontaneous CAPs and CAP states were utilized to measure the dynamic characteristics of brain networks. We used DPARSF and Dynamic Brain Connectome software to analyze each subregion-related CAP and CAP state for each group and then compared the between-group differences in the neural network biomarkers. We used Pearson’s correlation analysis to determine the associations between the neuroimaging measurements with between-group differences and the improvement in patients’ psychopathological symptoms.

Results: In the putamen-related CAPs, patients showed significantly increased intensity in the bilateral thalamus, bilateral supplementary motor areas, bilateral medial, and paracingulate gyrus, left paracentral lobule, left medial superior frontal gyrus, and left anterior cingulate gyrus compared with healthy controls. After treatment, thalamic signals in the putamen-related CAP 1 showed a significant increase, while the signals of the medial and paracingulate gyrus in the putamen-related CAP 3 revealed a significant decrease. The increase in thalamic signal intensity in the putamen-related CAP 1 was significantly and positively correlated with the percentage reduction in PANSS_P.

Conclusion: This study is the first to combine striatal CAPs and fMRI to explore treatment response-related biomarkers in the early phase of schizophrenia. Our findings suggest that dynamic changes in CAP states in the putamen-thalamus circuit may be potential biomarkers for predicting patients’ variation in the short-term treatment response of positive symptoms.

Introduction

Schizophrenia is an extremely debilitating psychiatric disorder that has a lifetime prevalence of approximately 0.5%. A considerable body of new evidence from functional magnetic resonance imaging (fMRI) studies supports that striatal functional connectivity (FC) is disrupted in this disease (1–3). Most antipsychotic agents, which compose first-line treatments for schizophrenia, essentially block dopamine D2 receptors in the striatum. It has been suggested that striatal dopamine is central to both schizophrenia pathology and antipsychotic action (1), and interindividual variation in striatal functional circuit is disrupted in schizophrenia and correlated with antipsychotic treatment response (4–6).

Pharmacoimaging studies based on fMRI data have preliminarily revealed antipsychotic agent treatment response-related alterations in FC within striatal subregions (4) as well as between striatal and extrastriatal regions such as the corticostriatal circuitry (7). The premise of these studies is that blood oxygenation level-dependent (BOLD) signals (8), which are obtained without explicit tasks or stimulation during fMRI scans, are temporally associated within functionally related brain areas. Despite a significant amount of existing research, the role of striatal FC as an antipsychotic treatment response biomarker is still being elucidated. The principal imaging approaches available for this goal have generally assumed stationarity that the temporal series do not alter their nature over time. The brain, however, is a complicated system characterized by dynamic FC (9). It has been found that interareal BOLD signal correlations demonstrated large variation over the timescale of tens of seconds to several minutes, confirming the non-stationary characteristic of FC (10), which provides an opportunity to increase the subtle information about the FC network computed during the analysis process.

Previous methods for depicting the non-stationary nature of FC have depended primarily on detecting temporal variability in BOLD signal coherences utilizing sliding window analysis or time–frequency analysis (10–13). These approaches are limited to a few regions of interest (ROIs) and their pairwise relationship. Another shortcoming of these time-domain approaches is the potential confounding effects of non-neurogenic signals such as instrumental drift, head motion, and thermal noise, which were found to be difficult to distinguish from the signals of interest, especially when utilizing shorter sliding windows (10). Recently, the spontaneous coactivation pattern (CAP) technique (14–16) has played a vital role in elucidating the dynamic nature of non-stationary resting state brain networks. As a data-driven approach, the CAP technique, which depends on few mathematical presumptions, is free of the confounding effect from the length of the sliding time window. Therefore, CAP analysis has been increasingly utilized to investigate aberrant brain network dynamics in schizophrenia (15) and other neuropsychiatric disorders such as depression (17), autism (18, 19), and Alzheimer’s disease (20). Yang et al. (15) found that patients with schizophrenia spend more time in CAP states involving the default mode and salience networks, but less time in CAP states involving the fronto-parietal network, and the abnormal dynamic characteristics of functional CAPs are correlated with patients’ symptom severity. However, to date, there is no direct evidence for the role of dynamic striatal FC in predicting antipsychotic treatment response.

This study therefore integrated striatal CAP measures to investigate treatment response- and schizophrenia-related biomarkers in the early phase of schizophrenia. All the patients in this study were treated with risperidone monotherapy for 8 weeks to control for the confounding effects of multidrug treatment on treatment response. Given the previous evidence (1, 4–6), we tested the hypothesis that striatal CAP abnormalities would be pathological traits characterizing schizophrenia and would relate to individual-level variation in antipsychotic response. To test the hypotheses, we applied the robust analytical approach CAPs (14, 15) and fMRI modalities to infer the aberrant state dynamics of the striatum. Given that the striatum consists of the neostriatum (caudate and putamen) and paleostriatum (pallidum) (21), we extracted the 3 regions of interest from the Automated Anatomical Labeling (AAL) template and analyzed their CAPs and dynamic states. The overlapping brain regions, of which activity in striatal CAPs showed significant alterations both at baseline and after treatment, were considered to be involved in both pathological and pharmacological mechanisms, and viewed as effective therapeutic targets. Then, Pearson’s correlation analysis was applied to compute the associations between neuroimaging measurements (including CAPs and dynamic CAP states) of the overlapping regions and the improvement of patients’ psychopathological symptoms to define the treatment response biomarkers.

Materials and methods

Participants

We recruited 42 treatment-naive patients with first-episode schizophrenia (Patients-0W) and 38 education-, gender-, and age- matched healthy controls in the Second Affiliated Hospital of Xinxiang Medical University in China from 2012.12 to 2014.01, China. This fMRI dataset has been used previously by our group (22, 23). Experienced psychiatrists diagnosed patients according to the Structured Clinical Interview for DSM-IV-TR. The illness duration of all the patients in this study was less than 12 months. SCID-non-patient edition was used to scan controls to confirm that they had no history of other psychiatric or neurological disorders. In this observational cohort study, patients were treated with risperidone, and then followed for 8 weeks.

The study was approved by the Ethics Committee of Henan Mental Hospital and Second Xiangya Hospital of Central South University (protocol code S088, 2012). It was conducted in accordance with the Declaration of Helsinki. Written informed consent was acquired from all the participants involved in the study.

Therapy and clinical assessments

The initial dose of risperidone was 2 or 1 mg per day within the first week. After 1 week, a slow titration was then utilized, and the risperidone dose was increased at 1-week intervals until patients showed satisfactory improvement. After 4 weeks, patients who did not exhibit obvious improvement were prescribed up to 6 mg/day of risperidone. All patients were treated with risperidone monotherapy for 8 weeks. Patients were also treated with benzhexol if they had extrapyramidal symptoms. During the 8-week intervention, some patients (n = 19) were given at least one dose (ranging from 1 to 2 mg) of lorazepam; however, this medication could not be taken within 1 week before the second scan. Patients were also not allowed to use antidepressants or mood stabilizers during the intervention period, as these medications are potential synergists. Treatment safety was evaluated weekly through clinical interviews with doctors.

Symptom severity was assessed both before and after 8 weeks of treatment on the day of MRI scanning by using the Positive and Negative Syndrome Scale (PANSS) (24), which primarily evaluates positive (PANSS-P, Items P1–P7) and negative (PANSS-N, Items N1–N7) symptom dimensions.

Imaging data acquisition

The participants underwent T1 weighted (T1W) imaging and fMRI on a 3.0T Siemens MRI scanner (Verio) at the Magnetic Imaging Centre of the Henan Mental Hospital. A standard 16-channel head coil was utilized. We scanned patients both at baseline and at the 8-week follow-up, while healthy controls only underwent the baseline scans. T1W images were obtained sagittally with a gradient echo pulse sequence. The detailed parameters were echo time (TE) = 2.52 ms, repetition time (TR) = 1,900 ms, field of view (FOV) = 250 mm × 250 mm, flip angle (FA) = 9°, slice gap = 0 mm, slice thickness (ST) = 1.0 mm, and number of slices = 176. Resting-state fMRI data were acquired with a gradient-echo echo-planar imaging sequence. The parameters were TE = 30 ms, TR = 2,000 ms, FOV = 240 mm × 240 mm, FA = 90°, time point = 240, matrix size = 64 × 64, ST = 3 mm, voxel size = 3.75 mm × 3.75 mm × 3 mm, and slice gap = 0. Patients were required to relax and keep their eyes closed while staying awake during the scanning.

Four patients withdrew from the follow-up scans. We therefore obtained imaging data from 42 patients at baseline (patients 0W) and 38 patients at the 8-week follow-up (patients-8W) as well as 38 healthy volunteers.

Image data preprocessing

We used the Data Processing Assistant for Resting-State fMRI (DPARSF,¹ DPARSF5.2) (25) to preprocess the fMRI images. The DPARSF is a toolbox for MATLAB (R2019b) and Data Processing and Analysis of Brain Imaging (DPABI)² (26). Slice timing and EPI undistortion correction were performed. Any subject with head motion >1.5° rotation or >1.5 mm translation in any direction was removed. We regressed out the motion parameters, white matter, and CSF from the time series. All subjects’ functional images were normalized to the Montreal Neurological Institute (MNI) space with the standard EPI template (resampled voxel size 3 × 3 × 3 mm³). We also conducted T1 coregistration analysis to preprocess the imaging data. The resampled images were then smoothed with the Gaussian kernel and detrended to remove the linear drift of BOLD signals.

The baseline fMRI data of four patients failed to meet the standards of quality control, and these four patients withdrew from the follow-up imaging scans.

Extraction of ROI time series

We imported the AAL template and extracted the time series of the 3 ROIs (putamen, pallidum, and caudate) in DPABI. The MNI coordinates and the position of the 3 ROIs in the Automated Anatomical Labeling (AAL) template are shown in Supplementary Figure 1. DPABI was used to extract the time series of the 3 ROIs and all other brain regions. We averaged the BOLD signal intensity of all the time points (fMRI time frames) of each ROI.

Extraction and normalization of CAP maps

With the goal of extracting striatum-related spontaneous CAPs, we first chose the top 15% of time points based on previous evidence (14), which demonstrated that the similarity (spatial associations) between the ROI-seeded association map utilizing all time points and the averaged signals of the selected time points rapidly increased while including more time points by decreasing the selection threshold; the spatial associations then reached a plateau after including ∼ 15% of all time points. Then, we selected the top 36 (240 × 15% = 36) time points for the 3 ROIs as the highly activated state. The similarity between our approach and a previous study (14) is the setting of the threshold (top 15% time points), while the difference lies in the selection of the seed-regions.

The three-dimensional brain maps of the 36 highly activated time points from the whole-brain voxel level were extracted and normalized by subtracting the mean BOLD signal intensity at each time point and dividing that value by the standard deviation. We averaged the normalized BOLD signal intensity of the highly activated time points on a regional level. The 3 ROI-related spontaneous CAPs in the control group are shown as an example in Supplementary Figure 2.

We compared the differences in striatum-related CAPs between Patients-0W and controls as well as between Patients-0W and Patients-8 by performing independent samples t-tests and paired t-tests, respectively. GRF correction was performed as well, with voxel level P < 0.01, cluster level P < 0.05, and voxel number >50.

Cluster analysis of dynamic CAP states

To extract further information about the dynamic variations in longitudinal ROI-related CAPs after treatment, after extracting striatum-related spontaneous CAPs, we then classified the 36 detected highly activated time points (in the ROI-related CAPs) according to their spatial similarity utilizing the k-means clustering approach. The above approach, i.e., clustering the top 15% of time points, has been applied previously (14). This process was performed on Dynamic Brain Connectome (Dynamic BC),³ a MATLAB toolbox (Supplementary Figure 3). The CAPs of both Patients-0W and Patients-8W were classified into three states. The brain maps in each state (CAP 1, CAP 2, and CAP 3) were extracted and averaged for Patients-0W and Patients-8W. We used paired sample t-tests to compare the longitudinal changes in BOLD signal intensity in each state, with GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50.

Investigating treatment response biomarkers

Overlapping areas between regions with case-control differences and those with longitudinal differences were viewed as effective treatment targets. We used Pearson’s correlation analysis to compute the correlations between patients’ clinical symptom improvement and longitudinal changes in imaging markers of the overlapping regions, including longitudinal alterations in CAPs and dynamic CAP states.

Results

Demographics and clinical symptoms

We detected no significant differences in the between-group comparisons of age, gender, handedness, and education data of patients and controls (Ps > 0.05; Supplementary Table 1). Patients showed significant clinical improvement in positive symptoms (Ps < 0.001; Supplementary Table 2), while their negative symptoms showed no significant alterations (P > 0.05; Supplementary Table 2).

Abnormal striatum-related spontaneous CAPs in schizophrenia

In the putamen-related CAPs, patients showed significantly increased intensity in the bilateral thalamus, bilateral supplementary motor area (SMA), bilateral medial, and paracingulate gyrus, left paracentral lobule, left medial superior frontal gyrus, and left anterior cingulate gyrus compared with healthy controls (GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50, Figure 1A and Supplementary Table 3). There were no significant between-group differences in the caudate- or pallidum-related CAPs.

FIGURE 1

Figure 1. Between-group comparisons of striatum-related spontaneous coactivation patterns (CAPs). (A) Patients showed significantly increased signal intensity in the putamen-related CAPs at baseline. (B) Patients had significantly increased signal intensity in the pallidum-related CAPs after treatment. (C) Patients had significantly decreased signal intensity in the putamen-related CAPs after treatment.

CAP alterations after treatment and associations with treatment response

In pallidum-related CAPs, patients had significantly increased intensity in the left lingual, fusiform, middle temporal, superior temporal, and angular gyrus after treatment (GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50, Figure 1B and Supplementary Table 4). In the putamen-related CAPs, patients had significantly decreased intensity in the left lingual gyrus (GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50, Figure 1C and Supplementary Table 4). Regarding the caudate-related CAPs, patients showed no significant longitudinal alterations.

There were no overlapping regions between areas with case-control differences and those with 8-week longitudinal differences in patients. We therefore did not further analyze the correlations between patients’ clinical symptom improvement and the longitudinal changes in neuroimaging measurements.

Longitudinal alterations of striatum-related CAP states and associations with treatment response

We detected significant longitudinal differences in the three subregion-related dynamic CAP states (GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50), except the putamen-related CAP 2 (see Table 1).

TABLE 1

Table 1. Longitudinal alterations of the three subregions-related CAP states.

By comparing the differential brain regions of Patients-0W vs. Controls (Figure 2A) and that of Patients-0W vs. Patients-8W (Table 1), we found that thalamus (detected in the putamen-related CAP 1, Table 1) and medial and paracingulate gyrus (detected in the putamen-related CAP 3, Table 1) were the overlapping brain regions. Thalamic signals in the putamen-related CAP 1 showed a significant increase (GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50, Figure 2A), while the signals of the medial and paracingulate gyrus in the putamen-related CAP 3 significantly decreased after treatment (GRF correction, voxel level P < 0.01, cluster level P < 0.05, and voxel number >50, Figure 2B).

FIGURE 2

Figure 2. Longitudinal alterations of striatum-related CAP states and associations with treatment response. (A) Thalamic signal in the Putamen-related CAP 1 showed significant increase after treatment. (B) Medial and paracingulate gyrus signal in the putamen-related CAP 3 showed significant decrease after treatment. (C) The increase of thalamic signal intensity demonstrated significantly positive correlations with the percentage reduction of patients’ positive symptoms (P = 0.0170 with Bonferroni correction). ***Represents P < 0.001, ****represents P < 0.0005.

The increase in thalamic signal intensity in the putamen-related CAP 1 showed significantly positive correlations with the percentage reduction (Patients-0W–Patients–8W/Patients-0W) of PANSS_P (r = 0.5302, P = 0.017 with Bonferroni correction, Figure 2C). We did not detect significant relations between the percentage reduction of PANSS_P and the signal intensity decrease of the medial and paracingulate gyrus in the putamen-related CAP 3 (P > 0.05 with Bonferroni correction).

Discussion

This study first integrates striatal CAP measures and fMRI modalities to explore treatment response biomarkers in the early phase of schizophrenia. As hypothesized, patients in the early phase of schizophrenia had striatal CAP abnormalities, parts of which were correlated with their variation in antipsychotic response of positive symptoms.

The diverse pathways and functional contributions of the striatum, coupled with known dopaminergic function differences in schizophrenia patients (1, 4), imply that the striatum may be a vital contributor to the etiology of this disease. Consistently, our current study identified that patients had increased signal intensity in the putamen-related CAPs, which involved the bilateral thalamus, SMA, medial and paracingulate gyrus, as well as the left paracentral lobule, medial superior frontal, and anterior cingulate gyrus. These findings emphasized the abnormal CAPs in the striatal-thalamic-cortico loop in individuals in the early phase of schizophrenia. This circuit plays an important role in modulating sustained attention, and functional and anatomical abnormalities in this circuitry have been confirmed in extensive MRI studies of schizophrenia (4, 27–30). We utilized a novel striatal CAP neuroimaging biomarker to demonstrate functional abnormalities associated with the pathologic mechanism of schizophrenia in the known striatal-thalamic-cortico circuit. In brief, our current study first identified schizophrenia pathology as resulting from striatal CAPs.

Our current study demonstrated that putamen-related thalamus signals were significantly increased after treatment, and the alterations were associated with patients’ positive symptom improvement. Notably, it is possible that aberrant activity of dopamine in mesolimbic neural circuits may be correlated with psychotic symptoms of schizophrenia (1, 31), and antipsychotic treatment may alter the mesolimbic neural pathways and thus improve patients’ psychotic symptoms. The striatum is an important mesolimbic structure. More importantly, the putamen, part of the dorsal striatum, is the gateway to basal nuclei (32); it receives excitatory afferents from the thalamus and cortical regions and forms the origin of the indirect and direct circuits (33–35). Our findings of the associations between the increase in putamen-related thalamic signal intensity and patients’ positive symptom alterations support the hypothesis that dynamic striatal CAP abnormalities are related to individual-level variation in antipsychotic response.

Approximately 1/3 of schizophrenia patients are defined as non-responders to first-line antipsychotic medications (36). Therefore, there is an urgent need for potential biomarkers to guide individualized precision therapy. In this regard, dynamic changes in CAP states in the putamen-thalamus circuit, which unveiled significant associations with antipsychotic response in the positive symptom dimension, could hold promise as biomarkers for guiding treatment choices. However, these results only apply to acute phase antipsychotic treatment. It remains unclear whether meaningful striatal CAP alterations in the brain can predict patients’ long-term treatment outcomes.

The findings in this study should be interpreted in light of several limitations. First, the sample size is small due to the challenges of recruiting drug-naive first-episode schizophrenia patients who would be treated with antipsychotic monotherapy. Future studies recruiting large samples are warranted. Second, the healthy volunteers were not scanned two times. Thus, the detected longitudinal striatal CAP differences in the patient group should be interpreted with caution, as it cannot be excluded that the differences are induced by the test-retest effect and longitudinal time alterations. Actually, it would be better to include patient control groups taking placebo, which could exclude (to a great extent) other causes of longitudinal alterations in imaging indices, including both test-retest effects, and progressive course of the illness (even over only 8 weeks). However, such a study design is clinically unethical given that known effective antipsychotic agents cannot be withheld from individuals with schizophrenia for 8 weeks. Third, there were no other patient control groups receiving other types of antipsychotic agents. We therefore cannot clarify whether the effects detected in our current study are specific to risperidone. Fourth, participants in this study might have been acclimated to the MRI scanning environment during the second scan, which may also influence our findings. Fifth, while our findings are consistent with a growing body of literature, replication in other exploratory cohorts is needed.

In conclusion, this study is the first to combine striatal CAPs and fMRI to explore schizophrenia- and treatment response-related biomarkers in the early phase of schizophrenia. Our findings suggest that patients in the early phase of schizophrenia had striatal CAP abnormalities, and dynamic changes in CAPs in the putamen-thalamus circuit may be potential biomarkers predicting patients’ variation in the short-term treatment response of positive symptoms.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Henan Mental Hospital and the Second Xiangya Hospital of Central South University. The patients/participants provided their written informed consent to participate in this study.

Author contributions

XZ, MH, HC, JJZ, and SW conceived and designed the study. KW, JBZ, SM, LK, NZ, LL, LXL, and DS performed the investigation and data analysis. XZ, JJZ, and MH contributed to the interpretation of the data. XZ, MH, and KW wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was funded by the National Natural Science Foundation of China (No. 81901357) and the National Key Research and Development Program of China (STI2030-Major Projects, 2021ZD0202000).

Acknowledgments

We would like to thank all participants for their time and excellent cooperation.

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.

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2023.1093030/full#supplementary-material

Footnotes

1. ^ http://www.rfmri.org/dparsf
2. ^ http://rfmri.org/DPABI
3. ^ http://restfmri.net/forum/DynamicBC

References

1. Howes O, Kapur S. The dopamine hypothesis of schizophrenia: version III–the final common pathway. Schizophr Bull. (2009) 35:549–62. doi: 10.1093/schbul/sbp006

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Martino M, Magioncalda P, Yu H, Li X, Wang Q, Meng Y, et al. Abnormal resting-state connectivity in a substantia nigra-related striato-thalamo-cortical network in a large sample of first-episode drug-naïve patients with schizophrenia. Schizophr Bull. (2018) 44:419–31. doi: 10.1093/schbul/sbx067

PubMed Abstract | CrossRef Full Text | Google Scholar

3. McCutcheon R, Abi-Dargham A, Howes O. Schizophrenia, dopamine and the striatum: from biology to symptoms. Trends Neurosci. (2019) 42:205–20. doi: 10.1016/j.tins.2018.12.004

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Li A, Zalesky A, Yue W, Howes O, Yan H, Liu Y, et al. A neuroimaging biomarker for striatal dysfunction in schizophrenia. Nat Med. (2020) 26:558–65. doi: 10.1038/s41591-020-0793-8

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Sarpal D, Argyelan M, Robinson D, Szeszko P, Karlsgodt K, John M, et al. Baseline striatal functional connectivity as a predictor of response to antipsychotic drug treatment. Am J Psychiatry. (2016) 173:69–77. doi: 10.1176/appi.ajp.2015.14121571

PubMed Abstract | CrossRef Full Text | Google Scholar

6. White T, Wigton R, Joyce D, Collier T, Fornito A, Shergill S. Dysfunctional striatal systems in treatment-resistant schizophrenia. Neuropsychopharmacology. (2016) 41:1274–85. doi: 10.1038/npp.2015.277

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Sarpal D, Robinson D, Lencz T, Argyelan M, Ikuta T, Karlsgodt K, et al. Antipsychotic treatment and functional connectivity of the striatum in first-episode schizophrenia. JAMA Psychiatry. (2015) 72:5–13. doi: 10.1001/jamapsychiatry.2014.1734

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Ogawa S, Lee T, Kay A, Tank D. Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proc Natl Acad Sci U.S.A. (1990) 87:9868–72.

Google Scholar

9. Zalesky A, Fornito A, Cocchi L, Gollo L, Breakspear M. Time-resolved resting-state brain networks. Proc Natl Acad Sci U.S.A. (2014) 111:10341–6. doi: 10.1073/pnas.1400181111

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chang C, Glover G. Time-frequency dynamics of resting-state brain connectivity measured with fMRI. Neuroimage. (2010) 50:81–98. doi: 10.1016/j.neuroimage.2009.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Majeed W, Magnuson M, Hasenkamp W, Schwarb H, Schumacher E, Barsalou L, et al. Spatiotemporal dynamics of low frequency BOLD fluctuations in rats and humans. Neuroimage. (2011) 54:1140–50. doi: 10.1016/j.neuroimage.2010.08.030

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Rack-Gomer A, Liu T. Caffeine increases the temporal variability of resting-state BOLD connectivity in the motor cortex. Neuroimage. (2012) 59:2994–3002. doi: 10.1016/j.neuroimage.2011.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Hutchison R, Womelsdorf T, Gati J, Everling S, Menon R. Resting-state networks show dynamic functional connectivity in awake humans and anesthetized macaques. Hum Brain Mapp. (2013) 34:2154–77. doi: 10.1002/hbm.22058

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Liu X, Duyn J. Time-varying functional network information extracted from brief instances of spontaneous brain activity. Proc Natl Acad Sci U.S.A. (2013) 110:4392–7. doi: 10.1073/pnas.1216856110

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Yang H, Zhang H, Di X, Wang S, Meng C, Tian L, et al. Reproducible coactivation patterns of functional brain networks reveal the aberrant dynamic state transition in schizophrenia. Neuroimage. (2021) 237:118193. doi: 10.1016/j.neuroimage.2021.118193

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Xu J, Lyu H, Li T, Xu Z, Fu X, Jia F, et al. Delineating functional segregations of the human middle temporal gyrus with resting-state functional connectivity and coactivation patterns. Hum Brain Mapp. (2019) 40:5159–71. doi: 10.1002/hbm.24763

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Kaiser R, Kang M, Lew Y, Van Der Feen J, Aguirre B, Clegg R, et al. Abnormal frontoinsular-default network dynamics in adolescent depression and rumination: a preliminary resting-state co-activation pattern analysis. Neuropsychopharmacology. (2019) 44:1604–12. doi: 10.1038/s41386-019-0399-3

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Paakki J, Rahko J, Kotila A, Mattila M, Miettunen H, Hurtig T, et al. Co-activation pattern alterations in autism spectrum disorder-A volume-wise hierarchical clustering fMRI study. Brain Behav. (2021) 11:e02174. doi: 10.1002/brb3.2174

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Marshall E, Nomi J, Dirks B, Romero C, Kupis L, Chang C, et al. Coactivation pattern analysis reveals altered salience network dynamics in children with autism spectrum disorder. Netw Neurosci. (2020) 4:1219–34. doi: 10.1162/netn_a_00163

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Ma X, Zhuo Z, Wei L, Ma Z, Li Z, Li H. Altered temporal organization of brief spontaneous brain activities in patients with Alzheimer’s disease. Neuroscience. (2020) 425:1–11. doi: 10.1016/j.neuroscience.2019.11.025

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Raz N, Rodrigue K, Kennedy K, Head D, Gunning-Dixon F, Acker J. Differential aging of the human striatum: longitudinal evidence. AJNR Am J Neuroradiol. (2003) 24:1849–56.

Google Scholar

22. Zong X, Hu M, Pantazatos S, Mann J, Wang G, Liao Y, et al. A dissociation in effects of risperidone monotherapy on functional and anatomical connectivity within the default mode network. Schizophr Bull. (2019) 45:1309–18. doi: 10.1093/schbul/sby175

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Zong X, He C, Huang X, Xiao J, Li L, Li M, et al. Predictive biomarkers for antipsychotic treatment response in early phase of schizophrenia: multi-omic measures linking subcortical covariant network, transcriptomic signatures, and peripheral epigenetics. Front Neurosci. (2022) 16:853186. doi: 10.3389/fnins.2022.853186

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Kay S, Fiszbein A, Opler L. The positive and negative syndrome scale (PANSS) for schizophrenia. Schizophr Bull. (1987) 13:261–76. doi: 10.1093/schbul/13.2.261

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Yan C, Wang X, Zuo X, Zang Y. DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics. (2016) 14:339–51. doi: 10.1007/s12021-016-9299-4

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Huang X, Pu W, Li X, Greenshaw A, Dursun S, Xue Z, et al. Decreased left putamen and thalamus volume correlates with delusions in first-episode schizophrenia patients. Front Psychiatry. (2017) 8:245. doi: 10.3389/fpsyt.2017.00245

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hazlett E, Buchsbaum M, Zhang J, Newmark R, Glanton C, Zelmanova Y, et al. Frontal-striatal-thalamic mediodorsal nucleus dysfunction in schizophrenia-spectrum patients during sensorimotor gating. Neuroimage. (2008) 42:1164–77. doi: 10.1016/j.neuroimage.2008.05.039

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Brickman A, Buchsbaum M, Ivanov Z, Borod J, Foldi N, Hahn E, et al. Internal capsule size in good-outcome and poor-outcome schizophrenia. J Neuropsychiatry Clin Neurosci. (2006) 18:364–76. doi: 10.1176/jnp.2006.18.3.364

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Buchsbaum M, Tang C, Peled S, Gudbjartsson H, Lu D, Hazlett E, et al. MRI white matter diffusion anisotropy and PET metabolic rate in schizophrenia. Neuroreport. (1998) 9:425–30. doi: 10.1097/00001756-199802160-00013

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Murawiec S. [Some questions about the essence of delusions in the light of recent neurobiological findings]. Psychiatr Pol. (2009) 43:403–10.

PubMed Abstract | Google Scholar

32. Kreitzer A, Malenka R. Striatal plasticity and basal ganglia circuit function. Neuron. (2008) 60:543–54. doi: 10.1016/j.neuron.2008.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Alexander G, DeLong M, Strick P. Parallel organization of functionally segregated circuits linking basal ganglia and cortex. Annu Rev Neurosci. (1986) 9:357–81. doi: 10.1146/annurev.ne.09.030186.002041

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Parent A, Hazrati L. Functional anatomy of the basal ganglia. I. The cortico-basal ganglia-thalamo-cortical loop. Brain Res Brain Res Rev. (1995) 20:91–127. doi: 10.1016/0165-0173(94)00007-c

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Haber S. The primate basal ganglia: parallel and integrative networks. J Chem Neuroanat. (2003) 26:317–30. doi: 10.1016/j.jchemneu.2003.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Howes O, McCutcheon R, Agid O, de Bartolomeis A, van Beveren N, Birnbaum M, et al. Treatment-resistant schizophrenia: treatment response and resistance in psychosis (TRRIP) working group consensus guidelines on diagnosis and terminology. Am J Psychiatry. (2017) 174:216–29. doi: 10.1176/appi.ajp.2016.16050503

PubMed Abstract | CrossRef Full Text | Google Scholar