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.

WONOEP 2022: Neurotechnology for the diagnosis of epilepsy.

Epilepsy's myriad causes and clinical presentations ensure that accurate diagnoses and targeted treatments remain a challenge. Advanced neurotechnologies are needed to better characterize individual patients across multiple modalities and analytical techniques. At the XVIth Workshop on Neurobiology of Epilepsy: Early Onset Epilepsies: Neurobiology and Novel Therapeutic Strategies (WONOEP 2022), the session on "advanced tools" highlighted a range of approaches, from molecular phenotyping of genetic epilepsy models and resected tissue samples to imaging-guided localization of epileptogenic tissue for surgical resection of focal malformations. These tools integrate cutting edge research, clinical data acquisition, and advanced computational methods to leverage the rich information contained within increasingly large datasets. A number of common challenges and opportunities emerged, including the need for multidisciplinary collaboration, multimodal integration, potential ethical challenges, and the multistage path to clinical translation. Despite these challenges, advanced epilepsy neurotechnologies offer the potential to improve our understanding of the underlying causes of epilepsy and our capacity to provide patient-specific treatment.

Cortical gene expression architecture links healthy neurodevelopment to the imaging, transcriptomics and genetics of autism and schizophrenia.

Human brain organization involves the coordinated expression of thousands of genes. For example, the first principal component (C1) of cortical transcription identifies a hierarchy from sensorimotor to association regions. In this study, optimized processing of the Allen Human Brain Atlas revealed two new components of cortical gene expression architecture, C2 and C3, which are distinctively enriched for neuronal, metabolic and immune processes, specific cell types and cytoarchitectonics, and genetic variants associated with intelligence. Using additional datasets (PsychENCODE, Allen Cell Atlas and BrainSpan), we found that C1-C3 represent generalizable transcriptional programs that are coordinated within cells and differentially phased during fetal and postnatal development. Autism spectrum disorder and schizophrenia were specifically associated with C1/C2 and C3, respectively, across neuroimaging, differential expression and genome-wide association studies. Evidence converged especially in support of C3 as a normative transcriptional program for adolescent brain development, which can lead to atypical supragranular cortical connectivity in people at high genetic risk for schizophrenia.

Long-term neuropsychological trajectories in children with epilepsy: does surgery halt decline?

Neuropsychological impairments are common in children with drug-resistant epilepsy. It has been proposed that epilepsy surgery may alleviate these impairments by providing seizure freedom; however, findings from prior studies have been inconsistent. We mapped long-term neuropsychological trajectories in children before and after undergoing epilepsy surgery, to measure the impact of disease course and surgery on functioning. We performed a retrospective cohort study of 882 children who had undergone epilepsy surgery at Great Ormond Street Hospital (1990-2018). We extracted patient information and neuropsychological functioning - obtained from IQ tests (domains: Full-Scale IQ, Verbal IQ, Performance IQ, Working Memory, and Processing Speed) and tests of academic attainment (Reading, Spelling and Numeracy) - and investigated changes in functioning using regression analyses. We identified 500 children (248 females) who had undergone epilepsy surgery (median age at surgery = 11.9 years, interquartile range = [7.8,15.0]) and neuropsychology assessment. These children showed declines in all domains of neuropsychological functioning in the time leading up to surgery (all p-values ≤ 0.001; e.g., ßFSIQ = -1.9, SEFSIQ = 0.3, pFSIQ < 0.001). Children lost on average one to four points per year, depending on the domain considered; 27-43% declined by 10 or more points from their first to their last preoperative assessment. At the time of presurgical evaluation, most children (46-60%) scored one or more standard deviations below the mean (<85) on the different neuropsychological domains; 37% of these met the threshold for intellectual disability (Full-Scale IQ < 70). On a group level, there was no change in performance from pre- to postoperative assessment on any of the domains (all p-values > 0.128). However, children who became seizure-free through surgery showed higher postoperative neuropsychological performance (e.g., rrb-FSIQ = 0.37, p < 0.001). These children continued to demonstrate improvements in neuropsychological functioning over the course of their long-term follow-up (e.g., ßFSIQ = 0.9, SEFSIQ = 0.3, pFSIQ = 0.004). Children who had discontinued antiseizure medication (ASM) treatment at one-year follow-up showed an eight-to-13-point advantage in postoperative Working Memory, Processing Speed, and Numeracy, and greater improvements in Verbal IQ, Working Memory, Reading, and Spelling (all p-values < 0.034) over the postoperative period compared to children who were seizure-free and still receiving ASMs. In conclusion, by providing seizure freedom and the opportunity for ASM cessation, epilepsy surgery may not only halt but reverse the downward trajectory that children with drug-resistant epilepsy display in neuropsychological functioning. To halt this decline as soon as possible, or potentially prevent it from occurring in the first place, children with focal epilepsy should be considered for epilepsy surgery as early as possible after diagnosis.

Combining OPM and lesion mapping data for epilepsy surgery planning: a simulation study.

When planning for epilepsy surgery, multiple potential sites for resection may be identified through anatomical imaging. Magnetoencephalography (MEG) using optically pumped sensors (OP-MEG) is a non-invasive functional neuroimaging technique which could be used to help identify the epileptogenic zone from these candidate regions. Here we test the utility of a-priori information from anatomical imaging for differentiating potential lesion sites with OP-MEG. We investigate a number of scenarios: whether to use rigid or flexible sensor arrays, with or without a-priori source information and with or without source modelling errors. We simulated OP-MEG recordings for 1309 potential lesion sites identified from anatomical images in the Multi-centre Epilepsy Lesion Detection (MELD) project. To localise the simulated data, we used three source inversion schemes: unconstrained, prior source locations at centre of the candidate sites, and prior source locations within a volume around the lesion location. We found that prior knowledge of the candidate lesion zones made the inversion robust to errors in sensor gain, orientation and even location. When the reconstruction was too highly restricted and the source assumptions were inaccurate, the utility of this a-priori information was undermined. Overall, we found that constraining the reconstruction to the region including and around the participant's potential lesion sites provided the best compromise of robustness against modelling or measurement error.

Transcriptional cartography integrates multiscale biology of the human cortex.

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Lesion detection in epilepsy surgery: Lessons from a prospective evaluation of a machine learning algorithm.

AIM: To evaluate a lesion detection algorithm designed to detect focal cortical dysplasia (FCD) in children undergoing stereoelectroencephalography (SEEG) as part of their presurgical evaluation for drug-resistant epilepsy. METHOD: This was a prospective, single-arm, interventional study (Idea, Development, Exploration, Assessment, and Long-Term Follow-Up phase 1/2a). After routine SEEG planning, structural magnetic resonance imaging sequences were run through an FCD lesion detection algorithm to identify putative clusters. If the top three clusters were not already sampled, up to three additional SEEG electrodes were added. The primary outcome measure was the proportion of patients who had additional electrode contacts in the SEEG-defined seizure-onset zone (SOZ). RESULTS: Twenty patients (median age 12 years, range 4-18 years) were enrolled, one of whom did not undergo SEEG. Additional electrode contacts were part of the SOZ in 1 out of 19 patients while 3 out of 19 patients had clusters that were part of the SOZ but they were already implanted. A total of 16 additional electrodes were implanted in nine patients and there were no adverse events from the additional electrodes. INTERPRETATION: We demonstrate early-stage prospective clinical validation of a machine learning lesion detection algorithm used to aid the identification of the SOZ in children undergoing SEEG. We share key lessons learnt from this evaluation and emphasize the importance of robust prospective evaluation before routine clinical adoption of such algorithms. WHAT THIS PAPER ADDS: The focal cortical dysplasia detection algorithm collocated with the seizure-onset zone (SOZ) in 4 out of 19 patients. The algorithm changed the resection boundaries in 1 of 19 patients undergoing stereoelectroencephalography for drug-resistant epilepsy. The patient with an altered resection due to the algorithm was seizure-free 1 year after resective surgery. Overall, the algorithm did not increase the proportion of patients in whom SOZ was identified.

The molecular genetic landscape of human brain size variation.

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity.

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.

Pediatric epilepsy surgery from 2000 to 2018: Changes in referral and surgical volumes, patient characteristics, genetic testing, and postsurgical outcomes.

OBJECTIVE: Neurosurgery is a safe and effective form of treatment for select children with drug-resistant epilepsy. Still, there is concern that it remains underutilized, and that seizure freedom rates have not improved over time. We investigated referral and surgical practices, patient characteristics, and postoperative outcomes over the past two decades. METHODS: We performed a retrospective cohort study of children referred for epilepsy surgery at a tertiary center between 2000 and 2018. We extracted information from medical records and analyzed temporal trends using regression analyses. RESULTS: A total of 1443 children were evaluated for surgery. Of these, 859 (402 females) underwent surgical resection or disconnection at a median age of 8.5 years (interquartile range [IQR] = 4.6-13.4). Excluding palliative procedures, 67% of patients were seizure-free and 15% were on no antiseizure medication (ASM) at 1-year follow-up. There was an annual increase in the number of referrals (7%, 95% confidence interval [CI] = 5.3-8.6; p < .001) and surgeries (4% [95% CI = 2.9-5.6], p < .001) over time. Duration of epilepsy and total number of different ASMs trialed from epilepsy onset to surgery were, however, unchanged, and continued to exceed guidelines. Seizure freedom rates were also unchanged overall but showed improvement (odds ratio [OR] 1.09, 95% CI = 1.01-1.18; p = .027) after adjustment for an observed increase in complex cases. Children who underwent surgery more recently were more likely to be off ASMs postoperatively (OR 1.04, 95% CI = 1.01-1.08; p = .013). There was a 17% annual increase (95% CI = 8.4-28.4, p < .001) in children identified to have a genetic cause of epilepsy, which was associated with poor outcome. SIGNIFICANCE: Children with drug-resistant epilepsy continue to be put forward for surgery late, despite national and international guidelines urging prompt referral. Seizure freedom rates have improved over the past decades, but only after adjustment for a concurrent increase in complex cases. Finally, genetic testing in epilepsy surgery patients has expanded considerably over time and shows promise in identifying patients in whom surgery is less likely to be successful.

Predicting seizure outcome after epilepsy surgery: Do we need more complex models, larger samples, or better data?

OBJECTIVE: The accurate prediction of seizure freedom after epilepsy surgery remains challenging. We investigated if (1) training more complex models, (2) recruiting larger sample sizes, or (3) using data-driven selection of clinical predictors would improve our ability to predict postoperative seizure outcome using clinical features. We also conducted the first substantial external validation of a machine learning model trained to predict postoperative seizure outcome. METHODS: We performed a retrospective cohort study of 797 children who had undergone resective or disconnective epilepsy surgery at a tertiary center. We extracted patient information from medical records and trained three models-a logistic regression, a multilayer perceptron, and an XGBoost model-to predict 1-year postoperative seizure outcome on our data set. We evaluated the performance of a recently published XGBoost model on the same patients. We further investigated the impact of sample size on model performance, using learning curve analysis to estimate performance at samples up to N = 2000. Finally, we examined the impact of predictor selection on model performance. RESULTS: Our logistic regression achieved an accuracy of 72% (95% confidence interval [CI] = 68%-75%, area under the curve [AUC] = .72), whereas our multilayer perceptron and XGBoost both achieved accuracies of 71% (95% CIMLP = 67%-74%, AUCMLP = .70; 95% CIXGBoost own = 68%-75%, AUCXGBoost own = .70). There was no significant difference in performance between our three models (all p > .4) and they all performed better than the external XGBoost, which achieved an accuracy of 63% (95% CI = 59%-67%, AUC = .62; pLR = .005, pMLP = .01, pXGBoost own = .01) on our data. All models showed improved performance with increasing sample size, but limited improvements beyond our current sample. The best model performance was achieved with data-driven feature selection. SIGNIFICANCE: We show that neither the deployment of complex machine learning models nor the assembly of thousands of patients alone is likely to generate significant improvements in our ability to predict postoperative seizure freedom. We instead propose that improved feature selection alongside collaboration, data standardization, and model sharing is required to advance the field.

Artificial intelligence for the detection of focal cortical dysplasia: Challenges in translating algorithms into clinical practice.

Focal cortical dysplasias (FCDs) are malformations of cortical development and one of the most common pathologies causing pharmacoresistant focal epilepsy. Resective neurosurgery yields high success rates, especially if the full extent of the lesion is correctly identified and completely removed. The visual assessment of magnetic resonance imaging does not pinpoint the FCD in 30%-50% of cases, and half of all patients with FCD are not amenable to epilepsy surgery, partly because the FCD could not be sufficiently localized. Computational approaches to FCD detection are an active area of research, benefitting from advancements in computer vision. Automatic FCD detection is a significant challenge and one of the first clinical grounds where the application of artificial intelligence may translate into an advance for patients' health. The emergence of new methods from the combination of health and computer sciences creates novel challenges. Imaging data need to be organized into structured, well-annotated datasets and combined with other clinical information, such as histopathological subtypes or neuroimaging characteristics. Algorithmic output, that is, model prediction, requires a technically correct evaluation with adequate metrics that are understandable and usable for clinicians. Publication of code and data is necessary to make research accessible and reproducible. This critical review introduces the field of automatic FCD detection, explaining underlying medical and technical concepts, highlighting its challenges and current limitations, and providing a perspective for a novel research environment.

A bipolar taxonomy of adult human brain sulcal morphology related to timing of fetal sulcation and trans-sulcal gene expression gradients.

We developed a computational pipeline (now provided as a resource) for measuring morphological similarity between cortical surface sulci to construct a sulcal phenotype network (SPN) from each magnetic resonance imaging (MRI) scan in an adult cohort (N=34,725; 45-82 years). Networks estimated from pairwise similarities of 40 sulci on 5 morphological metrics comprised two clusters of sulci, represented also by the bipolar distribution of sulci on a linear-to-complex dimension. Linear sulci were more heritable and typically located in unimodal cortex; complex sulci were less heritable and typically located in heteromodal cortex. Aligning these results with an independent fetal brain MRI cohort (N=228; 21-36 gestational weeks), we found that linear sulci formed earlier, and the earliest and latest-forming sulci had the least between-adult variation. Using high-resolution maps of cortical gene expression, we found that linear sulcation is mechanistically underpinned by trans-sulcal gene expression gradients enriched for developmental processes.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla.

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T2* has not been performed. Here we provide a dataset containing empirical T2* measurements acquired at 0.35 × 0.35 × 0.35 mm3 voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T2* (as well as R2*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T2*-weighted (f)MRI signal.

Networks Underlie Temporal Onset of Dysplasia-Related Epilepsy: A MELD Study.

OBJECTIVE: The purpose of this study was to evaluate if focal cortical dysplasia (FCD) co-localization to cortical functional networks is associated with the temporal distribution of epilepsy onset in FCD. METHODS: International (20 center), retrospective cohort from the Multi-Centre Epilepsy Lesion Detection (MELD) project. Patients included if >3 years old, had 3D pre-operative T1 magnetic resonance imaging (MRI; 1.5 or 3 T) with radiologic or histopathologic FCD after surgery. Images processed using the MELD protocol, masked with 3D regions-of-interest (ROI), and co-registered to fsaverage_sym (symmetric template). FCDs were then co-localized to 1 of 7 distributed functional cortical networks. Negative binomial regression evaluated effect of FCD size, network, histology, and sulcal depth on age of epilepsy onset. From this model, predictive age of epilepsy onset was calculated for each network. RESULTS: Three hundred eighty-eight patients had median age seizure onset 5 years (interquartile range [IQR] = 3-11 years), median age at pre-operative scan 18 years (IQR = 11-28 years). FCDs co-localized to the following networks: limbic (90), default mode (87), somatomotor (65), front parietal control (52), ventral attention (32), dorsal attention (31), and visual (31). Larger lesions were associated with younger age of onset (p = 0.01); age of epilepsy onset was associated with dominant network (p = 0.04) but not sulcal depth or histology. Sensorimotor networks had youngest onset; the limbic network had oldest age of onset (p values <0.05). INTERPRETATION: FCD co-localization to distributed functional cortical networks is associated with age of epilepsy onset: sensory neural networks (somatomotor and visual) with earlier onset, and limbic latest onset. These variations may reflect developmental differences in synaptic/white matter maturation or network activation and may provide a biological basis for age-dependent epilepsy onset expression. ANN NEUROL 2022;92:503-511.

IDEAL approach to the evaluation of machine learning technology in epilepsy surgery: protocol for the MAST trial.

Epilepsy and epilepsy surgery lend themselves well to the application of machine learning (ML) and artificial intelligence (AI) technologies. This is evidenced by the plethora of tools developed for applications such as seizure detection and analysis of imaging and electrophysiological data. However, few of these tools have been directly used to guide patient management. In recent years, the Idea, Development, Exploration, Assessment, Long-Term Follow-Up (IDEAL) collaboration has formalised stages for the evaluation of surgical innovation and medical devices, and, in many ways, this pragmatic framework is also applicable to ML/AI technology, balancing innovation and safety. In this protocol paper, we outline the preclinical (IDEAL stage 0) evaluation and the protocol for a prospective (IDEAL stage 1/2a) study to evaluate the utility of an ML lesion detection algorithm designed to detect focal cortical dysplasia from structural MRI, as an adjunct in the planning of stereoelectroencephalography trajectories in children undergoing intracranial evaluation for drug-resistant epilepsy.

Atlas of lesion locations and postsurgical seizure freedom in focal cortical dysplasia: A MELD study.

OBJECTIVE: 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. METHODS: The MELD (Multi-centre Epilepsy Lesion Detection) project collated a retrospective cohort of 580 patients with epilepsy attributed to FCD from 20 epilepsy centers worldwide. Magnetic resonance imaging-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. RESULTS: FCDs were nonuniformly distributed, concentrating in the superior frontal sulcus, frontal pole, and temporal pole. Epilepsy onset was typically before the age of 10 years. Earlier epilepsy onset was associated with lesions in primary sensory areas, whereas 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%. SIGNIFICANCE: FCD location is an important determinant of its size, the age at epilepsy onset, and the likelihood of seizure freedom postsurgery. 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 counseling in epilepsy.

The BigBrainWarp toolbox for integration of BigBrain 3D histology with multimodal neuroimaging.

Neuroimaging stands to benefit from emerging ultrahigh-resolution 3D histological atlases of the human brain; the first of which is 'BigBrain'. Here, we review recent methodological advances for the integration of BigBrain with multi-modal neuroimaging and introduce a toolbox, 'BigBrainWarp', that combines these developments. The aim of BigBrainWarp is to simplify workflows and support the adoption of best practices. This is accomplished with a simple wrapper function that allows users to easily map data between BigBrain and standard MRI spaces. The function automatically pulls specialised transformation procedures, based on ongoing research from a wide collaborative network of researchers. Additionally, the toolbox improves accessibility of histological information through dissemination of ready-to-use cytoarchitectural features. Finally, we demonstrate the utility of BigBrainWarp with three tutorials and discuss the potential of the toolbox to support multi-scale investigations of brain organisation.

Cortical patterning of morphometric similarity gradient reveals diverged hierarchical organization in sensory-motor cortices.

The topological organization of the cerebral cortex provides hierarchical axes, namely gradients, which reveal systematic variations of brain structure and function. However, the hierarchical organization of macroscopic brain morphology and how it constrains cortical function along the organizing axes remains unclear. We map the gradient of cortical morphometric similarity (MS) connectome, combining multiple features conceptualized as a "fingerprint" of an individual's brain. The principal MS gradient is anchored by motor and sensory cortices at two extreme ends, which are reliable and reproducible. Notably, divergences between motor and sensory hierarchies are consistent with the laminar histological thickness gradient but contrary to the canonical functional connectivity (FC) gradient. Moreover, the MS dissociates with FC gradients in the higher-order association cortices. The MS gradient recapitulates fundamental properties of cortical organization, from gene expression and cyto- and myeloarchitecture to evolutionary expansion. Collectively, our findings provide a heuristic hierarchical organization of cortical morphological neuromarkers.