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.

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.

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.

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.

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.

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.

Interpretable surface-based detection of focal cortical dysplasias: a Multi-centre Epilepsy Lesion Detection study

One outstanding challenge for machine learning in diagnostic biomedical imaging is algorithm interpretability. A key application is the identification of subtle epileptogenic focal cortical dysplasias (FCDs) from structural MRI. FCDs are difficult to visualize on structural MRI but are often amenable to surgical resection. We aimed to develop an open-source, interpretable, surface-based machine-learning algorithm to automatically identify FCDs on heterogeneous structural MRI data from epilepsy surgery centres worldwide. The Multi-centre Epilepsy Lesion Detection (MELD) Project collated and harmonized a retrospective MRI cohort of 1015 participants, 618 patients with focal FCD-related epilepsy and 397 controls, from 22 epilepsy centres worldwide. We created a neural network for FCD detection based on 33 surface-based features. The network was trained and cross-validated on 50% of the total cohort and tested on the remaining 50% as well as on 2 independent test sites. Multidimensional feature analysis and integrated gradient saliencies were used to interrogate network performance. Our pipeline outputs individual patient reports, which identify the location of predicted lesions, alongside their imaging features and relative saliency to the classifier. On a restricted 'gold-standard' subcohort of seizure-free patients with FCD type IIB who had T1 and fluid-attenuated inversion recovery MRI data, the MELD FCD surface-based algorithm had a sensitivity of 85%. Across the entire withheld test cohort the sensitivity was 59% and specificity was 54%. After including a border zone around lesions, to account for uncertainty around the borders of manually delineated lesion masks, the sensitivity was 67%. This multicentre, multinational study with open access protocols and code has developed a robust and interpretable machine-learning algorithm for automated detection of focal cortical dysplasias, giving physicians greater confidence in the identification of subtle MRI lesions in individuals with epilepsy.

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.

Brain charts for the human lifespan

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

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.

Relating quantitative 7T MRI across cortical depths to cytoarchitectonics, gene expression and connectomics

Ultra-high field MRI across the depth of the cortex has the potential to provide anatomically precise biomarkers and mechanistic insights into neurodegenerative disease like Huntington's disease that show layer-selective vulnerability. Here we compare multi-parametric mapping (MPM) measures across cortical depths for a 7T 500 μm whole brain acquisition to (a) layer-specific cell measures from the von Economo histology atlas, (b) layer-specific gene expression, using the Allen Human Brain atlas and (c) white matter connections using high-fidelity diffusion tractography, at a 1.3 mm isotropic voxel resolution, from a 300mT/m Connectom MRI system. We show that R2*, but not R1, across cortical depths is highly correlated with layer-specific cell number and layer-specific gene expression. R1- and R2*-weighted connectivity strength of cortico-striatal and intra-hemispheric cortical white matter connections was highly correlated with grey matter R1 and R2* across cortical depths. Limitations of the layer-specific relationships demonstrated are at least in part related to the high cross-correlations of von Economo atlas cell counts and layer-specific gene expression across cortical layers. These findings demonstrate the potential and limitations of combining 7T MPMs, gene expression and white matter connections to provide an anatomically precise framework for tracking neurodegenerative disease.

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.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

Quantitative In vivo MRI Assessment of Structural Asymmetries and Sexual Dimorphism of Transient Fetal Compartments in the Human Brain

Structural asymmetries and sexual dimorphism of the human cerebral cortex have been identified in newborns, infants, children, adolescents, and adults. Some of these findings were linked with cognitive and neuropsychiatric disorders, which have roots in altered prenatal brain development. However, little is known about structural asymmetries or sexual dimorphism of transient fetal compartments that arise in utero. Thus, we aimed to identify structural asymmetries and sexual dimorphism in the volume of transient fetal compartments (cortical plate [CP] and subplate [SP]) across 22 regions. For this purpose, we used in vivo structural T2-weighted MRIs of 42 healthy fetuses (16.43-36.86 gestational weeks old, 15 females). We found significant leftward asymmetry in the volume of the CP and SP in the inferior frontal gyrus. The orbitofrontal cortex showed significant rightward asymmetry in the volume of CP merged with SP. Males had significantly larger volumes in regions belonging to limbic, occipital, and frontal lobes, which were driven by a significantly larger SP. Lastly, we did not observe sexual dimorphism in the growth trajectories of the CP or SP. In conclusion, these results support the hypothesis that structural asymmetries and sexual dimorphism in relative volumes of cortical regions are present during prenatal brain development.

CIVET-Macaque: An automated pipeline for MRI-based cortical surface generation and cortical thickness in macaques

The MNI CIVET pipeline for automated extraction of cortical surfaces and evaluation of cortical thickness from in-vivo human 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.

Transcriptomic and cellular decoding of regional brain vulnerability to neurogenetic disorders

Neurodevelopmental disorders have a heritable component and are associated with region specific alterations in brain anatomy. However, it is unclear how genetic risks for neurodevelopmental disorders are translated into spatially patterned brain vulnerabilities. Here, we integrated cortical neuroimaging data from patients with neurodevelopmental disorders caused by genomic copy number variations (CNVs) and gene expression data from healthy subjects. For each of the six investigated disorders, we show that spatial patterns of cortical anatomy changes in youth are correlated with cortical spatial expression of CNV genes in neurotypical adults. By transforming normative bulk-tissue cortical expression data into cell-type expression maps, we link anatomical change maps in each analysed disorder to specific cell classes as well as the CNV-region genes they express. Our findings reveal organizing principles that regulate the mapping of genetic risks onto regional brain changes in neurogenetic disorders. Our findings will enable screening for candidate molecular mechanisms from readily available neuroimaging data.

Planning stereoelectroencephalography using automated lesion detection: Retrospective feasibility study

Objective

This retrospective, cross‐sectional study evaluated the feasibility and potential benefits of incorporating deep‐learning on structural magnetic resonance imaging (MRI) into planning stereoelectroencephalography (sEEG) implantation in pediatric patients with diagnostically complex drug‐resistant epilepsy. This study aimed to assess the degree of colocalization between automated lesion detection and the seizure onset zone (SOZ) as assessed by sEEG.

Methods

A neural network classifier was applied to cortical features from MRI data from three cohorts. (1) The network was trained and cross‐validated using 34 patients with visible focal cortical dysplasias (FCDs). (2) Specificity was assessed in 20 pediatric healthy controls. (3) Feasibility of incorporation into sEEG implantation plans was evaluated in 34 sEEG patients. Coordinates of sEEG contacts were coregistered with classifier‐predicted lesions. sEEG contacts in seizure onset and irritative tissue were identified by clinical neurophysiologists. A distance of <10 mm between SOZ contacts and classifier‐predicted lesions was considered colocalization.

Results

In patients with radiologically defined lesions, classifier sensitivity was 74% (25/34 lesions detected). No clusters were detected in the controls (specificity = 100%). Of the total 34 sEEG patients, 21 patients had a focal cortical SOZ, of whom eight were histopathologically confirmed as having an FCD. The algorithm correctly detected seven of eight of these FCDs (86%). In patients with histopathologically heterogeneous focal cortical lesions, there was colocalization between classifier output and SOZ contacts in 62%. In three patients, the electroclinical profile was indicative of focal epilepsy, but no SOZ was localized on sEEG. In these patients, the classifier identified additional abnormalities that had not been implanted.

Significance

There was a high degree of colocalization between automated lesion detection and sEEG. We have created a framework for incorporation of deep‐learning–based MRI lesion detection into sEEG implantation planning. Our findings support the prospective evaluation of automated MRI analysis to plan optimal electrode trajectories.

Estimates of cortical column orientation improve MEG source inversion

Determining the anatomical source of brain activity non-invasively measured from EEG or MEG sensors is challenging. In order to simplify the source localization problem, many techniques introduce the assumption that current sources lie on the cortical surface. Another common assumption is that this current flow is orthogonal to the cortical surface, thereby approximating the orientation of cortical columns. However, it is not clear which cortical surface to use to define the current source locations, and normal vectors computed from a single cortical surface may not be the best approximation to the orientation of cortical columns. We compared three different surface location priors and five different approaches for estimating dipole vector orientation, both in simulations and visual and motor evoked MEG responses. We show that models with source locations on the white matter surface and using methods based on establishing correspondences between white matter and pial cortical surfaces dramatically outperform models with source locations on the pial or combined pial/white surfaces and which use methods based on the geometry of a single cortical surface in fitting evoked visual and motor responses. These methods can be easily implemented and adopted in most M/EEG analysis pipelines, with the potential to significantly improve source localization of evoked responses.

MRI profiling of focal cortical dysplasia using multi-compartment diffusion models

Objective

Focal cortical dysplasia (FCD) lesion detection and subtyping remain challenging on conventional MRI. New diffusion models such as the spherical mean technique (SMT) and neurite orientation dispersion and density imaging (NODDI) provide measurements that potentially produce more specific maps of abnormal tissue microstructure. This study aims to assess the SMT and NODDI maps for computational and radiological lesion characterization compared to standard fractional anisotropy (FA) and mean diffusivity (MD).

Methods

SMT, NODDI, FA, and MD maps were calculated for 33 pediatric patients with suspected FCD (18 histologically confirmed). Two neuroradiologists scored lesion visibility on clinical images and diffusion maps. Signal profile changes within lesions and homologous regions were quantified using a surface‐based approach. Diffusion parameter changes at multiple cortical depths were statistically compared between FCD type IIa and type IIb.

Results

Compared to fluid‐attenuated inversion recovery (FLAIR) or T1‐weighted imaging, lesions conspicuity on NODDI intracellular volume fraction (ICVF) maps was better/equal/worse in 5/14/14 patients, respectively, while on SMT intra‐neurite volume fraction (INVF) in 3/3/27. Compared to FA or MD, lesion conspicuity on the ICVF was better/equal/worse in 27/4/2, while on the INVF in 20/7/6. Quantitative signal profiling demonstrated significant ICVF and INVF reductions in the lesions, whereas SMT microscopic mean, radial, and axial diffusivities were significantly increased. FCD type IIb exhibited greater changes than FCD type IIa. No changes were detected on FA or MD profiles.

Significance

FCD lesion‐specific signal changes were found in ICVF and INVF but not in FA and MD maps. ICVF and INVF showed greater contrast than FLAIR in some cases and had consistent signal changes specific to FCD, suggesting that they could improve current presurgical pediatric epilepsy imaging protocols and can provide features useful for automated lesion detection.

Shifts in myeloarchitecture characterise adolescent development of cortical gradients

We studied an accelerated longitudinal cohort of adolescents and young adults (n = 234, 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) overall increases in MT, and ii) 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. Enhanced MT in mid-to-deeper layers was found to spatially co-localise specifically with gene expression markers of oligodendrocytes. Interregional covariance analysis revealed that these intracortical changes contributed to a gradual differentiation of higher-order from 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 organisation.

Microstructural and functional gradients are increasingly dissociated in transmodal cortices

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.

Mapping Cortical Laminar Structure in the 3D BigBrain

Histological sections offer high spatial resolution to examine laminar architecture of the human cerebral cortex; however, they are restricted by being 2D, hence only regions with sufficiently optimal cutting planes can be analyzed. Conversely, noninvasive neuroimaging approaches are whole brain but have relatively low resolution. Consequently, correct 3D cross-cortical patterns of laminar architecture have never been mapped in histological sections. We developed an automated technique to identify and analyze laminar structure within the high-resolution 3D histological BigBrain. We extracted white matter and pial surfaces, from which we derived histologically verified surfaces at the layer I/II boundary and within layer IV. Layer IV depth was strongly predicted by cortical curvature but varied between areas. This fully automated 3D laminar analysis is an important requirement for bridging high-resolution 2D cytoarchitecture and in vivo 3D neuroimaging. It lays the foundation for in-depth, whole-brain analyses of cortical layering.

Automated detection of focal cortical dysplasia type II with surface-based magnetic resonance imaging postprocessing and machine learning

OBJECTIVE

Focal cortical dysplasia (FCD) is a major pathology in patients undergoing surgical resection to treat pharmacoresistant epilepsy. Magnetic resonance imaging (MRI) postprocessing methods may provide essential help for detection of FCD. In this study, we utilized surface-based MRI morphometry and machine learning for automated lesion detection in a mixed cohort of patients with FCD type II from 3 different epilepsy centers.

METHODS

Sixty-one patients with pharmacoresistant epilepsy and histologically proven FCD type II were included in the study. The patients had been evaluated at 3 different epilepsy centers using 3 different MRI scanners. T1-volumetric sequence was used for postprocessing. A normal database was constructed with 120 healthy controls. We also included 35 healthy test controls and 15 disease test controls with histologically confirmed hippocampal sclerosis to assess specificity. Features were calculated and incorporated into a nonlinear neural network classifier, which was trained to identify lesional cluster. We optimized the threshold of the output probability map from the classifier by performing receiver operating characteristic (ROC) analyses. Success of detection was defined by overlap between the final cluster and the manual labeling. Performance was evaluated using k-fold cross-validation.

RESULTS

The threshold of 0.9 showed optimal sensitivity of 73.7% and specificity of 90.0%. The area under the curve for the ROC analysis was 0.75, which suggests a discriminative classifier. Sensitivity and specificity were not significantly different for patients from different centers, suggesting robustness of performance. Correct detection rate was significantly lower in patients with initially normal MRI than patients with unequivocally positive MRI. Subgroup analysis showed the size of the training group and normal control database impacted classifier performance.

SIGNIFICANCE

Automated surface-based MRI morphometry equipped with machine learning showed robust performance across cohorts from different centers and scanners. The proposed method may be a valuable tool to improve FCD detection in presurgical evaluation for patients with pharmacoresistant epilepsy.

Morphometric Similarity Networks Detect Microscale Cortical Organization and Predict Inter-Individual Cognitive Variation

Macroscopic 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 organization 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 with 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.

Gene transcription profiles associated with inter-modular hubs and connection distance in human functional magnetic resonance imaging networks

Human functional magnetic resonance imaging (fMRI) brain networks have a complex topology comprising integrative components, e.g. long-distance inter-modular edges, that are theoretically associated with higher biological cost. Here, we estimated intra-modular degree, inter-modular degree and connection distance for each of 285 cortical nodes in multi-echo fMRI data from 38 healthy adults. We used the multivariate technique of partial least squares (PLS) to reduce the dimensionality of the relationships between these three nodal network parameters and prior microarray data on regional expression of 20 737 genes. The first PLS component defined a transcriptional profile associated with high intra-modular degree and short connection distance, whereas the second PLS component was associated with high inter-modular degree and long connection distance. Nodes in superior and lateral cortex with high inter-modular degree and long connection distance had local transcriptional profiles enriched for oxidative metabolism and mitochondria, and for genes specific to supragranular layers of human cortex. In contrast, primary and secondary sensory cortical nodes in posterior cortex with high intra-modular degree and short connection distance had transcriptional profiles enriched for RNA translation and nuclear components. We conclude that, as predicted, topologically integrative hubs, mediating long-distance connections between modules, are more costly in terms of mitochondrial glucose metabolism.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Obesity associated with increased brain age from midlife

Common mechanisms in aging and obesity are hypothesized to increase susceptibility to neurodegeneration, however, direct evidence in support of this hypothesis is lacking. We therefore performed a cross-sectional analysis of magnetic resonance image-based brain structure on a population-based cohort of healthy adults. Study participants were originally part of the Cambridge Centre for Ageing and Neuroscience (Cam-CAN) and included 527 individuals aged 20-87 years. Cortical reconstruction techniques were used to generate measures of whole-brain cerebral white-matter volume, cortical thickness, and surface area. Results indicated that cerebral white-matter volume in overweight and obese individuals was associated with a greater degree of atrophy, with maximal effects in middle-age corresponding to an estimated increase of brain age of 10 years. There were no similar body mass index-related changes in cortical parameters. This study suggests that at a population level, obesity may increase the risk of neurodegeneration.

Multiple markers of cortical morphology reveal evidence of supragranular thinning in schizophrenia

In vivo structural neuroimaging can reliably identify changes to cortical morphology and its regional variation but cannot yet relate these changes to specific cortical layers. We propose, however, that by synthesizing principles of cortical organization, including relative contributions of different layers to sulcal and gyral thickness, regional patterns of variation in thickness of different layers across the cortical sheet and profiles of layer variation across functional hierarchies, it is possible to develop indirect morphological measures as markers of more specific cytoarchitectural changes. We developed four indirect measures sensitive to changes specifically occurring in supragranular cortical layers, and applied these to test the hypothesis that supragranular layers are disproportionately affected in schizophrenia. Our findings from the four different measures converge to indicate a predominance of supragranular thinning in schizophrenia, independent of medication and illness duration. We propose that these indirect measures offer novel ways of identifying layer-specific cortical changes, offering complementary in vivo observations to existing post-mortem studies.