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.

The direct effect of genetic variations on clinical phenotypes within schizophrenia (SZ) remains elusive. We examined the previously identified association of reduced gray matter concentration in the insula - medial prefrontal cortex and a quantitative trait locus located in 12q24 in a SZ dataset. The main analysis was performed on 1461 SNPs and 830 participants. The highest contributing SNPs were localized in five genes including TMEM119, which encodes a microglial marker, that is associated with neuroinflammation and Alzheimer's disease. The gene set in 12q4 may partially explain brain alterations in SZ, but they may also relate to other psychiatric and developmental disorders.

Objective

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

Methods

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

Results

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

Conclusion

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

Significance

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

Background

Structural neuroimaging studies have identified similarities in the brains of individuals diagnosed with schizophrenia (SZ) and bipolar I disorder (BP), with overlap in regions of gray matter (GM) deficits between the two disorders. Recent studies have also shown that the symptom phenotypes associated with SZ and BP may allow for a more precise categorization than the current diagnostic criteria. In this study, we sought to identify GM alterations that were unique to each disorder and whether those alterations were also related to unique symptom profiles.

Materials and methods

We analyzed the GM patterns and clinical symptom presentations using independent component analysis (ICA), hierarchical clustering, and n-way biclustering in a large (N ∼ 3,000), merged dataset of neuroimaging data from healthy volunteers (HV), and individuals with either SZ or BP.

Results

Component A showed a SZ and BP < HV GM pattern in the bilateral insula and cingulate gyrus. Component B showed a SZ and BP < HV GM pattern in the cerebellum and vermis. There were no significant differences between diagnostic groups in these components. Component C showed a SZ < HV and BP GM pattern bilaterally in the temporal poles. Hierarchical clustering of the PANSS scores and the ICA components did not yield new subgroups. N-way biclustering identified three unique subgroups of individuals within the sample that mapped onto different combinations of ICA components and symptom profiles categorized by the PANSS but no distinct diagnostic group differences.

Conclusion

These multivariate results show that diagnostic boundaries are not clearly related to structural differences or distinct symptom profiles. Our findings add support that (1) BP tend to have less severe symptom profiles when compared to SZ on the PANSS without a clear distinction, and (2) all the gray matter alterations follow the pattern of SZ < BP < HV without a clear distinction between SZ and BP.

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

Background

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

Methods

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

Results

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

Conclusions

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