Source: https://www.nature.com/articles/s41467-018-05432-w?error=cookies_not_supported&code=d5905a0f-8e24-4320-9f47-a0f9fc93f205
Timestamp: 2019-04-23 22:02:57+00:00

Document:
Cognitive dysfunction is recognized as a core deficit in many psychotic disorders such as schizophrenia (SZ)1. Studies of psychopathology are now increasingly focusing on understanding how the human brain produces cognition, which may depend on knowledge of its large-scale organization. Mapping cognitive capability onto brain imaging2 has the potential to provide clues by exploiting links among enriched types of imaging and behavioral information for individuals3,4. Here a neuromarker is defined as a brain measure that is associated with a cognitive or behavioral outcome that can predict individual performance5. Neuroimaging techniques like structural and functional magnetic resonance imaging (MRI) have enabled identification of neuromarkers that bring psychiatry from subjective descriptive classification into objective and tangible brain-based measures6,7. For example, Rosenberg et al.8 demonstrated that whole-brain functional connectivity strength may serve as a neuromarker of sustained attention for both healthy and disease assessments. Finn et al.9 showed that functional connectivity (FC) profiles can predict levels of fluid intelligence. Drysdale et al.10 further used FC patterns to define novel subtypes of depression that may benefit treatment outcome prediction. These studies provided steps in developing biomarkers that will allow the field of imaging analysis and psychiatry to move forward to a new era, which motivates our work in identifying neuromarkers that are associated with cognitive composite ability11,12 and specific cognitive domains such as attention, working memory, and verbal learning.
The identified neuromarkers are most meaningful when they are replicable and can be used to predict new, previously unseen subjects. Though receiving increased attention13,14, current neurocognitive investigations are often focusing on one cognitive domain or analyzing one single modality, or performing correlations after separate unimodal analyses15. Therefore, the multi-domain and multimodal cross-information are either missing or not being fully leveraged to improve neuromarker identification, despite the evidence that such information is highly informative16,17,18. To address the above-mentioned issues, we performed successive studies in this paper to answer the following four questions using three independent subject cohorts (n = 294, n = 83, and n = 88), which both have three types of MRI data and are measured with similar cognitive metrics, respectively.
First, what combination of multimodal brain networks will be associated with the global cognitive ability, especially when impaired in schizophrenia? Can the identified multimodal neuromarker signatures be replicated in another independent dataset? Secondly, how are the identified neuromarker signatures associated with other cognitive domains? What are the most correlated cognitive domains? Furthermore, which brain networks are associated with specific representative cognitive domains? What are their commonality and differences across domains? Which imaging modality is more sensitive to cognitive domain discrepancy? Finally, are the identified multimodal brain networks able to predict cognitive performance for new individuals?
To this end, we searched for multimodal neuromarker signatures that can be used to quantify and predict cognitive performance, especially impaired in schizophrenia, by successive multivariate data mining and model generalization. According to the triple network model of major psychopathology proposed by Menon19, the aberrant intrinsic organization and interaction of the salience network (SAN), central executive network (CEN), and default-mode network (DMN) is characteristic of many psychiatric and neurological disorders. Considering both structural and functional dysfunction in these three networks have been found in schizophrenia20,21,22 and linked with cognitive deficits23,24, we hypothesize that the modality-specific SAN, CEN, and DMN would have pivotal roles in cognitive impairment in schizophrenia23, which may consist multimodal signatures that are defined as modality-specific brain regions25 exhibiting similar subject-wise covariation that can jointly predict cognitive performance for unseen individuals.
In order to address the first three questions, cognitive composite and multiple domain scores were used as references, respectively, to guide the three-way multimodal MRI data fusion in a discovery cohort: FBIRN (Function Biomedical Informatics Research Network, n = 294) by our proposed model26 (see Method for the analysis flowchart). Results were further replicated in an independent cohort (n = 83). After discovering the multimodal signatures particularly associated with each of the four cognitive domains (including Computerized Multiphasic Interactive Neuro-cognitive System [CMINDS]27 composite, CMINDS attention, CMINDS working memory, and CMINDS verbal learning), we compared them in two ways: firstly, across domains—to reveal the commonality and uniqueness of the domain-specific multimodal signatures; secondly, across modalities—to evaluate which imaging modality is more sensitive to cognitive domain discrepancy. Finally, by taking advantage of the above-extracted brain networks in one cohort, we built predictive models of cognitive metrics by performing linear regression on these identified neuromarkers. Remarkably, these models were successful in predicting the corresponding cognitive metrics for new individuals in another two independent cohorts: UNM (University of New Mexico, n = 83) and COBRE (Center for Biomedical Research Excellence, n = 88)28 This validates the generalizability of the identified modality-specific potential neuromarkers, which might be broadly applicable as predictors of multi-domain cognitive scores for new individuals.
We aim to identify multimodal co-varying and modality-specific brain networks associated with composite cognitive scores. Methods section presents the whole-analysis flowchart of data fusion and prediction. First, both the CMINDS composite score and the MCCB composite score were used as the reference for FBIRN and UNM cohorts, respectively, with MCCAR + jICA26,29 (multimodal canonical correlation analysis with reference plus joint independent component analysis). Then, to test the replicability of the identified multimodal networks, we calculated the permuted spatial correlations of their maps between cohorts and summarized the most affected cognitive domains related to schizophrenic deficit. Furthermore, cognitive domain scores of CMINDS attention, CMINDS working memory, and CMINDS verbal learning were used as reference to guide the three-way MRI fusion, respectively, aiming to identify the domain-common and modality-specific neuromarkers. Finally, after extracting the neuromarker features across multiple cognitive domains, we built linear regression models to predict individualized cognitive scores in FBIRN cohort, which were further generalized to predict corresponding cognitive measures of unseen subjects in two independent cohorts (UNM and COBRE).
The FBIRN Phase III study consisting of 147 SZs (mean 39.5 ± std11.7) and 147 HCs (37.4 ± 11), who were each measured by a neuropsychological battery that includes six neurocognitive domain tests called CMINDS27.
Three representative MRI features (fractional amplitude of low-frequency fluctuations (fALFF) from resting-state functional MRI (rs-fMRI), gray matter (GM) density from structural MRI (sMRI), and fractional anisotropy (FA) from diffusion MRI (dMRI)) were combined by a fusion with reference model26, in which CMINDS composite cognitive scores were adopted as the reference. MCCAR + jICA26 can simultaneously maximize the inter-modality covariation and correlations of certain imaging components with the referred clinical or cognitive measures (see more details in the Methods section). Based on the minimum description length (MDL)30 criterion, 24 components were estimated. After joint decomposition, independent components (ICs) and their subject-wise loadings were derived for each modality.
The supervised fusion model produced one joint IC showing significant correlations with the CMINDS cognitive composite scores for all modalities, and we denote this IC as ICref (Fig. 1). Figure 1a displays the spatial maps of FBIRN_ICref_composite for each modality. Figure 1b indicates that FBIRN_ICref_composite has positive correlations with the CMINDS composite scores for all modalities (r = 0.486*, 0.262*, 0.430* for sMRI, dMRI, and fMRI, respectively, * means FDR corrected for multiple comparisons, which represents the same meaning for all correlations and p values in this paper), namely, the higher the loadings, the better the cognitive performance. The p values derived from a permutation test for the correlation (details provided in Supplementary Note 1 and Supplementary Fig. 1) between FBIRN_ICref_composite and composite cognitive scores are pperm = 1.3×10−4*, 0.002, 1.0×10−4* for sMRI, dMRI, and fMRI, respectively (Fig. 1b). Moreover, the FA and GM components also show significant correlations with Positive and Negative Syndrome Scale (PANSS)31 negative scores (dMRI: r = −0.162, p = 0.05, sMRI: r = −0.285*, p = 5.3×10−4), which share many features with cognitive impairment32,33. No significant correlations were found with PANSS positive scores. The identified brain regions in FBIRN_ICref_composite are summarized in Supplementary Table 1 for fALFF components (Talairach labels), FA (WM tracts, from John Hopkins Atlas), and GM (Montreal Neurological Institute labels) respectively. Two sample t-tests were also performed on IC loadings between patients and controls. Notably, the FBIRN_ICref_composite also differs significantly in each feature with lower means in SZ (p = 7.5×10−17*, 2.7×10−7*, 5.3×10−12*), as shown in Fig. 1c, indicating a co-varying fALFF-FA-GM decrease in SZ. We denote this FBIRN_ICref_composite as a joint multimodal signature that is closely associated with the cognitive composite scores.
Due to the interferential effects of varying demographic distributions and measurement conditions, few neurocognitive studies have been strictly replicated for independent cohorts. Here, to validate the replicability of the above-extracted multimodal networks related to composite cognitive scores, we analyzed an independent cohort including 44 HCs and 39 SZs collected from UNM, which incorporated a different but related cognitive battery called the MATRICS (Measurement and Treatment Research to Improve Cognition in Schizophrenia) Consensus Cognitive Battery (MCCB)34. There was no subject overlap between the UNM and FBIRN and COBRE data. MCCB composite was used as reference to guide the 3-way fusion of UNM data. Twenty components were estimated based on MDL criterion30. A similar joint ICref (denoted as UNM_ ICref_composite) was identified to be associated with the MCCB composite scores (r = 0.271*, 0.281*, 0.311*) and also group-discriminating (p = 0.005*, 0.043, 0.002* for sMRI, dMRI, and fMRI, respectively), as displayed in Fig. 2. The identified areas in UNM_ ICref_composite are summarized in Supplementary Table 2. The p values from the permutation test for the correlation between UNM_ICref_composite and composite cognitive scores are 0.02, 0.01, 0.001 for sMRI, dMRI, and fMRI, respectively (Fig. 2b).
Notably, the identified multimodal networks exhibit substantial spatial overlap between discovery cohort (FBIRN) and replication cohort (UNM), while the cognitive performance was evaluated from two different systems (CMINDS and MCCB), which are similar but not identical in domain tasks27. Figure 3 summarizes the similarity of global cognition-associated brain maps between independent cohorts and modalities. The cross-cohort correlation between the identified components are GM: r = 0.51, pperm = 8.0×10−4, FA: r = 0.59, pperm = 2.0×10α; fALFF: r = 0.45, pperm = 0.002, where the significance pperm were resulted from 10000 permutations. Details provided in Supplementary Note 2 and Supplementary Table 3.
As shown in Fig. 3, sMRI and dMRI features show high spatial similarity (r > 0.5) between cohorts. Particularly, decreased GM volume in the salience network (SAN, including dorsal anterior cingulate cortex [ACC] and insula), dorsolateral prefrontal cortex (dlPFC, key node of CEN), and subcortical clusters (including striatum and thalamus) were detected in SZ for both cohorts. For dMRI, the FA reduction in SZ on white matter tracts (such as corpus callosum35, superior longitudinal fasciculus [SLF]36 and anterior thalamic radiation [ATR]37) were common for both cohorts. For fMRI, the prefrontal regions and posterior DMN were more consistent between cohorts with a lower fALFF in SZ.
Therefore, three key brain networks (SAN, CEN, pDMN) were identified in different modalities to be closely related to cognitive performance, which are also replicated in two independent cohorts (Figs. 1–3 and Supplementary Tables 1, 2). Note that these three networks have previously been associated with cognitive deficits in SZ in multiple reports38,39,40,41, but never in a multimodal analysis across multiple cohorts and domains. More importantly, the three networks are in accord with our initial hypothesis based on the triple network theory19 that aberrant organization and functioning within CEN, SAN, and DMN is closely associated with cognitive impairment in psychiatric disorders. Next, we investigated changes in the three networks in SZ in detail, as well as the associations of the identified multimodal components with other cognitive domains and different modalities.
Cognitive function encompasses a variety of cognitive domains, including but not limited to attention, working memory, visual/verbal learning, decision-making, reasoning, and social cognition12. We next investigate the association between the identified multimodal components with different cognitive domains by calculating loadings of FBIRN_ICref_composite in each modality with six CMINDS domain scores (Supplementary Table 4). Results indicate that the top three most correlated domains are attention (r = 0.441*, 0.269*, 0.333*), working memory (r = 0.402*, 0.240*, 0.341*), and verbal learning (r = 0.427*, 0.232*, 0.361*) for sMRI, dMRI, and fMRI, respectively, as displayed at the bottom of Fig. 4.
Figure 4 summarizes the key neurocognitive networks identified in FBIRN_ICref_composite and their associated cognitive domains with different colors. Specifically, there are higher fALFF values in pDMN (blue), including the ventromedial prefrontal cortex (vmPFC) and posterior cingulate cortex (PCC), but lower fALFF in the middle temporal gyrus (MTG) and primary visual cortex (V1) associated with a higher cognitive composite score42 (see also Fig. 1). Additional key compartments of CEN43 (purple), including dlPFC, superior temporal gyrus (STG), and posterior parietal cortex (PPC), were also identified in fALFF. For GM, higher GM volume in the salience network44 (cyan; ACC and bilateral insula) and subcortical regions (red; including hippocampus, caudate, and thalamus) are correlated with higher cognition, especially in the domain of attention. For dMRI, SZ shows lower FA values than HC in major white matter tracts (orange arrows) including forceps major (FMAJ), forceps minor (FMIN)35, SLF36, and ATR37, which can potentially connect the identified brain regions in GM or fALFF. For instance, the FMIN (genu of the corpus callosum) interconnect the left and right frontal lobe, linking dlPFC and vmPFC (both were detected in fALFF and GM), and provide evidence that disrupted anatomical connections in SZ in the anterior corpus callosum may relate with the PFC impairment in both function and structure, as well as cognitive deficits in multiple domains such as working memory and verbal learning45. Similarly, the identified regions in CEN in fALFF are also connected by the SLF, which has been implicated in poor executive function in SZ35. These results clearly suggest that alterations in one modality can be associated with correlated changes in distant, but connected regions in another modality (Figs. 3 and 4) and multimodal fusion proves to be a powerful tool to reveal this association.
Following the above findings, attention, working memory, and verbal learning are the top three domains that are most correlated with the multimodal networks associated with global cognition and all are impaired in SZ. Next, we investigated the multimodal brain networks that are specifically associated with these three domains, respectively; which were then further compared in two ways: firstly, across domains to reveal the commonality and uniqueness of the multimodal brain regions associated with particular cognitive domains and secondly, across modalities to evaluate which imaging modality is more sensitive to cognitive domain discrepancy, as illustrated in Fig. 5.
Toward this goal, CMINDS domain scores of working memory, attention or verbal learning were used as reference, respectively, to perform the same reference-directed fusion analysis as we did using CMINDS composite scores (see Methods section). For each cognitive domain, we extracted one joint ICref that is significantly correlated with specific domain scores and also group-discriminative in all modalities (all passed FDR correction at p < 0.05; see Supplementary Table 5). In order to quantify the robustness of the cognition-brain correlation across groups, we performed partial correlation to minimize the group effect. As shown in Supplementary Table 6, in the case of partial correlation, the cognition-imaging correlations remain significant (FDR corrected) after controlling for diagnosis in all four domains.
Figure 5a–d displays spatial maps of the identified FBIRN_ICref_composite (red), FBIRN_ICref_memory (green), FBIRN_ICref_attention (cyan), and FBIRN_ICref_learning (magenta), respectively. Figure 5e indicates the pairwise cross-domain similarity calculated by performing 3D correlations among the spatial maps of FBIRN_ICref_composite, FBIRN_ICref_memory, FBIRN_ICref_attention, and FBIRN_ICref_learning for GM, FA, and fALFF, respectively. Note that the darker blue and larger shading denotes a higher correlation. GM and FA maps demonstrate more consistent patterns across cognitive domains (r > 0.6), while fALFF maps exhibit more variance (r < 0.5 but still significant) and are more sensitive to domain differences. In particular, our results revealed a pattern of overlapping GM reduction in SAN and dlPFC across all domains, which is consistent with most recent meta-analyses in SZ patients. Hippocampus volume reduction is associated with the composite, working memory, and verbal learning domains. SZ patients also show decreased FA values in FMAJ, FMIN and ATR for all domains. For fALFF, the prefrontal cortex is consistently observed in all domains, while subcortical regions are more variable across domains.
In order to examine the domain-specific fALFF patterns more precisely, we display the representative slices of FBIRN_ICref_composite, FBIRN_ICref_memory, FBIRN_ICref_attention, and FBIRN_ ICref_learning. Figure 6a shows both positive and negative fALFF patterns. Figure 6b plots correlations between ICref loadings with its corresponding domain scores after controlling for diagnosis (Supplementary Table 6). It is clear that vmPFC and dlPFC are the domain-common fALFF regions that have been widely implicated in higher-order cognitive functions (e.g., attention, working memory, and verbal learning), as reflected in our results. Lower fALFF values in SZ were observed, corresponding to poorer cognitive performance than HCs46,47, which is also consistent with two related FBIRN fMRI studies48. In addition, pDMN was identified only for the cognitive composite and working memory domains, while other subcortical and cortical regions such as thalamus, hippocampus, STG, and visual cortex occur differently depending on the different domains. For instance, the increased fALFF in hippocampus in SZ is only detected for the verbal learning and composite scores, which is consistent with a recent result that suggests low-frequency hippocampal–cortical activity drives brain-wide resting-state fMRI connectivity and contributes to cognition49.
Our findings suggest that functional measures such as fALFF is more sensitive to differences between cognitive domains, whereas structural brain measures (GM/FA) are more similar and replicable across multi-domain cognitive impairment. Furthermore, in order to characterize the common core brain networks across the four cognitive domains (composite, working memory, attention, verbal learning), we extracted the overlapped brain regions of FBIRN_ICref_composite, FBIRN_ICref_memory, FBIRN_ICref_attention, and FBIRN_ICref_learning (see details in the Methods section on neuromarker extraction) and summarized the domain-common brain patterns for each modality as shown in Fig. 7a. Evidently, the SAN in GM (red), corpus callosum (CC) in FA (green), and PFC (dlPFC + vmPFC) in fALFF (yellow) are shared among all four domains, which may be treated as potential neuromarkers in each modality to characterize cognition quantitatively.
Based on the masks displayed in Fig. 7a, we further extracted the corresponding regions of interest (ROIs) in each modality to examine their original group differences and association with cognition (before fusion). The mean voxel values of each ROI were calculated for each subject and each modality. As expected, two sample t-tests between HC and SZ show significant group differences for GM_SAN (p = 6.8×10−9*), FA_CC (p = 0.001), and fALFF_PFC (p = 7.0×10−9*), respectively. The effect size for separating groups was also computed with Cohen’s d = 0.698, 0.487, and 0.623 for GM_SAN, FA_CC, and fALFF_PFC, respectively. In addition, the GM_SAN, FA_CC, and fALFF_PFC also show significant correlations with all CMINDS cognitive domain scores, as listed in Supplementary Table 7. All above results suggest that these extracted ROIs may have pivotal roles as neurocognitive substrates, i.e., so-called neuromarker signatures, as they are common core brain networks among multiple cognitive domains, and are widely impaired in mental disorders50.
Based on Eq. (1), correlation of r = 0.463* was achieved between the estimated CMINDS composite scores and its true values (Fig. 7b).
To test the generalizability of the identified four neuromarkers and the prediction model, we then extracted the same four ROIs in UNM and COBRE cohorts through masks (obtained from FBIRN cohort) and applied them to Eq. (1), to predict the unseen MCCB composite scores for both UNM and COBRE data (preprocessed using the same pipeline as in FBIRN). The model generalization from CMINDS to MCCB work well not only for composite score, but also for most cognitive domains (see Table 1 below), implicating the appreciable generalizability of the cognition-predictive models. CMINDS and MCCB are two similar but not identical cognitive measurement systems27; therefore, the cross-cohort generalization is a powerful evidence to validate the predictability of global cognition using the identified neuromarkers and model. As shown in Fig. 7d and c, Pearson correlations of r = 0.406 and r = 0.236 were achieved between the estimated MCCB composite scores and its true values for COBRE (42 HCs/46SZs) and UNM (41HCs/37SZs), respectively, suggesting good generalizability of the proposed cognition-prediction model. Note that the prediction models in Fig. 7b–d are the same, i.e., training in the FBIRN cohort to predict the CMINDS composite. Moreover, we also performed the prediction analysis within each group (HC or SZ) based on the four neuromarker signatures, i.e., using the group model trained by FBIRN to predict UNM + COBRE. The generalization in either case works well as shown in Supplementary Fig. 2 and Supplementary Note 3.
The same models are then used to predict the corresponding domain scores in MCCB with UNM and COBRE cohorts. Note that the beta weights β0, β1, β2, β3, and β4 are different for each domain. Table 1 lists the correlations between the estimated domains scores and true values for all three cohorts. It is clear that all CMINDS domain scores can be predicted by these four neuromarkers. More importantly, such a prediction can be generalized to most of the MCCB domains in two independent UNM and COBRE cohorts, suggesting that the identified neuromarker signatures provide a set of broadly applicable predictors on cognitive performance of either global or specific cognitive domains.
Furthermore, to examine the specific contribution power of each of the four neuromarkers on predicting multiple cognitive domains, we plotted the beta weights (β1, β2, β3, β4 as in Eq. (2)) for GM_SAN, FA_CC, fALFF_ PFC, and fALFF_pDMN, respectively, as in Fig. 8, which correspond to linear regression results in Table 1.
Evidently, GM_SAN and FA_CC have the higher beta weights for the majority of cognitive domains (p = 0.002, ANOVA test on |beta weights| among four neuromarkers). While structural neuromarkers were more predictive and less variable, functional MRI measures were more sensitive to domain differences. This can be reflected by the standard deviations of beta weights across cognitive domains, (i.e., GM_SAN:0.088;FA_CC:0.053;fALFF_PFC:0.130;fALFF_pDMN:0.134), indicating that the contributing power of fALFF shows more fluctuation than GM and FA across different cognitive domains.
In this study, we searched for multimodal neuromarker signatures that can be used to quantify and predict cognitive performance by successive data mining and model generalization. To the best of our knowledge, this is the first attempt to utilize cognition as a reference to guide the three-way multimodal MRI fusion, and to replicate findings in two independent cohorts and across multiple cognitive domains. Our goal was to answer the four challenging issues regarding the cognitive imaging biomarkers. By successfully predicting different measures of cognition in three independent datasets, GM_SAN, FA_CC, fALFF_PFC, and fALFF_pDMN demonstrated great potential as multimodal neuromarker signatures of generalized cognition.
One major finding was that CEN, SAN, and pDMN are three key networks that are especially important for uncovering multiple cognitive functions and their dysfunction in SZ20,21,22. More importantly, they are consistent with the triple network model of major psychopathology19. Particularly, in our study, bilateral GM volume reduction has been seen in SAN in schizophrenia patients and were closely associated with all cognitive domain scores. SAN deficits in schizophrenia have been linked to reality distortion, suggest that SAN abnormality may contribute to impaired salience that is associated with hallucinations and delusions in schizophrenia51. Together with other consistent findings51,52, we speculate that SAN in GM may have a crucial role as a structural substrate for neurocognition19,53. Moreover, CEN was detected in both fALFF and GM in our results, which is important for maintaining and manipulating information in decision-making54,55, working memory, and problem solving. In addition, pDMN was identified in fALFF maps related with cognitive composite and working memory domains. These modality-specific triple networks along with FA_CC covary subject-wisely, and were replicated in an independent cohort to be associated with global cognitive function (Fig. 3). Although they have previously been found to be associated with cognitive deficits in multiple reports38,39,40,41, but never in a multimodal joint analysis across multiple cohorts. Whereas our work provides evidence that a synthesis and interaction exist among the three intrinsically coupled networks19 in both brain function and structure, which are systematically engaged during cognition and impact multiple cognitive domains.
Another interesting finding was that the identified multimodal networks are most significantly correlated with three cognitive domains: attention, working memory, and verbal learning, which belong to key higher cognitive functions and are often severely impaired in SZ20,21,56. As seen in our results, working memory is correlated with fALFF and GM, which is in line with CEN (i.e., dlPFC, STG, and PPC) detected in fMRI, and partly in sMRI (dlPFC), where SZ patients showed decreased fALFF and GM values with lower domain scores. This is consistent with the fact that memory shows pronounced deficits in SZ, with working memory primarily affected22,47. In addition, attention encompasses a variety of functions that includes information selection, enhancement of selected information and inhibition of unselected information. We detected GM volume reduction in key nodes of the SAN (bilateral insula and ACC in blue frame as in Fig. 4) and dlPFC associated with attention deficits in SZ. Accordingly, SAN is known to be involved in a variety of monitoring, attention switching, and decision-making processes50. Furthermore, regions including AG and hippocampus (in green frame) were commonly detected in fMRI and sMRI, accompanying the domains of verbal learning that correlated with fALFF and GM. Specifically, the Broca’s area (BA 4457) and posterior Wernicke’s area (BA 2258) were, respectively, identified in fMRI and sMRI images and both are included in the language learning network. Therefore, the identified multimodal signatures represented a cortical–subcortical circuit in fALFF and GM, which may account for higher cognitive deficits in working memory, attention, verbal learning, and composite cognition in SZ59. Furthermore, FA map occurred in ATR (connecting the frontal–subcortical circuits) and SLF (linking frontal–parietal–temporal circuits), implying that alterations in one modality can be linked with changes in distant but connected regions in another modality. This strength in multimodal fusion reveals associations that cannot be discovered by separate multimodal analyses.
Third, our results indicate that structural features are more spatially consistent and replicable across multiple cognitive domains, while functional maps may better differentiate cognitive domains. Specifically, our results revealed one pattern of overlapped GM reduction in SAN extending across four domains (composite cognition, working memory, attention and learning) and in two cohorts. Similar findings were reported, indicating that GM lesions converged in ACC and dorsal-insula (key nodes of SAN)60 based on the voxel-based morphometry meta-analysis across six diagnostic groups (SZ, bipolar disorder, depression, addiction, obsessive-compulsive disorder, and anxiety), and affected most categories of psychiatric illness61. Coincidentally, reduced FA in corpus callosum that links the bilateral hemispheres has been reported in a recently ENIGMA large-scale coordinated study of white matter microstructural differences in schizophrenia62, while here we replicated it in a multimodal and multi-domain co-varying manner for the first time. By contrast, fALFF exhibit more spatial variance in cortical regions across a wide range of domains, except with vmPFC and dlPFC being commonly detected that closely associate with high-order cognitive functions, (i.e., working memory38, verbal learning, and attention63); whereas, subcortical regions such as hippocampus, thalamus are more task-specific. We provide new insights into fundamental cognitive imaging biomarkers by identifying co-varying multimodal neuromarker signatures in which functional features are more sensitive to cognitive domain differences, while brain structural and anatomical property show more unifying impaired patterns associated with multi-domain cognitive decline64.
Finally, the primary goal of neurocognitive imaging studies is to identify neuromarkers that can predict individual educational or health outcomes8,65. After discovering a set of stable, domain-shared, multimodal neuromarker signatures (Fig. 7a), we defined a linear regression model (Eq. (1)) that is able to predict individual cognitive composite scores, which can be generalized to two independent cohorts (UNM and COBRE) to predict new unseen subjects on similar but different cognitive measures (Figs. 7 and 8). We have validated the feasibility of the four neuromarkers as predictors of general cognitive performance, suggesting that they are descriptive in nature to cognitive function and may serve as potential biomarkers for cognitive impairments in SZ. Note we obtained high predictive power in this study by selecting four neuromarkers representing common patterns among four domains (Fig. 5). Future work could focus on incorporating domain unique patterns, which we expect to further increase predictive power by better leveraging the sensitivity of the fALFF feature, which showed more domain unique patterns.
A possible limitation of this work is that reference-directed fusion works on extracted features, rather than the original imaging data (e.g., using fALFF instead of 4D fMRI data). Although some of the temporal information was lost using this method, a “feature” tends to be more tractable than working with the large-dimensional original data66 and provides a simpler space in which to link the data67. In future work, we plan to incorporate features such as functional network connectivity matrices68, dynamic states69, and structural morphometric measures as fusion input to capture both temporal and spatial co-alterations. Furthermore, most participants were receiving antipsychotic medication at the time of scanning (medication information can be found in Supplementary Table 8 and Supplementary Note 4). In our current study, the correlation between medication dose and cognitive domain was not significant (Supplementary Table 9). In addition, no imaging voxels showed a significant correlation with medication dosage for any of the three modalities, demonstrating that medication has little or at best a very subtle effect on brain imaging in our current data. These results support our claim that the identified replicable multimodal co-varying patterns are associated with cognition but not medication exposure.
In summary, by jointly analyzing three MRI data in a supervised, cognition-guided multimodal fusion model, we successfully identified four multimodal neuromarker signatures that can be replicated in two independent cohorts. A triple brain network including GM_SAN, fALFF_CEN, and fALFF_pDMN may have a crucial role as anatomical substrates of neurocognition, which achieved preferable accuracy for predicting the multi-domain cognitive scores and can be generalized to predict previously unseen individuals. Our results suggest functional features are more sensitive to differentiate cognitive domains, while the salience network in GM and corpus callosum in FA are highly consistent and predictive on multiple cognitive domains. This work defines modality-specific brain networks that may be broadly applicable as neuromarkers of cognitive impairment, and our assessment can help better understand how functionally and anatomically connected brain systems both engender and constrain cognitive functions in SZ.
Based on supervised learning, we previously have developed a reference-guided fusion model called MCCAR + jICA26. Assume Xk represents multimodal dataset and each is a linear mixture of components Ck with a nonsingular mixing matrix Ak, k = 1, 2, 3, denoting the modality. Namely, Xk = AkCk, where Xk is a subjects-by-voxels feature matrix and Ak is a subjects by number of components (M) mixing matrix. MCCA with reference (MCCAR) imposes an additional constraint to maximize not only the covariations among loadings of each modality, but also the column-wise correlations between Ak and the reference signal, as shown in Eq. (3).
where ref is an N × 1 vector, denoting the referred measure, N is the subject number. corr(Ak,Aj) is the column-wise correlation between Ak and Aj, and corr (Ak,ref) is the column-wise correlation between Ak and ref. After optimization by MCCAR, we can obtain the potential target components Ci that are correlated with ref in each modality, as well as being most correlated across subjects between modalities. Then joint ICA is further applied to on the concatenated maps of [C1,…,CM], in order to keep the modality linkage of the potential target components and maximize the spatial independence. The final independent components (ICs) Sk, along with their mixing matrices Dk are obtained by linear source decomposition.
Among Sk, one or more joint independent components (ICs) that are specifically correlated with the ref will be identified as the target ICref. Therefore, by incorporating prior information, MCCAR + jICA enables identification of a joint multimodal component(s) that has robust correlations within referred measures and amongst themselves (inter-modality correlations), which may not be detected by a blind N-way multimodal fusion approach. For more details, please refer to our method paper on MCCAR + jICA26 and λ tuning for Supplementary Note 5.
Here the cognition-guided fusion analyses were implemented by feeding the preprocessed MRI features (i.e., GM density from sMRI, FA from dMRI and fALFF from resting-state fMRI) into MCCAR + jICA, as displayed in Fig. 9a and b. We sought to investigate the target joint independent components (ICref) that are not only significantly correlated with referred cognitive scores, such as CMINDS composite, MCCB composite or CMINDS domain scores, but also indicated significant linked functional–anatomical–structural alterations between schizophrenia and controls, i.e., group-discriminative.
After identifying multimodal networks of four CMINDS cognitive domains from the FBIRN data (Fig. 9c, d), we extracted the brain regions which were consistently involved within each modality as potential neuromarkers. Take the feature extraction of GM as an example. After converting component GM_ICref into Z scores and thresholding at |Z| > = 2, masks of GM_ICref for each of the four cognitive domains (composite, attention, working memory and verbal learning) were generated. A map of regions included in each of these four GM masks reveals the common ROIs across the four cognitive domains. The final GM mask of GM_SAN is shown in Fig. 7a. This mask of GM was then used to extract ROI features from every subject. The mean of the voxels within the obtained ROI was calculated for each subject, generating a N subj × 1 feature vector for GM_SAN. The other modalities (dMRI and fMRI) were processed in the same way to get the FA_CC and fALFF_PFC feature vectors. The fourth fALFF_pDMN feature was extracted from fALFF_ICref_composite from the FBIRN data. The resulting regions were those included in the triple network hypothesis19. Finally, we formed a feature matrix in dimension of N subj × 4 for the FBIRN data. For the UNM and COBRE cohorts, the ROI features of each modality were extracted by applying the four masks generated from the FBIRN to the UNM and COBRE data. Following this, the mean of each ROI was calculated for each subject, resulting in a N subj × 4 feature matrix for the UNM and COBRE cohorts, respectively.
After neuromarker extraction, for the FBIRN data, each of the four neuromarker vector was normalized to mean = 0, std = 1. These vectors were then treated as the linear regressors and the corresponding cognitive scores were treated as the targeted measures; together, they were input into the multiple linear regression model to obtain a linear equation for an estimate of the target measures. The same regression model (beta weights) achieved from the FBIRN prediction was then applied to the regressors obtained from UNM and COBRE data. Pearson correlations between true and predicted cognitive scores of each domain were assessed, as seen in Figs. 7, 8 and 9e. Based on the CMINDS regression model obtained from FBIRN cohort, if the estimated MCCB scores of UNM and COBRE cohorts are significantly correlated with the true MCCB scores, the generalization of predictability of cognition based on the four selected neuromarkers can be recognized to some degree. The beta weights of linear regression models for each predicted cognitive domain are displayed in Fig. 8.
Two independent data cohorts were used in this study. One is recruited from FBIRN phase III datasets including 147 SZs (39.5 × 11.7) and 147 HCs (37.4 × 11) that were matched for gender, age, handedness and race distributions. The demographic and cognitive information are summarized in Supplementary Table 10. All subjects were collected from seven FBIRN consortium sites (University of California Irvine, University of California Los Angeles, University of California San Francisco, Duke University, University of North Carolina, University of New Mexico, University of Iowa, and University of Minnesota). Each dataset included diagnosis, age at time of scan, gender, illness duration, symptom scores, and current medications when available, which were shared by each research group according to their site’s protocols. Inclusion criteria required all participants to be adults between the ages of 18 and 65 years. Diagnosis of SZ was confirmed by trained raters using the Structured Clinical Interview for DSM-IV (SCID)70. All patients were on a stable dose of antipsychotic medication either typical, atypical, or a combination for at least 2 months. Detailed medication information could be found in Supplementary Table 8. Symptom severity was rated using PANSS. All SZs were clinically stable at the time of scanning. In addition, HC participants were excluded for past or current psychiatric illness based on SCID assessment or for having a first-degree relative with a diagnosis of an Axis-I psychotic disorder. Written informed consent was obtained from all participants under protocols approved by the Institutional Review Boards at each study site. Demographic information for subjects of each site is provided in Supplementary Table 11.
The UNM cohort consisted of 39 SZs (35.6 ± 13.1) and 44 HCs (36.3 ± 12.5) who were collected from the University of New Mexico and recorded with MCCB scores. The third cohort include 42 patients with schizophrenia (39.3 ± 13.2) and 42 HCs (40 ± 11) who were collected from the COBRE project28. Details of the cognitive and clinical information of both UNM and COBRE data can be found in Supplementary Table 12 and Supplementary Table 13.
The resting-state fMRI data for the FBIRN cohort was collected on six 3T Siemens and one 3T General Electric (GE) scanner. The imaging protocol for the resting-state scans at all sites was the same: a T2∗-weighted AC-PC aligned echo planar imaging sequence (TR = 2 s, TE = 30 ms, flip angle = 77°, 3.4 × 3.4 × 4 mm with 1 mm gap, 162 frames, 5:38 min). For the resting scan, subjects were instructed to lie still with eyes closed.
DMRI data for the FBIRN cohort were also acquired on six 3T Siemens and one 3T GE. All parameters for these two scanners were the same except for TE (Siemens: 84 ms, GE: 81.7 ms). The rest of the parameters for both Siemens and GE were as follows: TR = 9 s; field of view (FOV) = 256 × 256 mm; slice thickness = 2 mm; number of slices = 72; slice gap = 2 mm; voxel resolution 2 × 2 × 2 mm; flip angle = 90°; number of diffusion gradient directions = 30, b = 800 seconds/mm2, and 5 measurements with b = 0. All images were registered to the first b = 0 image by FMRIB Linear Image Registration Tool (FLIRT: http://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FLIRT).
High-resolution T1-weighted brain imaging data for the FBIRN cohort were also acquired on six 3 T Siemens and one 3T GE. Siemens scan parameters were TR = 2.3 s, TE = 2.94 ms, flip angle = 9°, resolution = 256 × 256 × 160. GE scan parameters were TR = 5.95 s, TE = 1.99 ms, flip angle = 12°, resolution = 256 × 256 × 166. All scans covered the entire brain with FOV = 220 mm2, voxel size = 0.86 × 186 × 1.2 mm3.
The fMRI data were preprocessed using the automated analysis pipeline in SPM8 (http://www.fil.ion.ucl.ac.uk/spm) as following: motion correction, slice timing and normalization to MNI space, including reslicing to 3 × 3 × 3 mm voxels. We removed subjects who had framewise displacements (FD) exceeding 1.0 mm, and head motion exceeding 2.0 mm of maximal translation (in any direction of x, y, or z) or 1.0° of maximal rotation throughout the course of scanning. We also despiked the fMRI data, regressed out six head motion parameters, white matter, and cerebrospinal fluid in the denoising procedure. Results indicate FD (mean framewise displacements, mean of root of mean square frame-to-frame head motions assuming 50 mm head radius) for all subjects were <0.3 mm at every time point. There is no significant difference between patients and controls on mean FDs; namely, UNM, HC: mean = 0.22 ± 0.12 mm, SZ: 0.21 ± 0.11 mm, two sample t-test: p = 0.77, FBIRN, HC: mean = 0.25 ± 0.18 mm, SZ: 0.27 ± 0.21 mm, two sample t-test: p = 0.65. No imaging voxels showed a significant correlation with mean FD after the FDR multiple comparison correction for any of the three modalities, and also no significant correlations between mean FD and cognition. Detailed head motion correction could be found in the Supplementary Note 6. Data were then spatially smoothed with an 8 mm, full-width half-maximum (FWHM) Gaussian filter. The sum of the amplitude values in 0.01 to 0.08 Hz low-frequency range was divided by the sum of the amplitudes over the entire power spectrum48 to obtain fALFF. Finally, the fusion analysis was conducted on the spatial maps of fALFF. Considering there is no group difference in head motion, no significant correlations between mean FD and cognitive scores, and the correlations with cognitive scores are still significant after regressing out mean FD, we believe that micro-motion is not a major factor affecting the current results.
The dMRI data were preprocessed using the FMRIB Software Library (www.fmrib.ox.ac.uk/fsl) and consisted of the following steps. Firstly, quality check and remove gradient directions with excessive motion or vibration artifact. Secondly, motion and eddy correction. Thirdly, correction of gradient directions for image rotation due to the motion correction procedure. Finally, FA was calculated and smoothed using an 8 mm FWHM Gaussian filter.
The sMRI data were spatially normalized to MNI space using the unified segmentation method in SPM8, segmented into GM, white matter, and cerebral spinal fluid (CSF). Then, the GM volume were smoothed with an FWHM of 8 mm Gaussian filter. Subject outliers were further detected by spatial Pearson correlation with the template image to ensure that all subjects were segmented properly.
After preprocessing, the three-dimensional brain images of each subject were reshaped into a one-dimensional vector and stacked, forming a matrix (Nsubj × Nvoxel) for each of the three modalities. These three matrices were then normalized to have the same average sum of squares (computed across all subjects and all voxels for each modality) to ensure all modalities had the same ranges. To ensure that each of the modality features were not confounded by site-related differences in subject recruitment criteria or by other unidentified variables, multivariate analysis of covariance (MANCOVA) was performed on all feature matrices. Site, gender, age, and their interactions were all regressed out from GM, fALFF, and FA features, respectively, to minimize their impact on the brain imaging data. Thus, the resulting data were then ready for fusion analysis.
The code for the supervised fusion algorithm has been released and integrated in the Fusion ICA Toolbox (FIT, https://mialab.mrn.org/software/fit), which can be downloaded freely and used directly by users worldwide. The multimodal data used in the present study can be accessed upon request to the corresponding authors.
Lewis, D. A. Cortical circuit dysfunction and cognitive deficits in schizophrenia—implications for preemptive interventions. Eur. J. Neurosci. 35, 1871–1878 (2012).
Dolan, R. J. Neuroimaging of cognition: past, present, and future. Neuron 60, 496–502 (2008).
Sui, J. et al. In search of multimodal neuroimaging biomarkers of cognitive deficits in schizophrenia. Biol. Psychiatry 78, 794–804 (2015).
Woo, C. W., Chang, L. J., Lindquist, M. A. & Wager, T. D. Building better biomarkers: brain models in translational neuroimaging. Nat. Neurosci. 20, 365–377 (2017).
Abi-Dargham, A. & Horga, G. The search for imaging biomarkers in psychiatric disorders. Nat. Med. 22, 1248–1255 (2016).
Rosenberg, M. D. et al. A neuromarker of sustained attention from whole-brain functional connectivity. Nat. Neurosci. 19, 165–171 (2016).
Finn, E. S. et al. Functional connectome fingerprinting: identifying individuals using patterns of brain connectivity. Nat. Neurosci. 18, 1664–1671 (2015).
Petersen Steven, E. & Sporns, O. Brain networks and cognitive architectures. Neuron 88, 207–219 (2015).
Green, M. F., Kern, R. S. & Heaton, R. K. Longitudinal studies of cognition and functional outcome in schizophrenia: implications for MATRICS. Schizophr. Res. 72, 41–51 (2004).
Castellanos, F. X., Di Martino, A., Craddock, R. C., Mehta, A. D. & Milham, M. P. Clinical applications of the functional connectome. Neuroimage 80, 527–540 (2013).
Thomas Insel, B. C. et al. Research domain criteria (RDoC): toward a new classification framework for research on mental disorders. Am. J. Psychiatry 167, 748–751 (2010).
Sui, J., Huster, R., Yu, Q., Segall, J. M. & Calhoun, V. D. Function-structure associations of the brain: evidence from multimodal connectivity and covariance studies. Neuroimage 102(Pt 1), 11–23 (2014).
Meng, X. et al. Predicting individualized clinical measures by a generalized prediction framework and multimodal fusion of MRI data. Neuroimage 145, 218–229 (2017).
Liu S. et al. Linked 4-way multimodal brain differences in schizophrenia in a large Chinese Han population. Schizophrenia Bull. https://doi.org/10.1093/schbul/sby045 (2018).
Calhoun, V. D. & Sui, J. Multimodal fusion of brain imaging data: a key to finding the missing link(s) in complex mental illness. Biol. Psychiatry. 1, 230–244 (2016).
Menon, V. Large-scale brain networks and psychopathology: a unifying triple network model. Trends Cogn. Sci. 15, 483–506 (2011).
Premkumar, P. et al. Association between a longer duration of illness, age and lower frontal lobe grey matter volume in schizophrenia. Behav. Brain. Res. 193, 132–139 (2008).
Delawalla, Z., Csernansky, J. G. & Barch, D. M. Prefrontal cortex function in nonpsychotic siblings of individuals with schizophrenia. Biol. Psychiatry 63, 490–497 (2008).
Aleman, A., Hijman, R., De Haan, E. H. F. & Kahn, R. S. Memory impairment in schizophrenia: a meta-analysis. Am. J. Psychiatry 156, 1358–1366 (1999).
Palaniyappan, L., Mallikarjun, P., Joseph, V., White, T. P. & Liddle, P. F. Regional contraction of brain surface area involves three large-scale networks in schizophrenia. Schizophr. Res. 129, 163–168 (2011).
Williamson, P. Are anticorrelated networks in the brain relevant to schizophrenia? Schizophr. Bull. 33, 994–1003 (2007).
Erhardt, E. B., Allen, E. A., Damaraju, E. & Calhoun, V. D. On network derivation, classification, and visualization: a response to Habeck and Moeller. Brain Connect. 1, 1–19 (2011).
Qi, S. et al. Multimodal fusion with reference: searching for joint neuromarkers of working memory deficits in schizophrenia. IEEE Trans. Med. Imaging 37, 93–105 (2018).
van Erp, T. G. et al. Neuropsychological profile in adult schizophrenia measured with the CMINDS. Psychiatry Res. 230, 826–834 (2015).
Aine, C. J. et al. Multimodal neuroimaging in schizophrenia: description and dissemination. Neuroinformatics 15, 343–364 (2017).
Qi, S. et al. MicroRNA132 associated multimodal neuroimaging patterns in unmedicated major depressive disorder. Brain 141, 916–926 (2018).
Li, Y. O., Adali, T. & Calhoun, V. D. Estimating the number of independent components for functional magnetic resonance imaging data. Hum. Brain. Mapp. 28, 1251–1266 (2007).
Kay, S. R. F. A. & Opler, L. A. The positive and negative syndrome scale (PANSS) for schizophrenia. Schizophr. Bull. 13, 261–276 (1987).
Kirkpatrick, B., Fenton, W. S., Carpenter, W. T. Jr. & Marder, S. R. The NIMH-MATRICS consensus statement on negative symptoms. Schizophr. Bull. 32, 214–219 (2006).
Harvey, P. D., Koren, D., Reichenberg, A. & Bowie, C. R. Negative symptoms and cognitive deficits: what is the nature of their relationship? Schizophr. Bull. 32, 250–258 (2006).
Chan, R. C., Shum, D., Toulopoulou, T. & Chen, E. Y. Assessment of executive functions: review of instruments and identification of critical issues. Arch. Clin. Neuropsychol. 23, 201–216 (2008).
Bora, E. et al. Neuroanatomical abnormalities in schizophrenia: a multimodal voxelwise meta-analysis and meta-regression analysis. Schizophr. Res. 127, 46–57 (2011).
Szeszko, P. R. et al. Clinical and neuropsychological correlates of white matter abnormalities in recent onset schizophrenia. Neuropsychopharmacology 33, 976–984 (2008).
Ehrlich, S. et al. Associations of white matter integrity and cortical thickness in patients with schizophrenia and healthy controls. Schizophr. Bull. 40, 665–674 (2014).
Potkin, S. G. et al. Working memory and DLPFC inefficiency in schizophrenia: the FBIRN study. Schizophr. Bull. 35, 19–31 (2009).
Glahn, D. C. et al. Beyond hypofrontality: a quantitative meta-analysis of functional neuroimaging studies of working memory in schizophrenia. Hum. Brain. Mapp. 25, 60–69 (2005).
Michael, J., Minzenberg, A. R. L., Thelen, S., Carter, S. S., & Glahn, D. C. Meta-analysis of 41 functional neuroimaging studies of executive function in schizophrenia. Arch. Gen. Psychiatry 66, 811–822 (2009).
Xu, X., Yuan, H. & Lei, X. Activation and connectivity within the default mode network contribute independently to future-oriented thought. Sci. Rep. 6, 21001 (2016).
Manoliu, A. et al. Aberrant dependence of default mode/central executive network interactions on anterior insular salience network activity in schizophrenia. Schizophr. Bull. 40, 428–437 (2014).
Sun, H. et al. Two patterns of white matter abnormalities in medication-naive patients with first-episode schizophrenia revealed by diffusion tensor imaging and cluster analysis. JAMA Psychiatry 72, 678–686 (2015).
Ren, W. et al. Anatomical and functional brain abnormalities in drug-naive first-episode schizophrenia. Am. J. Psychiatry 170, 1308–1316 (2013).
Huang, X. Q. et al. Localization of cerebral functional deficits in treatment-naive, first-episode schizophrenia using resting-state fMRI. Neuroimage 49, 2901–2906 (2010).
Turner, J. A. et al. A multi-site resting state fMRI study on the amplitude of low frequency fluctuations in schizophrenia. Front. Neurosci. 7, 137 (2013).
Chan, R. W. et al. Low-frequency hippocampal–cortical activity drives brain-wide resting-state functional MRI connectivity. Proc. Natl Acad. Sci. USA 114, E6972–E6981 (2017).
Menon, V. & Uddin, L. Q. Saliency, switching, attention and control: a network model of insula function. Brain Struct. Funct. 214, 655–667 (2010).
White, T. P., Joseph, V., Francis, S. T. & Liddle, P. F. Aberrant salience network (bilateral insula and anterior cingulate cortex) connectivity during information processing in schizophrenia. Schizophr. Res. 123, 105–115 (2010).
Sridharan, D., Levitin, D. J. & Menon, V. A critical role for the right fronto-insular cortex in switching between central-executive and default-mode networks. Proc. Natl Acad. Sci. USA 105, 12569–12574 (2008).
Koechlin, E. & Summerfield, C. An information theoretical approach to prefrontal executive function. Trends Cogn. Sci. 11, 229–235 (2007).
Muller, N. G. & Knight, R. T. The functional neuroanatomy of working memory: contributions of human brain lesion studies. Neuroscience 139, 51–58 (2006).
Tregellas J. R. et al. Intrinsic hippocampal activity as a biomarker for cognition and symptoms in schizophrenia. Am. J. Psychiatry 171, 549–556 (2014).
Dronkers, N. F., Plaisant, O., Iba-Zizen, M. T. & Cabanis, E. A. Paul Broca’s historic cases: high resolution MR imaging of the brains of Leborgne and Lelong. Brain 130, 1432–1441 (2007).
Andreasen, N. C. P. S. & O’Leary, D. S. “Cognitive dysmetria” as an integrative theory of schizophrenia: a dysfunction in cortical subcortical-cerebellar circuitry. Schizophr. Bull. 24, 203–218 (1998).
Goodkind, M. et al. Identification of a common neurobiological substrate for mental illness. JAMA Psychiatry 72, 305–315 (2015).
Downar, J., Blumberger, D. M. & Daskalakis, Z. J. The Neural crossroads of psychiatric illness: an emerging target for brain stimulation. Trends Cogn. Sci. 20, 107–120 (2016).
Kelly S. et al. Widespread white matter microstructural differences in schizophrenia across 4322 individuals: results from the ENIGMA Schizophrenia DTI Working Group. Mol. Psychiatry 21, 1710–1716 (2017).
Dosenbach, N. U. et al. Distinct brain networks for adaptive and stable task control in humans. Proc. Natl Acad. Sci. USA 104, 11073–11078 (2007).
Ehrlich, S. et al. Associations of cortical thickness and cognition in patients with schizophrenia and healthy controls. Schizophr. Bull. 38, 1050–1062 (2012).
Jiang, R. et al. SMRI biomarkers predict electroconvulsive treatment outcomes: accuracy with independent data sets. Neuropsychopharmacology 43, 1078–1087 (2018).
Calhoun, V. D. & Adali, T. Feature-based fusion of medical imaging data. IEEE Trans. Inf. Technol. Biomed. 13, 711–720 (2009).
Smith, S. M. et al. Correspondence of the brain’s functional architecture during activation and rest. Proc. Natl Acad. Sci. USA 106, 13040–13045 (2009).
Calhoun, V. D., Miller, R., Pearlson, G. & Adali, T. The chronnectome: time-varying connectivity networks as the next frontier in fMRI data discovery. Neuron 84, 262–274 (2014).
Ma, S., Calhoun, V. D., Phlypo, R. & Adali, T. Dynamic changes of spatial functional network connectivity in healthy individuals and schizophrenia patients using independent vector analysis. Neuroimage 90, 196–206 (2014).
First MBS, R. L., Gibbon, M. G., & Williams, J. B. W. Structured Clinical Interview for DSM-IV-TR Axis I Disorders Patient Edn (Biometrics Research, New York, 2002).
This work was supported by the Chinese National Natural Science Foundation Nos. 81471367, 61773380, and 61703253, the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDB02060005), the National High-Tech Development Plan (863 program, No. 2015AA020513), “100 Talents Plan” of Chinese Academy of Sciences, and the NIH grants R01EB005846, 1R01MH094524, and P20GM103472 as well as NSF grant 1539067.
These authors contributed equally: Jing Sui, Shile Qi.
J.S. and V.D.C. designed the study; S.Q. performed the data analysis; J.S., S.Q., and V.D.C. wrote the paper. J.B., T.G.M.v.E., J.A.T., D.H.M., J.M.F., J.V., B.A.M., A.B., S.M., S.G.P., and A.P. contributed to the multimodal imaging and cognitive data. E.D., Z.F., and R.J. helped with data preprocessing. All authors contributed to the results interpretation and discussion.
Correspondence to Jing Sui or Vince D. Calhoun.

References: V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V. 
 V.