Search

Ischemic stroke results in a loss of tissue homeostasis and integrity, the underlying pathobiology of which stems primarily from the depletion of cellular energy stores and perturbation of available metabolites1. Hibernation in thirteen-lined ground squirrels (TLGS),Ictidomys tridecemlineatus, provides a natural model of ischemic tolerance as these mammals undergo prolonged periods of critically low cerebral blood flow without evidence of central nervous system (CNS) damage2. Studying the complex interplay of genes and metabolites that unfolds during hibernation may provide novel insights into key regulators of cellular homeostasis during brain ischemia. Herein, we interrogated the molecular profiles of TLGS brains at different time points within the hibernation cycle via RNA sequencing coupled with untargeted metabolomics. We demonstrate that hibernation in TLGS leads to major changes in the expression of genes involved in oxidative phosphorylation and this is correlated with an accumulation of the tricarboxylic acid (TCA) cycle intermediates citrate, cis-aconitate, and α-ketoglutarate-αKG. Integration of the gene expression and metabolomics datasets led to the identification of succinate dehydrogenase (SDH) as the critical enzyme during hibernation, uncovering a break in the TCA cycle at that level. Accordingly, the SDH inhibitor dimethyl malonate (DMM) was able to rescue the effects of hypoxia on human neuronal cellsin vitroand in mice subjected to permanent ischemic strokein vivo. Our findings indicate that studying the regulation of the controlled metabolic depression that occurs in hibernating mammals may lead to novel therapeutic approaches capable of increasing ischemic tolerance in the CNS.

BackgroundBrain MRIs acquired in clinical settings represent a valuable and underutilized scientific resource for investigating neurodevelopment. Utilization of these clinical scans has been limited because of their clinical acquisition and technical heterogeneity. These barriers have curtailed the interpretability and scientific value of retrospective studies of clinically acquired brain MRIs, compared to studies of prospectively acquired research quality brain MRIs.PurposeTo develop a scalable and rigorous approach to generate clinical brain growth chart models, to benchmark neuroanatomical differences in clinical MRIs, and to validate clinically-derived brain growth charts against those derived from large-scale research studies.Materials and MethodsWe curated a set of clinical MRIScans withLimitedImagingPathology (SLIP) – so-called “clinical controls” – from an urban pediatric healthcare system acquired between 2005 and 2020. The curation process included manual review of signed radiology reports, as well as automated and manual quality review of images without gross pathology. We measured global and regional volumetric imaging phenotypes in the SLIP sample using two alternative, advanced image processing pipelines, and quantitatively compared clinical brain growth charts to research brain growth charts derived from >123,000 MRIs.ResultsThe curated SLIP dataset included 372 patients scanned between the ages of 28 days post-birth and 22.2 years across nine 3T MRI scanners. Clinical brain growth charts were highly similar to growth charts derived from large-scale research datasets, in terms of the normative developmental trajectories predicted by the models. The clinical indication of the scans did not significantly bias the output of clinical brain charts. Tens of thousands of additional healthcare system scans meet inclusion criteria to be included in future brain growth charts.ConclusionBrain charts derived from clinical-controls are highly similar to brain charts from research-controls, suggesting that curated clinical scans could be used to supplement research datasets.Summary StatementBrain growth charts of pediatric clinical MRIs with limited imaging pathology (N=372) are highly correlated with charts from a large aggregated set of research controls (N>120,000).Key ResultsA cohort of brain MRI scans with limited reported imaging pathology (N=372, 186 female; ages 0.07 - 22.2 years, median = 10.2) were identified using signed radiology reports and processed using two segmentation pipelines. Growth charts generated from these scans are highly correlated with growth charts from a large aggregated set of research controls (r range 0.990 - 0.999). There was no evidence of bias due to the reason for each scan.

The relationship between obesity and brain structure is incompletely understood. Using diffusion-weighted MRI from ∼30,000 UK Biobank participants we test the hypothesis that obesity (waist-to-hip ratio, WHR) is associated with regional differences in two micro-structural MRI metrics: isotropic volume fraction (ISOVF), an index of free water, and intra-cellular volume fraction (ICVF), an index of neurite density. We observed significant associations with obesity in two coupled but distinct brain systems: a prefrontal-temporal-striatal system associated with ISOVF and a medial temporal-occipital-striatal system associated with ICVF. The ISOVF∼WHR system colocated with expression of genes enriched for innate immune functions, decreased glial density, and high mu opioid (MOR) and other neurotransmitter receptor density. Conversely, the ICVF∼WHR system co-located with expression of genes enriched for G-protein coupled receptors and decreased density of MOR and other receptors. To test whether these distinct brain phenotypes might differ in terms of their underlying shared genetics or relationship to maps of the inflammatory marker C-reactive Protein (CRP), we estimated the genetic correlations between WHR and ISOVF (rg= 0.026,P= 0.36) and ICVF (rg= 0.112,P <9 × 10−4) as well as comparing correlations between WHR maps and equivalent CRP maps for ISOVF and ICVF (p

Structural similarity networks are a central focus of magnetic resonance imaging (MRI) research into human brain connectomes in health and disease. We present Morphometric INverse Divergence (MIND), a robust method to estimate within-subject structural similarity between cortical areas based on the Kullback-Leibler divergence between the multivariate distributions of their structural features. Compared to the prior approach of morphometric similarity networks (MSNs) on N>10,000 data from the ABCD cohort, MIND networks were more consistent with known cortical symmetry, cytoarchitecture, and (in N=19 macaques) gold-standard tract-tracing connectivity, and were more invariant to cortical parcellation. Importantly, MIND networks were remarkably coupled with cortical gene co-expression, providing fresh evidence for the unified architecture of brain structure and transcription. Using kinship (N=1282) and genetic data (N=4085), we characterized the heritability of MIND phenotypes, identifying stronger genetic influence on the relationship between structurally divergent regions compared to structurally similar regions. Overall, MIND presents a biologically-validated lens for analyzing the structural organization of the cortex using readily-available MRI measurements.

The Allen Human Brain Atlas (AHBA) has catalysed many discoveries linking macroscale brain organisation to microscale neurobiology, including the existence of a transcriptional axis stretching from sensorimotor to association cortices. However, the complex spatial structure of the human cortex likely includes multiple gradients of gene expression that have yet to be characterised. Here, we reveal two additional transcriptional axes that are anatomically differentiated and have distinct profiles of metabolic and immune function, cytoarchitecture, and developmental trajectory. Moreover, we validate the three transcriptional axes in external data, show they explain more than half of spatial variation across key anatomical and functional neuroimaging maps, and demonstrate they are dissociably aligned with maps of autism, depression, and schizophrenia from one of the largest case-control imaging datasets. Together, these results reveal how human brain anatomy, development, and disorders are patterned along multiple cortical axes underpinned by the intrinsic structure of gene transcription.

Eating disorders (EDs) are complex psychiatric conditions that often emerge during adolescence, and affected individuals frequently demonstrate high rates of psychiatric comorbidity, particularly with depressive and anxiety disorders. Although risk for EDs reflects both genetic and neurobiological factors, knowledge of how genetic risk for EDs relates to neurobiology and psychiatric symptoms during critical developmental periods remains limited. We therefore implemented a novel multivariate framework, which sought to advance knowledge of the etiology of EDs by simultaneously estimating associations between genetic risk, brain structure and ED-related psychopathology symptoms in over 4,500 adolescents of European ancestry from the Adolescent Brain and Cognitive Development study (M(SD)age=119.29(7.49) months). Polygenic scores for anorexia nervosa (AN PGS) and body mass index (BMI PGS) were generated and related to three morphometric brain features— cortical thickness, surface area and subcortical grey matter volume—and to latent psychopathology factors using structural equation modeling. We identified a three-factor structure of ED-related psychopathology symptoms: eating, distress and fear factors. Increased BMI PGS were uniquely associated with greater eating factor scores, whereas AN PGS were unrelated to psychopathology factors. Moreover, genetic risk for high BMI and for AN had distinct neural correlates, where greater BMI PGS predicted widespread increases in cortical thickness and reductions in surface area while AN PGS were nominally related to reduced caudate volume. Altered default mode and visual network thickness was associated with greater eating factor scores, whereas distress and fear factor scores reflected a shared reduction in somatomotor network thickness. Our novel findings indicate that greater genetic risk for high BMI and altered cortical thickness of canonical brain networks underpin ED symptomatology in early adolescence. As neurobiological factors appear to shape disordered eating earlier in the life course than previously thought, these results underscore the need for early detection and intervention efforts for EDs.

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, how these gradients relate to adult brain function, and whether they 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 predicts the location of unseen brain tissue samples, ii) delineates known functional territories, and iii) explains the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical functional hierarchies differentiating primary sensory cortex from 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 non-human primate datasets, and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well known morphogens (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.

Animal studies of neurodevelopmental plasticity have shown that intrinsic brain activity evolves from high amplitude and globally synchronized to suppressed and sparse as plasticity declines and the cortex matures. Leveraging resting-state functional MRI data from 1033 individuals (8-23 years), we reveal that this stereotyped refinement of intrinsic activity occurs during human development and provides evidence for a cortical gradient of neurodevelopmental plasticity during childhood and adolescence. Specifically, we demonstrate that declines in the amplitude of intrinsic activity are initiated heterochronously across regions, coupled to the maturation of a plasticity-restricting structural feature, and temporally staggered along a hierarchical sensorimotor-association axis from ages 8 to 18. Youth from disadvantaged environments exhibit reduced intrinsic activity in regions further up the sensorimotor-association axis, suggestive of a reduced level of plasticity in late-maturing cortices. Our results uncover a hierarchical axis of neurodevelopment and offer insight into the temporal sequence of protracted neurodevelopmental plasticity in humans.

SummaryThere has been an immense research focus on the topic of cortical reorganization in human evolution, but much less is known regarding the reorganization of subcortical circuits which are intimate working partners of the cortex. Here, by combining advanced image analysis techniques with comparative neuroimaging data, we systematically map organizational differences in striatal, pallidal and thalamic anatomy between humans and chimpanzees. We relate interspecies differences, a proxy for evolutionary changes, to genetics and behavioral correlates. We show that highly heritable morphological measures are significantly expanded across species, in contrast to previous findings in the cortex. The identified morphological-cognitive latent variables were associated with striatal expansion, and affective latent variables were associated with more evolutionarily-conserved areas in the thalamus and globus pallidus. These findings provide new insight into the architecture of these subcortical hubs and can provide greater information on the role of these structures in health and illness.

Recent work in imaging genetics suggests high levels of genetic overlap within cortical regions for cortical thickness (CT) and surface area (SA). We model this multivariate system of genetic relationships by applying Genomic Structural Equation Modeling (Genomic SEM) and parsimoniously define five genomic brain factors underlying both CT and SA along with a general factor capturing genetic overlap across all brain regions. We validate these factors by demonstrating the generalizability of the model to a semi-independent sample and show that the factors align with biologically and functionally relevant parcellations of the cortex. We apply Stratified Genomic SEM to identify specific categories of genes (e.g., neuronal cell types) that are disproportionately associated with pleiotropy across specific subclusters of brain regions, as indexed by the genomic factors. Finally, we examine genetic associations with psychiatric and cognitive correlates, finding that broad aspects of cognitive function are associated with a general factor for SA and that psychiatric associations are null. These analyses provide key insights into the multivariate genomic architecture of two critical features of the cerebral cortex.

Imaging technologies are increasingly used to generate high-resolution reference maps of brain structure and function. Modern scientific discovery relies on making comparisons between new maps (e.g. task activations, group structural differences) and these reference maps. Although recent data sharing initiatives have increased the accessibility of such brain maps, data are often shared in disparate coordinate systems (or "spaces"), precluding systematic and accurate comparisons among them. Here we introduce the neuromaps toolbox, an open-access software package for accessing, transforming, and analyzing structural and functional brain annotations. We implement two registration frameworks to generate high-quality transformations between four standard coordinate systems commonly used in neuroimaging research. The initial release of the toolbox features >40 curated reference maps and biological ontologies of the human brain, including maps of gene expression, neurotransmitter receptors, metabolism, neurophysiological oscillations, developmental and evolutionary expansion, functional hierarchy, individual functional variability, and cognitive specialization. Robust quantitative assessment of map-to-map similarity is enabled via a suite of spatial autocorrelation-preserving null models. Finally, we demonstrate two examples of how neuromaps can be used to contextualize brain maps with respect to canonical annotations. By discovering novel associations with previously-established features of brain structure and function, neuromaps generates biological insight about new brain maps. Altogether, 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.

Diffuse gliomas have been hypothesized to originate from neural stem cells in the subventricular zone. Here, we evaluated this hypothesis by mapping independent sources of glioma localization and determining their relationships with neurogenic niches, genetic markers, and large-scale connectivity networks. Using lesion data from a total of 410 patients with glioma, we identified -- and replicated in an independent sample -- three lesion covariance networks (LCNs), which reflect clusters of frequent glioma co-localization. Each LCN overlapped with a distinct horn of the lateral ventricles. The first LCN, which overlapped with the anterior horn, was associated with low-grade, IDH-mutated/1p19q-codeleted tumors, as well as a neural transcriptomic signature and improved overall survival. Each LCN significantly corresponded with multiple brain networks, with LCN1 bearing an especially strong relationship with structural and functional connectivity, consistent with its neural transcriptomic profile. Finally, we identified subcortical, periventricular structures with functional connectivity patterns to the cortex that significantly matched each LCN. Cumulatively, our findings support a model wherein periventricular brain connectivity guides tumor development.

Current neuroimaging acquisition and processing approaches tend to be optimised for quality rather than speed. However, rapid acquisition and processing of neuroimaging data can lead to novel neuroimaging paradigms, such as adaptive acquisition, where rapidly processed data is used to inform subsequent image acquisition steps. Here we first evaluate the impact of several processing steps on the processing time and quality of registration of manually labelled T1-weighted MRI scans. Subsequently, we apply the selected rapid processing pipeline both to rapidly acquired multicontrast EPImix scans of 95 participants (which include T1-FLAIR, T2, T2*, T2-FLAIR, DWI & ADC contrasts, acquired in ∼1 minute), as well as to slower, more standard single-contrast T1-weighted scans of a subset of 66 participants. We quantify the correspondence between EPImix and single-contrast T1-weighted scans, using correlations between voxels and regions of interest across participants, measures of within- and between-participant identifiability as well as regional structural covariance networks. Furthermore, we explore the use of EPImix for the rapid construction of morphometric similarity networks. Finally, we quantify the reliability of EPImix-derived data using test-retest scans of 10 participants. Our results demonstrate that quantitative information can be derived from a neuroimaging scan acquired and processed within minutes, which could further be used to implement adaptive multimodal imaging and tailor neuroimaging examinations to individual patients.Abstract FigureGraphical abstract.

BackgroundRecent discovery of approximately 270 common genetic variants associated with schizophrenia has enabled polygenic risk scores (PRS) to be measured in the population. We hypothesized that normal variation in PRS would be associated with magnetic resonance imaging (MRI) phenotypes of brain morphometry and tissue composition.MethodsWe used the largest extant genome-wide association dataset (N = 69,369 cases and N = 236,642 healthy controls) to measure PRS for schizophrenia in a large sample of adults from the UK Biobank (Nmax = 29,878) who had multiple micro- and macro-structural MRI metrics measured at each of 180 cortical areas, seven subcortical structures, and 15 major white matter tracts. Linear mixed effect models were used to investigate associations between PRS and brain structure at global and regional scales, controlled for multiple comparisons.ResultsPolygenic risk was significantly associated with reduced neurite density index (NDI) at global brain scale, at 149 cortical regions, five subcortical structures and 14 white matter tracts. Other micro-structural parameters, e.g., fractional anisotropy, that were correlated with NDI were also significantly associated with PRS. Genetic effects on multiple MRI phenotypes were co-located in temporal, cingulate and prefrontal cortical areas, insula, and hippocampus. Post-hoc bidirectional Mendelian randomization analyses provided preliminary evidence in support of a causal relationship between (reduced) thalamic NDI and (increased) risk of schizophrenia.ConclusionsRisk-related reduction in NDI is plausibly indicative of reduced density of myelinated axons and dendritic arborization in large-scale cortico-subcortical networks. Cortical, subcortical and white matter micro-structure may be linked to the genetic mechanisms of schizophrenia.

Drug-resistant focal epilepsy is often caused by focal cortical dysplasias (FCDs). The distribution of these lesions across the cerebral cortex and the impact of lesion location on clinical presentation and surgical outcome are largely unknown. We created a neuroimaging cohort of patients with individually mapped FCDs to determine factors associated with lesion location and predictors of postsurgical outcome.The Multi-centre Epilepsy Lesion Detection (MELD) project collated a retrospective cohort of 580 patients with epilepsy attributed to FCD from 20 epilepsy centres worldwide. MRI-based maps of individual FCDs with accompanying demographic, clinical and surgical information were collected. We mapped the distribution of FCDs, examined for associations between clinical factors and lesion location, and developed a predictive model of postsurgical seizure freedom.FCDs were non-uniformly distributed, concentrating in the superior frontal sulcus, frontal pole and temporal pole. Epilepsy onset was typically before age 10. Earlier epilepsy onset was associated with lesions in primary sensory areas while later epilepsy onset was associated with lesions in association cortices. Lesions in temporal and occipital lobes tended to be larger than frontal lobe lesions. Seizure freedom rates varied with FCD location, from around 30% in visual, motor and premotor areas to 75% in superior temporal and frontal gyri. The predictive model of postsurgical seizure freedom had a positive predictive value of 70% and negative predictive value of 61%.FCD location is an important determinant of its size, the age of epilepsy onset and the likelihood of seizure freedom post-surgery. Our atlas of lesion locations can be used to guide the radiological search for subtle lesions in individual patients. Our atlas of regional seizure freedom rates and associated predictive model can be used to estimate individual likelihoods of postsurgical seizure freedom. Data-driven atlases and predictive models are essential for evidence-based, precision medicine and risk counselling in neurology.

Cortical regional identities develop through anterior-posterior (A-P) and dorsal-ventral (D-V) prenatal genomic patterning gradients. Here we find that A-P and D-V genomic patterning of cortical surface area (SA) and thickness (CT) is intact in typically developing and autistic toddlers with good language outcome, but is absent in autistic toddlers with poor early language outcome. Genes driving this effect are prominent in midgestational A-P and D-V gene expression gradients and prenatal cell types driving SA and CT variation (e.g., progenitor cells versus excitatory neurons). These genes are also important for vocal learning, human-specific evolution, and prenatal co-expression networks enriched for high-penetrance autism risk genes. Autism with poor early language outcome may be linked to atypical genomic cortical patterning starting in prenatal periods and which impacts later development of regional functional specialization and circuit formation.One Sentence SummaryGenomic patterning of the cortex is atypical in autistic toddlers with poor early language outcome.

Functional neuroimaging research in the non-human primate (NHP) has been advancing at a remarkable rate. The increase in available data establishes a need for robust analysis pipelines designed for NHP neuroimaging and accompanying template spaces to standardize the localization of neuroimaging results. Our group recently developed the NIMH Macaque Template (NMT), a high-resolution population average anatomical template and associated neuroimaging resources, providing researchers with a standard space for macaque neuroimaging (Seidlitz, Sponheim et al., 2018). Here, we release NMT v2, which includes both symmetric and asymmetric templates in stereotaxic orientation, with improvements in spatial contrast, processing efficiency, and segmentation. We also introduce the Cortical Hierarchy Atlas of the Rhesus Macaque (CHARM), a hierarchical parcellation of the macaque cerebral cortex with varying degrees of detail. These tools have been integrated into the neuroimaging analysis software AFNI (Cox, 1996) to provide a comprehensive and robust pipeline for fMRI processing, visualization and analysis of NHP data. AFNI’s new @animal_warper program can be used to efficiently align anatomical scans to the NMT v2 space, and afni_proc.py integrates these results with full fMRI processing using macaque-specific parameters: from motion correction through regression modeling. Taken together, the NMT v2 and AFNI represent an all-in-one package for macaque functional neuroimaging analysis, as demonstrated with available demos for both task and resting state fMRI.HighlightsThe NMT v2, a stereotaxically aligned symmetric macaque template, is introduced.A new atlas (CHARM), defined on NMT v2, parcellates the cortex at six spatial scales.AFNI’s @animal_warper aligns and maps data between monkey anatomicals and templates.AFNI’s afni_proc.py facilitates monkey fMRI analysis with automated scripting and QC.Demos of macaque task and resting state fMRI analysis with these tools are provided.

The MNI CIVET pipeline for automated extraction of cortical surfaces and evaluation of cortical thickness fromin-vivohuman MRI has been extended for processing macaque brains. Processing is performed based on the NIMH Macaque Template (NMT), as the reference template, with the anatomical parcellation of the surface following the D99 and CHARM atlases. The modifications needed to adapt CIVET to the macaque brain are detailed. Results have been obtained using CIVET-macaque to process the anatomical scans of the 31 macaques used to generate the NMT and another 95 macaques from the PRIME-DE initiative. It is anticipated that the open usage of CIVET-macaque will promote collaborative efforts in data collection and processing, sharing, and automated analyses from which the non-human primate brain imaging field will advance.

Neuroimaging non-human primates (NHPs) is a growing, yet highly specialized field of neuroscience. Resources that were primarily developed for human neuroimaging often need to be significantly adapted for use with NHPs or other animals, which has led to an abundance of custom, in-house solutions. In recent years, the global NHP neuroimaging community has made significant efforts to transform the field towards more open and collaborative practices. Here we present the PRIMatE Resource Exchange (PRIME-RE), a new collaborative online platform for NHP neuroimaging. PRIME-RE is a dynamic community-driven hub for the exchange of practical knowledge, specialized analytical tools, and open data repositories, specifically related to NHP neuroimaging. PRIME-RE caters to both researchers and developers who are either new to the field, looking to stay abreast of the latest developments, or seeking to collaboratively advance the field. HighlightsO_LIPRIME-RE is a bottom-up open science platform for non-human primate neuroimaging. C_LIO_LIThe website hosts a wiki and info about analytical tools and open data repositories. C_LIO_LISpecialized tools address structural, functional & diffusion-weighted analysis. C_LIO_LIMultiple templates & atlases for several non-human primate species are highlighted. C_LIO_LICommunication channels facilitate international conversation and collaboration C_LI

Regulation of gene expression drives protein interactions that govern synaptic wiring and neuronal activity. The resulting coordinated activity among neuronal populations supports complex psychological processes, yet how gene expression shapes cognition and emotion remains unknown. Here we directly bridge the microscale and macroscale by mapping gene expression patterns to functional activation patterns across the cortical sheet. Applying unsupervised learning to the Allen Human Brain Atlas and Neurosynth databases, we identify a ventromedial-dorsolateral gradient of gene assemblies that separate affective and cognitive domains. This topographic molecular-psychological signature reflects the hierarchical organization of the neocortex, including systematic variations in cell type, myeloarchitecture, laminar differentiation, and intrinsic network affiliation. In addition, this molecular-psychological signature is related to individual differences in cognitive performance, strengthens over neurodevelopment, and can be replicated in two independent repositories. Collectively, our results reveal spatially covarying transcriptomic and cognitive architectures, highlighting the influence that molecular mechanisms exert on psychological processes.

We hypothesized that there are sexual differences in human brain network development underlying the female > male divergence in adolescent depression. We tested for sex differences in parameters of brain network development (accelerated longitudinal fMRI, N=298 healthy adolescents, each scanned 1 to 3 times). Sexually divergent development of functional connectivity was located in default mode network (DMN), limbic cortex, and subcortical nuclei. Females had a more “disruptive” pattern of development, where weak functional connectivity at age 14 became stronger during adolescence. This fMRI-derived map of divergent adolescent development was co-located with (i) a map of functional dysconnectivity associated with adult major depressive disorder (MDD); and (ii) an adult brain gene expression pattern enriched for genes on the X chromosome, neurodevelopmental genes, and risk genes for MDD. Sexual divergence in disruptive development of DMN, limbic and subcortical functional networks is potentially relevant to the increased risk of depression in adolescent females.

Social-communication (SC) and restricted repetitive behaviors (RRB) are autism diagnostic symptom domains. SC and RRB severity can markedly differ within and between individuals and may be underpinned by different neural circuitry and genetic mechanisms. Modeling SC-RRB balance could help identify how neural circuitry and genetic mechanisms map onto such phenotypic heterogeneity. Here we developed a phenotypic stratification model that makes highly accurate (97-99%) out-of-sample SC=RRB, SC>RRB, and RRB>SC subtype predictions. Applying this model to resting state fMRI data from the EU-AIMS LEAP dataset (n=509), we find that while the phenotypic subtypes share many commonalities in terms of intrinsic functional connectivity, they also show replicable differences within some networks compared to a typically-developing group (TD). Specifically, the somatomotor network is hypoconnected with perisylvian circuitry in SC>RRB and visual association circuitry in SC=RRB. The SC=RRB subtype show hyperconnectivity between medial motor and anterior salience circuitry. Genes that are highly expressed within these networks show a differential enrichment pattern with known autism-associated genes, indicating that such circuits are affected by differing autism-associated genomic mechanisms. These results suggest that SC-RRB imbalance subtypes share many commonalities, but also express subtle differences in functional neural circuitry and the genomic underpinnings behind such circuitry.

Ageing is commonly associated with changes to segregation and integration of functional brain networks, but, in isolation, current network-based approaches struggle to elucidate changes across the many axes of functional organisation. However, the advent of gradient mapping techniques to neuroimaging provides a new means of studying functional organisation in a multi-dimensional connectivity space. Here, we studied ageing and behaviourally-relevant differences in a three-dimensional connectivity space using the Cambridge Centre for Ageing Neuroscience cohort (n=643). Building on gradient mapping techniques, we developed a set of measures to quantify the dispersion within and between functional communities. We detected a strong shift of the visual network across the lifespan from an extreme to a more central position in the 3D gradient space. In contrast, the dispersion distance of transmodal communities (dorsal attention, ventral attention, frontoparietal and default mode) did not change. However, these communities were increasingly dispersed with increasing age, reflecting more dissimilar functional connectivity profiles within each community. Increasing dispersion of frontoparietal and ventral attention networks, in particular, was associated negatively with cognition, measured by fluid intelligence. By using a technique that explicitly captures the ordering of functional systems in a multi-dimensional hierarchical framework, we identified behaviorally-relevant age-related differences of within and between network organisation. We propose that the study of functional gradients across the lifespan could provide insights that may facilitate the development of new strategies to maintain cognitive ability across the lifespan in health and disease.

Funder: Royal Children's Hospital Foundation; funder-id: http://dx.doi.org/10.13039/100014607, Funder: FP7 Ideas: European Research Council (); Grant(s): 319456, Interruption to gestation through preterm birth can significantly impact cortical development and have long-lasting adverse effects on neurodevelopmental outcome. We compared cortical morphology captured by high-resolution, multimodal magnetic resonance imaging (MRI) in n = 292 healthy newborn infants (mean age at birth = 39.9 weeks) with regional patterns of gene expression in the fetal cortex across gestation (n = 156 samples from 16 brains, aged 12 to 37 postconceptional weeks [pcw]). We tested the hypothesis that noninvasive measures of cortical structure at birth mirror areal differences in cortical gene expression across gestation, and in a cohort of n = 64 preterm infants (mean age at birth = 32.0 weeks), we tested whether cortical alterations observed after preterm birth were associated with altered gene expression in specific developmental cell populations. Neonatal cortical structure was aligned to differential patterns of cell-specific gene expression in the fetal cortex. Principal component analysis (PCA) of 6 measures of cortical morphology and microstructure showed that cortical regions were ordered along a principal axis, with primary cortex clearly separated from heteromodal cortex. This axis was correlated with estimated tissue maturity, indexed by differential expression of genes expressed by progenitor cells and neurons, and engaged in stem cell differentiation, neuron migration, and forebrain development. Preterm birth was associated with altered regional MRI metrics and patterns of differential gene expression in glial cell populations. The spatial patterning of gene expression in the developing cortex was thus mirrored by regional variation in cortical morphology and microstructure at term, and this was disrupted by preterm birth. This work provides a framework to link molecular mechanisms to noninvasive measures of cortical development in early life and highlights novel pathways to injury in neonatal populations at increased risk of neurodevelopmental disorder.

In the perinatal brain, regional cortical architecture and connectivity lay the foundations for functional circuits and emerging behaviour. Interruptions or atypical development during or before this period may therefore have long-lasting consequences. However, to be able to investigate these deviations, we need a measure of how this architecture evolves in the typically developing brain. To this end, in a large cohort of 241 term-born infants we used Magnetic Resonance Imaging to estimate cortical profiles based on morphometry and microstructure over the perinatal period (37-44 weeks post-menstrual age, PMA). Using the covariance of these profiles as a measure of inter-areal network similarity (Morphometric Similarity Networks; MSN), we clustered these networks into distinct modules. The resulting modules were consistent and symmetric, and corresponded to known functional distinctions, including sensory-motor, limbic and association regions and were spatially mapped onto known cytoarchitectonic tissue classes. Posterior (parietal, occipital) regions became more morphometrically similar with increasing PMA, while peri-cingulate and medial temporal regions became more dissimilar. Network strength was associated with PMA: Within-network similarity increased over PMA suggesting emerging network distinction. These changes in cortical network architecture over an eight-week period are consistent with, and likely underpin, the highly dynamic behavioural and cognitive development occurring during this critical period. The resulting cortical profiles might provide normative reference to investigate atypical early brain development.

The vast net of fibres within and underneath the cortex is optimised to support the convergence of different levels of brain organisation. Here we propose a novel coordinate system of the human cortex based on an advanced model of its connectivity. Our approach is inspired by seminal, but so far largely neglected models of cortico-cortical wiring established by post mortem anatomical studies and capitalizes on cutting-edge neuroimaging and machine learning. The new model expands the currently prevailing diffusion MRI tractography approach by incorporation of additional features of cortical microstructure and cortico-cortical proximity. Studying several datasets, we could show that our coordinate system robustly recapitulates established sensory-limbic and anterior-posterior dimensions of brain organisation. A series of validation experiments showed that the new wiring space reflects cortical microcircuit features (including pyramidal neuron depth and glial expression) and allowed for competitive simulations of functional connectivity and dynamics across a broad range contexts (based on resting-state fMRI, task-based fMRI, and human intracranial EEG coherence). Our results advance our understanding of how cell-specific neurobiological gradients produce a hierarchical cortical wiring scheme that is concordant with increasing functional sophistication of human brain organisation. Our evaluations demonstrate the cortical wiring space bridges across scales of neural organisation and can be easily translated to single individuals.

The biological processes underpinning adolescent brain maturation remain elusive. Expanding on previous work showing age-related changes in cortical morphology, we studied an accelerated longitudinal cohort of adolescents and young adults (n=223, two time points) to investigate dynamic reconfigurations in myeloarchitecture. Intracortical profiles were generated using magnetization transfer (MT) data, a myelin-sensitive magnetic resonance imaging contrast. Mixed-effect models of depth specific intracortical profiles demonstrated two separate processes i) related to overall increases in MT, and ii) showing a flattening of the MT profile related to enhanced signal in mid-to-deeper layers, especially in heteromodal and unimodal association cortices. This development was independent of morphological changes, and enhanced MT in mid-to-deeper layers was found to spatially co-localise specifically with gene expression markers of oligodendrocytes. Covariance analysis between all pairs of intracortical profiles revealed that these intracortical changes contributed to a gradual and dynamic differentiation from higher-order to lower-order systems. Depth-dependent trajectories of intracortical myeloarchitectural development contribute to the maturation of structural hierarchies in the human neocortex, providing a model for adolescent development that bridges microstructural and macroscopic scales of brain organization.eLife digestIntracortical myelin imposes a spatial structure on cortico-cortical connections, yet little is known about how myeloarchitecture develops throughout youth. We formulated a novel approach to study cortical myeloarchitecture in individual humans and leveraged an accelerated longitudinal design to track age-related changes from 14-27 years. We discovered two unique processes: one involving increasing mean myelin and another characterised by the preferential accumulation of myelin in mid-to-deeper cortical layers. Both processes contributed to an increasing segregation of lower-order from higher-order systems along the macroscale cortical hierarchy. These findings illustrate how layer specific microstructural changes contribute to the maturation of cortical organization and suggest adolescent fine tuning of hierarchical gradients of cortical networks.

Cortical development during childhood and adolescence has been characterised in recent years using metrics derived from Magnetic Resonance Imaging (MRI). Changes in cortical thickness are greatest in the first two decades of life and recapitulate the genetic organisation of the cortex, highlighting the potential early impact of gene expression on differences in cortical architecture over the lifespan. It is important to further our understanding of the possible neurobiological mechanisms that underlie these changes as differences in cortical thickness may act as a potential phenotypic marker of several common neurodevelopmental and psychiatric disorders.In this study, we combine MRI acquired from a large typically-developing childhood population (n=768) with comprehensive human gene expression databases to test the hypothesis that disrupted mechanisms common to neurodevelopmental disorders are encoded by genes expressed early in development and nested within those associated with typical cortical remodelling in childhood.We find that differential rates of thinning across the developing cortex are associated with spatially-varying gradients of gene expression. Genes that are expressed highly in regions of accelerated thinning are expressed predominantly in cortical neurons, involved in synaptic remodeling, and associated with common cognitive and neurodevelopmental disorders. Further, we identify subsets of genes that are highly expressed in the prenatal period and jointly associated with both developmental cortical morphology and neurodevelopmental disorders.

Adolescent changes in human brain function are not entirely understood. Here we used multi-echo functional magnetic resonance imaging (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in N=298 healthy adolescents. Participants were aged 14-26 years and were scanned on two or more occasions at least 6 months apart. We found two distinct modes of age-related change in FC: “conservative” and “disruptive”. Conservative development was characteristic of primary cortex, which was strongly connected at 14 years and became even more connected in the period 14-26 years. Disruptive development was characteristic of association cortex, hippocampus and amygdala, which were not strongly connected at 14 years but became more strongly connected during adolescence. We defined the maturational index (MI) as the signed coefficient of the linear relationship between baseline FC (at 14 years,FC14) and adolescent change in FC (∆FC14−26). Disruptive systems (with negative MI) were functionally specialised for social cognition and autobiographical memory and were significantly co-located with prior maps of aerobic glycolysis (AG), AG-related gene expression, post-natal expansion of cortical surface area, and adolescent shrinkage of cortical depth. We conclude that human brain organization is disrupted during adolescence by the emergence of strong functional connectivity of subcortical nuclei and association cortical areas, representing metabolically expensive re-modelling of synaptic connectivity between brain regions that were not strongly connected in childhood. We suggest that this re-modelling process may support emergence of social skills and self-awareness during healthy human adolescence.

Neurodevelopmental disorders are highly heritable and associated with spatially-selective disruptions of brain anatomy. The logic that translates genetic risks into spatially patterned brain vulnerabilities remains unclear but is a fundamental question in disease pathogenesis. Here, we approach this question by integrating (i)in vivoneuroimaging data from patient subgroups with known causal genomic copy number variations (CNVs), and (ii) bulk and single-cell gene expression data from healthy cortex. First, for each of six different CNV disorders, we show that spatial patterns of cortical anatomy change in youth are correlated with spatial patterns of expression for CNV region genes in bulk cortical tissue from typically-developing adults. Next, by transforming normative bulk-tissue cortical expression data into cell-type expression maps, we further link each disorder’s anatomical change map to specific cell classes and specific CNV-region genes that these cells express. Finally, we establish convergent validity of this “transcriptional vulnerability model” by inter-relating patient neuroimaging data with measures of altered gene expression in both brain and blood-derived patient tissue. Our work clarifies general biological principles that govern the mapping of genetic risks onto regional brain disruption in neurodevelopmental disorders. We present new methods that can harness these principles to screen for potential cellular and molecular determinants of disease from readily available patient neuroimaging data.

Schizophrenia has been conceived as a disorder of brain connectivity but it is unclear how this network phenotype is related to the emerging genetics. We used morphometric similarity analysis of magnetic resonance imaging (MRI) data as a marker of inter-areal cortical connectivity in three prior case-control studies of psychosis: in total, N=185 cases and N=227 controls. Psychosis was associated with globally reduced morphometric similarity (MS) in all 3 studies. There was also a replicable pattern of case-control differences in regional MS which was significantly reduced in patients in frontal and temporal cortical areas, but increased in parietal cortex. Using prior brain-wide gene expression data, we found that the cortical map of case-control differences in MS was spatially correlated with cortical expression of a weighted combination of genes enriched for neu-robiologically relevant ontology terms and pathways. In addition, genes that were normally over-expressed in cortical areas with reduced MS were significantly up-regulated in a prior post mortem study of schizophrenia. We propose that this combination of neuroimaging and transcriptional data provides new insight into how previously implicated genes and proteins, as well as a number of unreported proteins in their vicinity on the protein interaction network, may interact to drive structural brain network changes in schizophrenia.

While the role of cortical microstructure in organising neural function is well established, it remains unclear how structural constraints can give rise to more flexible elements of cognition. While nonhuman primate research has demonstrated a close structure–function correspondence, the relationship between microstructure and function remains poorly understood in humans, in part because of the reliance on post mortem analyses, which cannot be directly related to functional data. To overcome this barrier, we developed a novel approach to model the similarity of microstructural profiles sampled in the direction of cortical columns. Our approach was initially formulated based on an ultra-high–resolution 3D histological reconstruction of an entire human brain and then translated to myelin-sensitive magnetic resonance imaging (MRI) data in a large cohort of healthy adults. This novel method identified a system-level gradient of microstructural differentiation traversing from primary sensory to limbic regions that followed shifts in laminar differentiation and cytoarchitectural complexity. Importantly, while microstructural and functional gradients described a similar hierarchy, they became increasingly dissociated in transmodal default mode and fronto–parietal networks. Meta-analytic decoding of these topographic dissociations highlighted involvement in higher-level aspects of cognition, such as cognitive control and social cognition. Our findings demonstrate a relative decoupling of macroscale functional from microstructural gradients in transmodal regions, which likely contributes to the flexible role these regions play in human cognition., Modelling cortico-cortical networks of microstructural similarity in post mortem and in vivo brains reveals consistent gradients of microstructural differentiation running from sensory to transmodal cortices, and allows single-subject investigation of structure-function correspondence.

BackgroundGenetic risk is thought to drive clinical variation on a spectrum of schizophrenia-like traits but the underlying changes in brain structure that mechanistically link genomic variation to schizotypal experience and behaviour are unclear.MethodsWe assessed schizotypy using a self-reported questionnaire, and measured magnetization transfer (MT), as a putative micro-structural MRI marker of intra-cortical myelination, in 68 brain regions, in 248 healthy young people (aged 14-25 years). We used normative adult brain gene expression data, and partial least squares (PLS) analysis, to find the weighted gene expression pattern that was most co-located with the cortical map of schizotypy-related magnetization (SRM).ResultsMagnetization was significantly correlated with schizotypy in bilateral posterior cingulate cortex and precuneus (and for disorganized schizotypy also in medial prefrontal cortex; all FDR-corrected P < 0.05), which are regions of the default mode network specialized for social and memory functions. The genes most positively weighted on the whole genome expression map co-located with SRM were enriched for genes that were significantly down-regulated in two prior case-control histological studies of brain gene expression in schizophrenia. Conversely, the most negatively weighted genes were enriched for genes that were transcriptionally up-regulated in schizophrenia. Positively weighted (down-regulated) genes were enriched for neuronal, specifically inter-neuronal, affiliations and coded a network of proteins comprising a few highly interactive “hubs” such as parvalbumin and calmodulin.ConclusionsMicrostructural MRI maps of intracortical magnetization can be linked to both the behavioural traits of schizotypy and to prior histological data on dysregulated gene expression in schizophrenia.

Childhood socioeconomic status (SES) impacts cognitive development and mental health, but its association with structural brain development is not yet well-characterized. Here, we analyzed 1243 longitudinally-acquired structural MRI scans from 623 youth to investigate the relation between SES and cortical and subcortical morphology between ages 5 and 25 years. We found positive associations between SES and total volumes of the brain, cortical sheet, and four separate subcortical structures. These associations were developmentally fixed rather than age-dependent. Surface-based shape analysis revealed that higher SES is associated with areal expansion of (i) lateral prefrontal, anterior cingulate, lateral temporal, and superior parietal cortices and (ii) ventrolateral thalamic, and medial amygdalo-hippocampal sub-regions. Meta-analyses of functional imaging data indicate that cortical correlates of SES are centered on brain systems subserving sensorimotor functions, language, memory, and emotional processing. We further show that anatomical variation within a subset of these cortical regions partially mediates the positive association between SES and IQ. Finally, we identify neuroanatomical correlates of SES that exist above and beyond accompanying variation in IQ. Our findings clarify the spatiotemporal patterning of SES-related neuroanatomical variation and inform ongoing efforts to dissect the causal pathways underpinning observed associations between childhood SES and regional brain anatomy.

Sex chromosome aneuploidy (SCA) enhances risk for several psychiatric disorders associated with the limbic system, including mood and autism spectrum disorders. These patients provide a powerful genetics-first model for understanding the biological basis of psychopathology. Additionally, these disorders are frequently sex-biased in prevalence, further suggesting an etiological role for sex chromosomes. To clarify how limbic anatomy varies across sex and sex chromosome complement, we characterized amygdala and hippocampus structure in a uniquely large sample of patients carrying supernumerary sex chromosomes (n = 132) and typically developing controls (n=166). After correction for sex-differences in brain size, karyotypically normal males (XY) and females (XX) did not differ in volume or shape of either structure. In contrast, all SCAs were associated with lowered amygdala volume relative to gonadally-matched controls. This effect was robust to three different methods for total brain volume correction, including an allometric analysis that derived normative scaling rules for these structures in a separate, typically developing population (n = 79). Hippocampal volume was insensitive to SCA after correction for total brain volume. However, surface-based analysis revealed that SCA, regardless of specific karyotype, was consistently associated with a spatially specific pattern of shape change in both amygdala and hippocampus. In particular, SCA was accompanied by contraction around the basomedial nucleus of the amygdala and an area within the hippocampal surface that cuts across hippocampal subfields. These results demonstrate the power of SCA as a model to understand how copy number variation can precipitate changes in brain systems relevant to psychiatric disease.

Understanding heterogeneity in neural phenotypes is an important goal on the path to precision medicine for autism spectrum disorders (ASD). Age is a critically important variable in normal structural brain development and examining structural features with respect to age-related norms could help to explain ASD heterogeneity in neural phenotypes. Here we examined how cortical thickness (CT) in ASD can be parameterized as an individualized metric of deviance relative to typically-developing (TD) age-related norms. Across a large sample (n=870 per group) and wide age range (5-40 years), we applied a normative modelling approach that provides individualized whole-brain maps of age-related CT deviance in ASD. This approach isolates a subgroup of ASD individuals with highly age-deviant CT. The median prevalence of this ASD subgroup across all brain regions is 7.6%, and can reach as high as 10% for some brain regions. Testing age-normed CT scores also highlights on-average differentiation, and associations with behavioural symptomatology that is separate from insights gleaned from traditional case-control approaches. This work showcases a novel individualized approach for understanding ASD heterogeneity that could potentially further prioritize work on a subset of individuals with significant cortical pathophysiology represented in age-related CT deviance. Rather than cortical thickness pathology being widespread characteristic of most ASD patients, only a small subset of ASD individuals are actually highly deviant relative to age-norms. These individuals drive a large majority of small effect results from canonical case-control comparisons and should be prioritized in future research to better understand the mechanisms behind highly deviant CT patterns.

Non-human primate neuroimaging is a rapidly growing area of research that promises to transform and scale translational and cross-species comparative neuroscience.Unfortunately, the technological and methodological advances of the past two decades have outpaced the accrual of data, which is particularly challenging given the relatively few centers that have the necessary facilities and capabilities. The PRIMate Data Exchange (PRIME-DE) addresses this challenge by aggregating independently acquired non-human primate magnetic resonance imaging (MRI) datasets and openly sharing them via the International Neuroimaging Data-sharing Initiative (INDI). Here, we present the rationale, design and procedures for the PRIME-DE consortium, as well as the initial release, consisting of 13 independent data collections aggregated across 11 sites (total = 98 macaque monkeys). We also outline the unique pitfalls and challenges that should be considered in the analysis of the non-human primate MRI datasets, including providing automated quality assessment of the contributed datasets.

Evolutionary and developmental increases in primate brain size have been accompanied by systematic shifts in the proportionality of different primate brain systems. However, it remains unknown if and how brain patterning varies across the more than 2-fold inter-individual variation in brain size that occurs amongst typically-developing humans. Using in vivo neuroimaging data from 2 independent cohorts totaling nearly 3000 individuals, we find that larger-brained humans show preferential areal expansion within specific fronto-parietal cortical networks (default mode, dorsal attentional) and related subcortical regions, at the expense of primary sensory/motor systems. This targeted areal expansion recapitulates cortical remodeling across evolution, manifests by early childhood and is linked to molecular signatures of heightened metabolic cost. Our results define a new organizing principle in human brain patterning which governs the highly-coordinated remodeling of human brain shape as a function of naturally-occurring variations in brain size.One Sentence SummaryA hodologically and metabolically expensive brain network is preferentially expanded in larger-brained humans.

Complex network topology is characteristic of many biological systems, including anatomical and functional brain networks (connectomes). Here, we first constructed a structural covariance network (SCN) from MRI measures of cortical thickness on 296 healthy volunteers, aged 14-24 years. Next, we designed a new algorithm for matching sample locations from the Allen Brain Atlas to the nodes of the SCN. Subsequently we use this to define, transcriptomic brain networks (TBN) by estimating gene co-expression between pairs of cortical regions. Finally, we explore the hypothesis that TBN and the SCN are coupled.TBN and SCN were correlated across connection weights and showed qualitatively similar complex topological properties. There were differences between networks in degree and distance distributions. However, cortical areas connected to each other within modules of the SCN network had significantly higher levels of whole genome co-expression than expected by chance.Nodes connected in the SCN had significantly higher levels of expression and co-expression of a Human Supragranular Enriched (HSE) gene set that are known to be important for large-scale cortico-cortical connectivity. This coupling of brain transcriptome and connectome topologies was largely but not completely related to the common constraint of physical distance on both networks.

SummaryMacroscopic cortical networks are important for cognitive function, but it remains challenging to construct anatomically plausible individual structural connectomes from human neuroimaging. We introduce a new technique for cortical network mapping, based on inter-regional similarity of multiple morphometric parameters measured using multimodal MRI. In three cohorts (two human, one macaque), we find that the resulting morphometric similarity networks (MSNs) have a complex topological organisation comprising modules and high-degree hubs. Human MSN modules recapitulate known cortical cytoarchitectonic divisions, and greater inter-regional morphometric similarity was associated with stronger inter-regional co-expression of genes enriched for neuronal terms. Comparing macaque MSNs to tract-tracing data confirmed that morphometric similarity was related to axonal connectivity. Finally, variation in the degree of human MSN nodes accounted for about 40% of between-subject variability in IQ. Morphometric similarity mapping provides a novel, robust and biologically plausible approach to understanding how human cortical networks underpin individual differences in psychological functions.

Bulimia nervosa (BN) is a serious psychiatric illness defined by preoccupation with weight and shape, episodic binge-eating and compensatory behaviors. Although diagnosed BN has been associated with diffuse grey matter volume reductions, characterization of brain structure alterations in women with a range of BN symptoms has yet to be made. This study examined whether changes in cortical thickness (CT) scaled with BN symptom severity in a sample of 33 adult women (n = 10 BN; n = 5 EDNOS-BN). Our second objective was to assess global structural connectivity (SC) of CT and to determine if individual differences in global SC relate to BN symptom severity. We used the validated Eating Disorder Examination Questionnaire (EDE-Q; Fairburn & Beglin, 1994) as a continuous measure of BN symptom severity. Increased EDE-Q score was negatively related to global CT and local CT in the left middle frontal gyrus, right superior frontal gyrus and bilateral orbitofrontal cortex (OFC) and temporoparietal regions. Moreover, analysis of global SC indicated that BN-related cortical thinning preferentially occurred in regions with high global connectivity. Finally, we showed that individuals’ contribution to global SC at the group level were significantly related to EDE-Q score, where increased EDE-Q score correlated with reduced connectivity of the left OFC and middle temporal cortex and increased connectivity of the right superior parietal lobule. Our findings offer novel insight into CT alterations in BN and further suggest that the combination of CT and structural connectivity measures may be sensitive to individual differences in BN symptom severity.