Search

Formal thought disorder (FTD) is a clinical key factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, the relationship between FTD symptom dimensions and patterns of regional brain volume loss in schizophrenia remains to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles by enrolling a large multi-site cohort acquired by the ENIGMA Schizophrenia Working Group (752 schizophrenia patients and 1256 controls), to unravel the neuroanatomy of FTD in schizophrenia and using virtual histology tools on implicated brain regions to investigate the cellular basis. Based on the findings of previous clinical and neuroimaging studies, we decided to separately explore positive, negative and total formal thought disorder. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but positive and negative FTD demonstrated a dissociation: negative FTD showed a relative sparing of orbitofrontal cortical thickness, while positive FTD also affected lateral temporal cortices. Virtual histology identified distinct transcriptomic fingerprints associated for both symptom dimensions. Negative FTD was linked to neuronal and astrocyte fingerprints, while positive FTD also showed associations with microglial cell types. These results provide an important step towards linking FTD to brain structural changes and their cellular underpinnings, providing an avenue for a better mechanistic understanding of this syndrome.

Machine learning can be used to define subtypes of psychiatric conditions based on shared biological foundations of mental disorders. Here we analyzed cross-sectional brain images from 4,222 individuals with schizophrenia and 7038 healthy subjects pooled across 41 international cohorts from the ENIGMA, non-ENIGMA cohorts and public datasets. Using the Subtype and Stage Inference (SuStaIn) algorithm, we identify two distinct neurostructural subgroups by mapping the spatial and temporal trajectory of gray matter change in schizophrenia. Subgroup 1 was characterized by an early cortical-predominant loss with enlarged striatum, whereas subgroup 2 displayed an early subcortical-predominant loss in the hippocampus, striatum and other subcortical regions. We confirmed the reproducibility of the two neurostructural subtypes across various sample sites, including Europe, North America and East Asia. This imaging-based taxonomy holds the potential to identify individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Schizophrenia is a prototypical network disorder with widespread brain-morphological alterations, yet it remains unclear whether these distributed alterations robustly reflect the underlying network layout. We tested whether large-scale structural alterations in schizophrenia relate to normative structural and functional connectome architecture, and systematically evaluated robustness and generalizability of these network-level alterations. Leveraging anatomical MRI scans from 2439 adults with schizophrenia and 2867 healthy controls from 26 ENIGMA sites and normative data from the Human Connectome Project (n = 207), we evaluated structural alterations of schizophrenia against two network susceptibility models: (i) hub vulnerability, which examines associations between regional network centrality and magnitude of disease-related alterations; (ii) epicenter mapping, which identifies regions whose typical connectivity profile most closely resembles the disease-related morphological alterations. To assess generalizability and specificity, we contextualized the influence of site, disease stages, and individual clinical factors and compared network associations of schizophrenia with that found in affective disorders. Our findings show schizophrenia-related cortical thinning is spatially associated with functional and structural hubs, suggesting that highly interconnected regions are more vulnerable to morphological alterations. Predominantly temporo-paralimbic and frontal regions emerged as epicenters with connectivity profiles linked to schizophrenias alteration patterns. Findings were robust across sites, disease stages, and related to individual symptoms. Moreover, transdiagnostic comparisons revealed overlapping epicenters in schizophrenia and bipolar, but not major depressive disorder, suggestive of a pathophysiological continuity within the schizophrenia-bipolar-spectrum. In sum, cortical alterations over the course of schizophrenia robustly follow brain network architecture, emphasizing marked hub susceptibility and temporo-frontal epicenters at both the level of the group and the individual. Subtle variations of epicenters across disease stages suggest interacting pathological processes, while associations with patient-specific symptoms support additional inter-individual variability of hub vulnerability and epicenters in schizophrenia. Our work outlines potential pathways to better understand macroscale structural alterations, and inter- individual variability in schizophrenia.

BACKGROUND: Schizophrenia is associated with an increased risk of aggressive behaviour, which may partly be explained by illness-related changes in brain structure. However, previous studies have been limited by group-level analyses, small and selective samples of inpatients and long time lags between exposure and outcome. METHODS: This cross-sectional study pooled data from 20 sites participating in the international ENIGMA-Schizophrenia Working Group. Sites acquired T1-weighted and diffusion-weighted magnetic resonance imaging scans in a total of 2095 patients with schizophrenia and 2861 healthy controls. Measures of grey matter volume and white matter microstructural integrity were extracted from the scans using harmonised protocols. For each measure, normative modelling was used to calculate how much patients deviated (in z-scores) from healthy controls at the individual level. Ordinal regression models were used to estimate the associations of these deviations with concurrent aggressive behaviour (as odds ratios [ORs] with 99% confidence intervals [CIs]). Mediation analyses were performed for positive symptoms (i.e., delusions, hallucinations and disorganised thinking), impulse control and illness insight. Aggression and potential mediators were assessed with the Positive and Negative Syndrome Scale, Scale for the Assessment of Positive Symptoms or Brief Psychiatric Rating Scale. RESULTS: Aggressive behaviour was significantly associated with reductions in total cortical volume (OR [99% CI] = 0.88 [0.78, 0.98], p = .003) and global white matter integrity (OR [99% CI] = 0.72 [0.59, 0.88], p = 3.50 × 10-5) and additional reductions in dorsolateral prefrontal cortex volume (OR [99% CI] = 0.85 [0.74, 0.97], p =.002), inferior parietal lobule volume (OR [99% CI] = 0.76 [0.66, 0.87], p = 2.20 × 10-7) and internal capsule integrity (OR [99% CI] = 0.76 [0.63, 0.92], p = 2.90 × 10-4). Except for inferior parietal lobule volume, these associations were largely mediated by increased severity of positive symptoms and reduced impulse control. CONCLUSIONS: This study provides evidence that the co-occurrence of positive symptoms, poor impulse control and aggressive behaviour in schizophrenia has a neurobiological basis, which may inform the development of therapeutic interventions.

This work proposes a novel generative multimodal approach to jointly analyze multimodal data while linking the multimodal information to colors. We apply our proposed framework, which disentangles multimodal data into private and shared sets of features from pairs of structural (sMRI), functional (sFNC and ICA), and diffusion MRI data (FA maps). With our approach, we find that heterogeneity in schizophrenia is potentially a function of modality pairs. Results show (1) schizophrenia is highly multimodal and includes changes in specific networks, (2) non-linear relationships with schizophrenia are observed when interpolating among shared latent dimensions, and (3) we observe a decrease in the modularity of functional connectivity and decreased visual-sensorimotor connectivity for schizophrenia patients for the FA-sFNC and sMRI-sFNC modality pairs, respectively. Additionally, our results generally indicate decreased fractional corpus callosum anisotropy, and decreased spatial ICA map and voxel-based morphometry strength in the superior frontal lobe as found in the FA-sFNC, sMRI-FA, and sMRI-ICA modality pair clusters. In sum, we introduce a powerful new multimodal neuroimaging framework designed to provide a rich and intuitive understanding of the data which we hope challenges the reader to think differently about how modalities interact.

Machine learning can be used to define subtypes of psychiatric conditions based on shared clinical and biological foundations, presenting a crucial step toward establishing biologically based subtypes of mental disorders. With the goal of identifying subtypes of disease progression in schizophrenia, here we analyzed cross-sectional brain structural magnetic resonance imaging (MRI) data from 4,291 individuals with schizophrenia (1,709 females, age=32.5 years±11.9) and 7,078 healthy controls (3,461 females, age=33.0 years±12.7) pooled across 41 international cohorts from the ENIGMA Schizophrenia Working Group, non-ENIGMA cohorts and public datasets. Using a machine learning approach known as Subtype and Stage Inference (SuStaIn), we implemented a brain imaging-driven classification that identifies two distinct neurostructural subgroups by mapping the spatial and temporal trajectory of gray matter (GM) loss in schizophrenia. Subgroup 1 (n=2,622) was characterized by an early cortical-predominant loss (ECL) with enlarged striatum, whereas subgroup 2 (n=1,600) displayed an early subcortical-predominant loss (ESL) in the hippocampus, amygdala, thalamus, brain stem and striatum. These reconstructed trajectories suggest that the GM volume reduction originates in the Broca's area/adjacent fronto-insular cortex for ECL and in the hippocampus/adjacent medial temporal structures for ESL. With longer disease duration, the ECL subtype exhibited a gradual worsening of negative symptoms and depression/anxiety, and less of a decline in positive symptoms. We confirmed the reproducibility of these imaging-based subtypes across various sample sites, independent of macroeconomic and ethnic factors that differed across these geographic locations, which include Europe, North America and East Asia. These findings underscore the presence of distinct pathobiological foundations underlying schizophrenia. This new imaging-based taxonomy holds the potential to identify a more homogeneous sub-population of individuals with shared neurobiological attributes, thereby suggesting the viability of redefining existing disorder constructs based on biological factors.

Converging evidence suggests that schizophrenia (SZ) with primary, enduring negative symptoms (i.e., Deficit SZ (DSZ)) represents a distinct entity within the SZ spectrum while the neurobiological underpinnings remain undetermined. In the largest dataset of DSZ and Non-Deficit (NDSZ), we conducted a meta-analysis of data from 1560 individuals (168 DSZ, 373 NDSZ, 1019 Healthy Controls (HC)) and a mega-analysis of a subsampled data from 944 individuals (115 DSZ, 254 NDSZ, 575 HC) collected across 9 worldwide research centers of the ENIGMA SZ Working Group (8 in the mega-analysis), to clarify whether they differ in terms of cortical morphology. In the meta-analysis, sites computed effect sizes for differences in cortical thickness and surface area between SZ and control groups using a harmonized pipeline. In the mega-analysis, cortical values of individuals with schizophrenia and control participants were analyzed across sites using mixed-model ANCOVAs. The meta-analysis of cortical thickness showed a converging pattern of widespread thinner cortex in fronto-parietal regions of the left hemisphere in both DSZ and NDSZ, when compared to HC. However, DSZ have more pronounced thickness abnormalities than NDSZ, mostly involving the right fronto-parietal cortices. As for surface area, NDSZ showed differences in fronto-parietal-temporo-occipital cortices as compared to HC, and in temporo-occipital cortices as compared to DSZ. Although DSZ and NDSZ show widespread overlapping regions of thinner cortex as compared to HC, cortical thinning seems to better typify DSZ, being more extensive and bilateral, while surface area alterations are more evident in NDSZ. Our findings demonstrate for the first time that DSZ and NDSZ are characterized by different neuroimaging phenotypes, supporting a nosological distinction between DSZ and NDSZ and point toward the separate disease hypothesis.

Resting-state functional network connectivity (rsFNC) has shown utility for identifying characteristic functional brain patterns in individuals with psychiatric and mood disorders, providing a promising avenue for biomarker development. However, several factors have precluded widespread clinical adoption of rsFNC diagnostics, namely a lack of standardized approaches for capturing comparable and reproducible imaging markers across individuals, as well as the disagreement on the amount of data required to robustly detect intrinsic connectivity networks (ICNs) and diagnostically relevant patterns of rsFNC at the individual subject level. Recently, spatially constrained independent component analysis (scICA) has been proposed as an automated method for extracting ICNs standardized to a chosen network template while still preserving individual variation. Leveraging the scICA methodology, which solves the former challenge of standardized neuroimaging markers, we investigate the latter challenge of identifying a minimally sufficient data length for clinical applications of resting-state fMRI (rsfMRI). Using a dataset containing rsfMRI scans of individuals with schizophrenia and controls (M = 310) as well as simulated rsfMRI, we evaluated the robustness of ICN and rsFNC estimates at both the subject- and group-level, as well as the performance of diagnostic classification, with respect to the length of the rsfMRI time course. We found individual estimates of ICNs and rsFNC from the full-length (5 min) reference time course were sufficiently approximated with just 3-3.5 min of data (r = 0.85, 0.88, respectively), and significant differences in group-average rsFNC could be sufficiently approximated with even less data, just 2 min (r = 0.86). These results from the shorter clinical data were largely consistent with the results from validation experiments using longer time series from both simulated (30 min) and real-world (14 min) datasets, in which estimates of subject-level FNC were reliably estimated with 3-5 min of data. Moreover, in the real-world data we found rsFNC and ICN estimates generated across the full range of data lengths (0.5-14 min) more reliably matched those generated from the first 5 min of scan time than those generated from the last 5 min, suggesting increased influence of "late scan" noise factors such as fatigue or drowsiness may limit the reliability of FNC from data collected after 10+ min of scan time, further supporting the notion of shorter scans. Lastly, a diagnostic classification model trained on just 2 min of data retained 97%-98% classification accuracy relative to that of the full-length reference model. Our results suggest that, when decomposed with scICA, rsfMRI scans of just 2-5 min show good clinical utility without significant loss of individual FNC information of longer scan lengths.

Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

This study examined associations between resting-state amplitude of low frequency fluctuations (ALFF) and negative symptoms represented by total scores, second-order dimension (motivation and pleasure, expressivity), and first-order domain (anhedonia, avolition, asociality, alogia, blunted affect) factor scores in schizophrenia (n = 57). Total negative symptom scores showed positive associations with ALFF in temporal and frontal brain regions. Negative symptom domain scores showed predominantly stronger associations with regional ALFF compared to total scores, suggesting domain scores may better map to neural signatures than total scores. Improving our understanding of the neuropathology underlying negative symptoms may aid in addressing this unmet therapeutic need in schizophrenia.

Schizophrenia (SZ) is associated with an increased risk of life-long cognitive impairments, age-related chronic disease, and premature mortality. We investigated evidence for advanced brain ageing in adult SZ patients, and whether this was associated with clinical characteristics in a prospective meta-analytic study conducted by the ENIGMA Schizophrenia Working Group. The study included data from 26 cohorts worldwide, with a total of 2803 SZ patients (mean age 34.2 years; range 18-72 years; 67% male) and 2598 healthy controls (mean age 33.8 years, range 18-73 years, 55% male). Brain-predicted age was individually estimated using a model trained on independent data based on 68 measures of cortical thickness and surface area, 7 subcortical volumes, lateral ventricular volumes and total intracranial volume, all derived from T1-weighted brain magnetic resonance imaging (MRI) scans. Deviations from a healthy brain ageing trajectory were assessed by the difference between brain-predicted age and chronological age (brain-predicted age difference [brain-PAD]). On average, SZ patients showed a higher brain-PAD of +3.55 years (95% CI: 2.91, 4.19; I2 = 57.53%) compared to controls, after adjusting for age, sex and site (Cohen's d = 0.48). Among SZ patients, brain-PAD was not associated with specific clinical characteristics (age of onset, duration of illness, symptom severity, or antipsychotic use and dose). This large-scale collaborative study suggests advanced structural brain ageing in SZ. Longitudinal studies of SZ and a range of mental and somatic health outcomes will help to further evaluate the clinical implications of increased brain-PAD and its ability to be influenced by interventions.

IntroductionResting-state functional magnetic resonance imaging (rs-fMRI) is a powerful tool for assessing functional brain connectivity. Recent studies have focused on shorter-term connectivity and dynamics in the resting state. However, most of the prior work evaluates changes in time-series correlations. In this study, we propose a framework that focuses on time-resolved spectral coupling (assessed via the correlation between power spectra of the windowed time courses) among different brain circuits determined via independent component analysis (ICA).MethodsMotivated by earlier work suggesting significant spectral differences in people with schizophrenia, we developed an approach to evaluate time-resolved spectral coupling (trSC). To do this, we first calculated the correlation between the power spectra of windowed time-courses pairs of brain components. Then, we subgrouped each correlation map into four subgroups based on the connectivity strength utilizing quartiles and clustering techniques. Lastly, we examined clinical group differences by regression analysis for each averaged count and average cluster size matrices in each quartile. We evaluated the method by applying it to resting-state data collected from 151 (114 males, 37 females) people with schizophrenia (SZ) and 163 (117 males, 46 females) healthy controls (HC).ResultsOur proposed approach enables us to observe the change of connectivity strength within each quartile for different subgroups. People with schizophrenia showed highly modularized and significant differences in multiple network domains, whereas males and females showed less modular differences. Both cell count and average cluster size analysis for subgroups indicate a higher connectivity rate in the fourth quartile for the visual network in the control group. This indicates increased trSC in visual networks in the controls. In other words, this shows that the visual networks in people with schizophrenia have less mutually consistent spectra. It is also the case that the visual networks are less spectrally correlated on short timescales with networks of all other functional domains.ConclusionsThe results of this study reveal significant differences in the degree to which spectral power profiles are coupled over time. Importantly, there are significant but distinct differences both between males and females and between people with schizophrenia and controls. We observed a more significant coupling rate in the visual network for the healthy controls and males in the upper quartile. Fluctuations over time are complex, and focusing on only time-resolved coupling among time-courses is likely to miss important information. Also, people with schizophrenia are known to have impairments in visual processing but the underlying reasons for the impairment are still unknown. Therefore, the trSC approach can be a useful tool to explore the reasons for the impairments.

Brain functional networks identified from resting functional magnetic resonance imaging (fMRI) data have the potential to reveal biomarkers for brain disorders, but studies of complex mental illnesses such as schizophrenia (SZ) often yield mixed results across replication studies. This is likely due in part to the complexity of the disorder, the short data acquisition time, and the limited ability of the approaches for brain imaging data mining. Therefore, the use of analytic approaches which can both capture individual variability while offering comparability across analyses is highly preferred. Fully blind data-driven approaches such as independent component analysis (ICA) are hard to compare across studies, and approaches that use fixed atlas-based regions can have limited sensitivity to individual sensitivity. By contrast, spatially constrained ICA (scICA) provides a hybrid, fully automated solution that can incorporate spatial network priors while also adapting to new subjects. However, scICA has thus far only been used with a single spatial scale (ICA dimensionality, i.e., ICA model order). In this work, we present an approach using multi-objective optimization scICA with reference algorithm (MOO-ICAR) to extract subject-specific intrinsic connectivity networks (ICNs) from fMRI data at multiple spatial scales, which also enables us to study interactions across spatial scales. We evaluate this approach using a large N (N > 1,600) study of schizophrenia divided into separate validation and replication sets. A multi-scale ICN template was estimated and labeled, then used as input into scICA which was computed on an individual subject level. We then performed a subsequent analysis of multiscale functional network connectivity (msFNC) to evaluate the patient data, including group differences and classification. Results showed highly consistent group differences in msFNC in regions including cerebellum, thalamus, and motor/auditory networks. Importantly, multiple msFNC pairs linking different spatial scales were implicated. The classification model built on the msFNC features obtained up to 85% F1 score, 83% precision, and 88% recall, indicating the strength of the proposed framework in detecting group differences between schizophrenia and the control group. Finally, we evaluated the relationship of the identified patterns to positive symptoms and found consistent results across datasets. The results verified the robustness of our framework in evaluating brain functional connectivity of schizophrenia at multiple spatial scales, implicated consistent and replicable brain networks, and highlighted a promising approach for leveraging resting fMRI data for brain biomarker development.

In this work, we focus on explicitly nonlinear relationships in functional networks. We introduce a technique using normalized mutual information (NMI) that calculates the nonlinear relationship between different brain regions. We demonstrate our proposed approach using simulated data and then apply it to a dataset previously studied by Damaraju et al. This resting-state fMRI data included 151 schizophrenia patients and 163 age- and gender-matched healthy controls. We first decomposed these data using group independent component analysis (ICA) and yielded 47 functionally relevant intrinsic connectivity networks. Our analysis showed a modularized nonlinear relationship among brain functional networks that was particularly noticeable in the sensory and visual cortex. Interestingly, the modularity appears both meaningful and distinct from that revealed by the linear approach. Group analysis identified significant differences in explicitly nonlinear functional network connectivity (FNC) between schizophrenia patients and healthy controls, particularly in the visual cortex, with controls showing more nonlinearity (i.e., higher normalized mutual information between time courses with linear relationships removed) in most cases. Certain domains, including subcortical and auditory, showed relatively less nonlinear FNC (i.e., lower normalized mutual information), whereas links between the visual and other domains showed evidence of substantial nonlinear and modular properties. Overall, these results suggest that quantifying nonlinear dependencies of functional connectivity may provide a complementary and potentially important tool for studying brain function by exposing relevant variation that is typically ignored. Beyond this, we propose a method that captures both linear and nonlinear effects in a "boosted" approach. This method increases the sensitivity to group differences compared to the standard linear approach, at the cost of being unable to separate linear and nonlinear effects.

Background: While functional connectivity is widely studied, there has been little work studying functional connectivity at different spatial scales. Likewise, the relationship of functional connectivity between spatial scales is unknown. Methods: We proposed an independent component analysis (ICA)-based approach to capture information at multiple-model orders (component numbers), and to evaluate functional network connectivity (FNC) both within and between model orders. We evaluated the approach by studying group differences in the context of a study of resting-state functional magnetic resonance imaging (rsfMRI) data collected from schizophrenia (SZ) individuals and healthy controls (HC). The predictive ability of FNC at multiple spatial scales was assessed using support vector machine-based classification. Results: In addition to consistent predictive patterns at both multiple-model orders and single-model orders, unique predictive information was seen at multiple-model orders and in the interaction between model orders. We observed that the FNC between model orders 25 and 50 maintained the highest predictive information between HC and SZ. Results highlighted the predictive ability of the somatomotor and visual domains both within and between model orders compared with other functional domains. Also, subcortical-somatomotor, temporal-somatomotor, and temporal-subcortical FNCs had relatively high weights in predicting SZ. Conclusions: In sum, multimodel order ICA provides a more comprehensive way to study FNC, produces meaningful and interesting results, which are applicable to future studies. We shared the spatial templates from this work at different model orders to provide a reference for the community, which can be leveraged in regression-based or fully automated (spatially constrained) ICA approaches. Impact statement Multimodel order independent component analysis (ICA) provides a comprehensive way to study brain functional network connectivity within and between multiple spatial scales, highlighting findings that would have been ignored in single-model order analysis. This work expands upon and adds to the relatively new literature on resting functional magnetic resonance imaging-based classification and prediction. Results highlighted the differentiating power of specific intrinsic connectivity networks on classifying brain disorders of schizophrenia patients and healthy participants, at different spatial scales. The spatial templates from this work provide a reference for the community, which can be leveraged in regression-based or fully automated ICA approaches.

Revealing associations among various structural and functional patterns of the brain can yield highly informative results about the healthy and disordered brain. Studies using neuroimaging data have more recently begun to utilize the information within as well as across various functional and anatomical domains (i.e., groups of brain networks). However, most whole-brain approaches assume similar complexity of interactions throughout the brain. Here we investigate the hypothesis that interactions between brain networks capture varying amounts of complexity, and that we can better capture this information by varying the complexity of the model subspace structure based on available training data. To do this, we employ a Bayesian optimization-based framework known as the Tree Parzen Estimator (TPE) to identify, exploit and analyze patterns of variation in the information encoded by temporal information extracted from functional magnetic resonance imaging (fMRI) subdomains of the brain. Using a repeated cross-validation procedure on a schizophrenia classification task, we demonstrate evidence that interactions between specific functional subdomains are better characterized by more sophisticated model architectures compared to less complicated ones required by the others for optimally contributing towards classification and understanding the brain's functional interactions. We show that functional subdomains known to be involved in schizophrenia require more complex architectures to optimally unravel discriminatory information about the disorder. Our study points to the need for adaptive, hierarchical learning frameworks that cater differently to the features from different subdomains, not only for a better prediction but also for enabling the identification of features predicting the outcome of interest.

Graph-theoretical methods have been widely used to study human brain networks in psychiatric disorders. However, the focus has primarily been on global graphic metrics with little attention to the information contained in paths connecting brain regions. Details of disruption of these paths may be highly informative for understanding disease mechanisms. To detect the absence or addition of multistep paths in the patient group, we provide an algorithm estimating edges that contribute to these paths with reference to the control group. We next examine where pairs of nodes were connected through paths in both groups by using a covariance decomposition method. We apply our method to study resting-state fMRI data in schizophrenia versus controls. Results show several disconnectors in schizophrenia within and between functional domains, particularly within the default mode and cognitive control networks. Additionally, we identify new edges generating additional paths. Moreover, although paths exist in both groups, these paths take unique trajectories and have a significant contribution to the decomposition. The proposed path analysis provides a way to characterize individuals by evaluating changes in paths, rather than just focusing on the pairwise relationships. Our results show promise for identifying path-based metrics in neuroimaging data.

We introduce an extension of independent component analysis (ICA), called multiscale ICA, and design an approach to capture dynamic functional source interactions within and between multiple spatial scales. Multiscale ICA estimates functional sources at multiple spatial scales without imposing direct constraints on the size of functional sources, overcomes the limitation of using fixed anatomical locations, and eliminates the need for model-order selection in ICA analysis. We leveraged this approach to study sex-specific and sex-common connectivity patterns in schizophrenia. Results show dynamic reconfiguration and interaction within and between multi-spatial scales. Sex-specific differences occur (a) within the subcortical domain, (b) between the somatomotor and cerebellum domains, and (c) between the temporal domain and several others, including the subcortical, visual, and default mode domains. Most of the sex-specific differences belong to between-spatial-scale functional interactions and are associated with a dynamic state with strong functional interactions between the visual, somatomotor, and temporal domains and their anticorrelation patterns with the rest of the brain. We observed significant correlations between multi-spatial-scale functional interactions and symptom scores, highlighting the importance of multiscale analyses to identify potential biomarkers for schizophrenia. As such, we recommend such analyses as an important option for future functional connectivity studies.

N-methyl-D-aspartate receptor (NMDAR) hypofunction is a leading pathophysiological model of schizophrenia. Resting-state functional magnetic resonance imaging (rsfMRI) studies demonstrate a thalamic dysconnectivity pattern in schizophrenia involving excessive connectivity with sensory regions and deficient connectivity with frontal, cerebellar, and thalamic regions. The NMDAR antagonist ketamine, when administered at sub-anesthetic doses to healthy volunteers, induces transient schizophrenia-like symptoms and alters rsfMRI thalamic connectivity. However, the extent to which ketamine-induced thalamic dysconnectivity resembles schizophrenia thalamic dysconnectivity has not been directly tested. The current double-blind, placebo-controlled study derived an NMDAR hypofunction model of thalamic dysconnectivity from healthy volunteers undergoing ketamine infusions during rsfMRI. To assess whether ketamine-induced thalamic dysconnectivity was mediated by excess glutamate release, we tested whether pre-treatment with lamotrigine, a glutamate release inhibitor, attenuated ketamine's effects. Ketamine produced robust thalamo-cortical hyper-connectivity with sensory and motor regions that was not reduced by lamotrigine pre-treatment. To test whether the ketamine thalamic dysconnectivity pattern resembled the schizophrenia pattern, a whole-brain template representing ketamine's thalamic dysconnectivity effect was correlated with individual participant rsfMRI thalamic dysconnectivity maps, generating "ketamine similarity coefficients" for people with chronic (SZ) and early illness (ESZ) schizophrenia, individuals at clinical high-risk for psychosis (CHR-P), and healthy controls (HC). Similarity coefficients were higher in SZ and ESZ than in HC, with CHR-P showing an intermediate trend. Higher ketamine similarity coefficients correlated with greater hallucination severity in SZ. Thus, NMDAR hypofunction, modeled with ketamine, reproduces the thalamic hyper-connectivity observed in schizophrenia across its illness course, including the CHR-P period preceding psychosis onset, and may contribute to hallucination severity.

Schizophrenia is associated with widespread alterations in subcortical brain structure. While analytic methods have enabled more detailed morphometric characterization, findings are often equivocal. In this meta-analysis, we employed the harmonized ENIGMA shape analysis protocols to collaboratively investigate subcortical brain structure shape differences between individuals with schizophrenia and healthy control participants. The study analyzed data from 2,833 individuals with schizophrenia and 3,929 healthy control participants contributed by 21 worldwide research groups participating in the ENIGMA Schizophrenia Working Group. Harmonized shape analysis protocols were applied to each site's data independently for bilateral hippocampus, amygdala, caudate, accumbens, putamen, pallidum, and thalamus obtained from T1-weighted structural MRI scans. Mass univariate meta-analyses revealed more-concave-than-convex shape differences in the hippocampus, amygdala, accumbens, and thalamus in individuals with schizophrenia compared with control participants, more-convex-than-concave shape differences in the putamen and pallidum, and both concave and convex shape differences in the caudate. Patterns of exaggerated asymmetry were observed across the hippocampus, amygdala, and thalamus in individuals with schizophrenia compared to control participants, while diminished asymmetry encompassed ventral striatum and ventral and dorsal thalamus. Our analyses also revealed that higher chlorpromazine dose equivalents and increased positive symptom levels were associated with patterns of contiguous convex shape differences across multiple subcortical structures. Findings from our shape meta-analysis suggest that common neurobiological mechanisms may contribute to gray matter reduction across multiple subcortical regions, thus enhancing our understanding of the nature of network disorganization in schizophrenia.

The FAIR principles, as applied to clinical and neuroimaging data, reflect the goal of making research products Findable, Accessible, Interoperable, and Reusable. The use of the Collaborative Informatics and Neuroimaging Suite Toolkit for Anonymized Computation (COINSTAC) platform in the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) consortium combines the technological approach of decentralized analyses with the sociological approach of sharing data. In addition, ENIGMA + COINSTAC provides a platform to facilitate the use of machine-actionable data objects. We first present how ENIGMA and COINSTAC support the FAIR principles, and then showcase their integration with a decentralized meta-analysis of sex differences in negative symptom severity in schizophrenia, and finally present ongoing activities and plans to advance FAIR principles in ENIGMA + COINSTAC. ENIGMA and COINSTAC currently represent efforts toward improved Access, Interoperability, and Reusability. We highlight additional improvements needed in these areas, as well as future connections to other resources for expanded Findability.

BackgroundDysfunctional reward processing is implicated in multiple mental disorders. Novelty seeking (NS) assesses preference for seeking novel experiences, which is linked to sensitivity to reward environmental cues.MethodsA subset of 14-year-old adolescents (IMAGEN) with the top 20% ranked high-NS scores was used to identify high-NS-associated multimodal components by supervised fusion. These features were then used to longitudinally predict five different risk scales for the same and unseen subjects (an independent dataset of subjects at 19 years of age that was not used in predictive modeling training at 14 years of age) (within IMAGEN, n ≈1100) and even for the corresponding symptom scores of five types of patient cohorts (non-IMAGEN), including drinking (n = 313), smoking (n = 104), attention-deficit/hyperactivity disorder (n = 320), major depressive disorder (n = 81), and schizophrenia (n = 147), as well as to classify different patient groups with diagnostic labels.ResultsMultimodal biomarkers, including the prefrontal cortex, striatum, amygdala, and hippocampus, associated with high NS in 14-year-old adolescents were identified. The prediction models built on these features are able to longitudinally predict five different risk scales, including alcohol drinking, smoking, hyperactivity, depression, and psychosis for the same and unseen 19-year-old adolescents and even predict the corresponding symptom scores of five types of patient cohorts. Furthermore, the identified reward-related multimodal features can classify among attention-deficit/hyperactivity disorder, major depressive disorder, and schizophrenia with an accuracy of 87.2%.ConclusionsAdolescents with higher NS scores can be used to reveal brain alterations in the reward-related system, implicating potential higher risk for subsequent development of multiple disorders. The identified high-NS-associated multimodal reward-related signatures may serve as a transdiagnostic neuroimaging biomarker to predict disease risks or severity.

In this report, we present a case series involving four patients placed on the Clinical Institute Withdrawal Assessment for Alcohol, Revised (CIWA-Ar) protocol for alcohol or sedative-hypnotic withdrawal syndromes, who developed delirium on sustained or increasing symptom-triggered benzodiazepine dosages. In each of the four cases, delirium was not present on admission and resolved in the hospital itself with fixed benzodiazepine tapers. Cases were selected from an electronic medical record database of patients admitted to a United States-based university hospital and placed on CIWA-Ar between 2017 and 2018. This case series illustrates the major limitations of CIWA-Ar including its subjective nature, its susceptibility to inappropriate patient selection, and its requirement for providers to consider alternative etiologies to alcohol and benzodiazepine withdrawal syndromes. These cases demonstrate the necessity of considering other assessment and treatment options such as objective alcohol withdrawal scales, fixed benzodiazepine tapers, and even antiepileptics. An effective systems-based approach to overcoming these challenges may include setting time limits on CIWA-Ar orders within the electronic health record (EHR) system.

Background: Schizophrenia affects around 1% of the global population. Functional connectivity extracted from resting-state functional magnetic resonance imaging (rs-fMRI) has previously been used to study schizophrenia and has great potential to provide novel insights into the disorder. Some studies have shown abnormal functional connectivity in the default mode network (DMN) of individuals with schizophrenia, and more recent studies have shown abnormal dynamic functional connectivity (dFC) in individuals with schizophrenia. However, DMN dFC and the link between abnormal DMN dFC and symptom severity have not been well-characterized. Method: Resting-state fMRI data from subjects with schizophrenia (SZ) and healthy controls (HC) across two datasets were analyzed independently. We captured seven maximally independent subnodes in the DMN by applying group independent component analysis and estimated dFC between subnode time courses using a sliding window approach. A clustering method separated the dFCs into five reoccurring brain states. A feature selection method modeled the difference between SZs and HCs using the state-specific FC features. Finally, we used the transition probability of a hidden Markov model to characterize the link between symptom severity and dFC in SZ subjects. Results: We found decreases in the connectivity of the anterior cingulate cortex (ACC) and increases in the connectivity between the precuneus (PCu) and the posterior cingulate cortex (PCC) (i.e., PCu/PCC) of SZ subjects. In SZ, the transition probability from a state with weaker PCu/PCC and stronger ACC connectivity to a state with stronger PCu/PCC and weaker ACC connectivity increased with symptom severity. Conclusions: To our knowledge, this was the first study to investigate DMN dFC and its link to schizophrenia symptom severity. We identified reproducible neural states in a data-driven manner and demonstrated that the strength of connectivity within those states differed between SZs and HCs. Additionally, we identified a relationship between SZ symptom severity and the dynamics of DMN functional connectivity. We validated our results across two datasets. These results support the potential of dFC for use as a biomarker of schizophrenia and shed new light upon the relationship between schizophrenia and DMN dynamics.

BackgroundA number of studies in recent years have explored whole-brain dynamic connectivity using pairwise approaches. There has been less focus on trying to analyze brain dynamics in higher dimensions over time.MethodsWe introduce a new approach that analyzes time series trajectories to identify high traffic nodes in a high dimensional space. First, functional magnetic resonance imaging (fMRI) data are decomposed using spatial ICA to a set of maps and their associated time series. Next, density is calculated for each time point and high-density points are clustered to identify a small set of high traffic nodes. We validated our method using simulations and then implemented it on a real data set.ResultsWe present a novel approach that captures dynamics within a high dimensional space and also does not use any windowing in contrast to many existing approaches. The approach enables one to characterize and study the time series in a potentially high dimensional space, rather than looking at each component pair separately. Our results show that schizophrenia patients have a lower dynamism compared to healthy controls. In addition, we find patients spend more time in nodes associated with the default mode network and less time in components strongly correlated with auditory and sensorimotor regions. Interestingly, we also found that subjects oscillate between state pairs that show opposite spatial maps, suggesting an oscillatory pattern.ConclusionOur proposed method provides a novel approach to analyze the data in its native high dimensional space and can possibly provide new information that is undetectable using other methods.

BackgroundA number of studies in recent years have explored whole-brain dynamic connectivity using pairwise approaches. There has been less focus on trying to analyze brain dynamics in higher dimensions over time.MethodsWe introduce a new approach that analyzes time series trajectories to identify high traffic nodes in a high dimensional space. First, functional magnetic resonance imaging (fMRI) data are decomposed using spatial ICA to a set of maps and their associated time series. Next, density is calculated for each time point and high-density points are clustered to identify a small set of high traffic nodes. We validated our method using simulations and then implemented it on a real data set.ResultsWe present a novel approach that captures dynamics within a high dimensional space and also does not use any windowing in contrast to many existing approaches. The approach enables one to characterize and study the time series in a potentially high dimensional space, rather than looking at each component pair separately. Our results show that schizophrenia patients have a lower dynamism compared to healthy controls. In addition, we find patients spend more time in nodes associated with the default mode network and less time in components strongly correlated with auditory and sensorimotor regions. Interestingly, we also found that subjects oscillate between state pairs that show opposite spatial maps, suggesting an oscillatory pattern.ConclusionOur proposed method provides a novel approach to analyze the data in its native high dimensional space and can possibly provide new information that is undetectable using other methods.

Background: Schizophrenia affects around 1% of the global population. Functional connectivity extracted from resting-state functional magnetic resonance imaging (rs-fMRI) has previously been used to study schizophrenia and has great potential to provide novel insights into the disorder. Some studies have shown abnormal functional connectivity in the default mode network (DMN) of individuals with schizophrenia, and more recent studies have shown abnormal dynamic functional connectivity (dFC) in individuals with schizophrenia. However, DMN dFC and the link between abnormal DMN dFC and symptom severity have not been well-characterized. Method: Resting-state fMRI data from subjects with schizophrenia (SZ) and healthy controls (HC) across two datasets were analyzed independently. We captured seven maximally independent subnodes in the DMN by applying group independent component analysis and estimated dFC between subnode time courses using a sliding window approach. A clustering method separated the dFCs into five reoccurring brain states. A feature selection method modeled the difference between SZs and HCs using the state-specific FC features. Finally, we used the transition probability of a hidden Markov model to characterize the link between symptom severity and dFC in SZ subjects. Results: We found decreases in the connectivity of the anterior cingulate cortex (ACC) and increases in the connectivity between the precuneus (PCu) and the posterior cingulate cortex (PCC) (i.e., PCu/PCC) of SZ subjects. In SZ, the transition probability from a state with weaker PCu/PCC and stronger ACC connectivity to a state with stronger PCu/PCC and weaker ACC connectivity increased with symptom severity. Conclusions: To our knowledge, this was the first study to investigate DMN dFC and its link to schizophrenia symptom severity. We identified reproducible neural states in a data-driven manner and demonstrated that the strength of connectivity within those states differed between SZs and HCs. Additionally, we identified a relationship between SZ symptom severity and the dynamics of DMN functional connectivity. We validated our results across two datasets. These results support the potential of dFC for use as a biomarker of schizophrenia and shed new light upon the relationship between schizophrenia and DMN dynamics.

OBJECTIVE: Multimodal measurements of the same phenomena provide complementary information and highlight different perspectives, albeit each with their own limitations. A focus on a single modality may lead to incorrect inferences, which is especially important when a studied phenomenon is a disease. In this paper, we introduce a method that takes advantage of multimodal data in addressing the hypotheses of disconnectivity and dysfunction within schizophrenia (SZ). METHODS: We start with estimating and visualizing links within and among extracted multimodal data features using a Gaussian graphical model (GGM). We then propose a modularity-based method that can be applied to the GGM to identify links that are associated with mental illness across a multimodal data set. Through simulation and real data, we show our approach reveals important information about disease-related network disruptions that are missed with a focus on a single modality. We use functional MRI (fMRI), diffusion MRI (dMRI), and structural MRI (sMRI) to compute the fractional amplitude of low frequency fluctuations (fALFF), fractional anisotropy (FA), and gray matter (GM) concentration maps. These three modalities are analyzed using our modularity method. RESULTS: Our results show missing links that are only captured by the cross-modal information that may play an important role in disconnectivity between the components. CONCLUSION: We identified multimodal (fALFF, FA and GM) disconnectivity in the default mode network area in patients with SZ, which would not have been detectable in a single modality. SIGNIFICANCE: The proposed approach provides an important new tool for capturing information that is distributed among multiple imaging modalities.

ObjectiveMultimodal measurements of the same phenomena provide complementary information and highlight different perspectives, albeit each with their own limitations. A focus on a single modality may lead to incorrect inferences, which is especially important when a studied phenomenon is a disease. In this paper, we introduce a method that takes advantage of multimodal data in addressing the hypotheses of disconnectivity and dysfunction within schizophrenia (SZ).MethodsWe start with estimating and visualizing links within and among extracted multimodal data features using a Gaussian graphical model (GGM). We then propose a modularity-based method that can be applied to the GGM to identify links that are associated with mental illness across a multimodal data set. Through simulation and real data, we show our approach reveals important information about disease-related network disruptions that are missed with a focus on a single modality. We use functional MRI (fMRI), diffusion MRI (dMRI), and structural MRI (sMRI) to compute the fractional amplitude of low frequency fluctuations (fALFF), fractional anisotropy (FA), and gray matter (GM) concentration maps. These three modalities are analyzed using our modularity method.ResultsOur results show missing links that are only captured by the cross-modal information that may play an important role in disconnectivity between the components.ConclusionWe identified multimodal (fALFF, FA and GM) disconnectivity in the default mode network area in patients with SZ, which would not have been detectable in a single modality.SignificanceThe proposed approach provides an important new tool for capturing information that is distributed among multiple imaging modalities.

A common limitation of neuroimaging studies is their small sample sizes. To overcome this hurdle, the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Consortium combines neuroimaging data from many institutions worldwide. However, this introduces heterogeneity due to different scanning devices and sequences. ENIGMA projects commonly address this heterogeneity with random-effects meta-analysis or mixed-effects mega-analysis. Here we tested whether the batch adjustment method, ComBat, can further reduce site-related heterogeneity and thus increase statistical power. We conducted random-effects meta-analyses, mixed-effects mega-analyses and ComBat mega-analyses to compare cortical thickness, surface area and subcortical volumes between 2897 individuals with a diagnosis of schizophrenia and 3141 healthy controls from 33 sites. Specifically, we compared the imaging data between individuals with schizophrenia and healthy controls, covarying for age and sex. The use of ComBat substantially increased the statistical significance of the findings as compared to random-effects meta-analyses. The findings were more similar when comparing ComBat with mixed-effects mega-analysis, although ComBat still slightly increased the statistical significance. ComBat also showed increased statistical power when we repeated the analyses with fewer sites. Results were nearly identical when we applied the ComBat harmonization separately for cortical thickness, cortical surface area and subcortical volumes. Therefore, we recommend applying the ComBat function to attenuate potential effects of site in ENIGMA projects and other multi-site structural imaging work. We provide easy-to-use functions in R that work even if imaging data are partially missing in some brain regions, and they can be trained with one data set and then applied to another (a requirement for some analyses such as machine learning).

ObjectiveMultimodal measurements of the same phenomena provide complementary information and highlight different perspectives, albeit each with their own limitations. A focus on a single modality may lead to incorrect inferences, which is especially important when a studied phenomenon is a disease. In this paper, we introduce a method that takes advantage of multimodal data in addressing the hypotheses of disconnectivity and dysfunction within schizophrenia (SZ).MethodsWe start with estimating and visualizing links within and among extracted multimodal data features using a Gaussian graphical model (GGM). We then propose a modularity-based method that can be applied to the GGM to identify links that are associated with mental illness across a multimodal data set. Through simulation and real data, we show our approach reveals important information about disease-related network disruptions that are missed with a focus on a single modality. We use functional MRI (fMRI), diffusion MRI (dMRI), and structural MRI (sMRI) to compute the fractional amplitude of low frequency fluctuations (fALFF), fractional anisotropy (FA), and gray matter (GM) concentration maps. These three modalities are analyzed using our modularity method.ResultsOur results show missing links that are only captured by the cross-modal information that may play an important role in disconnectivity between the components.ConclusionWe identified multimodal (fALFF, FA and GM) disconnectivity in the default mode network area in patients with SZ, which would not have been detectable in a single modality.SignificanceThe proposed approach provides an important new tool for capturing information that is distributed among multiple imaging modalities.

A common limitation of neuroimaging studies is their small sample sizes. To overcome this hurdle, the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Consortium combines neuroimaging data from many institutions worldwide. However, this introduces heterogeneity due to different scanning devices and sequences. ENIGMA projects commonly address this heterogeneity with random-effects meta-analysis or mixed-effects mega-analysis. Here we tested whether the batch adjustment method, ComBat, can further reduce site-related heterogeneity and thus increase statistical power. We conducted random-effects meta-analyses, mixed-effects mega-analyses and ComBat mega-analyses to compare cortical thickness, surface area and subcortical volumes between 2897 individuals with a diagnosis of schizophrenia and 3141 healthy controls from 33 sites. Specifically, we compared the imaging data between individuals with schizophrenia and healthy controls, covarying for age and sex. The use of ComBat substantially increased the statistical significance of the findings as compared to random-effects meta-analyses. The findings were more similar when comparing ComBat with mixed-effects mega-analysis, although ComBat still slightly increased the statistical significance. ComBat also showed increased statistical power when we repeated the analyses with fewer sites. Results were nearly identical when we applied the ComBat harmonization separately for cortical thickness, cortical surface area and subcortical volumes. Therefore, we recommend applying the ComBat function to attenuate potential effects of site in ENIGMA projects and other multi-site structural imaging work. We provide easy-to-use functions in R that work even if imaging data are partially missing in some brain regions, and they can be trained with one data set and then applied to another (a requirement for some analyses such as machine learning).

The human brain is asymmetrically lateralized for certain functions (such as language processing) to regions in one hemisphere relative to the other. Asymmetries are measured with a laterality index (LI). However, traditional LI measures are limited by a lack of consensus on metrics used for its calculation. To address this limitation, source-based laterality (SBL) leverages an independent component analysis for the identification of laterality-specific alterations, identifying covarying components between hemispheres across subjects. SBL is successfully implemented with simulated data with inherent differences in laterality. SBL is then compared with a voxel-wise analysis utilizing structural data from a sample of patients with schizophrenia and controls without schizophrenia. SBL group comparisons identified three distinct temporal regions and one cerebellar region with significantly altered laterality in patients with schizophrenia relative to controls. Previous work highlights reductions in laterality (ie, reduced left gray matter volume) in patients with schizophrenia compared with controls without schizophrenia. Results from this pilot SBL project are the first, to our knowledge, to identify covarying laterality differences within discrete temporal brain regions. The authors argue SBL provides a unique focus to detect covarying laterality differences in patients with schizophrenia, facilitating the discovery of laterality aspects undetected in previous work.

The human brain is asymmetrically lateralized for certain functions (such as language processing) to regions in one hemisphere relative to the other. Asymmetries are measured with a laterality index (LI). However, traditional LI measures are limited by a lack of consensus on metrics used for its calculation. To address this limitation, source-based laterality (SBL) leverages an independent component analysis for the identification of laterality-specific alterations, identifying covarying components between hemispheres across subjects. SBL is successfully implemented with simulated data with inherent differences in laterality. SBL is then compared with a voxel-wise analysis utilizing structural data from a sample of patients with schizophrenia and controls without schizophrenia. SBL group comparisons identified three distinct temporal regions and one cerebellar region with significantly altered laterality in patients with schizophrenia relative to controls. Previous work highlights reductions in laterality (ie, reduced left gray matter volume) in patients with schizophrenia compared with controls without schizophrenia. Results from this pilot SBL project are the first, to our knowledge, to identify covarying laterality differences within discrete temporal brain regions. The authors argue SBL provides a unique focus to detect covarying laterality differences in patients with schizophrenia, facilitating the discovery of laterality aspects undetected in previous work.

The human brain is asymmetrically lateralized for certain functions (such as language processing) to regions in one hemisphere relative to the other. Asymmetries are measured with a laterality index (LI). However, traditional LI measures are limited by a lack of consensus on metrics used for its calculation. To address this limitation, source-based laterality (SBL) leverages an independent component analysis for the identification of laterality-specific alterations, identifying covarying components between hemispheres across subjects. SBL is successfully implemented with simulated data with inherent differences in laterality. SBL is then compared with a voxel-wise analysis utilizing structural data from a sample of patients with schizophrenia and controls without schizophrenia. SBL group comparisons identified three distinct temporal regions and one cerebellar region with significantly altered laterality in patients with schizophrenia relative to controls. Previous work highlights reductions in laterality (ie, reduced left gray matter volume) in patients with schizophrenia compared with controls without schizophrenia. Results from this pilot SBL project are the first, to our knowledge, to identify covarying laterality differences within discrete temporal brain regions. The authors argue SBL provides a unique focus to detect covarying laterality differences in patients with schizophrenia, facilitating the discovery of laterality aspects undetected in previous work.

Negative symptoms are core contributors to vocational and social deficits in schizophrenia (SZ). Available antipsychotic medications typically fail to reduce these symptoms. The neurohormone oxytocin (OT) is a promising treatment for negative symptoms, given its role in complex social behaviors mediated by the amygdala. In sample 1, we used a double-blind, placebo-controlled, crossover design to test the effects of a single dose of intranasal OT on amygdala resting-state functional connectivity (rsFC) in SZ (n = 22) and healthy controls (HC, n = 24) using a whole-brain corrected approach: we identified regions for which OT modulated SZ amygdala rsFC, assessed whether OT-modulated circuits were abnormal in SZ relative to HC on placebo, and evaluated whether connectivity on placebo and OT-induced connectivity changes correlated with baseline negative symptoms in SZ. Given our modest sample size, we used a second SZ (n = 183) and HC (n = 178) sample to replicate any symptom correlations. In sample 1, OT increased rsFC between the amygdala and left middle temporal gyrus, superior temporal sulcus, and angular gyrus (MTG/STS/AngG) in SZ compared to HC. Further, SZ had hypo-connectivity in this circuit compared to HC on placebo. More severe negative symptoms correlated with less amygdala-to-left-MTG/STS/AngG connectivity on placebo and with greater OT-induced connectivity increases. In sample 2, we replicated the correlation between amygdala-left-MTG/STS/AngG hypo-connectivity and negative symptoms, finding a specific association with expressive negative symptoms. These data suggest intranasal OT can normalize functional connectivity in an amygdala-to-left-MTG/STS/AngG circuit that contributes to negative symptoms in SZ.

BackgroundDynamic functional network connectivity (dFNC) of the brain has attracted considerable attention recently. Many approaches have been suggested to study dFNC with sliding window Pearson correlation (SWPC) being the most well-known. SWPC needs a relatively large sample size to reach a robust estimation but using large window sizes prevents us to detect rapid changes in dFNC.New methodHere we first calculate the gradients of each time series pair and use the magnitude of these gradients to calculate weighted average of shared trajectory (WAST) as a new estimator for dFNC.ResultsUsing WAST to compare healthy control and schizophrenia patients using a large dataset, we show disconnectivity between different regions associated with schizophrenia. In addition, WAST results reveals patients with schizophrenia stay longer in a connectivity state with negative connectivity between motor and sensory regions than do healthy controls.Comparison with existing methodsWe compare WAST with SWPC and multiplication of temporal derivatives (MTD) using different simulation scenarios. We show that WAST enables us to detect very rapid changes in dFNC (undetected by SWPC) while MTD performance is generally lower.ConclusionsAs large window sizes are unable to detect short states, using shorter window size is desirable if the estimator is robust enough. We provide evidence that WAST requires fewer samples (compared to SWPC) to reach a robust estimation. As a result, we were able to identify rapidly varying dFNC patterns undetected by SWPC while still being able to robustly estimate slower dFNC patterns.

There is growing evidence that rather than using a single brain imaging modality to study its association with physiological or symptomatic features, the field is paying more attention to fusion of multimodal information. However, most current multimodal fusion approaches that incorporate functional magnetic resonance imaging (fMRI) are restricted to second-level 3D features, rather than the original 4D fMRI data. This trade-off is that the valuable temporal information is not utilized during the fusion step. Here we are motivated to propose a novel approach called "parallel group ICA+ICA" that incorporates temporal fMRI information from group independent component analysis (GICA) into a parallel independent component analysis (ICA) framework, aiming to enable direct fusion of first-level fMRI features with other modalities (e.g., structural MRI), which thus can detect linked functional network variability and structural covariations. Simulation results show that the proposed method yields accurate intermodality linkage detection regardless of whether it is strong or weak. When applied to real data, we identified one pair of significantly associated fMRI-sMRI components that show group difference between schizophrenia and controls in both modalities, and this linkage can be replicated in an independent cohort. Finally, multiple cognitive domain scores can be predicted by the features identified in the linked component pair by our proposed method. We also show these multimodal brain features can predict multiple cognitive scores in an independent cohort. Overall, results demonstrate the ability of parallel GICA+ICA to estimate joint information from 4D and 3D data without discarding much of the available information up front, and the potential for using this approach to identify imaging biomarkers to study brain disorders.

There is growing evidence that rather than using a single brain imaging modality to study its association with physiological or symptomatic features, the field is paying more attention to fusion of multimodal information. However, most current multimodal fusion approaches that incorporate functional magnetic resonance imaging (fMRI) are restricted to second-level 3D features, rather than the original 4D fMRI data. This trade-off is that the valuable temporal information is not utilized during the fusion step. Here we are motivated to propose a novel approach called "parallel group ICA+ICA" that incorporates temporal fMRI information from group independent component analysis (GICA) into a parallel independent component analysis (ICA) framework, aiming to enable direct fusion of first-level fMRI features with other modalities (e.g., structural MRI), which thus can detect linked functional network variability and structural covariations. Simulation results show that the proposed method yields accurate intermodality linkage detection regardless of whether it is strong or weak. When applied to real data, we identified one pair of significantly associated fMRI-sMRI components that show group difference between schizophrenia and controls in both modalities, and this linkage can be replicated in an independent cohort. Finally, multiple cognitive domain scores can be predicted by the features identified in the linked component pair by our proposed method. We also show these multimodal brain features can predict multiple cognitive scores in an independent cohort. Overall, results demonstrate the ability of parallel GICA+ICA to estimate joint information from 4D and 3D data without discarding much of the available information up front, and the potential for using this approach to identify imaging biomarkers to study brain disorders.

The brain is highly dynamic, reorganizing its activity at different interacting spatial and temporal scales, including variation within and between brain networks. The chronnectome is a model of the brain in which nodal activity and connectivity patterns change in fundamental and recurring ways over time. Most literature assumes fixed spatial nodes/networks, ignoring the possibility that spatial nodes/networks may vary in time. Here, we introduce an approach to calculate a spatially fluid chronnectome (called the spatial chronnectome for clarity), which focuses on the variations of networks coupling at the voxel level, and identify a novel set of spatially dynamic features. Results reveal transient spatially fluid interactions between intra- and internetwork relationships in which brain networks transiently merge and separate, emphasizing dynamic segregation and integration. Brain networks also exhibit distinct spatial patterns with unique temporal characteristics, potentially explaining a broad spectrum of inconsistencies in previous studies that assumed static networks. Moreover, we show anticorrelative connections to brain networks are transient as opposed to constant across the entire scan. Preliminary assessments using a multi-site dataset reveal the ability of the approach to obtain new information and nuanced alterations that remain undetected during static analysis. Patients with schizophrenia (SZ) display transient decreases in voxel-wise network coupling within visual and auditory networks, and higher intradomain coupling variability. In summary, the spatial chronnectome represents a new direction of research enabling the study of functional networks which are transient at the voxel level, and the identification of mechanisms for within- and between-subject spatial variability.

Schizophrenia is a complex, debilitating mental disorder characterized by wide-ranging symptoms including delusions, hallucinations (so-called positive symptoms), and impaired motor and speech/language production (so-called negative symptoms). Salience-monitoring theorists propose that abnormal functional communication between the salience network (SN) and default mode network (DMN) begets positive and negative symptoms of schizophrenia, yet prior studies have predominately reported links between disrupted SN/DMN functional communication and positive symptoms. It remains unclear whether disrupted SN/DMN functional communication explains (1) solely positive symptoms or (2) both positive and negative symptoms of schizophrenia. To address this question, we incorporate time-lag-shifted functional network connectivity (FNC) analyses that explored coherence of the resting-state functional magnetic resonance imaging signal of 3 networks (anterior DMN, posterior DMN, and SN) with fixed time lags introduced between network time series (1 TR = 2 s; 2 TR = 4 s). Multivariate linear regression analysis revealed that severity of disordered thought and attentional deficits were negatively associated with 2 TR-shifted FNC between anterior DMN and posterior DMN. Meanwhile, severity of flat affect and bizarre behavior were positively associated with 1 TR-shifted FNC between anterior DMN and SN. These results provide support favoring the hypothesis that lagged SN/DMN functional communication is associated with both positive and negative symptoms of schizophrenia.

It has been widely shown that genomic factors influence both risk for schizophrenia and variation in functional brain connectivity. Moreover, schizophrenia is characterized by disrupted brain connectivity. In this work, we proposed a genome-connectome bipartite graph model to perform imaging genomic analysis. Functional network connectivity (FNC) was estimated after decomposing resting state functional magnetic resonance imaging data from both healthy controls (HC) and patients with schizophrenia (SZ) into spatial brain components using group independent component analysis (G-ICA). Then 83 FNC connections showing a group difference (HC vs SZ) were selected as fMRI nodes, and eighty-one schizophrenia-related single nucleotide polymorphisms (SNPs) were selected as genetic nodes respectively in the bipartite graph. Edges connecting pairs of genetic and fMRI nodes were defined based on the SNP-FNC associations across subjects evaluated by a general linear model. Results show that some SNP nodes in the bipartite graph have a high degree implying they are influential in modulating brain connectivity and may be more strongly associated with the risk of schizophrenia than other SNPs. A bi-clustering analysis detected a cluster with 15 SNPs interacting with 38 FNC connections, most of which were within or between somato-motor and visual brain areas. This suggests that the activity of these brain regions may be related to common SNPs and provides insights into the pathology of schizophrenia. The findings suggest that the SNP-FNC bipartite graph approach is a novel model to investigate genetic influences on functional brain connectivity in mental illness.

The analysis of time-varying activity and connectivity patterns (i.e., the chronnectome) using resting-state magnetic resonance imaging has become an important part of ongoing neuroscience discussions. The majority of previous work has focused on variations of temporal coupling among fixed spatial nodes or transition of the dominant activity/connectivity pattern over time. Here, we introduce an approach to capture spatial dynamics within functional domains (FDs), as well as temporal dynamics within and between FDs. The approach models the brain as a hierarchical functional architecture with different levels of granularity, where lower levels have higher functional homogeneity and less dynamic behavior and higher levels have less homogeneity and more dynamic behavior. First, a high-order spatial independent component analysis is used to approximate functional units. A functional unit is a pattern of regions with very similar functional activity over time. Next, functional units are used to construct FDs. Finally, functional modules (FMs) are calculated from FDs, providing an overall view of brain dynamics. Results highlight the spatial fluidity within FDs, including a broad spectrum of changes in regional associations, from strong coupling to complete decoupling. Moreover, FMs capture the dynamic interplay between FDs. Patients with schizophrenia show transient reductions in functional activity and state connectivity across several FDs, particularly the subcortical domain. Activity and connectivity differences convey unique information in many cases (e.g., the default mode) highlighting their complementarity information. The proposed hierarchical model to capture FD spatiotemporal variations provides new insight into the macroscale chronnectome and identifies changes hidden from existing approaches.

Functional connectivity is one of the most widely used tools for investigating brain changes due to schizophrenia. Previous studies have identified abnormal functional connectivity in schizophrenia patients at the resting state brain network level. This study tests the existence of functional connectivity effects at whole brain and domain levels. Domain level refers to the integration of data from several brain networks grouped by their functional relationship. Data integration provides more consistent and accurate information compared to an individual brain network. This work considers two domain level measures: functional connectivity strength and randomness. The first measure is simply an average of connectivities within the domain. The second measure assesses the unpredictability and lack of pattern of functional connectivity within the domain. Domains with less random connectivity have higher chance of exhibiting a biologically meaningful connectivity pattern. Consistent with prior observations, individuals with schizophrenia showed aberrant domain connectivity strength between subcortical, cerebellar, and sensorial brain areas. Compared to healthy volunteers, functional connectivity between cognitive and default mode domains showed less randomness, while connectivity between default mode-sensorial areas showed more randomness in schizophrenia patients. These differences in connectivity patterns suggest deleterious rewiring trade-offs among important brain networks.

Multimodal, imaging-genomics techniques offer a platform for understanding genetic influences on brain abnormalities in psychiatric disorders. Such approaches utilize the information available from both imaging and genomics data and identify their association. Particularly for complex disorders such as schizophrenia, the relationship between imaging and genomic features may be better understood by incorporating additional information provided by advanced multimodal modeling. In this study, we propose a novel framework to combine features corresponding to functional magnetic resonance imaging (functional) and single nucleotide polymorphism (SNP) data from 61 schizophrenia (SZ) patients and 87 healthy controls (HC). In particular, the features for the functional and genetic modalities include dynamic (i.e., time-varying) functional network connectivity (dFNC) features and the SNP data, respectively. The dFNC features are estimated from component time-courses, obtained using group independent component analysis (ICA), by computing sliding-window functional network connectivity, and then estimating subject specific states from this dFNC data using a k-means clustering approach. For each subject, both the functional (dFNC states) and SNP data are selected as features for a parallel ICA (pICA) based imaging-genomic framework. This analysis identified a significant association between a SNP component (defined by large clusters of functionally related SNPs statistically correlated with phenotype components) and time-varying or dFNC component (defined by clusters of related connectivity links among distant brain regions distributed across discrete dynamic states, and statistically correlated with genomic components) in schizophrenia. Importantly, the polygenetic risk score (PRS) for SZ (computed as a linearly weighted sum of the genotype profiles with weights derived from the odds ratios of the psychiatric genomics consortium (PGC)) was negatively correlated with the significant dFNC component, which were mostly present within a state that exhibited a lower occupancy rate in individuals with SZ compared with HC, hence identifying a potential dFNC imaging biomarker for schizophrenia. Taken together, the current findings provide preliminary evidence for a link between dFNC measures and genetic risk, suggesting the application of dFNC patterns as biomarkers in imaging genetic association study.

Brain functional connectivity has been shown to change over time during resting state fMRI experiments. Close examination of temporal changes have revealed a small set of whole-brain connectivity patterns called dynamic states. Dynamic functional network connectivity (dFNC) studies have demonstrated that it is possible to replicate the dynamic states across several resting state experiments. However, estimation of states and their temporal dynamicity still suffers from noisy and imperfect estimations. In regular dFNC implementations, states are estimated by comparing connectivity patterns through the data without considering time, in other words only zero order changes are examined. In this work we propose a method that includes first order variations of dFNC in the searching scheme of dynamic connectivity patterns. Our approach, referred to as temporal variation of functional network connectivity (tvFNC), estimates the derivative of dFNC, and then searches for reoccurring patterns of concurrent dFNC states and their derivatives. The tvFNC method is first validated using a simulated dataset and then applied to a resting-state fMRI sample including healthy controls (HC) and schizophrenia (SZ) patients and compared to the standard dFNC approach. Our dynamic approach reveals extra patterns in the connectivity derivatives complementing the already reported state patterns. State derivatives consist of additional information about increment and decrement of connectivity among brain networks not observed by the original dFNC method. The tvFNC shows more sensitivity than regular dFNC by uncovering additional FNC differences between the HC and SZ groups in each state. In summary, the tvFNC method provides a new and enhanced approach to examine time-varying functional connectivity.

Functional connectivity is one of the most widely used tools for investigating brain changes due to schizophrenia. Previous studies have identified abnormal functional connectivity in schizophrenia patients at the resting state brain network level. This study tests the existence of functional connectivity effects at whole brain and domain levels. Domain level refers to the integration of data from several brain networks grouped by their functional relationship. Data integration provides more consistent and accurate information compared to an individual brain network. This work considers two domain level measures: functional connectivity strength and randomness. The first measure is simply an average of connectivities within the domain. The second measure assesses the unpredictability and lack of pattern of functional connectivity within the domain. Domains with less random connectivity have higher chance of exhibiting a biologically meaningful connectivity pattern. Consistent with prior observations, individuals with schizophrenia showed aberrant domain connectivity strength between subcortical, cerebellar, and sensorial brain areas. Compared to healthy volunteers, functional connectivity between cognitive and default mode domains showed less randomness, while connectivity between default mode-sensorial areas showed more randomness in schizophrenia patients. These differences in connectivity patterns suggest deleterious rewiring trade-offs among important brain networks.

Brain functional connectivity has been shown to change over time during resting state fMRI experiments. Close examination of temporal changes have revealed a small set of whole-brain connectivity patterns called dynamic states. Dynamic functional network connectivity (dFNC) studies have demonstrated that it is possible to replicate the dynamic states across several resting state experiments. However, estimation of states and their temporal dynamicity still suffers from noisy and imperfect estimations. In regular dFNC implementations, states are estimated by comparing connectivity patterns through the data without considering time, in other words only zero order changes are examined. In this work we propose a method that includes first order variations of dFNC in the searching scheme of dynamic connectivity patterns. Our approach, referred to as temporal variation of functional network connectivity (tvFNC), estimates the derivative of dFNC, and then searches for reoccurring patterns of concurrent dFNC states and their derivatives. The tvFNC method is first validated using a simulated dataset and then applied to a resting-state fMRI sample including healthy controls (HC) and schizophrenia (SZ) patients and compared to the standard dFNC approach. Our dynamic approach reveals extra patterns in the connectivity derivatives complementing the already reported state patterns. State derivatives consist of additional information about increment and decrement of connectivity among brain networks not observed by the original dFNC method. The tvFNC shows more sensitivity than regular dFNC by uncovering additional FNC differences between the HC and SZ groups in each state. In summary, the tvFNC method provides a new and enhanced approach to examine time-varying functional connectivity.

AIMS:To provide quantitative conversions between commonly used scales for the assessment of negative symptoms in schizophrenia. METHOD:Linear regression analyses generated conversion equations between symptom scores from the Scale for the Assessment of Negative Symptoms (SANS), the Schedule for the Deficit Syndrome (SDS), the Positive and Negative Syndrome Scale (PANSS), or the Negative Symptoms Assessment (NSA) based on a cross sectional sample of 176 individuals with schizophrenia. Intraclass correlations assessed the rating conversion accuracy based on a separate sub-sample of 29 patients who took part in the initial study as well as an independent sample of 28 additional subjects with schizophrenia. RESULTS:Between-scale negative symptom ratings were moderately to highly correlated (r = 0.73-0.91). Intraclass correlations between the original negative symptom rating scores and those obtained via using the conversion equations were in the range of 0.61-0.79. CONCLUSIONS:While there is a degree of non-overlap, several negative symptoms scores reflect measures of similar constructs and may be reliably converted between some scales. The conversion equations are provided at http://www.converteasy.org and may be used for meta- and mega-analyses that examine negative symptoms.

BackgroundThe profile of cortical neuroanatomical abnormalities in schizophrenia is not fully understood, despite hundreds of published structural brain imaging studies. This study presents the first meta-analysis of cortical thickness and surface area abnormalities in schizophrenia conducted by the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) Schizophrenia Working Group.MethodsThe study included data from 4474 individuals with schizophrenia (mean age, 32.3 years; range, 11-78 years; 66% male) and 5098 healthy volunteers (mean age, 32.8 years; range, 10-87 years; 53% male) assessed with standardized methods at 39 centers worldwide.ResultsCompared with healthy volunteers, individuals with schizophrenia have widespread thinner cortex (left/right hemisphere: Cohen's d = -0.530/-0.516) and smaller surface area (left/right hemisphere: Cohen's d = -0.251/-0.254), with the largest effect sizes for both in frontal and temporal lobe regions. Regional group differences in cortical thickness remained significant when statistically controlling for global cortical thickness, suggesting regional specificity. In contrast, effects for cortical surface area appear global. Case-control, negative, cortical thickness effect sizes were two to three times larger in individuals receiving antipsychotic medication relative to unmedicated individuals. Negative correlations between age and bilateral temporal pole thickness were stronger in individuals with schizophrenia than in healthy volunteers. Regional cortical thickness showed significant negative correlations with normalized medication dose, symptom severity, and duration of illness and positive correlations with age at onset.ConclusionsThe findings indicate that the ENIGMA meta-analysis approach can achieve robust findings in clinical neuroscience studies; also, medication effects should be taken into account in future genetic association studies of cortical thickness in schizophrenia.