Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain.

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth.

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Large-scale analysis of structural brain asymmetries in schizophrenia via the ENIGMA consortium.

Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

Genomic loci influence patterns of structural covariance in the human brain.

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

neuromaps: structural and functional interpretation of brain maps.

Imaging technologies are increasingly used to generate high-resolution reference maps of brain structure and function. Comparing experimentally generated maps to these reference maps facilitates cross-disciplinary scientific discovery. Although recent data sharing initiatives increase the accessibility of brain maps, data are often shared in disparate coordinate systems, precluding systematic and accurate comparisons. Here we introduce neuromaps, a toolbox for accessing, transforming and analyzing structural and functional brain annotations. We implement functionalities for generating high-quality transformations between four standard coordinate systems. The toolbox includes curated reference maps and biological ontologies of the human brain, such as molecular, microstructural, electrophysiological, developmental and functional ontologies. Robust quantitative assessment of map-to-map similarity is enabled via a suite of spatial autocorrelation-preserving null models. neuromaps combines open-access data with transparent functionality for standardizing and comparing brain maps, providing a systematic workflow for comprehensive structural and functional annotation enrichment analysis of the human brain.

ASLPrep: a platform for processing of arterial spin labeled MRI and quantification of regional brain perfusion.

Arterial spin labeled (ASL) magnetic resonance imaging (MRI) is the primary method for noninvasively measuring regional brain perfusion in humans. We introduce ASLPrep, a suite of software pipelines that ensure the reproducible and generalizable processing of ASL MRI data.

Connectome architecture shapes large-scale cortical alterations in schizophrenia: a worldwide ENIGMA study.

Schizophrenia is a prototypical network disorder with widespread brain-morphological alterations, yet it remains unclear whether these distributed alterations robustly reflect the underlying network layout. We tested whether large-scale structural alterations in schizophrenia relate to normative structural and functional connectome architecture, and systematically evaluated robustness and generalizability of these network-level alterations. Leveraging anatomical MRI scans from 2439 adults with schizophrenia and 2867 healthy controls from 26 ENIGMA sites and normative data from the Human Connectome Project (n = 207), we evaluated structural alterations of schizophrenia against two network susceptibility models: (i) hub vulnerability, which examines associations between regional network centrality and magnitude of disease-related alterations; (ii) epicenter mapping, which identifies regions whose typical connectivity profile most closely resembles the disease-related morphological alterations. To assess generalizability and specificity, we contextualized the influence of site, disease stages, and individual clinical factors and compared network associations of schizophrenia with that found in affective disorders. Our findings show schizophrenia-related cortical thinning is spatially associated with functional and structural hubs, suggesting that highly interconnected regions are more vulnerable to morphological alterations. Predominantly temporo-paralimbic and frontal regions emerged as epicenters with connectivity profiles linked to schizophrenia's alteration patterns. Findings were robust across sites, disease stages, and related to individual symptoms. Moreover, transdiagnostic comparisons revealed overlapping epicenters in schizophrenia and bipolar, but not major depressive disorder, suggestive of a pathophysiological continuity within the schizophrenia-bipolar-spectrum. In sum, cortical alterations over the course of schizophrenia robustly follow brain network architecture, emphasizing marked hub susceptibility and temporo-frontal epicenters at both the level of the group and the individual. Subtle variations of epicenters across disease stages suggest interacting pathological processes, while associations with patient-specific symptoms support additional inter-individual variability of hub vulnerability and epicenters in schizophrenia. Our work outlines potential pathways to better understand macroscale structural alterations, and inter- individual variability in schizophrenia.

Sex differences in the functional topography of association networks in youth.

Prior work has shown that there is substantial interindividual variation in the spatial distribution of functional networks across the cerebral cortex, or functional topography. However, it remains unknown whether there are sex differences in the topography of individualized networks in youth. Here, we leveraged an advanced machine learning method (sparsity-regularized non-negative matrix factorization) to define individualized functional networks in 693 youth (ages 8 to 23 y) who underwent functional MRI as part of the Philadelphia Neurodevelopmental Cohort. Multivariate pattern analysis using support vector machines classified participant sex based on functional topography with 82.9% accuracy (P < 0.0001). Brain regions most effective in classifying participant sex belonged to association networks, including the ventral attention, default mode, and frontoparietal networks. Mass univariate analyses using generalized additive models with penalized splines provided convergent results. Furthermore, transcriptomic data from the Allen Human Brain Atlas revealed that sex differences in multivariate patterns of functional topography were spatially correlated with the expression of genes on the X chromosome. These results highlight the role of sex as a biological variable in shaping functional topography.

Diffusion MRI head motion correction methods are highly accurate but impacted by denoising and sampling scheme.

Head motion correction is particularly challenging in diffusion-weighted MRI (dMRI) scans due to the dramatic changes in image contrast at different gradient strengths and directions. Head motion correction is typically performed using a Gaussian Process model implemented in FSL's Eddy. Recently, the 3dSHORE-based SHORELine method was introduced that does not require shell-based acquisitions, but it has not been previously benchmarked. Here we perform a comprehensive evaluation of both methods on realistic simulations of a software fiber phantom that provides known ground-truth head motion. We demonstrate that both methods perform remarkably well, but that performance can be impacted by sampling scheme and the extent of head motion and the denoising strategy applied before head motion correction. Furthermore, we find Eddy benefits from denoising the data first with MP-PCA. In sum, we provide the most extensive known benchmarking of dMRI head motion correction, together with extensive simulation data and a reproducible workflow. PRACTITIONER POINTS: Both Eddy and SHORELine head motion correction methods performed quite well on a large variety of simulated data. Denoising with MP-PCA can improve head motion correction performance when Eddy is used. SHORELine effectively corrects motion in non-shelled diffusion spectrum imaging data.

A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging.

Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.

Network enrichment significance testing in brain-phenotype association studies.

Functional networks often guide our interpretation of spatial maps of brain-phenotype associations. However, methods for assessing enrichment of associations within networks of interest have varied in terms of both scientific rigor and underlying assumptions. While some approaches have relied on subjective interpretations, others have made unrealistic assumptions about spatial properties of imaging data, leading to inflated false positive rates. We seek to address this gap in existing methodology by borrowing insight from a method widely used in genetics research for testing enrichment of associations between a set of genes and a phenotype of interest. We propose network enrichment significance testing (NEST), a flexible framework for testing the specificity of brain-phenotype associations to functional networks or other sub-regions of the brain. We apply NEST to study enrichment of associations with structural and functional brain imaging data from a large-scale neurodevelopmental cohort study.

Similarity-based multimodal regression.

To better understand complex human phenotypes, large-scale studies have increasingly collected multiple data modalities across domains such as imaging, mobile health, and physical activity. The properties of each data type often differ substantially and require either separate analyses or extensive processing to obtain comparable features for a combined analysis. Multimodal data fusion enables certain analyses on matrix-valued and vector-valued data, but it generally cannot integrate modalities of different dimensions and data structures. For a single data modality, multivariate distance matrix regression provides a distance-based framework for regression accommodating a wide range of data types. However, no distance-based method exists to handle multiple complementary types of data. We propose a novel distance-based regression model, which we refer to as Similarity-based Multimodal Regression (SiMMR), that enables simultaneous regression of multiple modalities through their distance profiles. We demonstrate through simulation, imaging studies, and longitudinal mobile health analyses that our proposed method can detect associations between clinical variables and multimodal data of differing properties and dimensionalities, even with modest sample sizes. We perform experiments to evaluate several different test statistics and provide recommendations for applying our method across a broad range of scenarios.

QSIPrep: an integrative platform for preprocessing and reconstructing diffusion MRI data.

Diffusion-weighted magnetic resonance imaging (dMRI) is the primary method for noninvasively studying the organization of white matter in the human brain. Here we introduce QSIPrep, an integrative software platform for the processing of diffusion images that is compatible with nearly all dMRI sampling schemes. Drawing on a diverse set of software suites to capitalize on their complementary strengths, QSIPrep facilitates the implementation of best practices for processing of diffusion images.

Traumatic stress load and stressor reactivity score associated with accelerated gray matter maturation in youths indexed by normative models.

Understanding how traumatic stress affects typical brain development during adolescence is critical to elucidate underlying mechanisms related to both maladaptive functioning and resilience after traumatic exposures. The current study aimed to map deviations from normative ranges of brain gray matter for youths with traumatic exposures. For each cortical and subcortical gray matter region, normative percentiles of variations were established using structural MRI from typically developing youths without any traumatic exposure (n = 245; age range = 8-23) from the Philadelphia Neurodevelopmental Cohort (PNC). The remaining PNC participants with neuroimaging data (n = 1129) were classified as either within the normative range (5-95%), delayed (>95%) or accelerated (<5%) maturational ranges for each region using the normative model. An averaged quantile regression index was calculated across all regions. Mediation models revealed that high traumatic stress load was positively associated with poorer cognitive functioning and greater psychopathology, and these associations were mediated by accelerated gray matter maturation. Furthermore, higher stressor reactivity scores, which represent a less resilient response under traumatic stress, were positively correlated with greater acceleration of gray matter maturation (r = 0.224, 95% CI = [0.17, 0.28], p < 0.001), suggesting that more accelerated maturation was linked to greater stressor response regardless of traumatic stress load. We conclude that traumatic stress is a source of deviation from normative brain development associated with poorer cognitive functioning and more psychopathology in the long run.

Alprazolam modulates persistence energy during emotion processing in first-degree relatives of individuals with schizophrenia: a network control study.

Schizophrenia is marked by deficits in facial affect processing associated with abnormalities in GABAergic circuitry, deficits also found in first-degree relatives. Facial affect processing involves a distributed network of brain regions including limbic regions like amygdala and visual processing areas like fusiform cortex. Pharmacological modulation of GABAergic circuitry using benzodiazepines like alprazolam can be useful for studying this facial affect processing network and associated GABAergic abnormalities in schizophrenia. Here, we use pharmacological modulation and computational modeling to study the contribution of GABAergic abnormalities toward emotion processing deficits in schizophrenia. Specifically, we apply principles from network control theory to model persistence energy - the control energy required to maintain brain activation states - during emotion identification and recall tasks, with and without administration of alprazolam, in a sample of first-degree relatives and healthy controls. Here, persistence energy quantifies the magnitude of theoretical external inputs during the task. We find that alprazolam increases persistence energy in relatives but not in controls during threatening face processing, suggesting a compensatory mechanism given the relative absence of behavioral abnormalities in this sample of unaffected relatives. Further, we demonstrate that regions in the fusiform and occipital cortices are important for facilitating state transitions during facial affect processing. Finally, we uncover spatial relationships (i) between regional variation in differential control energy (alprazolam versus placebo) and (ii) both serotonin and dopamine neurotransmitter systems, indicating that alprazolam may exert its effects by altering neuromodulatory systems. Together, these findings provide a new perspective on the distributed emotion processing network and the effect of GABAergic modulation on this network, in addition to identifying an association between schizophrenia risk and abnormal GABAergic effects on persistence energy during threat processing.

Psychosis brain subtypes validated in first-episode cohorts and related to illness remission: results from the PHENOM consortium.

Using machine learning, we recently decomposed the neuroanatomical heterogeneity of established schizophrenia to discover two volumetric subgroups-a 'lower brain volume' subgroup (SG1) and an 'higher striatal volume' subgroup (SG2) with otherwise normal brain structure. In this study, we investigated whether the MRI signatures of these subgroups were also already present at the time of the first-episode of psychosis (FEP) and whether they were related to clinical presentation and clinical remission over 1-, 3-, and 5-years. We included 572 FEP and 424 healthy controls (HC) from 4 sites (Sao Paulo, Santander, London, Melbourne) of the PHENOM consortium. Our prior MRI subgrouping models (671 participants; USA, Germany, and China) were applied to both FEP and HC. Participants were assigned into 1 of 4 categories: subgroup 1 (SG1), subgroup 2 (SG2), no subgroup membership ('None'), and mixed SG1 + SG2 subgroups ('Mixed'). Voxel-wise analyses characterized SG1 and SG2 subgroups. Supervised machine learning analyses characterized baseline and remission signatures related to SG1 and SG2 membership. The two dominant patterns of 'lower brain volume' in SG1 and 'higher striatal volume' (with otherwise normal neuromorphology) in SG2 were identified already at the first episode of psychosis. SG1 had a significantly higher proportion of FEP (32%) vs. HC (19%) than SG2 (FEP, 21%; HC, 23%). Clinical multivariate signatures separated the SG1 and SG2 subgroups (balanced accuracy = 64%; p < 0.0001), with SG2 showing higher education but also greater positive psychosis symptoms at first presentation, and an association with symptom remission at 1-year, 5-year, and when timepoints were combined. Neuromorphological subtypes of schizophrenia are already evident at illness onset, separated by distinct clinical presentations, and differentially associated with subsequent remission. These results suggest that the subgroups may be underlying risk phenotypes that could be targeted in future treatment trials and are critical to consider when interpreting neuroimaging literature.

Associations between neighborhood socioeconomic status, parental education, and executive system activation in youth.

Socioeconomic status (SES) can impact cognitive performance, including working memory (WM). As executive systems that support WM undergo functional neurodevelopment during adolescence, environmental stressors at both individual and community levels may influence cognitive outcomes. Here, we sought to examine how SES at the neighborhood and family level impacts task-related activation of the executive system during adolescence and determine whether this effect mediates the relationship between SES and WM performance. To address these questions, we studied 1,150 youths (age 8-23) that completed a fractal n-back WM task during functional magnetic resonance imaging at 3T as part of the Philadelphia Neurodevelopmental Cohort. We found that both higher neighborhood SES and parental education were associated with greater activation of the executive system to WM load, including the bilateral dorsolateral prefrontal cortex, posterior parietal cortex, and precuneus. The association of neighborhood SES remained significant when controlling for task performance, or related factors like exposure to traumatic events. Furthermore, high-dimensional multivariate mediation analysis identified distinct patterns of brain activity within the executive system that significantly mediated the relationship between measures of SES and task performance. These findings underscore the importance of multilevel environmental factors in shaping executive system function and WM in youth.

Morphological integration of the human brain across adolescence and adulthood.

Brain structural covariance norms capture the coordination of neurodevelopmental programs between different brain regions. We develop and apply anatomical imbalance mapping (AIM), a method to measure and model individual deviations from these norms, to provide a lifespan map of morphological integration in the human cortex. In cross-sectional and longitudinal data, analysis of whole-brain average anatomical imbalance reveals a reproducible tightening of structural covariance by age 25 y, which loosens after the seventh decade of life. Anatomical imbalance change in development and in aging is greatest in the association cortex and least in the sensorimotor cortex. Finally, we show that interindividual variation in whole-brain average anatomical imbalance is positively correlated with a marker of human prenatal stress (birthweight disparity between monozygotic twins) and negatively correlated with general cognitive ability. This work provides methods and empirical insights to advance our understanding of coordinated anatomical organization of the human brain and its interindividual variation.

ALE meta-analyses of voxel-based morphometry studies: Parameter validation via large-scale simulations.

Activation likelihood estimation (ALE) meta-analysis has been applied to structural neuroimaging data since long, but up to now, any systematic assessment of the algorithm's behavior, power and sensitivity has been based on simulations using functional neuroimaging databases as their foundation. Here, we aimed to determine whether the guidelines offered by previous evaluations can be generalized to ALE meta-analyses of voxel-based morphometry (VBM) studies. We ran 365000 distinct ALE analyses filled with simulated experiments, randomly sampling parameters from BrainMap's VBM experiment database. We then examined the algorithm's sensitivity, its susceptibility to spurious convergence, and its susceptibility to excessive contributions by individual experiments. In general, the performance of the ALE algorithm was highly comparable between imaging modalities, with the algorithm's sensitivity and specificity reaching similar levels with structural data as previously observed with functional data. Because of the lower number of foci reported and the higher number of participants usually included in structural experiments, individual studies had, on average, a higher impact towards significant clusters. To prevent significant clusters from being driven by single experiments, we recommend that researchers include at least 23 experiments in a VBM ALE dataset, instead of the previously recommended minimum of n = 17. While these recommendations do not constitute hard borders, running ALE analyses on smaller datasets would require special diligence in assessing and reporting the contributions of experiments to individual clusters.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization.

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.