Structural white matter abnormalities in Schizophrenia and associations with neurocognitive performance and symptom severity.

Schizophrenia is associated with robust white matter (WM) abnormalities but influences of potentially confounding variables and relationships with cognitive performance and symptom severity remain to be fully determined. This study was designed to evaluate WM abnormalities based on diffusion tensor imaging (DTI) in individuals with schizophrenia, and their relationships with cognitive performance and symptom severity. Data from individuals with schizophrenia (SZ; n=138, mean age±SD=39.02±11.82; 105 males) and healthy controls (HC; n=143, mean age±SD=37.07±10.84; 102 males) were collected as part of the Function Biomedical Informatics Research Network Phase 3 study. Fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD), and mean diffusivity (MD) were compared between individuals with schizophrenia and healthy controls, and their relationships with neurocognitive performance and symptomatology assessed. Individuals with SZ had significantly lower FA in forceps minor and the left inferior fronto-occipital fasciculus compared to HC. FA in several tracts were associated with speed of processing and attention/vigilance and the severity of the negative symptom alogia. This study suggests that regional WM abnormalities are fundamentally involved in the pathophysiology of schizophrenia and may contribute to cognitive performance deficits and symptom expression observed in schizophrenia.

Genetic variants for head size share genes and pathways with cancer.

The size of the human head is highly heritable, but genetic drivers of its variation within the general population remain unmapped. We perform a genome-wide association study on head size (N = 80,890) and identify 67 genetic loci, of which 50 are novel. Neuroimaging studies show that 17 variants affect specific brain areas, but most have widespread effects. Gene set enrichment is observed for various cancers and the p53, Wnt, and ErbB signaling pathways. Genes harboring lead variants are enriched for macrocephaly syndrome genes (37-fold) and high-fidelity cancer genes (9-fold), which is not seen for human height variants. Head size variants are also near genes preferentially expressed in intermediate progenitor cells, neural cells linked to evolutionary brain expansion. Our results indicate that genes regulating early brain and cranial growth incline to neoplasia later in life, irrespective of height. This warrants investigation of clinical implications of the link between head size and cancer.

Differences in the neural correlates of schizophrenia with positive and negative formal thought disorder in patients with schizophrenia in the ENIGMA dataset.

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.

A method for estimating dynamic functional network connectivity gradients (dFNG) from ICA captures smooth inter-network modulation.

Dynamic functional network connectivity (dFNC) analysis is a widely used approach for studying brain function and offering insight into how brain networks evolve over time. Typically, dFNC studies utilized fixed spatial maps and evaluate transient changes in coupling among time courses estimated from independent component analysis (ICA). This manuscript presents a complementary approach that relaxes this assumption by spatially reordering the components dynamically at each timepoint to optimize for a smooth gradient in the FNC (i.e., a smooth gradient among ICA connectivity values). Several methods are presented to summarize dynamic FNC gradients (dFNGs) over time, starting with static FNC gradients (sFNGs), then exploring the reordering properties as well as the dynamics of the gradients themselves. We then apply this approach to a dataset of schizophrenia (SZ) patients and healthy controls (HC). Functional dysconnectivity between different brain regions has been reported in schizophrenia, yet the neural mechanisms behind it remain elusive. Using resting state fMRI and ICA on a dataset consisting of 151 schizophrenia patients and 160 age and gender-matched healthy controls, we extracted 53 intrinsic connectivity networks (ICNs) for each subject using a fully automated spatially constrained ICA approach. We develop several summaries of our functional network connectivity gradient analysis, both in a static sense, computed as the Pearson correlation coefficient between full time series, and a dynamic sense, computed using a sliding window approach followed by reordering based on the computed gradient, and evaluate group differences. Static connectivity analysis revealed significantly stronger connectivity between subcortical (SC), auditory (AUD) and visual (VIS) networks in patients, as well as hypoconnectivity in sensorimotor (SM) network relative to controls. sFNG analysis highlighted distinctive clustering patterns in patients and HCs along cognitive control (CC)/ default mode network (DMN), as well as SC/ AUD/ SM/ cerebellar (CB), and VIS gradients. Furthermore, we observed significant differences in the sFNGs between groups in SC and CB domains. dFNG analysis suggested that SZ patients spend significantly more time in a SC/ CB state based on the first gradient, while HCs favor the SM/DMN state. For the second gradient, however, patients exhibited significantly higher activity in CB domains, contrasting with HCs' DMN engagement. The gradient synchrony analysis conveyed more shifts between SM/ SC networks and transmodal CC/ DMN networks in patients. In addition, the dFNG coupling revealed distinct connectivity patterns between SC, SM and CB domains in SZ patients compared to HCs. To recap, our results advance our understanding of brain network modulation by examining smooth connectivity trajectories. This provides a more complete spatiotemporal summary of the data, contributing to the growing body of current literature regarding the functional dysconnectivity in schizophrenia patients. By employing dFNG, we highlight a new perspective to capture large scale fluctuations across the brain while maintaining the convenience of brain networks and low dimensional summary measures.

More reliable biomarkers and more accurate prediction for mental disorders using a label-noise filtering-based dimensional prediction method.

The integration of neuroimaging with artificial intelligence is crucial for advancing the diagnosis of mental disorders. However, challenges arise from incomplete matching between diagnostic labels and neuroimaging. Here, we propose a label-noise filtering-based dimensional prediction (LAMP) method to identify reliable biomarkers and achieve accurate prediction for mental disorders. Our method proposes to utilize a label-noise filtering model to automatically filter out unclear cases from a neuroimaging perspective, and then the typical subjects whose diagnostic labels align with neuroimaging measures are used to construct a dimensional prediction model to score independent subjects. Using fMRI data of schizophrenia patients and healthy controls (n = 1,245), our method yields consistent scores to independent subjects, leading to more distinguishable relabeled groups with an enhanced classification accuracy of 31.89%. Additionally, it enables the exploration of stable abnormalities in schizophrenia. In summary, our LAMP method facilitates the identification of reliable biomarkers and accurate diagnosis of mental disorders using neuroimages.

Linking enlarged choroid plexus with plasma analyte and structural phenotypes in clinical high risk for psychosis: A multisite neuroimaging study.

BACKGROUND: Choroid plexus (ChP) enlargement exists in first-episode and chronic psychosis, but whether enlargement occurs before psychosis onset is unknown. This study investigated whether ChP volume is enlarged in individuals with clinical high-risk (CHR) for psychosis and whether these changes are related to clinical, neuroanatomical, and plasma analytes. METHODS: Clinical and neuroimaging data from the North American Prodrome Longitudinal Study 2 (NAPLS2) was used for analysis. 509 participants (169 controls, 340 CHR) were recruited. Conversion status was determined after 2-years of follow-up, with 36 psychosis converters. The lateral ventricle ChP was manually segmented from baseline scans. A subsample of 31 controls and 53 CHR had plasma analyte and neuroimaging data. RESULTS: Compared to controls, CHR (d = 0.23, p = 0.017) and non-converters (d = 0.22, p = 0.03) demonstrated higher ChP volumes, but not in converters. In CHR, greater ChP volume correlated with lower cortical (r = -0.22, p < 0.001), subcortical gray matter (r = -0.21, p < 0.001), and total white matter volume (r = -0.28,p < 0.001), as well as larger lateral ventricle volume (r = 0.63,p < 0.001). Greater ChP volume correlated with makers functionally associated with the lateral ventricle ChP in CHR [CCL1 (r = -0.30, p = 0.035), ICAM1 (r = 0.33, p = 0.02)], converters [IL1ß (r = 0.66, p = 0.004)], and non-converters [BMP6 (r = -0.96, p < 0.001), CALB1 (r = -0.98, p < 0.001), ICAM1 (r = 0.80, p = 0.003), SELE (r = 0.59, p = 0.026), SHBG (r = 0.99, p < 0.001), TNFRSF10C (r = 0.78, p = 0.001)]. CONCLUSIONS: CHR and non-converters demonstrated significantly larger ChP volumes compared to controls. Enlarged ChP was associated with neuroanatomical alterations and analyte markers functionally associated with the ChP. These findings suggest that the ChP may be a key an important biomarker in CHR.

Spatial Dynamic Subspaces Encode Sex-Specific Schizophrenia Disruptions in Transient Network Overlap and Their Links to Genetic Risk.

BACKGROUND: Schizophrenia research reveals sex differences in incidence, symptoms, genetic risk factors, and brain function. However, a knowledge gap remains regarding sex-specific schizophrenia alterations in brain function. Schizophrenia is considered a dysconnectivity syndrome, but the dynamic integration and segregation of brain networks are poorly understood. Recent advances in resting-state functional magnetic resonance imaging allow us to study spatial dynamics, the phenomenon of brain networks spatially evolving over time. Nevertheless, estimating time-resolved networks remains challenging due to low signal-to-noise ratio, limited short-time information, and uncertain network identification. METHODS: We adapted a reference-informed network estimation technique to capture time-resolved networks and their dynamic spatial integration and segregation for 193 individuals with schizophrenia and 315 control participants. We focused on time-resolved spatial functional network connectivity, an estimate of network spatial coupling, to study sex-specific alterations in schizophrenia and their links to genomic data. RESULTS: Our findings are consistent with the dysconnectivity and neurodevelopment hypotheses and with the cerebello-thalamo-cortical, triple-network, and frontoparietal dysconnectivity models, helping to unify them. The potential unification offers a new understanding of the underlying mechanisms. Notably, the posterior default mode/salience spatial functional network connectivity exhibits sex-specific schizophrenia alteration during the state with the highest global network integration and is correlated with genetic risk for schizophrenia. This dysfunction is reflected in regions with weak functional connectivity to corresponding networks. CONCLUSIONS: Our method can effectively capture spatially dynamic networks, detect nuanced schizophrenia effects including sex-specific ones, and reveal the intricate relationship of dynamic information to genomic data. The results also underscore the clinical potential of dynamic spatial dependence and weak connectivity.

Networks extracted from nonlinear fMRI connectivity exhibit unique spatial variation and enhanced sensitivity to differences between individuals with schizophrenia and controls.

Functional magnetic resonance imaging (fMRI) studies often estimate brain intrinsic connectivity networks (ICNs) from temporal relationships between hemodynamic signals using approaches such as independent component analysis (ICA). While ICNs are thought to represent functional sources that play important roles in various psychological phenomena, current approaches have been tailored to identify ICNs that mainly reflect linear statistical relationships. However, the elements comprising neural systems often exhibit remarkably complex nonlinear interactions that may be involved in cognitive operations and altered in psychiatric conditions such as schizophrenia. Consequently, there is a need to develop methods capable of effectively capturing ICNs from measures that are sensitive to nonlinear relationships. Here, we advance a novel approach to estimate ICNs from explicitly nonlinear whole-brain functional connectivity (ENL-wFC) by transforming resting-state fMRI (rsfMRI) data into the connectivity domain, allowing us to capture unique information from distance correlation patterns that would be missed by linear whole-brain functional connectivity (LIN-wFC) analysis. Our findings provide evidence that ICNs commonly extracted from linear (LIN) relationships are also reflected in explicitly nonlinear (ENL) connectivity patterns. ENL ICN estimates exhibit higher reliability and stability, highlighting our approach's ability to effectively quantify ICNs from rsfMRI data. Additionally, we observed a consistent spatial gradient pattern between LIN and ENL ICNs with higher ENL weight in core ICN regions, suggesting that ICN function may be subserved by nonlinear processes concentrated within network centers. We also found that a uniquely identified ENL ICN distinguished individuals with schizophrenia from healthy controls while a uniquely identified LIN ICN did not, emphasizing the valuable complementary information that can be gained by incorporating measures that are sensitive to nonlinearity in future analyses. Moreover, the ENL estimates of ICNs associated with auditory, linguistic, sensorimotor, and self-referential processes exhibit heightened sensitivity towards differentiating between individuals with schizophrenia and controls compared to LIN counterparts, demonstrating the translational value of our approach and of the ENL estimates of ICNs that are frequently reported as disrupted in schizophrenia. In summary, our findings underscore the tremendous potential of connectivity domain ICA and nonlinear information in resolving complex brain phenomena and revolutionizing the landscape of clinical FC analysis.

4D DYNAMIC SPATIAL BRAIN NETWORKS AT REST LINKED TO COGNITION SHOW ATYPICAL VARIABILITY AND COUPLING IN SCHIZOPHRENIA.

Despite increasing interest in the dynamics of functional brain networks, most studies focus on the changing relationships over time between spatially static networks or regions. Here we propose an approach to study dynamic spatial brain net-works in human resting state functional magnetic resonance imaging (rsfMRI) data and evaluate the temporal changes in the volumes of these 4D networks. Our results show significant volumetric coupling (i.e., synchronized shrinkage and growth) between networks during the scan. We find that several features of such dynamic spatial brain networks are associated with cognition, with higher dynamic variability in these networks and higher volumetric coupling between network pairs positively associated with cognitive performance. We show that these networks are modulated differently in individuals with schizophrenia versus typical controls, resulting in network growth or shrinkage, as well as altered focus of activity within a network. Schizophrenia also shows lower spatial dynamical variability in several networks, and lower volumetric coupling between pairs of networks, thus upholding the role of dynamic spatial brain networks in cognitive impairment seen in schizophrenia. Our data show evidence for the importance of studying the typically overlooked voxelwise changes within and between brain networks.

Neural Correlates of Positive and Negative Formal Thought Disorder in Individuals with Schizophrenia: An ENIGMA Schizophrenia Working Group Study.

Formal thought disorder (FTD) is a key clinical factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, relationship between FTD symptom dimensions and patterns of regional brain volume deficiencies in schizophrenia remain to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles based on a large multi-site cohort through the ENIGMA Schizophrenia Working Group (752 individuals with schizophrenia and 1256 controls), to unravel the neuroanatomy of positive, negative and total FTD in schizophrenia and their cellular bases. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks for positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but 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 was also linked to microglial cell types. These findings relate different dimensions of FTD to distinct brain structural changes and their cellular underpinnings, improve our mechanistic understanding of these key psychotic symptoms.

Auditory oddball hypoactivation in schizophrenia.

Individuals with schizophrenia (SZ) show aberrant activations, assessed via functional magnetic resonance imaging (fMRI), during auditory oddball tasks. However, associations with cognitive performance and genetic contributions remain unknown. This study compares individuals with SZ to healthy volunteers (HVs) using two cross-sectional data sets from multi-center brain imaging studies. It examines brain activation to auditory oddball targets, and their associations with cognitive domain performance, schizophrenia polygenic risk scores (PRS), and genetic variation (loci). Both sample 1 (137 SZ vs. 147 HV) and sample 2 (91 SZ vs. 98 HV), showed hypoactivation in SZ in the left-frontal pole, and right frontal orbital, frontal pole, paracingulate, intracalcarine, precuneus, supramarginal and hippocampal cortices, and right thalamus. In SZ, precuneus activity was positively related to cognitive performance. Schizophrenia PRS showed a negative correlation with brain activity in the right-supramarginal cortex. GWA analyses revealed significant single-nucleotide polymorphisms associated with right-supramarginal gyrus activity. RPL36 also predicted right-supramarginal gyrus activity. In addition to replicating hypoactivation for oddball targets in SZ, this study identifies novel relationships between regional activity, cognitive performance, and genetic loci that warrant replication, emphasizing the need for continued data sharing and collaborative efforts.

Subcortical volumetric alterations in four major psychiatric disorders: a mega-analysis study of 5604 subjects and a volumetric data-driven approach for classification.

Differential diagnosis is sometimes difficult in practical psychiatric settings, in terms of using the current diagnostic system based on presenting symptoms and signs. The creation of a novel diagnostic system using objective biomarkers is expected to take place. Neuroimaging studies and others reported that subcortical brain structures are the hubs for various psycho-behavioral functions, while there are so far no neuroimaging data-driven clinical criteria overcoming limitations of the current diagnostic system, which would reflect cognitive/social functioning. Prior to the main analysis, we conducted a large-scale multisite study of subcortical volumetric and lateralization alterations in schizophrenia, bipolar disorder, major depressive disorder, and autism spectrum disorder using T1-weighted images of 5604 subjects (3078 controls and 2526 patients). We demonstrated larger lateral ventricles volume in schizophrenia, bipolar disorder, and major depressive disorder, smaller hippocampus volume in schizophrenia and bipolar disorder, and schizophrenia-specific smaller amygdala, thalamus, and accumbens volumes and larger caudate, putamen, and pallidum volumes. In addition, we observed a leftward alteration of lateralization for pallidum volume specifically in schizophrenia. Moreover, as our main objective, we clustered the 5,604 subjects based on subcortical volumes, and explored whether data-driven clustering results can explain cognitive/social functioning in the subcohorts. We showed a four-biotype classification, namely extremely (Brain Biotype [BB] 1) and moderately smaller limbic regions (BB2), larger basal ganglia (BB3), and normal volumes (BB4), being associated with cognitive/social functioning. Specifically, BB1 and BB2-3 were associated with severe and mild cognitive/social impairment, respectively, while BB4 was characterized by normal cognitive/social functioning. Our results may lead to the future creation of novel biological data-driven psychiatric diagnostic criteria, which may be expected to be useful for prediction or therapeutic selection.

NIDM-Terms: community-based terminology management for improved neuroimaging dataset descriptions and query.

The biomedical research community is motivated to share and reuse data from studies and projects by funding agencies and publishers. Effectively combining and reusing neuroimaging data from publicly available datasets, requires the capability to query across datasets in order to identify cohorts that match both neuroimaging and clinical/behavioral data criteria. Critical barriers to operationalizing such queries include, in part, the broad use of undefined study variables with limited or no annotations that make it difficult to understand the data available without significant interaction with the original authors. Using the Brain Imaging Data Structure (BIDS) to organize neuroimaging data has made querying across studies for specific image types possible at scale. However, in BIDS, beyond file naming and tightly controlled imaging directory structures, there are very few constraints on ancillary variable naming/meaning or experiment-specific metadata. In this work, we present NIDM-Terms, a set of user-friendly terminology management tools and associated software to better manage individual lab terminologies and help with annotating BIDS datasets. Using these tools to annotate BIDS data with a Neuroimaging Data Model (NIDM) semantic web representation, enables queries across datasets to identify cohorts with specific neuroimaging and clinical/behavioral measurements. This manuscript describes the overall informatics structures and demonstrates the use of tools to annotate BIDS datasets to perform integrated cross-cohort queries.

Cerebral cortical structural alteration patterns across four major psychiatric disorders in 5549 individuals.

According to the operational diagnostic criteria, psychiatric disorders such as schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and autism spectrum disorder (ASD) are classified based on symptoms. While its cluster of symptoms defines each of these psychiatric disorders, there is also an overlap in symptoms between the disorders. We hypothesized that there are also similarities and differences in cortical structural neuroimaging features among these psychiatric disorders. T1-weighted magnetic resonance imaging scans were performed for 5,549 subjects recruited from 14 sites. Effect sizes were determined using a linear regression model within each protocol, and these effect sizes were meta-analyzed. The similarity of the differences in cortical thickness and surface area of each disorder group was calculated using cosine similarity, which was calculated from the effect sizes of each cortical regions. The thinnest cortex was found in SZ, followed by BD and MDD. The cosine similarity values between disorders were 0.943 for SZ and BD, 0.959 for SZ and MDD, and 0.943 for BD and MDD, which indicated that a common pattern of cortical thickness alterations was found among SZ, BD, and MDD. Additionally, a generally smaller cortical surface area was found in SZ and MDD than in BD, and the effect was larger in SZ. The cosine similarity values between disorders were 0.945 for SZ and MDD, 0.867 for SZ and ASD, and 0.811 for MDD and ASD, which indicated a common pattern of cortical surface area alterations among SZ, MDD, and ASD. Patterns of alterations in cortical thickness and surface area were revealed in the four major psychiatric disorders. To our knowledge, this is the first report of a cross-disorder analysis conducted on four major psychiatric disorders. Cross-disorder brain imaging research can help to advance our understanding of the pathogenesis of psychiatric disorders and common symptoms.

Cortical morphology in patients with the deficit and non-deficit syndrome of schizophrenia: a worldwide meta- and mega-analyses.

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.

A deep residual model for characterization of 5D spatiotemporal network dynamics reveals widespread spatiodynamic changes in schizophrenia.

Schizophrenia is a severe brain disorder with serious symptoms including delusions, disorganized speech, and hallucinations that can have a long-term detrimental impact on different aspects of a patient's life. It is still unclear what the main cause of schizophrenia is, but a combination of altered brain connectivity and structure may play a role. Neuroimaging data has been useful in characterizing schizophrenia, but there has been very little work focused on voxel-wise changes in multiple brain networks over time, despite evidence that functional networks exhibit complex spatiotemporal changes over time within individual subjects. Recent studies have primarily focused on static (average) features of functional data or on temporal variations between fixed networks; however, such approaches are not able to capture multiple overlapping networks which change at the voxel level. In this work, we employ a deep residual convolutional neural network (CNN) model to extract 53 different spatiotemporal networks each of which captures dynamism within various domains including subcortical, cerebellar, visual, sensori-motor, auditory, cognitive control, and default mode. We apply this approach to study spatiotemporal brain dynamism at the voxel level within multiple functional networks extracted from a large functional magnetic resonance imaging (fMRI) dataset of individuals with schizophrenia (N = 708) and controls (N = 510). Our analysis reveals widespread group level differences across multiple networks and spatiotemporal features including voxel-wise variability, magnitude, and temporal functional network connectivity in widespread regions expected to be impacted by the disorder. We compare with typical average spatial amplitude and show highly structured and neuroanatomically relevant results are missed if one does not consider the voxel-wise spatial dynamics. Importantly, our approach can summarize static, temporal dynamic, spatial dynamic, and spatiotemporal dynamics features, thus proving a powerful approach to unify and compare these various perspectives. In sum, we show the proposed approach highlights the importance of accounting for both temporal and spatial dynamism in whole brain neuroimaging data generally, shows a high-level of sensitivity to schizophrenia highlighting global but spatially unique dynamics showing group differences, and may be especially important in studies focused on the development of brain-based biomarkers.

Alterations in grey matter structure linked to frequency-specific cortico-subcortical connectivity in schizophrenia via multimodal data fusion.

Schizophrenia (SZ) is a complex psychiatric disorder that is currently defined by symptomatic and behavioral, rather than biological, criteria. Neuroimaging is an appealing avenue for SZ biomarker development, as several neuroimaging-based studies comparing individuals with SZ to healthy controls (HC) have shown measurable group differences in brain structure, as well as functional brain alterations in both static and dynamic functional network connectivity (sFNC and dFNC, respectively). The recently proposed filter-banked connectivity (FBC) method extends the standard dFNC sliding-window approach to estimate FNC within an arbitrary number of distinct frequency bands. The initial implementation used a set of filters spanning the full connectivity spectral range, providing a unified approach to examine both sFNC and dFNC in a single analysis. Initial FBC results found that individuals with SZ spend more time in a less structured, more disconnected low-frequency (i.e., static) FNC state than HC, as well as preferential SZ occupancy in high-frequency connectivity states, suggesting a frequency-specific component underpinning the functional dysconnectivity observed in SZ. Building on these findings, we sought to link such frequency-specific patterns of FNC to covarying data-driven structural brain networks in the context of SZ. Specifically, we employ a multi-set canonical correlation analysis + joint independent components analysis (mCCA + jICA) data fusion framework to study the connection between grey matter volume (GMV) maps and FBC states across the full connectivity frequency spectrum. Our multimodal analysis identified two joint sources that captured co-varying patterns of frequency-specific functional connectivity and alterations in GMV with significant group differences in loading parameters between the SZ group and HC. The first joint source linked frequency-modulated connections between the subcortical and sensorimotor networks and GMV alterations in the frontal and temporal lobes, while the second joint source identified a relationship between low-frequency cerebellar-sensorimotor connectivity and structural changes in both the cerebellum and motor cortex. Together, these results show a strong connection between cortico-subcortical functional connectivity at both high and low frequencies and alterations in cortical GMV that may be relevant to the pathogenesis and pathophysiology of SZ.

Neural Correlates of Positive and Negative Formal Thought Disorder in Individuals with Schizophrenia: An ENIGMA Schizophrenia Working Group Study.

Formal thought disorder (FTD) is a key clinical factor in schizophrenia, but the neurobiological underpinnings remain unclear. In particular, relationship between FTD symptom dimensions and patterns of regional brain volume deficiencies in schizophrenia remain to be established in large cohorts. Even less is known about the cellular basis of FTD. Our study addresses these major obstacles based on a large multi-site cohort through the ENIGMA Schizophrenia Working Group (752 individuals with schizophrenia and 1256 controls), to unravel the neuroanatomy of positive, negative and total FTD in schizophrenia and their cellular bases. We used virtual histology tools to relate brain structural changes associated with FTD to cellular distributions in cortical regions. We identified distinct neural networks for positive and negative FTD. Both networks encompassed fronto-occipito-amygdalar brain regions, but 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 was also linked to microglial cell types. These findings relate different dimensions of FTD to distinct brain structural changes and their cellular underpinnings, improve our mechanistic understanding of these key psychotic symptoms.

Multi-model order spatially constrained ICA reveals highly replicable group differences and consistent predictive results from resting data: A large N fMRI schizophrenia study.

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.