Search

Background: Widely used psychotropic medications for obsessive-compulsive disorder (OCD) may change the volumes of subcortical brain structures, and differently in children vs. adults. We measured subcortical volumes cross-sectionally in patients finely stratified for age taking various common classes of OCD drugs. Methods: The ENIGMA-OCD consortium sample (1081 medicated/1159 unmedicated OCD patients and 2057 healthy controls aged 6–65) was divided into six successive 6–10-year age-groups. Individual structural MRIs were parcellated automatically using FreeSurfer into 8 regions-of-interest (ROIs). ROI volumes were compared between unmedicated and medicated patients and controls, and between patients taking serotonin reuptake inhibitors (SRIs), tricyclics (TCs), antipsychotics (APs), or benzodiazepines (BZs) and unmedicated patients. Results: Compared to unmedicated patients, volumes of accumbens, caudate, and/or putamen were lower in children aged 6–13 and adults aged 50–65 with OCD taking SRIs (Cohen's d = −0.24 to −0.74). Volumes of putamen, pallidum (d = 0.18–0.40), and ventricles (d = 0.31–0.66) were greater in patients aged 20–29 receiving APs. Hippocampal volumes were smaller in patients aged 20 and older taking TCs and/or BZs (d = −0.27 to −1.31). Conclusions: Results suggest that TCs and BZs could potentially aggravate hippocampal atrophy of normal aging in older adults with OCD, whereas SRIs may reduce striatal volumes in young children and older adults. Similar to patients with psychotic disorders, OCD patients aged 20–29 may experience subcortical nuclear and ventricular hypertrophy in relation to APs. Although cross-sectional, present results suggest that commonly prescribed agents exert macroscopic effects on subcortical nuclei of unknown relation to therapeutic response.

Objective: Recent work has shown that people with common epilepsies have characteristic patterns of cortical thinning, and that these changes may be progressive over time. Leveraging a large multicenter cross‐sectional cohort, we investigated whether regional morphometric changes occur in a sequential manner, and whether these changes in people with mesial temporal lobe epilepsy and hippocampal sclerosis (MTLE‐HS) correlate with clinical features. Methods: We extracted regional measures of cortical thickness, surface area, and subcortical brain volumes from T1‐weighted (T1W) magnetic resonance imaging (MRI) scans collected by the ENIGMA‐Epilepsy consortium, comprising 804 people with MTLE‐HS and 1625 healthy controls from 25 centers. Features with a moderate case–control effect size (Cohen d ≥ .5) were used to train an event‐based model (EBM), which estimates a sequence of disease‐specific biomarker changes from cross‐sectional data and assigns a biomarker‐based fine‐grained disease stage to individual patients. We tested for associations between EBM disease stage and duration of epilepsy, age at onset, and antiseizure medicine (ASM) resistance. Results: In MTLE‐HS, decrease in ipsilateral hippocampal volume along with increased asymmetry in hippocampal volume was followed by reduced thickness in neocortical regions, reduction in ipsilateral thalamus volume, and finally, increase in ipsilateral lateral ventricle volume. EBM stage was correlated with duration of illness (Spearman ρ = .293, p = 7.03 × 10−16), age at onset (ρ = −.18, p = 9.82 × 10−7), and ASM resistance (area under the curve = .59, p = .043, Mann–Whitney U test). However, associations were driven by cases assigned to EBM Stage 0, which represents MTLE‐HS with mild or nondetectable abnormality on T1W MRI. Significance: From cross‐sectional MRI, we reconstructed a disease progression model that highlights a sequence of MRI changes that aligns with previous longitudinal studies. This model could be used to stage MTLE‐HS subjects in other cohorts and help establish connections between imaging‐based progression staging and clinical features.

Deep neural networks show great promise for classifying brain diseases and making prognostic assessments based on neuroimaging data, but large, labeled training datasets are often required to achieve high predictive accuracy. Here we evaluated a range oftransfer learningor pre-training strategies to create useful MRI representations for downstream tasks that lack large amounts of training data, such as Alzheimer’s disease (AD) classification. To test our proposed pre-training strategies, we analyzed 4,098 3D T1-weighted brain MRI scans from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort and independently validated with an out-of-distribution test set of 600 scans from the Open Access Series of Imaging Studies (OASIS3) cohort for detecting AD. First, we trained 3D and 2D convolutional neural network (CNN) architectures. We tested combinations of multiple pre-training strategies based on (1) supervised, (2) contrastive learning, and (3) self-supervised learning - using pre-training data within versus outside the MRI domain. In our experiments, the 3D CNN pre-trained with contrastive learning provided the best overall results - when fine-tuned on T1-weighted scans for AD classification - outperformed the baseline by 2.8% when trained with all of the training data from ADNI. We also show test performance as a function of the training dataset size and the chosen pre-training method. Transfer learning offered significant benefits in low data regimes, with a performance boost of 7.7%. When the pre-trained model was used for AD classification, we were able to visualize an improved clustering of test subjects’ diagnostic groups, as illustrated via a uniform manifold approximation (UMAP) projection of the high-dimensional model embedding space. Further, saliency maps indicate the additional activation regions in the brain scan using pre-training, that then maximally contributed towards the final prediction score.

Emerging evidence suggests brain white matter alterations in adolescents with early-onset psychosis (EOP; age of onset d = 0.37), posterior corona radiata (d = 0.32), and superior fronto‐occipital fasciculus (d = 0.31). We also found widespread higher RD and more localized higher MD and AD. We detected significant effects of diagnostic subgroup, sex, and duration of illness, but not medication status. Using the largest EOP DTI sample to date, our findings suggest a profile of widespread white matter microstructure alterations in adolescents with EOP, most prominently in male individuals with early-onset schizophrenia and individuals with a shorter duration of illness.

This review compares the main brain abnormalities in schizophrenia (SZ), bipolar disorder (BD), major depressive disorder (MDD), and 22q11.2 Deletion Syndrome (22q11DS) determined by ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) consortium investigations. We obtained ranked effect sizes for subcortical volumes, regional cortical thickness, cortical surface area, and diffusion tensor imaging abnormalities, comparing each of these disorders relative to healthy controls. In addition, the studies report on significant associations between brain imaging metrics and disorder-related factors such as symptom severity and treatments. Visual comparison of effect size profiles shows that effect sizes are generally in the same direction and scale in severity with the disorders (in the order SZ > BD > MDD). The effect sizes for 22q11DS, a rare genetic syndrome that increases the risk for psychiatric disorders, appear to be much larger than for either of the complex psychiatric disorders. This is consistent with the idea of generally larger effects on the brain of rare compared to common genetic variants. Cortical thickness and surface area effect sizes for 22q11DS with psychosis compared to 22q11DS without psychosis are more similar to those of SZ and BD than those of MDD; a pattern not observed for subcortical brain structures and fractional anisotropy effect sizes. The observed similarities in effect size profiles for cortical measures across the psychiatric disorders mimic those observed for shared genetic variance between these disorders reported based on family and genetic studies and are consistent with shared genetic risk for SZ and BD and structural brain phenotypes.

Left–right asymmetry of the human brain is one of its cardinal features, and also a complex, multivariate trait. Decades of research have suggested that brain asymmetry may be altered in psychiatric disorders. However, findings have been inconsistent and often based on small sample sizes. There are also open questions surrounding which structures are asymmetrical on average in the healthy population, and how variability in brain asymmetry relates to basic biological variables such as age and sex. Over the last 4 years, the ENIGMA‐Laterality Working Group has published six studies of gray matter morphological asymmetry based on total sample sizes from roughly 3,500 to 17,000 individuals, which were between one and two orders of magnitude larger than those published in previous decades. A population‐level mapping of average asymmetry was achieved, including an intriguing fronto‐occipital gradient of cortical thickness asymmetry in healthy brains. ENIGMA's multi‐dataset approach also supported an empirical illustration of reproducibility of hemispheric differences across datasets. Effect sizes were estimated for gray matter asymmetry based on large, international, samples in relation to age, sex, handedness, and brain volume, as well as for three psychiatric disorders: autism spectrum disorder was associated with subtly reduced asymmetry of cortical thickness at regions spread widely over the cortex; pediatric obsessive–compulsive disorder was associated with altered subcortical asymmetry; major depressive disorder was not significantly associated with changes of asymmetry. Ongoing studies are examining brain asymmetry in other disorders. Moreover, a groundwork has been laid for possibly identifying shared genetic contributions to brain asymmetry and disorders., Left–right asymmetry of the human brain is one of its cardinal features, and also a complex, multivariate trait. Over the last four years, the ENIGMA‐Laterality Working Group has published six studies of grey matter morphological asymmetry in health and disease, based on total sample sizes from roughly 3,500 to 17,000 individuals, which were between one and two orders of magnitude larger than those published in previous decades. Here we review the findings from these six studies.

White matter tracts generated from whole brain tractography are often processed using automatic segmentation methods with standard atlases. Atlases are generated from hundreds of subjects, which becomes time-consuming to create and difficult to apply to all populations. In this study, we extended our prior work on using a deep generative model a Convolutional Variational Autoencoder - to map complex and data-intensive streamlines to a low-dimensional latent space given a limited sample size of 50 subjects from the ADNI3 dataset, to generate synthetic population-specific bundle templates using Kernel Density Estimation (KDE) on streamline embeddings. We conducted a quantitative shape analysis by calculating bundle shape metrics, and found that our bundle templates better capture the shape distribution of the bundles than the atlas data used in the original segmentation derived from young healthy adults. We further demonstrated the use of our framework for direct bundle segmentation from whole-brain tractograms.

Investigators in neuroscience have turned to Big Data to address replication and reliability issues by increasing sample sizes, statistical power, and representativeness of data. These efforts unveil new questions about integrating data arising from distinct sources and instruments. We focus on the most frequently assessed cognitive domain - memory testing - and demonstrate a process for reliable data harmonization across three common measures. We aggregated global raw data from 53 studies totaling N = 10,505 individuals. A mega-analysis was conducted using empirical bayes harmonization to remove site effects, followed by linear models adjusting for common covariates. A continuous item response theory (IRT) model estimated each individual’s latent verbal learning ability while accounting for item difficulties. Harmonization significantly reduced inter-site variance while preserving covariate effects, and our conversion tool is freely available online. This demonstrates that large-scale data sharing and harmonization initiatives can address reproducibility and integration challenges across the behavioral sciences.TeaserWe present a global effort to devise harmonization procedures necessary to meaningfully leverage big data.

Contains fulltext : 288062.pdf (Publisher’s version ) (Closed access) There is limited convergence in neuroimaging investigations into volumes of subcortical brain regions in social anxiety disorder (SAD). The inconsistent findings may arise from variations in methodological approaches across studies, including sample selection based on age and clinical characteristics. The ENIGMA-Anxiety Working Group initiated a global mega-analysis to determine whether differences in subcortical volumes can be detected in adults and adolescents with SAD relative to healthy controls. Volumetric data from 37 international samples with 1115 SAD patients and 2775 controls were obtained from ENIGMA-standardized protocols for image segmentation and quality assurance. Linear mixed-effects analyses were adjusted for comparisons across seven subcortical regions in each hemisphere using family-wise error (FWE)-correction. Mixed-effects d effect sizes were calculated. In the full sample, SAD patients showed smaller bilateral putamen volume than controls (left: d = -0.077, pFWE = 0.037; right: d = -0.104, pFWE = 0.001), and a significant interaction between SAD and age was found for the left putamen (r = -0.034, pFWE = 0.045). Smaller bilateral putamen volumes (left: d = -0.141, pFWE

Traumatic brain injury (TBI), a significant concern in military populations, is associated with alterations in brain structure and function, cognition, as well as physical and psychological dysfunction. Diffusion magnetic resonance imaging (dMRI) is particularly sensitive to changes in brain structure following TBI, as alterations in white matter (WM) microstructure are common. However, dMRI studies in mild TBI (mTBI) are conflicting, likely due to relatively small samples, sample heterogeneity (demographics, pre- and comorbidities) and injury characteristics (mechanism; chronicity). Furthermore, few studies account for brain asymmetry, which may impact cognitive functions subserved by WM tracts. Examining brain asymmetry in large samples may increase sensitivity to detect heterogeneous areas of subtle WM alteration in mTBI.Through the Enhancing Neuroimaging and Genetics through Meta-analysis (ENIGMA) Military-Relevant Brain Injury working group, we conducted a mega-analysis of neuroimaging and clinical data from 16 cohorts of Active Duty Service Members and Veterans (n=2,598; 2,321 males/277 females; age 19-85 years). 1,080 reported a deployment-related TBI, 480 had a history of only non-military-related TBI, 823 reported no history of TBI, and 215 did not differentiate between military and non-military TBI. dMRI data were processed in a harmonized manner along with harmonized demographic, injury, psychiatric, and cognitive measures. Hemispheric asymmetry of fractional anisotropy (FA, a common proxy for myelin organization) was calculated for 19 WM tracts and compared between those with and without TBI history.FA in the cingulum showed greater asymmetry in individuals with a history of deployment-related TBI; this effect was driven by greater left lateralization in the group with TBI. There was a trend towards lower FA of the right cingulum in the TBI group. These results remained significant after accounting for potentially confounding variables including posttraumatic stress disorder, depression, and handedness and were driven primarily by individuals who had sustained their worst TBI before age 40. We further found that alterations in the cingulum were associated with slower processing speed and poorer set shifting.The results indicate an enhancement of the previously reported natural left laterality, possibly due to vulnerability of the non-dominant hemisphere or compensatory mechanisms in the dominant hemisphere. The cingulum is one of the last WM tracts to mature, reaching peak FA around 42 years old. This effect was primarily detected in individuals whose worst injury occurred before age 40, suggesting that the protracted development of the cingulum may lead to increased vulnerability to insults, such as TBI.

BackgroundCurrent clinical assessments of Posttraumatic stress disorder (PTSD) rely solely on subjective symptoms and experiences reported by the patient, rather than objective biomarkers of the illness. Recent advances in data-driven computational approaches have been helpful in devising tools to objectively diagnose psychiatric disorders. Here we aimed to classify individuals with PTSD versus controls using heterogeneous brain datasets from the ENIGMA-PGC PTSD Working group.MethodsWe analyzed brain MRI data from 3,527 structural-MRI; 2,502 resting state-fMRI; and 1,953 diffusion-MRI. First, we identified the brain features that best distinguish individuals with PTSD from controls (TEHC and HC) using traditional machine learning methods. Second, we assessed the utility of the denoising variational autoencoder (DVAE) and evaluated its classification performance. Third, we assessed the generalizability and reproducibility of both models using leave-one-site-out cross-validation procedure for each modality.ResultsWe found lower performance in classifying PTSD vs. controls with data from over 20 sites (60% test AUC for s-MRI, 59% for rs-fMRI and 56% for d-MRI), as compared to other studies run on single-site data. The performance increased when classifying PTSD from HC without trauma history across all three modalities (75% AUC). The classification performance remained intact when applying the DVAE framework, which reduced the number of features. Finally, we found that the DVAE framework achieved better generalization to unseen datasets compared with the traditional machine learning frameworks, albeit performance was slightly above chance.ConclusionOur findings highlight the promise offered by machine learning methods for the diagnosis of patients with PTSD. The utility of brain biomarkers across three MRI modalities and the contribution of DVAE models for improving generalizability offers new insights into neural mechanisms involved in PTSD.Significance⍰Classifying PTSD from trauma-unexposed healthy controls (HC) using three imaging modalities performed well (∼75% AUC), but performance suffered markedly when classifying PTSD from trauma-exposed healthy controls (TEHC) using three imaging modalities (∼60% AUC).⍰Using deep learning for feature reduction (denoising variational auto-encoder; DVAE) dramatically reduced the number of features with no concomitant performance degradation.⍰Utilizing denoising variational autoencoder (DVAE) models improves generalizability across heterogeneous multi-site data compared with the traditional machine learning frameworks

Schizotypy represents an index of psychosis-proneness in the general population often associated with childhood trauma exposure. Both schizotypy and childhood trauma are linked to structural brain alterations, and it is possible that trauma exposure moderates the extent of brain morphological differences associated with schizotypy. We addressed this question using data from a total of 1,182 healthy adults (age range: 18-65 years old, 647 females/535 males), pooled from nine sites worldwide, contributing to the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Schizotypy working group. All participants completed both the Schizotypal Personality Questionnaire Brief version (SPQ-B), and the Childhood Trauma Questionnaire (CTQ), and underwent a 3D T1-weighted brain MRI scan from which regional indices of subcortical grey matter volume and cortical thickness were determined. A series of multiple linear regressions revealed that differences in cortical thickness in four regions-of-interest were significantly associated with interactions between schizotypy and trauma; subsequent moderation analyses indicated that increasing levels of schizotypy were associated with thicker left caudal anterior cingulate gyrus, right middle temporal gyrus and insula, and thinner left caudal middle frontal gyrus, in people exposed to higher (but not low or average) levels of childhood trauma. This was found in the context of thicker bilateral medial orbitofrontal gyri, right rostral anterior cingulate gyrus, left temporal pole, left insula, and thinner left paracentral lobule directly associated with increasing levels of schizotypy. In addition, thinner left postcentral, superior parietal and lingual gyri, as well as thicker left caudal middle frontal gyrus and smaller left thalamus and right caudate were associated with increasing levels of childhood trauma exposure. These results suggest that alterations in brain regions critical for higher cognitive and integrative processes that are associated with schizotypy may be enhanced in individuals exposed to high levels of trauma.

Neuroimaging studies point to neurostructural abnormalities in youth with anxiety disorders. Yet, findings are based on small-scale studies, often with small effect sizes, and have limited generalizability and clinical relevance. These issues have prompted a paradigm shift in the field towards highly powered (i.e., big data) individual-level inferences, which are data-driven, transdiagnostic, and neurobiologically informed. Here, we built and validated neurostructural machine learning (ML) models for individual-level inferences based on the largest-ever multi-site neuroimaging sample of youth with anxiety disorders (age: 10-25 years, N=3,343 individuals from 32 global sites), as compiled by three ENIGMA Anxiety Working Groups: Panic Disorder (PD), Generalized Anxiety Disorder (GAD), and Social Anxiety Disorder (SAD). ML classifiers were trained on MRI-derived regional measures of cortical thickness, surface area, and subcortical volumes to classify patients and healthy controls (HC) for each anxiety disorder separately and across disorders (transdiagnostic classification). Modest, yet robust, classification performance was achieved for PD vs. HC (AUC=0.62), but other disorder-specific and transdiagnostic classifications were not significantly different from chance. However, above chance-level transdiagnostic classifications were obtained in exploratory subgroup analyses of male patients vs. male HC, unmedicated patients vs. HC, and patients with low anxiety severity vs. HC (AUC 0.59-0.63). The above chance-level classifications were based on plausible and specific neuroanatomical features in fronto-striato-limbic and temporo-parietal regions. This study provides a realistic estimate of classification performance in a large, ecologically valid, multi-site sample of youth with anxiety disorders, and may as such serve as a benchmark.

Deep learning methods trained on brain MRI data from one scanner or imaging protocol can fail catastrophically when tested on data from other sites or protocols - a problem known asdomain shift. To address this, here we propose adomain adaptationmethod that trains a 3D CycleGAN (cycle-consistent generative adversarial network) to harmonize brain MRI data from 5 diverse sources (ADNI, WHIMS, OASIS, AIBL, and UK Biobank; total N=4,941 MRIs, age range: 46-96 years). The approach uses 2 generators and 2 discriminators to generate an image harmonized to a specific target dataset given an image from the source domain distribution andvice versa. We train the CycleGAN to jointly optimize an adversarial loss and cyclic consistency. We use a patch-based discriminator and impose identity loss to further regularize model training. To test the benefit of the harmonization, we show that brain age estimation - a common benchmarking task - is more accurate in GAN-harmonized versus raw data.t-SNE maps show the improved distributional overlap of the harmonized data in the latent space.

Current Alzheimer’s disease (AD) research has a major focus on validating and discovering noninvasive biomarkers that can detect AD, benchmark disease severity, and aid in testing the efficacy of interventions. Structural magnetic resonance imaging (sMRI) is a well-validated tool used in diagnosis and for monitoring disease progression in AD. Much of the sMRI literature centers around hippocampal and other medial temporal lobe structure atrophy, which are strongly associated with cognition and diagnosis. Because atrophy patterns are complex and vary by patient, researchers have made efforts to condense more brain information into validated metrics. Many of these methods use machine learning (ML), which can be difficult to interpret clinically, hampering clinical adoption. Here, we introduce a practical, clinically meaningful and interpretable index which we call an “AD-NeuroScore.” Our approach is automated and uses multiple regional brain volumes associated with cognitive decline. We used a modified Euclidean inspired distance function to calculate the differences between each participant and a cognitively normal (CN) older adult template, adjusting for intracranial volume, age, sex, and scanner model. Here we report validation results, including sensitivity to diagnosis (CN, mild cognitive impairment (MCI), and AD) and disease severity (Clinical Dementia Rating Scale Sum of Boxes (CDR-SB), Mini Mental State Exam (MMSE), and Alzheimer’s Disease Assessment Scale-Cognitive Subscale (ADAS-11) in 929 older adults (mean age=72.7 years, SD=6.3, Range=55.1-91.5, 50% Female) drawn from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) study. To determine if AD-NeuroScore might be predictive of disease progression, we assessed the relationship between the calculated AD-NeuroScore at baseline and change in both diagnosis and disease severity scores at 12, 24, 36, and 48-months. We performed additional validation in all analyses, benchmarking AD-NeuroScore against adjusted hippocampal volume (AHV). We found that AD-NeuroScore was significantly associated with diagnosis and all disease severity scores at baseline. Associations between AD-NeuroScore and disease severity (CDR-SB and ADAS-11) were significantly stronger than with AHV. Baseline AD-NeuroScore was also associated with change in diagnosis and changes in disease severity scores at all time points. Performance was equivalent, or in some cases superior, to AHV. These early validation results suggest that AD-NeuroScore has the potential to be a clinically meaningful biomarker for dementia.

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants’ demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LME INT), (2) LME that models both site-specific random intercepts and age-related random slopes (LME INT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2–81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3–85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ 2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ 2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ 2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects.

Automatic neuroimaging segmentation and parcellation tools provide convenient and systematic methods for extracting numerous features from brain MRI scans, and are becoming standard practice for large-scale coordinated studies. One such tool, FreeSurfer, provides an easy-to-use pipeline to extract metrics describing cortical and subcortical morphometry. Over the past two decades, there have been over 25 stable releases of FreeSurfer, and different versions are used across published works. Despite this, the reliability and compatibility of metrics derived from the most recent major version releases have yet to be assessed empirically. Here, we use test-retest data from three public brain MRI datasets to assess within-version reliability and between-version compatibility across 42 regional outputs from three versions of FreeSurfer: the latest, v7.1, and two previous stable releases - v5.3, and v6.0. We find v7.1 was less compatible with older versions for measuring cortical thickness. In particular, the thickness of the cingulate gyrus had low compatibility (intraclass correlation coefficient (ICC) between 0.37 and 0.61) between versions. Temporal and frontal poles, and the medial orbitofrontal surface area metrics, also showed low to moderate compatibility with v7.1. While our work compares all three versions, our sub-comparisons between the older versions (v5.3 and v6.0) replicates earlier findings of low compatibility of pallidum and putamen volumes. Low between-version compatibility was not always indicative of low within-version reliability – all versions showed good to excellent reliability across most regional measures (ICC>0.8). Age associations, quality control metrics, and Dice coefficients in an independent sample of 106 individual scans, processed with all three versions of FreeSurfer, revealed differences in results of downstream statistical analysis. As neuroimaging studies adopt more recently released software, we provide researchers with a reference to highlight the regions and metrics that may yield findings inconsistent with published works using older FreeSurfer software. An interactive viewer for the results is provided at http://data.brainescience.org/Freesurfer_Reliability/

BackgroundThe cerebellum critically contributes to higher-order cognitive and emotional functions such fear learning and memory. Prior research on cerebellar volume in PTSD is scant and has neglected neuroanatomical subdivisions of the cerebellum that differentially map on to motor, cognitive, and affective functions.MethodsWe quantified cerebellar lobule volumes using structural magnetic resonance imaging in 4,215 adults (PTSD n= 1640; Control n=2575) across 40 sites from the from the ENIGMA-PGC PTSD working group. Using a new state-of-the-art deep-learning based approach for automatic cerebellar parcellation, we obtained volumetric estimates for the total cerebellum and 28 subregions. Linear mixed effects models controlling for age, gender, intracranial volume, and site were used to compare cerebellum total and subregional volume in PTSD compared to healthy controls. The Benjamini-Hochberg procedure was used to control the false discovery rate (p-FDR< .05).ResultsPTSD was associated with significant grey and white matter reductions of the cerebellum. Compared to controls, people with PTSD demonstrated smaller total cerebellum volume. In addition, people with PTSD showed reduced volume in subregions primarily within the posterior lobe (lobule VIIB, crus II), but also the vermis (VI, VIII), flocculonodular lobe (lobule X), and cerebellar white matter (allp-FDR< 0.05). Effects of PTSD on volume were consistent, and generally more robust, when examining symptom severity rather than diagnostic status.ConclusionsThese findings implicate regionally specific cerebellar volumetric differences in the pathophysiology of PTSD. The cerebellum appears to play an important role in high-order cognitive and emotional processes, far beyond its historical association with vestibulomotor function. Further examination of the cerebellum in trauma-related psychopathology will help to clarify how cerebellar structure and function may disrupt cognitive and affective processes at the center of translational models for PTSD.

Movement disorders 38(1), 45-56 (2023). doi:10.1002/mds.29261, Published by Wiley, New York, NY

Machine learning methods have been used for over a decade for staging and subtyping a variety of brain diseases, offering fast and objective methods to classify neurodegenerative diseases such as Alzheimer’s disease (AD). Deep learning models based on convolutional neural networks (CNNs) have also been used to infer dementia severity and predict future clinical decline. Most CNN-based deep learning models use T1-weighted brain MRI scans to identify predictive features for these tasks. In contrast, we examine the added value of diffusion-weighted MRI (dMRI) - a variant of MRI, sensitive to microstructural tissue properties - as an additional input in CNN-based models of dementia severity. dMRI is sensitive to microstructural brain abnormalities not evident on standard anatomical MRI. By training CNNs on combined anatomical and diffusion MRI, we hypothesize that we could boost performance when predicting widely-used clinical assessments of dementia severity, such as individuals’ scores on the ADAS11, ADAS13, and MMSE (mini-mental state exam) clinical scales. For benchmarking, we evaluate CNNs that use T1-weighted MRI and dMRI to estimate “brain age” - the task of predicting a person’s chronological age from their neuroimaging data. To assess which dMRI-derived maps were most beneficial, we computed DWI-derived diffusion tensor imaging (DTI) maps of mean and radial diffusivity (MD/RD), axial diffusivity (AD) and fractional anisotropy (FA) for 1198 elderly subjects (age: 74.35 +/- 7.74 yrs.; 600 F/598 M, with a distribution of 636 CN/421 MCI/141 AD) from the Alzheimer’s Disease Neuroimaging Initiative (ADNI). We tested both 2D Slice CNN and 3D CNN neural network models for the above predictive tasks. Our results suggest that for at least some deep learning architectures, diffusion-weighted MRI may enhance performance for several AD-relevant deep learning tasks relative to using T1-weighted images alone.

Current knowledge about functional connectivity in obsessive-compulsive disorder (OCD) is based on small-scale studies, limiting the generalizability of results. Moreover, the majority of studies have focused only on predefined regions or functional networks rather than connectivity throughout the entire brain. Here, we investigated differences in resting-state functional connectivity between OCD patients and healthy controls (HC) using mega-analysis of data from 1,024 OCD patients and 1,028 HC from 28 independent samples of the ENIGMA-OCD consortium. We assessed group differences in whole-brain functional connectivity at both the regional and network level, and investigated whether functional connectivity could serve as biomarker to identify patient status at the individual level using machine learning analysis. The mega-analyses revealed widespread abnormalities in functional connectivity in OCD, with global hypo-connectivity (Cohen’s d: -0.27 to -0.13) and few hyper-connections, mainly with the thalamus (Cohen’s d: 0.19 to 0.22). Most hypo-connections were located within the sensorimotor network and no fronto-striatal abnormalities were found. Overall, classification performances were poor, with area-under-the-receiver-operating-characteristic curve (AUC) scores ranging between 0.567 and 0.673, with better classification for medicated (AUC=0.702) than unmedicated (AUC=0.608) patients versus healthy controls. These findings provide partial support for existing pathophysiological models of OCD and highlight the important role of the sensorimotor network in OCD. However, resting-state connectivity does not so far provide an accurate biomarker for identifying patients at the individual level.

Whole brain tractography is commonly used to study the brain’s white matter fiber pathways, but the large number of streamlines generated - up to one million per brain - can be challenging for large-scale population studies. We propose a robust dimensionality reduction framework for tractography, using a Convolutional Variational Autoencoder (ConvVAE) to learn low-dimensional embeddings from white matter bundles. The resulting embeddings can be used to facilitate downstream tasks such as outlier and abnormality detection, and mapping of disease effects on white matter tracts in individuals or groups. We design experiments to evaluate how well embeddings of different dimensions preserve distances from the original high-dimensional dataset, using distance correlation methods. We find that streamline distances and inter-bundle distances are well preserved in the latent space, with a 6-dimensional optimal embedding space. The generative ConvVAE model allows fast inference on new data, and the smooth latent space enables meaningful decodings that can be used for downstream tasks. We demonstrate the use of a ConvVAE model trained on control subjects’ data to detect structural anomalies in white matter tracts in patients with Alzheimer’s disease (AD). Using ConvVAEs to facilitate population analyses, we identified 6 tracts with statistically significant differences between AD and controls after controlling for age and sex effect, visualizing specific locations along the tracts with high anomalies despite large inter-subject variations in fiber bundle geometry.

Whole-brain tractograms generated from diffusion MRI digitally represent the white matter structure of the brain and are composed of millions of streamlines. Such tractograms can have false positive and anatomically implausible streamlines. To obtain anatomically relevant streamlines and tracts, supervised and unsupervised methods can be used for tractogram clustering and tract extraction. Here we propose FiberNeat, an unsupervised white matter tract filtering method. FiberNeat takes an input set of streamlines that could either be unlabeled clusters or labeled tracts. Individual clusters/tracts are projected into a latent space using nonlinear dimensionality reduction techniques, t-SNE and UMAP, to find spurious and outlier streamlines. In addition, outlier streamline clusters are detected using DBSCAN and then removed from the data in streamline space. We performed quantitative comparisons with expertly delineated tracts. We ran FiberNeat on 131 participants' data from the ADNI3 dataset. We show that applying FiberNeat as a filtering step after bundle segmentation improves the quality of extracted tracts and helps improve tractometry.

Large data initiatives and high-powered brain imaging analyses require the pooling of MR images acquired across multiple scanners, often using different protocols. Prospective cross-site harmonization often involves the use of a phantom or traveling subjects. However, as more datasets are becoming publicly available, there is a growing need for retrospective harmonization, pooling data from sites not originally coordinated together. Several retrospective harmonization techniques have shown promise in removing cross-site image variation. However, most unsupervised methods cannot distinguish between image-acquisition based variability and cross-site population variability, so they require that datasets contain subjects or patient groups with similar clinical or demographic information. To overcome this limitation, we consider cross-site MRI image harmonization as a style transfer problem rather than a domain transfer problem. Using a fully unsupervised deep-learning framework based on a generative adversarial network (GAN), we show that MR images can be harmonized by inserting the style information encoded from a reference image directly, without knowing their site/scanner labelsa priori. We trained our model using data from five large-scale multi-site datasets with varied demographics. Results demonstrated that our styleencoding model can harmonize MR images, and match intensity profiles, successfully, without relying on traveling subjects. This model also avoids the need to control for clinical, diagnostic, or demographic information. Moreover, we further demonstrated that if we included diverse enough images into the training set, our method successfully harmonized MR images collected from unseen scanners and protocols, suggesting a promising novel tool for ongoing collaborative studies.

Sex is a biological variable that contributes to individual variability in brain structure and behavior. Neuroimaging studies of population‐based samples have identified normative differences in brain structure between males and females, many of which are exacerbated in psychiatric and neurological conditions. Still, sex differences in MRI outcomes are understudied, particularly in clinical samples with known sex differences in disease risk, prevalence, and expression of clinical symptoms. Here we review the existing literature on sex differences in adult brain structure in normative samples and in 14 distinct psychiatric and neurological disorders. We discuss commonalities and sources of variance in study designs, analysis procedures, disease subtype effects, and the impact of these factors on MRI interpretation. Lastly, we identify key problems in the neuroimaging literature on sex differences and offer potential recommendations to address current barriers and optimize rigor and reproducibility. In particular, we emphasize the importance of large‐scale neuroimaging initiatives such as the Enhancing NeuroImaging Genetics through Meta‐Analyses consortium, the UK Biobank, Human Connectome Project, and others to provide unprecedented power to evaluate sex‐specific phenotypes in major brain diseases., Here, we review the human neuroimaging literature examining sex effects on adult brain structure using structural and diffusion MRI. We discuss normative sex differences based on population‐based studies as well as sex differences in 14 major brain diseases. Finally, we identify key barriers to advancing the science on neuroimaging sex effects and offer recommendations to mitigate these challenges, particularly through large‐scale neuroimaging.

BackgroundMajor depressive disorder (MDD) is associated with an increased risk of brain atrophy, aging-related diseases, and mortality. We examined potential advanced brain aging in MDD patients, and whether this process is associated with clinical characteristics in a large multi-center international dataset.MethodsWe performed a mega-analysis by pooling brain measures derived from T1-weighted MRI scans from 29 samples worldwide. Normative brain aging was estimated by predicting chronological age (10-75 years) from 7 subcortical volumes, 34 cortical thickness and 34 surface area, lateral ventricles and total intracranial volume measures separately in 1,147 male and 1,386 female controls from the ENIGMA MDD working group. The learned model parameters were applied to 1,089 male controls and 1,167 depressed males, and 1,326 female controls and 2,044 depressed females to obtain independent unbiased brain-based age predictions. The difference between predicted “brain age” and chronological age was calculated to indicate brain predicted age difference (brain-PAD).FindingsOn average, MDD patients showed a higher brain-PAD of +0.90 (SE 0.21) years (Cohen’s d=0.12, 95% CI 0.06-0.17) compared to controls. Relative to controls, first-episode and currently depressed patients showed higher brain-PAD (+1.2 [0.3] years), and the largest effect was observed in those with late-onset depression (+1.7 [0.7] years). In addition, higher brain-PAD was associated with higher self-reported depressive symptomatology (b=0.05, p=0.004).InterpretationThis highly powered collaborative effort showed subtle patterns of abnormal structural brain aging in MDD. Substantial within-group variance and overlap between groups were observed. Longitudinal studies of MDD and somatic health outcomes are needed to further assess the predictive value of these brain-PAD estimates.FundingThis work was supported, in part, by NIH grants U54 EB020403 and R01 MH116147.

Studies of posttraumatic stress disorder (PTSD) report volume abnormalities in multiple regions of the cerebral cortex. However, findings for many regions, particularly regions outside commonly studied emotion-related prefrontal, insular, and limbic regions, are inconsistent and tentative. Also, few studies address the possibility that PTSD abnormalities may be confounded by comorbid depression. A mega-analysis investigating all cortical regions in a large sample of PTSD and control subjects can potentially provide new insight into these issues. Given this perspective, our group aggregated regional volumes data of 68 cortical regions across both hemispheres from 1379 PTSD patients to 2192 controls without PTSD after data were processed by 32 international laboratories using ENIGMA standardized procedures. We examined whether regional cortical volumes were different in PTSD vs. controls, were associated with posttraumatic stress symptom (PTSS) severity, or were affected by comorbid depression. Volumes of left and right lateral orbitofrontal gyri (LOFG), left superior temporal gyrus, and right insular, lingual and superior parietal gyri were significantly smaller, on average, in PTSD patients than controls (standardized coefficients = −0.111 to −0.068, FDR corrected P values < 0.039) and were significantly negatively correlated with PTSS severity. After adjusting for depression symptoms, the PTSD findings in left and right LOFG remained significant. These findings indicate that cortical volumes in PTSD patients are smaller in prefrontal regulatory regions, as well as in broader emotion and sensory processing cortical regions.

ObjectiveSpinal cord damage is a hallmark of Friedreich ataxia (FRDA), but its progression and clinical correlates remain unclear. Here we performed a characterization of cervical spinal cord structural abnormalities in a large multisite FRDA cohort.MethodsWe performed a cross-sectional analysis of cervical spinal cord (C1 to C4) cross-sectional area (CSA) and eccentricity using MRI data from eight sites within the ENIGMA-Ataxia initiative, including 256 individuals with FRDA and 223 age- and sex-matched controls. Correlations and subgroup analyses within the FRDA cohort were undertaken based on disease duration, ataxia severity, and onset age.ResultsIndividuals with FRDA, relative to controls, had significantly reduced CSA at all examined levels, with large effect sizes (d>2.1) and significant correlations with disease severity (rd>1.2), but without significant clinical correlations. Subgroup analyses showed that CSA and eccentricity are abnormal at all disease stages. However, while CSA appears to decrease progressively, eccentricity remains stable over time.InterpretationPrevious research has shown that increased eccentricity reflects dorsal column (DC) damage, while decreased CSA reflects either DC or corticospinal tract (CST) damage or both. Hence, our data support the hypothesis that damage to DC and CST follow distinct courses in FRDA: developmental abnormalities likely define the DC, whereas CST alterations may be both developmental and degenerative. These results provide new insights about FRDA pathogenesis and indicate that CSA of the cervical spinal cord should be investigated further as a potential biomarker of disease progression.

Mechanisms underpinning age-related variations in cortical thickness in the human brain remain poorly understood. We investigated whether inter-regional age-related variations in cortical thinning (in a multicohort neuroimaging dataset from the ENIGMA Lifespan Working Group totalling 14,248 individuals, aged 4-89 years) depended on cell-specific marker gene expression levels. We found differences amidst early-life (60 years) in the patterns of association between inter-regional profiles of cortical thickness and expression profiles of marker genes for CA1 and S1 pyramidal cells, astrocytes, and microglia. Gene ontology and enrichment analyses indicated that each of the three life-stages was associated with different biological processes and cellular components: synaptic modeling in early life, neurotransmission in mid-life, and neurodegeneration in late-life. These findings provide mechanistic insights into age-related cortical thinning during typical development and aging.

Traumatic brain injury (TBI) is a major cause of death and disability in children in both developed and developing nations. Children and adolescents suffer from TBI at a higher rate than the general population, and specific developmental issues require a unique context since findings from adult research do not necessarily directly translate to children. Findings in pediatric cohorts tend to lag behind those in adult samples. This may be due, in part, both to the smaller number of investigators engaged in research with this population and may also be related to changes in safety laws and clinical practice that have altered length of hospital stays, treatment, and access to this population. The ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) Pediatric Moderate/Severe TBI (msTBI) group aims to advance research in this area through global collaborative meta-analysis of neuroimaging data. In this paper, we discuss important challenges in pediatric TBI research and opportunities that we believe the ENIGMA Pediatric msTBI group can provide to address them. With the paucity of research studies examining neuroimaging biomarkers in pediatric patients with TBI and the challenges of recruiting large numbers of participants, collaborating to improve statistical power and to address technical challenges like lesions will significantly advance the field. We conclude with recommendations for future research in this field of study.

BackgroundEmerging evidence suggests brain white matter alterations in adolescents with early-onset psychosis (EOP; age of onset MethodsOur sample included 321 adolescents with EOP (mean age: 16.3 ± 1.4 years, 46.4% females) and 265 adolescent CTR (mean age: 16.0 ± 1.7 years, 57.7% females) pooled from nine sites. All sites extracted mean fractional anisotropy (FA), mean diffusivity (MD), radial diffusivity (RD), and axial diffusivity (AD) for 25 white matter regions of interest per participant. ComBat harmonization was performed for all DTI measures to adjust for scanner differences. Multiple linear regression models were fitted to investigate case-control differences and associations with clinical variables in regional DTI measures.ResultsWe found widespread lower FA in EOP compared to CTR, with the largest effect sizes in the superior longitudinal fasciculus (Cohen’s d = 0.37), posterior corona radiata (d = 0.32), and superior fronto□occipital fasciculus (d = 0.31). We also found widespread higher RD and more localized higher MD and AD. We detected significant effects of diagnostic subgroup, sex, and duration of illness, but not medication status.ConclusionUsing the largest EOP DTI sample to date, our findings suggest a profile of widespread white matter microstructure alterations in adolescents with EOP, most prominently in male patients with early-onset schizophrenia and patients with a shorter duration of illness.

Numerous brain disorders demonstrate structural brain abnormalities, which are thought to arise from molecular perturbations or connectome miswiring. The unique and shared contributions of these molecular and connectomic vulnerabilities to brain disorders remain unknown, and has yet to be studied in a single multi-disorder framework. Using MRI morphometry from the ENIGMA consortium, we construct maps of cortical abnormalities for thirteen neurodevelopmental, neurological, and psychiatric disorders from N = 21 000 patients and N = 26 000 controls, collected using a harmonized processing protocol. We systematically compare cortical maps to multiple micro-architectural measures, including gene expression, neurotransmitter density, metabolism, and myelination (molecular vulnerability), as well as global connectomic measures including number of connections, centrality, and connection diversity (connectomic vulnerability). We find that regional molecular vulnerability and macroscale brain network architecture interact to drive the spatial patterning of cortical abnormalities in multiple disorders. Local attributes, particularly neurotransmitter receptor profiles, constitute the best predictors of both disorder-specific cortical morphology and cross-disorder similarity. Finally, we find that cross-disorder abnormalities are consistently subtended by a small subset of network epicentres in bilateral sensory-motor, medial temporal lobe, precuneus, and superior parietal cortex. Collectively, our results highlight how local biological attributes and global connectivity jointly shape cross-disorder cortical abnormalities.

Numerous brain disorders demonstrate structural brain abnormalities, which are thought to arise from molecular perturbations or connectome miswiring. The unique and shared contributions of these molecular and connectomic vulnerabilities to brain disorders remain unknown, and has yet to be studied in a single multi-disorder framework. Using MRI morphometry from the ENIGMA consortium, we construct maps of cortical abnormalities for thirteen neurodevelopmental, neurological, and psychiatric disorders from N = 21,000 participants and N = 26,000 controls, collected using a harmonised processing protocol. We systematically compare cortical maps to multiple micro-architectural measures, including gene expression, neurotransmitter density, metabolism, and myelination (molecular vulnerability), as well as global connectomic measures including number of connections, centrality, and connection diversity (connectomic vulnerability). We find a relationship between molecular vulnerability and white-matter architecture that drives cortical disorder profiles. Local attributes, particularly neurotransmitter receptor profiles, constitute the best predictors of both disorder-specific cortical morphology and cross-disorder similarity. Finally, we find that cross-disorder abnormalities are consistently subtended by a small subset of network epicentres in bilateral sensory-motor, inferior temporal lobe, precuneus, and superior parietal cortex. Collectively, our results highlight how local molecular attributes and global connectivity jointly shape cross-disorder cortical abnormalities.

Schizophrenia (SZ) is associated with an increased risk of life-long cognitive impairments, age-related chronic disease, and premature mortality. We investigated evidence for advanced brain ageing in adult SZ patients, and whether this was associated with clinical characteristics in a prospective meta-analytic study conducted by the ENIGMA Schizophrenia Working Group. The study included data from 26 cohorts worldwide, with a total of 2803 SZ patients (mean age 34.2 years; range 18-72 years; 67% male) and 2598 healthy controls (mean age 33.8 years, range 18-73 years, 55% male). Brain-predicted age was individually estimated using a model trained on independent data based on 68 measures of cortical thickness and surface area, 7 subcortical volumes, lateral ventricular volumes and total intracranial volume, all derived from T1-weighted brain magnetic resonance imaging (MRI) scans. Deviations from a healthy brain ageing trajectory were assessed by the difference between brain-predicted age and chronological age (brain-predicted age difference [brain-PAD]). On average, SZ patients showed a higher brain-PAD of +3.64 years (95% CI: 3.01, 4.26; I2 = 55.28%) compared to controls, after adjusting for age and sex (Cohen’s d = 0.50). Among SZ patients, brain-PAD was not associated with specific clinical characteristics (age of onset, duration of illness, symptom severity, or antipsychotic use and dose). This large-scale collaborative study suggests advanced structural brain ageing in SZ. Longitudinal studies of SZ and a range of mental and somatic health outcomes will help to further evaluate the clinical implications of increased brain-PAD and its ability to be influenced by interventions.

Contains fulltext : 281801.pdf (Publisher’s version ) (Open Access) It is increasingly recognized that multiple psychiatric conditions are underpinned by shared neural pathways, affecting similar brain systems. Here, we carried out a multiscale neural contextualization of shared alterations of cortical morphology across six major psychiatric conditions (autism spectrum disorder, attention deficit/hyperactivity disorder, major depression disorder, obsessive-compulsive disorder, bipolar disorder, and schizophrenia). Our framework cross-referenced shared morphological anomalies with respect to cortical myeloarchitecture and cytoarchitecture, as well as connectome and neurotransmitter organization. Pooling disease-related effects on MRI-based cortical thickness measures across six ENIGMA working groups, including a total of 28,546 participants (12,876 patients and 15,670 controls), we identified a cortex-wide dimension of morphological changes that described a sensory-fugal pattern, with paralimbic regions showing the most consistent alterations across conditions. The shared disease dimension was closely related to cortical gradients of microstructure as well as neurotransmitter axes, specifically cortex-wide variations in serotonin and dopamine. Multiple sensitivity analyses confirmed robustness with respect to slight variations in analytical choices. Our findings embed shared effects of common psychiatric conditions on brain structure in multiple scales of brain organization, and may provide insights into neural mechanisms of transdiagnostic vulnerability.