Search

Although specific risk factors for brain alterations in bipolar disorders (BD) are currently unknown, obesity impacts the brain and is highly prevalent in BD. Gray matter correlates of obesity in BD have been well documented, but we know much less about brain white matter abnormalities in people who have both obesity and BD. We obtained body mass index (BMI) and diffusion tensor imaging derived fractional anisotropy (FA) from 22 white matter tracts in 899 individuals with BD, and 1287 control individuals from 20 cohorts in the ENIGMA-BD working group. In a mega-analysis, we investigated the associations between BMI, diagnosis or medication and FA. Lower FA was associated with both BD and BMI in six white matter tracts, including the corpus callosum and thalamic radiation. Higher BMI or BD were uniquely associated with lower FA in three and six white matter tracts, respectively. People not receiving lithium treatment had a greater negative association between FA and BMI than people treated with lithium in the posterior thalamic radiation and sagittal stratum. In three tracts BMI accounted for 10.5 to 17% of the negative association between the number of medication classes other than lithium and FA. Both overweight/obesity and BD demonstrated lower FA in some of the same regions. People prescribed lithium had a weaker association between BMI and FA than people not on lithium. In contrast, greater weight contributed to the negative associations between medications and FA. Obesity may add to brain alterations in BD and may play a role in effects of medications on the brain., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Given the shared ectodermal origin and integrated development of the face and the brain, facial biomarkers emerge as potential candidates to assess vulnerability for disorders in which neurodevelopment is compromised, such as schizophrenia (SZ) and bipolar disorder (BD). The sample comprised 188 individuals (67 SZ patients, 46 BD patients and 75 healthy controls (HC)). Using a landmark-based approach on 3D facial reconstructions, we quantified global and local facial shape differences between SZ/BD patients and HC using geometric morphometrics. We also assessed correlations between facial and brain cortical measures. All analyses were performed separately by sex. Diagnosis explained 4.1 % - 5.9 % of global facial shape variance in males and females with SZ, and 4.5 % - 4.1 % in BD. Regarding local facial shape, we detected 43.2 % of significantly different distances in males and 47.4 % in females with SZ as compared to HC, whereas in BD the percentages decreased to 35.8 % and 26.8 %, respectively. We detected that brain area and volume significantly explained 2.2 % and 2 % of facial shape variance in the male SZ - HC sample. Our results support facial shape as a neurodevelopmental marker for SZ and BD and reveal sex-specific pathophysiological mechanisms modulating the interplay between the brain and the face., Competing Interests: Declaration of competing interest The authors declare no conflict of interest., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

The DDR1 locus is associated with the diagnosis of schizophrenia and with processing speed in patients with schizophrenia and first-episode psychosis. Here, we investigated whether DDR1 variants are associated with bipolar disorder (BD) features. First, we performed a case‒control association study comparing DDR1 variants between patients with BD and healthy controls. Second, we performed linear regression analyses to assess the associations of DDR1 variants with neurocognitive domains and psychosocial functioning. Third, we conducted a mediation analysis to explore whether neurocognitive impairment mediated the association between DDR1 variants and psychosocial functioning in patients with BD. Finally, we studied the association between DDR1 variants and white matter microstructure. We did not find any statistically significant associations in the case‒control association study; however, we found that the combined genotypes rs1264323AA-rs2267641AC/CC were associated with worse neurocognitive performance in patients with BD with psychotic symptoms. In addition, the combined genotypes rs1264323AA-rs2267641AC/CC were associated with worse psychosocial functioning through processing speed. We did not find correlations between white matter microstructure abnormalities and the neurocognitive domains associated with the combined genotypes rs1264323AA-rs2267641AC/CC. Overall, the results suggest that DDR1 may be a marker of worse neurocognitive performance and psychosocial functioning in patients with BD, specifically those with psychotic symptoms., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Multivariate techniques better fit the anatomy of complex neuropsychiatric disorders which are characterized not by alterations in a single region, but rather by variations across distributed brain networks. Here, we used principal component analysis (PCA) to identify patterns of covariance across brain regions and relate them to clinical and demographic variables in a large generalizable dataset of individuals with bipolar disorders and controls. We then compared performance of PCA and clustering on identical sample to identify which methodology was better in capturing links between brain and clinical measures. Using data from the ENIGMA-BD working group, we investigated T1-weighted structural MRI data from 2436 participants with BD and healthy controls, and applied PCA to cortical thickness and surface area measures. We then studied the association of principal components with clinical and demographic variables using mixed regression models. We compared the PCA model with our prior clustering analyses of the same data and also tested it in a replication sample of 327 participants with BD or schizophrenia and healthy controls. The first principal component, which indexed a greater cortical thickness across all 68 cortical regions, was negatively associated with BD, BMI, antipsychotic medications, and age and was positively associated with Li treatment. PCA demonstrated superior goodness of fit to clustering when predicting diagnosis and BMI. Moreover, applying the PCA model to the replication sample yielded significant differences in cortical thickness between healthy controls and individuals with BD or schizophrenia. Cortical thickness in the same widespread regional network as determined by PCA was negatively associated with different clinical and demographic variables, including diagnosis, age, BMI, and treatment with antipsychotic medications or lithium. PCA outperformed clustering and provided an easy-to-use and interpret method to study multivariate associations between brain structure and system-level variables. PRACTITIONER POINTS: In this study of 2770 Individuals, we confirmed that cortical thickness in widespread regional networks as determined by principal component analysis (PCA) was negatively associated with relevant clinical and demographic variables, including diagnosis, age, BMI, and treatment with antipsychotic medications or lithium. Significant associations of many different system-level variables with the same brain network suggest a lack of one-to-one mapping of individual clinical and demographic factors to specific patterns of brain changes. PCA outperformed clustering analysis in the same data set when predicting group or BMI, providing a superior method for studying multivariate associations between brain structure and system-level variables., (© 2024 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

The expression of Neuritin-1 (NRN1), a neurotrophic factor crucial for neurodevelopment and synaptic plasticity, is enhanced by the Brain Derived Neurotrophic Factor (BDNF). Although the receptor of NRN1 remains unclear, it is suggested that NRN1's activation of the insulin receptor (IR) pathway promotes the transcription of the calcium voltage-gated channel subunit alpha1 C (CACNA1C). These three genes have been independently associated with schizophrenia (SZ) risk, symptomatology, and brain differences. However, research on how they synergistically modulate these phenotypes is scarce. We aimed to study whether the genetic epistasis between these genes affects the risk and clinical presentation of the disorder via its effect on brain structure. First, we tested the epistatic effect of NRN1 and BDNF or CACNA1C on (i) the risk for SZ, (ii) clinical symptoms severity and functionality (onset, PANSS, CGI and GAF), and (iii) brain cortical structure (thickness, surface area and volume measures estimated using FreeSurfer) in a sample of 86 SZ patients and 89 healthy subjects. Second, we explored whether those brain clusters influenced by epistatic effects mediate the clinical profiles. Although we did not find a direct epistatic impact on the risk, our data unveiled significant effects on the disorder's clinical presentation. Specifically, the NRN1-rs10484320 x BDNF-rs6265 interplay influenced PANSS general psychopathology, and the NRN1-rs4960155 x CACNA1C-rs1006737 interaction affected GAF scores. Moreover, several interactions between NRN1 SNPs and BDNF-rs6265 significantly influenced the surface area and cortical volume of the frontal, parietal, and temporal brain regions within patients. The NRN1-rs10484320 x BDNF-rs6265 epistasis in the left lateral orbitofrontal cortex fully mediated the effect on PANSS general psychopathology. Our study not only adds clinical significance to the well-described molecular relationship between NRN1 and BDNF but also underscores the utility of deconstructing SZ into biologically validated brain-imaging markers to explore their mediation role in the path from genetics to complex clinical manifestation., (© 2024. The Author(s).)

Purpose: This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T 2 -based pore size estimation technique., Theory and Methods: A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images., Results: Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T 2 -based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes., Conclusion: The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T 2 -based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research., (© 2024 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.)

Schizophrenia may represent a trade-off in the evolution of human-specific ontogenetic mechanisms that guide neurodevelopment. Human Accelerated Regions (HARs) are evolutionary markers functioning as neurodevelopmental transcription enhancers that have been associated with brain configuration, neural information processing, and schizophrenia risk. Here, we have investigated the influence of HARs' polygenic load on neuroanatomical measures through a case-control approach (128 patients with schizophrenia and 115 controls). To this end, we have calculated the global schizophrenia Polygenic Risk Score (Global PRS SZ ) and that specific to HARs (HARs PRS SZ ). We have also estimated the polygenic burden restricted to the HARs linked to transcriptional regulatory elements active in the foetal brain (FB-HARs PRS SZ ) and the adult brain (AB-HARs PRS SZ ). We have explored the main effects of the PRSs and the PRSs x diagnosis interactions on brain regional cortical thickness (CT) and surface area (SA). The results indicate that a higher FB-HARs PRS SZ is associated with patients' lower SA in the lateral orbitofrontal cortex, the superior temporal cortex, the pars triangularis and the paracentral lobule. While noHARs-derived PRSs show an effect on the risk, our neuroanatomical findings suggest that the human-specific transcriptional regulation during the prenatal period underlies SA variability, highlighting the role of these evolutionary markers in the schizophrenia genomic architecture., (© 2023. Springer Nature Limited.)

Introduction: Estimating the risk of manic relapse could help the psychiatrist individually adjust the treatment to the risk. Some authors have attempted to estimate this risk from baseline clinical data. Still, no studies have assessed whether the estimation could improve by adding structural magnetic resonance imaging (MRI) data. We aimed to evaluate it., Material and Methods: We followed a cohort of 78 patients with a manic episode without mixed symptoms (bipolar type I or schizoaffective disorder) at 2-4-6-9-12-15-18 months and up to 10 years. Within a cross-validation scheme, we created and evaluated a Cox lasso model to estimate the risk of manic relapse using both clinical and MRI data., Results: The model successfully estimated the risk of manic relapse (Cox regression of the time to relapse as a function of the estimated risk: hazard ratio (HR)=2.35, p=0.027; area under the curve (AUC)=0.65, expected calibration error (ECE)<0.2). The most relevant variables included in the model were the diagnosis of schizoaffective disorder, poor impulse control, unusual thought content, and cerebellum volume decrease. The estimations were poorer when we used clinical or MRI data separately., Conclusion: Combining clinical and MRI data may improve the risk of manic relapse estimation after a manic episode. We provide a website that estimates the risk according to the model to facilitate replication by independent groups before translation to clinical settings., (Copyright © 2023 The Author(s). Published by Elsevier España S.L.U. All rights reserved.)

Purpose: Biophysical models of diffusion MRI have been developed to characterize microstructure in various tissues, but existing models are not suitable for tissue composed of permeable spherical cells. In this study we introduce Cellular Exchange Imaging (CEXI), a model tailored for permeable spherical cells, and compares its performance to a related Ball & Sphere (BS) model that neglects permeability., Methods: We generated DW-MRI signals using Monte-Carlo simulations with a PGSE sequence in numerical substrates made of spherical cells and their extracellular space for a range of membrane permeability. From these signals, the properties of the substrates were inferred using both BS and CEXI models., Results: CEXI outperformed the impermeable model by providing more stable estimates cell size and intracellular volume fraction that were diffusion time-independent. Notably, CEXI accurately estimated the exchange time for low to moderate permeability levels previously reported in other studies ( κ < 25 μ m / s $$ \kappa <25\kern0.3em \mu \mathrm{m}/\mathrm{s} $$ ). However, in highly permeable substrates ( κ = 50 μ m / s $$ \kappa =50\kern0.3em \mu \mathrm{m}/\mathrm{s} $$ ), the estimated parameters were less stable, particularly the diffusion coefficients., Conclusion: This study highlights the importance of modeling the exchange time to accurately quantify microstructure properties in permeable cellular substrates. Future studies should evaluate CEXI in clinical applications such as lymph nodes, investigate exchange time as a potential biomarker of tumor severity, and develop more appropriate tissue models that account for anisotropic diffusion and highly permeable membranes., (© 2023 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals LLC on behalf of International Society for Magnetic Resonance in Medicine.)

Despite diffusion tensor imaging (DTI) evidence for widespread fractional anisotropy (FA) reductions in the brain white matter of patients with bipolar disorder, questions remain regarding the specificity and sensitivity of FA abnormalities as opposed to other diffusion metrics in the disorder. We conducted a whole-brain voxel-based multicompartment diffusion MRI study on 316 participants (i.e., 158 patients and 158 matched healthy controls) employing four diffusion metrics: the mean diffusivity (MD) and FA estimated from DTI, and the intra-axonal signal fraction (IASF) and microscopic axonal parallel diffusivity (Dpar) derived from the spherical mean technique. Our findings provide novel evidence about widespread abnormalities in other diffusion metrics in BD. An extensive overlap between the FA and IASF results suggests that the lower FA in patients may be caused by a reduced intra-axonal volume fraction or a higher macromolecular content in the intra-axonal water. We also found a diffuse alteration in MD involving white and grey matter tissue and more localised changes in Dpar. A Machine Learning analysis revealed that FA, followed by IASF, were the most helpful metric for the automatic diagnosis of BD patients, reaching an accuracy of 72%. Number of mood episodes, age of onset/duration of illness, psychotic symptoms, and current treatment with lithium, antipsychotics, antidepressants, and antiepileptics were all significantly associated with microstructure abnormalities. Lithium treatment was associated with less microstructure abnormality., (© 2023 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods., Competing Interests: Declaration of Competing Interest The authors declare no competing interests., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T 2 and T 1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T 1 -T 2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm 2 ) and multiple inversion and echo times. A second reduced diffusion-T 2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement., Competing Interests: MB was an employee of Hays plc and a consultant for F. Hoffmann-La Roche Ltd. SM was an employee and shareholder of F. Hoffmann-La Roche Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Barakovic, Pizzolato, Tax, Rudrapatna, Magon, Dyrby, Granziera, Thiran, Jones and Canales-Rodríguez.)

Sentiment analysis is the automated coding of emotions expressed in text. Sentiment analysis and other types of analyses focusing on the automatic coding of textual documents are increasingly popular in psychology and computer science. However, the potential of treating automatically coded text collected with regular sampling intervals as a signal is currently overlooked. We use the phrase "text as signal" to refer to the application of signal processing techniques to coded textual documents sampled with regularity. In order to illustrate the potential of treating text as signal, we introduce the reader to a variety of such techniques in a tutorial with two case studies in the realm of social media analysis. First, we apply finite response impulse filtering to emotion-coded tweets posted during the US Election Week of 2020 and discuss the visualization of the resulting variation in the filtered signal. We use changepoint detection to highlight the important changes in the emotional signals. Then we examine data interpolation, analysis of periodicity via the fast Fourier transform (FFT), and FFT filtering to personal value-coded tweets from November 2019 to October 2020 and link the variation in the filtered signal to some of the epoch-defining events occurring during this period. Finally, we use block bootstrapping to estimate the variability/uncertainty in the resulting filtered signals. After working through the tutorial, the readers will understand the basics of signal processing to analyze regularly sampled coded text., (© 2022. The Author(s).)

Monte-Carlo diffusion simulations are a powerful tool for validating tissue microstructure models by generating synthetic diffusion-weighted magnetic resonance images (DW-MRI) in controlled environments. This is fundamental for understanding the link between micrometre-scale tissue properties and DW-MRI signals measured at the millimetre-scale, optimizing acquisition protocols to target microstructure properties of interest, and exploring the robustness and accuracy of estimation methods. However, accurate simulations require substrates that reflect the main microstructural features of the studied tissue. To address this challenge, we introduce a novel computational workflow, CACTUS (Computational Axonal Configurator for Tailored and Ultradense Substrates), for generating synthetic white matter substrates. Our approach allows constructing substrates with higher packing density than existing methods, up to 95% intra-axonal volume fraction, and larger voxel sizes of up to 500μm 3 with rich fibre complexity. CACTUS generates bundles with angular dispersion, bundle crossings, and variations along the fibres of their inner and outer radii and g-ratio. We achieve this by introducing a novel global cost function and a fibre radial growth approach that allows substrates to match predefined targeted characteristics and mirror those reported in histological studies. CACTUS improves the development of complex synthetic substrates, paving the way for future applications in microstructure imaging., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Villarreal-Haro, Gardier, Canales-Rodríguez, Fischi-Gomez, Girard, Thiran and Rafael-Patiño.)

We enable the estimation of the per-axon axial diffusivity from single encoding, strongly diffusion-weighted, pulsed gradient spin echo data. Additionally, we improve the estimation of the per-axon radial diffusivity compared to estimates based on spherical averaging. The use of strong diffusion weightings in magnetic resonance imaging (MRI) allows to approximate the signal in white matter as the sum of the contributions from only axons. At the same time, spherical averaging leads to a major simplification of the modeling by removing the need to explicitly account for the unknown distribution of axonal orientations. However, the spherically averaged signal acquired at strong diffusion weightings is not sensitive to the axial diffusivity, which cannot therefore be estimated although needed for modeling axons - especially in the context of multi-compartmental modeling. We introduce a new general method for the estimation of both the axial and radial axonal diffusivities at strong diffusion weightings based on kernel zonal modeling. The method could lead to estimates that are free from partial volume bias with gray matter or other isotropic compartments. The method is tested on publicly available data from the MGH Adult Diffusion Human Connectome project. We report reference values of axonal diffusivities based on 34 subjects, and derive estimates of axonal radii from only two shells. The estimation problem is also addressed from the angle of the required data preprocessing, the presence of biases related to modeling assumptions, current limitations, and future possibilities., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Marco Pizzolato reports financial support was provided by EUFramework Programme for Research and Innovation Marie Sklodowska-Curie Actions. Erick Jorge Canales-Rodríguez reports financial support was provided by Swiss National Science Foundation. Mariam Andersson reports financial support was provided by Capital Region of Denmark. Tim B. Dyrby reports financial support was provided by the European Research Council., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Objective: A major limitation of current suicide research is the lack of power to identify robust correlates of suicidal thoughts or behavior. Variation in suicide risk assessment instruments used across cohorts may represent a limitation to pooling data in international consortia., Method: Here, we examine this issue through two approaches: (a) an extensive literature search on the reliability and concurrent validity of the most commonly used instruments and (b) by pooling data (N ∼ 6,000 participants) from cohorts from the Enhancing NeuroImaging Genetics Through Meta-Analysis (ENIGMA) Major Depressive Disorder and ENIGMA-Suicidal Thoughts and Behaviour working groups, to assess the concurrent validity of instruments currently used for assessing suicidal thoughts or behavior., Results: We observed moderate-to-high correlations between measures, consistent with the wide range (κ range: 0.15-0.97; r range: 0.21-0.94) reported in the literature. Two common multi-item instruments, the Columbia Suicide Severity Rating Scale and the Beck Scale for Suicidal Ideation were highly correlated with each other (r = 0.83). Sensitivity analyses identified sources of heterogeneity such as the time frame of the instrument and whether it relies on self-report or a clinical interview. Finally, construct-specific analyses suggest that suicide ideation items from common psychiatric questionnaires are most concordant with the suicide ideation construct of multi-item instruments., Conclusions: Our findings suggest that multi-item instruments provide valuable information on different aspects of suicidal thoughts or behavior but share a modest core factor with single suicidal ideation items. Retrospective, multisite collaborations including distinct instruments should be feasible provided they harmonize across instruments or focus on specific constructs of suicidality. (PsycInfo Database Record (c) 2023 APA, all rights reserved).

Background: Obesity is highly prevalent and disabling, especially in individuals with severe mental illness including bipolar disorders (BD). The brain is a target organ for both obesity and BD. Yet, we do not understand how cortical brain alterations in BD and obesity interact., Methods: We obtained body mass index (BMI) and MRI-derived regional cortical thickness, surface area from 1231 BD and 1601 control individuals from 13 countries within the ENIGMA-BD Working Group. We jointly modeled the statistical effects of BD and BMI on brain structure using mixed effects and tested for interaction and mediation. We also investigated the impact of medications on the BMI-related associations., Results: BMI and BD additively impacted the structure of many of the same brain regions. Both BMI and BD were negatively associated with cortical thickness, but not surface area. In most regions the number of jointly used psychiatric medication classes remained associated with lower cortical thickness when controlling for BMI. In a single region, fusiform gyrus, about a third of the negative association between number of jointly used psychiatric medications and cortical thickness was mediated by association between the number of medications and higher BMI., Conclusions: We confirmed consistent associations between higher BMI and lower cortical thickness, but not surface area, across the cerebral mantle, in regions which were also associated with BD. Higher BMI in people with BD indicated more pronounced brain alterations. BMI is important for understanding the neuroanatomical changes in BD and the effects of psychiatric medications on the brain.

Tractography enables identifying and evaluating the healthy and diseased brain's white matter pathways from diffusion-weighted magnetic resonance imaging data. As previous evaluation studies have reported significant false-positive estimation biases, recent microstructure-informed tractography algorithms have been introduced to improve the trade-off between specificity and sensitivity. However, a major limitation for characterizing the performance of these techniques is the lack of ground truth brain data. In this study, we compared the performance of two relevant microstructure-informed tractography methods, SIFT2 and COMMIT, by assessing the subject specificity and reproducibility of their derived white matter pathways. Specifically, twenty healthy young subjects were scanned at eight different time points at two different sites. Subject specificity and reproducibility were evaluated using the whole-brain connectomes and a subset of 29 white matter bundles. Our results indicate that although the raw tractograms are more vulnerable to the presence of false-positive connections, they are highly reproducible, suggesting that the estimation bias is subject-specific. This high reproducibility was preserved when microstructure-informed tractography algorithms were used to filter the raw tractograms. Moreover, the resulting track-density images depicted a more uniform coverage of streamlines throughout the white matter, suggesting that these techniques could increase the biological meaning of the estimated fascicles. Notably, we observed an increased subject specificity by employing connectivity pre-processing techniques to reduce the underlaying noise and the data dimensionality (using principal component analysis), highlighting the importance of these tools for future studies. Finally, no strong bias from the scanner site or time between measurements was found. The largest intraindividual variance originated from the sole repetition of data measurements (inter-run)., Competing Interests: Conflict of interest G.F.P., T.H., and T.K. are employees of Siemens Healthcare AG, Switzerland., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022. Published by Elsevier Inc.)

Background: Morphology of the human cerebral cortex differs across psychiatric disorders, with neurobiology and developmental origins mostly undetermined. Deviations in the tangential growth of the cerebral cortex during pre/perinatal periods may be reflected in individual variations in cortical surface area later in life., Methods: Interregional profiles of group differences in surface area between cases and controls were generated using T1-weighted magnetic resonance imaging from 27,359 individuals including those with attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, schizophrenia, and high general psychopathology (through the Child Behavior Checklist). Similarity of interregional profiles of group differences in surface area and prenatal cell-specific gene expression was assessed., Results: Across the 11 cortical regions, group differences in cortical area for attention-deficit/hyperactivity disorder, schizophrenia, and Child Behavior Checklist were dominant in multimodal association cortices. The same interregional profiles were also associated with interregional profiles of (prenatal) gene expression specific to proliferative cells, namely radial glia and intermediate progenitor cells (greater expression, larger difference), as well as differentiated cells, namely excitatory neurons and endothelial and mural cells (greater expression, smaller difference). Finally, these cell types were implicated in known pre/perinatal risk factors for psychosis. Genes coexpressed with radial glia were enriched with genes implicated in congenital abnormalities, birth weight, hypoxia, and starvation. Genes coexpressed with endothelial and mural genes were enriched with genes associated with maternal hypertension and preterm birth., Conclusions: Our findings support a neurodevelopmental model of vulnerability to mental illness whereby prenatal risk factors acting through cell-specific processes lead to deviations from typical brain development during pregnancy., (Copyright © 2022 Society of Biological Psychiatry. Published by Elsevier Inc. All rights reserved.)

Aims: Rates of obesity have reached epidemic proportions, especially among people with psychiatric disorders. While the effects of obesity on the brain are of major interest in medicine, they remain markedly under-researched in psychiatry., Methods: We obtained body mass index (BMI) and magnetic resonance imaging-derived regional cortical thickness, surface area from 836 bipolar disorders (BD) and 1600 control individuals from 14 sites within the ENIGMA-BD Working Group. We identified regionally specific profiles of cortical thickness using K-means clustering and studied clinical characteristics associated with individual cortical profiles., Results: We detected two clusters based on similarities among participants in cortical thickness. The lower thickness cluster (46.8% of the sample) showed thinner cortex, especially in the frontal and temporal lobes and was associated with diagnosis of BD, higher BMI, and older age. BD individuals in the low thickness cluster were more likely to have the diagnosis of bipolar disorder I and less likely to be treated with lithium. In contrast, clustering based on similarities in the cortical surface area was unrelated to BD or BMI and only tracked age and sex., Conclusions: We provide evidence that both BD and obesity are associated with similar alterations in cortical thickness, but not surface area. The fact that obesity increased the chance of having low cortical thickness could explain differences in cortical measures among people with BD. The thinner cortex in individuals with higher BMI, which was additive and similar to the BD-associated alterations, may suggest that treating obesity could lower the extent of cortical thinning in BD., (© 2021 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Quantitative magnetic resonance imaging (qMRI) can increase the specificity and sensitivity of conventional weighted MRI to underlying pathology by comparing meaningful physical or chemical parameters, measured in physical units, with normative values acquired in a healthy population. This study focuses on multi-echo T 2 relaxometry, a qMRI technique that probes the complex tissue microstructure by differentiating compartment-specific T 2 relaxation times. However, estimation methods are still limited by their sensitivity to the underlying noise. Moreover, estimating the model's parameters is challenging because the resulting inverse problem is ill-posed, requiring advanced numerical regularization techniques. As a result, the estimates from distinct regularization strategies are different. In this work, we aimed to investigate the variability and reproducibility of different techniques for estimating the transverse relaxation time of the intra- and extra-cellular space ( T 2 I E ) in gray (GM) and white matter (WM) tissue in a clinical setting, using a multi-site, multi-session, and multi-run T 2 relaxometry dataset. To this end, we evaluated three different techniques for estimating the T 2 spectra (two regularized non-negative least squares methods and a machine learning approach). Two independent analyses were performed to study the effect of using raw and denoised data. For both the GM and WM regions, and the raw and denoised data, our results suggest that the principal source of variance is the inter-subject variability, showing a higher coefficient of variation (CoV) than those estimated for the inter-site, inter-session, and inter-run, respectively. For all reconstruction methods studied, the CoV ranged between 0.32 and 1.64%. Interestingly, the inter-session variability was close to the inter-scanner variability with no statistical differences, suggesting that T 2 I E is a robust parameter that could be employed in multi-site neuroimaging studies. Furthermore, the three tested methods showed consistent results and similar intra-class correlation (ICC), with values superior to 0.7 for most regions. Results from raw data were slightly more reproducible than those from denoised data. The regularized non-negative least squares method based on the L-curve technique produced the best results, with ICC values ranging from 0.72 to 0.92., Competing Interests: Authors GP, TH, TY, and TK are employees of Siemens Healthineers International AG, Switzerland. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Fischi-Gomez, Girard, Koch, Yu, Pizzolato, Brügger, Piredda, Hilbert, Cadic-Melchior, Beanato, Park, Morishita, Wessel, Schiavi, Daducci, Kober, Canales-Rodríguez, Hummel and Thiran.)

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T 1 and T 2 relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T 1 /T 2 data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions., (© 2021 John Wiley & Sons, Ltd.)

Previous cross-sectional research has largely associated binge drinking (BD) with changes in volume and thickness during adolescence and early adulthood. Nevertheless, the long-term alcohol-related effects on gray matter features in youths who had maintained a BD pattern over time have not yet been sufficiently explored. The present study aimed to assess group differences both cross-sectionally and longitudinally [using symmetric percent change (SPC)] on several structural measures (i.e., thickness, surface area, volume). For this purpose, magnetic resonance imaging was recorded twice within a 2-year interval; at baseline (18-19 years) and a follow-up (20-21 years). The sample included 44 university students who were classified as 16 stable binge drinkers (8 females) and 28 stable controls (13 females). Whole-brain analysis showed larger insular surface area in binge drinkers relative to controls at follow-up (cluster-wise p = 0.045). On the other hand, region of interest (ROI) analyses on thickness also revealed a group by sex interaction at follow-up ( p = 0.005), indicating that BD males had smaller right rostral middle frontal gyrus thickness than both control males ( p = 0.011) and BD females ( p = 0.029). Similarly, ROI-based analysis on longitudinal data showed a group by sex interaction in the right nucleus accumbens ( p = 0.009) which revealed a decreased volume across time in BD males than in control males ( p = 0.007). Overall, continued BD pattern during emerging adulthood appears to lead to gray matter abnormalities in regions intimately involved in reward processing, emotional regulation and executive functions. Notably, some anomalies varied significantly depending on sex, suggesting a sex-specific impact of BD on typical neurodevelopment processes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Pérez-García, Cadaveira, Canales-Rodríguez, Suárez-Suárez, Rodríguez Holguín, Corral, Blanco-Ramos and Doallo.)

Fetal brain diffusion magnetic resonance images (MRI) are often acquired with a lower through-plane than in-plane resolution. This anisotropy is often overcome by classical upsampling methods such as linear or cubic interpolation. In this work, we employ an unsupervised learning algorithm using an autoencoder neural network for single-image through-plane super-resolution by leveraging a large amount of data. Our framework, which can also be used for slice outliers replacement, overperformed conventional interpolations quantitatively and qualitatively on pre-term newborns of the developing Human Connectome Project. The evaluation was performed on both the original diffusion-weighted signal and the estimated diffusion tensor maps. A byproduct of our autoencoder was its ability to act as a denoiser. The network was able to generalize fetal data with different levels of motions and we qualitatively showed its consistency, hence supporting the relevance of pre-term datasets to improve the processing of fetal brain images., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Kebiri, Canales-Rodríguez, Lajous, de Dumast, Girard, Alemán-Gómez, Koob, Jakab and Bach Cuadra.)

Background: Bipolar disorder (BD) is associated with cortical and subcortical structural brain abnormalities. It is unclear whether such alterations progressively change over time, and how this is related to the number of mood episodes. To address this question, we analyzed a large and diverse international sample with longitudinal magnetic resonance imaging (MRI) and clinical data to examine structural brain changes over time in BD., Methods: Longitudinal structural MRI and clinical data from the ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) BD Working Group, including 307 patients with BD and 925 healthy control subjects, were collected from 14 sites worldwide. Male and female participants, aged 40 ± 17 years, underwent MRI at 2 time points. Cortical thickness, surface area, and subcortical volumes were estimated using FreeSurfer. Annualized change rates for each imaging phenotype were compared between patients with BD and healthy control subjects. Within patients, we related brain change rates to the number of mood episodes between time points and tested for effects of demographic and clinical variables., Results: Compared with healthy control subjects, patients with BD showed faster enlargement of ventricular volumes and slower thinning of the fusiform and parahippocampal cortex (0.18

In magnetic resonance imaging, the application of a strong diffusion weighting suppresses the signal contributions from the less diffusion-restricted constituents of the brain's white matter, thus enabling the estimation of the transverse relaxation time T 2 that arises from the more diffusion-restricted constituents such as the axons. However, the presence of cell nuclei and vacuoles can confound the estimation of the axonal T 2 , as diffusion within those structures is also restricted, causing the corresponding signal to survive the strong diffusion weighting. We devise an estimator of the axonal T 2 based on the directional spherical variance of the strongly diffusion-weighted signal. The spherical variance T 2 estimates are insensitive to the presence of isotropic contributions to the signal like those provided by cell nuclei and vacuoles. We show that with a strong diffusion weighting these estimates differ from those obtained using the directional spherical mean of the signal which contains both axonal and isotropically-restricted contributions. Our findings hint at the presence of an MRI-visible isotropically-restricted contribution to the signal in the white matter ex vivo fixed tissue (monkey) at 7T, and do not allow us to discard such a possibility also for in vivo human data collected with a clinical 3T system., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Delineating the association of age and cortical thickness in healthy individuals is critical given the association of cortical thickness with cognition and behavior. Previous research has shown that robust estimates of the association between age and brain morphometry require large-scale studies. In response, we used cross-sectional data from 17,075 individuals aged 3-90 years from the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to infer age-related changes in cortical thickness. We used fractional polynomial (FP) regression to quantify the association between age and cortical thickness, and we computed normalized growth centiles using the parametric Lambda, Mu, and Sigma method. Interindividual variability was estimated using meta-analysis and one-way analysis of variance. For most regions, their highest cortical thickness value was observed in childhood. Age and cortical thickness showed a negative association; the slope was steeper up to the third decade of life and more gradual thereafter; notable exceptions to this general pattern were entorhinal, temporopolar, and anterior cingulate cortices. Interindividual variability was largest in temporal and frontal regions across the lifespan. Age and its FP combinations explained up to 59% variance in cortical thickness. These results may form the basis of further investigation on normative deviation in cortical thickness and its significance for behavioral and cognitive outcomes., (© 2021 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

For many traits, males show greater variability than females, with possible implications for understanding sex differences in health and disease. Here, the ENIGMA (Enhancing Neuro Imaging Genetics through Meta-Analysis) Consortium presents the largest-ever mega-analysis of sex differences in variability of brain structure, based on international data spanning nine decades of life. Subcortical volumes, cortical surface area and cortical thickness were assessed in MRI data of 16,683 healthy individuals 1-90 years old (47% females). We observed significant patterns of greater male than female between-subject variance for all subcortical volumetric measures, all cortical surface area measures, and 60% of cortical thickness measures. This pattern was stable across the lifespan for 50% of the subcortical structures, 70% of the regional area measures, and nearly all regions for thickness. Our findings that these sex differences are present in childhood implicate early life genetic or gene-environment interaction mechanisms. The findings highlight the importance of individual differences within the sexes, that may underpin sex-specific vulnerability to disorders., (© 2020 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

The hippocampus consists of anatomically and functionally distinct subfields that may be differentially involved in the pathophysiology of bipolar disorder (BD). Here we, the Enhancing NeuroImaging Genetics through Meta-Analysis Bipolar Disorder workinggroup, study hippocampal subfield volumetry in BD. T1-weighted magnetic resonance imaging scans from 4,698 individuals (BD = 1,472, healthy controls [HC] = 3,226) from 23 sites worldwide were processed with FreeSurfer. We used linear mixed-effects models and mega-analysis to investigate differences in hippocampal subfield volumes between BD and HC, followed by analyses of clinical characteristics and medication use. BD showed significantly smaller volumes of the whole hippocampus (Cohen's d = -0.20), cornu ammonis (CA)1 (d = -0.18), CA2/3 (d = -0.11), CA4 (d = -0.19), molecular layer (d = -0.21), granule cell layer of dentate gyrus (d = -0.21), hippocampal tail (d = -0.10), subiculum (d = -0.15), presubiculum (d = -0.18), and hippocampal amygdala transition area (d = -0.17) compared to HC. Lithium users did not show volume differences compared to HC, while non-users did. Antipsychotics or antiepileptic use was associated with smaller volumes. In this largest study of hippocampal subfields in BD to date, we show widespread reductions in nine of 12 subfields studied. The associations were modulated by medication use and specifically the lack of differences between lithium users and HC supports a possible protective role of lithium in BD., (© 2020 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

First-degree relatives of patients diagnosed with schizophrenia (SZ-FDRs) show similar patterns of brain abnormalities and cognitive alterations to patients, albeit with smaller effect sizes. First-degree relatives of patients diagnosed with bipolar disorder (BD-FDRs) show divergent patterns; on average, intracranial volume is larger compared to controls, and findings on cognitive alterations in BD-FDRs are inconsistent. Here, we performed a meta-analysis of global and regional brain measures (cortical and subcortical), current IQ, and educational attainment in 5,795 individuals (1,103 SZ-FDRs, 867 BD-FDRs, 2,190 controls, 942 schizophrenia patients, 693 bipolar patients) from 36 schizophrenia and/or bipolar disorder family cohorts, with standardized methods. Compared to controls, SZ-FDRs showed a pattern of widespread thinner cortex, while BD-FDRs had widespread larger cortical surface area. IQ was lower in SZ-FDRs (d = -0.42, p = 3 × 10 -5 ), with weak evidence of IQ reductions among BD-FDRs (d = -0.23, p = .045). Both relative groups had similar educational attainment compared to controls. When adjusting for IQ or educational attainment, the group-effects on brain measures changed, albeit modestly. Changes were in the expected direction, with less pronounced brain abnormalities in SZ-FDRs and more pronounced effects in BD-FDRs. To conclude, SZ-FDRs and BD-FDRs show a differential pattern of structural brain abnormalities. In contrast, both had lower IQ scores and similar school achievements compared to controls. Given that brain differences between SZ-FDRs and BD-FDRs remain after adjusting for IQ or educational attainment, we suggest that differential brain developmental processes underlying predisposition for schizophrenia or bipolar disorder are likely independent of general cognitive impairment., (© 2020 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

Age has a major effect on brain volume. However, the normative studies available are constrained by small sample sizes, restricted age coverage and significant methodological variability. These limitations introduce inconsistencies and may obscure or distort the lifespan trajectories of brain morphometry. In response, we capitalized on the resources of the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) Consortium to examine age-related trajectories inferred from cross-sectional measures of the ventricles, the basal ganglia (caudate, putamen, pallidum, and nucleus accumbens), the thalamus, hippocampus and amygdala using magnetic resonance imaging data obtained from 18,605 individuals aged 3-90 years. All subcortical structure volumes were at their maximum value early in life. The volume of the basal ganglia showed a monotonic negative association with age thereafter; there was no significant association between age and the volumes of the thalamus, amygdala and the hippocampus (with some degree of decline in thalamus) until the sixth decade of life after which they also showed a steep negative association with age. The lateral ventricles showed continuous enlargement throughout the lifespan. Age was positively associated with inter-individual variability in the hippocampus and amygdala and the lateral ventricles. These results were robust to potential confounders and could be used to examine the functional significance of deviations from typical age-related morphometric patterns., (© 2021 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

MRI-derived brain measures offer a link between genes, the environment and behavior and have been widely studied in bipolar disorder (BD). However, many neuroimaging studies of BD have been underpowered, leading to varied results and uncertainty regarding effects. The Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Bipolar Disorder Working Group was formed in 2012 to empower discoveries, generate consensus findings and inform future hypothesis-driven studies of BD. Through this effort, over 150 researchers from 20 countries and 55 institutions pool data and resources to produce the largest neuroimaging studies of BD ever conducted. The ENIGMA Bipolar Disorder Working Group applies standardized processing and analysis techniques to empower large-scale meta- and mega-analyses of multimodal brain MRI and improve the replicability of studies relating brain variation to clinical and genetic data. Initial BD Working Group studies reveal widespread patterns of lower cortical thickness, subcortical volume and disrupted white matter integrity associated with BD. Findings also include mapping brain alterations of common medications like lithium, symptom patterns and clinical risk profiles and have provided further insights into the pathophysiological mechanisms of BD. Here we discuss key findings from the BD working group, its ongoing projects and future directions for large-scale, collaborative studies of mental illness., (© 2020 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)

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., (© 2021 The Authors. Human Brain Mapping published by Wiley Periodicals LLC.)