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.

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.

Objective

Multimodal measurements of the same phenomena provide complementary information and highlight different perspectives, albeit each with their own limitations. A focus on a single modality may lead to incorrect inferences, which is especially important when a studied phenomenon is a disease. In this paper, we introduce a method that takes advantage of multimodal data in addressing the hypotheses of disconnectivity and dysfunction within schizophrenia (SZ).

Methods

We start with estimating and visualizing links within and among extracted multimodal data features using a Gaussian graphical model (GGM). We then propose a modularity-based method that can be applied to the GGM to identify links that are associated with mental illness across a multimodal data set. Through simulation and real data, we show our approach reveals important information about disease-related network disruptions that are missed with a focus on a single modality. We use functional MRI (fMRI), diffusion MRI (dMRI), and structural MRI (sMRI) to compute the fractional amplitude of low frequency fluctuations (fALFF), fractional anisotropy (FA), and gray matter (GM) concentration maps. These three modalities are analyzed using our modularity method.

Results

Our results show missing links that are only captured by the cross-modal information that may play an important role in disconnectivity between the components.

Conclusion

We identified multimodal (fALFF, FA and GM) disconnectivity in the default mode network area in patients with SZ, which would not have been detectable in a single modality.

Significance

The proposed approach provides an important new tool for capturing information that is distributed among multiple imaging modalities.

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.

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

Background

A 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.

Methods

We 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.

Results

We 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.

Conclusion

Our 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.

Objective

Previous 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.

Methods

Volumetric 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 I_ATT), and 60 typically developing age-matched controls.

Results

The 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 I_ATT.

Conclusions

We 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 I_ATT group, suggesting that disease effects contribute substantially to the clinical attentional difficulties seen in this population.

Background

Comprehensive 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.

Methods

The 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.

Discussion

This 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 registration

ClinicalTrials.gov NCT04578912 . Registered on October 8, 2020.

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.

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.