Search

BackgroundComprehensive Behavioral Intervention for Tics (CBIT) is a first-line treatment for tic disorders that aims to improve controllability over tics that an individual finds distressing or impairing. However, it is only effective for approximately half of patients. Supplementary motor area (SMA)-directed neurocircuitry plays a strong role in motor inhibition, and activity in this region is thought to contribute to tic expression. Targeted modulation of SMA using transcranial magnetic stimulation (TMS) may increase CBIT efficacy by improving patients' ability to implement tic controllability behaviors.MethodsThe CBIT + TMS trial is a two-phase, milestone-driven early-stage randomized controlled trial. The trial will test whether augmenting CBIT with inhibitory, non-invasive stimulation of SMA with TMS modifies activity in SMA-mediated circuits and enhances tic controllability in youth ages 12-21 years with chronic tics. Phase 1 will directly compare two rTMS augmentation strategies (1 Hz rTMS vs. cTBS) vs. sham in N = 60 participants. Quantifiable, a priori "Go/No Go Criteria" guide the decision to proceed to phase 2 and the selection of the optimal TMS regimen. Phase 2 will compare the optimal regimen vs. sham and test the link between neural target engagement and clinical outcomes in a new sample of N = 60 participants.DiscussionThis clinical trial is one of few to date testing TMS augmentation of therapy in a pediatric sample. The results will provide insight into whether TMS is a potentially viable strategy for enhancing CBIT efficacy and reveal potential neural and behavioral mechanisms of change.Trial registrationClinicalTrials.gov NCT04578912 . Registered on October 8, 2020.

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.

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.

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.

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.

ObjectiveTo assess our hypothesis that brain macrostructure is different in individuals with mucopolysaccharidosis type I (MPS I) and healthy controls (HC), we conducted a comprehensive multicenter study using a uniform quantitative magnetic resonance imaging (qMRI) protocol, with analyses that account for the effects of disease phenotype, age, and cognition.MethodsBrain MRIs in 23 individuals with attenuated (MPS IA) and 38 with severe MPS I (MPS IH), aged 4-25 years, enrolled under the study protocol NCT01870375, were compared to 98 healthy controls.ResultsCortical and subcortical gray matter, white matter, corpus callosum, ventricular and choroid plexus volumes in MPS I significantly differed from HC. Thicker cortex, lower white matter and corpus callosum volumes were already present at the youngest MPS I participants aged 4-5 years. Age-related differences were observed in both MPS I groups, but most markedly in MPS IH, particularly in cortical gray matter metrics. IQ scores were inversely associated with ventricular volume in both MPS I groups and were positively associated with cortical thickness only in MPS IA.ConclusionsQuantitatively-derived MRI measures distinguished MPS I participants from HC as well as severe from attenuated forms. Age-related neurodevelopmental trajectories in both MPS I forms differed from HC. The extent to which brain structure is altered by disease, potentially spared by treatment, and how it relates to neurocognitive dysfunction needs further exploration.

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.

Major depressive disorder (MDD) is associated with an increased risk of brain atrophy, aging-related diseases, and mortality. We examined potential advanced brain aging in adult MDD patients, and whether this process is associated with clinical characteristics in a large multicenter international dataset. We performed a mega-analysis by pooling brain measures derived from T1-weighted MRI scans from 19 samples worldwide. Healthy brain aging was estimated by predicting chronological age (18-75 years) from 7 subcortical volumes, 34 cortical thickness and 34 surface area, lateral ventricles and total intracranial volume measures separately in 952 male and 1236 female controls from the ENIGMA MDD working group. The learned model coefficients were applied to 927 male controls and 986 depressed males, and 1199 female controls and 1689 depressed females to obtain independent unbiased brain-based age predictions. The difference between predicted "brain age" and chronological age was calculated to indicate brain-predicted age difference (brain-PAD). On average, MDD patients showed a higher brain-PAD of +1.08 (SE 0.22) years (Cohen's d = 0.14, 95% CI: 0.08-0.20) compared with controls. However, this difference did not seem to be driven by specific clinical characteristics (recurrent status, remission status, antidepressant medication use, age of onset, or symptom severity). This highly powered collaborative effort showed subtle patterns of age-related structural brain abnormalities in MDD. Substantial within-group variance and overlap between groups were observed. Longitudinal studies of MDD and somatic health outcomes are needed to further assess the clinical value of these brain-PAD estimates.

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.

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.

BackgroundSchizophrenia is associated with robust hippocampal volume deficits but subregion volume deficits, their associations with cognition, and contributing genes remain to be determined.MethodsHippocampal formation (HF) subregion volumes were obtained using FreeSurfer 6.0 from individuals with schizophrenia (n = 176, mean age ± s.d. = 39.0 ± 11.5, 132 males) and healthy volunteers (n = 173, mean age ± s.d. = 37.6 ± 11.3, 123 males) with similar mean age, gender, handedness, and race distributions. Relationships between the HF subregion volume with the largest between group difference, neuropsychological performance, and single-nucleotide polymorphisms were assessed.ResultsThis study found a significant group by region interaction on hippocampal subregion volumes. Compared to healthy volunteers, individuals with schizophrenia had significantly smaller dentate gyrus (DG) (Cohen's d = -0.57), Cornu Ammonis (CA) 4, molecular layer of the hippocampus, hippocampal tail, and CA 1 volumes, when statistically controlling for intracranial volume; DG (d = -0.43) and CA 4 volumes remained significantly smaller when statistically controlling for mean hippocampal volume. DG volume showed the largest between group difference and significant positive associations with visual memory and speed of processing in the overall sample. Genome-wide association analysis with DG volume as the quantitative phenotype identified rs56055643 (β = 10.8, p < 5 × 10-8, 95% CI 7.0-14.5) on chromosome 3 in high linkage disequilibrium with MOBP. Gene-based analyses identified associations between SLC25A38 and RPSA and DG volume.ConclusionsThis study suggests that DG dysfunction is fundamentally involved in schizophrenia pathophysiology, that it may contribute to cognitive abnormalities in schizophrenia, and that underlying biological mechanisms may involve contributions from MOBP, SLC25A38, and RPSA.

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.

ObjectiveWe propose and develop a novel biclustering (N-BiC) approach for performing N-way biclustering of neuroimaging data. Our approach is applicable to an arbitrary number of features from both imaging and behavioral data (e.g., symptoms). We applied it to structural MRI data from patients with schizophrenia.MethodsIt uses a source-based morphometry approach [i.e., independent component analysis of gray matter segmentation maps] to decompose the data into a set of spatial maps, each of which includes regions that covary among individuals. Then, the loading parameters for components of interest are entered to an exhaustive search, which incorporates a modified depth-first search technique to carry out the biclustering, with the goal of obtaining submatrices where the selected rows (individuals) show homogeneity in their expressions of selected columns (components) and vice versa.ResultsFindings demonstrate that multiple biclusters have an evident association with distinct brain networks for the different types of symptoms in schizophrenia. The study identifies two components: inferior temporal gyrus (16) and brainstem (7), which are related to positive (distortion/excess of normal function) and negative (diminution/loss of normal function) symptoms in schizophrenia, respectively.ConclusionN-BiC is a data-driven method of biclustering MRI data that can exhaustively explore relationships/substructures from a dataset without any prior information with a higher degree of robustness than earlier biclustering applications.SignificanceThe use of such approaches is important to investigate the underlying biological substrates of mental illness by grouping patients into homogeneous subjects, as the schizophrenia diagnosis is known to be relatively nonspecific and heterogeneous.

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.

Schizophrenia is a complex disorder that affects perception, cognition, and emotion causing symptoms such as delusions, hallucinations, and suspiciousness. Schizophrenia is also associated with structural cortical abnormalities including lower gray matter (GM) concentration, GM volume, and cortical thickness relative to healthy control individuals. However, the association between GM measures and symptom dimensions in schizophrenia is still not well understood. Here, we applied parallel independent component analysis (pICA), a higher-order statistical approach that identifies covarying patterns within two (or more) data modalities simultaneously, to link covarying brain networks of GM concentration with covarying linear combinations of the positive and negative syndrome scale (PANSS). In a large sample of patients with schizophrenia (n = 337) the association between these two data modalities was investigated. The pICA revealed a distinct PANSS profile characterized by increased delusional symptoms, suspiciousness, hallucinations, and anxiety, that was associated with a pattern of lower GM concentration in inferior temporal gyri and fusiform gyri and higher GM concentration in the sensorimotor cortex. GM alterations replicate previous findings; additionally, applying a multivariate technique, we were able to map a very specific symptom profile onto these GM alterations extending our understanding of cortical abnormalities associated with schizophrenia. Techniques like parallel ICA can reveal linked patterns of alterations across different data modalities that can help to identify biologically-informed phenotypes which might help to improve future treatment targets.

ObjectivePrevious research suggests attention and white matter (WM) abnormalities in individuals with mucopolysaccharidosis type I (MPS I); this cross-sectional comparison is one of the first to examine the relationship of WM structural abnormalities as measured by corpus callosum (CC) volumes with attention scores to evaluate this relationship in a larger sample of patients with MPS I.MethodsVolumetric MRI data and performance on a computerized measure of sustained attention were compared for 18 participants with the severe form of MPS I (MPS IH), 18 participants with the attenuated form of MPS I (MPS IATT), and 60 typically developing age-matched controls.ResultsThe MPS I groups showed below-average mean attention scores (p < 0.001) and smaller CC volumes (p < 0.001) than controls. No significant associations were found between attention performance and CC volume for controls. Attention was associated with posterior CC volumes in the participants with MPS IH (p = 0.053) and total (p = 0.007) and anterior (p < 0.001) CC volumes in participants with MPS IATT.ConclusionsWe found that attention and CC volumes were reduced in participants with MPS I compared to typically developing controls. Smaller CC volumes in participants with MPS I were associated with decreased attention; such an association was not seen in controls. While hematopoietic cell transplantation used to treat MPS IH may compound these effects, attention difficulties were also seen in the MPS IATT group, suggesting that disease effects contribute substantially to the clinical attentional difficulties seen in this population.

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.

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.

This study assessed genetic contributions to six cognitive domains, identified by the MATRICS Cognitive Consensus Battery as relevant for schizophrenia, cognition-enhancing, clinical trials. Psychiatric Genomics Consortium Schizophrenia polygenic risk scores showed significant negative correlations with each cognitive domain. Genome-wide association analyses identified loci associated with attention/vigilance (rs830786 within HNF4G), verbal memory (rs67017972 near NDUFS4), and reasoning/problem solving (rs76872642 within HDAC9). Gene set analysis identified unique and shared genes across cognitive domains. These findings suggest involvement of common and unique mechanisms across cognitive domains and may contribute to the discovery of new therapeutic targets to treat cognitive deficits in schizophrenia.

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.

Hallucinations characterize schizophrenia, with approximately 59% of patients reporting auditory hallucinations and 27% reporting visual hallucinations. Prior neuroimaging studies suggest that hallucinations are linked to disrupted communication across distributed (sensory, salience-monitoring and subcortical) networks. Yet, our understanding of the neurophysiological mechanisms that underlie auditory and visual hallucinations in schizophrenia remains limited. This study integrates two resting-state functional magnetic resonance imaging (fMRI) analysis methods - amplitudes of low-frequency fluctuations (ALFF) and functional network connectivity (FNC) - to explore the hypotheses that (1) abnormal FNC between salience and sensory (visual/auditory) networks underlies hallucinations in schizophrenia, and (2) disrupted hippocampal oscillations (as measured by hippocampal ALFF) beget changes in FNC linked to hallucinations. Our first hypothesis was supported by the finding that schizophrenia patients reporting hallucinations have higher FNC between the salience network and an associative auditory network relative to healthy controls. Hippocampal ALFF was negatively associated with FNC between primary auditory cortex and the salience network in healthy subjects, but was positively associated with FNC between these networks in patients reporting hallucinations. These findings provide indirect support favoring our second hypothesis. We suggest future studies integrate fMRI with electroencephalogram (EEG) and/or magnetoencephalogram (MEG) methods to directly probe the temporal relation between altered hippocampal oscillations and changes in cross-network functional communication.

Cognitive impairment is a feature of many psychiatric diseases, including schizophrenia. Here we aim to identify multimodal biomarkers for quantifying and predicting cognitive performance in individuals with schizophrenia and healthy controls. A supervised learning strategy is used to guide three-way multimodal magnetic resonance imaging (MRI) fusion in two independent cohorts including both healthy individuals and individuals with schizophrenia using multiple cognitive domain scores. Results highlight the salience network (gray matter, GM), corpus callosum (fractional anisotropy, FA), central executive and default-mode networks (fractional amplitude of low-frequency fluctuation, fALFF) as modality-specific biomarkers of generalized cognition. FALFF features are found to be more sensitive to cognitive domain differences, while the salience network in GM and corpus callosum in FA are highly consistent and predictive of multiple cognitive domains. These modality-specific brain regions define-in three separate cohorts-promising co-varying multimodal signatures that can be used as predictors of multi-domain cognition.