White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

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.

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 PRSSZ) and that specific to HARs (HARs PRSSZ). We have also estimated the polygenic burden restricted to the HARs linked to transcriptional regulatory elements active in the foetal brain (FB-HARs PRSSZ) and the adult brain (AB-HARs PRSSZ). 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 PRSSZ 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.

A relatively large number of studies have investigated the power of structural magnetic resonance imaging (sMRI) data to discriminate patients with schizophrenia from healthy controls. However, very few of them have also included patients with bipolar disorder, allowing the clinically relevant discrimination between both psychotic diagnostics. To assess the efficacy of sMRI data for diagnostic prediction in psychosis we objectively evaluated the discriminative power of a wide range of commonly used machine learning algorithms (ridge, lasso, elastic net and L0 norm regularized logistic regressions, a support vector classifier, regularized discriminant analysis, random forests and a Gaussian process classifier) on main sMRI features including grey and white matter voxel-based morphometry (VBM), vertex-based cortical thickness and volume, region of interest volumetric measures and wavelet-based morphometry (WBM) maps. All possible combinations of algorithms and data features were considered in pairwise classifications of matched samples of healthy controls (N = 127), patients with schizophrenia (N = 128) and patients with bipolar disorder (N = 128). Results show that the selection of feature type is important, with grey matter VBM (without data reduction) delivering the best diagnostic prediction rates (averaging over classifiers: schizophrenia vs. healthy 75%, bipolar disorder vs. healthy 63% and schizophrenia vs. bipolar disorder 62%) whereas algorithms usually yielded very similar results. Indeed, those grey matter VBM accuracy rates were not even improved by combining all feature types in a single prediction model. Further multi-class classifications considering the three groups simultaneously made evident a lack of predictive power for the bipolar group, probably due to its intermediate anatomical features, located between those observed in healthy controls and those found in patients with schizophrenia. Finally, we provide MRIPredict (https://www.nitrc.org/projects/mripredict/), a free tool for SPM, FSL and R, to easily carry out voxelwise predictions based on VBM images.

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.

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.

The aim of this study was to identify differential global and local brain structural patterns in Dravet Syndrome (DS) patients as compared with a control subject group, using brain morphometry techniques which provide a quantitative whole-brain structural analysis that allows for specific patterns to be generalized across series of individuals. Nine patients with the diagnosis of DS that tested positive for mutation in the SCN1A gene and nine well-matched healthy controls were investigated using voxel brain morphometry (VBM), cortical thickness and cortical gyrification measurements. Global volume reductions of gray matter (GM) and white matter (WM) were related to DS. Local volume reductions corresponding to several white matter regions in brainstem, cerebellum, corpus callosum, corticospinal tracts and association fibers (left inferior fronto-occipital fasciculus and left uncinate fasciculus) were also found. Furthermore, DS showed a reduced cortical folding in the right precentral gyrus. The present findings describe DS-related brain structure abnormalities probably linked to the expression of the SCN1A mutation.

Our goal is to study the human brain anatomical network. For this, the anatomical connection probabilities (ACP) between 90 cortical and subcortical brain gray matter areas are estimated from diffusion-weighted Magnetic Resonance Imaging (DW-MRI) techniques. The ACP between any two areas gives the probability that those areas are connected at least by a single nervous fiber. Then, the brain is modeled as a non-directed weighted graph with continuous arc weights given by the ACP matrix. Based on this approach, complex networks properties such as small-world attributes, efficiency, degree distribution, vulnerability, betweenness centrality and motifs composition are studied. The analysis was carried out for 20 right-handed healthy subjects (mean age: 31.10, S.D.: 7.43). According to the results, all networks have small-world and broad-scale characteristics. Additionally, human brain anatomical networks present bigger local efficiency and smaller global efficiency than the corresponding random networks. In a vulnerability and betweenness centrality analysis, the most indispensable and critical anatomical areas were identified: putamens, precuneus, insulas, superior parietals and superior frontals. Interestingly, some areas have a negative vulnerability (e.g. superior temporal poles, pallidums, supramarginals and hechls), which suggest that even at the cost of losing in global anatomical efficiency, these structures were maintained through the evolutionary processes due to their important functions. Finally, symmetrical characteristic building blocks (motifs) of size 3 and 4 were calculated, obtaining that motifs of size 4 are the expanded version of motif of size 3. These results are in agreement with previous anatomical studies in the cat and macaque cerebral cortex.

Individuals with bipolar disorders (BD) frequently suffer from obesity, which is often associated with neurostructural alterations. Yet, the effects of obesity on brain structure in BD are under-researched. We obtained MRI-derived brain subcortical volumes and body mass index (BMI) from 1134 BD and 1601 control individuals from 17 independent research sites within the ENIGMA-BD Working Group. We jointly modeled the effects of BD and BMI on subcortical volumes using mixed-effects modeling and tested for mediation of group differences by obesity using nonparametric bootstrapping. All models controlled for age, sex, hemisphere, total intracranial volume, and data collection site. Relative to controls, individuals with BD had significantly higher BMI, larger lateral ventricular volume, and smaller volumes of amygdala, hippocampus, pallidum, caudate, and thalamus. BMI was positively associated with ventricular and amygdala and negatively with pallidal volumes. When analyzed jointly, both BD and BMI remained associated with volumes of lateral ventricles  and amygdala. Adjusting for BMI decreased the BD vs control differences in ventricular volume. Specifically, 18.41% of the association between BD and ventricular volume was mediated by BMI (Z = 2.73, p = 0.006). BMI was associated with similar regional brain volumes as BD, including lateral ventricles, amygdala, and pallidum. Higher BMI may in part account for larger ventricles, one of the most replicated findings in BD. Comorbidity with obesity could explain why neurostructural alterations are more pronounced in some individuals with BD. Future prospective brain imaging studies should investigate whether obesity could be a modifiable risk factor for neuroprogression.

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.

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 profiles of cortical abnormalities in schizophrenia and bipolar disorder, and how far they resemble each other, have only been studied to a limited extent. The aim of this study was to identify and compare the changes in cortical morphology associated with these pathologies.

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

Microstructure imaging from diffusion magnetic resonance (MR) data represents an invaluable tool to study non-invasively the morphology of tissues and to provide a biological insight into their microstructural organization. In recent years, a variety of biophysical models have been proposed to associate particular patterns observed in the measured signal with specific microstructural properties of the neuronal tissue, such as axon diameter and fiber density. Despite very appealing results showing that the estimated microstructure indices agree very well with histological examinations, existing techniques require computationally very expensive non-linear procedures to fit the models to the data which, in practice, demand the use of powerful computer clusters for large-scale applications. In this work, we present a general framework for Accelerated Microstructure Imaging via Convex Optimization (AMICO) and show how to re-formulate this class of techniques as convenient linear systems which, then, can be efficiently solved using very fast algorithms. We demonstrate this linearization of the fitting problem for two specific models, i.e. ActiveAx and NODDI, providing a very attractive alternative for parameter estimation in those techniques; however, the AMICO framework is general and flexible enough to work also for the wider space of microstructure imaging methods. Results demonstrate that AMICO represents an effective means to accelerate the fit of existing techniques drastically (up to four orders of magnitude faster) while preserving accuracy and precision in the estimated model parameters (correlation above 0.9). We believe that the availability of such ultrafast algorithms will help to accelerate the spread of microstructure imaging to larger cohorts of patients and to study a wider spectrum of neurological disorders.

Fronto-limbic white matter (WM) abnormalities are assumed to lie at the heart of the pathophysiology of bipolar disorder (BD); however, diffusion tensor imaging (DTI) studies have reported heterogeneous results and it is not clear how the clinical heterogeneity is related to the observed differences. This study aimed to identify WM abnormalities that differentiate patients with BD from healthy controls (HC) in the largest DTI dataset of patients with BD to date, collected via the ENIGMA network. We gathered individual tensor-derived regional metrics from 26 cohorts leading to a sample size of N = 3033 (1482 BD and 1551 HC). Mean fractional anisotropy (FA) from 43 regions of interest (ROI) and average whole-brain FA were entered into univariate mega- and meta-analyses to differentiate patients with BD from HC. Mega-analysis revealed significantly lower FA in patients with BD compared with HC in 29 regions, with the highest effect sizes observed within the corpus callosum (R2 = 0.041, Pcorr < 0.001) and cingulum (right: R2 = 0.041, left: R2 = 0.040, Pcorr < 0.001). Lithium medication, later onset and short disease duration were related to higher FA along multiple ROIs. Results of the meta-analysis showed similar effects. We demonstrated widespread WM abnormalities in BD and highlighted that altered WM connectivity within the corpus callosum and the cingulum are strongly associated with BD. These brain abnormalities could represent a biomarker for use in the diagnosis of BD. Interactive three-dimensional visualization of the results is available at www.enigma-viewer.org.

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.

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.

Spherical deconvolution (SD) methods are widely used to estimate the intra-voxel white-matter fiber orientations from diffusion MRI data. However, while some of these methods assume a zero-mean Gaussian distribution for the underlying noise, its real distribution is known to be non-Gaussian and to depend on many factors such as the number of coils and the methodology used to combine multichannel MRI signals. Indeed, the two prevailing methods for multichannel signal combination lead to noise patterns better described by Rician and noncentral Chi distributions. Here we develop a Robust and Unbiased Model-BAsed Spherical Deconvolution (RUMBA-SD) technique, intended to deal with realistic MRI noise, based on a Richardson-Lucy (RL) algorithm adapted to Rician and noncentral Chi likelihood models. To quantify the benefits of using proper noise models, RUMBA-SD was compared with dRL-SD, a well-established method based on the RL algorithm for Gaussian noise. Another aim of the study was to quantify the impact of including a total variation (TV) spatial regularization term in the estimation framework. To do this, we developed TV spatially-regularized versions of both RUMBA-SD and dRL-SD algorithms. The evaluation was performed by comparing various quality metrics on 132 three-dimensional synthetic phantoms involving different inter-fiber angles and volume fractions, which were contaminated with noise mimicking patterns generated by data processing in multichannel scanners. The results demonstrate that the inclusion of proper likelihood models leads to an increased ability to resolve fiber crossings with smaller inter-fiber angles and to better detect non-dominant fibers. The inclusion of TV regularization dramatically improved the resolution power of both techniques. The above findings were also verified in human brain data.

Bipolar disorders (BDs) are among the leading causes of morbidity and disability. Objective biological markers, such as those based on brain imaging, could aid in clinical management of BD. Machine learning (ML) brings neuroimaging analyses to individual subject level and may potentially allow for their diagnostic use. However, fair and optimal application of ML requires large, multi-site datasets. We applied ML (support vector machines) to MRI data (regional cortical thickness, surface area, subcortical volumes) from 853 BD and 2167 control participants from 13 cohorts in the ENIGMA consortium. We attempted to differentiate BD from control participants, investigated different data handling strategies and studied the neuroimaging/clinical features most important for classification. Individual site accuracies ranged from 45.23% to 81.07%. Aggregate subject-level analyses yielded the highest accuracy (65.23%, 95% CI = 63.47-67.00, ROC-AUC = 71.49%, 95% CI = 69.39-73.59), followed by leave-one-site-out cross-validation (accuracy = 58.67%, 95% CI = 56.70-60.63). Meta-analysis of individual site accuracies did not provide above chance results. There was substantial agreement between the regions that contributed to identification of BD participants in the best performing site and in the aggregate dataset (Cohen's Kappa = 0.83, 95% CI = 0.829-0.831). Treatment with anticonvulsants and age were associated with greater odds of correct classification. Although short of the 80% clinically relevant accuracy threshold, the results are promising and provide a fair and realistic estimate of classification performance, which can be achieved in a large, ecologically valid, multi-site sample of BD participants based on regional neurostructural measures. Furthermore, the significant classification in different samples was based on plausible and similar neuroanatomical features. Future multi-site studies should move towards sharing of raw/voxelwise neuroimaging data.

Background: Structural imaging studies of cannabis users have found evidence of both cortical and subcortical volume reductions, especially in cannabinoid receptor-rich regions such as the hippocampus and amygdala. However, the findings have not been consistent. In the present study, we examined a sample of adult heavy cannabis users without other substance abuse to determine whether long-term use is associated with brain structural changes, especially in the subcortical regions. Method: We compared the gray matter volume of 14 long-term, heavy cannabis users with non-using controls. To provide robust findings, we conducted two separate studies using two different MRI techniques. Each study used the same sample of cannabis users and a different control group, respectively. Both control groups were independent of each other. First, whole-brain voxel-based morphometry (VBM) was used to compare the cannabis users against 28 matched controls (HC1 group). Second, a volumetric analysis of subcortical regions was performed to assess differences between the cannabis users and a sample of 100 matched controls (HC2 group) obtained from a local database of healthy volunteers. Results: The VBM study revealed that, compared to the control group HC1, the cannabis users did not show cortical differences nor smaller volume in any subcortical structure but showed a cluster (p < 0.001) of larger GM volume in the basal ganglia, involving the caudate, putamen, pallidum, and nucleus accumbens, bilaterally. The subcortical volumetric analysis revealed that, compared to the control group HC2, the cannabis users showed significantly larger volumes in the putamen (p = 0.001) and pallidum (p = 0.0015). Subtle trends, only significant at the uncorrected level, were also found in the caudate (p = 0.05) and nucleus accumbens (p = 0.047). Conclusions: This study does not support previous findings of hippocampal and/or amygdala structural changes in long-term, heavy cannabis users. It does, however, provide evidence of basal ganglia volume increases.